Papers for Monday, Jul 25 2022

We detail the follow up and characterization of a transiting exo-Venus identified by TESS, GJ 3929b, (TOI-2013b) and its non-transiting companion planet, GJ 3929c (TOI-2013c). GJ 3929b is an Earth-sized exoplanet in its star's Venus-zone (P$_{b}$ = 2.616272 $\pm$ 0.000005 days; S$_{b}$ = 17.3$^{+0.8}_{-0.7}$ S$_{\oplus}$) orbiting a nearby M dwarf. GJ 3929c is most likely a non-transiting sub-Neptune. Using the new, ultra-precise NEID spectrometer on the WIYN 3.5 m Telescope at Kitt Peak National Observatory, we are able to modify the mass constraints of planet b reported in previous works and consequently improve the significance of the mass measurement to almost 4$\sigma$ confidence (M$_{b}$ = 1.75 $\pm$ 0.45 M$_{\oplus}$). We further adjust the orbital period of planet c from its alias at 14.30 $\pm$ 0.03 days to the likely true period of 15.04 $\pm$ 0.03 days, and adjust its minimum mass to m$\sin i$ = 5.71 $\pm$ 0.92 M$_{\oplus}$. Using the diffuser-assisted ARCTIC imager on the ARC 3.5 m telescope at Apache Point Observatory, in addition to publicly available TESS and LCOGT photometry, we are able to constrain the radius of planet b to R$_{p}$ = 1.09 $\pm$ 0.04 R$_{\oplus}$. GJ 3929b is a top candidate for transmission spectroscopy in its size regime (TSM = 14 $\pm$ 4), and future atmospheric studies of GJ 3929b stand to shed light on the nature of small planets orbiting M dwarfs.

In the next decade, the Laser Interferometer Space Antenna (LISA) will detect the coalescence of massive black hole binaries (MBHBs) in the range $[10^4, 10^8] \, \rm M_{\odot}$, up to $z\sim10$. Their gravitational wave (GW) signal is expected to be accompanied by an electromagnetic counterpart (EMcp), generated by the gas accreting on the binary or on the remnant BH. In this work, we present the number and characteristics (such as redshift and mass distribution, apparent magnitudes or fluxes) of EMcps detectable jointly by LISA and some representative EM telescopes. We combine state-of-the-art astrophysical models for the galaxies formation and evolution to build the MBHBs catalogues, with Bayesian tools to estimate the binary sky position uncertainty from the GW signal. Exploiting additional information from the astrophysical models, such as the amount of accreted gas and the BH spins, we evaluate the expected EM emission in the soft X-ray, optical and radio bands. Overall, we predict between 7 and 21 EMcps in 4 yrs of joint observations by LISA and the considered EM facilities, depending on the astrophysical model. We also explore the impact of the hydrogen and dust obscuration of the optical and X-ray emissions, as well as of the collimation of the radio emission: these effects reduce the number to EMcps to 2 or 3, depending on the astrophysical model, again in 4 yrs of observations. Most of the EMcps are characterised by faint EM emission, challenging the observational capabilities of future telescopes. Finally, we also find that systems with multi-modal sky position posterior distributions represent only a minority of cases and do not affect significantly the number of EMcps.

We examine the effect of neutrino flavor transformation by the fast flavor instability (FFI) on long-term mass ejection from accretion disks formed after neutron star mergers. Neutrino emission and absorption in the disk set the composition of the disk ejecta, which subsequently undergoes $r$-process nucleosynthesis upon expansion and cooling. Here we perform 28 time-dependent, axisymmetric, viscous-hydrodynamic simulations of accretion disks around hypermassive neutron stars (HMNSs) of variable lifetime, using a 3-species neutrino leakage scheme for emission and an annular-lightbulb scheme for absorption. We include neutrino flavor transformation due the FFI in a parametric way, by modifying the absorbed neutrino fluxes and temperatures, allowing for flavor mixing at various levels of flavor equilibration, and also in a way that aims to respect the lepton-number preserving symmetry of the neutrino self-interaction Hamiltonian. We find that for a promptly-formed black hole (BH), the FFI lowers the average electron fraction of the disk outflow due to a decrease in neutrino absorption, driven primarily by a drop in electron neutrino/antineutrino flux upon flavor mixing. For a long-lived HMNS, the disk emits more heavy lepton neutrinos and reabsorbs more electron neutrinos than for a BH, with a smaller drop in flux compensated by a higher neutrino temperature upon flavor mixing. The resulting outflow has a broader electron fraction distribution, a more proton-rich peak, and undergoes stronger radiative driving. Disks with intermediate HMNS lifetimes show results that fall in between these two limits. In most cases, the impact of the FFI on the outflow is moderate, with changes in mass ejection, average velocity, and average electron fraction of order $\sim 10\%$, and changes in the lanthanide/actinide mass fraction of up to a factor $\sim 2$.

We investigate whether spatial-kinematic substructure in young star-forming regions can be quantified using Moran's $I$ statistic. Its presence in young star clusters would provide an indication that the system formed from initially substructured conditions, as expected by the hierarchical model of star cluster formation, even if the cluster were spatially smooth and centrally concentrated. Its absence, on the other hand, would be evidence that star clusters form monolithically. The Moran's $I$ statistic is applied to $N$-body simulations of star clusters with different primordial spatial-velocity structures, and its evolution over time is studied. It is found that this statistic can be used to reliably quantify spatial-kinematic substructure, and can be used to provide evidence as to whether the spatial-kinematic structure of regions with ages $\lesssim$ 6 Myr is best reproduced by the hierarchical or monolithic models of star formation. Moran's $I$ statistic is also able to conclusively say whether the data are $not$ consistent with initial conditions that lack kinematic substructure, such as the monolithic model, in regions with ages up to, and potentially beyond, 10 Myrs. This can therefore provide a kinematic signature of the star cluster formation process that is observable for many Myr after any initial spatial structure has been erased.

Based on the Sloan Digital Sky Survey Data Release 16, we have detected the large-scale structure of Ly$\alpha$ emission in the Universe at redshifts $z = 2$--3.5 by cross-correlating quasar positions and Ly$\alpha$ emission imprinted in the residual spectra of luminous red galaxies. We apply an analytical model to fit the corresponding Ly$\alpha$ surface brightness profile and multipoles of the redshift-space quasar-Ly$\alpha$ emission cross-correlation function. The model suggests an average cosmic Ly$\alpha$ luminosity density of ${6.6_{-3.1}^{+3.3}}\times 10^{40} {\rm erg\, s^{-1} cMpc^{-3}}$, a $\sim 2\sigma$ detection with a median value about 8--9 times those estimated from deep narrowband surveys of Ly$\alpha$ emitters at similar redshifts. Although the low signal-to-noise ratio prevents us from a significant detection of the Ly$\alpha$ forest-Ly$\alpha$ emission cross-correlation, the measurement is consistent with the prediction of our best-fit model from quasar-Ly$\alpha$ emission cross-correlation within current uncertainties. We rule out the scenario that these Ly$\alpha$ photons mainly originate from quasars. We find that Ly$\alpha$ emission from star-forming galaxies, including contributions from that concentrated around the galaxy centers and that in the diffuse Ly$\alpha$ emitting halos, is able to explain the bulk of the the Ly$\alpha$ luminosity density inferred from our measurements. Ongoing and future surveys can further improve the measurements and advance our understanding of the cosmic Ly$\alpha$ emission field.

