Papers for Thursday, May 27 2021

The detection of dark matter subhalos without a stellar component in the Galactic halo remains a challenge. We use supervised machine learning to identify high-latitude gamma-ray sources with dark matter-like spectra among unassociated gamma-ray sources in the 4FGL-DR2. Out of 843 4FGL-DR2 unassociated sources at $|b| \geq 10\mathrm{^\circ}$, we select 73 dark matter subhalo candidates. Of the 69 covered by the Neil Gehrels Swift Observatory (Swift), 17 show at least one X-ray source within the 95% LAT error ellipse and 52 where we identify no new sources. This latest inventory of dark subhalos candidates allows us to investigate the possible dark matter substructure responsible for the perturbation in the GD-1 stellar stream. In particular, we examine the possibility that the alleged GD-1 dark subhalo may appear as a 4FGL-DR2 gamma-ray source from dark matter annihilation into Standard Model particles.

In star-forming galaxies, the far-infrared (FIR) and radio-continuum luminosities obey a tight empirical relation over a large range of star-formation rates (SFR). To understand the physics, we examine magneto-hydrodynamic galaxy simulations, which follow the genesis of cosmic ray (CR) protons at supernovae and their advective and anisotropic diffusive transport. We show that gravitational collapse of the proto-galaxy drives a small-scale dynamo, which exponentially amplifies weak seed magnetic fields. After saturation at small scales, they grow in scale to reach equipartition with thermal and CR energy densities in Milky Way-mass galaxies. In small galaxies, the magnetic energy saturates at the turbulent energy while it fails to reach equipartition with thermal and CR energies. We solve for steady-state spectra of CR protons, secondary electrons/positrons from hadronic CR-proton interactions with the interstellar medium, and primary shock-accelerated electrons at supernovae. The modeled radio-synchrotron emission is dominated by primary electrons, irradiates the magnetised disc and bulge of our simulated Milky Way-mass galaxy and weakly traces bubble-shaped magnetically-loaded outflows. Our star-forming and star-bursting galaxies with saturated magnetic fields match the global FIR-radio correlation (FRC) across four orders of magnitude. Its intrinsic scatter arises due to (i) different magnetic saturation levels that result from different seed magnetic fields, (ii) different radio synchrotron luminosities for different specific SFRs at fixed SFR and (iii) a varying radio intensity with galactic inclination. In agreement with observations, several 100-pc-sized regions within star-forming galaxies also obey the FRC, while the centres of starbursts exceed the FRC by a substantial amount.

We make use of ALMA continuum observations of $15$ luminous Lyman-break galaxies at $z$$\sim$$7$$-$$8$ to probe their dust-obscured star-formation. These observations are sensitive enough to probe to obscured SFRs of $20$ $M_{\odot}$$/$$yr$ ($3\sigma$). Six of the targeted galaxies show significant ($\geq$$3$$\sigma$) dust continuum detections, more than doubling the number of known dust-detected galaxies at $z$$>$$6.5$. Their IR luminosities range from $2.7$$\times$$10^{11}$ $L_{\odot}$ to $1.1$$\times$$10^{12}$ $L_{\odot}$, equivalent to obscured SFRs of $20$ to $105$ $M_{\odot}$$/$$yr$. We use our results to quantify the correlation of the infrared excess IRX on the UV-continuum slope $\beta_{UV}$ and stellar mass. Our results are most consistent with an SMC attenuation curve for intrinsic $UV$-slopes $\beta_{UV,intr}$ of $-2.63$ and most consistent with an attenuation curve in-between SMC and Calzetti for $\beta_{UV,intr}$ slopes of $-2.23$, assuming a dust temperature $T_d$ of $50$ K. Our fiducial IRX-stellar mass results at $z$$\sim$$7$$-$$8$ are consistent with marginal evolution from $z$$\sim$$0$. We then show how both results depend on $T_d$. For our six dust-detected sources, we estimate their dust masses and find that they are consistent with dust production from SNe if the dust destruction is low ($<$$90$%). Finally we determine the contribution of dust-obscured star formation to the star formation rate density for $UV$ luminous ($<$$-$$21.5$ mag: $\gtrsim$$1.7$$L_{UV} ^*$) $z$$\sim$$7$$-$$8$ galaxies, finding that the total SFR density at $z$$\sim$$7$ and $z$$\sim$$8$ from bright galaxies is $0.18_{-0.10}^{+0.08}$ dex and $0.20_{-0.09}^{+0.05}$ dex higher, respectively, i.e. $\sim$$\frac{1}{3}$ of the star formation in $\gtrsim$$1.7$$L_{UV} ^*$ galaxies at $z$$\sim$$7$$-$$8$ is obscured by dust.

An extinction-free estimator of the star-formation rate (SFR) of galaxies is critical for understanding the high-redshift universe. To this end, the nearly linear, tight correlation of far-infrared (FIR) and radio luminosity of star-forming galaxies is widely used. While the FIR is linked to massive star formation, which also generates shock-accelerated cosmic ray (CR) electrons and radio synchrotron emission, a detailed understanding of the underlying physics is still lacking. Hence, we perform three-dimensional magneto-hydrodynamical (MHD) simulations of isolated galaxies over a broad range of halo masses and SFRs using the moving-mesh code AREPO, and evolve the CR proton energy density self-consistently. In post-processing, we calculate the steady-state spectra of primary, shock-accelerated and secondary CR electrons, which result from hadronic CR proton interactions with the interstellar medium. The resulting total radio luminosities correlate with the FIR luminosities as observed and are dominated by primary CR electrons if we account for anisotropic CR diffusion. The increasing contribution of secondary emission up to 30 per cent in starbursts is compensated by the larger bremsstrahlung and Coulomb losses. CR electrons are in the calorimetric limit and lose most of their energy through inverse Compton interactions with star-light and cosmic microwave background (CMB) photons while less energy is converted to synchrotron emission. This implies steep steady-state synchrotron spectra in starbursts. Interestingly, we find that thermal free-free emission hardens the total radio spectra at high radio frequencies and reconciles calorimetric theory with observations while free-free absorption explains the observed low-frequency flattening towards the central regions of starbursts.

[abridged] We analyse four transits of WASP-33b observed with the optical high-resolution HARPS-N spectrograph to confirm its nodal precession, study its atmosphere and investigate the presence of star-planet interactions.We extract the mean line profiles of the spectra by using the LSD method, and analyse the Doppler shadow and the RVs. We also derive the transmission spectrum of the planet, correcting it for the stellar contamination due to rotation, CLV and pulsations. We confirm the previously discovered nodal precession of WASP-33b, almost doubling the time coverage of the inclination and projected spin-orbit angle variation. We find that the projected obliquity reached a minimum in 2011 and use this constraint to derive the geometry of the system, in particular its obliquity at that epoch ($\epsilon=113.99^{\circ}\pm 0.22^{\circ}$) and the inclination of the stellar spin axis ($i_{\rm s}=90.11^{\circ}\pm 0.12^{\circ}$), as well as the gravitational quadrupole moment of the star $J_2=(6.73\pm 0.22)\times 10^{-5}$. We present detections of H$\alpha$ and H$\beta$ absorption in the atmosphere of the planet with a contrast almost twice smaller than previously detected in the literature. We also find evidence for the presence of a pre-transit signal, which repeats in all four analysed transits. The most likely explanation lies in a possible excitation of a stellar pulsation mode by the presence of the planetary companion. Future common analysis of all available datasets in the literature will help shedding light on the possibility that the observed Balmer lines transit depth variations are related to stellar activity and/or pulsation, and to set constraints on the energetics possibly driving atmospheric escape. A complete orbital phase coverage of WASP-33b with high-resolution spectroscopic (spectro-polarimetric) observations could help understanding the nature of the pre-transit signal.

