Papers for Thursday, Feb 17 2022

The Milky Way halo was predominantly formed by the merging of numerous progenitor galaxies. However, our knowledge of this process is still incomplete, especially in regard to the total number of mergers, their global dynamical properties and their contribution to the stellar population of the Galactic halo. Here, we uncover the Milky Way mergers by detecting groupings of globular clusters, stellar streams and satellite galaxies in action ($\mathbf{J}$) space. While actions fully characterize the orbits, we additionally use the redundant information on their energy ($\textit{E}$) to enhance the contrast between groupings. For this endeavour, we use $\textit{Gaia}$ EDR3 based measurements of $170$ globular clusters, $41$ streams and $46$ satellites to derive their $\mathbf{J}$ and $\textit{E}$. To detect groups, we use the $\texttt{ENLINK}$ software, coupled with a statistical procedure that accounts for the observed phase-space uncertainties of these objects. We detect a total of $N=6$ groups, including the previously known mergers $\textit{Sagittarius}$, $\textit{Cetus}$, $\textit{Gaia-Sausage/Enceladus}$, $\textit{LMS-1/Wukong}$, $\textit{Arjuna/Sequoia/I'itoi}$ and one new merger that we call $\textit{Pontus}$. All of these mergers, together, comprise $62$ objects ($\approx 25\%$ of our sample). We discuss their members, orbital properties and metallicity distributions. We find that the three most metal-poor streams of our Galaxy -- "C-19" ([Fe/H]$=-3.4$ dex), "Sylgr" ([Fe/H]$=-2.9$ dex) and "Phoenix" ([Fe/H]$=-2.7$ dex) -- are associated with $\textit{LMS-1/Wukong}$; showing it to be the most metal-poor merger. The global dynamical atlas of Milky Way mergers that we present here provides a present-day reference for galaxy formation models.

A mild correlation exists in active galaxies between the mean black hole accretion, as traced by the mean X-ray luminosity $\left<{\rm L}_{\rm X}\right>$, and the host galaxy stellar mass M$_*$, characterised by a normalisation steadily decreasing with cosmic time and lower in more quiescent galaxies. We create comprehensive semi-empirical mock catalogues of active black holes to pin down which parameters control the shape and evolution of the $\left<{\rm L}_{\rm X}\right>-{\rm M}_*$ relation of X-ray detected active galaxies. We find that the normalisation of the $\left<{\rm L}_{\rm X}\right>-{\rm M}_*$ relation is largely independent of the fraction of active galaxies (the duty cycle), but strongly dependent on the mean Eddington ratio, when adopting a constant underlying M$_{\rm BH}-{\rm M}_*$ relation as suggested by observational studies. The data point to a decreasing mean Eddington ratio with cosmic time and with galaxy stellar mass at fixed redshift. Our data can be reproduced by black holes and galaxies evolving on similar M$_{\rm BH}-{\rm M}_*$ relations but progressively decreasing their average Eddington ratios, mean X-ray luminosities, and specific star formation rates, when moving from the starburst to the quiescent phase. Models consistent with the observed $\left<{\rm L}_{\rm X}\right>-{\rm M}_*$ relation and independent measurements of the mean Eddington ratios, are characterised by M$_{\rm BH}-{\rm M}_*$ relations lower than those derived from dynamically measured local black holes. Our results point to the $\left<{\rm L}_{\rm X}\right>-{\rm M}_*$ relation as a powerful diagnostic to: 1) probe black hole-galaxy scaling relations and the level of accretion onto black holes; 2) efficiently break the degeneracies between duty cycles and accretion rates in cosmological models of black holes.

Recent observations of the stellar halo have uncovered the debris of an ancient merger, Gaia-Sausage-Enceladus, estimated to have occurred ~8 Gyr ago. Follow-up studies have associated GSE with a large-scale tilt in the stellar halo that links two well-known stellar over-densities in diagonally opposing octants of the Galaxy (the Hercules-Aquila Cloud and Virgo Overdensity; HAC and VOD). In this paper, we study the plausibility of such unmixed merger debris persisting over several Gyr in the Galactic halo. We employ the simulated stellar halo from Naidu et al. (2021), which reproduces several key properties of the merger remnant, including the large-scale tilt. By integrating the orbits of these simulated stellar halo particles, we show that adoption of a spherical halo potential results in rapid phase mixing of the asymmetry. However, adopting a tilted halo potential preserves the initial asymmetry in the stellar halo for many Gyr. The asymmetry is preserved even when a realistic growing disk is added to the potential. These results suggest that HAC and VOD are long-lived structures that are associated with GSE and that the dark matter halo of the Galaxy is tilted with respect to the disk and aligned in the direction of HAC-VOD. Such halo-disk misalignment is common in modern cosmological simulations. Lastly, we study the relationship between the local and global stellar halo in light of a tilted global halo comprised of highly radial orbits. We find that the local halo offers a dynamically biased view of the global halo due to its displacement from the Galactic Center.

We present GIGA-Lens: a gradient-informed, GPU-accelerated Bayesian framework for modeling strong gravitational lensing systems, implemented in TensorFlow and JAX. The three components, optimization using multi-start gradient descent, posterior covariance estimation with variational inference, and sampling via Hamiltonian Monte Carlo, all take advantage of gradient information through automatic differentiation and massive parallelization on graphics processing units (GPUs). We test our pipeline on a large set of simulated systems and demonstrate in detail its high level of performance. The average time to model a single system on four Nvidia A100 GPUs is 105 seconds. The robustness, speed, and scalability offered by this framework make it possible to model the large number of strong lenses found in current surveys and present a very promising prospect for the modeling of $\mathcal{O}(10^5)$ lensing systems expected to be discovered in the era of the Vera C. Rubin Observatory, Euclid, and the Nancy Grace Roman Space Telescope.

Vector resonant relaxation (VRR) is known to be the fastest gravitational process that shapes the geometry of stellar orbits in nuclear star clusters. This leads to the realignment of the orbital planes on the corresponding VRR time scale $t_{\rm VRR}$ of a few million years, while the eccentricity $e$ and semimajor axis $a$ of the individual orbits are approximately conserved. The distribution of orbital inclinations reaches an internal equilibrium characterised by two conserved quantities, the total potential energy among stellar orbits, $E_{\rm tot}$, and the total angular momentum, $L_{\rm tot}$. On timescales longer than $t_\mathrm{VRR}$, the eccentricities and semimajor axes change slowly and the distribution of orbital inclinations are expected to evolve through a series of VRR equilibria. Using a Monte Carlo Markov Chain method, we determine the equilibrium distribution of orbital inclinations in the microcanonical ensemble with fixed $E_{\rm tot}$ and $L_{\rm tot}$ for isolated nuclear star clusters with a power-law distribution of $a$, $e$, and $m$, where $m$ is the stellar mass. We explore the possible equilibria for $9$ representative $E_{\rm tot}$--$L_{\rm tot}$ pairs that cover the possible parameter space. For all cases, the equilibria show anisotropic mass segregation where the distribution of more massive objects is more flattened than that for lighter objects. Given that stellar black holes are more massive than the average main sequence stars, these findings suggest that black holes reside in disc-like structures within nuclear star clusters for a wide range of initial conditions.

We analyse FIR dust continuum measurements for 14 galaxies ($z\approx 7$) in the ALMA REBELS LP to derive their physical properties. Our model uses three input data: (a) the UV spectral slope, $\beta$, (b) the observed UV continuum flux at $1500$A, $F_{\rm UV}$, (c) the observed continuum flux at $\approx 158\mu$m, $F_{158}$, and considers Milky Way (MW) and SMC extinction curves, along with different dust geometries. We find that REBELS galaxies have (28-90.5)% of their star formation obscured; the total (UV+IR) star formation rates are in the range $31.5 < {\rm SFR}/ (M_\odot {\rm yr}^{-1}) < 129.5$. The sample-averaged dust mass and temperature are $(1.3\pm 1.1)\times 10^7 M_\odot$ and $52 \pm 11$ K, respectively. In some galaxies dust is abundant (REBELS-14, $M'_d \approx 3.4 \times 10^7 M_\odot$), or hot (REBELS-18, $T'_d \approx 67$ K). The dust distribution is compact ($<0.3$ kpc for 70% of the galaxies). The dust yield per supernova is $0.1 \le y_d/M_\odot \le 3.3$, with 70% of the galaxies requiring $y_d < 0.25 M_\odot$. Three galaxies (REBELS-12, 14, 39) require $y_d > 1 M_\odot$. With the SFR predicted by the model and a MW extinction curve, REBELS galaxies detected in [CII] nicely follow the local $L_{\rm CII}-$SFR relation, and are approximately located on the Kennicutt-Schmidt relation. The sample-averaged gas depletion time is of $0.11\, y_P^{-2}$ Gyr, where $y_P$ is the ratio of the gas-to-stellar distribution radius. For some systems a solution simultaneously matching the observed ($\beta, F_{\rm UV}, F_{158}$) values cannot be found. This occurs when the index $I_m = (F_{158}/F_{\rm UV})/(\beta-\beta_{\rm int})$, where $\beta_{\rm int}$ is the intrinsic UV slope, exceeds $I_m^*\approx 1120$ for a MW curve. For these objects we argue that the FIR and UV emitting regions are not co-spatial, questioning the use of the IRX-$\beta$ relation.