LISA will extend the search for gravitational waves (GWs) at $0.1\,{-}\,100$ mHz where loud signals from coalescing binary black holes of $ 10^4 \,{-}\,10^7\,\rm M_{\odot}$ are expected. Depending on their mass and luminosity distance, the uncertainty in the LISA sky-localization decreases from hundreds of deg$^2$ during the inspiral phase to fractions of a deg$^2$ after the merger. By using the semi-analytical model L-Galaxies applied to the Millennium-I merger trees, we generate a simulated Universe to identify the hosts of $z\,{\leq}\,3$ coalescing binaries with total mass of $3\,{\times}\,10^{5}$, $3\,{\times}\,10^6$ and $3\,{\times}\,10^7\rm M_{\odot}$, and varying mass ratio. We find that, even at the time of merger, the number of galaxies around the LISA sources is too large (${\gtrsim}\,10^2$) to allow direct host identification. However, if an X-ray counterpart is associated to the GW sources at $z\,{<}\,1$, all LISA fields at merger are populated by ${\lesssim}\,10$ AGNs emitting above ${\sim}\, 10^{-17} \, \rm erg\,cm^{-2}\,s^{-1}$. For sources at higher redshifts, the poorer sky-localization causes this number to increase up to ${\sim}\, 10^3$. Archival data from eRosita will allow discarding ${\sim}\, 10\%$ of these AGNs, being too shallow to detect the dim X-ray luminosity of the GW sources. Inspiralling binaries in an active phase with masses ${\lesssim}\,10^6\rm M_{\odot}$ at $z\,{\leq}\,0.3$ can be detected, as early as $10$ hours before the merger, by future X-ray observatories in less than a few minutes. For these systems, ${\lesssim}\,10$ AGNs are within the LISA sky-localization area. Finally, the LISA-Taiji network would guarantee the identification of an X-ray counterpart $10$ hours before merger for all binaries at $z\,{\lesssim}\,1$.

HR\,8799 is a young planetary system composed of 4 planets and a double debris belt. Being the first multi-planetary system discovered with the direct imaging technique, it has been observed extensively since 1998. This wide baseline of astrometric measurements, counting over 50 observations in 20 years, permits a detailed orbital and dynamical analysis of the system. To explore the orbital parameters of the planets, their dynamical history, and the planet-to-disk interaction, we made follow-up observations of the system during the VLT/SPHERE GTO program. We obtained 21 observations, most of them in favorable conditions. In addition, we observed HR\,8799 with the instrument LBT/LUCI. All the observations were reduced with state-of-the-art algorithms implemented to apply the spectral and angular differential imaging method. We re-reduced the SPHERE data obtained during the commissioning of the instrument and in 3 open-time programs to have homogeneous astrometry. The precise position of the 4 planets with respect to the host star was calculated by exploiting the fake negative companions method. To improve the orbital fitting, we also took into account all of the astrometric data available in the literature. From the photometric measurements obtained in different wavelengths, we estimated the planets' masses following the evolutionary models. We obtained updated parameters for the orbits with the assumption of coplanarity, relatively small eccentricities, and periods very close to the 2:1 resonance. We also refined the dynamical mass of each planet and the parallax of the system (24.49 $\pm$ 0.07 mas). We also conducted detailed $N$-body simulations indicating possible positions of a~putative fifth innermost planet with a mass below the present detection limits of $\simeq 3$~\MJup.

We present a study of the narrow Fe K$\alpha$ line in seven bright, nearby AGN that have been observed extensively with the Chandra High Energy Transmission Grating (HETG). The HETG data reveal a wider Fe K$\alpha$ line in the first order spectrum than in the second and third order spectra, which we interpret as the result of spatially extended Fe K$\alpha$ emission. We utilize these differences in narrow Fe K$\alpha$ line widths in the multi-order Chandra HETG spectra to determine the spatial extent and intrinsic velocity width of the emitting material in each object. We find that there is modest evidence for spatially extended emission in each object, corresponding to extension of $r\sim5-100$ pc. These distances are significantly larger than those inferred from velocity widths assuming gravitational motions, which give $r\sim0.01-1$ pc. This implies that either the gas is emitting at a range of radii, with smaller radii dominating the velocity width and larger radii dominating the spatial extent, or that the gas is exhibiting non-gravitational motions, which we suggest would be outflows due to slight excess redshift in the line and velocities that exceed the freefall velocity. We also use the spatial extent information to estimate the mass of the emitting gas by counting fluorescing iron atoms, finding masses on the order of $M_\mathrm{gas}\sim10^5-10^8\,M_\odot$. Future work with observatories like XRISM will be able to extend this study to a larger number of AGN and decrease uncertainties that arise due to the low signal-to-noise of the higher order HETG data.

A commonly-used, simplifying assumption when modeling the sources of ultra-high energy cosmic rays (UHECRs) is that all of them accelerate particles to the same maximum energy. Motivated by the fact that candidate astrophysical accelerators exhibit a vast diversity in terms of their relevant properties such as luminosity, Lorentz factor, and magnetic field strength, we study the compatibility of a population of sources with non-identical maximum cosmic-ray energies with the observed energy spectrum and composition of UHECRs at Earth. For this purpose, we compute the UHECR spectrum emerging from a population of sources with a power-law distribution of maximum energies applicable to a broad range of astrophysical scenarios. We find that the allowed source-to-source variance of the maximum energy must be small to describe the data. Even in the most extreme scenario, with a very sharp cutoff of individual source spectra and negative redshift evolution of the accelerators, the maximum energies of 90% of sources must be identical within a factor of three -- in contrast to the variance expected for astrophysical sources.

We present a novel machine learning based approach for detecting galaxy-scale gravitational lenses from interferometric data, specifically those taken with the International LOFAR Telescope (ILT), which is observing the northern radio sky at a frequency of 150 MHz, an angular resolution of 350 mas and a sensitivity of 90 uJy beam-1 (1 sigma). We develop and test several Convolutional Neural Networks to determine the probability and uncertainty of a given sample being classified as a lensed or non-lensed event. By training and testing on a simulated interferometric imaging data set that includes realistic lensed and non-lensed radio sources, we find that it is possible to recover 95.3 per cent of the lensed samples (true positive rate), with a contamination of just 0.008 per cent from non-lensed samples (false positive rate). Taking the expected lensing probability into account results in a predicted sample purity for lensed events of 92.2 per cent. We find that the network structure is most robust when the maximum image separation between the lensed images is greater than 3 times the synthesized beam size, and the lensed images have a total flux density that is equivalent to at least a 20 sigma (point-source) detection. For the ILT, this corresponds to a lens sample with Einstein radii greater than 0.5 arcsec and a radio source population with 150 MHz flux densities more than 2 mJy. By applying these criteria and our lens detection algorithm we expect to discover the vast majority of galaxy-scale gravitational lens systems contained within the LOFAR Two Metre Sky Survey.

