Papers for Wednesday, Mar 27 2024

The Voids-within-Voids-within-Voids (VVV) project used dark-matter-only simulations to study the abundance and structure of dark matter haloes over the full mass range populated in the standard $\Lambda\mathrm{CDM}$ cosmology. Here we explore how baryonic effects modify these results for $z=0$ halo masses in the range $10^4$ to $10^7~\mathrm{M_\odot}$, below the threshold for galaxy formation. Our main study focuses on three simulations from identical initial conditions at $z=127$, one following dark matter only, one including non-radiative gas, and one additionally including the baryonic physics relevant in this halo mass range (cooling and photoheating). In the non-radiative simulation, above $10^{5.5}~\mathrm{M_\odot}$, halo abundance and internal structure are very similar to the dark-matter-only simulation, and the baryon to dark matter ratio is everywhere close to the cosmic value. At lower mass, this ratio drops and haloes are less concentrated and less massive in the non-radiative case. Test simulations at higher resolution show this to be mainly a resolution effect; the expected drop in baryon content due to residual pressure effects only becomes substantial for $z=0$ haloes below $\sim 10^{2.7}~\mathrm{M_\odot}$. However, gas is heated by reionization at $z=6$ in our ``full physics'' run, and this results in almost complete expulsion of gas from all haloes in our simulated mass range. This suppresses the halo mass function by $\sim 30 \%$, lowers halo concentration, and consequently weakens the dark matter annihilation signal by $\sim 40-60 \%$.

Using more than 100 galaxies in the MUSE Ultra Deep Field with spectroscopy from the Hubble Space Telescope's Wide Field Camera 3 and the Very Large Telescope's Multi Unit Spectroscopic Explorer, we extend the gas-phase mass-metallicity relation (MZR) at $z\approx\,$1$\,$-$\,$2 down to stellar masses of M$_{\star}$ $\approx$ 10$^{7.5}$ M$_{\odot}$. The sample reaches six times lower in stellar mass and star formation rate (SFR) than previous HST studies at these redshifts, and we find that galaxy metallicities decrease to log(O/H) + 12 $\approx$ 7.8 $\pm$ 0.1 (15% solar) at log(M$_{\star}$/M$_{\odot}$) $\approx$ 7.5, without evidence of a turnover in the shape of the MZR at low masses. We validate our strong-line metallicities using the direct method for sources with [O III] $\lambda$4363 and [O III] $\lambda$1666 detections, and find excellent agreement between the techniques. The [O III] $\lambda$1666-based metallicities double existing measurements with S/N $\geq$ 5 for unlensed sources at $z~>$ 1, validating the strong-line calibrations up to $z \sim$2.5. We confirm that the MZR resides $\sim$0.3 dex lower in metallicity than local galaxies and is consistent with the fundamental metallicity relation (FMR) if the low mass slope varies with SFR. At lower redshifts ($z\sim$0.5) our sample reaches $\sim$0.5 dex lower in SFR than current calibrations and we find enhanced metallicities that are consistent with extrapolating the MZR to lower SFRs. Finally, we detect only a $\sim$0.1 dex difference in the metallicities of galaxies in groups versus isolated environments. These results are based on robust calibrations and reach the lowest masses and SFRs that are accessible with HST, providing a critical foundation for studies with the Webb and Roman Space Telescopes.

While the tip of the red giant branch (TRGB) has been used as a distance indicator since the early 1990's, its application to measure the Hubble Constant as a primary distance indicator occurred only recently. The TRGB is also currently at an interesting crossroads as results from the James Webb Space Telescope (JWST) are beginning to emerge. In this chapter, we provide a review of the TRGB as it is used to measure the Hubble constant. First, we provide an essential review of the physical and observational basis of the TRGB as well as providing a summary for its use for measuring the Hubble Constant. More attention is then given is then given to recent, but still pre-JWST, developments, including new calibrations and developments with algorithms. We also address challenges that arise while measuring a TRGB-based Hubble Constant. We close by looking forward to the exciting prospects from telescopes such as JWST and Gaia.

In this study we use MUSE Integral Field Spectroscopy (IFS), along with multi-line diagnostics, for the optical identification of Supernova Remnants (SNRs) in the galaxy NGC 7793. We find in total 238 SNR candidates, 225 of them new identifications, increasing significantly the number of known SNRs in this galaxy. The velocity dispersion of the candidate SNRs was calculated, giving a mean value of 27 $ km\,s^{ -1} $. We construct the H$\rm \alpha$, [S II], [O III], and [S II] - H$\rm \alpha$ luminosity functions, and for the first time, the [N II], [N II] - H$\rm \alpha$, [N II] - [S II], [O III] - [S II], and [O III] - [N II] luminosity functions of the candidate SNRs. Shock models, along with the observed multi-line information were used, in order to estimate shock velocities. The $\sim$ 65% of the SNRs present velocities < 300 $ km\,s^{ -1} $. There is a clear correlation between shock velocity and [O III]/H$\rm \beta$ ratio, and a less clear but still evident correlation in the relation between shock velocity and the [S II]/H$\rm \alpha$, [N II]/H$\rm \alpha$ ratios. We also use the [S II]6716/31 ratio of the SNR candidates to calculate their post-shock density, assuming different temperatures. The median value of the density of our sample is $\sim$ 80 $ cm^{ -3} $, for a temperature of T = $10^4\,K$. No correlation between shock velocity and density, or density and SNRs with [S II]/H$\rm \alpha$ > 0.4 and [S II]/H$\rm \alpha$ < 0.4 is observed.

Forthcoming measurements of the Sunyaev-Zeldovich (SZ) effect in galaxy clusters will dramatically improve our understanding of the main intra-cluster medium (ICM) properties and how they depend on the particular thermal and dynamical state of the associated clusters. Using a sample of simulated galaxy clusters we assess the impact of the ICM internal dynamics on both the thermal and kinetic SZ effects (tSZ and kSZ, respectively). We produce synthetic maps of the SZ effect, for the simulated clusters. For each galaxy cluster in the sample, its dynamical state is estimated by using a combination of well-established indicators. We use the correlations between SZ maps and cluster dynamical state, to look for the imprints of the evolutionary events, mainly mergers, on the SZ signals. The kinetic effect shows a remarkable correlation with the dynamical state: unrelaxed clusters present a higher radial profile and an overall stronger signal at all masses and radii. Furthermore, the kSZ signal is correlated with rotation for relaxed clusters, while for the disturbed systems the effect is dominated by other motions such as bulk flows, turbulence, etc. The kSZ effect shows a dipolar pattern when averaging over cluster dynamical classes, especially for the relaxed population. This feature can be exploited to stack multiple kSZ maps in order to recover a stronger dipole signal that would be correlated with the global rotation properties of the sample. The SZ effect can be used as a tool to estimate the dynamical state of galaxy clusters, especially to segregate those clusters with a quiescent evolution from those with a rich record of recent merger events. Our results suggest that the forthcoming observational data measuring the SZ signal in clusters could be used as a complementary strategy to classify the evolutionary history of galaxy clusters.

We revisit the picture of X-ray emission of groups through the study of systematic differences in the optical properties of groups with and without X-ray emission and study the effect of large-scale density field on scaling relations. We present the identification of X-ray galaxy groups using a combination of RASS and SDSS data. We include new X-ray reanalysis of RASS, to include very extended (up to a size of half a degree) sources and account for differences in the limiting sensitivity towards compact and very extended X-ray emission. X-ray groups exhibit less scatter in the scaling relations and selecting the groups based on the extended X-ray emission leads to an additional scatter reduction. Most of the scatter for the optical groups is associated with a small (6%) fraction of outliers, primarily associated with low optical luminosity groups found in dense regions of the cosmic web. These groups are the primary candidates for being the contaminants in the optical group catalogues. Removing those groups from the optical group sample using optically measured properties only, leads to a substantial reduction in the scatter in the most scaling relations of the optical groups. We find a density dependence of both the X-ray and optical luminosity of groups, which we associate with the assembly bias. Abridged.