Evolved Wolf-Rayet stars form a key aspect of massive star evolution, and their strong outflows determine their final fates. In this study, we calculate grids of stellar models for a wide range of initial masses at five metallicities (ranging from solar down to just 2% solar). We compare a recent hydrodynamically-consistent wind prescription with two earlier frequently-used wind recipes in stellar evolution and population synthesis modelling, and we present the ranges of maximum final masses at core He-exhaustion for each wind prescription and metallicity Z. Our model grids reveal qualitative differences in mass-loss behaviour of the wind prescriptions in terms of "convergence". Using the prescription from Nugis & Lamers the maximum stellar black hole is found to converge to a value of 20-30Msun, independent of host metallicity, however when utilising the new physically-motivated prescription from Sander & Vink there is no convergence to a maximum black hole mass value. The final mass is simply larger for larger initial He-star mass, which implies that the upper black hole limit for He-stars below the pair-instability gap is set by prior evolution with mass loss, or the pair instability itself. Quantitatively, we find the critical Z for pair-instability (Z_PI) to be as high as 50% Zsolar, corresponding to the host metallicity of the LMC. Moreover, while the Nugis & Lamers prescription would not predict any black holes above the approx 130Msun pair-instability limit, with Sander & Vink winds included, we demonstrate a potential channel for very massive helium stars to form such massive black holes at ~2% Zsolar or below.

Identifying stars found in the Milky Way as having formed in situ or accreted can be a complex and uncertain undertaking. We use Gaia kinematics and APOGEE elemental abundances to select stars belonging to the Gaia-Sausage-Enceladus (GSE) and Sequoia accretion events. These samples are used to characterize the GSE and Sequoia population metallicity distribution functions, elemental abundance patterns, age distributions, and progenitor masses. We find that the GSE population has a mean [Fe/H] $\sim -1.15$ and a mean age of 10-12 Gyr. GSE has a single sequence in [Mg/Fe] vs [Fe/H] consistent with the onset of SN Ia Fe contributions and uniformly low [Al/Fe] of $\sim -0.25$ dex. The derived properties of the Sequoia population are strongly dependent on the kinematic selection. We argue the selection with the least contamination is $J_{\phi}/J_{\mbox{tot}} < -0.6$ and $(J_z - J_R)/J_{\mbox{tot}} < 0.1$. This results in a mean [Fe/H] $\sim -1.3$ and a mean age of 12-14 Gyr. The Sequoia population has a complex elemental abundance distribution with mainly high [Mg/Fe] stars. We use the GSE [Al/Fe] vs [Mg/H] abundance distribution to inform a chemically-based selection of accreted stars, which is used to remove possible contaminant stars from the GSE and Sequoia samples.

We use the Cosmic Assembly Deep Near-infrared Extragalactic Legacy Survey (CANDELS) data to study the relationship between quenching and the stellar mass surface density within the central radius of 1 kpc ($\Sigma_1$) of low-mass galaxies (stellar mass $M_* \lesssim 10^{9.5} M_\odot$) at $0.5 \leq z < 1.5$. Our sample is mass complete down to $\sim 10^9 M_\odot$ at $0.5 \leq z < 1.0$. We compare the mean $\Sigma_1$ of star-forming galaxies (SFGs) and quenched galaxies (QGs) at the same redshift and $M_*$. We find that low-mass QGs have higher $\Sigma_1$ than low-mass SFGs, similar to galaxies above $10^{10} M_\odot$. The difference of $\Sigma_1$ between QGs and SFGs increases slightly with $M_*$ at $M_* \lesssim 10^{10} M_\odot$ and decreases with $M_*$ at $M_* \gtrsim 10^{10} M_\odot$. The turnover mass is consistent with the mass where quenching mechanisms transition from internal to environmental quenching. At $0.5 \leq z < 1.0$, we find that the $\Sigma_1$ of galaxies increases by about 0.25 dex in the green valley (i.e., the transitioning region from star forming to fully quenched), regardless of their $M_*$. Using the observed specific star formation rate (sSFR) gradient in the literature as a constraint, we estimate that the quenching timescale (i.e., time spent in the transition) of low-mass galaxies is a few ($\sim4$) Gyrs at $0.5 \leq z < 1.0$. The mechanisms responsible for quenching need to gradually quench star formation in an outside-in way, i.e., preferentially ceasing star formation in outskirts of galaxies while maintaining their central star formation to increase $\Sigma_1$. An interesting and intriguing result is the similarity of the growth of $\Sigma_1$ in the green valley between low-mass and massive galaxies, which suggests that the role of internal processes in quenching low-mass galaxies is a question worthy of further investigation.

Recent observations reveal that, at a given stellar mass, blue galaxies tend to live in haloes with lower mass while red galaxies live in more massive host haloes. The physical driver behind this is still unclear because theoretical models predict that, at the same halo mass, galaxies with high stellar masses tend to live in early-formed haloes which naively leads to an opposite trend. Here, we show that the {\sc Simba} simulation quantitatively reproduces the colour bimodality in SHMR and reveals an inverse relationship between halo formation time and galaxy transition time. It suggests that the origin of this bimodality is rooted in the intrinsic variations of the cold gas content due to halo assembly bias. {\sc Simba}'s SHMR bimodality quantitatively relies on two aspects of its input physics: (1) Jet-mode AGN feedback, which quenches galaxies and sets the qualitative trend; and (2) X-ray AGN feedback, which fully quenches galaxies and yields better agreement with observations. The interplay between the growth of cold gas and the AGN quenching in {\sc Simba} results in the observed SHMR bimodality.

We forecast the ability of future-generation experiments to detect the fine-structure lines of the carbon and oxygen ions, [CII] and [OIII] in intensity mapping (IM) from the Epoch of Reionization ($z \sim 6-8$). Combining the latest empirically derived constraints relating the luminosity of the [OIII] line to the ambient star-formation rate, and using them in conjunction with previously derived estimates for the abundance of [CII] in haloes, we predict the expected auto-correlation IM signal to be observed using next-generation facilities based on the Fred Young Submillimetre Telescope (FYST) and the balloon-borne facility, Experiment for Cryogenic Large-Aperture Intensity Mapping (EXCLAIM) over $z \sim 5.3 - 7$. We describe how improvements to both the ground-based and balloon-based surveys in the future will enable a cross-correlation signal to be detected at $\sim$ 10-40 $\sigma$ over $z \sim 5.3 - 7$. Finally, we propose a space-based mission targeting the [OIII] 88 and 52 $\mu$m lines along with the [CII] 158 $\mu$m line, configured to enhance the signal-to-noise ratio of cross-correlation measurements. We find that such a configuration can achieve a high-significance detection (hundreds to thousands of $\sigma$) in both auto- and cross-correlation modes.

Super-Jupiters, brown dwarfs, and stars can form from the collapse of self-gravitating discs. Such discs are turbulent, with flocculent spiral arms accelerating gas to transonic speeds horizontally and vertically. Objects that fragment from gravito-turbulent discs should spin with a wide range of directions, reflecting the random orientations of their parent eddies. We show by direct numerical simulation that obliquities of newly collapsed fragments can range up to 45$^\circ$. Subsequent collisions between fragments can further broaden the obliquity distribution up to 90$^\circ$. The large obliquities of newly discovered super-Jupiters on wide orbits around young stars may thus be gravito-turbulent in origin. Obliquely spinning fragments are born on orbits that may be inclined relative to their parent discs by up to 20$^\circ$, and gravitationally stir leftover material to many times the pre-fragmentation disc thickness.

When two Neutron Stars collide a multi-band electromagnetic emission, known as Kilonova (KN), follows being powered by the radioactive decay of ejecta products. In this contribution we discuss how the equation of state of dense matter, impacts the mass and velocity in the KN ejecta and thus its light curve. Using this information encoded in the stellar mass-radius relationship, we ellaborate on how the future experimental observations in photon channels, in addition to complementary multimessenger probes, could provide a new and more detailed insight into the equation of state of nuclear matter.