Blue Compact Dwarfs (BCDs) are low-luminosity (M$_{K} > -21$ mag), metal-poor ($\frac{1}{50}$ $\le Z/Z_{\odot} \le\frac{1}{2}$), centrally concentrated galaxies with bright clumps of star-formation. Cosmological surface brightness dimming and small size limit their detection at high redshifts, making their formation process difficult to observe. Observations of BCDs are needed at intermediate redshifts, where they are still young enough to show their formative stages, particularly in the outer regions where cosmic gas accretion should drive evolution. Here, we report the discovery of excess far-ultraviolet (FUV) emission in the outer regions of 11 BCDs in the GOODS-South field at redshifts between 0.1 and 0.24, corresponding to look back times of 1.3 - 2.8 Gyr in standard cosmology. These observations were made by the Ultra-Violet Imaging Telescope (UVIT) on AstroSat. For ten BCDs, the radial profiles of intrinsic FUV emission, corrected for the instrument point spread function, have larger scale-lengths than their optical counterparts observed with the Hubble Space Telescope. Such shallow FUV profiles suggest extended star-formation in cosmically accreting disks. Clumpy structure in the FUV also suggests the outer FUV disks are gravitationally unstable. Dynamical friction on the clumps drives them inward at an average rate exceeding $10^6~M_{\odot}$Gyr$^{-1}$.

Understanding the chemical processes during starless core and prestellar core evolution is an important step in understanding the initial stages of star and disk formation. This project is a study of deuterated ammonia, o-NH$_2$D, in the L1251 star-forming region toward Cepheus. Twenty-two dense cores (twenty of which are starless or prestellar, and two of which have a protostar), previously identified by p-NH$_3$ (1,1) observations, were targeted with the 12m Arizona Radio Observatory telescope on Kitt Peak. o-NH$_2$D J$_{\rm{K_a} \rm{K_c}}^{\pm} =$ $1_{11}^{+} \rightarrow 1_{01}^{-}$ was detected in 13 (59\%) of the NH$_3$-detected cores with a median sensitivity of $\sigma_{T_{mb}} = 17$ mK. All cores detected in o-NH$_2$D at this sensitivity have p-NH$_3$ column densities $> 10^{14}$ cm$^{-2}$. The o-NH$_2$D column densities were calculated using the constant excitation temperature (CTEX) approximation while correcting for the filling fraction of the NH$_3$ source size. The median deuterium fraction was found to be 0.11 (including 3$\sigma$ upper limits). However, there are no strong, discernible trends in plots of deuterium fraction with any physical or evolutionary variables. If the cores in L1251 have similar initial chemical conditions, then this result is evidence of the cores physically evolving at different rates.

Galaxies and their central supermassive black holes are known to coevolve, but the physical background for this is unknown as of yet. The High-$z$ Universe probed via Lensing by QSOs (HULQ) project aims to investigate this coevolution by using quasi-stellar object (QSO) host galaxies acting as gravitational lenses (QSO lenses). We present the results of the spectroscopic observation of the first QSO lens candidate from the HULQ project, HULQ J0002+0239, which consists of a QSO host galaxy at $z_{\rm d} = 1.455$ and four seemingly lensed objects in a cross-like configuration. Deep optical spectra of two of the possibly lensed objects with $z \sim 24.5$ mag were obtained with the Gemini Multi-Object Spectrograph on the Gemini North Telescope. Their spectra reveal that the objects are newly discovered galaxies at $z=0.29$ and $z=1.11$, and we conclude that HULQ J0002+0239 is not a QSO lens. Our QSO lens search results are so far in agreement with the predicted number of QSO lenses, and we discuss how the future investigation of additional QSO lens candidates could tell us more about the evolution of the black hole mass and host galaxy scaling relations.

Two independent groups reported the discovery of an isolated dark stellar remnant in the microlensing event OGLE-2011-BLG-0462 based on photometric ground-based observations coupled with astrometric measurements taken with the Hubble Space Telescope. These two analyses yielded discrepant mass measurements, with the first group reporting that the lensing object is a black hole of 7.1 +/- 1.3 solar masses whereas the other concluded that the microlensing event was caused by either a neutron star or a low-mass black hole (1.6-4.4 solar masses). Here, we scrutinize the available photometric and astrometric data and conclude that systematic errors are a cause of the discrepant measurements. We find that the lens is an isolated black hole with a mass of 7.88 +/- 0.82 solar masses located at a distance of 1.49 +/- 0.12 kpc. We also study the impact of blending on the accuracy of astrometric microlensing measurements. We find that low-level blending by source companions is a major, previously unrecognized, challenge to astrometric microlensing measurements of black hole masses.

[Abridged] Far-infrared observations with Herschel revealed a surprisingly low abundance of cold-water reservoirs in protoplanetary discs. On the other hand, a handful of discs show emission of hot water transitions excited at temperatures above a few hundred Kelvin. In particular, the protoplanetary discs around the Herbig Ae stars HD 100546 and HD 163296 show opposite trends in terms of cold versus hot water emission: in the first case, the ground-state transitions are detected and the high-J lines are undetected, while the trend is opposite in HD 163296. We performed a spectral analysis using the thermo-chemical model DALI. We find that HD 163296 is characterised by a water-rich (abundance $\gtrsim 10^{-5}$) hot inner disc (within the snowline) and a water-poor ($< 10^{-10}$) outer disc: the relative abundance may be due to the thermal desorption of icy grains that have migrated inward. Remarkably, the size of the H$_2$O emitting region corresponds to a narrow dust gap visible in the millimeter continuum at $r=10\,$au with ALMA. The low-J lines detected in HD 100546 instead imply an abundance of a few $10^{-9}$ in the cold outer disc ($> 40$ au). The emitting region of the cold H$_2$O transitions is spatially coincident with that of the H$_2$O ice previously seen in the near-infrared. Notably, millimetre observations with ALMA reveal the presence of a large dust gap between nearly 40 and 150 au, likely opened by a massive embedded protoplanet. In both discs, we find that the warm molecular layer in the outer region (beyond the snow line) is highly depleted of water molecules, implying an oxygen-poor chemical composition of the gas. We speculate that gas-phase oxygen in the outer disc is readily depleted and its distribution in the disc is tightly coupled to the dynamics of the dust grains.

We report our search for quasi-periodic signals in long-term optical and $\gamma$-ray data for the blazar PKS~1222+216, where the data are from the Steward Observatory blazar monitoring program and the all-sky survey with the Large Area Telescope onboard the {\it Fermi Gamma-ray Space Telescope}, respectively. A quasi-periodic signal, with a period of $\simeq$420\,days and a significance of $>5\sigma$, is found in the measurements of the optical linear polarization degree for the source, while no similar signals are found in the optical and $\gamma$-ray light curves covering approximately the same time period of $\sim$10\,yr. We study the quasi-periodic variations by applying a helical jet model and find that the model can provide a good explanation. This work shows that polarimetry can be a powerful tool for revealing the physical properties, in particular the configuration of the magnetic fields of jets from galactic supermassive black holes.

Pulsar wind nebulae (PWNe) are one of the most energetic galactic sources with bright emissions from radio waves to very high-energy gamma-rays. We perform wideband X-ray spectroscopy of four energetic PWNe, N157B, PSR J1813-1749, PSR J1400-6325, and G21.5-0.9, with the Suzaku, Chandra, NuSTAR, and Hitomi observatories. A significant spectral break or cutoff feature is found in the hard X-ray band for all the samples, except for N157B. The break energies in the broken power-law fitting are in the range of 4--14 keV, whereas the cutoff energies in the cutoff power-law fitting are at 22 keV or higher. The break or cutoff energy does not show a significant correlation with either the spin-down energy or characteristic age of the hosting pulsars. A possible correlation is found between the photon index change in the broken power-law fitting and the X-ray emitting efficiency of the pulsars, although its significance is not high enough to be conclusive. We discuss what determines the break parameters based on simple models.