Multi-wavelength extragalactic nuclear transients, particularly those detectable as multi-messengers, are among the primary drivers for the next-generation observatories. X-ray quasi-periodic eruptions (QPEs) are the most recent and perhaps most peculiar addition to this group. Here, we report a first estimate of the volumetric rate of QPEs based on the first four discoveries with the eROSITA X-ray telescope onboard the Spectrum Roentgen Gamma observatory. Under the assumption, supported by a suite of simulated light curves, that these four sources sample the intrinsic population somewhat homogeneously, we correct for their detection efficiency and compute a QPE abundance of $\mathscr{R}_{{\rm vol}} = 0.60_{-0.43}^{+4.73} \times 10^{-6}\,$Mpc$^{-3}$ above an intrinsic average $\log L_{\rm 0.5-2.0\,keV}^{\rm peak} > 41.7$. Since the exact lifetime of QPEs ($\tau_{\rm life}$) is currently not better defined than between a few years or few decades, we convert this to a formation rate of $\mathscr{R}_{\rm vol}/\tau_{\rm life}\approx 0.6 \times 10^{-7} (\tau_{\rm life}/10\,\mathrm{y})^{-1}\,$Mpc$^{-3}\,$year$^{-1}$. As a comparison, this value is a factor $\sim10\,\tau_{\rm life}$ times smaller than the formation rate of tidal disruption events. The origin of QPEs is still debated, although lately most models suggest that they are the electromagnetic counterpart of extreme mass ratio inspirals (EMRIs). In this scenario, the QPE rate would thus be the first-ever constraint (i.e. a lower limit) to the EMRI rate from observations alone. Future discoveries of QPEs and advances in their theoretical modeling will consolidate or rule out their use for constraining the number of EMRIs detectable by the LISA mission.

We present a new all-sky catalogue of X-ray detected groups (AXES-2MRS), based on the identification of large X-ray sources found in the ROSAT All-Sky Survey (RASS) with the Two Micron Redshift Survey (2MRS) Bayesian Group Catalogue. We study the basic properties of these galaxy groups to gain insights into the effect of different group selections on the properties. In addition to X-ray luminosity coming from shallow survey data of RASS, we have obtained detailed X-ray properties of the groups by matching the AXES-2MRS catalogue to archival X-ray observations by XMM-Newton and complemented this by adding the published XMM-Newton results on galaxy clusters in our catalogue. We analyse temperature and density to the lowest overdensity accessible by the data, obtaining hydrostatic mass estimates and comparing them to the velocity dispersions of the groups. We find a large spread in the central mass to virial mass ratios for galaxy groups in the XMM-Newton subsample. This can either indicate large non-thermal pressure of galaxy groups affecting our X-ray mass measurements, or the effect of a diversity of halo concentrations on X-ray properties of galaxy groups. Previous catalogues, based on detecting the peak of the X-ray emission preferentially sample the high-concentration groups, while our new catalogue includes many low-concentration groups. Abridged.

The discovery of planetary systems beyond the solar system has revealed a diversity of architectures, most of which differ significantly from our system. The initial detection of an exoplanet is often followed by subsequent discoveries within the same system as observations continue, measurement precision is improved, or additional techniques are employed. The HD 104067 system is known to consist of a bright K dwarf host star and a giant planet in a $\sim$55 day period eccentric orbit. Here we report the discovery of an additional planet within the HD 104067 system, detected through the combined analysis of radial velocity data from the HIRES and HARPS instruments. The new planet has a mass similar to Uranus and is in an eccentric $\sim$14 day orbit. Our injection-recovery analysis of the radial velocity data exclude Saturn-mass and Jupiter-mass planets out to 3 AU and 8 AU, respectively. We further present TESS observations that reveal a terrestrial planet candidate ($R_p = 1.30\pm0.12$ $R_\oplus$) in a $\sim$2.2~day period orbit. Our dynamical analysis of the three planet model shows that the two outer planets produce significant eccentricity excitation of the inner planet, resulting in tidally induced surface temperatures as high as $\sim$2600 K for an emissivity of unity. The terrestrial planet candidate may therefore be caught in a tidal storm, potentially resulting in its surface radiating at optical wavelengths.

Context. White-light stellar flares are proxies for some of the most energetic types of flares, but their triggering mechanism is still poorly understood. As they are associated with strong X and UV emission, their study is particularly relevant to estimate the amount of high-energy irradiation onto the atmospheres of exoplanets, especially those in their stars' habitable zone. Aims. We used the high-cadence, high-photometric capabilities of the CHEOPS and TESS space telescopes to study the detailed morphology of white-light flares occurring in a sample of 130 late-K and M stars, and compared our findings with results obtained at a lower cadence. We developed dedicated software for this purpose. Results. Multi-peak flares represent a significant percentage ($\gtrsim 30$\%) of the detected outburst events. Our findings suggest that high-impulse flares are more frequent than suspected from lower-cadence data, so that the most impactful flux levels that hit close-in exoplanets might be more time-limited than expected. We found significant differences in the duration distributions of single-peak and complex flare components, but not in their peak luminosity. A statistical analysis of the flare parameter distributions provides marginal support for their description with a log-normal instead of a power-law function, leaving the door open to several flare formation scenarios. We tentatively confirmed previous results about quasi-periodic pulsations in high-cadence photometry, report the possible detection of a pre-flare dip, and did not find hints of photometric variability due to an undetected flare background. Conclusions. The high-cadence study of stellar hosts might be crucial to evaluate the impact of their flares on close-in exoplanets, as their impulsive phase emission might otherwise be incorrectly estimated. Future telescopes such as PLATO and Ariel will help in this respect.

Sigmoids produce strong eruptive events. Earlier studies have shown that the ICME axial magnetic field Bz can be predicted with some credibility by observing the corresponding filament or the polarity inversion line in the region of eruption and deriving the magnetic field direction from that. Sigmoids are coronal structures often associated with filaments in the sigmoidal region. In this study, firstly we compare filament chirality with sigmoid handedness to observe their correlation. Secondly, we perform non-linear force-free approximations of the coronal magnetic connectivity using photospheric vector magnetograms underneath sigmoids to obtain a weighted-average value of the force-free parameter and to correlate it with filament chirality and the observed coronal sigmoid handedness. Importantly, we find that the sigmoids and their filament counterparts do not always have the same helicity signs. Production of eruptive events by regions that do not have the same signs of helicities is $\sim$3.5 times higher than when they do. A case study of magnetic energy/ helicity evolution in NOAA AR 12473 is also presented.

The presence of dual active galactic nuclei (AGN) on scales of a few tens of kpc can be used to study merger-induced accretion on supermassive black holes (SMBHs) and offer insights about SMBH mergers, using dual AGNs as merger precursors. In this study, we use the Romulus25 cosmological simulation to investigate the properties and evolution of dual AGNs. We first analyze the properties of AGNs ($L_{bol} > 10^{43} \rm erg/s$) and their neighboring SMBHs (any SMBHs closer than 30 pkpc to an AGN) at $z \leq 2$. This is our underlying population. Subsequently, we applied the luminosity threshold of $L_{bol} > 10^{43} erg/s$ to the neighboring SMBHs thereby identifying dual and multiple AGNs in our simulation. We examined the properties and statistics of dual AGNs in comparison to single AGNs. Our findings indicate an increase in the number of both single and dual AGNs from lower to higher redshifts. All dual AGNs in our sample resulted from major mergers. Compared with the single AGN population, duals are characterized by a lower black hole to halo mass ratio. We found that the properties of dual AGN host halos, including halo mass, stellar mass, star formation rate (SFR), and gas mass, are generally consistent with those of single AGN halos, albeit tending towards the higher end of their respective property ranges. Our analysis uncovered a diverse array of evolutionary patterns among dual AGNs, including rapidly evolving systems, slower ones, and instances where SMBH mergers are ineffective.

This study presents a semi-analytical model that models the effects that a static magnetized medium with a tangled field produces in relativistic collimated jets. The model is a first approximation that addresses the magnetic field present in the medium and is based on pressure equilibrium principles between the jet, cocoon, and external medium. A fraction of the ambient medium field is allowed to be entrained in the cocoon. We find that the jet and cocoon properties are affected at high magnetic fields and that they are sensitive to the mixing in the cocoon. For low-mixing a slower-broader jet with a broader and more energetic cocoon would be produced. On the other hand, high-mixing would produce a faster-narrower jet with a narrow and less-energetic cocoon. Two-dimensional hydrodynamical simulations are used to validate the model and to constrain the mixing parameter. We found good qualitative agreement between the model and the simulations. For high magnetization, the results were found to be consistent with the low mixing case of our semi-analytical model.

