Abstract visibility
Abstract text size

Papers for Wednesday, Feb 23 2022

Papers with local authors

Yi Feng, Di Li, Yuan-Pei Yang, Yongkun Zhang, Weiwei Zhu, Bing Zhang, Wenbin Lu, Pei Wang, Shi Dai, Ryan S. Lynch, Jumei Yao, Jinchen Jiang, Jiarui Niu, Dejiang Zhou, Heng Xu, Chenchen Miao, Chenhui Niu, Lingqi Meng, Lei Qian, Chao-Wei Tsai, Bojun Wang, Mengyao Xue, Youling Yue, Mao Yuan, Songbo Zhang, Lei Zhang
0 votes
Paper 19 — arXiv:2202.09601
0 votes
Paper 19 — arXiv:2202.09601

The polarization of fast radio bursts (FRBs), bright astronomical transients, contains crucial information about their environments. We report polarization measurements of five repeating FRBs, the abundant signals of which enable wide-band observations with two telescopes. A clear trend of lower polarization at lower frequencies was found, which can be well characterized by a single parameter rotation-measure-scatter ($\sigma_{\mathrm{RM}}$) and modeled by multi-path scatter. Sources with higher $\sigma_{\mathrm{RM}}$ have higher RM magnitude and scattering timescales. The two sources with the most substantial $\sigma_{\mathrm{RM}}$, FRB 20121102A and FRB 20190520B, are associated with a compact persistent radio source. These properties indicate a complex environment near the repeating FRBs, such as a supernova remnant or a pulsar wind nebula, consistent with their arising from young populations.

Yuan-Pei Yang, Wenbin Lu, Yi Feng, Bing Zhang, Di Li

10 pages, 1 figures, submitted to AAS journals on Feb. 7, 2022, comments welcome!

0 votes
Paper 20 — arXiv:2202.09602
0 votes
Paper 20 — arXiv:2202.09602

Some repeating fast radio burst (FRB) sources exhibit complex polarization behaviors, including frequency-dependent depolarization, variation of rotation measure (RM), and oscillating spectral structures of polarized components. Very recently, Feng et al. (2022) reported that active repeaters exhibit conspicuous frequency-dependent depolarization and a strong correlation between RM scatter ($\sigma_{\rm RM}$) and the temporal scattering time ($\tau_{\rm s}$), $\sigma_{\rm RM}\propto\tau_{\rm s}^{1.0\pm0.2}$, both of which can be well described by multi-path propagation through a magnetized inhomogeneous plasma screen. This observation strongly suggests that the temporal scattering and RM scatter originate from the same region. Besides, a particular finding of note in Feng et al. (2022) is that the FRBs with compact persistent radio sources (PRS) tend to have extreme $\sigma_{\rm RM}$. In this work, we analyze the temporal scattering, RM scatter and the PRS emission contributed by the magnetized inhomogeneous plasma screen near an FRB source. The behaviors of the RM scatter imply that the magnetized plasma environment is consistent with a supernova remnant or a pulsar wind nebula, and the predicted $\sigma_{\rm RM}$-$\tau_{\rm s}$ relation is $\sigma_{\rm RM}\propto\tau_{\rm s}^{(0.54-0.83)}$ for different astrophysical scenarios. We also show that the specific luminosity of a PRS should have a positive correlation with the RM contributed by the plasma screen. This is consistent with the observations of FRB 121102 and FRB 190520B.

Thomas Mikal-Evans, David K. Sing, Joanna K. Barstow, Tiffany Kataria, Jayesh Goyal, Nikole Lewis, Jake Taylor, Nathan. J. Mayne, Tansu Daylan, Hannah R. Wakeford, Mark S. Marley, Jessica J. Spake

Published in Nature Astronomy (publisher version is Open Access)

0 votes
Paper 32 — arXiv:2202.09884
0 votes
Paper 32 — arXiv:2202.09884

The temperature profile of a planetary atmosphere is a key diagnostic of radiative and dynamical processes governing the absorption, redistribution, and emission of energy. Observations have revealed dayside stratospheres that either cool or warm with altitude for a small number of gas giant exoplanets, while other dayside stratospheres are consistent with constant temperatures. Here we report spectroscopic phase curve measurements for the gas giant WASP-121b, which constrain stratospheric temperatures throughout the diurnal cycle. Variations measured for a water vapour spectral feature reveal a temperature profile that transitions from warming with altitude on the dayside hemisphere to cooling with altitude on the nightside hemisphere. The data are well explained by models assuming chemical equilibrium, with water molecules thermally dissociating at low pressures on the dayside and recombining on the nightside. Nightside temperatures are low enough for perovskite (CaTiO3) to condense, which could deplete titanium from the gas phase and explain recent non-detections at the day-night terminator. Nightside temperatures are also consistent with the condensation of refractory species such as magnesium, iron, and vanadium. Detections of these metals at the day-night terminator suggest, however, that if they do form nightside clouds, cold trapping does not efficiently remove them from the upper atmosphere. Horizontal winds and vertical mixing could keep these refractory condensates aloft in the upper atmosphere of the nightside hemisphere until they are recirculated to the hotter dayside hemisphere and vaporised.

Changseok Kim, Jong-Hak Woo, Yashashree Jadhav, Aeree Chung, Junhyun Baek, Jeong Ae Lee, Jaejin Shin, Ho Seong Hwang, Rongxin Luo, Donghoon Son, HyunGi Kim, Hyuk Woo

Accepted for publication in ApJ (Feb. 9th, 2022). 23 pages, 15 figures

0 votes
Paper 49 — arXiv:2202.10044
0 votes
Paper 49 — arXiv:2202.10044

We present the star formation rate (SFR) measurements based on the spectral energy distribution (SED) analysis with new sub-mm fluxes combined with archival multi-wavelength data for a sample of 52 AGN host galaxies at z $<0.2$. We carried out sub-mm observations using the SCUBA-2 camera at the James Clerk Maxwell Telescope, and obtained flux or an upper limit at 450 and 850 $\mu$m for each target. By experimenting the effect of the AGN dust component in the SED fit, we find that dust luminosity can be overestimated if AGN contribution is ignored. While the SFR based on 4000{\AA} break shows a significant offset compared to dust luminosity based SFR, the SFR obtained by the artificial neural network (Ellison et al. 2016) generally shows consistency albeit with a large scatter. We find that SFR correlates with AGN outflow strength manifested by the [OIII] $\lambda 5007$ emission line, and that AGNs with higher Eddington ratios and stronger outflows are in general hosted by galaxies with higher SFR, which is consistent with the correlation reported by Woo et al. (2020). This suggests no instantaneous quenching of star formation due to AGN feedback.

Nitsan Bar, Shany Danieli, Kfir Blum

12 pages, 6 figures + appendices

0 votes
Paper 59 — arXiv:2202.10179
0 votes
Paper 59 — arXiv:2202.10179

Ultra-diffuse galaxies that contain a large sample of globular clusters (GCs) offer an opportunity to test the predictions of galactic dynamics theory. NGC5846-UDG1 is an excellent example, with a high-quality sample of dozens of GC candidates. We show that the observed distribution of GCs in NGC5846-UDG1 is suggestive of mass segregation induced by gravitational dynamical friction. We present simple analytic calculations, backed by a series of numerical simulations, that naturally explain the observed present-day pattern of GC masses and radial positions. Subject to some assumptions on the GC population at birth, the analysis supports the possibility that NGC5846-UDG1 resides in a massive dark matter halo. This is an example for the use of GC-rich systems as dynamical (in addition to kinematical) tracers of dark matter.

All other papers

Zac Johnston, Sheldon Wasik, Rachel Titus, MacKenzie L. Warren, Evan P. O'Connor, Remco Zegers, Sean M. Couch

14 pages, 8 figures, 1 table. Submitted to ApJ

Recent studies have highlighted the sensitivity of core-collapse supernovae (CCSNe) models to electron-capture (EC) rates on neutron-rich nuclei near the N=50 closed-shell region. In this work, we perform a large suite of one-dimensional (1D) CCSN simulations for 200 stellar progenitors using recently-updated EC rates in this region. For comparison, we repeat the simulations using two previous implementations of EC rates: a microphysical library with parametrized N=50 rates (LMP), and an older single-nucleus approximation (SNA). We follow the simulations through shock revival up to several seconds post-bounce, and show that the EC rates produce a consistent imprint on CCSN properties, often surpassing the role of the progenitor itself. Notable impacts include the timescale of core collapse, the electron fraction and mass of the inner core at bounce, the accretion rate through the shock, the success or failure of revival, and the properties of the central compact remnant. We also compare the observable neutrino signal of the neutronization burst in a DUNE-like detector, and find consistent impacts to the counts and mean energies. Overall, the updated rates result in properties that are intermediate between LMP and SNA, and yet slightly more favorable to explosion than both.

Andrea V. Macciò, Mohamad Ali-Dib, Pavle Vulanović, Hind Al Noori, Fabian Walter, Nico Krieger, Tobias Buck

8 pages, 5 figures, accepted for publication on MNRAS

Cosmological galaxy formation simulations are still limited by their spatial/mass resolution and cannot model from first principles some of the processes, like star formation, that are key in driving galaxy evolution. As a consequence they still rely on a set of 'effective parameters' that try to capture the scales and the physical processes that cannot be directly resolved in the simulation. In this study we show that it is possible to use Machine Learning techniques applied to real and simulated images of galaxies to discriminate between different values of these parameters by making use of the full information content of an astronomical image instead of collapsing it into a limited set of values like size, or stellar/ gas masses. In this work we apply our method to the NIHAO simulations and the THINGS and VLA-ANGST observations of HI maps in nearby galaxies to test the ability of different values of the star formation density threshold $n$ to reproduce observed HI maps. We show that observations indicate the need for a high value of $n \gtrsim 80$ ,cm$^{-3}$ (although the exact numerical value is model-dependent), which has important consequences for the dark matter distribution in galaxies. Our study shows that with innovative methods it is possible to take full advantage of the information content of galaxy images and compare simulations and observations in an interpretable, non-parametric and quantitative manner.

U. Meštrić, E. Vanzella, A. Zanella, M. Castellano, F. Calura, P. Rosati, P. Bergamini, A. Mercurio, M. Meneghetti, C. Grillo, G.B. Caminha, M. Nonino, E. Merlin, G. Cupani, E. Sani

22 pages, 18 figures, 1 table, submitted to MNRAS. Comments welcome

We study the physical properties (size, stellar mass, luminosity, star formation rate) and scaling relations for a sample of 166 star-forming clumps with redshift $z \sim 2-6.2$. They are magnified by the Hubble Frontier Field galaxy cluster MACS~J0416 and have robust lensing amplification ($2\lesssim \mu \lesssim 82$) computed by using our high-precision lens model, based on 182 multiple images. Our sample extends by $\sim 3$ times the number of spectroscopically-confirmed lensed clumps at $z \gtrsim 2$. We identify clumps in ultraviolet continuum images and find that, whenever the effective spatial resolution (enhanced by gravitational lensing) increases, they fragment into smaller entities, likely reflecting the hierarchically-organized nature of star formation. Kpc-scale clumps, most commonly observed in non-lensed fields, are not found in our sample. The physical properties of our sample extend the parameter space typically probed by $z \gtrsim 1$ non-lensed observations and simulations, by populating the low mass (M$_\star \lesssim 10^7$ M$_\odot$), low star formation rate (SFR $\lesssim 0.5$ M$_\odot$ yr$^{-1}$), and small size (R$_\mathrm{eff} \lesssim 100$ pc) regime. The new domain probed by our study approaches the regime of compact stellar complexes and star clusters. In the mass-size plane our sample span the region between galaxies and globular clusters, with a few clumps laying in the region populated by young star clusters and globular-clusters. For the bulk of our sample, we measure star-formation rates which are higher than those observed locally in compact stellar systems, indicating different conditions for star formation at high redshift and in the local Universe.

Reconstructing lens potentials and lensed sources can easily become an underconstrained problem, even when the degrees of freedom are low, due to degeneracies, particularly when potential perturbations superimposed on a smooth lens are included. Regularization has traditionally been used to constrain the solutions where the data failed to do so, e.g. in unlensed parts of the source. In this exploratory work, we go beyond the usual choices of regularization and adopt observationally motivated priors for the source brightness. We also perform a similar comparison when reconstructing lens potential perturbations, which are assumed to be stationary, i.e. permeate the entire field of view. We find that physically motivated priors lead to lower residuals, avoid overfitting, and are decisively preferred within a Bayesian quantitative framework in all the examples considered. For the perturbations, choosing the wrong regularization can have a detrimental effect that even high-quality data cannot correct for, while using a purely smooth lens model can absorb them to a very high degree and lead to biased solutions. Finally, our new implementation of the semi-linear inversion technique provides the first quantitative framework for measuring degeneracies between the source and the potential perturbations.

Tristan L. Smith, Matteo Lucca, Vivian Poulin, Guillermo F. Abellan, Lennart Balkenhol, Karim Benabed, Silvia Galli, Riccardo Murgia

19 pages, 9 figures