The European Solar Telescope (EST) is a project aimed at studying the magnetic connectivity of the solar atmosphere, from the deep photosphere to the upper chromosphere. Its design combines the knowledge and expertise gathered by the European solar physics community during the construction and operation of state-of-the-art solar telescopes operating in visible and near-infrared wavelengths: the Swedish 1m Solar Telescope (SST), the German Vacuum Tower Telescope (VTT) and GREGOR, the French T\'elescope H\'eliographique pour l'\'Etude du Magn\'etisme et des Instabilit\'es Solaires (TH\'EMIS), and the Dutch Open Telescope (DOT). With its 4.2 m primary mirror and an open configuration, EST will become the most powerful European ground-based facility to study the Sun in the coming decades in the visible and near-infrared bands. EST uses the most innovative technological advances: the first adaptive secondary mirror ever used in a solar telescope, a complex multi-conjugate adaptive optics with deformable mirrors that form part of the optical design in a natural way, a polarimetrically compensated telescope design that eliminates the complex temporal variation and wavelength dependence of the telescope Mueller matrix, and an instrument suite containing several (etalon-based) tunable imaging spectropolarimeters and several integral field unit spectropolarimeters. This publication summarises some fundamental science questions that can be addressed with the telescope, together with a complete description of its major subsystems.

With its exquisite sensitivity, wavelength coverage, and spatial and spectral resolution, the James Webb Space Telescope is poised to revolutionise our view of the distant, high-redshift ($z>5$) Universe. While Webb's spectroscopic observations will be transformative for the field, photometric observations play a key role in identifying distant objects and providing more comprehensive samples than accessible to spectroscopy alone. In addition to identifying objects, photometric observations can also be used to infer physical properties and thus be used to constrain galaxy formation models. However, inferred physical properties from broadband photometric observations, particularly in the absence of spectroscopic redshifts, often have large uncertainties. With the development of new tools for forward modelling simulations it is now routinely possible to predict observational quantities, enabling a direct comparison with observations. With this in mind, in this work, we make predictions for the colour evolution of galaxies at $z=5-15$ using the FLARES: First Light And Reionisation Epoch Simulations cosmological hydrodynamical simulation suite. We predict a complex evolution, driven predominantly by strong nebular line emission passing through individual bands. These predictions are in good agreement with existing constraints from Hubble and Spitzer as well as some of the first results from Webb. We also contrast our predictions with other models in the literature: while the general trends are similar we find key differences, particularly in the strength of features associated with strong nebular line emission. This suggests photometric observations alone should provide useful discriminating power between different models.

Young massive stellar clusters are extreme environments and potentially provide the means for efficient particle acceleration. Indeed, they are increasingly considered as being responsible for a significant fraction of cosmic rays (CRs) accelerated within the Milky Way. Westerlund 1, the most massive known young stellar cluster in our Galaxy is a prime candidate for studying this hypothesis. While the very-high-energy $\gamma$-ray source HESS J1646-458 has been detected in the vicinity of Westerlund 1 in the past, its association could not be firmly identified. We aim to identify the physical processes responsible for the $\gamma$-ray emission around Westerlund 1 and thus to better understand the role of massive stellar clusters in the acceleration of Galactic CRs. Using 164 hours of data recorded with the High Energy Stereoscopic System (H.E.S.S.), we carried out a deep spectromorphological study of the $\gamma$-ray emission of HESS J1646-458. We furthermore employed H I and CO observations of the region to infer the presence of gas that could serve as target material for interactions of accelerated CRs. We detected large-scale ($\sim 2^\circ$ diameter) $\gamma$-ray emission with a complex morphology, exhibiting a shell-like structure and showing no significant variation with $\gamma$-ray energy. The combined energy spectrum of the emission extends to several tens of TeV, and is uniform across the entire source region. We did not find a clear correlation of the $\gamma$-ray emission with gas clouds as identified through H I and CO observations. We conclude that, of the known objects within the region, only Westerlund 1 can explain the bulk of the $\gamma$-ray emission. Several CR acceleration sites and mechanisms are conceivable, and discussed in detail. (abridged)

The recently released Type Ia supernovae (SNe Ia) sample, Pantheon+, is an updated version of Pantheon and has very important cosmological implications. To explore the origin of the enhanced constraining power and internal correlations of datasets in different redshifts, we perform a comprehensively tomographic analysis of the Pantheon+ sample. Using the Pantheon+ data alone, we give the $2\,\sigma$ lower bound on the Hubble constant $H_0>45.7$ km s$^{-1}$ Mpc$^{-1}$ and the matter fraction $\Omega_m=0.367\pm0.030$, which shows the evidence of dark energy at the $21\,\sigma$ confidence level but is in a $1.7\,\sigma$ tension with that from the Planck-2018 measurement. Combining the Pantheon+ sample with cosmic microwave background, baryon acoustic oscillations, cosmic chronometers, galaxy clustering and weak lensing data, we give the strongest constraint $H_0=67.88\pm0.42$ km s$^{-1}$ Mpc$^{-1}$ at the $1\,\sigma$ confidence level. After dividing the full sample to 10 bins, we find that the first bin in the redshift range $z\in[0.00122, \, 0.227235]$ dominates the constraining power of the whole sample. We also investigate the effects of low-z and high-z subsamples of Pantheon+ on $H_0$ and $\Omega_{m}$, and find that low-z SNe Ia do not have enough constraining power until $z\sim0.1$. Interestingly, high-z SNe Ia in the redshift range $z>0.227235$ can give an competitive constraint on $\Lambda$CDM when compared to three low-z bins. We expect that future high-precision SNe Ia data can independently determine both $H_0$ and $\Omega_{m}$.

The VRI light curves were measured for the low-mass eclipsing binary V608 Cam as a part of our long-term observational project for studying of eclipsing binaries with a short orbital period. The Tess light curve solution in Phoebe results to the detached configuration, where the temperature of primary component was fixed to $T_1 = 5300$ K according to Gaia results, which gives us $T_2 = 4110 \pm 50$ K for the secondary. The spectral type of the primary component was derived to be K0 and the photometric mass ratio was estimated $q = 0.92 \pm 0.07$. Characteristics and temporal variation of the cold region on the surface of the secondary component were estimated and are attributed to apparent period changes of this eclipsing binary with a cycle of about 2.4 yr.

Context: While star formation on large molecular cloud scales and on small core and disk scales has been investigated intensely over the past decades, the connection of the large-scale interstellar material with the densest small-scale cores has been a largely neglected field. Methods: Using NOEMA and the IRAM 30\,m telescope, we mapped large areas (640\,arcmin$^2$) of the archetypical star formation complex Cygnus X at 3.6\,mm wavelengths in line and continuum emission. Results: The scope and outline of The Cygnus Allscale Survey of Chemistry and Dynamical Environments (CASCADE) is presented. We then focus on the first observed subregion in Cygnus X, namely the DR20 star formation site, which comprises sources in a range of evolutionary stages from cold pristine gas clumps to more evolved ultracompact H{\sc ii} regions. The data covering cloud to cores scales at a linear spatial resolution of $<5000$\,au reveal several kinematic cloud components that are likely part of several large-scale flows onto the central cores. The temperature structure of the region is investigated by means of the HCN/HNC intensity ratio and compared to dust-derived temperatures. We find that the deuterated DCO$^+$ emission is almost exclusively located toward regions at low temperatures below 20\,K. Investigating the slopes of spatial power spectra of dense gas tracer intensity distributions (HCO$^+$, H$^{13}$CO$^+$, and N$_2$H$^+$), we find comparatively flat slopes between $-2.9$ and $-2.6$, consistent with high Mach numbers and/or active star formation in DR20.