Ices are an important constituent of protoplanetary disks. New observational facilities, notably JWST, will greatly enhance our view of disk ices by measuring their infrared spectral features. We present a suite of models to complement these upcoming observations. Our models use a kinetics-based gas-grain chemical evolution code to simulate the distribution of ices in a disk, followed by a radiative transfer code using a subset of key ice species to simulate the observations. We present models reflecting both molecular inheritance and chemical reset initial conditions. We find that H$_2$O, CO$_2$, and CH$_3$OH near-to-mid-IR absorption features are readily observable in disk-integrated spectra of highly-inclined disks while CO, NH$_3$, and CH$_4$ ice do not show prominent features. CH$_3$OH ice has low abundance and is not observable in the reset model, making this species an excellent diagnostic of initial chemical conditions. CO$_2$ ice features exhibit the greatest change over disk lifetime: decreasing and increasing for the inheritance and reset models, respectively. Spatially-resolved spectra of edge-on disks, possible with JWST's integral field unit observing modes, are ideal for constraining the vertical distribution of ices and may be able to isolate features from ices closer to the midplane (e.g., CO) given sufficient sensitivity. Spatially-resolved spectra of face-on disks can trace scattered-light features from H$_2$O, CO$_2$, and CH$_3$OH, plus CO and CH$_4$ from the outermost regions. We additionally simulate far-IR H$_2$O ice emission features and find they are strongest for disks viewed face-on.

We present a dynamical study of the narrow-line regions in two nearby QSO2s. We construct dynamical models based on detailed photoionization models of the emission-line gas, including the effects of internal dust, to apply to observations of large-scale outflows from these AGNs. We use Mrk 477 and Mrk 34 in order to test our models against recent HST STIS observations of [O III] emission-line kinematics since these AGNs possess more energetic outflows than found in Seyfert galaxies. We find that the outflows within 500 pc are consistent with radiative acceleration of dusty gas, however, the outflows in Mrk 34 are significantly more extended and may not be directly accelerated by radiation. We characterize the properties of X-ray winds found from the expansion of [O III]-emitting gas close to the black hole. We show that such winds possess the kinetic energy density to disturb [O III] gas at 1.8 kpc, and have sufficient energy to entrain the [O III] clouds at 1.2 kpc. Assuming that the X-ray wind possesses the same radial mass distribution as the [O III] gas, we find that the peak kinetic luminosity for this wind is 2% of Mrk 34's bolometric luminosity, which is in the 0.5% - 5% range required by some models for efficient feedback. Our work shows that, although the kinetic luminosity as measured from [O III]-emitting gas is frequently low, X-ray winds may provide more than one order of magnitude higher kinetic power.

We present a new measurement of the Ly$\alpha$ luminosity function at redshift $z=6.9$, finding moderate evolution from $z=5.7$ that is consistent with a fully or largely ionized $z\sim7$ intergalactic medium. Our result is based on four fields of the LAGER (Lyman Alpha Galaxies in the Epoch of Reionization) project. Our survey volume of $6.1\times10^{6}$ Mpc$^{3}$ is double that of the next largest $z\sim 7$ survey. We combine two new LAGER fields (WIDE12 and GAMA15A) with two previously reported LAGER fields (COSMOS and CDFS). In the new fields, we identify $N=95$ new $z=6.9$ Ly$\alpha$ emitters (LAEs); characterize our survey's completeness and reliability; and compute Ly$\alpha$ luminosity functions. The best-fit Schechter luminosity function parameters for all four LAGER fields are in good general agreement. Two fields (COSMOS and WIDE12) show evidence for a bright-end excess above the Schechter function fit. We find that the Ly$\alpha$ luminosity density declines at the same rate as the UV continuum LF from $z=5.7$ to $z=6.9$. This is consistent with an intergalactic medium that was fully ionized as early as redshift $z\sim 7$, or with a volume-averaged neutral hydrogen fraction of $x_{HI} < 0.33$ at $1\sigma$.

We aim to investigate the IR excess of 45 MW Cepheids combining different observables to constrain the presence of CSE. We used the SPIPS algorithm, a robust implementation of the parallax-of-pulsation method that combines photometry, angular diameter, stellar effective temperature, and RV measurements in a global modelling of the pulsation. We obtained new photometric measurements at mid-IR with the VLT/VISIR complemented with literature data. We then compared the mean magnitude from 0.5$\mu$m to 70$\mu$m with stellar atmosphere models to infer the IR excess, which we attribute to the presence of a CSE. We report that at least 29% of our sample have a detected IR excess. We estimated a mean excess of 0.08mag in K and 0.13mag in N. Other Cepheids possibly have IR excess, but they were rejected due to their low detection level compared to a single-star model. We do not see any correlation between the IR excess and the pulsation period as previously suspected for MW Cepheids, but a rather constant trend. We also do not find any correlation between the CO absorption and the presence of a CSE, but rather with the stellar effective temperature, which confirms that the CO features previously reported are mostly photospheric. No bias caused by the presence of the CSE is detected on the average distance estimates from a SPIPS analysis with a fitted colour excess. We also do not find correlation between the presence of IR excess and the evolution stage. We report a fraction of 29% of Cepheids with an IR excess likely produced by the CSE. Longer period Cepheids do not exhibit greater excess than short periods as previously suspected from observations and theoretical dusty-wind models. Other mechanisms such as free-free emission, among others, may be at the origin of their formation. We also show that not fitting the colour excess leads to a bias on the distance estimates in our Galaxy.

We argue that the magnetically closed corona evolves primarily quasi-statically, punctuated by many localized bursts of activity associated with magnetic reconnection at a myriad of small current sheets. The sheets form by various processes that do not involve a traditional turbulent cascade whereby energy flows losslessly through a continuum of spatial scales starting from the large scale of the photospheric driving. If such an inertial range is a defining characteristic of turbulence, then the magnetically closed corona is not a turbulent system. It nonetheless has a complex structure that bears no direct relationship to the pattern of driving.

We search for ultra-luminous QSOs at high redshift using photometry from the SkyMapper Southern Survey DR3, in combination with 2MASS, VHS DR6, VIKING DR5, AllWISE, and CatWISE2020, as well as parallaxes and proper motions from Gaia DR2 and eDR3. We report 119 newly discovered Southern QSOs, of which 97 with -29 < M_145 < -27 and 4 < z < 5.5 are found an effective search area of 14,120 deg^2. In combination with already known QSOs, we construct a sample that is >90% complete for M_145 < -27.5 at z=4.4 and for M_145 < -28 at z=5.4. This Southern sample has a surface density that is over 3 times higher than in previous searches in the Northern hemisphere, which is only partly due to a more inclusive selection. We derive the bright end of the QSO luminosity function at restframe 145nm and measure its slope as $\beta \approx -3.84$ at z~5. We also present the first z~5 QSO luminosity function at restframe 300nm.

We analyze space-based time series photometry of Sun-like stars, mostly in the Pleiades, but also field stars and the Sun itself. We focus on timescales between roughly 1 hour and 1 day. In the corresponding frequency band these stars display brightness fluctuations with a decreasing power-law continuous spectrum. K2 and Kepler observations show that the RMS flicker due to this Mid-Frequency Continuum (MFC) can reach almost 1%, approaching the modulation amplitude from active regions. The MFC amplitude varies by a factor up to 40 among Pleiades members with similar Teff, depending mainly on the stellar Rossby number Ro. For Ro<0.04, the mean amplitude is roughly constant at about 0.4%; at larger Ro the amplitude decreases rapidly, shrinking by about two orders of magnitude for Ro~1. Among stars, the MFC amplitude correlates poorly with that of modulation from rotating active regions. Among field stars observed for 3 years by Kepler, the quarterly average modulation amplitudes from active regions are much more time-variable than the quarterly MFC amplitudes. We argue that the process causing the MFC is largely magnetic in nature, and that its power-law spectrum comes from magnetic processes distinct from the star's global dynamo, with shorter timescales. By analogy with solar phenomena, we hypothesize that the MFC arises from a (sometimes energetic) variant of the solar magnetic network, perhaps combined with rotation-related changes in the morphology of supergranules.