We investigate constraints on early dark energy (EDE) using ACT DR4, SPT-3G 2018, Planck polarization, and restricted Planck temperature data (at $\ell<650$), finding a $3.3\sigma$ preference ($\Delta\chi^2=-16.2$ for 3 additional degrees of freedom) for EDE over $\Lambda$CDM. The EDE contributes a maximum fractional energy density of $f_{\rm EDE}(z_c)=0.163^{+0.047}_{-0.04}$ at a redshift $z_c=3357\pm200$ and leads to a CMB inferred value of the Hubble constant $H_0=74.2^{+1.9}_{-2.1}$ km/s/Mpc. We find that Planck and ACT DR4 data provide the majority of the improvement in $\chi^2$, and that the inclusion of SPT-3G pulls the posterior of $f_{\rm EDE}(z_c)$ away from $\Lambda$CDM. This is the first time that a moderate preference for EDE has been reported for these three combined CMB data sets. We find that including measurements of supernovae luminosity distances and the baryon acoustic oscillation standard ruler only minimally affects the preference ($3.0\sigma$), while measurements that probe the clustering of matter at late times - the lensing potential power spectrum from Planck and $f \sigma_8$ from BOSS - decrease the significance of the preference to 2.6$\sigma$. Conversely, adding a prior on the $H_0$ value as reported by the SH0ES collaboration increases the preference to the $4-5\sigma$ level. In the absence of this prior, the inclusion of Planck TT data at $\ell>1300$ reduces the preference from $3.0\sigma$ to $2.3\sigma$ and the constraint on $f_{\rm EDE}(z_c)$ becomes compatible with $\Lambda$CDM at $1\sigma$. We explore whether systematic errors in the Planck polarization data may affect our conclusions and find that changing the TE polarization efficiencies significantly reduces the Planck preference for EDE. More work will be necessary to establish whether these hints for EDE within CMB data alone are the sole results of systematic errors or an opening to new physics.

Adam Coogan, Thomas D. P. Edwards, Horng Sheng Chia, Richard N. George, Katherine Freese, Cody Messick, Christian N. Setzer, Christoph Weniger, Aaron Zimmerman

15 pages, 6 figures. Comments welcome! Code can be found at this https URL

The most sensitive search pipelines for gravitational waves from compact binary mergers use matched filters to extract signals from the noisy data stream coming from gravitational wave detectors. Matched-filter searches require banks of template waveforms covering the physical parameter space of the binary system. Unfortunately, template bank construction can be a time-consuming task. Here we present a new method for efficiently generating template banks that utilizes automatic differentiation to calculate the parameter space metric. Principally, we demonstrate that automatic differentiation enables accurate computation of the metric for waveforms currently used in search pipelines, whilst being computationally cheap. Additionally, by combining random template placement and a Monte Carlo method for evaluating the fraction of the parameter space that is currently covered, we show that search-ready template banks for frequency-domain waveforms can be rapidly generated. Finally, we argue that differentiable waveforms offer a pathway to accelerating stochastic placement algorithms. We implement all our methods into an easy-to-use Python package based on the jax framework, diffbank, to allow the community to easily take advantage of differentiable waveforms for future searches.

Benjamin Gregg, Daniela Calzetti, Mark Heyer

26 pages, 11 figures

We present an extensive archival analysis of a sample of local galaxies, combining multi-wavelength data from GALEX, Spitzer and Herschel to investigate "blue-side" mid-infrared (MIR) and "red-side" far-infrared (FIR) color-color correlations within the observed infrared spectral energy distributions (IR SEDs). Our sample largely consists of the KINGFISH galaxies, with the important addition of a select few including NGC5236 (M83) and NGC4449. With data from the far-ultraviolet FUV (${\sim}0.15$ $\mu$m) through 500 $\mu$m convolved to common angular resolution, we measure photometry of $kpc$-scale star-forming regions 36$"\times$36$"$ in size. Star formation rates (SFRs), stellar masses and metallicity distributions are derived throughout our sample. Focusing on the $f_{70}/f_{500}$ "FIR" and $f_{8}/f_{24}$ "MIR" flux density ratios (colors), we find that a sub-sample of galaxies demonstrate a strong IR color-color correlation within their star-forming regions, while others demonstrate uncorrelated colors. This division is driven by two main effects: 1) the local strength of star formation (SF) and 2) the metal content of the interstellar medium (ISM). Galaxies uniformly dominated by high surface densities of SF (e.g. NGC5236) demonstrate strong IR color-color correlations, while galaxies that exhibit lower levels of SF and mixed environments (e.g. NGC5457) demonstrate weaker or no correlation--explained by the increasing effect of varying ISM heating and metal content on the IR colors, specifically in the MIR. We find large dispersion in the SFR-$L_{8}$ (8 $\mu$m luminosity) relation that is traced by the metallicity distributions, consistent with extant studies, highlighting its problematic use as a SFR indicator across diverse systems/samples.

D. Deras, A. Arellano Ferro, I. Bustos Fierro, M. A. Yepez

13 pages, 8 figures, To be published by Revista Mexicana de Astronomia y Astrofisica

We present a $VI$ CCD photometric study and a membership analysis of the globular cluster M56 (NGC 6779). This produced a decontaminated CMD from field stars, which enabled a better confrontation with theoretical isochrones; Zero-Age Horizontal Branches (ZAHB) and post-ZAHB evolutionary tracks. Post He-flash evolutionary models with a He-core mass of 0.5 M$_\odot$ and envelopes of 0.04 - 0.18 M$_\odot$, cover the complete Horizontal Branch. Models with total mass $\sim$ 0.68 M$_\odot$ explain the RR Lyrae, while those with a mass $\sim$ 0.56 M$_\odot$ and a very subtle envelope, explain the Pop II cepheids with a progenitor in the blue tail of the HB. %Using the mean $(V-I)$ at minimum color curve for the RRab stars, we derived a value for the reddening of $E(B-V)$ = 0.26 $\pm$ 0.02. Based on the Fourier decomposition of the $V$ light curve of a single cluster member RRc star, we determined a metallicity of $\rm [Fe/H]_{ ZW}$ = -1.96 $\pm$ 0.09. Several independent distance determination approaches lead to a mean distance to M56 of $\big< d \big>$ = 9.4 $\pm$ 0.4 kpc. Finally, we report 5 new variables: 1 SX Phe, 3 EB, and 1 RRc.

Rishabh Singh Teja, Avinash Singh, D.K. Sahu, G.C. Anupama, Brajesh Kumar, A.J. Nayana

25 pages, 20 figures, and 7 tables. Submitted to the AAS Journals

We present spectroscopic and photometric observations of the Type IIP supernova, SN$~$2020jfo, in ultraviolet and optical wavelengths. SN$~$2020jfo occurred in the spiral galaxy M61 (NGC$\,$4303), with eight observed supernovae in the past 100 years. SN$~$2020jfo exhibited a short plateau lasting $58\pm2\,d$, and achieved a maximum brightness in $V$-band of $\rm M_V=-17.4\pm0.4\,mag$ at about $\rm 7.7\pm0.5\,d$ since explosion. From the bolometric light curve, we have estimated the mass of $\rm ^{56}Ni$ synthesised in the explosion to be $\rm 0.033\pm0.006\,M_\odot$. The observed spectral features are typical for a type IIP supernova except for shallow H$\alpha$ absorption throughout the evolution and the presence of stable $\rm ^{58}$Ni feature at 7378$\,$\r{A}, in the nebular phase. Using hydrodynamical modelling in the MESA$\,+\,$STELLA framework, we have estimated an ejecta mass of $\rm \sim 5\,M_\odot$. Models also indicate SN$~$2020jfo could be the result of a Red Super Giant progenitor with $\rm M_{ZAMS}\,\sim\,12\,M_\odot$. Bolometric light curve modelling revealed the presence of a secondary radiation source for initial $\rm \sim 20\,d$, which has been attributed to interaction with a circumstellar material of mass $\rm \sim0.2\,M_\odot$, which most likely was ejected due to enhanced mass loss about 20 years prior to the supernova explosion.

Kendall Sullivan, Adam L. Kraus

28 pages, 11 figures, accepted for publication in ApJ

Accretion is one of the defining characteristics of classical T Tauri stars, fueled by the presence of a circumstellar disk comprised of dust and gas. Accretion produces a UV and optical excess, while re-radiated emission at the inner edge of the dust component of the disk produces a near-infrared (NIR) excess. The interplay between stars and their disks helps regulate protoplanetary disk evolution and dispersal, which is key to a full understanding of planet formation. To investigate the relations between NIR excess and optical excess in both single and binary stars, we used an archival sample of spectroscopically characterized members of the Taurus star-forming region ($\tau \sim$ 1-2 Myr) with measured luminosities, spectral types, and optical veiling. We combined the archival sample with 2MASS and WISE NIR photometry and high-resolution imaging surveys. We found that NIR and optical excesses are correlated in multiple NIR photometric bands, suggesting that they are closely related, likely because more massive disks have higher inner dust disk walls and are also associated with higher accretion rates. We also found that multiplicity has no impact on accretion or inner disk properties in a sample with a wide range of separations, but the sample was too small to specifically investigate close binaries, where the effects of multiplicity on disk properties should be most significant.

Demet Kirmizibayrak, Jeremy Heyl

Accepted for publication in International Astronomical Union Proceedings Series No. 363, 2021 : Neutron Star Astrophysics at the Crossroads: Magnetars and the Multimessenger Revolution

We present our findings on magnetar spectral line analysis in the context of upcoming high resolution, high effective area, high throughput X-ray telescopes for two cases: persistent magnetar emission and magnetar bursts. For magnetars in quiescence, we present our preliminary work on modelling for phase-resolved emission. Our results reveal the necessity of constraining line depth and width concurrently with line energy to conclusively determine hotspot emission and corresponding magnetic field geometry. We then present the results of our simulations using effective area and response of various current and upcoming X-ray telescopes for magnetar spectral line detection and expand on the exciting opportunities upcoming telescopes provide to probe quiescent and burst emission region geometry and propagation in the extreme magnetic field of a magnetar

Teng Liu, Andrea Merloni, Julien Wolf, Mara Salvato, Thomas Reiprich, Riccardo Arcodia, Georg Lamer, Antonis Georgakakis, Tom Dwelly, Jeremy Sanders, Johannes Buchner, Frank Haberl, Miriam Ramos-Ceja, Joern Wilms, Kirpal Nandra, Hermann Brunner, Marcella Brusa, Axel Schwope, Jan Robrade, Michael J. Freyberg, Thomas Boller, Chandreyee Maitra, Angie Veronica, Adam Malyali

Submitted to A&A

Context: The eROSITA X-ray telescope on board the SRG observatory performed Calibration and Performance Verification (CalPV) observations between September 2019 and December 2019, ahead of the planned four-year all-sky surveys. Most of them were deep, pointing-mode observations. Aims: We present here the X-ray catalog detected from the set of extra-galactic CalPV observations released to the public by the German eROSITA consortium, and the multi-band counterparts of these X-ray sources. Methods: We develop a source detection method optimized for point-like X-ray sources by including extended X-ray emission in the background measurement. The multi-band counterparts are identified using a Bayesian method from the CatWISE catalog. Results: Combining 11 CalPV fields, we present a catalog containing 9522 X-ray sources, whose X-ray fluxes are measured through spectral fitting. CatWISE counterparts are presented for 77% of the sources. Significant variabilities are found in 179 of the sources, which are also presented with this paper. Most of these fields show similar number counts of point sources as typical extragalactic fields, and a few harbor particular stellar populations.

Vikram V. Dwarkadas (University of Chicago)

3 figures, 1 table. 13 pages. Author version of published paper, with reduced size figures 2 and 3. For full resolution figures see the published paper at this https URL Accepted to Galaxies. This article belongs to the special issue of the conference Asymmetrical Post Main-Sequence Nebulae 8

Using a code that employs a self-consistent method for computing the effects of photo-ionization on circumstellar gas dynamics, we model the formation of wind-driven nebulae around massive stars. We take into account changes in stellar properties and mass-loss over the star's evolution. Our simulations show how various properties, such as the density and ionization fraction, change throughout the evolution of the star. The multi-dimensional simulations reveal the presence of strong ionization front instabilities in the main-sequence phase, similar to those seen in galactic ionization fronts. Hydrodynamic instabilities at the interfaces lead to the formation of filaments and clumps that are continually being stripped off and mixed with the low density interior. Even though the winds start out as completely radial, the spherical symmetry is quickly destroyed, and the shocked wind region is manifestly asymmetrical. The simulations demonstrate that it is important to include the effects of the photoionizing photons from the star, and simulations that do not include this may fail to reproduce the observed density profile and ionization structure of wind-blown bubbles around massive stars.

Bachelet Etienne, Specht David, Penny Matthew, Hundertmark Markus, Awiphan Supachai, Beaulieu Jean-Philippe, Dominik Martin, Kerins Eamonn, Maoz Dan, Meade Evan, Nucita Achille, Poleski Radek, Ranc Clement, Rhodes Jason, Robin Annie

Accepted in A&A

As the Kepler mission has done for hot exoplanets, the ESA Euclid and NASA Roman missions have the potential to create a breakthrough in our understanding of the demographics of cool exoplanets, including unbound, or "free-floating", planets (FFPs). In this study, we demonstrate the complementarity of the two missions and propose two joint-surveys to better constrain the mass and distance of microlensing events. We first demonstrate that an early brief Euclid survey (7 h) of the Roman microlensing fields will allow the measurement of a large fraction of events relative proper motions and lens magnitudes. Then, we study the potential of simultaneous observations by Roman and Euclid to enable the measurement of the microlensing parallax for the shortest microlensing events. Using detailed simulations of the joint detection yield we show that within one year Roman-Euclid observations will be at least an order of magnitude more sensitive than current ground-based measurements. Depending on the exact distribution of FFP, a joint Roman-Euclid campaign should detect around 130 FFP events within a year, including 110 with measured parallax that strongly constrain the FFP mass, and around 30 FFP events with direct mass and distance measurements. The ability of the joint survey to completely break the microlens mass-distance-velocity degeneracy for a significant subset of events provides a unique opportunity to verify unambiguously the FFP hypothesis or else place abundance limits for FFPs between Earth and Jupiter masses that are up to two orders of magnitude stronger than provided by ground-based surveys. Finally, we study the capabilities of the joint survey to enhance the detection and charcterization of exomoons, and found that it could lead to the detection of the first exomoon.

David Stevenson, Peter Bodenheimer, Jack J. Lissauer, Gennaro D'Angelo

20 pages, 1 table, 8 figures. Submitted to The Planetary Science Journal of the American Astronomical Society