We present VLA observations toward the Class 0 protostar L1157 MMS at 6.8 mm and 9 mm with a resolution of ~0.04" (14 au). We detect two sources within L1157 MMS and interpret these sources as a binary protostar with a separation of ~16 au. The material directly surrounding the binary system within the inner 50 au radius of the system has an estimated mass of 0.11 M_sun, calculated from the observed dust emission. We interpret the observed binary system in the context of previous observations of its flattened envelope structure, low rates of envelope rotation from 5000 to 200 au scales, and an ordered, poloidal magnetic field aligned with the outflow. Thus, L1157 MMS is a prototype system for magnetically-regulated collapse and the presence of a compact binary within L1157 MMS demonstrates that multiple star formation can still occur within envelopes that likely have dynamically important magnetic fields.

Stars can be tidally destroyed or swallowed by supermassive black hole binaries. Using a large number of accurate few-body simulations, we investigate the enhancement and suppression of full and partial disruption and direct capture events by hard supermassive black hole binaries with a wide ranges of key parameters, i.e., the primary black hole mass ($10^{5}-10^{8}M_{\odot}$), the binary mass ratio ($10^{-3}-1$), the ratio of the binary semimajor axis to the hardening radius ($10^{-4}-1$), the binary eccentricity ($0.0-0.9$) and the stellar mass ($0.3-3M_{\odot}$). This is a significant extension of the parameter space compared to previous work. We show that the encounter probabilities of all three events are well-described by the encounter cross section, which is proportional to the pericenter distance. The probability of full disruptions by supermassive black hole binaries can be enhanced by up to a factor of $40-50$ or suppressed by up to a factor of $10$, relative to that by single black holes, depending on the binary parameters. The Post-Newtonian effects are not important for the primary black hole mass $\leq 10^{7}M_{\odot}$, but can provide an additional enhancement of the full disruption probability by less than a factor of $2-3$ for higher primary black hole masses. We provide a fitting formula for the full disruption probability by the hard supermassive black hole binaries that works for a wide range of parameters. We also find that partial disruption events can occur multiple times before full disruptions or direct captures, and their probabilities can be greater than that of full disruption events by a factor of three. Because partial disruption events can induce stellar spins and mass loss and change the orbits, it can significantly affect the overall full disruption event rate and the shape of the light curves.

Using a penalised maximum likelihood we estimate, for the first time, the velocity distribution of white dwarfs in the Solar neighbourhood. Our sample consists of 129 675 white dwarfs within 500 pc in Gaia Early Data Release 3 The white dwarf velocity distributions reveal a similar structure to the rest of the Solar neighbourhood stars, reflecting that white dwarfs are subjected to the same dynamical processes. In the velocity distribution for three magnitude-binned subsamples we however find a novel structure at $(U, V) = (7, -19)$ km s$^{-1}$ in fainter samples, potentially related to the Coma Berenices stream. We also see a double-peaked feature in $U$-$W$ at $U \approx -30$ km s$^{-1}$ and in $V$-$W$ at $V \approx -20$ km s$^{-1}$ for fainter samples. We determine the velocity distribution and velocity moments as a function of absolute magnitude for two samples based on the bifurcation identified in Gaia Data Release 2 in the colour-magnitude diagram. The brighter, redder sequence has a larger velocity dispersion than the fainter, bluer sequence across all magnitudes. It is hard to reconcile this kinematic difference with a bifurcation caused purely by atmospheric composition, while it fits neatly with a significant age difference between the two sequences. Our results provide novel insights into the kinematic properties of white dwarfs and demonstrate the power of analytical techniques that work for the large fraction of stars that do not have measured radial velocities in the current era of large-scale astrometric surveys.

X-ray observations of the hot intra-cluster medium (ICM) in galaxy groups and clusters provide quantities such as their gas mass, X-ray luminosity, and temperature. The analysis of the scaling relations between these observable properties gives considerable insight into the physical processes taking place in the ICM. Furthermore, an understanding of the scaling relations between ICM properties and the total cluster mass is essential for cosmological studies with clusters. For these reasons, the X-ray scaling relations of groups and clusters have been a major focus of research over the past several decades. In this Chapter, after presenting the expectations from the self-similar model, based on the assumption that only gravity drives the evolution of the ICM, we discuss how the processes of gas cooling and non-gravitational heating are believed to be responsible for the observed deviations from the self-similar scenario. We also describe important complications that must be considered when measuring and interpreting the scaling relations.

Single-degenerate (SD) binary systems composed of a white dwarf and a non-degenerate helium (He)-star companion have been proposed as the potential progenitors of Type Ia supernovae (SNe Ia). The He-star companions are expected to survive the SN Ia explosion in this SD progenitor model. In the present work, we map the surviving He-star companion models computed from our previous three-dimensional hydrodynamical simulations of ejecta-companion interaction into the one-dimensional stellar evolution code MESA to follow their long-term evolution to make predictions on their post-impact observational properties, which can be helpful for searches of such surviving He-star companions in future observations. By comparing with the very late-epoch light curve of the best observed SN Ia, SN 2011fe, we find that our surviving He-star companions become significantly more luminous than SN 2011fe about 1000d after the maximum light. This suggests that a He star is very unlikely to be a companion to the progenitor of SN 2011fe.

The level of central segregation of Blue Straggler stars proved to be an excellent tracer of the dynamical evolution of old star clusters (the so-called "dynamical clock"), both in the Milky Way and in the Large Magellanic Cloud. The $A^{+}$ parameter, used to measure the Blue Stragglers degree of segregation, has in fact been found to strongly correlate with the parent cluster central relaxation time. Here we studied the Blue-Straggler population of two young stellar systems in the Small Magellanic Cloud, namely NGC 339 (which is 6 Gyr old) and NGC 419 (with an age of only 1.5 Gyr), in order to study their dynamical state. Thanks to multi-epoch, high angular resolution Hubble Space Telescope observations available for both clusters, we took advantage of the stellar proper motions measured in the regions of the two systems and we selected a population of likely cluster members, removing the strong contamination from Small Magellanic Cloud stars. This enabled us to study, with unprecedented accuracy, the radial distribution of Blue Stragglers in these two extragalactic clusters and to measure their dynamical age. As expected for such young clusters, we found that both systems are poorly evolved from the dynamical point of view, also fully confirming that the $A^{+}$ parameter is a sensitive "clock hand" even in the dynamically-young regime.

We present an algorithm to detect the outer edges of Coronal Mass Ejection (CME) events as seen in differences of Heliospheric Imager STEREO SECCHI HI-1 images from either A or B spacecraft, as well as its implementation in Python.

The power spectrum has been a workhorse for cosmological studies of large-scale structure. However, the present-day matter distribution is highly non-Gaussian and significant cosmological information is also contained in higher-order correlation functions. Meanwhile, baryon physics (particularly AGN feedback) has previously been shown to strongly affect the two-point statistics but there has been limited exploration of its effects on higher-order functions to date. Here we use the BAHAMAS suite of cosmological hydrodynamical simulations to explore the effects of baryon physics and massive neutrinos on the halo bispectrum. In contrast to matter clustering which is suppressed by baryon physics, we find that the halo clustering is typically enhanced. The strength of the effect and the scale over which it extends depends on how haloes are selected. On small scales (k > 1 $h$ Mpc$^{-1}$, dominated by satellites of groups/clusters), we find that the bispectrum is highly sensitive to the efficiency of star formation and feedback, making it an excellent testing ground for galaxy formation models. We show that the effects of feedback and the effects of massive neutrinos are largely separable (independent of each other) and that massive neutrinos strongly suppress the halo bispectrum on virtually all scales up to the free-streaming length (apart from the smallest scales, where baryon physics dominates). The strong sensitivity of the bispectrum to neutrinos on the largest scales and galaxy formation physics on the smallest scales bodes well for upcoming precision measurements from the next generation of wide-field surveys.