Previous detections of an X-ray emission line near 3.5 keV in galaxy clusters and other dark matter-dominated objects have been interpreted as observational evidence for the decay of sterile neutrino dark matter. Motivated by this, we report on a search for a 3.5 keV emission line from the Milky Way's galactic dark matter halo with HaloSat. As a single pixel, collimated instrument, HaloSat observations are impervious to potential systematic effects due to grazing incidence reflection and CCD pixelization, and thus may offer a check on possible instrumental systematic errors in previous analyses. We report non-detections of a $\sim$3.5 keV emission line in four HaloSat observations near the Galactic Center. In the context of the sterile neutrino decay interpretation of the putative line feature, we provide 90% confidence level upper limits on the 3.5 keV line flux and 7.1 keV sterile neutrino mixing angle: $F \leq 0.077$ ph cm$^{-2}$ s$^{-1}$ sr$^{-1}$ and $\sin^2(2\theta) \leq 4.25 \times 10^{-11}$. The HaloSat mixing angle upper limit was calculated using a modern parameterization of the Milky Way's dark matter distribution, and in order to compare with previous limits, we also report the limit calculated using a common historical model. The HaloSat mixing angle upper limit places constraints on a number of previous mixing angle estimates derived from observations of the Milky Way's dark matter halo and galaxy clusters, and excludes several previous detections of the line. The upper limits cannot, however, entirely rule out the sterile neutrino decay interpretation of the 3.5 keV line feature.

We conducted a long-slit spectrophotometry analysis to obtain the chemical abundances of seven metal-poor HII regions in three galaxies: UM 160, UM 420, and TOL 0513-393. The data have been taken with the Focal Reducer Low Dispersion Spectrograph 1 (FORS1) at the 8.2-m Very Large Telescope. We derived the physical conditions and the chemical abundances of N, O, Ne, S, Ar, and Cl. We also performed a detailed analysis that involves abundance determinations using the $t^2$ formalism. Based on HeI recombination line intensity ratios, together with the Helio14 code, we derived the abundance of He. In addition, for a value $\Delta Y/\Delta Z_O =3.3\pm 0.7$, we have estimated that the primordial helium abundance by mass is $Y_{\rm P}=0.2448\pm0.0033$. This value agrees with values derived from Standard Big Bang Nucleosynthesis and with other recent determinations of $Y_{\rm P}$.

We explore whether assumptions about dust grain shape affect resulting estimates of the composition and grain size distribution of the AU Microscopii (AU Mic) debris disk from scattered light data collected by Lomax et al. (2018). The near edge-on orientation of the AU Mic debris disk makes it ideal for studying the effect of the scattering phase function (SPF) on the measured flux ratios as a function of wavelength and projected distance. Previous efforts to model the AU Mic debris disk have invoked a variety of dust grain compositions and explored the effect of porosity, but did not undertake a systematic effort to explore a full range of size distributions and compositions to understand possible degeneracies in fitting the data. The degree to which modelling dust grains with more realistic shapes compounds these degeneracies has also not previously been explored. We find differences in the grain properties retrieved depending on the grain shape model used. We also present here our calculations of porous grains of size parameters x = 0.1 to 48 and complex refractive indices (m = n+ik) ranging from n = 1.1 to 2.43 and k = 0 to 1.0, covering multiple compositions at visible and near infrared wavelengths such as ice, silicates, amorphous carbon, and tholins.

Kilonovae are ultraviolet, optical, and infrared transients powered by the radioactive decay of heavy elements following a neutron star merger. Joint observations of kilonovae and gravitational waves can offer key constraints on the source of Galactic $r$-process enrichment, among other astrophysical topics. However, robust constraints on heavy element production requires rapid kilonova detection (within $\sim 1$ day of merger) as well as multi-wavelength observations across multiple epochs. In this study, we quantify the ability of 13 wide field-of-view instruments to detect kilonovae, leveraging a large grid of over 900 radiative transfer simulations with 54 viewing angles per simulation. We consider both current and upcoming instruments, collectively spanning the full kilonova spectrum. The Roman Space Telescope has the highest redshift reach of any instrument in the study, observing kilonovae out to $z \sim 1$ within the first day post-merger. We demonstrate that BlackGEM, DECam, GOTO, the Vera C. Rubin Observatory's LSST, ULTRASAT, and VISTA can observe some kilonovae out to $z \sim 0.1$ ($\sim$475 Mpc), while DDOTI, MeerLICHT, PRIME, $Swift$/UVOT, and ZTF are confined to more nearby observations. Furthermore, we provide a framework to infer kilonova ejecta properties following non-detections and explore variation in detectability with these ejecta parameters.

We extend previous works by considering two additional radio frequencies (K band and X/Ka band) with the aim to study the frequency dependence of the source positions and its potential connection with the physical properties of the underlying AGN. We compared the absolute source positions measured at four different wavelengths, that is, the optical position from the Gaia Early Data Release 3 (EDR3) and the radio positions at the dual S/X, X/Ka combinations and at K band, as available from the third realization of the International Celestial Reference Frame (ICRF3), for 512 common sources. We first aligned the three ICRF3 individual catalogs onto the Gaia EDR3 frame and compare the optical-to-radio offsets before and after the alignment. Then we studied the correlation of optical-to-radio offsets with the observing (radio) frequency, source morphology, magnitude, redshift, and source type. The deviation among optical-to-radio offsets determined in the different radio bands is less than 0.5 mas, but there is statistical evidence that the optical-to-radio offset is smaller at K band compared to S/X band for sources showing extended structures. The optical-to-radio offset was found to statistically correlate with the structure index. Large optical-to-radio offsets appear to favor faint sources but are well explained by positional uncertainty, which is also larger for these sources. We did not detect any statistically significant correlation between the optical-to-radio offset and the redshift. The radio source structure might also be a major cause for the radio-to-optical offset. For the alignment of with the Gaia celestial reference frame, the S/X band frame remains the preferred choice at present.

We present multiwavelength observations of the Type II SN 2020pni. Classified at $\sim 1.3$ days after explosion, the object showed narrow (FWHM $<250\,\textrm{km}\,\textrm{s}^{-1}$) recombination lines of ionized helium, nitrogen, and carbon, as typically seen in flash-spectroscopy events. Using the non-LTE radiative transfer code CMFGEN to model our first high resolution spectrum, we infer a progenitor mass-loss rate of $\dot{M}=(3.5-5.3)\times10^{-3}$ M$_{\odot}$ yr$^{-1}$ (assuming a wind velocity of $v_w=200\,\textrm{km}\,\textrm{s}^{-1}$), estimated at a radius of $R_{\rm in}=2.5\times10^{14}\,\rm{cm}$. In addition, we find that the progenitor of SN 2020pni was enriched in helium and nitrogen (relative abundances in mass fractions of 0.30$-$0.40, and $8.2\times10^{-3}$, respectively). Radio upper limits are also consistent with a dense CSM, and a mass-loss rate of $\dot M>5 \times 10^{-4}\,\rm{M_{\odot}\,yr^{-1}}$. During the first 4 days after first light, we also observe an increase in velocity of the hydrogen lines (from $\sim 250\,\textrm{km}\,\textrm{s}^{-1}$ to $\sim 1000\,\textrm{km}\,\textrm{s}^{-1}$), which suggests a complex CSM. The presence of dense and confined CSM, as well as its inhomogeneous structure, suggest a phase of enhanced mass loss of the progenitor of SN 2020pni during the last year before explosion. Finally, we compare SN 2020pni to a sample of other shock-photoionization events. We find no evidence of correlations among the physical parameters of the explosions and the characteristics of the CSM surrounding the progenitors of these events. This favors the idea that the mass-loss experienced by massive stars during their final years could be governed by stochastic phenomena, and that, at the same time, the physical mechanisms responsible for this mass-loss must be common to a variety of different progenitors.

We investigate the $H_0$ tension in a range of extended model frameworks beyond the standard $\Lambda$CDM without the data from cosmic microwave background (CMB). Specifically, we adopt the data from baryon acoustic oscillation, big bang nucleosynthesis and type Ia supernovae as indirect measurements of $H_0$ to study the tension. We show that the estimated value of $H_0$ from indirect measurements is overall lower than that from direct local ones regardless of the data sets and a range of extended models to be analyzed, which indicates that, although the significance of the tension varies depending on models, the $H_0$ tension persists in a broad framework beyond the standard $\Lambda$CDM model even without CMB data.