Dust grains are important in various astrophysical processes and serve as indicators of interstellar medium structures, density, and mass. Understanding their physical properties and chemical composition is crucial in astrophysics. Dust polarization is a valuable tool for studying these properties. The Radiative Torque (RAT) paradigm, which includes Radiative Torque Alignment (RAT-A) and Radiative Torque Disruption (RAT-D), is essential to interpret the dust polarization data and constrain the fundamental properties of dust grains. However, it has been used primarily to interpret observations at a single wavelength. In this study, we analyze the thermal dust polarization spectrum obtained from observations with SOFIA/HAWC+ and JCMT/POL-2 in the OMC-1 region and compare the observational data with our numerical results using the RAT paradigm. We find that the dense gas exhibits a positive spectral slope, whereas the warm regions show a negative one. We demonstrate that a one-layer dust (one-phase) model can only reproduce the observed spectra at certain locations and cannot match those with prominent V-shape spectra (for which the degree of polarization initially decreases with wavelength from 54 to $\sim$ 300$\,\mu$m and then increases at longer wavelengths). To address this, we improve our model by incorporating two dust components (warm and cold) along the line of sight, resulting in a two-phase model. This improved model successfully reproduces the V-shaped spectra. The best model corresponds to a mixture composition of silicate and carbonaceous grains in the cold medium. Finally, by assuming the plausible model of grain alignment, we infer the inclination angle of the magnetic fields in OMC-1. This approach represents an important step toward better understanding the physics of grain alignment and constraining 3D magnetic fields using dust polarization spectrum.

Protoplanetary disks are often assumed to change slowly and smoothly during planet formation. Here, we investigate the time evolution of isolated disks subject to viscosity and a disk wind. The viscosity is assumed to increase rapidly at around 900 K due to thermal ionization of alkali metals, or thermionic and ion emission from dust, and the onset of magneto-rotational instability (MRI). The disks generally undergo large, rapid fluctuations for a wide range of time-averaged mass accretion rates. Fluctuations involve coupled waves in temperature and surface density that move radially in either direction through the inner 1.5 AU of the disk. Two types of wave are seen with radial speeds of roughly 50 and 1000 cm/s respectively. The pattern of waves repeats with a period of roughly 10,000 years that depends weakly on the average mass accretion rate. Viscous transport due to MRI is confined to the inner disk. This region is resupplied by mass flux from the outer disk driven by the disk wind. Interior to 1 AU, the temperature and surface density can vary by a factor of 2--10 on timescales of years to ky. The stellar mass accretion rate varies by 3 orders of magnitude on a similar timescale. This behavior lasts for at least 1 My for initial disks comparable to the minimum-mass solar nebula.

We present JWST NIRCam imaging of B14-65666 ("Big Three Dragons"), a bright Lyman-break galaxy system ($M_\text{UV}=-22.5$ mag) at $z=7.15$. The high angular resolution of NIRCam reveals the complex morphology of two galaxy components: galaxy E has a compact core (E-core), surrounded by diffuse, extended, rest-frame optical emission, which is likely to be tidal tails; and galaxy W has a clumpy and elongated morphology with a blue UV slope ($\beta_\text{UV}=-2.2\pm0.1$). The flux excess, F356W$-$F444W, peaks at the E-core ($1.05^{+0.08}_{-0.09}$ mag), tracing the presence of strong [OIII] 4960,5008 \r{A} emission. ALMA archival data show that the bluer galaxy W is brighter in dust continua than the redder galaxy E, while the tails are bright in [OIII] 88 $\mathrm{\mu m}$. The UV/optical and sub-mm SED fitting confirms that B14-65666 is a major merger in a starburst phase as derived from the stellar mass ratio (3:1 to 2:1) and the star-formation rate, $\simeq1$ dex higher than the star-formation main sequence at the same redshift. The galaxy E is a dusty ($A_\text{V}=1.2\pm0.1$ mag) starburst with a possible high dust temperature ($\ge63$-$68$ K). The galaxy W would have a low dust temperature ($\le27$-$33$ K) or patchy stellar-and-dust geometry, as suggested from the infrared excess (IRX) and $\beta_\text{UV}$ diagram. The high optical-to-FIR [OIII] line ratio of the E-core shows its lower gas-phase metallicity ($\simeq0.2$ Z$_{\odot}$) than the galaxy W. These results agree with a scenario where major mergers disturb morphology and induce nuclear dusty starbursts triggered by less-enriched inflows. B14-65666 shows a picture of complex stellar buildup processes during major mergers in the epoch of reionization.

The co-evolution between supermassive black holes and their environment is most directly traced by the hot atmospheres of dark matter halos. Cooling of the hot atmosphere supplies the central regions with fresh gas, igniting active galactic nuclei (AGN) with long duty cycles. Outflows from the central engine tightly couple with the surrounding gaseous medium and provide the dominant heating source preventing runaway cooling. Every major modern hydrodynamical simulation suite now includes a prescription for AGN feedback to reproduce realistic populations of galaxies. However, the mechanisms governing the feeding/feedback cycle between the central black holes and their surrounding galaxies and halos are still poorly understood. Galaxy groups are uniquely suited to constrain the mechanisms governing the cooling-heating balance, as the energy supplied by the central AGN can exceed the gravitational binding energy of halo gas particles. Here we provide a brief overview of our knowledge of the impact of AGN on the hot atmospheres of galaxy groups, with a specific focus on the thermodynamic profiles of groups. We then present our on-going efforts to improve on the implementation of AGN feedback in galaxy evolution models by providing precise benchmarks on the properties of galaxy groups. We introduce the \XMM~ Group AGN Project (X-GAP), a large program on \XMM~ targeting a sample of 49 galaxy groups out to $R_{500c}$.

The connection between the binary black hole (BBH) mergers observed by LIGO-Virgo-KAGRA (LVK) and their stellar progenitors remains uncertain. Specifically, the fraction $\epsilon$ of stellar mass that ends up in BBH mergers and the delay time $\tau$ between star formation and BBH merger carry information about the astrophysical processes that give rise to merging BBHs. We model the BBH merger rate in terms of the cosmic star formation history, coupled with a metallicity-dependent efficiency $\epsilon$ and a distribution of delay times $\tau$, and infer these parameters with data from the Third Gravitational-Wave Transient Catalog (GWTC-3). We find that the progenitors to merging BBHs preferentially form in low metallicity environments with a low metallicity efficiency of $\log_{10}\epsilon_{<Z_t}=-3.99^{+0.68}_{-0.87}$ and a high metallicity efficiency of $\log_{10}\epsilon_{<Z_t}=-4.60^{+0.30}_{-0.34}$ at the 90% credible level. The data also prefer short delay times. For a power-law distribution $p(\tau)\propto \tau^\alpha$, we find $\tau_\text{min}<1.9 $ Gyr and $\alpha<-1.32$ at 90% credibility. Our model allows us to extrapolate the mass density in BBHs out to high redshifts. We cumulatively integrate our modelled density rate over cosmic time to get the total mass density of merging stellar mass BBHs as a function of redshift. Today, stellar-mass BBH mergers make up only $\sim 0.01\%$ of the total stellar mass density created by high-mass ($>10\,M_\odot$) progenitors. However, because massive stars are so short-lived, there may be more mass in merging BBHs than in living massive stars as early as $\sim 2.5$ Gyr ago. We also compare to the mass in supermassive BHs, finding that the mass densities were comparable $\sim 12.5$ Gyr ago, but the mass density in SMBHs quickly increased to $\sim 75$ times the mass density in merging stellar mass BBHs by $z\sim 1$.

We use the kinematics of field OB stars to estimate the frequencies of runaway stars generated by the dynamical ejection scenario (DES), the binary supernova scenario (BSS), and the combined two-step mechanism. We update the proper motions for field OB and OBe stars in the Small Magellanic Cloud (SMC) using Gaia DR3. Our sample now contains 336 stars from the Runaways and Isolated O-Type Star Spectroscopic Survey of the SMC (RIOTS4), and we update our algorithm to calculate more accurate velocities compared to those obtained previously from DR2. We find a decrease in median velocity from 39 to 29 km/s, implying that the proper motions from our previous work were systematically overestimated. We present the velocity distribution for OBe stars and quantitatively compare it to those of non-compact binaries and high-mass X-ray binaries. We confirm that OBe stars appear to be dominated by the BSS and are likely post-SN binary systems, further supporting the mass-transfer model to explain the origin of their emission-line disks. In contrast, normal OB stars may show a bimodal velocity distribution, as may be expected from different processes that occur with dynamical ejections. The kinematics of fast-rotating OB stars are similar to those of normal OB stars rather than OBe stars, suggesting that the origin of their high v_r*sin(i) is different from that of OBe stars. We update our model parameters describing the kinematic origins of the SMC field population, still confirming that for runaway stars, the DES mechanism dominates, and two-step ejections seem comparable in frequency to pure BSS ejections.