Simulations of Jupiter's formation are presented that incorporate mixing of H-He with denser material entering the planet as solids. Heavy compounds and gas mix substantially when the planet becomes roughly as massive as Earth, because incoming planetesimals can fully vaporize. Supersaturation of vaporized silicates causes the excess to sink as droplets, but water remains at higher altitudes. Because the mean molecular weight decreases rapidly outward, some of the compositional inhomogeneities produced during formation can survive for billions of years. After 4.57 Gyr, our Jupiter model retains compositional gradients; proceeding outwards one finds: i) an inner heavy-element core, the outer part derived from hot supersaturated rain-out; ii) a composition-gradient region, containing most of the heavy elements, where H-He abundance increases outward, reaching about 0.9 mass fraction at 0.3 of the radius, with silicates enhanced relative to water in the lower parts and depleted in the upper parts; iii) a uniform composition region (neglecting He immiscibility) that is enriched over protosolar and contains most of the planet's mass; and iv) an outer region where cloud formation (condensation) of heavy constituents occurs. This radial compositional profile has heavy elements more broadly distributed than predicted by classical formation models, but less diluted than suggested by Juno-constrained gravity models. The compositional gradients in the region containing the bulk of the heavy elements prevent convection, in both our models and those fitting current gravity, resulting in a hot interior where much of the accretion energy remains trapped.

Thomas J. Maccarone (Texas Tech University), Nathalie Degenaar (Amsterdam), Bailey E. Tetarenko (Texas Tech), Craig Heinke (Alberta), Rudy Wijnands (Amsterdam), Gregory R. Sivakoff (Alberta)

7 pages, accepted to MNRAS

The presence of some X-ray sources in the Galactic Centre region which show variability, but do not show outbursts in over a decade of monitoring has been used to argue for the presence of a large population of stellar mass black holes in this region. A core element of the arguments that these objects are accreting black holes is the claim that neutron stars (NSs) in low mass X-ray binaries (LMXBs) do not have long transient recurrence times. We demonstrate in this paper that about half of the known transient LMXBs with clear signatures for NS primaries have recurrence times in excess of a decade for outbursts at the sensitivity of MAXI. We furthermore show that, in order to reconcile the expected total population of NS LMXBs with the observed one and with the millisecond radio pulsar (MSRP) population of the Galaxy, systems with recurrence times well in excess of a century for outbursts detectable by instruments like MAXI must be the dominant population of NS LMXBs, and that few of these systems have yet to be discovered.

Christian Fruck, Markus Gaug, Alexander Hahn, Victor Acciari, Jürgen Besenrieder, Dijana Dominis Prester, Daniela Dorner, David Fink, Lluís Font, Saša Mićanović, Razmik Mirzoyan, Dominik Müller, Lovro Pavletić, Felix Schmuckermaier, Martin Will

40 pages, submitted to MNRAS

We present a new elastic LIDAR concept, based on a bi-axially mounted Nd:YAG laser and a telescope with HPD readout, combined with fast FADC signal digitization and offline pulse analysis. The LIDAR return signals have been extensively quality checked and absolutely calibrated. We analyze seven years of quasi-continuous LIDAR data taken during those nights when the MAGIC telescopes were operating. Characterization of the nocturnal ground layer yields zenith and azimuth angle dependent aerosol extinction scale heights for clear nights. We derive aerosol transmission statistics for light emitted from various altitudes throughout the year and separated by seasons. We find further seasonal dependencies of cloud base and top altitudes, but none for the LIDAR ratios of clouds. Finally, the night sky background light is characterized using the LIDAR photon backgrounds. abstract.txt

Jais Kumar (1), Prasun Dutta (1), Samir Choudhuri (2), Nirupam Roy (3) ((1) Department of Physics, Indian Institute of Technology (Banaras Hindu University), Varanasi, India, (2) Astronomy Unit, Queen Mary University of London, London, United Kingdom, (3) Department of Physics, Indian Institute of Science, Bangalore, India)

15 Pages, 9 Figures, Accepted for Publication in MNRAS

Observation of redshifted 21-cm signals from neutral hydrogen holds the key to understanding the structure formation and its evolution during the reionization and post-reionization era. Apart from the presence of orders of magnitude larger foregrounds in the observed frequency range, the instrumental effects of the interferometers combined with the ionospheric effects present a considerable challenge in the extraction of 21-cm signals from strong foregrounds. The systematic effects of time and frequency correlated residual gain errors originating from the measurement process introduce a bias and enhance the variance of the power spectrum measurements. In this work, we study the effect of time-correlated residual gain errors in the presence of strong foreground. We present a method to produce analytic estimates of the bias and vari ance in the power spectrum. We use simulated observations to confirm the efficacy of this method and then use it to understand various effects of the gain errors. We find that as the standard deviation in the residual gain errors increases, the bias in the estimation supersedes the variance. It is observed that an optimal choice of the time over which the gain solutions are estimated minimizes the risk. We also find that the interferometers with higher baseline densities are preferred instruments for these studies.

S. Geier, M. Dorsch, I. Pelisoli, N. Reindl, U. Heber, A. Irrgang

28 pages, A&A accepted

Hot subdwarf stars represent a late and peculiar stage in the evolution of low-mass stars, because they are likely formed by close binary interactions. Here we performed a radial velocity (RV) variability study of a sample of 646 hot subdwarfs with multi-epoch radial velocities from Sloan Digital Sky Survey (SDSS) and Large Sky Area Multi-Object Fibre Spectroscopic Telescope (LAMOST) spectra [...] As diagnostics we used the fraction of RV-variable stars and the distribution of the maximum RV variations $\Delta RV_{\rm max}$. Both indicators turned out to be quite inhomogeneous across the studied parameter ranges. A striking feature is the completely different behaviour of He-poor and He-rich hot subdwarfs. While the former have a high fraction of close binaries, almost no significant RV variations could be detected for the latter. This led us to the conclusion that there likely is no evolutionary connection between these subtypes. Intermediate He-rich- and extreme He-rich sdOB/Os on the other hand are likely related. We conclude further that the vast majority of this population is formed via one or several binary merger channels. Hot subdwarfs with temperatures cooler than $\sim24\,000\,{\rm K}$ tend to show less and smaller RV-variations. These objects might constitute a new subpopulation of binaries with longer periods and late-type or compact companions. The RV-variability properties of the extreme horizontal branch (EHB) and corresponding post-EHB populations of the He-poor hot subdwarfs match and confirm the predicted evolutionary connection between them. Stars found below the canonical EHB at somewhat higher surface gravities show large RV-variations and a high RV-variability fraction, which is consistent with most of them being low-mass EHB stars or progenitors of low-mass helium white dwarfs in close binaries.

Qing-Wen Deng, Feng Wang, Hui Deng, Ying Mei, Jing Li, Oleg Smirnov, Shao-Guang Guo

11 pages, 10 figures, accepted by RAA

The Square Kilometre Array (SKA) is the largest radio interferometer under construction in the world. Its immense amount of visibility data poses a considerable challenge to the subsequent processing by the science data processor (SDP). Baseline dependent averaging (BDA), which reduces the amount of visibility data based on the baseline distribution of the radio interferometer, has become a focus of SKA SDP development. This paper developed and implemented a full-featured BDA module based on Radio Astronomy Simulation, Calibration and Imaging Library (RASCIL). Simulated observations were then performed with RASCIL based on a full-scale SKA1-LOW configuration. The performance of the BDA was systematically investigated and evaluated based on the simulated data. The experimental results presented that the amount of visibility data is reduced by about 50\% to 85\% for different time intervals ($\Delta t_{max}$). In addition, different $\Delta t_{max}$ have a significant effect on the imaging quality. The smaller the $\Delta t_{max}$, the smaller the degradation of the imaging quality.

Liang Zhang, Ruisheng Zheng, Changhui Rao, Bing Wang, Huadong Chen, Libo Zhong, Yao Chen

It is generally believed that filament formation involves a process of the accumulation of magnetic energy. However, in this paper we discuss the idea that filaments will not erupt and will only deform when the stored magnetic energy is released gradually. Combining high-quality observations from Solar Dynamics Observatory and other instruments, we present the formation and immediate deformation of a small filament (F1) in the active region (AR) 12760 on 28-30 April 2020. Before the filament formation, three successive dipoles quickly emerged with separation motions in the center of AR 12760. Due to the magnetic interaction between magnetic dipoles and pre-existing positive polarities, coronal brightenings consequently appeared in the overlying atmosphere. Subsequently, because of the continuous cancellation of magnetic flux that happened around the adjacent ends of F1 and another nearby filament (F2), the magnetic reconections occurred intermittently occurred between F1 and F2. Finally, F1 lessened in the shear, and F2 became shorter. All the results show that the formation of F1 was closely associated with intermittent interactions between the sequence of emerging dipoles and pre-existing magnetic polarities, and the immediate deformation of F1 was intimately related to intermittent interactions between F1 and F2. We also suggest that the intermittent magnetic interactions driven by the continuous magnetic activities (magnetic-flux emergence, cancellation, and convergence) play an important role in the formation and deformation of filaments.

Niels F.W. Ligterink, Aida Ahmadi, Bijaya. Luitel, Audrey Coutens, Hannah Calcutt, Łukasz Tychoniec, Harold Linnartz, Jes K. Jørgensen, Robin T. Garrod, Jordy Bouwman

Accepted for publication in ACS Earth and Space Chemistry; 59 pages, 10 figures

This work aims to constrain the abundances of interstellar amides, by searching for this group of prebiotic molecules in the intermediate-mass protostar Serpens SMM1-a. ALMA observations are conducted toward Serpens SMM1. A spectrum is extracted toward the SMM1-a position and analyzed with the CASSIS line analysis software for the presence of characteristic rotational lines of a number of amides and other molecules. NH$_{2}$CHO, NH$_{2}$CHO $\nu_{12}$=1, NH$_{2}^{13}$CHO, CH$_{3}$C(O)NH$_{2}$ $\nu$=0,1, CH$_{2}$DOH, CH$_{3}$CHO, and CH$_{3}$C(O)CH$_{3}$ are securely detected, while trans-NHDCHO, NH$_{2}$CDO, CH$_{3}$NHCHO $\nu$=0,1, CH$_{3}$COOH, and HOCH$_{2}$CHO are tentatively identified. The results of this work are compared with detections presented in the literature. A uniform CH$_{3}$C(O)NH$_{2}$/NH$_{2}$CHO ratio is found for a group of interstellar sources with vast physical differences. A similar ratio is seen for CH$_{3}$NHCHO, based on a smaller data sample. The D/H ratio of NH$_{2}$CHO is about 1--3\% and is close to values found in the low-mass source IRAS~16293--2422B. The formation of CH$_{3}$C(O)NH$_{2}$ and NH$_{2}$CHO is likely linked. Formation of these molecules on grain surfaces during the dark cloud stage is a likely scenario. The high D/H ratio of NH$_{2}$CHO is also seen as an indication that these molecules are formed on icy dust grains. As a direct consequence, amides are expected to be present in the most pristine material from which planetary systems form, thus providing a reservoir of prebiotic material.

Sándor Frey, Ingrid Tar, Krisztina Perger

5 pages, 2 figures. To appear in the proceedings of the European VLBI Network Mini-Symposium and Users' Meeting 2021, Proceedings of Science, PoS(EVN2021)006

RX J1456.0+5048 is a prominent X-ray source detected by ROSAT. There is ~100-mJy level radio emission associated with the X-ray source. However, interferometric observations with increasing angular resolution revealed that three distinct objects located within 2 arcmin are responsible for the measured total flux density. Whether these radio sources lining up in the sky are physically associated or just seen close to each other in projection is not immediately clear. In fact, incorrect cross-identification of the X-ray, optical and radio sources can already be found in the literature. Here we summarise the current knowledge about this intriguing group of objects, where two of the three sources show compact radio emission detected with very long baseline interferometry (VLBI). We present a VLBI image of one of them for the first time, based on archival European VLBI Network (EVN) data taken at 5 GHz.

Tomohisa Kawashima, Seiji Ishiguro, Toseo Moritaka, Ritoku Horiuchi, Kohji Tomisaka

13 pages, 11 figures, accepted for publication in ApJ

We study the kinetic plasma dynamics in collisionless relativistic jets with velocity shear, by carrying out particle-in-cell simulations in the transverse plane of a jet. It is discovered that intermittent magnetic reconnections (MRs) are driven by Mushroom instability (MI), which is an important kinetic-scale plasma instability in the plasma shear-flows with relativistic bulk speed. We refer to this sequence of kinetic plasma phenomena as "MI-driven MR". The MI-driven MRs intermittently occur with moving the location of the reconnection points from the vicinity of the initial velocity-shear surface towards the center of the jet. As a consequence, the number density of high energy electrons, which are accelerated by MI-driven MRs, increases with time in the region inside the initial velocity-shear surface with accompanying the generation and subsequent amplification of magnetic fields by MI. The maximum Lorentz factor of electrons increases with initial bulk Lorentz factor of the jet. A possible relation of MI-driven MR to the bright synchrotron emission in jet-spine of active galactic nucleus jets is also discussed.