Classical novae are cataclysmic binary star systems in which the matter of a companion star is accreted on a white dwarf (WD). Accumulation of the matter in a layer eventually causes a thermonuclear explosion on the surface of the WD, brightening the WD to ~ 10 5 solar luminosities and triggering ejection of the accumulated matter. They provide extreme conditions required to accelerate particles, electrons or protons, to high energies. Here we present the detection of gamma rays by the MAGIC telescopes from the 2021 outburst of RS Ophiuchi (RS Oph), a recurrent symbiotic nova, that allowed us, for the first time, to accurately characterize the emission from a nova in the 60 GeV to 250 GeV energy range. The theoretical interpretation of the combined Fermi -LAT and MAGIC data suggests that protons are accelerated to hundreds of GeV in the nova shock. Such protons should create bubbles of enhanced Cosmic Ray density up to about 13 pc from the recurrent novae.

We conduct three-dimensional hydrodynamical simulations of weak jets that we launch into a core collapse supernovae (CCSNe) ejecta half an hour after the explosion and find that the interaction of the fast jets with the CCSN ejecta creates high pressure zones that induce a backflow that results in mass accretion onto the newly born neutron star. In cases of weak jets, a total power of 10^45-10^46 erg, the backflow mass accretion might power up to an order of magnitude more energetic jets. In total, the jets of the two post-explosion jet-launching episodes have enough energy to influence the morphology of the very inner ejecta, a mass of 0.1 M_O. Our results imply that in some, probably a minority of, CCSN remnants the very inner regions might display a bipolar structure that results from post-explosion weak jets. The regions outside this part might display the morphology of jittering jets.

In this paper we present the third data release (DR3) of the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS). We identify likely near-infrared counterparts to submillimetre sources in the South Galactic Pole (SGP) field using the VISTA VIKING survey. We search for the most probable counterparts within 15 arcsec of each Herschel source using a probability measure based on the ratio between the likelihood the true counterpart is found close to the submillimetre source and the likelihood that an unrelated object is found in the same location. For 110 374 (57.0$\%$) sources we find galaxies on the near-infrared images where the probability that the galaxy is associated to the source is greater than 0.8. We estimate the false identification rate to be 4.8$\%$, with a probability that the source has an associated counterpart on the VIKING images of 0.835$\pm$0.009. We investigate the effects of gravitational lensing and present 41 (0.14 deg$^{-2}$) candidate lensed systems with observed flux densities > 100 mJy at 500 $\mu$m. We include in the data release a probability that each source is gravitationally lensed and discover an additional 5 923 sources below 100 mJy that have a probability greater than 0.94 of being gravitationally lensed. We estimate that $\sim$ 400 - 1 000 sources have multiple true identifications in VIKING based on the similarity of redshift estimates for multiple counterparts close to a Herschel source. The data described in this paper can be found at the H-ATLAS website.

We studied the spectral changes of the high-mass X-ray binary system LMC X-4 to understand the origin and mechanisms beyond its super-orbital modulation (30.4 days). To this aim, we obtained a monitoring campaign with Swift/XRT (0.3-10 keV) and complemented these data with the years-long Swift/BAT survey data (15-60 keV). We found a self-consistent, physically motivated, description of the broadband X-ray spectrum using a Swift/XRT and a NuSTAR observation at the epoch of maximum flux. We decomposed the spectrum into the sum of a bulk+thermal Comptonization, a disk-reflection component and a soft contribution from a standard Shakura-Sunyaev accretion disk. We applied this model to 20 phase-selected Swift spectra along the super-orbital period. We found a phase-dependent flux ratio of the different components, whereas the absorption column does not significantly vary. The disk emission is decoupled with respect to the hard flux. We interpret this as a geometrical effect in which the inner parts of the disk are tilted with respect to the obscuring outer regions.

We investigate the model where electrons and dark matter interact with dark energy through the rolling of a scalar field which comes from extra dimensional theories such as the braneworld theory and Brans-Dicke theory. In this model, dark energy couples to dark matter and electrons, which leads to larger values of the mass energies of dark matter and electrons in the early universe. We also fit our model to the cosmological data. By analyzing the data from Planck, baryon acoustic oscillation (BAO), light curves (Pantheon), and type-Ia supernovae (SH0ES), it can be seen that the Hubble tension is relieved in our model and the coupling parameter prefers a non-zero value with a significance of over 2{\sigma}.

We present 12 new transit light curves and 16 new out-of-transit radial velocity measurements for the XO-3 system. By modelling our newly collected measurements together with archival photometric and Doppler velocimetric data, we confirmed the unusual configuration of the XO-3 system, which contains a massive planet ($M_P=11.92^{+0.59}_{-0.63} M_J$) on a relatively eccentric ($e=0.2853^{+0.0027}_{-0.0026}$) and short-period ($3.19152 \pm 0.00145\,$day) orbit around a massive star ($M_*=1.219^{+0.090}_{-0.095} M_{\odot}$). Furthermore, we find no strong evidence for a temporal change of either $V\sin i_{*}$ (and by extension, the stellar spin vector of XO-3), or the transit profile (and thus orbital angular momentum vector of XO-3b). We conclude that the discrepancy in previous Rossiter-McLaughlin measurements ($70.0^{\circ} \pm 15.0^{\circ}$ (Hebrard et al. 2008); $37.3^{\circ} \pm 3.7^{\circ}$ (Winn et al. 2009); $37.3^{\circ} \pm 3.0^{\circ}$ (Hirano et al. 2011)) may have stemmed from systematic noise sources.

We present optical ccd images of the large supernova remnant (SNR) G132.7$+$1.3 (HB3) covering its full extent for the first time, in the emission lines of H$\alpha+$[N II], [S II] and [O III], where new and known filamentary and diffuse structures are detected. These observations are supplemented by new low-resolution long-slit spectra and higher-resolution images in the same emission lines. Both the flux-calibrated images and spectra confirm that the optical emission originates from shock-heated gas since the [S II]/H$\alpha$ $>$ 0.4. Our findings are also consistent with the recently developed emission line ratio diagnostics for distinguishing SNRs from H II regions. A multi-wavelength comparison among our optical data and relevant observations in radio, X-rays, $\gamma$-rays and CO bands, provided additional evidence on the interaction of HB3 with the surrounding clouds and clarified the borders of the SNR and the adjacent cloud. We discuss the supernova (SN) properties and evolution that led to the current observables of HB3 and we show that the remnant has most likely passed at the pressure driven snowplow phase. The estimated SN energy was found to be $\left(3.7 \pm 1.5\right) \times 10^{51}$ erg and the current SNR age $\left(5.1 \pm 2.1\right) \times 10^4$ yrs. We present an alternative scenario according to which the SNR evolved in the wind bubble cavity excavated by the progenitor star and currently is interacting with its density walls. We show that the overall mixed morphology properties of HB3 can be explained if the SN resulted by a Wolf-Rayet progenitor star with mass $\sim 34 \rm~M_{\rm\odot}$.

Using the Space Telescope Imaging Spectrograph, we have obtained ultraviolet (UV) spectra from $\sim 1200$ to 2000 \AA\ of known Lyman continuum (LyC) emitting galaxies at low redshift ($z \sim 0.3-0.4$) with varying absolute LyC escape fractions (fesc $\sim 0.01 - 0.72$). Our observations include in particular the galaxy J1243+4646, which has the highest known LyC escape fraction at low redshift. While all galaxies are known Lyman alpha emitters, we consistently detect an inventory of additional emission lines, including CIV 1550, HeII 1640, OIII] 1666, and CIII] 1909, whose origin is presumably essentially nebular. CIV 1550 emission is detected above 4 $\sigma$ in six out of eight galaxies, with equivalent widths of EW(CIV)$=12-15$ Ang for two galaxies, which exceeds the previously reported maximum emission in low-$z$ star-forming galaxies. We detect CIV 1550 emission in all LyC emitters with escape fractions fesc $> 0.1$ and find a tentative increase in the flux ratio CIV 1550/ CIII] 1909 with fesc. Based on the data, we propose a new criterion to select and classify strong leakers (galaxies with fesc $> 0.1$): CIV 1550/ CIII] 1909 $> 0.75$. Finally, we also find HeII 1640 emission in all the strong leakers with equivalent widths from 3 to 8 Ang rest frame. These are among the highest values observed in star-forming galaxies and are primarily due to a high rate of ionizing photon production. The nebular HeII 1640 emission of the strong LyC emitters does not require harder ionizing spectra at $>54$ eV compared to those of typical star-forming galaxies at similarly low metallicity.