The multi-wavelength outburst activity in the BL Lacertae source OJ 287 has sparked a lot of controversy about whether the source contains one or two black holes (BHs) and what characteristics of this black hole binary would be. In this article we present the results of analysis of the X-ray flaring activity of OJ 287 using the data of Swift/XRT observations. We discovered that the energy spectra in all spectral states can be adequately fit with the XSPEC BMC model (the Comptonization one). As a result we found that the X-ray photon index of the BMC model, $\Gamma$ correlates with the mass accretion rate, $\dot M$. We found the photon index $\Gamma$ to increase monotonically with accretion rate $\dot M$ from $\Gamma\sim 2$ in the intermediate state (IS) to $\Gamma\sim2.5$ the high/soft state (HSS) with subsequent saturation at $\Gamma\sim$ 2.6 level at higher luminosities. This type of behavior of the spectral index is remarkably similar to the pattern observed in a number of established stellar-mass black hole candidates. Assuming the universality of the observed pattern of the correlation between the photon index and mass accretion rate, we estimate the BH mass in OJ 287 to be around $2\times10^8$ solar masses, using the well studied BH binaries GX 339-4 and XTE J1859-226 to calibrate the model.

The quantum kinetic evolution of neutrinos in dense environments, such as the core-collapse supernovae or the neutron star mergers, can result in fast flavor conversion (FFC), presenting a significant challenge to achieving robust astrophysical modeling of these systems. Recent works that directly simulate the quantum kinetic transport of neutrinos in localized domains have suggested that the asymptotic outcome of FFCs can be modeled by simple analytical prescriptions. In this Letter, we incorporate the analytical prescriptions into global simulations that solve the classical neutrino transport equation including collisions and advection under spherical symmetry. We demonstrate excellent agreement between results obtained using this approach and those directly from the corresponding global quantum kinetic simulations. In particular, this effective method can also precisely capture the collisional feedback effect for cases where the FFC happens inside the neutrinosphere. Our work highlights that a robust integration of FFCs in classical neutrino transport used in astrophysical simulation can be feasible.

The unprecedented imaging power of JWST provides new abilities to observe the shapes of objects in the early Universe in a way that has not been possible before. Recently, JWST acquired a deep field image inside the same field imaged in the past as the HST Ultra Deep Field. Computer-based quantitative analysis of spiral galaxies in that field shows that among 34 galaxies for which their rotation of direction can be determined by the shapes of the arms, 24 rotate clockwise, and just 10 rotate counterclockwise. The one-tailed binomial distribution probability to have asymmetry equal or stronger than the observed asymmetry by chance is $\sim$0.012. While the analysis is limited by the small size of the data, the observed asymmetry is aligned with all relevant previous large-scale analyses from all premier digital sky surveys, all show a higher number of galaxies rotating clockwise in that part of the sky, and the magnitude of the asymmetry increases as the redshift gets higher. This paper also provides data and analysis to reproduce previous experiments suggesting that the distribution of galaxy rotation in the Universe is random, to show that the exact same data used in these studies in fact show non-random distribution, and in excellent agreement with the results shown here. These findings reinforce consideration of the possibility that the directions of rotation of spiral galaxies as observed from Earth are not necessarily randomly distributed. The explanation can be related to the large-scale structure of the Universe, but can also be related to a possible anomaly in the physics of galaxy rotation.

In this paper, we investigate the benefits of teaming up data from the radio to the far- 1 infrared (FIR) regime for the characterization of Dusty Star-Forming Galaxies (DSFGs). These galaxies 2 are thought to be the star-forming progenitors of local massive quiescent galaxies, and play a pivotal 3 role in the reconstruction of the cosmic star formation rate density up to high redshift. Due to their 4 dust-enshrouded nature, DSFGs are often invisible in the near-infrared/Optical/UV bands. Therefore, 5 they necessitate observations at longer wavelengths, primarily the FIR, where dust emission occurs, 6 and radio, which is not affected by dust absorption. Combining data from these two spectral windows 7 makes it possible to characterize even the dustiest objects, enabling the retrieval of information about 8 their age, dust temperature, and star-formation status, and facilitates the differentiation between 9 various galaxy populations that evolve throughout cosmic history. Despite the detection of faint radio 10 sources being a challenging task, this study demonstrates that an effective strategy to build statistically 11 relevant samples of DSFGs would be reaching deep sensitivities in the radio band, even restricted to 12 smaller areas, and then combining these radio observations with FIR/submm data. Additionally, 13 the paper quantifies the improvement in the Spectral Energy Distribution (SED) reconstruction of 14 DSFGs by incorporating ALMA band measurements, in particular, in its upgraded status thanks to 15 the anticipated Wideband Sensitivity Upgrade.

Motivated by measurements of the rotation speed of accretor stars in post-mass-transfer (post-MT) systems, we investigate how magnetic braking affects the spin-down of individual stars during binary evolution with the \texttt{MESAbinary} module. Unlike the conventional assumption of tidal synchronization coupled with magnetic braking in binaries, we first calculate whether tides are strong enough to synchronize the orbit. Subsequently, this influences the spin-down of stars and the orbital separation. In this study, we apply four magnetic braking prescriptions to reduce the spin angular momentum of the two stars throughout the entire binary evolution simulation. Our findings reveal that despite magnetic braking causing continuous spin-down of the accretor, when the donor begins to transfer material onto the accretor, the accretor can rapidly spin up to its critical rotation rate. After MT, magnetic braking becomes more important in affecting the angular momentum evolution of the stars. Post-MT accretor stars thus serve as a valuable testbed for observing how the magnetic braking prescriptions operate in spinning down stars from their critical rotation, including the saturation regimes of the magnetic braking. The rotation rate of the accretor star, combined with its mass, could provide age information since the cessation of MT. By comparing the models against observation, the magnetic braking prescription by Garraffo et al. (2018b) is found to better align with the rotation data of post-MT accretors.

We present a novel maximum a posteriori estimator to jointly estimate band-powers and the covariance of the three-dimensional power spectrum (P3D) of Lyman-alpha forest flux fluctuations, called MAPLE. Our Wiener-filter based algorithm reconstructs a window-deconvolved P3D in the presence of complex survey geometries typical for Lyman-alpha surveys that are sparsely sampled transverse to and densely sampled along the line-of-sight. We demonstrate our method on idealized Gaussian random fields with two selection functions: (i) a sparse sampling of 30 background sources per square degree designed to emulate the currently observing the Dark Energy Spectroscopic Instrument (DESI); (ii) a dense sampling of 900 background sources per square degree emulating the upcoming Prime Focus Spectrograph Galaxy Evolution Survey. Our proof-of-principle shows promise, especially since the algorithm can be extended to marginalize jointly over nuisance parameters and contaminants, i.e.offsets introduced by continuum fitting. Our code is implemented in JAX and is publicly available on GitHub.

Vortex fiber nulling (VFN) is a technique for detecting and characterizing faint companions at small separations from their host star. A near-infrared ($\sim2.3 \mu$m) VFN demonstrator mode was deployed on the Keck Planet Imager and Characterizer (KPIC) instrument at the Keck Observatory and presented earlier. In this paper, we present the first VFN companion detections. Three targets, HIP 21543 Ab, HIP 94666 Ab, and HIP 50319 B, were detected with host-companion flux ratios between 70 and 430 at and within one diffraction beamwidth ($\lambda/D$). We complement the spectra from KPIC VFN with flux ratio and position measurements from the CHARA Array to validate the VFN results and provide a more complete characterization of the targets. This paper reports the first direct detection of these three M dwarf companions, yielding their first spectra and flux ratios. Our observations provide measurements of bulk properties such as effective temperatures, radial velocities, and v$\sin{i}$, and verify the accuracy of the published orbits. These detections corroborate earlier predictions of the KPIC VFN performance, demonstrating that the instrument mode is ready for science observations.

An emission line at ~6.7 keV is attributable to a He-like iron K-shell transition, which indicates existence of a thin thermal plasma with a temperature of several keV. Using Suzaku archival data, we searched for the iron K-line from the spiral galaxy NGC 6946, and found the iron K-line at 6.68+/-0.07 keV at the 3.1 sigma level in the central r<2.'5 region. The iron line luminosity from the central region was estimated to be (2.3+/-1.2)x10^{37} erg s^{-1} at a distance of 5.5 Mpc. The origin of the iron emission line is discussed.