In the next decades, it is necessary to forge new late-universe cosmological probes to precisely constrain the Hubble constant and the equation of state of dark energy simultaneously. In this work, we show that the four typical late-universe cosmological probes, the 21 cm intensity mapping (IM), fast radio burst (FRB), gravitational wave (GW) standard siren, and strong gravitational lensing (SGL), are expected to be forged into useful tools in solving the Hubble tension and exploring dark energy. We propose that the synergy of them is rather important in cosmology. We simulate the 21 cm IM, FRB, GW, and SGL data based on the hypothetical observations of the Hydrogen Intensity and Real-time Analysis eXperiment (HIRAX), the Square Kilometre Array (SKA), the Einstein Telescope (ET), and the Large Synoptic Survey Telescope (LSST), respectively. We find that the four probes show obviously different parameter degeneracy orientations in cosmological constraints, so any combination of them can break the parameter degeneracies and thus significantly improve the constraint precision. The joint 21 cm IM+FRB+GW+SGL data can provide the constraint errors of $\sigma(\Omega_{\rm m})=0.0022$ and $\sigma(H_0)=0.16\ \rm km\ s^{-1}\ Mpc^{-1}$ in the $\Lambda$CDM model, which meet the standard of precision cosmology, i.e., the constraint precision of parameters is better than 1%. In addition, the joint data give $\sigma(w)=0.020$ in the $w$CDM model, and $\sigma(w_0)=0.066$ and $\sigma(w_a)=0.25$ in the $w_0w_a$CDM model, which are all better than the constraints obtained by the CMB+BAO+SN data. We show that the synergy between the four late-universe cosmological probes has magnificent prospects.

K.H. Lee, I. Bartos, A. Cook, A. Corsi, Z. Marka, G.C. Privon, S. Marka

6 pages, 1 figure

The heaviest elements in the Universe are synthesized through rapid neutron capture ($r$-process) in extremely neutron rich outflows. Neutron star mergers were established as an important $r$-process source through the multi-messenger observation of GW170817. Collapsars were also proposed as a potentially major source of heavy elements; however, this is difficult to probe through optical observations due to contamination by other emission mechanisms. Here we present observational constraints on $r$-process nucleosynthesis by collapsars based on radio follow-up observations of nearby long gamma-ray bursts. We make the hypothesis that late-time radio emission arises from the collapsar wind ejecta responsible for forging $r$-process elements, and consider the constraints that can be set on this scenario using radio observations of a sample of Swift/BAT GRBs located within 2 Gpc. No radio counterpart was identified in excess of the radio afterglow of the GRBs in our sample, limiting the collapsar $r$-process contribution to $\lesssim0.2$ M$_\odot$ under the models we considered, with constant circum-merger densities giving more stringent constraints. While our results are in tension with collapsars being the majority $r$-process production sites, the ejecta mass and velocity profile of collapsar winds is not yet well modeled. As such, our results are currently subject to large uncertainties, but further theoretical work could greatly improve them.

V. Aivazyan, M. Almualla, S. Antier, A. Baransky, K. Barynova, S. Basa, F. Bayard, S. Beradze, D. Berezin, M. Blazek, D. Boutigny, D. Boust, E. Broens, O. Burkhonov, A. Cailleau, N. Christensen, D. Cejudo, A. Coleiro, M. W. Coughlin, D. Datashvili, T. Dietrich, F. Dolon, J.-G. Ducoin, P.-A. Duverne, G. Marchal-Duval, C. Galdies, L. Granier, V. Godunova, P. Gokuldass, H. B. Eggenstein, M. Freeberg, P. Hello, R. Inasaridze, E. O. Ishida, P. Jaquiery, D. A. Kann, G. Kapanadze, S. Karpov, R. W. Kiendrebeogo, A. Klotz, R. Kneip, N. Kochiashvili, W. Kou, F. Kugel, C. Lachaud, S. Leonini, A. Leroy, N. Leroy, A. Le Van Su, D. Marchais, M. Masek, T. Midavaine, A. Moller, D. Morris, R. Natsvlishvili, F. Navarete, K. Noysena, S. Nissanke, K. Noonan, N. B. Orange, J. Peloton, A. Popowicz, T. Pradier, et al. (20 additional authors not shown)

22 pages, 9 figures

We present our follow-up observations with GRANDMA of transient sources revealed by the Zwicky Transient Facility (ZTF). Over a period of six months, all ZTF triggers were examined in real time by a dedicated science module implemented in the Fink broker, which will be used for the data processing of the Vera C. Rubin Observatory. In this article, we present three selection methods to identify kilonova candidates. Out of more than 35 million candidates, a hundred sources have passed our selection criteria. Six were then followed-up by GRANDMA (by both professional and amateur astronomers). The majority were finally classified either as asteroids or as supernovae events. We mobilized 37 telescopes, bringing together a large sample of images, taken under various conditions and quality. To complement the orphan kilonova candidates (those without associated gamma-ray bursts, which were all), we included three additional supernovae alerts to conduct further observations of during summer 2021. We demonstrate the importance of the amateur astronomer community that contributed images for scientific analyzes of new sources discovered in a magnitude range r'=17-19 mag. We based our rapid kilonova classification on the decay rate of the optical source that should exceed 0.3 mag/day. GRANDMA's follow-up determined the fading rate within 1.5+/-1.2 days post-discovery, without waiting for further observations from ZTF. No confirmed kilonovae were discovered during our observing campaign. This work will be continued in the coming months in the view of preparing for kilonova searches in the next gravitational-wave observing run O4.

Ondrej Jaura, Simon C. O. Glover, Katharina M. J. Wollenberg, Ralf S. Klessen, Sam Geen, Lionel Haemmerlé

22 pages, 20 figures, 4 tables. Accepted by MNRAS

Radiative feedback from massive Population III (Pop III) stars in the form of ionising and photodissociating photons is widely believed to play a central role in shutting off accretion onto these stars. Understanding whether and how this occurs is vital for predicting the final masses reached by these stars and the form of the Pop III stellar initial mass function. To help us better understand the impact of UV radiation from massive Pop III stars on the gas surrounding them, we carry out high resolution simulations of the formation and early evolution of these stars, using the AREPO moving-mesh code coupled with the innovative radiative transfer module SPRAI. Contrary to most previous results, we find that the ionising radiation from these stars is trapped in the dense accretion disk surrounding them. Consequently, the inclusion of radiative feedback has no significant impact on either the number or the total mass of protostars formed during the 20 kyr period that we simulate. We show that the reason that we obtain qualitatively different results from previous studies of Pop III stellar feedback lies in how the radiation is injected into the simulation. HII region trapping only occurs if the photons are injected on scales smaller than the local scale height of the accretion disk, a criterion not fulfilled in previous 3D simulations of this process. Finally, we speculate as to whether outflows driven by the magnetic field or by Lyman-alpha radiation pressure may be able to clear enough gas away from the star to allow the HII region to escape from the disk.

Arka Banerjee, Subinoy Das, Anshuman Maharana, Ravi Kumar Sharma

6 figures

Cosmologies with Light Massive Relics (LiMRs) as a subdominant component of the dark sector are well-motivated from a particle physics perspective, and can also have implications for the $\sigma_8$ tension between early and late time probes of matter clustering. The effects of LiMRs on the Cosmic Microwave Background (CMB) and structure formation on large (linear) scales have been investigated extensively. In this paper, we initiate a systematic study of the effects of LiMRs on smaller, nonlinear scales using cosmological $N$-body simulations; focusing on quantities relevant for photometric galaxy surveys. For most of our study, we use a particular model of nonthermal LiMRs but the methods developed easily generalize to a large class of models of LiMRs -- we explicitly demonstrate this by considering the Dodelson-Widrow form of the velocity distribution. We find that, in general, the effects of LiMR on small scales are distinct from those of a $\Lambda$CDM universe, even when the value of $\sigma_8$ is matched between the models. We show that weak lensing measurements around massive clusters, between $\sim 0.1 h^{-1}$Mpc and $\sim 10 h^{-1}$Mpc, should have sufficient signal-to-noise in future surveys to distinguish between $\Lambda$CDM and LiMR models that are tuned to fit both CMB data and large (linear) scale structure data at late times. Furthermore, we find that different LiMR cosmologies which are indistinguishable by conventional linear probes can be distinguished by these probes if their velocity distributions are sufficiently different. LiMR models can, therefore, be best tested and constrained by jointly analyzing data from CMB and late-time structure formation on both large \textit{and} small scales.

Arya Farahi, Dhayaa Anbajagane, August Evrard

The code is publicly available at this https URL and can be installed through `pip install kllr`

The underlying physics of astronomical systems governs the relation between their measurable properties. Consequently, quantifying the statistical relationships between system-level observable properties of a population offers insights into the astrophysical drivers of that class of systems. While purely linear models capture behavior over a limited range of system scale, the fact that astrophysics is ultimately scale-dependent implies the need for a more flexible approach to describing population statistics over a wide dynamic range. For such applications, we introduce and implement a class of Kernel-Localized Linear Regression (KLLR) models. KLLR is a natural extension to the commonly-used linear models that allows the parameters of the linear model -- normalization, slope, and covariance matrix -- to be scale-dependent. KLLR performs inference in two steps: (1) it estimates the mean relation between a set of independent variables and a dependent variable and; (2) it estimates the conditional covariance of the dependent variables given a set of independent variables. We demonstrate the model's performance in a simulated setting and showcase an application of the proposed model in analyzing the baryonic content of dark matter halos. As a part of this work, we publicly release a Python implementation of the KLLR method.

M. J. Boschini, S. Della Torre, M. Gervasi, D. Grandi, G. Johannesson, G. La Vacca, N. Masi, I. V. Moskalenko, S. Pensotti, T. A. Porter, L. Quadrani, P. G. Rancoita, D. Rozza, M. Tacconi

17 pages, 12 figures, 10 tables; Submitted to ApJ

Since its launch, the Alpha Magnetic Spectrometer-02 (AMS-02) has delivered outstanding quality measurements of the spectra of cosmic-ray (CR) species, $\bar{p}$, $e^{\pm}$, and nuclei (H-Si, Fe), which resulted in a number of breakthroughs. The most recent AMS-02 result is the measurement of the spectra of CR sodium and aluminum up to $\sim$2 TV. Given their low solar system abundances, a significant fraction of each element is produced in fragmentations of heavier species, predominantly Ne, Mg, and Si. In this paper, we use precise measurements of the sodium and aluminum spectra by AMS-02 together with ACE-CRIS and Voyager 1 data to test their origin. We show that the sodium spectrum agrees well with the predictions made with the GalProp-HelMod framework, while aluminum spectrum shows a significant excess in the rigidity range from 2-7 GV. In this context, we discuss the origin of other low-energy excesses in Li, F, and Fe found earlier. The observed excesses in Li, F, and Al appear to be consistent with the local Wolf-Rayet (WR) stars hypothesis, invoked to reproduce anomalous $^{22}$Ne/$^{20}$Ne, $^{12}$C/$^{16}$O, and $^{58}$Fe/$^{56}$Fe ratios in CRs, while excess in Fe is likely connected with a past SN activity in the solar neighborhood. We also provide updated local interstellar spectra (LIS) of sodium and aluminum in the rigidity range from few MV to $\sim$2 TV. Our calculations employ the self-consistent GalProp-HelMod framework that has proved to be a reliable tool in deriving the LIS of CR $\bar{p}$, $e^{-}$, and nuclei $Z\le28$.

Knox S. Long, William P. Blair, P. Frank Winkler, Lorenza Della Bruna, Angela Adamo, Anna F. McLeod, Phillippe Amram

28 page, including 14 figures and 4 tables; accepted for publication in ApJ

Here we describe a new study of the SNRs and SNR candidates in nearby face-on spiral galaxy M83, based primarily on MUSE integral field spectroscopy. Our revised catalog of SNR candidates in M83 has 366 objects, 81 of which are reported here for the first time. Of these, 229 lie within the MUSE observation region, 160 of which have spectra with [SII]:Halpha ratios exceeding 0.4, the value generally accepted as confirmation that an emission nebula is shock-heated. Combined with 51 SNR candidates outside the MUSE region with high [SII]:Halpha ratios, there are 211 spectroscopically-confirmed SNRs in M83, the largest number of confirmed SNRs in any external galaxy. MUSE's combination of relatively high spectral resolution and broad wavelength coverage has allowed us to explore two other properties of SNRs that could serve as the basis of future SNR searches. Specifically, most of the objects identified as SNRs on the basis of [SII]:Halpha ratios exhibit more velocity broadening and lower ratios of [SIII]:[SII] emission than HII regions. A search for nebulae with the very broad emission lines expected from young, rapidly expanding remnants revealed none, except for the previously identified B12-174a. The SNRs identified in M83 are, with few exceptions, middle-aged ISM-dominated ones. Smaller diameter candidates show a larger range of velocity broadening and a larger range of gas densities than the larger diameter objects, as expected if the SNRs expanding into denser gas brighten and then fade from view at smaller diameters than those expanding into a more tenuous ISM

Ken-ichi Tadaki, Akiyoshi Tsujita, Yoichi Tamura, Kotaro Kohno, Bunyo Hatsukade, Daisuke Iono, Minju M. Lee, Yuichi Matsuda, Tomonari Michiyama, Tohru Nagao, Kouichiro Nakanishi, Yuri Nishimura, Toshiki Saito, Hideki Umehata, Jorge Zavala

8 pages, 5 figures, accepted for publication in PASJ

We present observations of [NII] 205 $\mu$m, [OIII] 88 $\mu$m and dust emission in a strongly-lensed, submillimeter galaxy (SMG) at $z=6.0$, G09.83808, with the Atacama Large Millimeter/submillimeter Array (ALMA). Both [NII] and [OIII] line emissions are detected at $>12\sigma$ in the 0.8$"$-resolution maps. Lens modeling indicates that the spatial distribution of the dust continuum emission is well characterized by a compact disk with an effective radius of 0.64$\pm$0.02 kpc and a high infrared surface brightness of $\Sigma_\mathrm{IR}=(1.8\pm0.3)\times10^{12}~L_\odot$ kpc$^{-2}$. This result supports that G09.83808 is the progenitors of compact quiescent galaxies at $z\sim4$, where the majority of its stars are expected to be formed through a strong and short burst of star formation. G09.83808 and other lensed SMGs show a decreasing trend of the [NII] line to infrared luminosity ratio with increasing continuum flux density ratio between 63 $\mu$m and 158 $\mu$m, as seen in local luminous infrared galaxies (LIRGs). The decreasing trend can be reproduced by photoionization models with increasing ionization parameters. Furthermore, by combining the [NII]/[OIII] luminosity ratio with far-infrared continuum flux density ratio in G09.83808, we infer that the gas phase metallicity is already $Z\approx 0.5-0.7~Z_\odot$. G09.83808 is likely one of the earliest galaxies that has been chemically enriched at the end of reionization.