Globular clusters (GCs) are considered strong candidates for hosting rogue (free-floating) planets. Since they are not bound to a star, they are undetectable by any traditional detection methods: transit, radial velocity, or direct imaging. Gravitational microlensing (ML), which causes transient brightening of background stars by passing foreground masses, is, on the other hand, an established method of detecting planets and proves promising for application in GCs. By employing the image subtraction technique, differential photometry on the time-series images of GCs could extract variability events, build light curves and inspect them for the presence of microlensing. However, instrumental anomalies and varying observing conditions over a long observational campaign period result in the distortion of stellar Point Spread Function (PSF), which affects the subtraction quality and leads to false-positive transient detection and large-scale noise structure in the subtracted images. We propose an iterative image reconstruction method as a modification to the Scaled Gradient Projection (SGP) algorithm, called the Flux-Conserving Scaled Gradient Projection (FC-SGP), to restore the shapes of stars while preserving their flux well within the photometrically accepted tolerance. We perform an extensive empirical comparative study of FC-SGP with different image restoration algorithms like the Richardson-Lucy (RL) and the original SGP algorithms, using several physically motivated metrics and experimental convergence analysis. We find that FC-SGP could be a promising approach for astronomical image restoration. In the future, we aim to extend its application to different image formats while maintaining the performance of the proposed algorithm.

Results of 111-MHz monitoring observations carried out on the Big Scanning Antenna of the Pushchino Radio Astronomy Observatory during September 1-28, 2015 are presented. Fifty-four pulsating sources were detected at declinations $-9^o < \delta < +42^o$. Forty-seven of these are known pulsars, five are new sources, and two are previously discovered transients. Estimates of the peak flux densities and dispersion measures are presented or all these sources.

In this paper, we study stellar light curves from the TESS satellite (Transiting Exoplanet Survey Satellite) for the presence of stellar flares. The main aim is to detect stellar flares using two-minutes cadence data and to perform statistical analysis. To find and analyze stellar flares we prepared automatic software WARPFINDER. We implemented three methods described in this paper: trend, difference, and profile fitting. Automated search for flares was accompanied by visual inspection. Using our software we analyzed two-minute cadence light curves of 330,000 stars located in the first 39 sectors of TESS observations. As a result, we detected over 25,000 stars showing flare activity with the total number of more than 140,000 flares. This means that about 7.7% of all the analyzed objects are flaring stars. The estimated flare energies range between $10^{31}$ and $10^{36}$ erg. We prepared a preliminary preview of the statistical distribution of parameters such as a flare duration, amplitudes and energy, and compared it with previous results. The relationship between stellar activity and its spectral type, temperature and mass was also statistically analyzed. Based on the scaling laws, we estimated the average values of the magnetic field strength and length of the flare loops. In our work, we used both single (about 60%), and double (about 40%) flare profiles to fit the observational data. The components of the double profile are supposed to be related to the direct heating of the photosphere by non-thermal electrons and back warming processes.

Abell 1033 is a merging galaxy cluster of moderate mass (M500 = 3.2e14 Msun) which hosts a broad variety of diffuse radio sources linked to different astrophysical phenomena. The most peculiar one is an ultra-steep spectrum elongated feature which is the prototype of the category of gently re-energized tails (GReET). In addition, the cluster hosts sources previously classified as a radio phoenix and a radio halo. In this work, we aim to improve the understanding of the cosmic ray acceleration mechanisms in galaxy clusters in a frequency and mass range poorly explored so far. To investigate the ultra-steep synchrotron emission in the cluster, we perform a full calibration of a LOFAR observation centered at 54 MHz. We analyze this observation together with re-calibrated data of the LOFAR Two-meter Sky Survey at 144 MHz and an archival GMRT observation at 323 MHz. We perform a spectral study of the radio galaxy tail connected to the GReET to test if the current interpretation of the source is in agreement with observational evidence below 100 MHz. Additionally, we study the radio halo at different frequencies. We report an extreme spectral curvature for the GReET, the spectral index flattens from $\alpha_{144}^{323} = -4$ to $\alpha_{54}^{144} = -2$. This indicates the presence of a cut-off in the electron energy spectrum. At the cluster center, we detect the radio halo at 54, 144 and at lower significance at 323 MHz. We categorize it as an ultra-steep spectrum radio halo with a spectral index $\alpha = -1.65 \pm 0.17$. Additionally, it is found to be significantly above the radio power-to-cluster mass correlations reported in the literature. Furthermore, the synchrotron spectrum of the halo is found to further steepen between 144 and 323 MHz, in agreement with the presence of a break in the electron spectrum, which is a prediction of homogeneous re-acceleration models.

The gravitational potential of initially Poisson distributed primordial black holes (PBH) can induce a stochastic gravitational-wave background (SGWB) at second order in cosmological perturbation theory. This SGWB was previously studied in the context of general relativity (GR) and modified gravity setups by assuming a monochromatic PBH mass function. Here we extend the previous analysis in the context of GR by studying the aforementioned SGWB within more physically realistic regimes where PBHs have different masses. In particular, starting from a power-law cosmologically motivated primordial curvature power spectrum we extract the extended PBH mass function and the associated to it PBH gravitational potential which acts as the source of the scalar induced SGWB. At the end, by taking into account the dynamical evolution of the PBH gravitational potential during the transition from the matter era driven by PBHs to the radiation era we extract the respective GW signal today. Interestingly, in order to trigger an early PBH-dominated era and avoid the GW constraints at BBN we find that the spectral index $n_\mathrm{s}$ of our primordial curvature power spectrum should be within the narrow range $n_\mathrm{s}\in[1.1482,1,1493]$ while at the same time the GW signal is found to be potentially detectable by LISA.

Despite image reconstruction becoming more widespread when interpreting OIFITS Data, model fitting in u,v space often remains the best way to interpret data, either because of the sparsity of the data, or because a quantitative measurement needs to be done. PMOIRED, is a flexible Python library to visualize, manipulate and model OIFITS data using simple geometric models. The strength of PMOIRED resides in its capability to combine linearly various simple components to create complex scenes, while linking, constraining, and adding priors to fitted parameters. The code also enables grid search to find global minima, as well as data resampling to better evaluate uncertainties. In addition to analytical functions, arbitrary radial profiles, azimuthal variations or sparse wavelet modelling of spectra are implemented.