We present a systematic spectral-timing analysis of a fast appearance/disappearance of a type-B quasi-periodic oscillation (QPO), observed in four NICER observations of MAXI J1348-630. By comparing the spectra of the period with and without the type-B QPO, we found that the main difference appears at energy bands above ~2 keV, suggesting that the QPO emission is dominated by the hard Comptonised component. During the transition, a change in the relative contribution of the disk and Comptonised emission was observed. The disk flux decreased while the Comptonised flux increased from non-QPO to type-B QPO. However, the total flux did not change too much in the NICER band. Our results reveal that the type-B QPO is associated with a redistribution of accretion power between the disk and Comptonised emission. When the type-B QPO appears, more accretion power is dissipated into the Comptonised region than in the disk. Our spectral fits give a hint that the increased Comptonised emission may come from an additional component that is related to the base of the jet.

Using the short cadence data (1-minute interval) of the Kepler space telescope, we conducted a statistical analysis for the characteristic time of stellar flares on Sun-like stars (SLS). Akin to solar flares, stellar flares show rise and decay light curve profile, which reflects two distinct phases (rise phase and decay phase) of flare process. We derived the characteristic times of the two phases for the stellar flares of SLS, resulting the median rise time of about 5.9 minutes and the median decay time of 22.6 minutes. It is found that both the rise time and the decay time of the stellar flares follow the log-normal distribution. The peak positions of the log-normal distributions for flare rise time and decay time are 3.5 minutes and 14.8 minutes, respectively. These time values of stellar flares are similar to the timescale of solar flares, which supports that stellar flares and solar flares have the same physical mechanism. The statistical results obtained in this work for SLS can be a benchmark of flare characteristic times when comparing with other types of stars.

Context. Luminous Blue Variables (LBVs) are thought to be in a transitory phase between O stars on the main-sequence and the Wolf-Rayet stage. Recent studies suggest that they might be formed through binary interaction. Only a few are known in binary systems but their multiplicity fraction is uncertain. Aims. This study aims at deriving the binary fraction among the Galactic (confirmed and candidate) LBV population. We combine multi-epoch spectroscopy and long-baseline interferometry. Methods. We use cross-correlation to measure their radial velocities. We identify spectroscopic binaries through significant RV variability (larger than 35 km/s). We investigate the observational biases to establish the intrinsic binary fraction. We use CANDID to detect interferometric companions, derive their parameters and positions. Results. We derive an observed spectroscopic binary fraction of 26 %. Considering period and mass ratio ranges from Porb=1 to 1000 days, and q = 0.1-1.0, and a representative set of orbital parameter distributions, we find a bias-corrected binary fraction of 62%. From interferometry, we detect 14 companions out of 18 objects, providing a binary fraction of 78% at projected separations between 1 and 120 mas. From the derived primary diameters, and the distances of these objects, we measure for the first time the exact radii of Galactic LBVs to be between 100 and 650 Rsun, making unlikely to have short-period systems. Conclusions. This analysis shows that the binary fraction among the Galactic LBV population is large. If they form through single-star evolution, their orbit must be initially large. If they form through binary channel that implies that either massive stars in short binary systems must undergo a phase of fully non-conservative mass transfer to be able to sufficiently widen the orbit or that LBVs form through merging in initially binary or triple systems.

The guest star of AD 1181 is the only historical supernova of the last millennium that is without a definite counterpart. The previously proposed association with supernova remnant 3C58 is in strong doubt because of the inferred age of this remnant. Here we report a new identification of SN 1181 with our codiscovery of the hottest known Wolf Rayet star of the Oxygen sequence (dubbed Parkers star) and its surrounding nebula Pa 30. Our spectroscopy of the nebula shows a fast shock with extreme velocities of approx. 1,100kms. The derived expansion age of the nebula implies an explosive event approx 1,000 years ago which agrees with the 1181 AD event. The on sky location also fits the historical Chinese and Japanese reports of SN 1181 to 3.5degrees. Pa 30 and Parkers star have previously been proposed to be the result of a double-degenerate merger, leading to a rare Type Iax supernova. The likely historical magnitude and the distance suggest the event was subluminous for normal supernova. This agrees with the proposed Type Iax association which would also be the first of its kind in the Galaxy. Taken together, the age, location, event magnitude and duration elevate Pa 30 to prime position as the counterpart of SN 1181. This source is the only Type Iax supernova where detailed studies of the remnant star and nebula are possible. It provides strong observational support for the double degenerate merger scenario for Type Iax supernovae.

The Coronal Multi-channel Polarimeter (CoMP) has opened up exciting opportunities to probe transverse MHD waves in the Sun's corona. The archive of CoMP data is utilised to generate a catalogue of quiescent coronal loops that can be used for studying propagating kink waves. The catalogue contains 120 loops observed between 2012-2014. This catalogue is further used to undertake a statistical study of propagating kink waves in the quiet regions of the solar corona, investigating phase speeds, loop lengths, footpoint power ratio and equilibrium parameter values. The statistical study enables us to establish the presence of a relationship between the rate of damping and the length of the coronal loop, with longer coronal loops displaying weaker wave damping. We suggest the reason for this behaviour is related to a decreasing average density contrast between the loop and ambient plasma as loop length increases. The catalogue presented here will provide the community with the foundation for the further study of propagating kink waves in the quiet solar corona.

We review the main physical processes that lead to the formation of stellar binary black holes (BBHs) and to their merger. BBHs can form from the isolated evolution of massive binary stars. The physics of core-collapse supernovae and the process of common envelope are two of the main sources of uncertainty about this formation channel. Alternatively, two black holes can form a binary by dynamical encounters in a dense star cluster. The dynamical formation channel leaves several imprints on the mass, spin and orbital properties of BBHs.

Understanding oscillations and waves in astrophysical fluid bodies helps to elucidate their observed variability and the underlying physical mechanisms. Indeed, global oscillations and bending modes of accretion discs or tori may be relevant to quasi-periodicity and warped structures around compact objects. While most studies rely on linear theory, observationally significant, nonlinear dynamics is still poorly understood, especially in Keplerian discs for which resonances typically demand a separate treatment. In this work we introduce a novel analytical model which exactly solves the ideal, compressible fluid equations for a non-self-gravitating elliptical cylinder within a local shearing sheet. The aspect ratio of the ring is an adjustable parameter, allowing a continuum of models ranging from a torus of circular cross-section to a thin ring. We restrict attention to flow fields which are a linear function of the coordinates, capturing the lowest order global motions and reducing the dynamics to a set of coupled ordinary differential equations (ODEs). This system acts as a framework for exploring a rich range of hydrodynamic phenomena in both the large amplitude and Keplerian regimes. We demonstrate the connection between tilting tori and warped discs within this model, showing that the linear modes of the ring correspond to oppositely precessing global bending modes. These are further confirmed within a numerical grid based simulation. Crucially, the ODE system developed here allows for a more tractable investigation of nonlinear dynamics. This will be demonstrated in a subsequent paper which evidences mode coupling between warping and vertical motions in thin tilted rings.

Various attempts have been made in the literature at describing the origin and the physical mechanisms behind flaring events in blazars with radiative emission models, but detailed properties of multi-wavelength (MWL) light curves still remain difficult to reproduce. We have developed a versatile radiative code, based on a time-dependent treatment of particle acceleration, escape and radiative cooling, allowing us to test different scenarios to connect the continuous low-state emission self-consistently with that during flaring states. We consider flares as weak perturbations of the quiescent state and apply this description to the February 2010 MWL flare of Mrk 421, the brightest Very High Energy (VHE) flare ever detected from this archetypal blazar, focusing on interpretations with a minimum number of free parameters. A general criterion is obtained, which disfavours a one-zone model connecting low and high state under our assumptions. A two-zone model combining physically connected acceleration and emission regions yields a satisfactory interpretation of the available time-dependent MWL light curves and spectra of Mrk 421, although certain details remain difficult to reproduce. The two-zone scenario finally proposed for the complex quiescent and flaring VHE emitting region involves both Fermi-I and Fermi-II acceleration mechanisms, respectively at the origin of the quiescent and flaring emission.