The solar wind upstream of Mars's bow shock can be described in terms of Alfv\'enic turbulence, with an incompressible energy cascade rate of $10^{-17}$ J m$^{-3}$ s$^{-1}$ at magnetohydrodynamics (MHD) scales. The solar wind has more Alfv\'en waves propagating outwards from the Sun (than inwards) and a median Alfv\'en ratio of $\sim0.33$. Newly ionized planetary protons associated with the extended hydrogen corona generate waves at the local proton cyclotron frequency. These 'proton cyclotron waves' (PCW) mostly correspond to fast magnetosonic waves, although the ion cyclotron (Alfv\'enic) wave mode is possible for large Interplanetary Magnetic Field cone angles. PCW do not show significant effects on the solar wind energy cascade rates at MHD scales but could affect smaller scales. The magnetosheath displays high amplitude wave activity, with high occurrence rate of Alfv\'en waves. Turbulence appears not fully developed in the magnetosheath, suggesting fluctuations do not have enough time to interact in this small-size region. Some studies suggest PCW affect turbulence in the magnetosheath. Overall, wave activity is reduced inside the magnetic pile-up region and the Martian ionosphere. However, under certain conditions, upstream waves can reach the upper ionosphere. So far, there have not been conclusive observations of Alfv\'en waves in the ionosphere or along crustal magnetic fields, which could be due to the lack of adequate observations.

Utilizing a kinematic decomposition of simulated galaxies, we focus on galaxies with tiny kinematically inferred stellar halos, indicative of weak external influences. We investigate the intricate interplay between internal (natural) and external (nurture) processes in shaping the scaling relationships of specific angular momentum ($j_\star$), stellar mass ($M_\star$), and size of disk galaxies within the IllustrisTNG simulation. The correlation among mass, size, and angular momentum of galaxies is examined by comparing simulations with observations and the theoretical predictions of the exponential hypothesis. Galaxies with tiny stellar halos exhibit a large scatter in the $j_\star$-$M_\star$ relation, which suggests that it is inherently present in their initial conditions. The analysis reveals that the disks of these galaxies adhere to the exponential hypothesis, resulting in a tight fiducial $j_\star$-$M_\star$-scale length (size) relation that is qualitatively consistent with observations. The inherent scatter in $j_\star$ provides a robust explanation for the mass-size relation and its substantial variability. Notably, galaxies that are moderately influenced by external processes closely adhere to a scaling relation akin to that of galaxies with tiny stellar halos. This result underscores the dominant role of internal processes in shaping the overall $j_\star$-$M_\star$ and mass-size relation, with external effects playing a relatively minor role in disk galaxies. Furthermore, the correlation between galaxy size and the virial radius of the dark matter halo exists but fails to provide strong evidence of the connection between galaxies and their parent dark matter halos.

This paper presents a hydrodynamic simulation that couples detailed non-local thermodynamic equilibrium (NLTE) calculations of the hydrogen and helium level populations to model the H$\alpha$ and He 10830 transmission spectra of the hot Jupiter HAT-P-32b. A Monte Carlo simulation is applied to calculate the number of Ly$\alpha$ resonance scatterings, which is the main process for populating H(2). In the examined parameter space, only the models with H/He $\geq$ 99.5/0.5, $(0.5 \sim 3.0)$ times the fiducial value of $F_{\rm XUV}$, $\beta_m = 0.16\sim 0.3$, can explain the H$\alpha$ and He 10830 lines simultaneously. We find a mass-loss rate of $\sim (1.0\sim 3.1) \times 10^{13}$ g s$^{-1}$, consistent with previous studies. Moreover, we find that the stellar Ly$\alpha$ flux should be as high as $4 \times 10^{5}$ erg cm$^{-2}$ s$^{-1}$, indicating high stellar activity during the observation epoch of the two absorption lines. Despite the fact that the metallicity in the lower atmosphere of HAT-P-32b may be super-solar, our simulations tentatively suggest it is close to solar in the upper atmosphere. The difference in metallicity between the lower and upper atmospheres is essential for future atmospheric characterisations.

We investigate the attainable maximum energy of particles accelerated in the core-collapse supernova remnant (SNR) shock propagating in the free wind region with the Parker-spiral magnetic field, current sheet, and the wind termination shock (WTS) by using test particle simulations. This work focuses on Wolf-Rayet stars as progenitors. The magnetic field amplification in the free wind region (shock upstream region) is not considered in this work. Test particle simulations show that particles escaped from the core-collapse SNR reach and move along the WTS, and eventually return to the SNR shock from the poles or equator of the WTS. The particle attainable energy can be boosted by this cyclic motion between the SNR shock and WTS and can be larger than the particle energy that is limited by escape from the SNR shock. The particle energy limited by the cyclic motion between the SNR shock and WTS is about $10-100~{\rm TeV}$. Thus, the core-collapse SNR without upstream magnetic field amplification can be the origin of the break around $10~{\rm TeV}$ of the energy spectrum of observed cosmic ray protons and helium.

We use both chirped-pulse Fourier transform and frequency modulated absorption spectroscopy to study the rotational spectrum of 2-methoxyethanol in several frequency regions ranging from 8.7-500 GHz. The resulting rotational parameters permitted a search for this molecule in Atacama Large Millimeter/submillimeter Array (ALMA) observations toward the massive protocluster NGC 6334I as well as source B of the low-mass protostellar system IRAS 16293-2422. 25 rotational transitions are observed in the ALMA Band 4 data toward NGC 6334I, resulting in the first interstellar detection of 2-methoxyethanol. A column density of $1.3_{-0.9}^{+1.4} \times 10^{17}$ cm$^{-2}$ is derived at an excitation temperature of $143_{-39}^{+31}$ K. However, molecular signal is not observed in the Band 7 data toward IRAS 16293-2422B and an upper limit column density of $2.5 \times 10^{15}$ cm$^{-2}$ is determined. Various possible formation pathways--including radical recombination and insertion reactions--are discussed. We also investigate physical differences between the two interstellar sources that could result in the observed abundance variations.

The extragalactic background light (EBL) is the cumulative radiation outside the Milky Way. The determination of its corresponding primary emitting sources as well as its total energy level across the entire electromagnetic spectrum has profound implications for both cosmology and galaxy formation. However, the detailed origin of the EBL at far-infrared wavelengths, particularly those close to the peak of the cosmic infrared background, remains unclear. Here we report the results of our ongoing SCUBA-2 450 $\mu$m survey of 10 massive galaxy cluster fields. By exploiting the strong gravitational lensing offered by these clusters, we obtain significant counts down to an unprecedented depth of $\sim$0.1 mJy at this wavelength, about ten times deeper than that reached by any other previous survey. The cumulative energy density based on the counts is 138.1$^{+23.9}_{-19.3}$ Jy degree$^{-2}$, or 0.45$^{+0.08}_{-0.06}$ MJy sr$^{-1}$. Comparing our measurements to those made by the COBE and Planck satellites, we find that at this flux density level, the 450 $\mu$m EBL is entirely resolved by our SCUBA-2 observations. Thus, we find for the first time that discrete sources produce fully to the 450 $\mu$m EBL, and that about half of it comes from sources with sub-mJy flux densities. Our deep number counts provide strong constraints on galaxy formation models.

Solar energetic particle (SEP) events are one of the most crucial aspects of space weather that require continuous monitoring and forecasting. Their prediction depends on various factors including source eruptions. In the present work, we use the Geostationary Solar Energetic Particle (GSEP) data set covering solar cycles 22, 23, and 24. We develop a framework using time series-based machine learning (ML) models with the aim of developing robust short-term forecasts by classifying SEP events. For this purpose, we introduce an ensemble learning approach that merges the results from univariate time series of three proton channels (10 MeV, 50 MeV, and 100 MeV) and the long band X-ray flux channel from the Geostationary Operational Environmental Satellite (GOES) missions and analyze their performance. We consider three models, namely, time series forest (TSF), supervised time series forest (STSF) and bag of SFA symbols (BOSS). Our study also focuses on understanding and developing confidence in the predictive capabilities of our models. Therefore, we utilize multiple evaluation techniques and metrics. Based on that, we find STSF to perform well in all scenarios. The summary of metrics for the STSF model is as follows: AUC = 0.981; F1-score = 0.960; TSS = 0.919; HSS = 0.920; GSS = 0.852; and MCC = 0.920. The Brier score loss of the STSF model is 0.077. This work lays the foundation for building near-real-time (NRT) short-term SEP event predictions using robust ML methods.