Ilaria Ermolli, Fabrizio Giorgi, Mariarita Murabito, Marco Stangalini, Vincenzo Guido, Marco Molinaro, Paolo Romano, Salvatore L. Guglielmino, Giorgio Viavattene, Gianna Cauzzi, Serena Criscuoli, Kevin P. Reardon, Alexandra Tritschler

22 pages, 13 figures, accepted for publication in Astronomy and Astrophysics

The IBIS data Archive (IBIS-A) stores data acquired with the Interferometric BIdimensional Spectropolarimeter (IBIS), which was operated at the Dunn Solar Telescope of the US National Solar Observatory from 2003 to 2019. The instrument provided series of high-resolution narrowband spectropolarimetric imaging observations of the photosphere and chromosphere in the range 5800$-$8600 \AA~ and co-temporal broadband observations in the same spectral range and with the same field of view of the polarimetric data. We present the data currently stored in IBIS-A, as well as the interface utilized to explore such data and facilitate its scientific exploitation. To this purpose we also describe the use of IBIS-A data in recent and undergoing studies relevant to solar physics and space weather research. IBIS-A includes raw and calibrated observations, as well as science-ready data. The latter comprise maps of the circular, linear, and net circular polarization, and of the magnetic and velocity fields derived for a significant fraction of the series available in the archive. IBIS-A furthermore contains links to observations complementary to the IBIS data, such as co-temporal high-resolution observations of the solar atmosphere available from the instruments onboard the Hinode and IRIS satellites, and full-disc multiband images from INAF solar telescopes. IBIS-A currently consists of 30 TB of data taken with IBIS during 28 observing campaigns performed in 2008 and from 2012 to 2019 on 159 days. Metadata and movies of each calibrated and science-ready series are also available to help users evaluating observing conditions. IBIS-A represents a unique resource for investigating the plasma processes in the solar atmosphere and the solar origin of space weather events.

There is much debate about the channels for astrophysical origins of the stellar-mass binary black hole (BBH) mergers detected by LIGO and Virgo. Active galactic nuclei (AGNs) are promising sites for the efficient formation and rapid mergers of BBHs due to migration traps in high-density gas disks within an inner radii. In this paper, we carry out Monte Carlo simulations to explore the mass properties of mergers over time and hierarchical mergers|with one of the black holes (BHs) being the remnant of a previous merger. We find that the predicted merger rate is $\sim 27-37~{\rm Gpc^{-3}~yr^{-1}}$ and the detection rate of LIGO/Virgo accompanying with extreme mass-ratio mergers will increase in the late stage of AGNs. The fraction of hierarchical mergers is $\sim 24\%$, and its mass-ratio peak is $\sim 0.15-0.35$. Compared with the low-generation mergers in hierarchical mergers, the mass ratio of high-generation mergers has a flatter distribution around its peak. These reveal that the BBH mergers detected by LIGO/Virgo with the extreme mass ratio or heavy component masses can be well explained in AGN channel.

Yoshifusa Ita (1), Takashi Ichikawa (1), Hironori Tsutsui (2), Takumi Hanaue (1), Takahiro Komiyama (1), Hiroki Onozato (3, 1), Atsushi Iwamatsu (1), Ryosuke Morita (1), Yuta Habasaki (1), Ryuto Amemiya (1), Miho Hanawa (1), Kenshi Yanagisawa (3, 2), Hideyuki Izumiura (2), Yoshikazu Nakada (4) ((1) Astronomical Institute, Graduate School of Science, Tohoku University, 6-3 Aramaki Aoba, Aoba-ku, Sendai, Miyagi 980-8578, Japan, (2) Subaru Telescope Okayama Branch Office, National Astronomical Observatory of Japan, 3037-5 Honjo, Kamogata, Asakuchi, Okayama, 719-0232, Japan, (3) National Astronomical Observatory of Japan, 2-21-1 Osawa, Mitaka, Tokyo 181-8588, Japan, (4) Institute of Astronomy, Graduate School of Science, The University of Tokyo, 2-21-1 Osawa, Mitaka, Tokyo, 181-0015, Japan)

Accepted for publication in Publications of the Astronomical Society of Japan

A near-infrared telescope with an effective aperture diameter of thirty millimeters has been developed. The primary objective of the development is to observe northern bright stars in the $J$, $H$, and $K_{\rm s}$ bands and provide accurate photometric data on those stars. The second objective is to repeatedly observe a belt-like region along the northern Galactic plane ($|b| \le 5^\circ$ and $\delta \ge -30^\circ$) to monitor bright variable stars there. The telescope has been in use since December 2016. The purpose of this paper is to describe the design and operational performances of the telescope, photometric calibration methods, and our scientific goals. We show that the telescope has the ability to provide photometry with an uncertainty of less than 5\% for stars brighter than 7, 6.5, and 6~mag in the $J$, $H$, and $K_{\rm s}$ bands, respectively. The repeatability of the photometric measurements for the same star is better than 1\% for bright stars. Our observations will provide accurate photometry on bright stars that are lacking in the Two Micron Sky Survey and the Two Micron All-Sky Survey. Repeated observations at a good cadence will also reveal their nature of the variability in the near-infrared.

Q. Lin, D. Fouchez, J. Pasquet, M. Treyer, R. Ait Ouahmed, S. Arnouts, O. Ilbert

29 pages, 12+11 figures, 2+3 tables; accepted in Astronomy & Astrophysics

Deep Learning models have been increasingly exploited in astrophysical studies, yet such data-driven algorithms are prone to producing biased outputs detrimental for subsequent analyses. In this work, we investigate two major forms of biases, i.e., class-dependent residuals and mode collapse, in a case study of estimating photometric redshifts as a classification problem using Convolutional Neural Networks (CNNs) and galaxy images with spectroscopic redshifts. We focus on point estimates and propose a set of consecutive steps for resolving the two biases based on CNN models, involving representation learning with multi-channel outputs, balancing the training data and leveraging soft labels. The residuals can be viewed as a function of spectroscopic redshifts or photometric redshifts, and the biases with respect to these two definitions are incompatible and should be treated in a split way. We suggest that resolving biases in the spectroscopic space is a prerequisite for resolving biases in the photometric space. Experiments show that our methods possess a better capability in controlling biases compared to benchmark methods, and exhibit robustness under varying implementing and training conditions provided with high-quality data. Our methods have promises for future cosmological surveys that require a good constraint of biases, and may be applied to regression problems and other studies that make use of data-driven models. Nonetheless, the bias-variance trade-off and the demand on sufficient statistics suggest the need for developing better methodologies and optimizing data usage strategies.

Ting Li, Xudong Sun, Yijun Hou, Anqin Chen, Shuhong Yang, Jun Zhang

7 pages, 4 figures, 1 table and 1 online database

With the aim of investigating how the magnetic field in solar active regions (ARs) controls flare activity, i.e., whether a confined or eruptive flare occurs, we analyze 106 flares of Geostationary Operational Environmental Satellite (GOES) class $\geq$M1.0 during 2010$-$2019. We calculate mean characteristic twist parameters $\alpha$$_{FPIL}$ within the "flaring polarity inversion line" region and $\alpha$$_\mathrm{HFED}$ within the area of high photospheric magnetic free energy density, which both provide measures of the nonpotentiality of AR core region. Magnetic twist is thought to be related to the driving force of electric current-driven instabilities, such as the helical kink instability. We also calculate total unsigned magnetic flux ($\Phi$$_\mathrm{AR}$) of ARs producing the flare, which describes the strength of the background field confinement. By considering both the constraining effect of background magnetic fields and the magnetic non-potentiality of ARs, we propose a new parameter $\alpha$/$\Phi$$_\mathrm{AR}$ to measure the probability for a large flare to be associated with a coronal mass ejection (CME). We find that in about 90\% of eruptive flares, $\alpha$$_\mathrm{FPIL}$/$\Phi$$_\mathrm{AR}$ and $\alpha$$_\mathrm{HFED}$/$\Phi$$_\mathrm{AR}$ are beyond critical values (2.2$\times$$10^{-24}$ and 3.2$\times$$10^{-24}$ Mm$^{-1}$ Mx$^{-1}$), whereas they are less than critical values in $\sim$ 80\% of confined flares. This indicates that the new parameter $\alpha$/$\Phi$$_\mathrm{AR}$ is well able to distinguish eruptive flares from confined flares. Our investigation suggests that the relative measure of magnetic nonpotentiality within the AR core over the restriction of the background field largely controls the capability of ARs to produce eruptive flares.

This work presents a first attempt to apply fuzzy cluster analysis (FCA) to analyzing stellar spectra. FCA is adopted to categorize line indices measured from LAMOST low-resolution spectra, and automatically remove the least metallicity-sensitive indices. The FCA-processed indices are then transferred to the artificial neural network (ANN) to derive metallicities for 147 very metal-poor (VMP) stars that have been analyzed by high-resolution spectroscopy. The FCA-ANN method could derive robust metallicities for VMP stars, with a precision of 0.2 dex compared with high-resolution analysis. The recommended FCA threshold value \lambda for this test is between 0.9965 and 0.9975. After reducing the dimension of the line indices through FCA, the derived metallicities are still robust, with no loss of accuracy, and the FCA-ANN method performs stably for different spectral quality from [Fe/H] from -1.8 down to -3.5. Compared with traditional classification methods, FCA considers ambiguity in groupings and noncontinuity of data, and is thus more suitable for observational data analysis. Though this early test uses FCA to analyze low-resolution spectra, and feeds the input to the ANN method to derive metallicities, FCA should be able to, in the large data era, also analyze slitless spectroscopy and multiband photometry, and prepare the input for methods not limited to ANN, in the field of stellar physics for other studies, e.g., stellar classification, identification of peculiar objects. The literature-collected high-resolution sample can help improve pipelines to derive stellar metallicities, and systematic offsets in metallicities for VMP stars for three published LAMOST catalogs have been discussed.

The next Galactic core-collapse supernova (SN) should yield a large number of observed neutrinos. Using Bayesian techniques, we show that with an SN at a known distance up to 25 kpc, the neutrino events in a water Cherenkov detector similar to Super-Kamiokande (SK) could be used to distinguish between seven one-dimensional neutrino emission models assuming no flavor oscillations or the standard Mikheyev-Smirnov-Wolfenstein effect. Some of these models could still be differentiated with an SN at a known distance of 50 kpc. We also consider just the relative distributions of neutrino energy and arrival time predicted by the models and find that a detector like SK meets the requirement to distinguish between these distributions with an SN at an unknown distance up to $\sim 10$ kpc.

The positive-to-negative transition of spectral lag is an uncommon feature reported in a small number of GRBs. An application of such a feature has been made to constrain the critical quantum gravity energy ($E_{\rm QG}$) of the light photons under the hypothesis that the Lorentz invariance might be violated. Motivated by previous case studies, this paper systematically examined the up-to-date Fermi/GBM GRB sample for the lag transition feature to establish a comprehensive physical limit on the Lorentz Invariance Violation (LIV). This search resulted in 32 GRBs with redshift available, which exhibit the lag-transition phenomenon. We first fit each of the lag-E relations of the 32 GRBs with an empirical smoothly broken power law function, and found that the lag transition occurs typically at about 400 keV. We then implemented the LIV effect into the fit, which enabled us to constrain the lower limit of the linear and quadratic values of $E_{\rm QG}$, which are typically distributed at $1.5\times 10^{14}$ GeV and $8\times 10^{5}$ GeV, respectively.

Tianjun Gan, Abderahmane Soubkiou, Sharon X. Wang, Zouhair Benkhaldoun, Shude Mao, Étienne Artigau, Pascal Fouqué, Steven Giacalone, Christopher A. Theissen, Christian Aganze, Karen A. Collins, Avi Shporer, Khalid Barkaoui, Mourad Ghachoui, Steve B. Howell, Claire Lamman, Olivier D. S. Demangeon, Artem Burdanov, Charles Cadieux, Jamila Chouqar, Kevin I. Collins, Neil J. Cook, Laetitia Delrez, Brice-Olivier Demory, René Doyon, Georgina Dransfield, Courtney D. Dressing, Elsa Ducrot, Jiahao Fan, Lionel Garcia, Holden Gill, Michaël Gillon, Crystal L. Gnilka, Yilen Gómez Maqueo Chew, Maximilian N. Günther, Christopher E. Henze, Chelsea X. Huang, Emmanuel Jehin, Eric L. N. Jensen, Zitao Lin, James McCormac, Catriona A. Murray, Prajwal Niraula, Peter P. Pedersen, Francisco J. Pozuelos, et al. (18 additional authors not shown)

19 pages, 15 figures, submitted to MNRAS

We present the discovery of TOI-2136b, a sub-Neptune planet transiting every 7.85 days a nearby M4.5V-type star, identified through photometric measurements from the TESS mission. The host star is located $33$ pc away with a radius of $R_{\ast} = 0.34\pm0.02\ R_{\odot}$, a mass of $0.34\pm0.02\ M_{\odot}$ and an effective temperature of $\rm 3342\pm100\ K$. We estimate its stellar rotation period to be $75\pm5$ days based on archival long-term photometry. We confirm and characterize the planet based on a series of ground-based multi-wavelength photometry, high-angular-resolution imaging observations, and precise radial velocities from CFHT/SPIRou. Our joint analysis reveals that the planet has a radius of $2.19\pm0.17\ R_{\oplus}$, and a mass measurement of $6.4\pm2.4\ M_{\oplus}$. The mass and radius of TOI2136b is consistent with a broad range of compositions, from water-ice to gas-dominated worlds. TOI-2136b falls close to the radius valley for low-mass stars predicted by the thermally driven atmospheric mass loss models, making it an interesting target for future studies of its interior structure and atmospheric properties.