The volume available on small satellites restricts the size of optical apertures to a few centimetres, limiting the Ground-Sampling Distance (GSD) in the visible to typically 3 m at 500 km. We present in this paper the latest development of a laboratory demonstrator of a segmented deployable telescope that will triple the achievable ground resolution and improve photometric capability of CubeSat imagers. Each mirror segment is folded for launch and unfolds in space. We demonstrate through laboratory validation very high deployment repeatability of the mirrors <{\pm}5 {\mu}m. To enable diffraction-limited imaging, segments are controlled in piston, tip, and tilt. This is achieved by an initial coarse alignment of the mirrors followed by a fine phasing step. Finally, we investigate the impact of the thermal environment on high-order wavefront error and the conceptual design of a deployable secondary fitting inside 1U.

We study the molecular gas content of 24 star-forming galaxies at $z=3-4$, with a median stellar mass of $10^{9.1}$ M$_{\odot}$, from the MUSE Hubble Ultra Deep Field (HUDF) Survey. Selected by their Lyman-alpha-emission and H-band magnitude, the galaxies show an average EW $\approx 20$ angstrom, below the typical selection threshold for Lyman Alpha Emitters (EW $> 25$ angstrom), and a rest-frame UV spectrum similar to Lyman Break Galaxies. We use rest-frame optical spectroscopy from KMOS and MOSFIRE, and the UV features observed with MUSE, to determine the systemic redshifts, which are offset from Lyman alpha by 346 km s$^{-1}$, with a 100 to 600 km s$^{-1}$ range. Stacking CO(4-3) and [CI](1-0) (and higher-$J$ CO lines) from the ALMA Spectroscopic Survey of the HUDF (ASPECS), we determine $3\sigma$ upper limits on the line luminosities of $4.0\times10^{8}$ K km s$^{-1}$pc$^{2}$ and $5.6\times10^{8}$ K km s$^{-1}$pc$^{2}$, respectively (for a 300 km s$^{-1}$ linewidth). Stacking the 1.2 mm and 3 mm dust continuum flux densities, we find a $3\sigma$ upper limits of 9 $\mu$Jy and $1.2$ $\mu$Jy, respectively. The inferred gas fractions, under the assumption of a 'Galactic' CO-to-H$_{2}$ conversion factor and gas-to-dust ratio, are in tension with previously determined scaling relations. This implies a substantially higher $\alpha_{\rm CO} \ge 10$ and $\delta_{\rm GDR} \ge 1200$, consistent with the sub-solar metallicity estimated for these galaxies ($12 + \log(O/H) \approx 7.8 \pm 0.2$). The low metallicity of $z \ge 3$ star-forming galaxies may thus make it very challenging to unveil their cold gas through CO or dust emission, warranting further exploration of alternative tracers, such as [CII].

We briefly present the history of technical solutions aimed at improving the efficiency of spectroscopy on small- and moderate-diameter telescopes. We assess the current state of spectroscopy techniques and some of the perspectives.

The gamma-ray emission of RX J1713.7-3946, despite being extensively studied in the GeV and TeV domain, remains poorly understood. This is mostly because in this range, two competing mechanisms can efficiently produce gamma rays: the inverse Compton scattering of accelerated electrons, and interactions of accelerated protons with nuclei of the ISM. In addition to the acceleration of particles from the thermal pool, the reacceleration of pre-existing CRs is often overlooked, and shall in fact also been taken into account. Especially, because of the distance to the SNR ($\sim 1$ kpc), and the low density in which the shock is currently expanding ($\sim 10^{-2}$ cm$^{-3}$), the re-acceleration of CR electrons pre-existing in the ISM, can account for a significant fraction of the observed gamma-ray emission, and contribute to the shaping of the spectrum in the GeV-TeV range. Remarkably, this emission of leptonic origin is found to be close to the level of the gamma-ray signal in the TeV range, provided that the spectrum of pre-exisiting cosmic ray electrons is similar to the one observed in the local interstellar medium. The overall gamma-ray spectrum of RX J1713.7-3946 is naturally produced as the sum of a leptonic emission from reaccelerated CR electrons, and a subdominant hadronic emission from accelerated protons. We also argue that neutrino observations with next-generation detectors might lead to a detection even in the case of a lepto-hadronic origin of the gamma-ray emission.

Coronal loops form the basic building blocks of the magnetically closed solar corona yet much is still to be determined concerning their possible fine-scale structuring and the rate of heat deposition within them. Using an improved multi-stranded loop model to better approximate the numerically challenging transition region, this paper examines synthetic NASA Solar Dynamics Observatory's (SDO) Atmospheric Imaging Assembly (AIA) emission simulated in response to a series of prescribed spatially and temporally random, impulsive and localised heating events across numerous sub-loop elements with a strong weighting towards the base of the structure; the nanoflare heating scenario. The total number of strands and nanoflare repetition times are varied systematically in such a way that the total energy content remains approximately constant across all the cases analysed. Repeated time lag detection during an emission time series provides a good approximation for the nanoflare repetition time for low-frequency heating. Furthermore, using a combination of AIA 171/193 and 193/211 channel ratios in combination with spectroscopic determination of the standard deviation of the loop apex temperature over several hours alongside simulations from the outlined multi-stranded loop model, it is demonstrated that both the imposed heating rate and number of strands can be realised.

Some generalizations of the relation between high-energy astrophysical neutrino and cosmic ray fluxes are obtained, taking into account present results on the cosmic ray spectrum and composition as well as a more realistic modeling of the Galactic and extragalactic cosmic ray components down to PeV energies. It is found that the level of neutrino fluxes measured by IceCube can be consistent with sources that are thin to escaping protons. This could also make it easier for heavier nuclei to be emitted from the sources without suffering excessive disintegration processes.

Extensions have been made recently to the coronal approximation for the purpose of modelling line emission from carbon and oxygen in the lower solar atmosphere. The same modelling is used here for other elements routinely observed in the solar transition region: N, Ne, Mg, Si and S. The modelling includes the effects of higher densities suppressing dielectronic recombination and populating long-lived, metastable levels; the presence of metastable levels typically causes effective ionisation rates to increase and recombination rates to decrease. Processes induced by the radiation field, namely photo-ionisation and photo-excitation, have been included, along with charge transfer, which occurs when electrons are exchanged during atom-ion and ion-ion collisions. The resulting ion balances are shown, and indicate significant changes compared to the frequently-employed coronal approximation. The effect on level populations within ions caused by photo-excitation is also assessed. To give an illustration of how line emission could be altered by these processes, selected line contribution functions are presented at the end.

The NASA Kepler and follow-on K2 missions (2009-2018) left a legacy of data and discoveries, finding thousands of exoplanets, and also obtaining high-precision long time-series data for hundreds of thousands of stars, including many types of pulsating variables. Here we highlight a few of the ongoing discoveries from Kepler data on $\delta$ Scuti pulsating variables, which are core hydrogen-burning stars of about twice the mass of the Sun. We discuss many unsolved problems surrounding the properties of the variability in these stars, and the progress enabled by Kepler data in using pulsations to infer their interior structure, a field of research known as asteroseismology.

We propose a novel mechanism for enhancing the primordial gravitational waves without significantly affecting the curvature perturbations produced during inflation. This is achieved due to non-linear sourcing of resonantly amplified scalar field fluctuations. Our result is an explicit scale-dependent counter-example of the famous Lyth bound, which opens up a promising perspective of producing detectable inflationary tensor modes with low-scale inflation and a sub-Planckian field excursion. We explicitly demonstrate the testability of our mechanism with upcoming Cosmic Microwave Background B-mode observations.

We re-examine the XUV luminosity evolution of TRAPPIST-1 utilizing new observational constraints (XUV and bolometric luminosity) from multi-epoch X-ray/UV photometry. Following the formalism presented in Fleming et al. (2020), we infer that TRAPPIST-1 maintained a saturated XUV luminosity, relative to the bolometric luminosity, of $\log_{10}$(L$_{\rm XUV}$/L$_{\rm bol}$) $= -3.03_{-0.23}^{+0.25}$ at early times for a period of $t_{sat} = 3.14_{-1.46}^{+2.22}$ Gyr. After the saturation phase, we find L$_{\rm XUV}$ decayed over time by an exponential rate of $\beta_{\rm XUV} = -1.17_{-0.28}^{+0.27}$. Compared to our inferred age of the system, ${\rm age} = 7.96_{-1.87}^{+1.78}$ Gyr, our result for $t_{sat}$ suggests that there is only a $\sim4\%$~chance that TRAPPIST-1 still remains in the saturated phase today, which is significantly lower than the previous estimate of 40\%. Despite this reduction in $t_{sat}$, our results remain consistent in the conclusion that the TRAPPIST-1 planets likely received an extreme amount XUV energy -- an estimated integrated XUV energy of $\sim10^{30}-10^{32}$ erg over the star's lifetime.