We present the results of a comparison between different methods to estimate the power of relativistic jets from active galactic nuclei (AGN). We selected a sample of 32 objects (21 flat-spectrum radio quasars, 7 BL Lacertae Objects, 2 misaligned AGN, and 2 changing-look AGN) from the Very Large Baseline Array (VLBA) observations at 43 GHz of the Boston University blazar program. We then calculated the total, radiative, and kinetic jet power from both radio and high-energy gamma-ray observations, and compare the values. We found an excellent agreement between the radiative power calculated by using the Blandford and K\"onigl model with 37 or 43 GHz data, and the values derived from the high-energy $\gamma-$ray luminosity. The agreement is still acceptable if 15 GHz data are used, although with a larger dispersion, but it improves if we use a constant fraction of the $\gamma-$ray luminosity. We found a good agreement also for the kinetic power calculated with Blandford and K\"onigl model with 15 GHz data, and the value from the extended radio emission. We also propose some easy-to-use equations to estimate the jet power.

No matter how fascinating and exotic other terrestrial planets are revealed to be, nothing generates more excitement than announcements regarding their habitability. From the observation of Mars to present-day efforts toward Venus and the characterization of exoplanets, the search for life, or at least environments that could accommodate life, has been a major drive for space exploration. So far, we have found no other unquestionably habitable world besides Earth. The conditions of the habitability of terrestrial planets have proved elusive, as surface conditions depend on the complex interplay of many processes throughout the evolution of a planet. Here, we review how the interior of a rocky planet can drive the evolution of surface conditions and the atmosphere. Instead of listing criteria assumed to be critical for life, we discuss how the bulk-silicate planet can affect the onset, continuation and cessation of habitability. We then consider how it can be observed and current efforts towards this end.