Detailed knowledge about stellar winds and evolution at different metallicities is crucial for understanding stellar populations and feedback in the Local Group of galaxies and beyond. Despite efforts in the literature, we still lack a comprehensive, empirical view of the dependence of wind properties on metallicity ($Z$). Here, we investigate the winds of O and B stars in the Milky Way (MW) and Small Magellanic Cloud (SMC). We gathered a sample of 96 stars analyzed by means of the NLTE code CMFGEN. We explored their wind strengths and terminal velocities to address the $Z$ dependence, over a large luminosity range. The empirical wind-luminosity relation (WLR) obtained updates and extends previous results in the literature. It reveals a luminosity and $Z$ dependence, in agreement with the radiatively driven wind theory. For bright objects ($\log L/L_\odot \gtrsim 5.4$), we infer that $\dot{M} \sim Z^{0.5-0.8}$. However, this dependence seems to get weaker or vanish at lower luminosities. The analysis of the terminal velocities suggests a shallow $Z^n$ dependence, with $n \sim 0.1-0.2$, but it should be confirmed with a larger sample and more accurate $V_{\infty}$ determinations. Recent results on SMC stars based on the PoWR code support our inferred WLR. On the other hand, recent bow-shocks measurements stand mostly above our derived WLR. Theoretical calculations of the WLR are not precise, specially at low $L$, where the results scatter. Deviations between our results and recent predictions are identified to be due to the weak wind problem and the extreme terminal velocities predicted by the models. The Z dependence suggested by our analysis deserves further investigations, given its astrophysical implications.

These workshop notes present an introduction to the concepts and mathematical foundations of polarimetry. One of the main goals of this workshop is to develop an understanding of the relationships between a physical description of the signal path (e.g. gain, delay, rotation, coupling, etc.), the corresponding transformations of the electric field vector, and the equivalent transformations of the Stokes parameters. The adopted algebraic/geometric approach is either directly copied from or heavily inspired by the work of Britton (2000) and Hamaker (2000), and some justification is provided for preferring these over other approaches and parameterizations.

We present the orbital period variations of 14,127 contact eclipsing binaries (CEBs) based on the OGLE-III\&IV observations in the Galactic bulge. New times of minimum lights for the CEBs were derived by binary modeling for the full seasonal light curves, which were made from survey observations at an interval of 1 year. The orbital period changes of the systems were classified based on the statistical inference, multiple-hypothesis testing error measure, and visual inspection of the eclipse timing diagrams. As the results, we identified 13,716 CEBs with a parabola, 307 CEBs with a sinusoid, and 104 CEBs with the two variations. The period distributions of the inner close binaries and the outer companions were in the ranges of $0.235-0.990$ days and $5.0-14.0$ years, respectively. In our sample of 13,820 CEBs showing a parabolic variation, the highest decreasing and increasing period rates were determined to be $\dot P=-1.38\pm0.06\times10^{-5}$ day year$^{-1}$ for OGLE-BLG-ECL-169991 and $\dot P=+8.99\pm0.44\times10^{-6}$ day year$^{-1}$ for OGLE-BLG-ECL-189805, respectively. The secular period change rates were distributed almost symmetrically around zero, and most of them lie within $\dot P=\pm5.0\times10^{-6}$ day year$^{-1}$.

Einstein's theory of general relativity (GR) has been precisely tested on solar system scales, but extragalactic tests are still poorly performed. In this work, we use a newly compiled sample of galaxy-scale strong gravitational lenses to test the validity of GR on kiloparsec scales. In order to solve the circularity problem caused by the pre-assumption of a specific cosmological model based on GR, we employ the distance sum rule in the Friedmann-Lema\^{\i}tre-Robertson-Walker metric to directly estimate the parameterized post-Newtonian (PPN) parameter $\gamma_{\rm PPN}$ and the cosmic curvature $\Omega_k$ by combining observations of strong lensing and type Ia supernovae. This is the first simultaneous measurement of $\gamma_{\rm PPN}$ and $\Omega_k$ without any assumptions about the contents of the universe or the theory of gravity. Our results show that $\gamma_{\rm PPN}=1.11^{+0.11}_{-0.09}$ and $\Omega_{k}=0.48^{+1.09}_{-0.71}$, indicating a strong degeneracy between the two quantities. The measured $\gamma_{\rm PPN}$, which is consistent with the prediction of 1 from GR, provides a precise extragalactic test of GR with a fractional accuracy better than 9.0\%. If a prior of spatial flatness (i.e., $\Omega_{k}=0$) is adopted, the PPN parameter constraint can be further improved to $\gamma_{\rm PPN}=1.07^{+0.07}_{-0.07}$, representing a precision of 6.5\%. On the other hand, in the framework of a GR (i.e., $\gamma_{\rm PPN}=1$), our results are still marginally compatible with zero curvature ($\Omega_k=-0.12^{+0.48}_{-0.36}$), supporting no significant deviation from a flat universe.

We use the age measurements of 114 old astrophysical objects (OAO) in the redshift range $0\lesssim z\lesssim 8$ to explore the Hubble tension. The age of the Universe at any $z$ is inversely proportional to the Hubble constant, $H_0$, so requiring the Universe to be older than the OAO it contains at any $z$ will lead to an upper limit on $H_0$. Assuming flat $\Lambda$CDM and setting a Gaussian prior on the matter density parameter $\Omega_{\rm m}=0.315\pm0.007$ informed by {\it Planck}, we obtain a 95\% confidence-level upper limit of $H_0<70.6 \rm{~km} \rm{~s}^{-1} \rm{~Mpc}^{-1}$, representing a $2\sigma$ tension with the measurement using the local distance ladder. We find, however, that the inferred upper limit on $H_{0}$ depends quite sensitively on the prior for $\Omega_{\rm m}$, and the Hubble tension between early-time and local measurements of $H_{0}$ may be due in part to the inference of both $\Omega_{\rm m}$ and $H_0$ in {\it Planck}, while the local measurement uses only $H_{0}$. The age-redshift data may also be used for cosmological model comparisons. We find that the $R_{\rm h}=ct$ universe accounts well for the data, with a reasonable upper limit on $H_{0}$, while Einstein-de Sitter fails to pass the cosmic-age test. Finally, we present a model-independent estimate of the spatial curvature using the ages of 61 galaxies and the luminosity distances of 1,048 Pantheon Type Ia supernovae. This analysis suggests that the geometry of the Universe is marginally consistent with spatial flatness at a confidence level of $1.6\sigma$, characterized as $\Omega_{k}=0.43^{+0.27}_{-0.27}$.

The existence of magnetized turbulence in the interstellar HI is well accepted. A number of techniques to obtain turbulence spectrum and magnetic field direction and strength have been developed and successfully applied to HI spectroscopic data. To better separate the imprints of density and velocity fluctuations to the channel maps, a new theory-based technique, the Velocity Decomposition Algorithm (VDA,Yuen et.al 2021), has been created. The technique demonstrates that the intensity fluctuations are separated into a component pv that mostly arises from velocity fluctuations and pd that mostly arise from density fluctuations. The VDA helps to clarify the nature of the filamentary structure observed in channel maps. A recent publication (Kalberla et.al 2022,K22) claims that the application of VDA to HI4PI data provides negative correlation of pv and pd,which according to the authors invalidates the technique since it requires that pv and pd have zero correlation. However, the quantities pv and pd given by VDA are naturally orthogonal which can be trivially checked analytically or numerically. That means the correct application of the VDA to any data must provide zero correlation. This is the point that we clarify in this paper and search for the cause of the mistake in the application of the VDA in K22 that resulted in the erroneous conclusion. We prove analytically that by construction pv and pd are not correlated. We identify the likely mistake in the VDA expression that K22 used and reproduce their figures with the incorrect expression. We find that 14 out of 15 figures in K22 are invalid, and thus their criticism of the VDA is ill-founded and arises from their use of incorrect expressions. We conclude that the detrimental mistake that K22 made at their analysis completely invalidate their scientific claim that Y21 isnt compatible to observation.