Modifications were made to the Sloan Digital Sky Survey's Baryonic Oscillations Spectroscopic Survey (SDSS/BOSS) optical fibres assigned to quasar targets in order to improve the signal-to-noise ratio in the Lyman-alpha forest. However, the penalty for these modifications is that quasars observed in this way require additional flux correction procedures in order to recover the correct spectral shapes. In this paper we describe such a procedure, based on the geometry of the problem and other observational parameters. Applying several correction methods to four SDSS quasars with multiple observations permits a detailed check on the relative performances of the different flux correction procedures. We contrast our method (which takes into account a wavelength dependent seeing profile) with the BOSS pipeline approach (which does not). Our results provide independent confirmation that the geometric approach employed in the SDSS pipeline works well, although with room for improvement. By separating the contributions from four effects we are able to quantify their relative importance. Most importantly, we demonstrate that wavelength dependence has a significant impact on the derived spectral shapes and thus should not be ignored.

We have performed Monte Carlo simulations of the trajectories of several runaway stars using their parallaxes and proper motions from the Gaia EDR3 catalogue. We have confirmed the hypothesis that the stars AE Aur and $\mu$Col are a product of the multiple system breakup $\sim$2.5 Myr ago and the Orion Trapezium may be the parent cluster for this pair of stars. We show that the data from the Gaia EDR3 catalogue for the star $\iota$Ori, mainly the parallax, do not allow us to talk about the breakup of the multiple system of AE Aur, $\mu$Col, and $\iota$Ori. The existence of close pair encounters between the stars HD 30112 and HD 43112 $\sim$1 Myr ago has been confirmed. Close triple encounters confirm the hypothesis that the stars HD 30112 and HD 43112 escaped from the parent cluster Col 69. We show that the stars HIP 28133 and TYC 5368-1541-1 have a nonzero probability of escape from the region within 10 pc of the center of the Orion Trapezium cluster and a fairly high probability (about 8\%) that they were both at distances less than 20 pc from the center of the Orion Trapezium $\sim$2.5 Myr ago. It has been established for the first time that the stars Gaia EDR3 3021115184676332288 and Gaia EDR3 2983790269606043648 have a probability of about 0.5\% that they broke up as a binary system $\sim$1.1 Myr ago. The star Gaia EDR3 3021115184676332288 has a probability of about 16\% that it escaped from the region within 10 pc of the center of the Orion Trapezium cluster $\sim$1 Myr ago.

The cosmological model-independent method Gaussian process (GP) has been widely used in the reconstruction of Hubble constant $H_0$, and the hyperparameters inside GP influence the reconstructed result derived from GP. Different hyperparameters inside GP are used in the constraint of $H_0$ derived from GP with observational Hubble parameter $H(z)$ data (OHD), and the influence of the hyperparameters inside GP on the reconstruction of $H_0$ with GP is discussed. The discussion about the hyperparameters inside GP and the forecasts for future data show that the consideration of the lower and upper bounds on the GP's hyperparameters are necessary in order to get an extrapolated result of $H_0$ from GP reliably and robustly.

High-redshift luminous quasars powered by accreting supermassive black holes (SMBHs) with mass $\gtrsim 10^9 M_\odot$ constrain their formation pathways. We investigate the formation of heavy seeds of SMBHs through gas collapse in the quasar host progenitors, using merger trees to trace the halo growth in highly-biased, overdense regions of the universe. The progenitor halos are likely irradiated by intense H$_2$-photodissociating radiation from nearby star-forming galaxies and heat the interior gas by successive mergers. The kinetic energy of the gas originating from mergers as well as baryonic streaming motion prevents gas collapse and delays prior star formation. With a streaming velocity higher than the root-mean-square value, gas clouds in nearly all $10^4$ realizations of merger trees enter the atomic-cooling stage and begin to collapse isothermally with $T \simeq 8000 K$ via Ly$\alpha$ cooling. The fraction of trees which host isothermal gas collapse is $14\%$ and increases with streaming velocity, while the rest form H$_2$-cooled cores after short isothermal phases. If the collapsing gas is enriched to $Z_{crit}\sim 2\times 10^{-3} Z_\odot$, requiring efficient metal mixing, this fraction could be reduced by additional cooling via metal fine-structure lines. In the massive collapsing gas, the accretion rate onto a newly-born protostar ranges between $3 \times 10^{-3}-5 M_\odot yr^{-1}$, among which a large fraction exceeds the critical rate suppressing stellar radiative feedback. As a result, we expect a distribution of stellar mass (presumably BH mass) ranging from several hundred to above $10^5 M_\odot$, potentially forming massive BH binary mergers and yielding gravitational wave events.

We apply for the first time Gaussian Process Regression (GPR) as a foreground removal technique in the context of single-dish, low redshift HI intensity mapping, and present an open-source python toolkit for doing so. We use MeerKAT and SKA1-MID-like simulations of 21cm foregrounds (including polarisation leakage), HI cosmological signal and instrumental noise. We find that it is possible to use GPR as a foreground removal technique in this context, and that it is better suited in some cases to recover the HI power spectrum than Principal Component Analysis (PCA), especially on small scales. GPR is especially good at recovering the radial power spectrum, outperforming PCA when considering the full bandwidth of our data. Both methods are worse at recovering the transverse power spectrum, since they rely on frequency-only covariance information. When halving our data along frequency, we find that GPR performs better in the low frequency range, where foregrounds are brighter. It performs worse than PCA when frequency channels are missing, to emulate RFI flagging. We conclude that GPR is an excellent foreground removal option for the case of single-dish, low redshift HI intensity mapping. Our python toolkit gpr4im and the data used in this analysis are publicly available on GitHub. The GitHub symbol in the caption of each figure links to a jupyter notebook showing how the figure was produced.

We study eight different gamma-ray burst (GRB) data sets to examine whether current GRB measurements -- that probe a largely unexplored part of cosmological redshift ($z$) space -- can be used to reliably constrain cosmological model parameters. We use three Amati-correlation samples and five Combo-correlation samples to simultaneously derive correlation and cosmological model parameter constraints. The intrinsic dispersion of each GRB data set is taken as a goodness measurement. We examine the consistency between the cosmological bounds from GRBs with those determined from better-established cosmological probes, such as baryonic acoustic oscillation (BAO) and Hubble parameter $H(z)$ measurements. We use the Markov chain Monte Carlo method implemented in \textsc{MontePython} to find best-fit correlation and cosmological parameters, in six different cosmological models, for the eight GRB samples, alone or in conjunction with BAO and $H(z)$ data. For the Amati correlation case, we compile a data set of 118 bursts, the A118 sample, which is the largest -- about half of the total Amati-correlation GRBs -- current collection of GRBs suitable for constraining cosmological parameters. This updated GRB compilation has the smallest intrinsic dispersion of the three Amati-correlation GRB data sets we examined. We are unable to define a collection of reliable bursts for current Combo-correlation GRB data. Cosmological constraints determined from the A118 sample are consistent with -- but significantly weaker than -- those from BAO and $H(z)$ data. They also are consistent with the spatially-flat $\Lambda$CDM model as well as with dynamical dark energy models and non-spatially-flat models. Since GRBs probe a largely unexplored region of $z$, it is well worth acquiring more and better-quality burst data which will give a more definitive answer to the question of the title.