The reflection of X-ray radiation produced near a compact object from its stellar companion contributes to the orbital variability of polarization in X-ray binaries. The X-rays are reflected mainly via Thomson scattering resulting in a high polarization. The orbital variability of the polarization strongly depends on the inclination and the orbital parameters allowing us to constrain them. To explore this phenomenon, we present analytical single-scattering models for the polarized reflection. We find that while diluted by the direct emission, the reflection can produce a polarization degree of about 1$\%$ in the case of a large reflection albedo. We fitted the orbital variations of the X-ray polarization observed by the Imaging X-ray Polarimetry Explorer from an accreting weakly magnetized neutron star `clocked burster' GS 1826$-$238 and found that the amplitude of the variations is too large to be primarily caused by the companion star. The polarized reflection is more significant if the compact object is obscured from the observer, and thus it should be more easily observable in certain high-inclination targets.

We investigate the Wavelet Scattering Transform (WST) as a tool for the study of Primordial non-Gaussianity (PNG) in Large Scale Structure (LSS), and compare its performance with that achievable via a joint analysis with power spectrum and bispectrum (P+B). We consider the three main primordial bispectrum shapes - local, equilateral and orthogonal - and produce Fisher forecast for the corresponding fNL amplitude parameters, jointly with standard cosmological parameters. We analyze simulations from the publicly available "Quijote" and "Quijote-png" N-body suites, studying both the dark matter and halo fields. We find that the WST outperforms the power spectrum alone on all parameters, both on the fNL's and on cosmological ones. In particular, on fNL_loc for halos, the improvement is about 27%. When B is combined with P, halo constraints from WST are weaker for fNL_loc (at ~ 15% level), but stronger for fNL_eq (~ 25%) and fNL_ortho (~ 28%). Our results show that WST, both alone and in combination with P+B, can improve the extraction of information on PNG from LSS data over the one attainable by a standard P+B analysis. Moreover, we identify a class of WST in which the origin of the extra information on PNG can be cleanly isolated.

Shock waves are excited by coronal mass ejections (CMEs) and large-scale extreme-ultraviolet (EUV) wave fronts and can result in low-frequency radio emission under certain coronal conditions. In this work, we investigate a moving source of low-frequency radio emission as a CME and an associated EUV wave front move along a channel of a lower density, magnetic field, and Alfv\'en speed in the solar corona. Observations from the Atmospheric Imaging Assembly on board the Solar Dynamics Observatory, the Nan\c{c}ay Radio Heliograph (NRH), and the Irish Low Frequency Array(I-LOFAR) were analysed. Differential emission measure maps were generated to determine densities and Alfv\'en maps, and the kinematics of the EUV wave front was tracked using CorPITA. The radio sources' positions and velocity were calculated from NRH images and I-LOFAR dynamic spectra. The EUV wave expanded radially with a uniform velocity of $\sim$ 500 km s$^{-1}$. However, the radio source was observed to be deflected and appeared to move along a channel of a lower Alfv\'en speed, abruptly slowing from 1700 km s$^{-1}$ to 250 km s$^{-1}$ as it entered a quiet-Sun region. A shock wave with an apparent radial velocity of > 420 km s$^{-1}$ was determined from the drift rate of the associated Type II radio burst. The apparent motion of the radio source may have resulted from a wave front moving along a coronal wave guide or by different points along the wave front emitting at locations with favourable conditions for shock formation.

One X-ray short burst accompanied by a Galactic fast radio burst has been detected so far, and its X-ray cut-off energy was significantly higher than that of other X-ray short bursts. Such X-ray bursts are thought to be emitted from a fireball in the magnetosphere of magnetars. If a fireball is formed around a magnetic pole, it expands and is accelerated under its radiation pressure, later producing photon emission and plasma outflow, which can be responsible for radio bursts. We numerically study the radiative acceleration of this outflowing fireball, which consists of electron-positron pairs and radiation, and obtain the spectrum of the escaped X-ray photons. We consistently take into account cyclotron resonant scattering, which enhances the scattering cross section resulting in a strong radiative force and high optical depth. Our results show that similar spectra to the observed X-ray spectrum in the Galactic fast radio burst are realized in outflowing fireballs and that the plasma outflow is simultaneously accelerated to a high Lorentz factor owing to cyclotron resonant scattering.

We update the ephemerides of 16 transiting exoplanets using our ground-based observations, new TESS data, and previously published observations including those of amateur astronomers. All these light curves were modeled by making use of a set of quantitative criteria with the EXOFAST code to obtain mid-transit times. We searched for statistically significant secular and/or periodic trends in the mid-transit times. We found that the timing data are well modeled by a linear ephemeris for all systems except for XO-2 b, for which we detect an orbital decay with the rate of -12.95 $\pm$ 1.85 ms/yr that can be confirmed with future observations. We also detect a hint of potential periodic variations in the TTV data of HAT-P-13 b which also requires confirmation with further precise observations.

The Unified Model of Active Galactic Nuclei (UMAGN) is a comprehensive theoretical framework aimed at elucidating the diverse observations of AGN, encompassing quasars and Seyfert galaxies. The Model attributes the observed variations to different orientations of a surrounding matter disk around a supermassive black hole (SMBH), with the primary factor influencing observational diversity being the alignment of the AGN with the observer's line of sight. We present a comprehensive overview of the observational evidence, empirical and theoretical research, tracing key milestones that led to a unified perspective. We encapsulate the scientific journey culminating in the proposal of the UMAGN, including insights into the accretion disk, torus, and relativistic jet, and emphasize the properties of objects within UMAGN. Additionally, we underscore recent progress in multimessenger research involving electromagnetic waves, gravitational waves, astroparticles, and neutrinos, notably through collaborations such as the Event Horizon Telescope, LIGO/Virgo, IceCube, Pierre Auger, and KM3Net. We argue that these advancements present opportunities to enhance and refine UMAGN, contributing to a deeper understanding of AGNs and their implications for the formation and evolution of galaxies. The convergence of observational and theoretical research, coupled with emerging multimessenger techniques, paves the way for further strides in comprehending these enigmatic cosmic phenomena.

In previous work we have shown that the dipole in the low redshift supernovae of the Pantheon+SH0ES data does not agree with the one inferred from the velocity of the solar system as obtained from CMB data. We interpreted this as the presence of significant bulk velocities. In this paper we study the monopole, dipole and quadrupole in the Pantheon+SH0ES data. We find that in addition to the dipole also both, the monopole and the quadrupole are detected with high significance. They are of similar amplitudes as the bulk flow. While the monopole is only significant at very low redshift, the quadrupole even increases with redshift.

The formation history of Jupiter has been of interest due to its ability to shape the solar system's history. Yet little attention has been paid to the formation and growth of Saturn and the other giant planets. Here, we explore the implications of the simplest disc and pebble accretion model with steady-state accretion on the formation of giant planets in the solar system through N-body simulations. We conducted a statistical survey of different disc parameters and initial conditions of the protoplanetary disc to establish which combination best reproduces the present outer solar system. We examined the effect of the initial planetesimal disc mass, the number of planetesimals and their size-frequency distribution slope, pebble accretion prescription, and sticking efficiency on the likelihood of forming gas giants and their orbital distribution. The results reveal that the accretion sticking efficiency is the most sensitive parameter for controlling the final masses and number of giant planets. We have been unable to replicate the formation of all three types of giant planets in the solar system in a single simulation. The probability distribution of the final location of the giant planets is approximately constant in $\log r$, suggesting there is a slight preference for formation closer to the Sun but no preference for more massive planets to form closer. The eccentricity distribution has a higher mean for more massive planets, indicating that systems with more massive planets are more violent. The formation timescales of the cores of the gas giants are distinct, suggesting that they formed sequentially.

We report the discovery of ten new pulsars in the globular cluster Terzan 5 as part of the Transients and Pulsars with MeerKAT (TRAPUM) Large Survey Project. We observed Terzan 5 at L-band (856--1712 MHz) with the MeerKAT radio telescope for four hours on two epochs, and performed acceleration searches of 45 out of 288 tied-array beams covering the core of the cluster. We obtained phase-connected timing solutions for nine discoveries, covering nearly two decades of archival observations from the Green Bank Telescope for all but one. Highlights include PSR J1748$-$2446ao which is an eccentric ($e = 0.32$) wide-orbit (orbital period $P_{\rm b} = 57.55$ d) system. We were able to measure the rate of advance of periastron ($\dot{\omega}$) for this system allowing us to determine a total mass of $3.17 \pm \, 0.02\, \rm M_{\odot}$. With a minimum companion mass ($M_{\rm c}$) of $\sim 0.8\, \rm M_{\odot}$, PSR J1748$-$2446ao is a candidate double neutron star (DNS) system. If confirmed to be a DNS, it would be the fastest spinning pulsar ($P = 2.27$ ms) and the longest orbital period measured for any known DNS system. PSR J1748$-$2446ap has the second highest eccentricity for any recycled pulsar ($e \sim 0.905$) and for this system we can measure the total mass ($1.997 \pm 0.006\, \rm M_{\odot}$) and also estimate the individual pulsar and companion masses. PSR J1748$-$2446ar is an eclipsing redback (minimum $M_{\rm c} \sim 0.34\, \rm M_{\odot}$) system whose properties confirm it to be the counterpart to a previously published source identified in radio and X-ray imaging. With these discoveries, the total number of confirmed pulsars in Terzan 5 is 49, the highest for any globular cluster so far. These discoveries further enhance the rich set of pulsars known in Terzan 5 and provide scope for a deeper understanding of binary stellar evolution, cluster dynamics and ensemble population studies.

We observe systematic profile changes in the visible pulsar of the compact double neutron star system PSR~J1946+2052 using observations with the Five-hundred-meter Aperture Spherical radio Telescope (FAST). The interpulse of PSR~J1946+2052 changed from single-peak to double-peak shape from 2018 to 2021. We attribute this evolution as the result of the relativistic spin precession of the pulsar. With the high sensitivity of FAST, we also measure significant polarization for the first time, allowing us to model this with the precessional rotating vector model. Assuming, to the first order, a circular hollow-cone-like emission beam pattern and taking the validity of general relativity, we derive the binary's orbital inclination angle (${63^\circ}^{+5^\circ}_{-3^\circ}$) and pulsar's spin geometry. Pulsar's spin vector and the orbital angular momentum vector are found to be only slightly misaligned (${0.21^\circ}^{+0.28^\circ}_{-0.10^\circ}$).The quoted uncertainties do not reflect the systematic uncertainties introduced by our model assumptions. By simulating future observations of profile and polarization evolution, we estimate that we could constrain the precession rate within a $43\%$ uncertainty in 9 years. Hence, we suggest that the system's profile evolution could be combined with precise pulsar timing to test general relativity in the future.

In the fall of 2019, during the in-flight calibration phase of the SRG observatory, the onboard eROSITA and Mikhail Pavlinsky ART-XC telescopes carried out a series of observations of PG 1634+706 - one of the most luminous (an X-ray luminosity $\sim 10^{46}$ erg/s) quasars in the Universe at $z<2$. Approximately at the same dates this quasar was also observed by the XMM-Newton observatory. Although the object had already been repeatedly studied in X-rays previously, its new observations allowed its energy spectrum to be measured more accurately in the wide range $1-30$ keV (in the quasar rest frame). Its spectrum can be described by a two-component model that consists of a power-law continuum with a slope $\Gamma\approx 1.9$ and a broadened iron emission line at an energy of about 6.4 keV. The X-ray variability of the quasar was also investigated. On time scales of the order of several hours (here and below, in the source rest frame) the X-ray luminosity does not exhibit a statistically significant variability. However, it changed noticeably from observation to observation in the fall of 2019, having increased approximately by a factor of 1.5 in 25 days. A comparison of the new SRG and XMM-Newton measurements with the previous measurements of other X-ray observatories has shown that in the entire 17-year history of observations of the quasar PG 1634+706 its X-ray luminosity has varied by no more than a factor of 2.5, while the variations on time scales of several weeks and several years are comparable in amplitude.

This study presents bounds on transient variations of fundamental constants, with typical timescales ranging from minutes to months, using clocks in space. The underlying phenomenology describing such transient variations relies on models for Dark Matter (DM) which suggest possible encounters of macroscopic compact objects with the Earth, due to the motion of the solar system in the galactic halo. If such compact objects possess an effective feeble interaction with the ordinary matter beyond the gravitational one, it may result in effective transient variations of fundamental constants. Such variations leave signatures on clocks onboard GNSS satellites. In this paper, we introduce a phenomenological study dedicated to the search for such DM transient objects using the network of passive hydrogen masers (H-Masers) onboard Galileo satellites. We first model the signature of transient variations of fundamental constants as a frequency modulation in the difference between two satellite clocks, considering the satellite trajectories relative to the transient event. Then, we present first results based on a fast analysis method, the maximum reach analysis. The main result is a significant extension of the discovery range for DM transients, with a sensitivity never achieved before. We investigate indeed the range of transient sizes from $10^5$ to $10^9$ kilometres, which, apart from indirect and model-dependent non-transient effects, has never been explored previously.

Mutual gravitational interactions between the five major Uranian satellites raise small quasi-periodic fluctuations on their orbital elements. At the same time, tidal interactions between the satellites and the planet induce a slow outward drift of the orbits, while damping the eccentricities and the inclinations. In this paper, we revisit the current and near past evolution of this system using a N-body integrator, including spin evolution and tidal dissipation with the weak friction model. We update the secular eigenmodes of the system and show that it is unlikely that any of the main satellites were recently captured into a high obliquity Cassini state. We rather expect that the Uranian satellites are in a low obliquity Cassini state and compute their values. We also show that the current eccentricities of the satellites are not forced, and estimate the free eccentricities and inclinations. We constrain the quality factor of Uranus to be $Q_U = (8.6 \pm 2.9)\times10^3$, and that of the satellites to be $Q_S \sim 500$. We find that the system most likely encountered the 5/3 mean motion resonance between Ariel and Umbriel in the past, at about ($0.7\pm0.2$) Gyr ago. We additionally determine the eccentricities and inclinations of all satellites just after the resonance passage that comply with the current system. We finally show that, from the crossing of the 5/3 MMR to the present, the evolution of the system is mostly peaceful and dominated by tides raised on Uranus by the satellites.

At present, the main satellites of Uranus are not involved in any low order mean motion resonance (MMR). However, owing to tides raised in the planet, Ariel and Umbriel most likely crossed the 5/3 MMR in the past. Previous studies on this resonance passage relied on limited time-consuming N-body simulations or simplified models focusing solely on the effects of the eccentricity or the inclination. In this paper, we aim to provide a more comprehensive view on how the system evaded capture in the 5/3 MMR. For that purpose, we developed a secular resonant two-satellite model with low eccentricities and low inclinations, including tides using the weak friction model. By performing a large number of numerical simulations, we show that capture in the 5/3 MMR is certain if the initial eccentricities of Ariel, $e_1$, and Umbriel, $e_2$, are related through $(e_1^2 + e_2^2)^{1/2} < 0.007$. Moreover, we observe that the eccentricity of Ariel is the key variable to evade the 5/3 MMR with a high probability. We determine that for $e_1 > 0.015$ and $e_2 < 0.01$, the system avoids capture in at least 60\% of the cases. We also show that, to replicate the currently observed system, the initial inclinations of Ariel and Umbriel must lay within $I_1 \leq 0.05^{\circ}$ and $0.06^{\circ} \leq I_2 \leq 0.11^{\circ}$, respectively. We checked these results using a complete N-body model with the five main satellites and did not observe any significant differences.

Of the hundreds of $z\gtrsim6$ quasars discovered to date, only one is known to be gravitationally lensed, despite the high lensing optical depth expected at $z\gtrsim6$. High-redshift quasars are typically identified in large-scale surveys by applying strict photometric selection criteria, in particular by imposing non-detections in bands blueward of the Lyman-$\alpha$ line. Such procedures by design prohibit the discovery of lensed quasars, as the lensing foreground galaxy would contaminate the photometry of the quasar. We present a novel quasar selection methodology, applying contrastive learning (an unsupervised machine learning technique) to Dark Energy Survey imaging data. We describe the use of this technique to train a neural network which isolates an 'island' of 11 sources, of which 7 are known $z\sim6$ quasars. Of the remaining four, three are newly discovered quasars (J0109-5424, $z=6.07$; J0122-4609, $z=5.99$; J0603-3923, $z=5.94$), as confirmed by follow-up Gemini-South/GMOS and archival NTT/EFOSC2 spectroscopy, implying a 91 per cent efficiency for our novel selection method; the final object on the island is a brown dwarf. In one case (J0109-5424), emission below the Lyman limit unambiguously indicates the presence of a foreground source, though high-resolution optical/near-infrared imaging is still needed to confirm the quasar's lensed (multiply-imaged) nature. Detection in the g band has led this quasar to escape selection by traditional colour cuts. Our findings demonstrate that machine learning techniques can thus play a key role in unveiling populations of quasars missed by traditional methods.

We forecast constraints on cosmological parameters enabled by three surveys conducted with SPT-3G, the third-generation camera on the South Pole Telescope. The surveys cover separate regions of 1500, 2650, and 6000 ${\rm deg}^{2}$ to different depths, in total observing 25% of the sky. These regions will be measured to white noise levels of roughly 2.5, 9, and 12 $\mu{\rm K-arcmin}$, respectively, in CMB temperature units at 150 GHz by the end of 2024. The survey also includes measurements at 95 and 220 GHz, which have noise levels a factor of ~1.2 and 3.5 times higher than 150 GHz, respectively, with each band having a polarization noise level ~$\sqrt{\text{2}}$ times higher than the temperature noise. We use a novel approach to obtain the covariance matrices for jointly and optimally estimated gravitational lensing potential bandpowers and unlensed CMB temperature and polarization bandpowers. We demonstrate the ability to test the $\Lambda{\rm CDM}$ model via the consistency of cosmological parameters constrained independently from SPT-3G and Planck data, and consider the improvement in constraints on $\Lambda{\rm CDM}$ extension parameters from a joint analysis of SPT-3G and Planck data. The $\Lambda{\rm CDM}$ cosmological parameters are typically constrained with uncertainties up to ~2 times smaller with SPT-3G data, compared to Planck, with the two data sets measuring significantly different angular scales and polarization levels, providing additional tests of the standard cosmological model.

Light rays passing very close to black holes may wind several times before escaping. For any given electromagnetic source around the black hole, a distant observer would thus observe two infinite sequences of images on either side of the black hole. These images are generated by light rays performing an increasing numbers of loops. The strong deflection limit provides a simple analytic formalism to describe such higher order images for spherically symmetric metrics, while for axially symmetric black holes one typically resorts to numerical approaches. Here we present the leading order perturbation to higher order images when the black hole spin is turned on. We show that the images slide around the black hole shadow as an effect of space-time dragging. We derive analytical formulae for their shifts and the perturbation of their time delays. We also discuss how such simple analytical formulae for images by Kerr black holes can be of great help in many applications.

Unimodular gravity became an object of increasing interest in the late $80$-ties and was recently used in primordial Universe modeling with cosmological constant, in the context of the Brans-Dicke gravity including scalar field. In the present article we investigate the possibility of imposing the unimodular condition within the $5$-dimensional Kaluza-Klein theory including the scalar field. The variational principle is formulated in $5$ dimensions first, and dimensional reduction is applied to the resulting set of equations. A cosmological model based on these equations is then presented and discussed.

In this work we analyze plasma and magnetic field data provided by the Parker Solar Probe (\emph{PSP}) and Solar Orbiter (\emph{SO}) missions to investigate the radial evolution of the heating of Alfv\'enic slow wind (ASW) by imbalanced Alfv\'en-Wave (AW) turbulent fluctuations from 0.06 au to 1 au. in our analysis we focus on slow solar-wind intervals with highly imbalanced and incompressible turbulence (i.e., magnetic compressibility $C_B=\delta B/B\leq 0.25$, plasma compressibility $C_n=\delta n/n\leq 0.25$ and normalized cross-helicity $\sigma_c\geq 0.65$). First, we estimate the AW turbulent dissipation rate from the wave energy equation and find that the radial profile trend is similar to the proton heating rate. Second, we find that the scaling of the empirical AW turbulent dissipation rate $Q_W$ obtained from the wave energy equation matches the scaling from the phenomenological AW turbulent dissipation rate $Q_{\rm CH09}$ (with $Q_{\rm CH09}\simeq 1.55 Q_W$) derived by~\cite{chandran09} based on the model of reflection-driven turbulence. Our results suggest that, as in the fast solar wind, AW turbulence plays a major role in the ion heating that occurs in incompressible slow-wind streams.

The first direct detection of gravitational waves from binary neutron stars on the 17th of August, 2017, (GW170817) heralded the arrival of a new messenger for probing neutron star astrophysics and provided the first constraints on neutron star equation of state from gravitational wave observations. Significant computational effort was expended to obtain these first results and therefore, as observations of binary neutron star coalescence become more routine in the coming observing runs, there is a need to improve the analysis speed and flexibility. Here, we present a rapid approach for inferring the neutron star equation of state based on Normalising Flows. As a demonstration, using the same input data, our approach, ASTREOS, produces results consistent with those presented by the LIGO-Virgo collaboration but requires < 1 sec to generate neutron star equation of state confidence intervals. Furthermore, ASTREOS allows for non-parametric equation of state inference. This rapid analysis will not only facilitate neutron star equation of state studies but can potentially enhance future alerts for electromagnetic follow-up observations of binary neutron star mergers.

We present a new numerical tool designed to probe the dense layers of neutron star crusts. It is based on the Time-Dependent Hartree-Fock-Bogoliubov theory with generalized Skyrme nuclear energy density functionals, such as the Brussels-Montreal ones. We use it to study the time evolution of a nucleus accelerating through superfluid neutron medium in the inner crust of a neutron star. We extract an effective mass in the low velocity limit. We observe a threshold velocity and specify mechanisms of dissipation: phonon emission, Cooper pairs breaking, and vortex rings creation. The microscopic effects we study have impact on neutron star. Moreover, the mechanisms, we described, are general and apply also to other fermionic superfluid mixtures like liquid helium, or ultracold gases.

Vacuum decay and symmetry breaking play an important role in the fundamental structure of the matter and the evolution of the universe. In this work we study how the purely classical effect of accretion of fundamental fields onto black holes can lead to shells of symmetry restoration in the midst of a symmetry broken phase. We also show how it can catalyze vacuum decay, forming a bubble that expands asymptotically at the speed of light. These effects offer an alternative, purely classical mechanism to quantum tunnelling for seeding phase transitions in the universe.

I present an introduction and topical review on axions as a dark matter candidate. Emphasis is placed on issues surrounding the cosmology of axion dark matter that are relevant for present-day searches, including: early-Universe production mechanisms, predictions of the axion mass, bounds on axion properties derived from cosmological data, as well as the direct and indirect detection of relic axion populations.

We consider a double polytropic cosmological fluid and demonstrate that, when one constituent resembles a bare cosmological constant while the other emulates a generalized Chaplygin gas, a good description of the Universe's large-scale dynamics is obtained. In particular, our double polytropic reduces to the Murnaghan equation of state, whose applications are already well established in solid state physics and classical thermodynamics. Intriguingly, our model approximates the conventional $\Lambda$CDM paradigm while reproducing the collective effects of logotropic and generalized Chaplygin fluids across different regimes. To check the goodness of our fluid description, we analyze first order density perturbations, refining our model through various orders of approximation, utilizing $\sigma_8$ data alongside other cosmological data sets. Encouraging results suggest that our model, based on the Murnaghan equation of state, outperforms the standard cosmological background within specific approximate regimes and, on the whole, surpasses the standard phenomenological reconstruction of dark energy.