HD 129929 is a slowly-rotating $\beta$ Cephei pulsator with a rich spectrum of detected oscillations, including two rotational multiplets. The asteroseismic interpretation revealed the presence of radial differential rotation in this massive star of $\sim$9.35 M . The stellar core is indeed estimated to spin $\sim$3.6 times faster than the surface. The surface rotation was consequently derived as $\sim$2 km/s. This massive star represents an ideal counter-part to the wealth of space-based photometry results for main-sequence and evolved low-mass stars. Those latter have revealed a new, and often unexpected, picture of the angular momentum transport processes acting in stellar interiors. We investigate in a new way the constraints on the internal rotation of HD 129929, focusing on their interpretation for the evolution of the internal rotation during the main sequence of a massive star. We test separately hydrodynamic and magnetic instability transport processes of angular momentum. We used the best asteroseismic model obtained in an earlier work. We calibrated stellar models including rotation, with different transport processes, to reproduce that reference model. We then looked whether one process is favoured to reproduce the rotation profile of HD 129929, based on the fit of the asteroseismic multiplets. The impact of the Tayler magnetic instability on the angular momentum transport predicts a ratio of the core-to-surface rotation rate of only 1.6, while the recently revised prescription of this mechanism predicts solid-body rotation. Both are too low in comparison with the asteroseismic inference. The models with only hydrodynamic processes are in good agreement with the asteroseismic measurements. Strikingly, we can also get a constraint on the profile of rotation on the zero age main sequence: likely, the ratio between the core and surface rotation was at least $\sim$1.7.

Context: Broadband optical constants of astrophysical ice analogues in the infrared (IR) and terahertz (THz) ranges are required for modeling the dust continuum emission and radiative transfer in dense and cold regions, where thick icy mantles are formed on the surface of dust grains. Aims: In this paper, the THz time-domain spectroscopy (TDS) and the Fourier-transform IR spectroscopy (FTIR) are combined to study optical constants of CO and CO$_2$ ices in the broad THz-IR spectral range. Methods: The measured ices are grown at cryogenic temperatures by gas deposition on a cold Si window. A method to quantify the broadband THz-IR optical constants of ices is developed based on the direct reconstruction of the complex refractive index of ices in the THz range from the TDS data, and the use of the Kramers-Kronig relation in the IR range for the reconstruction from the FTIR data. Uncertainties of the Kramers-Kronig relation are eliminated by merging the THz and IR spectra. The reconstructed THz-IR response is then analyzed using classical models of complex dielectric permittivity. Results: The complex refractive index of CO and CO$_2$ ices deposited at the temperature of $28$ K is obtained in the range of 0.3--12.0 THz. Based on the measured dielectric constants, opacities of the astrophysical dust with CO and CO$_2$ icy mantles are computed. Conclusions: The developed method can be used for a model-independent reconstruction of optical constants of various astrophysical ice analogs in a broad THz-IR range. Such data can provide important benchmarks to interpret the broadband observations from the existing and future ground-based facilities and space telescopes. The reported results will be useful to model sources that show a drastic molecular freeze-out, such as central regions of prestellar cores and mid-planes of protoplanetary disks, as well as CO and CO$_2$ snow lines in disks.

Open clusters serve as a useful tool for calibrating models of the relationship between mass, rotation, and age for stars with an outer convection zone due to the homogeneity of the stars within the cluster. Cluster to cluster comparisons are essential to determine whether the universality of spin down relations holds. NGC 6709 is selected as a young open cluster for which no rotation periods of members have previously been obtained. This cluster is at a distance of over 1 kpc and has two red giant members. Isochrones place the age of the cluster at around 150 Myr, or approximately the same age as the Pleiades. Photometry is obtained over a multi-month observing season at the robotic observatory STELLA. After basic processing, PSF photometry was derived using Daophot II, and a suite of related software allowed us to create time series of relative magnitude changes for each star. Four time series analysis methods are then applied to these light curves to obtain rotation periods for members stars. We obtain for the first time rotation periods for 45 FGK cluster members of NGC 6709. We compare our rotation periods to Gaia EDR3 colors and find a slow-rotating sequence with increasing rotation periods towards redder stars and a smaller clump of rapid rotators that have not yet joined this sequence. NGC 6709 has rotation periods very similar to that of another Pleiades-age open cluster, NGC 2516.

Gamma-ray bursts (GRBs) are known to have the most relativistic jets, with initial Lorentz factors in the order of a few hundreds. Many GRBs display an early X-ray light-curve plateau, which was not theoretically expected and therefore puzzled the community for many years. Here, we show that this observed signal is naturally obtained within the classical GRB "fireball" model, provided that the initial Lorentz factor is rather a few tens, and the expansion occurs into a medium-low density "wind". The range of Lorentz factors in GRB jets is thus much wider than previously thought and bridges an observational gap between mildly relativistic jets inferred in active galactic nuclei, to highly relativistic jets deduced in few extreme GRBs. Furthermore, long GRB progenitors are either not Wolf-Rayet stars, or the wind properties during the final stellar evolution phase are different than at earlier times. We discuss several testable predictions of this model.

Mass discrepancy in galaxies invokes dark matter, or alternatively modification of gravity or inertia. These theoretical possibilities may be distinguished by the statistical relation between the centripetal acceleration of particles in orbital motion and the expected Newtonian acceleration for the observed distribution of baryons in galaxies. Here predictions of cold dark matter halos, modified gravity, and modified inertia are compared and tested by a statistical sample of rotation curves of galaxies. Modified gravity under an estimated mean external field correctly predicts the observed statistical relation of accelerations. Cold dark matter halos predict systematically deviating relations and modified inertia is inconsistent with the apparently seen difference between the inner and outer parts. All aspects of rotation curves are most naturally explained by modified gravity.

Characterization of directly imaged exoplanets is one of the most eagerly anticipated science functions of the James Webb Space Telescope. MIRI, the mid-IR instrument has the capability to provide unique spatially resolved photometric data points in a spectral range never achieved so far for such objects. We aim to present the very first on-sky contrast measurements of the MIRI's coronagraphs. In addition to a classical Lyot coronagraph at the longest wavelength, this observing mode implements the concept of the four quadrant phase mask for the very first time in a space telescope. We observed single stars together with a series of reference stars to measure raw contrasts as they are delivered on the detector, as well as reference subtracted contrasts. MIRI's coronagraphs achieve raw contrasts greater than $10^3$ at the smallest angular separations (within $1''$) and about $10^5$ further out (beyond $5\sim6''$). Subtracting the residual diffracted light left unattenuated by the coronagraph has the potential to bring the final contrast down to the background and detector limited noise floor at most angular separations (a few times $10^4$ at less than $1''$). MIRI coronagraphs behave as expected from simulations. In particular the raw contrasts for all four coronagraphs are fully consistent with the diffractive model. Contrasts obtained with subtracting reference stars also meet expectations and are fully demonstrated for two four quadrant phase masks (F1065C and F1140C). The worst contrast, measured at F1550C, is very likely due to a variation of the phase aberrations at the primary mirror during the observations, and not an issue of the coronagraph itself. We did not perform reference star subtraction with the Lyot mask at F2300C, but we anticipate that it would bring the contrast down to the noise floor.