(16) Psyche is the largest M-type asteroid in the main belt and the target of the NASA Discovery-class Psyche mission. Despite gaining considerable interest in the scientific community, Psyche's composition and formation remain unconstrained. Originally, Psyche was considered to be almost entirely composed of metal due to its high radar albedo and spectral similarities to iron meteorites. More recent telescopic observations suggest the additional presence of low-Fe pyroxene and exogenic carbonaceous chondrites on the asteroid's surface. To better understand the abundances of these additional materials, we investigated visible near-infrared (0.35 - 2.5 micron) spectral properties of three-component laboratory mixtures of metal, low-Fe pyroxene, and carbonaceous chondrite. We compared the band depths and spectral slopes of these mixtures to the telescopic spectrum of (16) Psyche to constrain material abundances. We find that the best matching mixture to Psyche consists of 82.5% metal, 7% low-Fe pyroxene, and 10.5% carbonaceous by weight, suggesting that the asteroid is less metallic than originally estimated (~94%). The relatively high abundance of carbonaceous chondrite material estimated from our laboratory experiments implies the delivery of this exogenic material through low velocity collisions to Psyche's surface. Assuming that Psyche's surface is representative of its bulk material content, our results suggest a porosity of 35% to match recent density estimates.

To evaluate the role of Total Solar Irradiance (TSI) on Northern Hemisphere (NH) surface air temperature trends it is important to have reliable estimates of both quantities. 16 different TSI estimates were compiled from the literature. 1/2 of these estimates are low variability and 1/2 are high variability. 5 largely-independent methods for estimating NH temperature trends were evaluated using: 1) only rural weather stations; 2) all available stations whether urban or rural (the standard approach); 3) only sea surface temperatures; 4) tree-ring temperature proxies; 5) glacier length temperature proxies. The standard estimates using urban as well as rural stations were anomalous as they implied a much greater warming in recent decades than the other estimates. This suggests urbanization bias might still be a problem in current global temperature datasets despite the conclusions of some earlier studies. Still, all 5 estimates confirm it is currently warmer than the late 19th century, i.e., there has been some global warming since 1850. For the 5 estimates of NH temperatures, the contribution from direct solar forcing for all 16 estimates of TSI was evaluated using simple linear least-squares fitting. The role of human activity in recent warming was then calculated by fitting the residuals to the UN IPCC's recommended anthropogenic forcings time series. For all 5 NH temperature series, different TSI estimates implied everything from recent global warming being mostly human-caused to it being mostly natural. It seems previous studies (including the most recent IPCC reports) that had prematurely concluded the former failed to adequately consider all the relevant estimates of TSI and/or to satisfactorily address the uncertainties still associated with NH temperature trend estimates. Several recommendations are provided on how future research could more satisfactorily resolve these issues.

In order to solve the hierarchy problem, the relaxion must remain trapped in the correct minimum, even if the electroweak symmetry is restored after reheating. In this scenario, the relaxion starts rolling again until the back-reaction potential, with its set of local minima, reappears. Depending on the time of barrier-reappearance, Hubble friction alone may be insufficient to re-trap the relaxion in a large portion of the parameter space. Thus, an additional source of friction is required, which might be provided by coupling to a dark photon. The dark photon experiences a tachyonic instability as the relaxion rolls, which slows down the relaxion by backreacting to its motion, and efficiently creates anisotropies in the dark photon energy-momentum tensor, sourcing gravitational waves. We calculate the spectrum of the resulting gravitational wave background from this new mechanism, and evaluate its observability by current and future experiments. We further investigate the possibility that the coherently oscillating relaxion constitutes dark matter and present the corresponding constraints from gravitational waves.

Solutions to the two-body problem in general relativity allow us to predict the mass, spin and recoil velocity of a black-hole merger remnant given the masses and spins of its binary progenitors. In this paper we address the inverse problem: given a binary black-hole merger, can we use the parameters measured by gravitational-wave interferometers to tell if the binary components are of hierarchical origin, i.e. if they are themselves remnants of previous mergers? If so, can we determine at least some of the properties of their parents? This inverse problem is in general overdetermined. We show that hierarchical mergers occupy a characteristic region in the plane composed of the effective spin parameters $\chi_{\rm eff}$ and $\chi_{\rm p}$, and therefore a measurement of these parameters can add weight to the hierarchical-merger interpretation of some gravitational-wave events, including GW190521. If one of the binary components has hierarchical origin and its spin magnitude is well measured, we derive exclusion regions on the properties of its parents: for example we infer that the parents of GW190412 (if hierarchical) must have had unequal masses and low spins. Our formalism is quite general, and it can be used to infer constraints on the astrophysical environment producing hierarchical mergers.

Atom-in-jellium calculations of the Einstein frequency in condensed matter and of the equation of state were used to predict the variation of shear modulus from zero pressure to $\sim 10^7$ g/cm$^3$, for several elements relevant to white dwarf (WD) stars and other self-gravitating systems. This is by far the widest range reported electronic structure calculation of shear modulus, spanning from ambient through the one-component plasma to extreme relativistic conditions. The predictions were based on a relationship between Debye temperature and shear modulus which we assess to be accurate at the $o(10\%)$ level, and is the first known use of atom-in-jellium theory to calculate a shear modulus. We assessed the overall accuracy of the method by comparing with experimental measurements and more detailed electronic structure calculations at lower pressures.

The recent discovery of binary neutron star mergers has opened a new and exciting venue of research into hot and dense strongly interacting matter. For the first time this elusive state of matter, described by the theory of quantum chromo dynamics, can be studied in two very different environments. On the macroscopic scale in the collisions of neutron stars and on the microscopic scale in collisions of heavy ions at particle collider facilities. We will discuss the conditions that are created in these mergers and the corresponding high energy nuclear collisions. This includes the properties of QCD matter, i.e. the expected equation of state as well as expected chemical and thermodynamic properties of this exotic matter. To explore this matter in the laboratory - a new research prospect is available at the Facility for Antiproton and Ion Research, FAIR. The new facility is being constructed adjacent to the existing accelerator complex of the GSI Helmholtz Centre for Heavy Ion Research at Darmstadt/Germany, expanding the research goals and technical possibilities substantially. The worldwide unique accelerator and experimental facilities of FAIR will open the way for a broad spectrum of unprecedented research supplying a variety of experiments in hadron, nuclear, atomic and plasma physics as well as biomedical and material science which will be briefly described.

Mike Lockwood and Mathew Owens discuss how eclipse observations are aiding the development of a climatology of near-Earth space

We explore the multi-faceted important features of turbulence (e.g., anisotropy, dispersion, diffusion) in the three-dimensional (3D) wavenumber domain ($k_\parallel$, $k_{\perp,1}$, $k_{\perp,2}$), by employing the k-filtering technique to the high-quality measurements of fields and particles from the MMS multi-spacecraft constellation. We compute the 3D power spectral densities (PSDs) of magnetic and electric fluctuations (marked as $\rm{PSD}(\delta \mathbf{B}(\mathbf{k}))$ and $\rm{PSD}(\delta \mathbf{E}'_{\langle\mathbf{v}_\mathrm{i}\rangle}(\mathbf{k}))$), both of which show a prominent spectral anisotropy in the sub-ion range. We give the first 3D image of the bifurcation between power spectra of the electric and magnetic fluctuations, by calculating the ratio between $\rm{PSD}(\delta \mathbf{E}'_{ \langle\mathbf{v}_\mathrm{i}\rangle}(\mathbf{k}))$ and $\rm{PSD}(\delta \mathbf{B}(\mathbf{k}))$, the distribution of which is related to the non-linear dispersion relation. We also compute the ratio between electric spectra in different reference frames defined by the ion bulk velocity, that is $\mathrm{PSD}(\delta{\mathbf{E}'_{\mathrm{local}\ \mathbf{v}_\mathrm{i}}})/\mathrm{PSD}(\delta{\mathbf{E}'_{ \langle\mathbf{v}_\mathrm{i}\rangle}})$, to visualize the turbulence ion diffusion region (T-IDR) in wavenumber space. The T-IDR has an anisotropy and a preferential direction of wavevectors, which is generally consistent with the plasma wave theory prediction based on the dominance of kinetic Alfv\'en waves (KAW). This work manifests the worth of the k-filtering technique in diagnosing turbulence comprehensively, especially when the electric field is involved.

Quantum coherence is one of the most striking features of quantum mechanics rooted in the superposition principle. Recently it has been demonstrated that it is possible to harvest the quantum coherence from a coherent scalar field. In order to explore a new method of detecting axion dark matter, we quantify a coherent measure of a detector and show that the detector can harvest the quantum coherence from the axion dark matter.