The cosmic curvature density parameter has been constrained in the present work independent of any background cosmological model. The reconstruction is performed adopting the non-parametric Gaussian Processes (GP). The constraints on $\Omega_{k0}$ are obtained via a Markov Chain Monte Carlo (MCMC) analysis. Late-time cosmological probes viz., the Supernova (SN) distance modulus data, the Cosmic Chronometer (CC) and the radial Baryon Acoustic Oscillations ($r$BAO) measurements of the Hubble data have been utilized for this purpose. The results are further combined with the data from redshift space distortions (RSD) which studies the growth of large scale structure in the universe. The only \textit{a priori} assumption is that the universe is homogeneous and isotropic, described by the FLRW metric. Results indicate that a spatially flat universe is well consistent in 2$\sigma$ within the domain of reconstruction $0<z<2$ for the background data. On combining the RSD data we find that the results obtained are consistent with spatial flatness mostly within 2$\sigma$ and always within 3$\sigma$ in the domain of reconstruction $0<z<2$.

Aims. The first relativistic solar proton event of solar cycle 25 (SC25) was detected on 28 October 2021 by neutron monitors (NMs) on the ground and particle detectors onboard spacecraft in the near-Earth space. This is the first ground level enhancement (GLE) of the current cycle. A detailed reconstruction of the NM response together with the identification of the solar eruption that generated these particles is investigated based on in-situ and remote-sensing measurements. Methods. In-situ proton observations from a few MeV to $\sim$500 MeV were combined with the detection of a solar flare in soft X-rays (SXRs), a coronal mass ejection (CME), radio bursts and extreme ultraviolet (EUV) observations to identify the solar origin of the GLE. Timing analysis was performed and a relation to the solar sources was outlined. Results. GLE73 reached a maximum particle rigidity of $\sim$2.4 GV and is associated with type III, type II, type IV radio bursts and an EUV wave. A diversity of time profiles recorded by NMs was observed. This points to an anisotropic nature of the event. The peak flux at E$>$10 MeV was only $\sim$30 pfu and remained at this level for several days. The release time of $\geq$1 GV particles was found to be $\sim$15:40 UT. GLE73 had a moderately hard rigidity spectrum at very high energies ($\gamma$ $\sim$5.5). Comparison of GLE73 to previous GLEs with similar solar drivers is performed

In this work, starting from the well-accepted relations in literature, we introduce a new formalism to compute the astrometric membership probabilities for sources in star clusters, and we provide an application to the case of the open cluster M 37. The novelty of our approach is a refined -- and magnitude-dependent -- modelling of the parallax distribution of the field stars. We employ the here-derived list of members to estimate the cluster's mean systemic astrometric parameters, which are based on the most recent Gaia's catalog (EDR3).

Photochemistry has the potential to substantially impact the atmospheric composition of exoplanets with consequences on the radiative transfer, thermal structure and dynamics of the atmospheres, particularly in UV-rich stellar environments. Here, we present the results of a first laboratory experimental simulation of photochemistry in carbon-rich exoplanet atmospheres at elevated temperatures. Evolution of gas-phase molecular composition was quantitatively monitored with infrared spectroscopy and mass spectrometry. We found that H2/CO gas compositions can change significantly from thermal equilibria compositions when irradiated with Lyman-alpha photons at temperatures ranging from 600 K to 1500 K. Carbon dioxide and water were found to be the main products caused by photolysis, while formation of methane was also observed to a lesser extent. We find that photochemistry efficiency is strongly correlated with increasing temperature. Our finding that water is efficiently produced by photochemistry in a super Solar C/O=1 environment, representing C enhancement relative to solar values C/O ratio = 0.54, has significant implications for the interpretation of many exoplanet transmission spectra. We also find the formation of an organic solid condensate at 1500 K and under Lyman-alpha UV-radiation, confirming the possibility of forming photochemical hazes in hot-Jupiter exoplanet atmospheres with an enhanced C/O ratio compared to Solar.

Classical Cepheids are the most popular distance indicators and tracers of young stellar populations. The key advantage is that they are bright and they can be easily identified in Local Group and Local Volume galaxies. Their evolutionary and pulsation properties depend on their chemical abundances. The main aim of this investigation is to perform a new and accurate abundance analysis of two tens of calibrating Galactic Cepheids using high spectral resolution (R$\sim$40,000-115,000) and high S/N spectra ($\sim$400) covering the entire pulsation cycle. We focus our attention on possible systematics affecting the estimate of atmospheric parameters and elemental abundances along the pulsation cycle. We cleaned the line list by using atomic transition parameters based on laboratory measurements and by removing lines that are either blended or display abundance variations along the pulsation cycle. The spectroscopic approach that we developed brings forward small dispersions in the variation of the atmospheric parameters ($\sigma$($T_{\rm eff}$)$\sim$50 K, $\sigma$($\log{g}$)$\sim$0.2 dex, and $\sigma$($\xi$)$\sim$0.2 km/s) and in the abundance of both iron ($\lesssim$ 0.05 dex) and alpha elements ($\lesssim$0.10 dex) over the entire pulsation cycle. We also provide new and accurate effective temperature templates by splitting the calibrating Cepheids into four different period bins, ranging from short to long periods. For each period bin, we performed an analytical fit with Fourier series providing $\theta = 5040/{T_{\rm eff}}$ as a function of the pulsation phase. The current findings are a good viaticum to trace the chemical enrichment of the Galactic thin disk by using classical Cepheids and a fundamental stepping stone for further investigations into the more metal-poor regime typical of Magellanic Cepheids.

Two close-in planets have been recently found around the M-dwarf flare star AU Microscopii (AU Mic). These Neptune-sized planets (AU Mic b and c) seem to be located very close to the so-called "evaporation valley" in the exoplanet population, making this system an important target for studying atmospheric loss on exoplanets. This process, while mainly driven by the high-energy stellar radiation, will be strongly mediated by the space environment surrounding the planets. Here we present an investigation on this last area, performing 3D numerical modeling of the quiescent stellar wind from AU Mic, as well as time-dependent simulations describing the evolution of a highly energetic Coronal Mass Ejection (CME) event in this system. Observational constraints on the stellar magnetic field and properties of the eruption are incorporated in our models. We carry out qualitative and quantitative characterizations of the stellar wind, the emerging CMEs, as well as the expected steady and transient conditions along the orbit of both exoplanets. Our results predict an extreme space weather for AU Mic and its planets. This includes sub-Alfv\'enic regions for the large majority of the exoplanet orbits, very high dynamic and magnetic pressure values in quiescence (varying within $10^{2} - 10^{5}$ times the dynamic pressure experienced by the Earth), and an even harsher environment during the passage of any escaping CME associated with the frequent flaring observed in AU Mic. These space weather conditions alone pose an immense challenge for the survival of the exoplanetary atmospheres (if any) in this system.

Longitudinal oscillations in prominences are common phenomena on the Sun. These oscillations can be used to infer the geometry and intensity of the filament magnetic field. Previous theoretical studies of longitudinal oscillations made two simplifying assumptions: uniform gravity and semi-circular dips on the supporting flux tubes. However, the gravity is not uniform and realistic dips are not semi-circular. To understand the effects of including the nonuniform solar gravity on longitudinal oscillations, and explore the validity of the pendulum model with different flux-tube geometries. We first derive the equation describing the motion of the plasma along the flux tube including the effects of nonuniform gravity, yielding corrections to the original pendulum model. We also compute the full numerical solutions for the normal modes, and compare them with the new pendulum approximation. We have found that the nonuniform gravity introduces a significant modification in the pendulum model. We have also found a cut-off period, i.e. the longitudinal oscillations cannot have a period longer than 167 minutes. In addition, considering different tube geometries, the period depends almost exclusively on the radius of curvature at the bottom of the dip. We conclude that nonuniform gravity significantly modifies the pendulum model. These corrections are important for prominence seismology, because the inferred values of the radius of curvature and minimum magnetic-field strength differ substantially from those of the old model. However, we find that the corrected pendulum model is quite robust and is still valid for non-circular dips.