We report the results from new observations from a long-term radial velocity monitoring campaign complemented with high resolution spectroscopy, as well as new astrometry and seismology of a sample of 41 red giants from the third version of APOKASC, which includes young alpha rich (YAR) stars. The aim is to better characterize the YAR stars in terms of binarity fraction, mass, abundance trends and kinematic properties. The radial velocities of HERMES, APOGEE and Gaia were combined to determine the binary fraction among YAR stars. In combination with their mass estimate, their evolutionary status, chemical composition and kinematic properties, it allows to better constrain the nature of these objects. We find that the frequency of binaries among over-massive stars is not significantly different than that of the other stars in our sample, but that the most massive YAR stars are indeed single, which has been predicted by population synthesis models. Studying their [C/N], [C/Fe] and [N/Fe] trends with mass, many over-massive stars do not follow the APOKASC stars, favouring the scenario that most of them are product of mass transfer. Our sample further includes two under-massive stars, with sufficiently low masses so that these stars could not have reached the red giant phase without significant mass loss. Both over-massive and under-massive stars might show some anomalous APOGEE abundances such as N, Na, P, K and Cr, although higher resolution optical spectroscopy might be needed to confirm these findings. Considering the significant fraction of stars that are formed in pairs and the variety of ways that make mass transfer possible, the diversity in properties in terms of binarity and chemistry of the over-massive and under-massive stars studied here implies that it is not safe to directly relate the mass of the YAR stars with age and that most of these objects are likely not young.

We present extended Lyman-{\alpha} (Ly{\alpha}) emission out to 800 kpc of 1034 [O III]-selected galaxies at redshifts 1.9<z<2.35 using the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX). The locations and redshifts of the galaxies are taken from the 3D-HST survey. The median-stacked surface brightness profile of Ly{\alpha} emission of the [O III]-selected galaxies agrees well with that of 968 bright Ly{\alpha}-emitting galaxies (LAEs) at r>40 kpc from the galaxy centers. The surface brightness in the inner parts (r<10 kpc) around the [O III]-selected galaxies, however, is ten times fainter than that of the LAEs. Our results are consistent with the notion that photons dominating the outer regions of the Ly{\alpha} halos are not produced in the central galaxies but originate outside of them.

Building on previous Bayesian approaches, we introduce a novel formulation of probabilistic cross-identification, where detections are directly associated to (hypothesized) astronomical objects in a globally optimal way. We show that this new method scales better for processing multiple catalogs than enumerating all possible candidates, especially in the limit of crowded fields, which is the most challenging observational regime for new-generation astronomy experiments such as the Rubin Observatory Legacy Survey of Space and Time (LSST). Here we study simulated catalogs where the ground-truth is known and report on the statistical and computational performance of the method. The paper is accompanied by a public software tool to perform globally optimal catalog matching based on directional data.

We present the first James Webb Space Telescope/NIRCam-led determination of $7<z<9$ galaxy properties based on broadband imaging from 0.8 to 5 microns as part of the GLASS-JWST Early Release Science program. This is the deepest dataset acquired at these wavelengths to date, with an angular resolution $\lesssim0.14$ arcsec. We robustly identify 14 galaxies with S/N>8 in F444W from 8 arcmin$^2$ of data at $m_{AB}\leq 28$ from a combination of dropout and photometric redshift selection. From simulated data modeling, we estimate the dropout sample purity to be $\gtrsim90\%$. We find that the number density of these sources is broadly consistent with expectations from the UV luminosity function determined from Hubble Space Telescope data. We characterize galaxy physical properties using a Bayesian Spectral Energy Distribution fitting method, finding median stellar mass $10^{8.7}M_\odot$ and age 130 Myr, indicating they started ionizing their surroundings at redshift $z>9.5$. Their star formation main sequence is consistent with predictions from simulations. Lastly, we introduce an analytical framework to constrain main-sequence evolution at $z>7$ based on galaxy ages and basic assumptions, through which we find results consistent with expectations from cosmological simulations. While this work only gives a glimpse of the properties of typical galaxies that are thought to drive the reionization of the universe, it clearly shows the potential of JWST to unveil unprecedented details on galaxy formation in the first billion years.

Alpha$^2$ Canum Venaticorum (AM CVn) is a strongly magnetic star with peculiar chemical signatures and periodic variability that have been long attributed to the diffusion of magnetic elements through the photosphere, leading to chemical spots across the stellar surface. However, recent studies of other magnetic hot stars are consistent with magnetospheric clouds above the surface. Here we take a renewed approach to modeling AM CVn with a simplified dynamical magnetosphere (DM) and a tilted, offset magnetic dipole to reproduce its Transiting Exoplanet Survey Satellite (TESS) variability. Our dipole model also reproduces well the magnetic surface map of AM CVn from Silvester, Kochukhov, & Wade (2014). Its ultraviolet variability, from IUE archival spectra, is also consistent with traditional reddening models. We further discuss the implications of a magnetosphere on other observable quantities from the system to conclude that it is unlikely to be present in AM CVn.

The discovery of extended gamma-ray emission toward a number of middle-aged pulsars suggests the possibility of long-lived particle confinement beyond the classical pulsar wind nebula (PWN) stage. How this emerging source class can be extrapolated to a Galactic population remains unclear. We aim to evaluate how pulsar halos fit in existing TeV observations, under the assumption that all middle-aged pulsars develop halos similar to those observed toward the J0633+1746 or B0656+14 pulsars. We modeled the populations of supernova remnants, PWNe, and pulsar halos in the Milky Way. The PWN-halo evolutionary sequence is described in a simple yet coherent framework, and both kinds of objects are assumed to share the same particle injection properties. We then assessed the contribution of the different source classes to the very-high-energy emission from the Galaxy. The synthetic population can be made consistent with the flux distribution of all known objects, including unidentified objects, for a reasonable set of parameters. The fraction of the populations predicted to be detectable in surveys of the Galactic plane with HESS. and HAWC is then found to be in good agreement with their actual outcome, with a number of detectable halos ranging from 30 to 80% of the number of detectable PWNe. Prospects for CTA involve the detection of 250-300 sources in the Galactic Plane Survey, including 170 PWNe and up to 100 halos. The extent of diffusion suppression in halos has a limited impact on such prospects but its magnitude has a strong influence. The level of diffuse emission from unresolved populations in each survey is found to be dominated by halos and comparable to large-scale interstellar radiation powered by cosmic rays above 0.1-1TeV. Pulsar halos are shown to be viable counterparts to a fraction of the currently unidentified sources if they develop around most middle-aged pulsars (abridged).

We present a reduction and analysis of the \textit{James Webb Space Telescope} (JWST) SMACS~0723 field to conduct a search for ultra-high-redshift galaxies ($9<z<12$) present within the Epoch of Reionization. We use a combination of photometric redshifts and spectral energy distribution (SED) modelling-based selection criteria to optimise sample completeness while minimising contamination. We find four $z > 9$ candidate galaxies which have not previously been identified, with one object at $z = 11.5$, and another which is possibly a close pair of galaxies. These sources are fairly bright with $m_{277} \sim 26 - 28$. A significant fraction of these sources show evidence for Balmer-breaks or extreme emission lines from H$\beta$ and [OIII], demonstrating that the stellar populations could be quite advanced in age or very young depending on the exact cause of the F444W excess. We discuss the stellar masses and resolved structures of these early galaxies and find that the S\'{e}rsic indices reveal a mixture of light concentration levels, but that the sizes of all our systems are exceptionally small in general with $< 0.5$~kpc. These systems have stellar masses M$_{*} \sim 10^{9.5}$ M$_{\odot}$ with our $z \sim 11.5$ candidate a dwarf galaxy with a stellar mass M$_{*} \sim 10^{7.9}$ M$_{\odot}$ These candidate ultra high-redshift galaxies are excellent targets for future NIRSpec observations to understand better their physical nature.