Hachem Dhouib, Stéphane Mathis, Lisa Bugnet, Timothy Van Reeth, Conny Aerts

21 pages, 15 figures, 3 tables, abstract shortened for arXiv. Accepted for publication in A&A

Asteroseismology has revealed small core-to-the-surface rotation contrasts in stars in the whole HR diagram. This is the signature of strong transport of angular momentum (AM) in stellar interiors. One of the plausible candidates to efficiently carry AM is magnetic fields with various topologies that could be present in stellar radiative zones. Among them, strong axisymmetric azimuthal magnetic fields have received a lot of interest. Indeed, if they are subject to the so-called Tayler instability, the accompanying triggered Maxwell stresses can transport AM efficiently. In addition, the electromotive force induced by the fluctuations of magnetic and velocity fields could potentially sustain a dynamo action that leads to the regeneration of the initial strong axisymmetric azimuthal magnetic field. The key question we aim to answer is: can we detect signatures of these deep strong azimuthal magnetic fields? The only way to answer this question is asteroseismology and the best laboratories of study are intermediate-mass and massive stars. Most of these are rapid rotators during their main-sequence. Therefore, we have to study stellar pulsations propagating in stably stratified, rotating, and potentially strongly magnetised radiative zones. We generalise the traditional approximation of rotation by taking simultaneously general axisymmetric differential rotation and azimuthal magnetic fields into account in a non-perturbative way. Using this new formalism, we derive the asymptotic properties of magneto-gravito-inertial (MGI) waves and their period spacings. We find that toroidal magnetic fields induce a shift in the period spacings of MGI modes. An equatorial azimuthal magnetic field with an amplitude of the order of $10^5\,\rm G$ leads to signatures that can be detectable thanks to modern space photometry. More complex hemispheric configurations are more difficult to observe.

Fan Yang, Wei Wang, Xing Wei, Hui Zhang, Ji-Lin Zhou, Su-Su Shan, Jie Zheng, Wei-Kai Zong, Ming Yang, Yu Bai, Song Wang, Jia-Chen Zheng, Yu-Ru Xu, Yu-Feng Li, You-Jun Lu, Ji-Feng Liu

Accept for publication to RAA

Transit Timing Variation (TTV) of hot Jupiters provides direct observational evidence of planet tidal dissipation. Detecting tidal dissipation through TTV needs high precision transit timings and long timing baselines. In this work, we predict and discuss the potential scientific contribution of SiTian Survey in detecting and analyzing exoplanet TTV. We develop a tidal dissipation detection pipeline for SiTian Survey that aims at time-domain astronomy with 72 1-meter optical telescopes. The pipeline includes the modules of light curve deblending, transit timing obtaining, and TTV modeling. SiTian is capable to detect more than 25,000 exoplanets among which we expect $\sim$50 sources showing evidence of tidal dissipation. We present detection and analysis of tidal dissipating targets, based on simulated SiTian light curves of XO-3b and WASP-161b. The transit light curve modeling gives consistent results within 1$\sigma$ to input values of simulated light curves. Also, the parameter uncertainties predicted by Monte-Carlo Markov Chain are consistent with the distribution obtained from simulating and modeling the light curve 1000 times. The timing precision of SiTian observations is $\sim$ 0.5 minutes with one transit visit. We show that differences between TTV origins, e.g., tidal dissipation, apsidal precession, multiple planets, would be significant, considering the timing precision and baseline. The detection rate of tidal dissipating hot Jupiters would answer a crucial question of whether the planet migrates at an early formation stage or random stages due to perturbations, e.g., planet scattering, secular interaction. SiTian identified targets would be constructive given that the sample would extend tenfold.

Ravit Helled, David J. Stevenson, Jonathan I. Lunine, Scott J. Bolton, Nadine Nettelmann, Sushil Atreya, Tristan Guillot, Burkhard Militzer, Yamila Miguel, William B. Hubbard

invited review for Icarus, accepted for publication

The Juno mission has revolutionized and challenged our understanding of Jupiter. As Juno transitioned to its extended mission, we review the major findings of Jupiter's internal structure relevant to understanding Jupiter's formation and evolution. Results from Juno's investigation of Jupiter's interior structure imply that the planet has compositional gradients and is accordingly non-adiabatic, with a complex internal structure. These new results imply that current models of Jupiter's formation and evolution require a revision. In this paper, we discuss potential formation and evolution paths that can lead to an internal structure model consistent with Juno data, and the constraints they provide. We note that standard core accretion formation models, including the heavy-element enrichment during planetary growth is consistent with an interior that is inhomogeneous with composition gradients in its deep interior. However, such formation models typically predict that this region, which could be interpreted as a primordial dilute core, is confined to about 10% of Jupiter's total mass. In contrast, structure models that fit Juno data imply that this region contains 30% of the mass or more. One way to explain the origin of this extended region is by invoking a relatively long (about 2 Myrs) formation phase where the growing planet accretes gas and planetesimals delaying the runaway gas accretion. Alternatively, Jupiter's fuzzy core could be a result of a giant impact or convection post-formation. These novel scenarios require somewhat special and specific conditions. Clarity on the plausibility of such conditions could come from future high-resolution observations of planet-forming regions around other stars, from the observed and modeled architectures of extrasolar systems with giant planets, and future Juno data obtained during its extended mission.

Ravit Helled, Naor Movshovitz, Nadine Nettelmann

a review chapter, to appear in "Planetary systems now" edited by Luisa M. Lara and David C. Jewitt

Revealing the true nature of the gas giant planets in our Solar System is challenging. The masses of Jupiter and Saturn are about 318 and 95 Earth masses, respectively. While they mostly consist of hydrogen and helium, the total mass and distribution of the heavier elements, which reveal information on their origin, are still unknown. Recent accurate measurements of the gravitational fields of Jupiter and Saturn together with knowledge of the behavior of planetary materials at high pressures allow us to better constrain their interiors. Updated structure models of Jupiter and Saturn suggest that both planets have complex interiors that include composition inhomogeneities, non-convective regions, and fuzzy cores. In addition, it is clear that there are significant differences between Jupiter and Saturn and that each giant planet is unique. This has direct implications for giant exoplanet characterization and for our understanding of gaseous planets as a class of astronomical objects. In this review we summarize the methods used to model giant planet interiors and recent developments in giant planet structure models.

The Hubble tension that the standard $\Lambda$CDM model is suffering from can be resolved with pre-recombination early dark energy. We present the first constraint on the tensor-to-scalar ratio $r$ in corresponding Hubble-tension-free cosmologies using the most recent BICEP/Keck cosmic microwave background (CMB) B-mode polarization data. We find, combining BICEP/Keck with Planck18 CMB and baryon acoustic oscillation data, that the models with larger Hubble constant $H_0$ will have tighter upper bound on $r$, and resolution $H_0\sim73$ km/s/Mpc of the Hubble tension tightens the upper bound to $r<0.028\ (95\%\text{C.L.})$, $25\%$ tighter than the $\Lambda$CDM constraint $r<0.036$. We clarify the origin of this tightening bound.

Anastasia Chaikova, Dmitriy Kostunin, Sergei B. Popov

4 pages, submitted to RNAAS

The CHIME/FRB collaboration has recently published a catalog containing about half a thousand fast radio bursts (FRBs) including their spectra and several reconstructed properties, like signal widths, amplitudes, etc. We have developed a model-independent approach for a classification of these bursts using cross-correlation and clustering algorithms applied to one-dimensional intensity profiles of the bursts (i.e., to amplitudes as a function of time averaged over the frequency). Using this algorithm we identified two major classes of FRBs featuring different waveform morphology, and, simultaneously, different distribution of brightness temperature. Bursts from one of the identified cluster have lower brightness temperatures and larger widths than events from the other cluster. Both groups include bursts from the repeaters and one-off events.

G. Buldgen, M. Farnir, P. Eggenberger, J. Bétrisey, C. Pezzotti, C. Pinçon, M. Deal, S. J. A. J. Salmon

Accepted for publication in Astronomy and Astrophysics

The advent of space-based photometry observations provided high-quality asteroseismic data for a large number of stars. These observations enabled the adaptation of advanced techniques, until then restricted to helioseismology, to study the best asteroseismic targets. Amongst these, the 16Cyg binary system holds a special place, being the brightest solar twins observed by Kepler. For this system, modellers have access to high-quality asteroseismic, spectroscopic and interferometric data, making it the perfect testbed for the limitations of stellar models. We aim to further constrain the internal structure and fundamental parameters of 16CygA&B using linear seismic inversion techniques of both global indicators and localised corrections of the internal structure. We start from the models defined by detailed modelling in our previous paper and extend our analysis by applying variational inversions to these models. We carried out inversions of so-called seismic indicators and provided local corrections of the internal structure of the two stars. Our results indicate that linear seismic inversions alone are not able to discriminate between standard and non-standard models for 16CygA&B. We confirm the results of our previous studies that used linear inversion techniques, but consider that the differences could be linked to small fundamental parameters variations rather than to a missing process in the models. We confirm the robustness and reliability of the results of the modelling performed in our previous paper. We conclude that non-linear inversions are likely required to further investigate the properties of 16CygA&B from a seismic point of view, but that these inversions should be coupled to analyses of the depletion of light elements such as lithium and beryllium to constrain the macroscopic transport of chemicals and potential non-standard evolutionary paths.

G. Buldgen, P. Eggenberger

Proceedings of the Sixteenth Marcel Grossman meeting (July 5 - 10 2021, Session on "Rotation in Stellar Evolution")

The possibility of measuring the internal rotation of the Sun and stars thanks to helio- and asteroseismology offers tremendous constraints on hydro- and magnetohydrodynamical processes acting in stellar interiors. Understanding the processes responsible for the transport of angular momentum in stellar interiors is crucial as they will also influence the transport of chemicals and thus the evolution of stars. Here we present some of the key results obtained in both fields and how detailed seismic analyses can provide stringent constraints on the physics of angular momentum transport in the interior of low mass stars and potentially rule out some candidates.

Andreas J. Weiss, Teresa Nieves-Chinchilla, Christian Möstl, Martin A. Reiss, Tanja Amerstorfer, Rachel L. Bailey

14 pages, 6 figures

We present a new analytical approach with the aim to describe generally curved and twisted magnetic flux ropestructures, that are embedded within interplanetary coronal mass ejections, under the constraint of invariant axialflux. In this paper we showcase the simplest case of a generally curved flux rope with a circular cross-sectionwhich can be described in terms of the curvature and the torsion of the Frenet-Serret equations. The magneticfield configuration, for the axial and poloidal field components, are described in terms of a radial expansion usinga Legendre basis. We further derive equations that allow us to configure our model for any arbitrary magnetictwist and also evaluate the force distribution. We show the effects and differences of our proposed modelcompared to a purely cylindrical or toroidal geometry using an arbitrarily twisted exemplary flux rope structurewith an uniformly twisted magnetic field configuration. In order to indirectly compare our model with realin-situ measurements we generate two synthetic in-situ profiles using virtual spacecraft trajectories, realisticallysimulating apex and flank encounters of an interplanetary coronal mass ejection. This proposed model presentsan intermediate steps towards describing more complex flux rope structures with arbitrary cross-section shapes.

G. Tautvaišienė, Š. Mikolaitis, A. Drazdauskas, E. Stonkutė, R. Minkevičiūtė, E. Pakštienė, H. Kjeldsen, K. Brogaard, Y. Chorniy, C. von Essen, F. Grundahl, M. Ambrosch, V. Bagdonas, A. Sharma, C. Viscasillas Vázquez

17 pages, 8 figures. arXiv admin note: text overlap with arXiv:2005.07526

In fulfilling the aims of the planetary and asteroseismic research missions, such as that of the NASA Transiting Exoplanet Survey Satellite (TESS) space telescope, accurate stellar atmospheric parameters and a detailed chemical composition are required as input. We have observed high-resolution spectra for all 848 bright (V<8 mag) stars that are cooler than F5 spectral class in the area up to 12 deg surrounding the northern TESS continuous viewing zone and uniformly determined the main atmospheric parameters, ages, orbital parameters, velocity components, and precise abundances of up to 24 chemical species (C(C2), N(CN), [O I], Na I, Mg I, Al I, Si I, Si I, Ca I, Ca II, Sc I, Sc II, Ti I, Ti II, V I, Cr I, Cr II, Mn I, Fe I, Fe II, Co I, Ni I, Cu I, and Zn I) for 740 slowly rotating stars. The analysis of 25 planet-hosting stars in our sample drove us to the following conclusions: the dwarf stars hosting high-mass planets are more metal rich than those with low-mass planets. We find slightly negative C/O and Mg/Si slopes toward the stars with high-mass planets. All the low-mass planet hosts in our sample show positive $\Delta$[El/Fe] versus condensation temperature slopes, in particular, the star with the large number of various planets. The high-mass planet hosts have a diversity of slopes, but in more metal rich, older, and cooler stars, the positive elemental abundance slopes are more common.

David Modiano, Rudy Wijnands, Aastha Parikh, Jari van Opijnen, Sill Verberne, Marieke van Etten

31 pages, 39 figures. Accepted for publication in Astronomy & Astrophysics (date of acceptance: 09/02/2022)