The impact of the magnetic field on the postbounce supernova dynamics of non-rotating stellar cores is studied by performing three-dimensional magnetohydrodynamics simulations with spectral neutrino transport. We compare the explodability between initially strongly and weakly magnetized models of $20$ and $27$ $M_{\odot}$ pre-supernova progenitors. We find that although the efficiency for the conversion of the neutrino heating into the turbulent energy including magnetic fields in the gain region is not significantly different between the strong and weak field models, the amplified magnetic field due to the neutrino-driven convection on large hot bubbles just behind a stalled shock leads to the faster and more energetic explosion in the strongly magnetized models. In addition, by comparing the difference between 2nd- and 5th-order spatial accuracy of the simulation in the strong field model for $27$ $M_{\odot}$ progenitor, we also find that the higher-order accuracy in space is positive for the explosion because it enhances the growth of neutrino-driven convection in the gain region. A new possibility of the origin of the magnetic field of the proto-neutron star (PNS) is proposed based on our results of core-collapse supernova simulations for the non-rotating model. The magnetic field is accumulated and amplified to the magnetar level, that is, $\mathcal{O}(10^{14})$ G, in the convectively stable shell in the vicinity of the PNS surface.

Fast radio bursts (FRBs) are extragalactic radio transients of extraordinary luminosity. Studying the diverse temporal and spectral behaviour recently observed in a number of FRBs may help determine the nature of the entire class. For example, a fast spinning or highly magnetised neutron star might generate the rotation-powered acceleration required to explain the bright emission. Periodic, sub-second components, suggesting such rotation, were recently reported in one FRB, and potentially in two more. Here we report the discovery of FRB 20201020A with Apertif, an FRB showing five components regularly spaced by 0.415 ms. This sub-millisecond structure in FRB 20201020A carries important clues about the progenitor of this FRB specifically, and potentially about that of FRBs in general. We thus contrast its features to the predictions of the main FRB source models. We perform a timing analysis of the FRB 20201020A components to determine the significance of the periodicity. We compare these against the timing properties of the previously reported CHIME FRBs with sub-second quasi-periodic components, and against two Apertif bursts from repeating FRB 20180916B that show complex time-frequency structure. We find the periodicity of FRB 20201020A to be marginally significant at 2.5$\sigma$. Its repeating subcomponents cannot be explained as a pulsar rotation since the required spin rate of over 2 kHz exceeds the limits set by typical neutron star equations of state and observations. The fast periodicity is also in conflict with a compact object merger scenario. These quasi-periodic components could, however, be caused by equidistant emitting regions in the magnetosphere of a magnetar. The sub-millisecond spacing of the components in FRB 20201020A, the smallest observed so far in a one-off FRB, may rule out both neutron-star rotation and binary mergers as the direct source of quasi-periodic FRBs.

Multiple stellar populations (MPs) representing star-to-star light-element abundance variations are common in nearly all ancient Galactic globular clusters. Here we provide the strongest evidence yet that the populous, ~ 1.7 Gyr-old Large Magellanic Cloud cluster NGC 2173 also exhibits light-element abundance variations. Thus, our results suggest that NGC 2173 is the youngest cluster for which MPs have been confirmed to date. Our conclusion is based on the distinct bifurcation at the tip of its red-giant branch in high-quality color--magnitude diagrams generated from Hubble Space Telescope imaging observations. Our results are further supported by a detailed analysis of 'pseudo-$UBI$' maps, which reveal clear evidence of a bimodality in the cluster's red-giant-branch color distribution. Young clusters in the Magellanic Clouds can provide critical insights into galaxy evolution histories. Our discovery of MPs in NGC 2173 suggests that ancient Galactic globular clusters and young massive clusters might share a common formation process.

We present a new version of the FIT3D and Pipe3D codes, two packages to derive properties of the stellar populations and the ionized emission lines from optical spectroscopy and integral field spectroscopy data respectively. The new codes have been fully transcribed to Python from the original Perl and C versions, modifying the algorithms when needed to make use of the unique capabilities of this language with the main goals of (1) respecting as much as possible the original philosophy of the algorithms, (2) maintaining a full compatibility with the original version in terms of the format of the required input and produced output files, and (3) improving the efficiency and accuracy of the algorithms, and solving known (and newly discovered) bugs. The complete package is freely distributed, with an available repository online. pyFIT3D and pyPipe3D are fully tested with data of the most recent IFS data surveys and compilations (e.g. CALIFA, MaNGA, SAMI and AMUSING++), and confronted with simulations. We describe here the code, its new implementation, its accuracy in recovering the parameters based on simulations, and a showcase of its implementation on a particular dataset.

We study the small-scale asymptotic behaviour of the cold dark matter density fluctuation power spectrum in the Zel'dovich approximation, without introducing an ultraviolet cut-off. Assuming an initially correlated Gaussian random field and spectral index $0 < n_s < 1$, we derive the small-scale asymptotic behaviour of the initial momentum-momentum correlations. This result is then used to derive the asymptotics of the power spectrum in the Zel'dovich approximation. Our main result is an asymptotic series, dominated by a $k^{-3}$ tail at large wave-numbers, containing higher-order terms that differ by integer powers of $k^{n_s-1}$ and logarithms of $k$. Furthermore, we show that dark matter power spectra with an ultraviolet cut-off develop an intermediate range of scales, where the power spectrum is accurately described by the asymptotics of dark matter without a cut-off. These results reveal information about the mathematical structure that underlies the perturbative terms in kinetic field theory and thus the non-linear power spectrum. We also discuss the sensitivity of the small-scale asymptotics to the spectral index $n_s$.

The properties of young massive clusters (YMCs) are key to understanding the star formation mechanism in starburst systems, especially mergers. We present ALMA high-resolution ($\sim$10 pc) continuum (100 and 345 GHz) data of YMCs in the overlap region of the Antennae galaxy. We identify 6 sources in the overlap region, including two sources that lie in the same giant molecular cloud (GMC). These YMCs correspond well with radio sources in lower resolution continuum (100 and 220 GHz) images at GMC scales ($\sim$60 pc). We find most of these YMCs are bound clusters through virial analysis. We estimate their ages to be $\sim$1 Myr and to be either embedded or just beginning to emerge from their parent cloud. We also compare each radio source with Pa$\beta$ source and find they have consistent total ionizing photon numbers, which indicates they are tracing the same physical source. By comparing the free-free emission at $\sim$10 pc scale and $\sim$60 pc scale, we find that $\sim$50% of the free-free emission in GMCs actually comes from these YMCs. This indicates that roughly half of the stars in massive GMCs are formed in bound clusters. We further explore the mass correlation between YMCs and GMCs in the Antennae and find it generally agrees with the predictions of the star cluster simulations. The most massive YMC has a stellar mass that is 1% - 5% of its host GMC mass.

The spatial range for feedback from star formation varies from molecular cloud disruption on parsec scales to supershells and disk blowout on kiloparsec scales. The relative amounts of energy and momentum given to these scales is important for understanding the termination of star formation in any one region and the origin of interstellar turbulence and disk stability in galaxies as a whole. Here we measure for eleven THINGS galaxies the excess kinetic energy, velocity dispersion and surface density of HI gas associated with regions of excess star formation, where the excess is determined from the difference between the observed local value and the azimuthal average. We find small decreases in the excess kinetic energy and velocity dispersion in regions of excess star formation rate density, suggesting that most of the feedback energy does not go into local HI motion. Most likely it disrupts molecular clouds and dissipates rapidly at high gas density. Some could also be distributed over larger regions, filling in spaces between the peaks of star formation and contributing to other energy sources from self-gravity and spiral arm shocks.

We discuss a problem about magnetic collapse as a possible process for singularity formation of the magnetic field in a finite time within ideal magneto-hydrodynamics for incompressible fluids. This process is very important from the point of view of various astrophysical applications, in particular, as a mechanism of magnetic filaments formation in the convective zone of the Sun. The collapse possibility is connected with compressibility of continuously distributed magnetic field lines. A well-known example of the formation of magnetic filaments in the kinematic dynamo approximation with a given velocity field, first considered by Parker in 1963, rather indicates that the increase in the magnetic field is exponential in time. In the case of the kinematic approximation for the induction equation, the magnetic filaments formation is shown to occur in areas with a hyperbolic velocity profile.