Theories of dark energy that affect the speed of gravitational waves $c_{\rm GW}$ on cosmological scales naturally lead to a frequency-dependent transition of that speed close to the LIGO/Virgo/KAGRA (LVK) band. While observations such as GW170817 assure us that $c_{\rm GW}$ is extremely close to the speed of light in the LVK band, a frequency-dependent transition below the LVK band is a smoking-gun signal for large classes of dynamical dark energy theories. Here we discuss 1) how the remnants of such a transition can be constrained with observations in the LVK band, 2) what signatures are associated with such a transition in the LISA band, and 3) how joint observations in the LVK and LISA bands allow us to place tight constraints on this transition and the underlying theories. We find that deviations of $c_{\rm GW}$ can be constrained down to a level of $\sim 10^{-17}$ in the LVK ${\textit and}$ LISA bands even for mild frequency-dependence, much stronger than existing bounds for frequency-independent $c_{\rm GW} \neq c$. We use the strain data from GW170817 to bound the deviation of $c_{\rm GW}$ to be less than $10^{-17}$ at 100 Hz and less than $10^{-18}$ at 500 Hz. We also identify a particularly interesting type of transition in between the LVK and LISA bands and show how multi-band observations can constrain this further. Finally, we discuss what these current and forecasted constraints imply for the underlying dark energy theories.

We report on a study of the natural warm inflationary paradigm (WNI). We show two important new results arise in this model. One is that the observational constraints on the primordial power spectrum from the cosmic microwave background (CMB) can be satisfied without going beyond the Planck scale of the effective field theory. The second is that WNI can inevitably provide perfect conditions for the production of primordial black holes (PBHs) in the golden window of black-hole mass range ($10^{-16} -10^{-11}M_{\odot}$) where it can account for all of the the dark matter content of the universe while satisfying observational constraints.

We analyze the role of sterile neutrinos in the framework of the $\mu\nu$SSM, where the presence of right-handed neutrinos provides a simultaneous solution to $\mu$- and $\nu$-problems in supersymmetry. We adopt a minimalistic approach, reproducing light neutrino masses and mixing angles at tree-level using just two right-handed neutrinos as part of the seesaw mechanism. A third right-handed neutrino does not contribute significantly to the mass of the three active ones, behaving as a sterile neutrino with a mass in the range keV$-$MeV. Furthermore, a keV sterile neutrino can be a good candidate for dark matter with a lifetime larger than the age of the Universe. In particular, the tree-body decay to active neutrinos gives the dominant contribution to its lifetime. The one-loop decay to gamma and active neutrino is subdominant, but relevant for observations such as astrophysical X-rays. We find regions of the parameter space of the $\mu\nu$SSM, with different values of the sterile neutrino mass, fulfilling not only these constraints but also collider constraints from the Higgs sector. Interestingly, the claimed $3.5$~keV line detection can also potentially be explained in the $\mu\nu$SSM with a sterile neutrino of $7$~keV mass.

We explore the possibility of a scalar field driving ekpyrotic contraction through a non-canonical kinetic energy density rather than a negative potential. We find that this kinetically-driven ekpyrosis ("k-ekpyrosis") can be achieved in a variety of models, including scalar field theories with power-law, polynomial, or DBI-like kinetic terms in the action. Of these examples, the ekpyrotic phase is best sustained in power-law models, which can generate large and constant equation-of-state parameters, followed by DBI-like models, which can exhibit dynamical attractors toward similarly large equations of state. We show that for a broad class of theories including these examples, phases of k-ekpyrosis are accompanied by a preceding or concurrent phase of superluminality.

Motivated by the exciting features and a recent proposed general form of the function of non-metricity scalar Q, we investigate the cosmological implications in $f(Q)$ gravity, through the resulting effective dark energy sector, extracting analytical expressions for the dark energy density, equation-of-state parameter and the deceleration parameters. We show that even in the absence of a cosmological constant, the universe exhibits the usual thermal history, with the sequence of matter and dark energy eras, and the dark-energy equation-of-state parameter always lie in the phantom regime. Additionally, calculating the age of the universe, through the extracted analytical equations of the scenario at hand, we show that the result coincide with the value corresponding to $\Lambda$CDM scenario within 1$\sigma$. Moreover, we show the excellent agreement of the scenario at hand with Supernovae type Ia observational data. Lastly, comparing the cosmological behavior in the case of the absence of an explicit cosmological constant, with the one of the presence of a cosmological constant we show that $f(Q)$ gravity can mimic the cosmological constant in a very efficient way, providing very similar behavior, revealing the advantages and capabilitites of the scenario at hand.

We introduce a time-delay function in bulk viscosity cosmology. Even for bulk viscosity functions where closed-form solutions are known, because of the time-delay term the exact solutions are lost. Therefore in order to study the cosmological evolution of the resulting models we perform a detail analysis of the stability of the critical points, which describe de Sitter solutions, by using Lindstedt's method. We find that for the stability of the critical points it depends also on the time-delay parameter, where a critical time-delay value is found which play the role of a bifurcation point. For time-delay values near to the critical value, the cosmological evolution has a periodic evolution, this oscillating behaviour is because of the time-delay function. We find a new behaviour near the exponential expansion point, which can be seen also as an alternative way to exit the exponential inflation.

In this paper, we construct \textit{Dirac-boson stars} (DBSs) model composed of a scalar field and two Dirac fields. The scalar field and both Dirac fields are in the ground state. We consider the solution families of the DBSs for the synchronized frequency $\tilde{\omega}$ and the nonsynchronized frequency $\tilde{\omega}_D$ cases, respectively. We find several different solutions when the Dirac mass $\tilde{\mu}_D$ and scalar field frequency $\tilde{\omega}_S$ are taken in some particular ranges. In contrast, no similar case has been found in previous studies of multistate boson stars. Moreover, we discuss the characteristics of each type of solution family of the DBSs and present the relationship between the ADM mass $M$ of the DBSs and the synchronized frequency $\tilde{\omega}$ or the nonsynchronized frequency $\tilde{\omega}_D$. Finally, we calculate the binding energy $E_B$ of the DBSs and investigate the relationship of $E_B$ with the synchronized frequency $\tilde{\omega}$ or the nonsynchronized frequency $\tilde{\omega}_D$.

Dark matter is estimated to make up ~60% of all normal/baryonic matter but cannot be directly imaged. Despite the fact that dark matter cannot be directly observed yet, its influence on the motion of stars and gas in spiral galaxies has been detected. One way to show motion in galaxies is to produce rotation curves that are plots of velocity measurements of how fast stars and gas move in a galaxy around the center of mass. According to Newton's Law of Gravitation, the rotational velocity is an indication of the amount of visible and non-visible mass in the galaxy. Given that the visible matter can be estimated using photometry, dark matter mass in galaxies can be calculated. In order to gain a greater appreciation of the research scientists' findings about dark matter, their method should be easily reproduced by any curious individual. Our interactive workshop is an excellent educational tool to investigate how dark matter impacts the rotation of visible matter by providing a guide to produce galactic rotation curves. The Python-based notebooks are set up to walk you through the steps of producing rotation curves and to allow you to learn about each component of the galaxy. The three steps of the rotation curve building process is plotting the measured velocity data, constructing the rotation curves for each component, and fitting the total velocity to the measured values. After completing all the modules of the workshop, reading and understanding scientific journals of dark matter findings should be considerably easier.