Despite the prevalence of transient-searching facilities operating across most wavelengths, the ultraviolet (UV) transient sky remains to be systematically studied. We have recently initiated the Transient Ultraviolet Objects (TUVO) project, with which we search for serendipitous UV transients in data from currently available UV instruments, with a focus on the UV/Optical (UVOT) telescope aboard the Neil Gehrels Swift Observatory (an overview of the TUVO project is described in a companion paper). Here we describe TUVOpipe, the pipeline we constructed in order to find such transients in the UVOT data, using difference image analysis. The pipeline is run daily on all new public UVOT data (which are available 6-8 hours after the observations are performed), so we discover transients in near real-time. This allows for follow-up observations to be performed. From October 1, 2020, to the time of submission, we have processed 111,330 individual UVOT images and we currently detect an average rate of ~100 transient candidates per day. Of these daily candidates, on average ~30 are real transients, separated by human vetting from the remaining `bogus' transients which were not discarded automatically within the pipeline. Most of the real transients we detect are known variable stars, though we also detect many known active galactic nuclei and accreting white dwarfs. TUVOpipe can additionally run in archival mode, whereby all archival UVOT data of a given field is scoured for `historical' transients; in this mode we also mostly find variable stars. However, some of the transients we find (in particular in the real-time mode) represent previously unreported new transients, or undiscovered outbursts of known transients, predominantly outbursts from cataclysmic variables. In this paper we describe the operation of (both modes of) TUVOpipe and some of the initial results we have so far obtained.

Yijun Hou, Ting Li, Shuhong Yang, Shin Toriumi, Yilin Guo, Jun Zhang

23 pages, 11 figures, 2 tables, accepted for publication in ApJ

Light bridges (LBs) are among the most striking sub-structures in sunspots, where various activities have been revealed by recent high-resolution observations from the Interface Region Imaging Spectrograph (IRIS). According to the variety of physical properties, we classified these activities into four distinct categories: transient brightening (TB), intermittent jet (IJ), type-I light wall (LW-I), and type-II light wall (LW-II). In IRIS 1400/1330 {\AA} observations, TBs are characterized by abrupt emission enhancements, and IJs appear as collimated plasma ejections with a width of 1-2 Mm at some LB sites. Most observed TBs are associated with IJs and show superpositions of some chromosphere absorption lines on enhanced and broadened wings of C II and Si IV lines, which could be driven by intermittent magnetic reconnection in the lower atmosphere. LW-I and LW-II are wall-shaped structures with bright fronts above the whole LB. An LW-I has a continuous oscillating front with a typical height of several Mm and an almost stationary period of 4-5 minutes. On the contrary, an LW-II has a indented front with a height of over 10 Mm, which has no stable period and is accompanied by recurrent TBs in the entire LB. These results support that LW-IIs are driven by frequent reconnection occurring along the whole LB due to large-scale magnetic flux emergence or intrusion, rather than the leakage of waves producing LW-Is. Our observations reveal a highly dynamical scenario of activities above LBs driven by different basic physical processes, including magneto-convection, magnetic reconnection, and wave leakage.

K. Kawauchi, F. Murgas, E. Palle, N. Narita, A. Fukui, T. Hirano, H. Parviainen, H. T. Ishikawa, N. Watanabe, E. Esparaza-Borges, M. Kuzuhara, J. Orell-Miquel, V. Krishnamurthy, M. Mori, T. Kagetani, Y. Zou, K. Isogai, J. H. Livingston, S. B. Howell, N. Crouzet, J. P. de Leon, T. Kimura, T. Kodama, J. Korth, S. Kurita, A. Laza-Ramos, R. Luque, A. Madrigal-Aguado, K. Miyakawa, G. Morello, T. Nishiumi, G. E. F. Rodríguez, M. Sánchez-Benavente, M. Stangret, H. Teng, Y. Terada, C. L. Gnilka, N. Guerrero, H. Harakawa, K. Hodapp, Y. Hori, M. Ikoma, S. Jacobson, M. Konishi, T. Kotani, T. Kudo, T. Kurokowa, N. Kusakabe, J. Nishikawa, M. Omiya, T. Serizawa, M. Tamura, A. Ueda, S. Vievard

16 pages, 14 figures, submitted for publication in Astronomy & Astrophysics

The NASA space telescope $TESS$ is currently in the extended mission of its all-sky search for new transiting planets. Of the thousands of candidates that TESS is expected to deliver, transiting planets orbiting nearby M dwarfs are particularly interesting targets since they provide a great opportunity to characterize their atmospheres by transmission spectroscopy. We aim to validate and characterize the new sub-Neptune sized planet candidate TOI-2136.01 orbiting a nearby M dwarf ($d = 33.36 \pm 0.02$ pc, $T_{eff} = 3373 \pm 108$ K) with an orbital period of 7.852 days. We use TESS data, ground-based multi-color photometry, and radial velocity measurements with the InfraRed Doppler (IRD) instrument on the Subaru Telescope to validate the planetary nature of TOI-2136.01 and estimate the stellar and planetary parameters. We also conduct high-resolution transmission spectroscopy to search for helium in its atmosphere. We confirmed that TOI-2136.01 (now named as TOI-2136b) is a bona fide planet with a planetary radius of $R_p = 2.2 \pm 0.07$ $R_{Earth}$ and a mass of $M_p = 4.7^{+3.1}_{-2.6}$ $M_{Earth}$. We also search for helium 10830 \r{A} absorption lines and place an upper limit on the equivalent width of $<$ 7.8 m\r{A} (95 % confidence) and on the absorption signal of $<$ 1.44 % (95 % confidence). TOI-2136b is a sub-Neptune transiting a nearby and bright star (J=10.8) and is a potential hycean planet, making it an excellent target for atmospheric studies to understand the formation, evolution, and habitability of the small planets.

Qiuhan He, James Nightingale, Richard Massey, Andrew Robertson, Aristeidis Amvrosiadis, Shaun Cole, Carlos S. Frenk, Ran Li, Nicola C. Amorisco, R. Benton Metcalf, Xiaoyue Cao, Amy Etherington

Submitted to MNRAS. Comments welcome!

Strong gravitational lensing offers a compelling test of the cold dark matter paradigm, as it allows for subhaloes with masses of $\sim10^{9}$ M$_\odot$ and below to be detected. We test commonly-used techniques for detecting dark matter subhaloes superposed in images of strongly lensed galaxies. For the lens we take a simulated galaxy in a $\sim10^{13}$ M$_\odot$ halo grown in a high-resolution cosmological hydrodynamics simulation, which we view from two different directions. To remove particle noise, we represent the projected galaxy mass distribution by a series of analytic Gaussian profiles which precisely capture the features of the projected galaxy. We first model the lens mass as a broken power-law density profile (the parameterization of which includes a core ensuring the results are unaffected by a numerical central core in the simulated galaxy) and then search for small haloes. Of the two projections, one has a regular elliptical shape, while the other has distinct deviations from an elliptical shape. For the former, the broken power-law model gives no false positives and correctly recovers the mass of the superposed small halo, but for the latter we do find false positives and the inferred halo mass is overestimated by a factor of $\sim4-5$. We then replace the broken power-law model by a more complex model in which the lens mass is decomposed into separate stellar and dark matter components. In this case, we show that we can precisely capture the complex projected structures of the simulated galaxy and correctly infer the input parameters of the superposed small halo. By improving the lens galaxy mass model, we argue this enables more robust subhalo inference that can detect lower mass dark matter subhaloes, strengthening the technique's ability to test alternative models of dark matter which do not form dark matter subhaloes below a cut-off mass.

Andrew G. Fullard, John T. O'Brien, Wolfgang E. Kerzendorf, Manisha Shrestha, Jennifer L. Hoffman, Richard Ignace, Patrick van der Smagt

24 pages, 13 figures, accepted to ApJ

Understanding the evolution of massive binary stars requires accurate estimates of their masses. This understanding is critically important because massive star evolution can potentially lead to gravitational wave sources such as binary black holes or neutron stars. For Wolf-Rayet stars with optically thick stellar winds, their masses can only be determined with accurate inclination angle estimates from binary systems which have spectroscopic $M \sin i$ measurements. Orbitally-phased polarization signals can encode the inclination angle of binary systems, where the Wolf-Rayet winds act as scattering regions. We investigated four Wolf-Rayet + O star binary systems, WR 42, WR 79, WR 127, and WR 153, with publicly available phased polarization data to estimate their masses. To avoid the biases present in analytic models of polarization while retaining computational expediency, we used a Monte Carlo radiative transfer model accurately emulated by a neural network. We used the emulated model to investigate the posterior distribution of parameters of our four systems. Our mass estimates calculated from the estimated inclination angles put strong constraints on existing mass estimates for three of the systems, and disagrees with the existing mass estimates for WR 153. We recommend a concerted effort to obtain polarization observations that can be used to estimate the masses of Wolf-Rayet binary systems and increase our understanding of their evolutionary paths.

A. Bemporad, V. Andretta, R. Susino, S. Mancuso, D. Spadaro, M. Mierla, D. Berghmans, E. D'Huys, A. N. Zhukov, D.-C. Talpeanu, R. Colaninno, P. Hess, J. Koza, S. Jejcic, P. Heinzel, E. Antonucci, V. Da Deppo, S. Fineschi, F. Frassati, G. Jerse, F. Landini, G. Naletto, G. Nicolini, M. Pancrazzi, M. Romoli, C. Sasso, A. Slemer, M. Stangalini, L. Teriaca

15 pages, 14 figures, accepted for publication on A&A

On February 12, 2021 two subsequent eruptions occurred above the West limb, as seen along the Sun-Earth line. The first event was a typical slow Coronal Mass Ejection (CME), followed $\sim 7$ hours later by a smaller and collimated prominence eruption, originating Southward with respect to the CME, followed by a plasma blob. These events were observed not only by SOHO and STEREO-A missions, but also by the suite of remote sensing instruments on-board Solar Orbiter (SolO). This work shows how data acquired by the Full Sun Imager (FSI), Metis coronagraph, and Heliospheric Imager (SoloHI) from the SolO perspective can be combined to study the eruptions and the different source regions. Moreover, we show how Metis data can be analyzed to provide new information about solar eruptions. Different 3D reconstruction methods were applied to the data acquired by different spacecraft including remote sensing instruments on-board SolO. Images acquired by both Metis channels in the Visible Light (VL) and H I Lyman$-\alpha$ line (UV) were combined to derive physical information on the expanding plasma. The polarization ratio technique was also applied for the first time to the Metis images acquired in the VL channel. The two eruptions were followed in 3D from their source region to their expansion in the intermediate corona. Thanks to the combination of VL and UV Metis data, the formation of a post-CME Current Sheet (CS) was followed for the first time in the intermediate corona. The plasma temperature gradient across a post-CME blob propagating along the CS was also measured for the first time. Application of the polarization ratio technique to Metis data shows that, thanks to the combination of four different polarization measurements, the errors are reduced by $\sim 5-7$\%, thus better constraining the 3D distribution of plasma.

Ofek Bengyat, Avishay Gal-Yam

6 pages of text, 7 figures. Submitted to ApJ. Code available in GitHub under this https URL

A quantitative data-driven comparison among supernovae (SNe) based on their spectral time series combined with multi-band photometry is presented. We use an unsupervised Random Forest algorithm as a metric on a set of 82 well-documented SNe representing all the main spectroscopic types, in order to embed these in an abstract metric space reflecting shared correlations between the objects. We visualize the resulting metric space in 3d, revealing strong agreement with the current spectroscopic classification scheme. The embedding splits Type Ib supernovae into two groups, with one subgroup exhibiting broader, less prominent, higher-velocity lines than the other, possibly suggesting a new SN Ib subclass is required. The method could be to classify newly discovered SNe according to their distance from known event groups, or ultimately to devise a new, spectral-temporal classification scheme. Such an embedding could also depend on hidden parameters which may perhaps be physically interpretable.

Angeles I. Diaz, S. Zamora

17 pages, 18 figures, accepted for publication in MNRAS

We present a methodology for the use of sulphur as global metallicity tracer in galaxies, allowing performing a complete abundance analysis using only the red-to-near infrared spectral region. We have applied it to a compilation of high-quality data split into two samples: HII regions (DHR) in spiral and irregular galaxies, and dwarf galaxies dominated by a strong starburst (HIIGal). Sulphur abundances have been derived by direct methods under the assumption of an ionisation structure composed of two zones: an intermediate one where S{++} is originated and a low ionisation one where S{+} is formed. Ionisation correction factors (ICF) have been calculated from the Ar{3+}/Ar{3+} ratio and are shown to correlate with the hardness of the radiation field. Only about 10% of the objects show S{3+} contributions to the total abundance larger than 30%. A good correlation exists between sulphur abundance and ionising temperature with low metallicity objects being ionised by hotter stars. No correlation is found between ionisation parameter and total S/H abundance. Most of the HIIGal objects show S/O ratios below the solar value and a trend for increasing S/O ratios with increasing sulphur abundances while DHR objects show S/O ratios larger than solar and a tendency for lower S/O ratios for higher metallicities. Finally, we present a calibration of the sulphur abundance through the S{23} parameter that remains single valued up to sulphur abundances well beyond the solar value. S{23} is independent of the ionisation parameter and only weakly dependent on ionising temperature.

Joshua Krissansen-Totton, Maggie Thompson, Max L. Galloway, Jonathan J. Fortney

16 pages, 5 figures

The search for life on exoplanets is motivated by the universal ways in which life could modify its planetary environment. Atmospheric gases such as oxygen and methane are promising candidates for such environmental modification due to the evolutionary benefits their production would confer. However, confirming that these gases are produced by life, rather than by geochemical or astrophysical processes, will require a thorough understanding of planetary context, including the expected counterfactual atmospheric evolution for lifeless planets. Here, we evaluate current understanding of planetary context for several candidate biosignatures and their upcoming observability. We review the contextual framework for oxygen and describe how conjectured abiotic oxygen scenarios may be testable. In contrast to oxygen, current understanding of how planetary context controls non-biological methane (CH$_4$) production is limited, even though CH$_4$ biosignatures in anoxic atmospheres may be readily detectable with the James Webb Space Telescope. We assess environmental context for CH$_4$ biosignatures and conclude that abundant atmospheric CH$_4$ coexisting with CO$_2$, and CO:CH$_4$ << 1 is suggestive of biological production, although precise thresholds are dependent on stellar context and sparsely characterized abiotic CH$_4$ scenarios. A planetary context framework is also considered for alternative or agnostic biosignatures. Whatever the distribution of life in the Universe, observations of terrestrial exoplanets in coming decades will provide a quantitative understanding of the atmospheric evolution of lifeless worlds. This knowledge will inform future instrument requirements to either corroborate the presence of life elsewhere or confirm its apparent absence.