Radioactive nuclei are the key to understanding the circumstances of the birth of our Sun because meteoritic analysis has proven that many of them were present at that time. Their origin, however, has been so far elusive. The ERC-CoG-2016 RADIOSTAR project is dedicated to investigating the production of radioactive nuclei by nuclear reactions inside stars, their evolution in the Milky Way Galaxy, and their presence in molecular clouds. So far, we have discovered that: (i) radioactive nuclei produced by slow ($^{107}$Pd and $^{182}$Hf) and rapid ($^{129}$I and $^{247}$Cm) neutron captures originated from stellar sources - asymptotic giant branch (AGB) stars and compact binary mergers, respectively - within the galactic environment that predated the formation of the molecular cloud where the Sun was born; (ii) the time that elapsed from the birth of the cloud to the birth of the Sun was of the order of 10$^7$ years, and (iii) the abundances of the very short-lived nuclei $^{26}$Al, $^{36}$Cl, and $^{41}$Ca can be explained by massive star winds in single or binary systems, if these winds directly polluted the early Solar System. Our current and future work, as required to finalise the picture of the origin of radioactive nuclei in the Solar System, involves studying the possible origin of radioactive nuclei in the early Solar System from core-collapse supernovae, investigating the production of $^{107}$Pd in massive star winds, modelling the transport and mixing of radioactive nuclei in the galactic and molecular cloud medium, and calculating the galactic chemical evolution of $^{53}$Mn and $^{60}$Fe and of the p-process isotopes $^{92}$Nb and $^{146}$Sm.

We present the observations of three sympathetic filament eruptions occurring on 19 July 2015 namely F1, F2, and F3. The events were observed in UV/EUV wavelengths by Atmospheric Imaging Assembly onboard the Solar Dynamics Observatory and by Global Oscillation Network Group telescope in H{\alpha} line. As filament F1 starts to erupt, a part of it falls close to the location of the F2 and F3 filaments. This causes the eruption of F2 and F3 during which the two filaments merge together and trigger a medium-class CME and a long-duration GOES C2.1 class flare. We discuss the dynamics and kinematics of these three filament eruptions and related phenomena.

The Cepheid variables in SMC, LMC, the Milky Way, M33 and M31 are used to examine the dependence of pulsation mode on metallicity which was previously found in red supergiants. The initial samples of Cepheids are collected from the Cepheid catalogs identified from the OGLE, PS1, DIRECT, WISE and ZTF surveys. The contaminants are removed with the help of the Gaia/EDR3 astrometric information for extra galaxies or by comparing the geometric distance and the distance from the P-L relation for the Milky Way. The division of fundamental and first-overtone mode is refined according to the gap between the two modes in the P-L diagram of the objects in each galaxy. The ratio of FU/(FU+1O) is found to be 0.59, 0.60, 0.69, 0.83 and 0.85 for SMC, LMC, the Milky Way, M33 and M31 respectively in order of metallicity, which confirms that the pulsation mode depends on metallicity in the way that the ratio of FU/(FU+1O) increases with metallicity. This dependence is not changed if the incompleteness of the samples is taken into account.

The taxonomy of galaxy morphology is critical in astrophysics as the morphological properties are powerful tracers of galaxy evolution. With the upcoming Large-scale Imaging Surveys, billions of galaxy images challenge astronomers to accomplish the classification task by applying traditional methods or human inspection. Consequently, machine learning, in particular supervised deep learning, has been widely employed to classify galaxy morphologies recently due to its exceptional automation, efficiency, and accuracy. However, supervised deep learning requires extensive training sets, which causes considerable workloads; also, the results are strongly dependent on the characteristics of training sets, which leads to biased outcomes potentially. In this study, we attempt Few-shot Learning to bypass the two issues. Our research adopts the dataset from Galaxy Zoo Challenge Project on Kaggle, and we divide it into five categories according to the corresponding truth table. By classifying the above dataset utilizing few-shot learning based on Siamese Networks and supervised deep learning based on AlexNet, VGG_16, and ResNet_50 trained with different volumes of training sets separately, we find that few-shot learning achieves the highest accuracy in most cases, and the most significant improvement is $21\%$ compared to AlexNet when the training sets contain 1000 images. In addition, to guarantee the accuracy is no less than 90\%, few-shot learning needs $\sim$6300 images for training, while ResNet_50 requires 13000 images. Considering the advantages stated above, foreseeably, few-shot learning is suitable for the taxonomy of galaxy morphology and even for identifying rare astrophysical objects, despite limited training sets consisting of observational data only.

In this paper we present the deep, wide-field and multi-band imaging of the LEO I pair NGC~3379-NGC~3384, from the VST Early-type GAlaxy Survey (VEGAS). The main goal of this study is to map the intra-group baryons in the pair, in the form of diffuse light and globular clusters (GCs). Taking advantage from the large covered area, which extends for $\sim$ 3.9 square degrees around the pair, and the long integration time, we can map the light distribution out to $\sim$ 63 kpc and down to $\sim$ 30 mag/arcsec$^2$ in the $g$ band and $\sim$ 29 mag/arcsec$^2$ in the $r$ band, deeper than previous data available for this target. The map of the intra-group light (IGL) presents two very faint ($\mu_g \sim$ 28-29 mag/arcsec$^2$) streams protruding from the brightest group member NGC~3379 and elongated toward North-West and South. We estimate that the fraction of the stellar halo around NGC~3379 plus the IGL is $\sim 17 \pm 2\%$ in both $g$ and $r$ bands, with an average color $g$-$r$= $0.75 \pm 0.04$~mag. The color distribution of the GCs appears multi-modal, with two dominant peaks at (u-r) = 1.8 mag and (u-r) = 2.1 mag, respectively. The GC population stretches from North-East to South-West and from North-West to South of the pair, in the last case overlapping with the streams of IGL, as well as the PNe distribution found by \citet{Hartke2020} and \citet{Hartke2022arXiv220108710H}. Since these structures are elongated in the direction of the two nearby galaxies M96 and NGC~3338, they could be the remnant of a past gravitational interactions with the pair.

Recurrent Novae result from thermonuclear explosions in the outer layers of White Dwarfs, as they accrete from their Red Giant companions. Ejected material drives an expanding shock into the companion star's wind, accelerating particles to relativistic energies. We report the H.E.S.S. detection of very-high-energy gamma rays from the recurrent Nova RS Ophiuchi up to a month after the 2021 outburst. A common origin of the H.E.S.S. emission and the high-energy emission detected with Fermi-LAT is favoured, due to their similar decay profiles, $\propto t^{-1.7}$. The peak flux in very-high-energies is delayed by two days with respect to Fermi-LAT. These observations reveal time-dependent particle energization, and provide a real-time window on an efficient cosmic accelerator. With this measurement, we establish recurrent Novae as multi-TeV Galactic transient sources.

Counts of galaxy clusters offer a high-precision probe of cosmology, but control of systematic errors will determine the accuracy of this measurement. Using Buzzard simulations, we quantify one such systematic, the triaxiality distribution of clusters identified with the redMaPPer optical cluster finding algorithm, which was used in the Dark Energy Survey Year-1 (DES Y1) cluster cosmology analysis. We test whether redMaPPer selection biases the clusters' shape and orientation and find that it only biases orientation, preferentially selecting clusters with their major axes oriented along the line of sight. Modeling the richness-mass relation as a log-linear relation, we find that the log-richness amplitude $\ln(A)$ is boosted from the lowest to highest orientation bin with a significance of $14\sigma$, while the orientation dependence of the richness-mass slope and intrinsic scatter is minimal. We also find that the weak lensing shear-profile ratios of cluster-associated dark halos in different orientation bins resemble a "bottleneck" shape that can be quantified with a Cauchy function. We test the correlation of orientation with two other leading systematics in cluster cosmology -- miscentering and projection -- and find a null correlation. Analytic templates for the triaxiality bias of observed-richness and lensing profiles are mapped as corrections to the observable of richness-binned lensing profiles for redMaPPer clusters. The resulting mass bias confirms the DES Y1 finding that triaxiality is a leading source of bias in cluster cosmology. However, the richness-dependence of the bias confirms that triaxiality does not fully resolve the tension at low-richness between DES Y1 cluster cosmology and other probes. Our model can be used for quantifying the impact of triaxiality bias on cosmological constraints for upcoming weak lensing surveys of galaxy clusters.

How do we appropriately fit a model based on an idealised Friedmann-Lema{\i}tre Robertson-Walker (FLRW) spacetime to observations made from a single location in a lumpy Universe? We address this question for surveys that measure the imprints of the Baryon Acoustic Oscillation (BAO) in galaxy distribution and the peak apparent magnitude of the type 1A supernova (SN1A). These observables are related to the cosmological model through the Alcock-Paczynski parameters and the distance-redshift relation. Using the corresponding inhomogeneous spacetime expressions of these as observed data, we perform a parameter inference assuming that the background FLRW model is the correct model of the universe. This process allows us to estimate the best fit Hubble rate and the deceleration parameter. We find that the inferred Hubble rate from the monopole of the Alcock-Paczynski parameters is in tension with the Hubble rate determined using the distance-redshift relation. The latter gives the best fit Hubble rate for the cosmological expansion. The constraint on the Hubble rate from the Alcock-Paczynski parameters is contaminated by the environment. When the environmental contribution is restricted to modes in the Hubble flow, we find about (9-12)\% discrepancy in the Hubble rate. Finally, we comment on the insufficiency of the method of cosmography in constraining the deceleration parameter.

Dwarf galaxies are by far the most numerous galaxies in the Universe, showing properties that are quite different from those of their larger and more luminous cousins. This review focuses on the physical and chemical properties of the interstellar medium of those dwarfs that are known to host significant amounts of gas and dust. The neutral and ionized gas components and the impact of the dust will be discussed, as well as first indications for the existence of active nuclei in these sources. Cosmological implications are also addressed, considering the primordial helium abundance and the similarity of local Green Pea galaxies with young, sometimes proto-galactic sources in the early Universe.

Recent cosmological analyses rely on the ability to accurately sample from high-dimensional posterior distributions. A variety of algorithms have been applied in the field, but justification of the particular sampler choice and settings is often lacking. Here we investigate three such samplers to motivate and validate the algorithm and settings used for the Dark Energy Survey (DES) analyses of the first 3 years (Y3) of data from combined measurements of weak lensing and galaxy clustering. We employ the full DES Year 1 likelihood alongside a much faster approximate likelihood, which enables us to assess the outcomes from each sampler choice and demonstrate the robustness of our full results. We find that the ellipsoidal nested sampling algorithm $\texttt{MultiNest}$ reports inconsistent estimates of the Bayesian evidence and somewhat narrower parameter credible intervals than the sliced nested sampling implemented in $\texttt{PolyChord}$. We compare the findings from $\texttt{MultiNest}$ and $\texttt{PolyChord}$ with parameter inference from the Metropolis-Hastings algorithm, finding good agreement. We determine that $\texttt{PolyChord}$ provides a good balance of speed and robustness, and recommend different settings for testing purposes and final chains for analyses with DES Y3 data. Our methodology can readily be reproduced to obtain suitable sampler settings for future surveys.

How do we appropriately fit a model based on an idealised Friedmann-Lema\^itre Robertson-Walker (FLRW) spacetime to observations made from one location in a lumpy Universe? We address this question for the apparent magnitude of Type 1A supernova (SN1A) and the acoustic peaks in the density field. We show that distance in a perturbed universe is given by the FLRW spacetime on average, except in the immediate neighbourhood of an observer. We demonstrate how this distinction impacts the determination of the Hubble rate and the SN1A absolute magnitude.

We explore the use of the mass spectrum of neutron stars and black holes in gravitational-wave compact binary sources as a cosmological probe. These standard siren sources provide direct measurements of luminosity distance. In addition, features in the mass distribution, such as mass gaps or peaks, will redshift, and thus provide independent constraints on their redshift distribution. We argue that the mass spectrum of LIGO/Virgo/KAGRA events introduces at least five independent mass "features": the upper and lower edges of the pair instability supernova (PISN) gap, the upper and lower edges of the neutron star-black hole gap, and the minimum neutron star mass. We find that the PISN gap dominates the cosmological inference with current detectors (2G), as shown in previous work. We further argue that the lower mass gap will provide the most powerful constraints in the era of Cosmic Explorer and Einstein Telescope (3G). We demonstrate that degeneracies between redshift evolution of the source masses and cosmology can be broken, unless an astrophysical conspiracy shifts all features of the full mass distribution simultaneously following the (non-trivial) Hubble diagram evolution. We find that this "spectral siren" method has the potential to constrain both cosmology and the evolution of the mass distribution, with 2G achieving better than $10\%$ precision on $H(z)$ at $z\lesssim1$ within a year, and 3G reaching $\lesssim1\%$ at $z\gtrsim2$ within one month.

Polycyclic aromatic hydrocarbons (PAHs) play a key role in the chemical and hydrodynamical evolution of the atmospheres of exoplanets and planet-forming discs. If they can survive the planet formation process, PAHs are likely to be involved in pre-biotic chemical reactions eventually leading to more complex molecules such as amino acids and nucleotides, which form the basis for life as we know it. However, the abundance and specific role of PAHs in these environments is largely unknown due to limitations in sensitivity and range of wavelength of current and previous space-borne facilities. Upcoming infrared space spectroscopy missions, such as Twinkle and Ariel, present a unique opportunity to detect PAHs in the atmospheres of exoplanets and planet-forming discs. In this work we present synthetic observations based on conservative numerical modeling of typical planet-forming discs and a transiting hot Saturnian planet around solar type star. Our models show that Twinkle and Ariel might both be able to detect the 3.3 micron PAH feature within reasonable observing time in discs and transiting planets, assuming that PAHs are present with an abundance of at least one tenth of the interstellar medium value.

Artificial skyglow is a form of light pollution with wide ranging implications on the environment. The extent, intensity and color of skyglow depends on the artificial light sources and weather conditions. Skyglow can be best determined with ground based instruments. We mapped the skyglow of Berlin, Germany, for clear sky and overcast sky conditions inside and outside of the city limits. We conducted observations using a transect from the city center of Berlin towards a rural place more than 58 km south of Berlin using all-sky photometry with a calibrated commercial digital camera and a fisheye lens. From the multispectral imaging data, we processed luminance and correlated color temperature maps. We extracted the night sky brightness and correlated color temperature at zenith, as well as horizontal and scalar illuminance simultaneously. We calculated cloud amplification factors at each site and investigated the changes of brightness and color with distance, particularly showing differences inside and outside of the city limits. We found high values for illuminance above full moon light levels and amplification factors as high as 25 in the city center and a gradient towards the city limit and outside of the city limit. We further observed that clouds decrease the correlated color temperature in almost all cases. We discuss advantages and weaknesses of our method, compare the results with modeled night sky brightness data and provide recommendations for future work.

At low Reynolds numbers, the wind flow in the wake of a single wind turbine is generally not turbulent. However, turbines in wind farms affect each other's wakes so that a turbulent flow can arise. In the present work, an analogue of this effect for the massless charge carrier flow around obstacles in graphene is outlined. We use a relativistic hydrodynamic simulation to analyze the flow in a sample containing impurities. Depending on the density of impurities in the sample, we indeed find evidence for potentially turbulent flow and discuss experimental consequences.

The pyramids of the Giza plateau have fascinated visitors since ancient times and are the last of the Seven Wonders of the ancient world still standing. It has been half a century since Luiz Alvarez and his team used cosmic-ray muon imaging to look for hidden chambers in Khafres Pyramid. Advances in instrumentation for High-Energy Physics (HEP) allowed a new survey, ScanPyramids, to make important new discoveries at the Great Pyramid (Khufu) utilizing the same basic technique that the Alvarez team used, but now with modern instrumentation. The Exploring the Great Pyramid Mission plans to field a very-large muon telescope system that will be transformational with respect to the field of cosmic-ray muon imaging. We plan to field a telescope system that has upwards of 100 times the sensitivity of the equipment that has recently been used at the Great Pyramid, will image muons from nearly all angles and will, for the first time, produce a true tomographic image of such a large structure.

In the current study, we investigated a specific model of anisotropic strange stars specially Her X-1, in the background of modified f(R,T) gravity by choosing f(R,T) = R+2{\xi}T, where R is Ricci scalar, T is the trace of the energy-momentum tensor and {\xi} is a coupling constant. To obtained the solution for the modified field equations, we apply Buchdahl metric to our equations. We consider the case, when the matter is governed by MIT bag model equation of state as Pr =1/3({\rho} -4B), where B is bag constant. We calculate the values of unknown parameters using Schwarzschild interior space-time followed by choosing the appropriate values of parameter {\xi} , K, {\beta} and also tabulated different values of bag constant for different {\xi} . We examine the physical validity of our model by performing tests such as energy conditions, equilibrium of the forces, the adiabatic index, redshift and some more. The observational results showed that the proposed f(R,T) model satisfies all these tests and are quite acceptable.