Forrest Mozer, Stuart Bale, Cynthia Cattell, Jasper Halekas, Ivan Vasko, Jae Verniero, Paul Kellogg

12 pages, 6 figures. arXiv admin note: substantial text overlap with arXiv:2111.07161

Perihelion passes on Parker Solar Probe orbits six through nine have been studied to show that solar wind core electrons emerged from 15 solar radii with a temperature of 55 plus or minus 5 eV, independent of the solar wind speed which varied from 300 to 800 km/sec. After leaving 15 solar radii and in the absence of triggered ion acoustic waves at greater distances, the core electron temperature varied with radial distance, R, in solar radii, as 1900R-4/3 electron volts because of cooling produced by the adiabatic expansion. The coefficient, 1900, reproduces the minimum core electron perpendicular temperature observed during the 25 days of observation. In the presence of triggered ion acoustic waves, the core electrons were isotropically heated as much as a factor of two above the minimum temperature, 1900R-4/3 eV. Triggered ion acoustic waves were the only waves observed in coincidence with the electron core heating. They are the dominant wave mode at frequencies greater than 100 Hz at solar distances between 15 and 30 solar radii.

S.E. Harper, C. Dickinson, A. Barr, R. Cepeda-Arroita, R.D.P. Grumitt, H.M. Heilgendorff, L. Jew, J.L. Jonas, M.E. Jones, J.P. Leahy, J. Leech, T.J. Pearson, M.W. Peel, A.C.S. Readhead, A.C. Taylor

20 pages, 13 figures, to be submitted to MNRAS

The C-Band All-Sky Survey (C-BASS) has observed the Galaxy at 4.76GHz with an angular resolution of $0.73^\circ$ full-width half-maximum, and detected Galactic synchrotron emission with high signal-to-noise ratio over the entire northern sky ($\delta > -15^{\circ}$). We present the results of a spatial correlation analysis of Galactic foregrounds at mid-to-high ($b > 10^\circ$) Galactic latitudes using a preliminary version of the C-BASS map. We jointly fit for synchrotron, dust, and free-free components between $20$ and $1000$GHz and look for differences in the Galactic synchrotron spectrum, and the emissivity of anomalous microwave emission (AME) when using either the C-BASS map or the 408MHz all-sky map to trace synchrotron emission. We find marginal evidence for a steepening ($\left<\Delta\beta\right> = -0.06\pm0.02$) of the Galactic synchrotron spectrum at high frequencies resulting in a mean spectral index of $\left<\beta\right> = -3.10\pm0.02$ over $4.76-22.8$GHz. Therefore, the synchrotron emission can be well modelled by a single power-law up to tens of GHz. Due to this, we find that the AME emissivity is not sensitive to changing the synchrotron tracer from the 408MHz map to the 4.76GHz map. We interpret this as strong evidence for the origin of AME being spinning dust emission.

Andreia Carrillo, Keith Hawkins, Paula Jofré, Danielle de Brito Silva, Payel Das, Madeline Lucey

Accepted for publication in MNRAS. 25 pages, 14 figures, and 7 tables. Includes appendix A on line selection and B on NLTE corrections. Table of abundances will be available online but can be obtained through emailing the primary author

Understanding the assembly of our Galaxy requires us to also characterize the systems that helped build it. In this work, we accomplish this by exploring the chemistry of accreted halo stars from the Gaia-Enceladus/Gaia-Sausage (GES) selected in the infrared from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) Data Release 16. We use high resolution optical spectra for 62 GES stars to measure abundances in 20 elements spanning the $\alpha$, Fe-peak, light, odd-Z, and notably, the neutron-capture groups of elements to understand their trends in the context of and in contrast to the Milky Way and other stellar populations. Using these derived abundances we find that the optical and the infrared abundances agree to within 0.15 dex except for O, Co, Na, Cu, and Ce. These stars have enhanced neutron-capture abundance trends compared to the Milky Way, and their [Eu/Mg] and neutron-capture abundance ratios (e.g., [Y/Eu], [Ba/Eu], [Zr/Ba], [La/Ba], and [Nd/Ba]) point to r-process enhancement and a delay in s-process enrichment. Their [$\alpha$/Fe] trend is lower than the Milky Way trend for [Fe/H]$>$-1.5 dex, similar to previous studies of GES stars and consistent with the picture that these stars formed in a system with a lower rate of star formation. This is further supported by their depleted abundances in Ni, Na, and Cu abundances, again, similar to previous studies of low-$\alpha$ stars with accreted origins.

Marcella Wijngaarden, Katerina Chatziioannou, Andreas Bauswein, James A. Clark, Neil J. Cornish

18 pages, 17 figures

The gravitational wave signal emitted during the coalescence of two neutron stars carries information about the stars' internal structure. During the long inspiral phase the main matter observable is the tidal interaction between the binary components, an effect that can be parametrically modeled with compact-binary solutions to General Relativity. After the binary merger the main observable is frequency modes of the remnant, most commonly giving rise to a short-duration signal accessible only through numerical simulations. The complicated morphology and the decreasing detector sensitivity in the relevant frequencies currently hinder detection of the post-merger signal and motivate separate analyses for the pre-merger and post-merger data. However, planned and ongoing detector improvements could soon put the post-merger signal within reach. In this study we target the whole pre-merger and post-merger signal without an artificial separation at the binary merger. We construct a hybrid analysis that models the inspiral with templates based on analytical calculations and calibrated to numerical relativity and the post-merger signal with a flexible morphology-independent analysis. Applying this analysis to GW170817 we find, as expected, that the post-merger signal remains undetected. We further study simulated signals and find that we can reconstruct the full signal and simultaneously estimate both the pre-merger tidal deformation and the post-merger signal frequency content. Our analysis allows us to study neutron star physics using all the data available and directly test the pre-merger and post-merger signal for consistency thus probing effects such as the onset of the hadron-quark phase transition.

Fast magnetic reconnection is defined by the topology of the magnetic field lines changing on a time scale that is approximately an order of magnitude longer than the topology-conserving ideal-evolution time scale. Fast reconnection is an intrinsic property of Faraday's law when the evolving magnetic field depends non-trivially on all three spatial coordinates and is commonly observed -- even when the effects that allow topology breaking are arbitrarily small. The associated current density need only be enhanced by a factor of approximately ten but lies in thin but broad ribbons along the magnetic field. These results follow from the variation in the separation of neighboring pairs of magnetic field lines, which in an ideal evolution typically increases exponentially with time, and the existence of a spatial scale below which magnetic field lines feely change their identities due to non-ideal effects such a resistivity. Traditional reconnection theory ignores exponentially large variations and relies on the current density reaching a magnitude that is exponentially larger than is actually required. Here, an analysis of the behavior of magnetic field lines in the neighborhood of an arbitrarily chosen line is used to obtain more precise and rigorous results on intrinsic reconnection.

Subhajit Saha, Abdulla Al Mamon, Somnath Saha

9 pages, 1 figure, 3 tables, Comments welcome

This paper deals with the study of the evolution of PBHs in an FLRW universe with dissipation due to bulk viscosity which is considered to be in the form of isentropic gravitational particle creation. Assuming that the process of evaporation is quite suppressed during the early radiation era, we obtain an analytic solution for the evolution of PBH mass by accretion during this era, subject to an initial condition. We also obtain an upper bound on the accretion efficiency $\epsilon$ for $a \sim a_r$, where $a_r$ is the point of transition from the early de Sitter to the radiation era. Furthermore, aided by Maple 13, we obtain numerical solutions for the mass of a hypothetical PBH with intial mass 100 g assumed to be formed at an epoch when the value of the Hubble parameter was, say, 1 km/s/Mpc. We consider three values of the accretion efficiency, $\epsilon=0.23,0.5$, and $0.89$ for our study. The analysis reveals that the mass of the PBH increases rapidly due to the accretion of radiation in the early stages of its evolution. The accretion continues but its rate decreases gradually with the evolution of the Universe. Finally, Hawking radiation comes into play and the rate of evaporation surpasses the rate of accretion so that the PBH mass starts to decrease. As the Universe grows, evaporation becomes the dominant phenomenon, and the mass of the PBH decreases at a faster rate.

George Savvidy

79 pages, 17 figures. Based on lectures at the International Bogolyubov Conference "Problems of Theoretical and MathematicalPhysics" at the Steklov Mathematical Institute, Moscow-Dubna, September 9-13, 2019 [ arXiv:2001.01785 ] and seminars at the Niels Bohr Institute, at the CERN Theory Department and A. Alikhanian National Laboratory in Yerevan

We give a general review on the application of Ergodic theory to the investigation of dynamics of the Yang-Mills gauge fields and of the gravitational systems, as well as its application in the Monte Carlo method and fluid dynamics. In ergodic theory the maximally chaotic dynamical systems (MCDS) can be defined as dynamical systems that have nonzero Kolmogorov entropy. The hyperbolic dynamical systems that fulfil the Anosov C-condition belong to the MCDS insofar as they have exponential instability of their phase trajectories and positive Kolmogorov entropy. It follows that the C-condition defines a rich class of MCDS that span over an open set in the space of all dynamical systems. The large class of Anosov-Kolmogorov MCDS is realised on Riemannian manifolds of negative sectional curvatures and on high-dimensional tori. The interest in MCDS is rooted in the attempts to understand the relaxation phenomena, the foundations of the statistical mechanics, the appearance of turbulence in fluid dynamics, the non-linear dynamics of Yang-Mills field and gravitating N-body systems as well as black hole thermodynamics. Our aim is to investigate classical- and quantum-mechanical properties of MCDS and their role in the theory of fundamental interactions.

Pierre Dussarrat, Bertrand Theodore, Dorothee Coppens, Carsten Standfuss, Bernard Tournier

EUMETSAT is developing the on-ground processing chain of the infrared Fourier transform spectrometers (IRS) on-board of the Meteosat Third Generation sounding satellites (MTG-S). In this context, the authors have investigated the impact of a particular type of radiometric error, called the calibration ringing. It arises in Fourier transform spectrometers when the instrument optical transmission varies within the domain of the spectral response function. The expected radiometric errors are simulated in the context of the MTG-IRS instrument in the long wave infrared (LWIR) band. A software correction methodology is designed and its performance is assessed.

A. Lovato, I. Bombaci, D. Logoteta, M. Piarulli, R. B. Wiringa

15 pages, 3 figures. arXiv admin note: text overlap with arXiv:1908.04426

We present the equation of state of infinite neutron matter as obtained from highly-realistic Hamiltonians that include nucleon-nucleon and three-nucleon coordinate-space potentials. We benchmark three independent many-body methods: Brueckner-Bethe-Goldstone (BBG), Fermi hypernetted chain/single-operator chain (FHNC/SOC), and auxiliary-field diffusion Monte Carlo (AFDMC). We find them to provide similar equations of state when the Argonne $v_{18}$ and the Argonne $v_{6}^\prime$ nucleon-nucleon potentials are used in combination with the Urbana IX three-body force. Only at densities larger than about 1.5 the nuclear saturation density ($\rho_0 = 0.16\,\rm{fm}^{-3}$) the FHNC/SOC energies are appreciably lower than the other two approaches. The AFDMC calculations carried out with all of the Norfolk potentials fitted to reproduce the experimental trinucleon ground-state energies and $nd$ doublet scattering length yield unphysically bound neutron matter, associated with the formation of neutron droplets. Including tritium $\beta$-decay in the fitting procedure, as in the second family of Norfolk potentials, mitigates but does not completely resolve this problem. An excellent agreement between the BBG and AFDMC results is found for the subset of Norfolk interactions that do not make neutron-matter collapse, while the FHNC/SOC equations of state are moderately softer.

Axion-like particles (ALPs) are very light, neutral, spin zero bosons predicted by superstring theory. ALPs interact primarily with two photons and in the presence of an external magnetic field they generate photon-ALP oscillations and the change of the polarization state of photons. While well motivated from a theoretical point of view, hints on ALP existence come from astrophysics. In this paper, we demonstrate a strict relationship between initial photon polarization and photon-ALP conversion probability - which can be extrapolated by observed astrophysical spectra - so that, in the presence of ALPs, flux-measuring observatories become also porarimeters.

"El Cielo de Salamanca" ("The Sky of Salamanca") is a quarter-sphere-shaped vault 8.70 metres in diameter. It was painted sometime between 1480 and 1493 and shows five zodiacal constellations, three boreal and six austral. The Sun and Mercury are also represented. It formed part of a three times larger vault depicting the 48 Ptolemaic constellations and the rest of the planets known at the time. This was a splendid work of art that covered the ceiling of the first library of the University of Salamanca, one of the oldest in Europe having obtained its royal charter in 1218. But it was also a pioneering scientific work: a planetarium used to teach astronomy, the first of its kind in the history of Astronomy that has come down to us in the preserved part that we now call "The Sky of Salamanca". We describe the scientific context surrounding the chair of Astrology founded around 1460 at the University of Salamanca, which led to the production of this unique scientific work of art and to the flourishing of Astronomy in Salamanca. We analyse the possible dates compatible with it, showing that they are extremely infrequent. In the period of 1100 years from 1200 to 2300 that we studied there are only 23 years that have feasible days. We conclude that the information contained in "El Cielo de Salamanca" is not sufficient to assign it to a specific date but rather to an interval of several days that circumstantial evidence seems to place in August 1475. The same configuration of the sky will be observable, for the first time in 141 years, from the 22nd to the 25th of August 2022. The next occasion to observe it live will be in 2060.