Papers for Thursday, Sep 02 2021

The epoch of Cosmic Dawn, when the first stars and galaxies were born, is widely considered the final frontier of observational cosmology today. Mapping the period between Cosmic Dawn and the present-day provides access to more than 90% of the baryonic (normal) matter in the Universe, and unlocks several thousand times more Fourier modes of information than available in today's cosmological surveys. We review the progress in modelling baryonic gas observations as tracers of the cosmological large-scale structure from Cosmic Dawn to the present day. We illustrate how the description of dark matter haloes can be extended to describe baryonic gas abundances and clustering. This innovative approach allows us to fully utilize our current knowledge of astrophysics to constrain cosmological parameters from future observations. Combined with the information content of multi-messenger probes, this will also elucidate the properties of the first supermassive black holes at Cosmic Dawn. We present a host of fascinating implications for constraining physics beyond the $\Lambda$CDM model, including tests of the theories of inflation and the cosmological principle, the effects of non-standard dark matter, and possible deviations from Einstein's general relativity on the largest scales.

Out-of-focus (OOF) holography can be used to determine aperture deformations of radio telescopes that lead to errors in the phase of the complex aperture distribution. In contrast to traditional holography, OOF can be done without a reference antenna, which has a number of practical advantages. The aim of this work is to develop a standard procedure for OOF holography at the Effelsberg telescope. This includes performing OOF holography observations and the development of a software, the pyoof package, to compute aberrations of the telescope's optical system. Based on the OOF holography method developed by Nikolic et al. (2007a), we adapted the approach to the Effelsberg 100-m telescope in order to determine the aberrations of the aperture phase distribution (phase-error maps). The developed OOF holography software is presented as well as the results from observations performed at Effelsberg. Early results reveal a possible gravitationally-caused residual deformation not contained in the previously existing aperture and pointing model, and hence we propose to make changes to the model to counteract aberrations in the telescope's surface. The OOF holography method (observations and pyoof package) works as expected at the Effelsberg 100-m telescope and is able to validate the good performance of the existing finite element model. Test measurements show that slight improvements of the aperture efficiency and gain elevation dependence are possible with a more extensive OOF holography campaign.

Context. Time delay lensing is a powerful tool to measure the Hubble constant $H_0$. In order to obtain an accurate estimate of $H_0$ from a sample of time delay strong lenses, however, it is necessary to have a very good knowledge of the mass structure of the lens galaxies. Strong lensing data on their own are not sufficient to break the degeneracy between $H_0$ and the lens model parameters, on a single object basis. Aims. The goal of this study is to determine whether it is possible to break the $H_0$-lens structure degeneracy with the statistical combination of a large sample of time-delay lenses, relying purely on strong lensing data (that is, with no stellar kinematics information). Methods. I simulated a set of 100 lenses with doubly imaged quasars and related time delay measurements. I fitted these data with a Bayesian hierarchical method and a flexible model for the lens population, emulating the lens modelling step. Results. The sample of 100 lenses, on its own, provides a measurement of $H_0$ with $3\%$ precision, but with a $-4\%$ bias. However, the addition of prior information on the lens structural parameters from a large sample of lenses with no time delays, such as that considered in Paper I, allows for a $1\%$-level inference. Conclusions. Breaking the $H_0$-lens model degeneracy with lensing data alone is possible, but $1\%$ measurements of $H_0$ require either a number of time delay lenses much larger than 100, or the knowledge of the structural parameter distribution of the lens population from a separate sample of lenses.

Long $\rm \gamma$-ray bursts (GRBs) are produced by the dissipation of ultra-relativistic jets launched by newly-born black holes after a collapse of massive stars. Right after the luminous and highly variable $\gamma$-ray emission, the multi-wavelength afterglow is released by the external dissipation of the jet in circumburst medium. We report the discovery of very bright ($\rm \sim 10$ mag) optical emission $\rm \sim 28$ s after the explosion of the extremely luminous and energetic GRB 210619B located at redshift 1.937. Early multi-filter observations allowed us to witness the end of the shock wave propagation into the GRB ejecta. We observed the spectral transition from a bright reverse to the forward shock emission, demonstrating that the GRB multi-wavelength emission is originated from a narrow and highly magnetised jet propagating into a rarefied interstellar medium. We also find evidence of an additional component of radiation, coming from the jet wings which is able explain the uncorrelated optical/X-ray emission.

We quantify the impact of galaxy formation on dark matter halo shapes using cosmological simulations at redshift $z=0$. The haloes are drawn from the IllustrisTNG project, a suite of magneto-hydrodynamic simulations of galaxies. We focus on haloes of mass $10^{10-14} M_\odot$ from the 50-Mpc (TNG50) and 100-Mpc (TNG100) boxes, and compare them to dark matter-only (DMO) analogues and other simulations e.g. NIHAO and Eagle. We further quantify the prediction uncertainty by varying the baryonic feedback models in a series of smaller 25 Mpc $h^{-1}$ boxes. We find that: (i) galaxy formation results in rounder haloes compared to the DMO simulations, in qualitative agreement with past hydrodynamic models. Haloes of mass $\approx 2\times 10^{12} M_\odot$ are most spherical, with an average minor-to-major axis ratio of $\left< s \right> \approx 0.75$ in the inner halo, an increase of 40 per cent compared to their DMO counterparts. No significant change in halo shape is found for low-mass $10^{10} M_\odot$ haloes; (ii) stronger feedback, e.g. increasing galactic wind speed, reduces the impact of baryons; (iii) the inner halo shape correlates with the stellar mass fraction, which can explain the dependence of halo shapes on different feedback models; (iv) the fiducial and weaker feedback models are most consistent with observational estimates of the Milky Way halo shape. Yet, at fixed halo mass, very diverse and possibly unrealistic feedback models all predict inner halo shapes that are closer to one another than to the DMO results. This implies that a larger observational sample would be required to statistically distinguish between different baryonic prescriptions due to large halo-to-halo variation in halo shapes.

Spectroscopically, TDEs are characterized by broad ( 10$^{4}$ km/s) emission lines and show large diversity as well as different line profiles. After carefully and consistently performing a series of data reduction tasks including host galaxy light subtraction, we present here the first detailed, spectroscopic population study of 16 optical/UV TDEs. We report a time lag between the peaks of the optical light-curves and the peak luminosity of H$\alpha$ spanning between 7 - 45 days. If interpreted as light-echoes, these lags correspond to distances of 2 - 12 x 10$^{16}$ cm, one to two orders of magnitudes larger than the estimated blackbody radii (R$_{\rm BB}$) of the same TDEs and we discuss the possible origin of this surprisingly large discrepancy. We also report time lags for the peak luminosity of He I $\lambda$5876 line; smaller than the ones of H$\alpha$ for H TDEs and similar or larger for N III Bowen TDEs. We report that N III Bowen TDEs have lower H$\alpha$ velocity widths compared to the rest of the TDEs in our sample and we also find that a strong X-ray to optical ratio might imply weakening of the line widths. Furthermore, we study the evolution of line luminosities and ratios with respect to their radii (R$_{\rm BB}$) and temperatures (T$_{\rm BB}$). We find a linear relationship between H$\alpha$ luminosity and the R$_{\rm BB}$ and potentially an inverse power-law relation with T$_{\rm BB}$ leading to weaker H$\alpha$ emission for T$_{\rm BB}$ $\geq$ 25000 K. The He II/He I ratio becomes large at the same temperatures possibly pointing to an ionization effect. The He II/H$\alpha$ ratio becomes larger as the photospheric radius recedes, implying a stratified photosphere where Helium lies deeper than Hydrogen. We suggest that the large diversity of the spectroscopic features seen in TDEs along with their X-ray properties, can potentially be attributed to viewing angle effects.

A defining prediction of the cold dark matter (CDM) cosmological model is the existence of a very large population of low-mass haloes. This population is absent in models in which the dark matter particle is warm (WDM). These alternatives can, in principle, be distinguished observationally because halos along the line-of-sight can perturb galaxy-galaxy strong gravitational lenses. Furthermore, the WDM particle mass could be deduced because the cut-off in their halo mass function depends on the mass of the particle. We systematically explore the detectability of low-mass haloes in WDM models by simulating and fitting mock lensed images. Contrary to previous studies, we find that halos are harder to detect when they are either behind or in front of the lens. Furthermore, we find that the perturbing effect of haloes increases with their concentration: detectable haloes are systematically high-concentration haloes, and accounting for the scatter in the mass-concentration relation boosts the expected number of detections by as much as an order of magnitude. Haloes have lower concentration for lower particle masses and this further suppresses the number of detectable haloes beyond the reduction arising from the lower halo abundances alone. Taking these effects into account can make lensing constraints on the value of the mass function cut-off at least an order of magnitude more stringent than previously appreciated.

Big bang nucleosynthesis (BBN) is the standard model theory for the production of the light nuclides during the early stages of the universe, taking place for a period of about 20 minutes after the big bang. Deuterium production, in particular, is highly sensitive to the primordial baryon density and the number of neutrino species, and its abundance serves as a sensitive test for the conditions in the early universe. The comparison of observed deuterium abundances with predicted ones requires reliable knowledge of the relevant thermonuclear reaction rates, and their corresponding uncertainties. Recent observations reported the primordial deuterium abundance with percent accuracy, but some theoretical predictions based on BBN are at tension with the measured values because of uncertainties in the cross section of the deuterium-burning reactions. In this work, we analyze the S-factor of the D(p,$\gamma$)$^3$He reaction using a hierarchical Bayesian model. We take into account the results of eleven experiments, spanning the period of 1955--2021; more than any other study. We also present results for two different fitting functions, a two-parameter function based on microscopic nuclear theory and a four-parameter polynomial. Our recommended reaction rates have a 2.2\% uncertainty at $0.8$~GK, which is the temperature most important for deuterium BBN. Differences between our rates and previous results are discussed.

The efficiency of regolith production is key in understanding the properties of airless surfaces. Debris aprons, of fillets, around rocks are an ubiquitous morphology on many surfaces without atmosphere, which origin and evolution are largely unknown. Here we show that fillet originates from the juxtaposed rock under abrasion and that rocks of different cohesion have fillets with distinct morphological evolution. Thus, a fillet around a rock allows to disentangle rock cohesion from its surface exposure age. By combing topographic diffusion modeling with images of blocks of known age on the Moon we find abrasion rates for cm-sized boulders similar to regional rates (0.2 mm/Myr), whereas for 10-m sized blocks the rate is two order of magnitude higher (20 mm/Myr). Rates for instances of rocks of higher strength are reduced by ~50%. Fillets around lunar rocks are consistent with abrasion by isotropic micrometeoroid bombardment.

The irradiance of the Sun is modulated on all time scales. Small-scale solar magnetic elements composed of quiet-Sun network and active region plages contribute to this modulation on solar cycle time scales. The evaluation of their contrast as a function of their magnetic field strength is an important constraint for models of solar irradiance variation. In this thesis, we improve and extend the results of earlier contrast studies taking advantage of the high resolution data delivered by the balloon-borne solar observatory Sunrise during its science flights in 2009 and 2013. In addition to the high quality of the Sunrise data, access to the UV wavelengths makes these data unique in this kind of study. In the first two parts of this thesis, we find that the contrast variation with the magnetic field strength agrees qualitatively with radiative magnetohydrodynamic simulations of the quiet Sun. In particular, at a continuum wavelength around 500\,nm, the contrast does not decrease at large field strengths, indicative of having resolved these features. In addition, the contrast at all wavelengths in the quiet-Sun network is found to be higher than in solar plage, and the difference decreases with atmospheric height in accordance with empirical models of flux tubes. In the last part of this thesis, we make use of solar limb data observed by Sunrise to evaluate the contribution of stray light to the total point spread function of the telescope.

Ultra-hot Jupiters are highly irradiated gas giants with equilibrium temperatures typically higher than 2000K. Atmospheric studies of these planets have shown that their transmission spectra are rich in metal lines, with some of these metals being ionised due to the extreme temperatures. Here, we use two transit observations of WASP-76b obtained with the CARMENES spectrograph to study the atmosphere of this planet using high-resolution transmission spectroscopy. Taking advantage of the two channels and the coverage of the red and near-infrared wavelength ranges by CARMENES, we focus our analysis on the study of the CaII infrared triplet (IRT) at 8500A and the HeI triplet at 10830A. We present the discovery of the CaII IRT at 7$\sigma$ in the atmosphere of WASP-76b using the cross-correlation technique, which is consistent with previous detections of the CaII H&K lines in the same planet, and with the atmospheric studies of other ultra-hot Jupiters reported to date. The low mass density of the planet, and our calculations of the XUV irradiation received by the exoplanet, show that this planet is a potential candidate to have a HeI evaporating envelope and, therefore, we performed further investigations focussed on this aspect. The transmission spectrum around the HeI triplet shows a broad and red-shifted absorption signal in both transit observations. However, due to the strong telluric contamination around the HeI lines and the relatively low signal-to-noise ratio of the observations, we are not able to unambiguously conclude if the absorption is due to the presence of helium in the atmosphere of WASP-76b, and we consider the result to be only an upper limit. Finally, we revisit the transmission spectrum around other lines such as NaI, LiI, H$\alpha$, and KI. The upper limits reported here for these lines are consistent with previous studies.

Despite extensive searches and the relative proximity of solar system objects (SSOS) to Earth, many remain undiscovered and there is still much to learn about their properties and interactions. This work is the first in a series dedicated to detecting and analyzing SSOs in the all-sky NOIRLab Source Catalog (NSC). We search the first data release of the NSC with CANFind, a Computationally Automated NSC tracklet Finder. NSC DR1 contains 34 billion measurements of 2.9 billion unique objects, which CANFind categorizes as belonging to "stationary" (distant stars, galaxies) or moving (SSOs) objects via an iterative clustering method. Detections of stationary bodies for proper motion (mu) less than 2.5"/hr (0.017 degrees/day) are identified and analyzed separately. Remaining detections belonging to hi-mu objects are clustered together over single nights to form "tracklets". Each tracklet contains detections of an individual moving object, and is validated based on spatial linearity and motion through time. Proper motions are then calculated and used to connect tracklets and other unassociated measurements over multiple nights by predicting their locations at common times forming "tracks". This method extracted 527,055 tracklets from NSC DR1 in an area covering 29,971 square degrees of the sky. The data show distinct groups of objects with similar observed mu in ecliptic coordinates, namely Main Belt Asteroids, Jupiter Trojans, and Kuiper Belt Objects. Apparent magnitudes range from 10-25 mag in the ugrizY and VR bands. Color-color diagrams show a bimodality of tracklets between primarily carbonaceous and siliceous groups, supporting prior studies.

In this paper, we have selected a sample of 64 teraelectronvolt blazars, with redshift, from those classified in the fourth Fermi Large Area Telescope source catalog\footnote{\url{https://fermi.gsfc.nasa.gov/ssc/data/access/lat/8yr_catalog/}}. We have obtained the values of the relevant physical parameters by performing a log-parabolic fitting of the average-state multiwavelength spectral energy distributions. We estimate the range of the radiation zone parameters, such as the Doppler factor (${D}$), the magnetic field strength ($B$), the radiative zone radius ($R$) and the peak Lorentz factor (${\gamma _{\rm p}}$) of nonthermal electrons. Here, we show that (1) there is a strong linear positive correlation between the intrinsic synchrotron peak frequency and the intrinsic inverse Compton scattering (ICs) peak frequency among different types of blazars; (2) if radio bands are excluded, the spectral index of each band is negatively correlated with the intrinsic peak frequency; (3) there is a strong linear negative correlation between the curvature at the peak and the intrinsic peak frequency of the synchrotron bump, and a weak positive correlation between the curvature at the peak and the intrinsic peak frequency of the ICs bump; (4) there is a strong linear positive correlation between the intrinsic ICs peak luminosity and intrinsic $\gamma$-ray luminosity and between the intrinsic ICs peak frequency and peak Lorentz factor; (5) there is a strong negative linear correlation between $\rm log{\;B}$ and $\rm log{\;\gamma_{p}}$; and (6) there is no correlation between $\rm log{\;R}$ and $\rm log{\;\gamma_{p}}$.

In this study we demonstrate that stellar masses of galaxies (Mstar) are universally correlated through a double power law function with the product of the dynamical velocities (Ve) and sizes to one-fourth power (Re^0.25) of galaxies, both measured at the effective radii. The product VeRe^0.25 represents the fourth root of the total binding energies within effective radii of galaxies. This stellar mass-binding energy correlation has an observed scatter of 0.14 dex in log(VeRe^0.25) and 0.46 dex in log(Mstar). It holds for a variety of galaxy types over a stellar mass range of nine orders of magnitude, with little evolution over cosmic time. A toy model of self-regulation between binding energies and supernovae feedback is shown to be able to reproduce the observed slopes, but the underlying physical mechanisms are still unclear. The correlation can be a potential distance estimator with an uncertainty of 0.2 dex independent of the galaxy type.

We explore a possible time variation of the fine structure constant ($\alpha \equiv e^2/\hbar c$) using the Sunyaev-Zel'dovich effect measurements of galaxy clusters along with their X-ray observations. Specifically, the ratio of the integrated Compto-ionization parameter $Y_{SZ}D_A^2$ and its X-ray counterpart $Y_X$ is used as an observable to constrain the bounds on the variation of $\alpha$. Considering the violation of cosmic distance duality relation, this ratio depends on the fine structure constant as $\sim \alpha^3$. We use the quintessence model to provide the origin of $\alpha$ time variation. In order to give a robust test on $\alpha$ variation, two galaxy cluster samples, the 61 clusters provided by the Planck collaboration and the 58 clusters detected by the South Pole Telescope, are collected for analysis. Their X-ray observations are given by the XMM-Newton survey. Our results give $\zeta=-0.203^{+0.101}_{-0.099}$ for the Planck sample and $\zeta=-0.043^{+0.165}_{-0.148}$ for the SPT sample, indicating that $\alpha$ is constant with redshift within $3\sigma$ and $1\sigma$ for the two samples, respectively.

We present a method to determine sodium-abundance ratios ([Na/Fe]) using the \ion{Na}{1} D doublet lines in low-resolution ($R \sim$ 2000) stellar spectra. As stellar \ion{Na}{1} D lines are blended with those produced by the interstellar medium (ISM), we developed a technique for removing the interstellar \ion{Na}{1} D lines using the relationship between extinction, which is proportional to $E(B-V)$, and the equivalent width (EW) of the interstellar \ion{Na}{1} D absorption lines. When measuring [Na/Fe], we also considered corrections for non-local thermodynamic equilibrium (NLTE) effects. Comparisons with data from high-resolution spectroscopic surveys suggest that the expected precision of [Na/Fe] from low-resolution spectra is better than 0.3 dex for stars with [Fe/H] $>$ $-$3.0. We also present a simple application employing the estimated [Na/Fe] values for a large number of stellar spectra from the Sloan Digital Sky Survey (SDSS). After classifying the SDSS stars into Na-normal, Na-high, and Na-extreme, we explore their relation to stars in Galactic globular clusters (GCs). We find that, while the Na-high SDSS stars exhibit a similar metallicity distribution function (MDF) to that of the GCs, indicating that the majority of such stars may have originated from GC debris, the MDF of the Na-normal SDSS stars follows that of typical disk and halo stars. As there is a high fraction of carbon-enhanced metal-poor stars among the Na-extreme stars, they may have a non-GC origin, perhaps due to mass-transfer events from evolved binary companions.

Weak gravitational lensing flexions are a kind of weak lensing distortion which are defined as the spin 1 and spin 3 combinations of the third order derivatives of gravitational lensing potential. Since the shear has spin 2 combination of the second order derivative, the flexion signal gives a partly independent information from shear signal and is more sensitive to the local mass distribution than shear signal. Thus its measurement is expected to play important roles in observational cosmology. However, since the weakness of the flexion signal as well as the complicatedness of its intrinsic noise made its accurate observation very difficult. We propose a new method of measuring the flexion signal using ERA method which is a method to measure weak lensing shear without any approximation. We find two particular combinations of the flexions which provide the quantities with only lensing information and free of intrinsic noise when taken average. It is confirmed by simple numerical simulation that the statistical average of these combinations do not in fact depend on the strength of the intrinsic distortion.

High-resolution spectra are unique indicators of three-dimensional processes in exoplanetary atmospheres. For instance, in 2020, Ehrenreich et al. reported transmission spectra from the ESPRESSO spectrograph yielding an anomalously large Doppler blueshift from the ultra-hot Jupiter WASP-76b. Interpretations of these observations invoke toy model depictions of gas-phase iron condensation in lower-temperature regions of the planet's atmosphere. In this work, we forward model the atmosphere of WASP-76b with double-gray general circulation models (GCMs) and ray-striking radiative transfer to diagnose the planet's high-resolution transmission spectrum. We confirm that a physical mechanism driving strong east-west asymmetries across the terminator must exist to reproduce large Doppler blueshifts in WASP-76b's transmission spectrum. We identify low atmospheric drag and a deep radiative-convective boundary as necessary components of our GCM to produce this asymmetry (the latter is consistent with existing Spitzer phase curves). However, we cannot reproduce either the magnitude or the time-dependence of the WASP-76b Doppler signature with gas-phase iron condensation alone. Instead, we find that high-altitude, optically thick clouds composed of $\rm Al_2O_3$, Fe, or $\rm Mg_2SiO_4$ provide reasonable fits to the Ehrenreich et al. observations -- with marginal contributions from condensation. This fit is further improved by allowing a small orbital eccentricity ($e \approx 0.01$), consistent with prior WASP-76b orbital constraints. We additionally validate our forward-modeled spectra by reproducing lines of nearly all species detected in WASP-76b by Tabernero et al. 2021. Our procedure's success in diagnosing phase-resolved Doppler shifts demonstrates the benefits of physical, self-consistent, three-dimensional simulations in modeling high-resolution spectra of exoplanet atmospheres.

We re-examine the extremely metal-poor (XMP) dwarf galaxy AGC 198691 using a high quality spectrum obtained by the LBT's MODS instrument. Previous spectral observations obtained from KOSMOS on the Mayall 4-m and the Blue Channel spectrograph on the MMT 6.5-m telescope did not allow for the determination of sulfur, argon, or helium abundances. We report an updated and full chemical abundance analysis for AGC 198691, including confirmation of the extremely low "direct" oxygen abundance with a value of 12 + log(O/H) = 7.06 $\pm$ 0.03. AGC 198691's low metallicity potentially makes it a high value target for helping determine the primordial helium abundance ($Y_p$). Though complicated by a Na I night sky line partially overlaying the He I $\lambda$5876 emission line, the LBT/MODS spectrum proved adequate for determining AGC 198691's helium abundance. We employ the recently expanded and improved model of Aver et al. (2021), incorporating higher Balmer and Paschen lines, augmented by the observation of the infrared helium emission line He I $\lambda$10830 obtained by Hsyu et al. (2020). Applying our full model produced a reliable helium abundance determination, consistent with the expectation for its metallicity. Although this is the lowest metallicity object with a detailed helium abundance, unfortunately, due to its faintness (EW(H$\beta$) $<$ 100 AA) and the compromised He I $\lambda$5876, the resultant uncertainty on the helium abundance is too large to allow a significant improvement on the measurement of $Y_p$.

Apsidal motion is a gradual shift in the position of periastron. The impact of dynamic tides on apsidal motion has long been debated, because the contribution could not be quantified due to the lack of high quality observations. KIC 4544587 with tidally excited oscillations has been observed by \textit{Kepler} high-precision photometric data based on long time baseline and short-cadence schema. In this paper, we compute the rate of apsidal motion that arises from the dynamic tides as $19.05\pm 1.70$ mrad yr$^{-1}$ via tracking the orbital phase shifts of tidally excited oscillations. We also calculate the procession rate of the orbit due to the Newtonian and general relativistic contribution as $21.49 \pm 2.8$ and $2.4 \pm 0.06$ mrad yr$^{-1}$, respectively. The sum of these three factors is in excellent agreement with the total observational rate of apsidal motion $42.97 \pm 0.18$ mrad yr$^{-1}$ measured by eclipse timing variations. The tidal effect accounts for about 44\% of the overall observed apsidal motion and is comparable to that of the Newtonian term. Dynamic tides have a significant contribution to the apsidal motion. The analysis method mentioned in this paper presents an alternative approach to measuring the contribution of the dynamic tides quantitatively.

X-ray binaries in outburst typically show two canonical X-ray spectral states, i.e. hard and soft states, in which the physical properties of the accretion flow and of the jet are known to change. Recently, the JED-SAD paradigm has been proposed for black hole X-ray binaries, aimed to address the accretion-ejection interplay in these systems. According to this model, the accretion flow is composed by an outer standard Shakura-Sunyaev disk (SAD) and an inner hot Jet Emitting Disk (JED). The JED produces both the hard X-ray emission, effectively playing the role of the hot corona, and the radio jets. In this paper, we use the JED-SAD model to describe the evolution of the accretion flow in the black hole transient MAXI J1820+070 during its hard and hard-intermediate states. Contrarily to the previous applications of this model, the Compton reflection component has been taken into account. We use eight broadband X-rays spectra, including NuSTAR, NICER and Swift data, providing a total spectral coverage of 0.8-190 keV. The data were directly fitted with the JED-SAD model. Our results suggest that the optically thick disk (i.e. the SAD) does not extend down to the ISCO in any of the considered epochs. In particular, as the system evolves towards the hard/intermediate state, we find that the inner radius decreases from $\sim$60 R$_{\rm G}$ in the first observation down to $\sim$30 R$_{\rm G}$ in the last one. This trend is accompanied by an increase of the mass-accretion rate. In all hard-intermediate state observations, two reflection components, characterized by different values of ionization, are required to adequately explain the data. These components likely originate from different regions of the SAD. We show that a flared outer disk could, in principle, explain the double reflection component.

Using the DIANOGA hydrodynamical zoom-in simulation set of galaxy clusters, we analyze the dynamics traced by stars belonging to the Brightest Cluster Galaxies (BCGs) and their surrounding diffuse component, forming the intracluster light (ICL), and compare it to the dynamics traced by dark matter and galaxies identified in the simulations. We compute scaling relations between the BCG and cluster velocity dispersions and their corresponding masses (i.e. $M_\mathrm{BCG}^{\star}$- $\sigma_\mathrm{BCG}^{\star}$, $M_{200}$- $\sigma_{200}$, $M_\mathrm{BCG}^{\star}$- $M_{200}$, $\sigma_\mathrm{BCG}^{\star}$- $\sigma_{200}$), we find in general a good agreement with observational results. Our simulations also predict $\sigma_\mathrm{BCG}^{\star}$- $\sigma_{200}$ relation to not change significantly up to redshift $z=1$, in line with a relatively slow accretion of the BCG stellar mass at late times. We analyze the main features of the velocity dispersion profiles, as traced by stars, dark matter, and galaxies. As a result, we discuss that observed stellar velocity dispersion profiles in the inner cluster regions are in excellent agreement with simulations. We also report that the slopes of the BCG velocity dispersion profile from simulations agree with what is measured in observations, confirming the existence of a robust correlation between the stellar velocity dispersion slope and the cluster velocity dispersion (thus, cluster mass) when the former is computed within $0.1 R_{500}$. Our results demonstrate that simulations can correctly describe the dynamics of BCGs and their surrounding stellar envelope, as determined by the past star-formation and assembly histories of the most massive galaxies of the Universe.

The ESPRESSO spectrograph is a new powerful tool to detect and characterize extrasolar planets. Its design allows unprecedented radial velocity precision (down to a few tens of cm/s) and long-term thermo-mechanical stability. We present the first standalone detection of an extrasolar planet by blind radial velocity search using ESPRESSO and aim at showing the power of the instrument in characterizing planetary signals at different periodicities in long time spans. We use 41 ESPRESSO measurements of HD\,22496 within a time span of 895 days with a median photon noise of 18 cm/s. A radial velocity analysis is performed to test the presence of planets in the system and to account for the stellar activity of this K5-K7 main sequence star. For benchmarking and comparison, we attempt the detection with 43 archive HARPS measurements and compare the results yielded by the two datasets. We also use four TESS sectors to search for transits. We find radial velocity variations compatible with a close-in planet with an orbital period of $P=5.09071\pm0.00026$ days when simultaneously accounting for the effects of stellar activity at longer time scales ($P_{\rm rot}=34.99^{+0.58}_{-0.53}$ days). We characterize the physical and orbital properties of the planet and find a minimum mass of $5.57^{+0.73}_{-0.68}$ $\mathrm{M}_{\oplus}$, right in the dichotomic regime between rocky and gaseous planets. Although not transiting according to TESS data, if aligned with the stellar spin axis, the absolute mass of the planet must be below 16 $\mathrm{M}_{\oplus}$. We find no significant evidence for additional signals with semi-amplitudes above 56 cm/s at 95% confidence. With a modest set of radial velocity measurements, ESPRESSO is capable of detecting and characterizing low-mass planets and constrain the presence of planets in the habitable zone of K-dwarfs down to the rocky-mass regime.

We present a photometric and spectroscopic study of HD 50526, an ellipsoidal binary member of the group Double Periodic Variable stars. Performing data-mining in photometric surveys and conducting new spectroscopic observations with several spectrographs during 2008 to 2015, we obtained orbital and stellar parameters of the system. The radial velocities were analyzed with the genetic PIKAIA algorithm, whereas Doppler tomography maps for the H$\alpha$ and H$\beta$ lines were constructed with the Total Variation Minimization code. An optimized simplex-algorithm was used to solve the inverse-problem adjusting the light curve with the best stellar parameters for the system. We find an orbital period of $6.701 \pm 0.001 ~\mathrm{d}$ and a long photometric cycle of $191 \pm 2 ~\mathrm{d}$. We detected the spectral features of the coldest star, and modeled it with a $\log{g} = 2.79 \pm 0.02 ~\mathrm{dex}$ giant of mass $1.13 \pm 0.02 ~\mathrm{M_{\odot}}$ and effective temperature $10500 \pm 125 ~\mathrm{K}$. In addition, we determine a mass ratio $q= 0.206 \pm 0.033$ and that the hot star is a B-type dwarf of mass $5.48 \pm 0.02 ~\mathrm{M_{\odot}}$. The $V$-band orbital light curve can be modeled including the presence of an accretion disk around the hotter star. This fills the Roche lobe of the hotter star, and has a radius $14.74 \pm 0.02 ~\mathrm{R_{\odot}}$ and temperature at the outer edge $9400 ~\mathrm{K}$. Two bright spots located in the disk account for the global morphology of the light curve. The Doppler tomography maps of H$\alpha$ and H$\beta$, reveal complex structures of mass fluxes in the system.

We present a detailed photometric and spectroscopic analysis of DD CMa, based on published survey photometry and new spectroscopic data. We find an improved orbital period of $P_\mathrm{o}= 2.0084530 \pm 0.0000006 ~\mathrm{d}$. Our spectra reveal H$\beta$ and H$\alpha$ absorptions with weak emission shoulders and we also find color excess in the WISE multiband photometry, interpreted as signatures of circumstellar matter. We model the $V$-band orbital light curve derived from the ASAS and ASAS-SN surveys, assuming a semidetached configuration and using the mass ratio and temperature of the hotter star derived from our spectroscopic analysis. Our model indicates that the system consists of a B 2.5 dwarf and a B 9 giant of radii 3.2 and 3.7 $\mathrm{R_{\odot}}$, respectively, orbiting in a circular orbit of radius 6.75 $\mathrm{R_{\odot}}$. We also found $M_{\mathrm{c}} = 1.7 \pm 0.1 ~\mathrm{M_{\odot}}$, $T_{\mathrm{c}} = 11350 \pm 100 ~\mathrm{K}$ and $M_{\mathrm{h}} = 6.4 \pm 0.1 ~\mathrm{M_{\odot}}$, $T_{\mathrm{h}} = 20000 \pm 500 ~\mathrm{K}$, for the cooler and hotter star, respectively. We find broad single emission peaks in H$\alpha$ and H$\beta$ after subtracting the synthetic stellar spectra. Our results are consistent with mass exchange between the stars, and suggest the existence of a stream of gas being accreted onto the early B-type star.

In this paper we further elaborate on the Bose-Einstein condensate (BEC) dark matter model extended in our preceding work [Phys. Rev. D 102, 083510 (2020)] by the inclusion of 6th order (or three-particle) repulsive self-interaction term. Herein, our goal is to complete the picture through adding to the model the 4th order attractive self-interaction. The results of our analysis confirm the following: while in the preceding work the two-phase structure and the possibility of first-order phase transition was established, here we demonstrate that with the two competing self-interactions involved, the nontrivial phase structure of the enriched model remains intact. For this to hold, we study the conditions which the parameters of the model, including the interaction parameters, should satisfy. As a by-product and in order to provide some illustration, we obtain the rotation curves and the (bipartite) entanglement entropy for the case of particular dwarf galaxy.

We investigate the spatial, temporal, and spectral properties of 10 microflares from AR12721 on 2018 September 9 and 10 observed in X-rays using the Nuclear Spectroscopic Telescope ARray (NuSTAR) and the Solar Dynamic Observatory's Atmospheric Imaging Assembly and Helioseismic and Magnetic Imager (SDO/AIA and HMI). We find GOES sub-A class equivalent microflare energies of 10$^{26}$-10$^{28}$ erg reaching temperatures up to 10 MK with consistent quiescent or hot active region core plasma temperatures of 3-4 MK. One microflare (SOL2018-09-09T10:33), with an equivalent GOES class of A0.1, has non-thermal HXR emission during its impulsive phase (of non-thermal power $\sim$7$\times$10$^{24}$ erg s$^{-1}$) making it one of the faintest X-ray microflares to have direct evidence for accelerated electrons. In 4 of the 10 microflares, we find that the X-ray time profile matches fainter and more transient sources in the EUV, highlighting the need for observations sensitive to only the hottest material that reaches temperatures higher than those of the active region core ($>$5 MK). Evidence for corresponding photospheric magnetic flux cancellation/emergence present at the footpoints of 8 microflares is also observed.

We present radio observations of the most slowly rotating known radio pulsar PSR J0250+5854. With a 23.5 s period, it is close, or even beyond, the $P$-$\dot{P}$ diagram region thought to be occupied by active pulsars. The simultaneous observations with FAST, the Chilbolton and Effelsberg LOFAR international stations, and NenuFAR represent a five-fold increase in the spectral coverage of this object, with the detections at 1250 MHz (FAST) and 57 MHz (NenuFAR) being the highest- and lowest-frequency published respectively to date. We measure a flux density of $4\pm2$ $\mu$Jy at 1250 MHz and an exceptionally steep spectral index of $-3.5^{+0.2}_{-1.5}$, with a turnover below $\sim$95 MHz. In conjunction with observations of this pulsar with the GBT and the LOFAR Core, we show that the intrinsic profile width increases drastically towards higher frequencies, contrary to the predictions of conventional radius-to-frequency mapping. We examine polarimetric data from FAST and the LOFAR Core and conclude that its polar cap radio emission is produced at an absolute height of several hundreds of kilometres around 1.5 GHz, similar to other rotation-powered pulsars across the population. Its beam is significantly underfilled at lower frequencies, or it narrows because of the disappearance of conal outriders. Finally, the results for PSR J0250+5854 and other slowly spinning rotation-powered pulsars are contrasted with the radio-detected magnetars. We conclude that magnetars have intrinsically wider radio beams than the slow rotation-powered pulsars, and that consequently the latter's lower beaming fraction is what makes objects such as PSR J0250+5854 so scarce.

Context. Theoretical results showed the possibility that neutron capture elements were produced in the early Universe by two different sources: a frequent s-process source hosted by rotating massive stars, and a rare r-process source hosted most likely by neutron star mergers. The two sources produce barium with different isotopic compositions. Aims. We aim to investigate the lines of barium in two halo stars, HD 6268 and HD 4306. The spectra present an exquisite quality, both in terms of resolution (R > 100'000) and signal-to-noise (400). Due to hyperfine splitting (hfs) effects, barium lines are expected to show slightly different profiles depending on the barium isotopic fraction. Methods. We applied a standard local thermodynamic equilibrium synthesis of the barium lines. We compared the synthetic results assuming an s-process isotopic pattern or an r-process isotopic pattern for the two barium lines for each star that exhibited hfs. We also applied a methodology, less dependent on the accuracy of the theoretical Ba hfs structure, that transforms the lines of HD 4306 into those we would observe if its atmospheric parameter values (i.e. Teff, log g, micro- and macro-turbulence, Vsin i, and Ba abundance) were the same as those of HD 6268. Results. With both methods, our results show that the barium lines with hfs effects of HD 4306 are in agreement with an s-process composition and the lines in HD 6268 have a different profile, which is most likely linked to the presence of an r-process isotopic pattern. Conclusions. Two lines of barium of HD 6268 and HD 4306 seem to confirm the theoretical expectation that both r-process events and also s-process contribution by rotating massive stars have polluted the ancient halo of our Galaxy.

We search for isotropic stochastic gravitational-wave background (SGWB) in the International Pulsar Timing Array second data release. By modeling the SGWB as a power-law, we find very strong Bayesian evidence for a common-spectrum process, and further this process has scalar transverse (ST) correlations allowed in general metric theory of gravity as the Bayes factor in favor of the ST-correlated process versus the spatially uncorrelated common-spectrum process is $30\pm 2$. The median and the $90\%$ equal-tail amplitudes of ST mode are $\mathcal{A}_{\mathrm{ST}}= 1.29^{+0.51}_{-0.44} \times 10^{-15}$, or equivalently the energy density parameter per logarithm frequency is $\Omega_{\mathrm{GW}}^{\mathrm{ST}} = 2.31^{+2.19}_{-1.30} \times 10^{-9}$, at frequency of 1/year. However, we do not find any statistically significant evidence for the tensor transverse (TT) mode and then place the $95\%$ upper limits as $\mathcal{A}_{\mathrm{TT}}< 3.95 \times 10^{-15}$, or equivalently $\Omega_{\mathrm{GW}}^{\mathrm{TT}}< 2.16 \times 10^{-9}$, at frequency of 1/year.

Massive cosmological neutrinos suppress the Large-Scale Structure (LSS) in the Universe by smoothing the cosmic over-densities, and hence structure formation is delayed relative to that in the standard Lambda-Cold Dark Matter ($\Lambda$CDM) model. We characterize the merger and mass accretion history of dark matter halos with the halo formation time $a_{1/2}$, tree entropy $s$ and halo leaf function $\ell(X)$ and measure them using neutrino-involved N-body simulations. We show that a non-zero sum of neutrino masses $M_\nu$ delays the $a_{1/2}$ for halos with virial mass between $10^{13} M_\odot$ and $3\times 10^{13} M_\odot$, whereas a non-zero neutrino asymmetry parameter $\eta^2$ has the opposite effect. While the mean tree entropy $\bar s$ does not depend significantly on either $M_\nu$ or $\eta^2$, the halo leaf function does. Furthermore, the dependencies of $\ell$ on $M_\nu$ and $\eta^2$ have significant evolution in redshift $z$, with the relative contributions of $M_\nu$ and $\eta^2$ showing a sigmoid-like transition as a function of $z$ around $z \approx 0.6$. Together with the matter power spectrum, these observables allow us to break the parameter degeneracy between $M_\nu$ and $\eta^2$ so that they can both be constrained in principle.

Several supersonic runaway pulsar wind nebulae (sPWNe) with jet-like extended structures have been recently discovered in X-rays. If these structures are the product of electrons escaping the system and diffusing into the surrounding interstellar medium, they can produce a radio halo extending for several arcmin around the source. We model the expected radio emission in this scenario in the Lighthouse Nebula sPWN. We assume a constant particle injection rate during the source lifetime, and isotropic diffusion into the surrounding medium. Our predictions strongly depend on the low- and high-energy cutoffs given in the particle distribution. Our results indicate that extended radio emission can be detected from the Lighthouse Nebula without the need to invoke extreme values for the model parameters. We provide synthetic synchrotron maps that can be used to constrain these results with observations by current highly sensitive radio instruments.

Carbon-enhanced metal-poor (CEMP) r/s-stars show surface-abundance distributions characteristic of the so-called intermediate neutron capture process (i-process) of nucleosynthesis. We previously showed that the ingestion of protons in the convective helium-burning region of a low-mass low-metallicity star can explain the surface abundance distribution observed in CEMP r/s stars relatively well. Such an i-process requires detailed reaction network calculations involving hundreds of nuclei for which reaction rates have not yet been determined experimentally. We investigate the nuclear physics uncertainties affecting the i-process during the AGB phase of low-metallicity low-mass stars by propagating the theoretical uncertainties in the radiative neutron capture cross sections, as well as the 13C(a,n)16O reaction rate, and estimating their impact on the surface-abundance distribution. It is found that considering systematic uncertainties on the various nuclear ingredients affecting the radiative neutron capture rates, surface elemental abundances are typically predicted within +/-0.4 dex. The 56 < Z < 59 region of the spectroscopically relevant heavy-s elements of Ba-La-Ce-Pr as well as the r-dominated Eu element remain relatively unaffected by nuclear uncertainties. In contrast, the inclusion of the direct capture contribution impacts the rates in the neutron-rich A~45, 100, 160, and 200 regions, and the i-process production of the Z~45 and 65-70 elements. Uncertainties in the photon strength function also impact the overabundance factors by typically 0.2-0.4 dex. Nuclear level densities tend to affect abundance predictions mainly in the Z=74-79 regions. The uncertainties associated with the neutron-producing reaction 13C(a,n)16O and the unknown beta-decay rates are found to have a low impact on the overall surface enrichment

Rapid advance has been made recently in accurate distance measurements for nearby ($D < 11$ Mpc) galaxies based on the magnitude of the tip of red giant branch stars resolved with the Hubble Space Telescope. We use observational properties of galaxies presented in the last version of Updated Nearby Galaxy Catalog to derive a halo mass of luminous galaxies via orbital motion of their companions. Our sample contains 298 assumed satellites with known radial velocities around 25 Milky Way-like massive galaxies and 65 assumed satellites around 47 fainter dominant galaxies. The average total mass-to-$K$-band luminosity ratio is $31\pm6 M_\odot/L_\odot$ for the luminous galaxies, increasing up to $\sim200 M_\odot/L_\odot$ toward dwarfs. The bulge-dominated luminous galaxies are characterized with $\langle{}M_T/L_K\rangle = 73\pm15 M_\odot/L_\odot$, while the disc-dominated spirals have $\langle{}M_T/L_K\rangle = 17.4\pm2.8 M_\odot/L_\odot$. We draw attention to a particular subsample of luminous spiral galaxies with signs of declining rotation curve, which have a radial velocity dispersion of satellites less than 55 km/s and a poor dark matter halo with $\langle{}M_T/L_K\rangle = 5.5\pm1.1 M_\odot/L_\odot$. We note that a fraction of quenched (dSph, dE) companions around Milky Way-like galaxies decreases with their linear projected separation as $0.75 \exp(-R_p/350\,\mathrm{kpc})$.

Detecting continuous nanohertz gravitational waves (GWs) generated by individual close binaries of supermassive black holes (CB-SMBHs) is one of the primary objectives of pulsar timing arrays (PTAs). The detection sensitivity is slated to increase significantly as the number of well-timed millisecond pulsars will increase by more than an order of magnitude with the advent of next-generation radio telescopes. Currently, the Bayesian analysis pipeline using parallel tempering Markov chain Monte Carlo has been applied in multiple studies for CB-SMBH searches, but it may be challenged by the high dimensionality of the parameter space for future large-scale PTAs. One solution is to reduce the dimensionality by maximizing or marginalizing over uninformative parameters semi-analytically, but it is not clear whether this approach can be extended to more complex signal models without making overly simplified assumptions. Recently, the method of diffusive nested (DNest) sampling shown the capability of coping with high dimensionality and multimodality effectively in Bayesian analysis. In this paper, we apply DNest to search for continuous GWs in simulated pulsar timing residuals and find that it performs well in terms of accuracy, robustness, and efficiency for a PTA including $\mathcal{O}(10^2)$ pulsars. DNest also allows a simultaneous search of multiple sources elegantly, which demonstrates its scalability and general applicability. Our results show that it is convenient and also high beneficial to include DNest in current toolboxes of PTA analysis.

The rapid accumulation of observed Fast Radio Bursts (FRBs) originating from cosmological distances makes it likely that some will be strongly lensed by intervening matter along the line of sight. Detection of lensed FRB repeaters, which account for a noteworthy fraction of the total population, will allow not only an accurate measurement of the lensing time delay, but also follow-up high-resolution observations to pinpoint the location of the lensed images. Recent works proposed to use such strongly-lensed FRBs to derive constraints on the current expansion rate $ H_{0} $ as well as on cosmic curvature. Here we study the prospects for placing constraints on departures from general relativity via such systems. Using an ensemble of simulated events, we focus on the gravitational slip parameter $\gamma_{\rm PN}$ in screened modified gravity models and show that FRB time-delay measurements can yield constraints as tight as $ \left| \gamma_{\rm PN}-1\right| \lesssim 0.04\times(\Lambda/100\rm kpc)\times[N/10]^{-1/2} $ at $1\sigma$ with $10$ detections.

As catalogs of gravitational-wave transients grow, new records are set for the most extreme systems observed to date. The most massive observed black holes probe the physics of pair instability supernovae while providing clues about the environments in which binary black hole systems are assembled. The least massive black holes, meanwhile, allow us to investigate the purported neutron star-black hole mass gap, and binaries with unusually asymmetric mass ratios or large spins inform our understanding of binary and stellar evolution. Existing outlier tests generally implement leave-one-out analyses, but these do not account for the fact that the event being left out was by definition an extreme member of the population. This results in a bias in the evaluation of outliers. We correct for this bias by introducing a coarse-graining framework to investigate whether these extremal events are true outliers or whether they are consistent with the rest of the observed population. Our method enables us to study extremal events while testing for population model misspecification. We show that this ameliorates biases present in the leave-one-out analyses commonly used within the gravitational-wave community. Applying our method to results from the second LIGO--Virgo transient catalog, we find qualitative agreement with the conclusions of Abbott et al, ApJL 913 L7 (2021). GW190814 is an outlier because of its small secondary mass. We find that neither GW190412 nor GW190521 are outliers.

GRS 1915+105 has been in a bright flux state for more than 2 decades, but in 2018 a significant drop in flux was observed, partly due to changes in the central engine along with increased X-ray absorption. The aim of this work is to explore how X-ray spectro-polarimetry can be used to derive the basic geometrical properties of the absorbing and reflecting matter. In particular, the expected polarisation of the radiation reflected off the disc and the putative outflow is calculated. We use \textit{NuSTAR} data collected after the flux drop to derive the parameters of the system from hard X-ray spectroscopy. The spectroscopic parameters are then used to derive the expected polarimetric signal, using results from a MonteCarlo radiative transfer code both in the case of neutral and fully ionised matter. From the spectral analysis, we find that the continuum emission becomes softer with increasing flux, and that in all flux levels the obscuring matter is highly ionised. This analysis, on the other hand, confirms that spectroscopy alone is unable to put constraints on the geometry of the reflectors. Simulations show that X-ray polarimetric observations, like those that will be provided soon by the Imaging X-ray Polarimetry Explorer (IXPE), will help to determine the geometrical parameters which are left unconstrained by the spectroscopic analysis.

Returning humans to the Moon presents an unprecedented opportunity to determine the origin of volatiles stored in the permanently shaded regions (PSRs), which trace the history of lunar volcanic activity, solar wind surface chemistry, and volatile delivery to the Earth and Moon through impacts of comets, asteroids, and micrometeoroids. So far, the source of the volatiles sampled by the Lunar Crater Observation and Sensing Satellite (LCROSS) plume (Colaprete et al. 2010; Gladstone et al. 2010) has remained undetermined. We show here that the source could not be volcanic outgassing and the composition is best explained by cometary impacts. Ruling out a volcanic source means that volatiles in the top 1-3 meters of the Cabeus PSR regolith may be younger than the latest volcanic outgassing event (~1 billion years ago; Gya) (Needham et al. 2017).

For gravitational wave (GW) detected neutron star mergers, one of the leading candidates for electromagnetic (EM) counterparts is the afterglow from an ultra-relativistic jet. Where this afterglow is observed, it will likely be viewed off-axis, such as the afterglow following GW170817/GRB 170817A. The temporal behaviour of an off-axis observed GRB afterglow can be used to reveal the lateral jet structure, and statistical model fits can put constraints on the various model free-parameters. Amongst these parameters is the inclination of the system to the line of sight. Along with the GW detection, the afterglow modelling provides the best constraint on the inclination to the line-of-sight and can improve the estimates of cosmological parameters e.g. the Hubble constant, from GW-EM events. However, modelling of the afterglow depends on the assumed jet structure and, often overlooked, the effects of lateral spreading. Here we show how the inclusion of lateral spreading in the afterglow models can affect the estimated inclination of GW-EM events.

The Kennicutt-Schmidt law is an empirical relation between the star formation rate surface density ($\Sigma_{SFR}$) and the gas surface density ($\Sigma_{gas}$) in disc galaxies. The relation has a power-law form $\Sigma_{SFR} \propto \Sigma_{gas}^n$. Assuming that star formation results from gravitational collapse of the interstellar medium, $\Sigma_{SFR}$ can be determined by dividing $\Sigma_{gas}$ by the local free-fall time $t_{ff}$. The formulation of $t_{ff}$ yields the relation between $\Sigma_{SFR}$ and $\Sigma_{gas}$, assuming that a constant fraction ($\varepsilon_{SFE}$) of gas is converted into stars every $t_{ff}$. This is done here for the first time using Milgromian dynamics (MOND). Using linear stability analysis of a uniformly rotating thin disc, it is possible to determine the size of a collapsing perturbation within it. This lets us evaluate the sizes and masses of clouds (and their $t_{ff}$) as a function of $\Sigma_{gas}$ and the rotation curve. We analytically derive the relation $\Sigma_{SFR} \propto \Sigma_{gas}^{n}$ both in Newtonian and Milgromian dynamics, finding that $n=1.4$. The difference between the two cases is a change only to the constant pre-factor, resulting in increased $\Sigma_{SFR}$ of up to 25\% using MOND in the central regions of dwarf galaxies. Due to the enhanced role of disk self-gravity, star formation extends out to larger galactocentric radii than in Newtonian gravity, with the clouds being larger. In MOND, a nearly exact representation of the present-day main sequence of galaxies is obtained if $\epsilon_{SFE} = \text{constant} \approx 1.1\%$. We also show that empirically found correction terms to the Kennicutt-Schmidt law are included in the here presented relations. Furthermore, we determine that if star formation is possible, then the temperature only affects $\Sigma_{SFR}$ by at most a factor of $\sqrt{2}$.

Magnetospheric accretion models predict that matter from protoplanetary disks accretes onto the star via funnel flows which follow the stellar field lines and shock on the stellar surface leaving a hot spot with a density gradient. Previous work has inferred different densities in the hot spot, but has not been sensitive to the radial density distribution. Attempts have been made to measure this with X-ray observations, but X-ray emission only traces a fraction of the hot spot and also coronal emission. Here we report periodic ultraviolet and optical light curves of the accreting star GM Aur that display a time lag of about 1 day between their peaks. The periodicity arises as the source of the ultraviolet and optical emission moves into and out of view as it rotates along with the star. The time-lag indicates a difference in the spatial distribution of ultraviolet and optical brightness over the stellar surface. Within the framework of the magnetospheric accretion model, this indicates a radial density gradient in a hot spot on the stellar surface since different density parts of the hot spot are expected to emit radiation at different wavelengths. These results are the first observational confirmation of the magnetospheric accretion model's prediction of a density gradient in the hot spot and demonstrate the insights gained from focusing on the wavelengths where the bulk of the accretion energy can be observed.

We investigate the evolution of orbital shapes at the Inner Lindblad Resonance region of a rotating three-dimensional bar, the mass of which is growing with time. We evaluate in time-dependent models, during a 5 Gyr period, the importance of orbits with initial conditions known to play a significant role in supporting peanut-like structures in autonomous systems. These orbits are the central family of periodic orbits (x1) and vertical perturbations of it, orbits of its standard three-dimensional bifurcations at the region (x1v1 and x1v2), as well as orbits in their neighbourhood. The knowledge of the regular or chaotic character of these orbits is essential as well, because it allows us to estimate their contribution to the support of a rotating bar and, more importantly, the dynamical mechanisms that make it possible. This is calculated by means of the GALI2 index. We find that orbital patterns existing in the autonomous case, persist for longer times in the more massive bar models, and even more so in a model in which the central spheroid component of our adopted galactic potential becomes rather insignificant. The peanut-supporting orbits which we find, have a regular or, in most cases, a weakly chaotic character. There are cases in which orbits starting close to unstable periodic orbits in an autonomous model behave as regular and support the bar when its mass increases with time. An approximation for the orbital dynamics of our non-autonomous models at a certain time, can be considered that of the corresponding frozen system around that time.

The existence of cosmological compact objects with very strong gravity is a prediction of General Relativity and an exact solution of the Einstein equations. These objects are called black holes and recently we had the first observations of them. However, the theory of black hole formation has some disadvantages. In order to avoid these, some scientists suggest the existence of gravastars (gravitation vacuum stars), an alternative stellar model which seems to solve the problems of the black hole theory. In this work we compare black holes and gravastars using a wide range of the literature and we emphasize the properties of gravastars, which are consistent with the current cosmological observations. Also, we propose gravastars as the solution of the "pair-instability" effect and a possible explanation for the observed masses of the compact objects, before the collapse, from the gravitational signal GW190521, since in the formation of a gravastar there aren't mass restrictions.

In this paper, we analyze the cosmological evolution, allowed parameter space, and observational prospects for a dark sector consisting of thermally produced pseudo-Dirac fermions with a small mass splitting, coupled to the Standard Model through a dark photon. This scenario is particularly notable in the context of sub-GeV dark matter, where the mass-off-diagonal leading interaction limits applicability of both CMB energy injection constraints and standard direct detection searches. We present the first general study of the thermal history of pseudo-Dirac DM with splittings from 100 eV to MeV, focusing on the depletion of the heavier "excited" state abundance via scatterings and decays, and on the distinctive signals arising from its small surviving abundance. We analyze CMB energy injection bounds on both DM annihilation and decay, accelerator-based probes, and new line-like direct-detection signals from the excited DM down-scattering on either nuclei or electrons, as well as future search prospects in each channel. We also comment on the relevance of this signal to the few-keV Xenon1T electron excess and on possible diurnal modulation of this signal, and introduce a signal-strength parametrization to facilitate the comparison of future experimental results to theoretical expectations.

We introduce a family of solutions of Einstein's gravity minimally coupled to an anisotropic fluid, describing asymptotically flat black holes with hair and a regular horizon. These spacetimes can describe the geometry of galaxies harboring supermassive black holes, and are extensions of Einstein clusters to include horizons. They are useful to constrain the environment surrounding astrophysical black holes, using electromagnetic or gravitational-wave observations. We compute the main properties of the geometry, including the corrections to the ringdown stage induced by the external matter and fluxes by orbiting particles. The leading order effect to these corrections is a gravitational-redshift, but gravitational-wave propagation is affected by the galactic potential in a nontrivial way, and may be characterized with future observatories.

We investigate toy models for spatial and temporal instabilities in collective neutrino oscillations induced by neutrino self-interactions, with special emphasis on inhomogeneous systems with densities following a profile. Simulations are based on a mathematica program that solves the Liouville equation with or without vacuum terms, refractive terms from a background medium, and neutrino-neutrino forward scattering, in one space dimension and in time. A discrete number of momentum modes are characterized by the neutrino velocity projection on the spatial direction. We also consider the effects of charged current interaction source terms and neutral current scattering contributions. We find that refractive effects from the medium, in particular for density distributions with a profile, and neutral current non-forward scattering off the background medium can strongly influence fast collective flavor transformations. Specifically we find that if both are present, fast flavor conversions can be strongly suppressed or at least delayed.

The ultra-slow-roll (USR) inflation represents a class of single-field models with sharp deceleration of the rolling dynamics on small scales, leading to a significantly enhanced power spectrum of the curvature perturbations and primordial black hole (PBH) formation. Such a sharp transition of the inflationary background can trigger the coherent motion of scalar condensates with effective potentials governed by the rolling rate of the inflaton field. We show that a scalar condensate carrying (a combination of) baryon or lepton number can achieve successful baryogenesis through the Affleck-Dine mechanism from unconventional initial conditions excited by the USR transition. Viable parameter space for creating the correct baryon asymmetry of the Universe naturally incorporates the specific limit for PBHs to contribute significantly to dark matter, shedding light on the cosmic coincidence problem between the baryon and dark matter densities today.

It is shown that in the context of the standard model of particles, weak interaction of cosmic microwave background (CMB) and cosmic neutrino background (C$\nu$B), in the order of one loop forward scattering, can generate non-vanishing TB and EB power spectra in the presence of scalar perturbation which is in contrast with the standard scenario cosmology. By comparing our results with the current experimental data, significant information can be extracted about the nature of C$\nu$B: the contribution of CMB-C$\nu$B forward scattering for TB, TE, and EB power spectra in the case of C$\nu$B as Majorana particles is twice larger than that in the case of Dirac one. Finally, we have discussed the evidence of cosmic neutrino background in the Planck 2018 polarization data reported on a new measurement of the cosmic birefringence (CB) angle $\beta$. It is shown that the mean opacity due to cosmic neutrino background can behave as an anisotropic birefringent medium changing the linear polarization rotation angle. Despite the good agreement of our results with the reported data (for EB power spectrum), CB angle from the average value of $\beta|_\nu\simeq1/2\bar\kappa$ is about two (four) times larger than what reported in \cite{minami}, in the case of C$\nu$B as Dirac (Majorana) particles. This point will be discussed more precisely.

We investigate the nonrotating neutron stars in $f(T)$ gravity with $f(T)=T+\alpha T^2$, where $T$ is the torsion scalar in the teleparallel formalism of gravity. In particular, we utilize the SLy and BSk family of equations of state for perfect fluid to describe the neutron stellar matter and search for the effects of the $f(T)$ modification on the models of neutron stars. For positive $\alpha$, the modification results in a stronger gravitation exerted on the stellar matter, leading to a smaller stellar mass in comparison to general relativity. Moreover, there seems to be an upper limit for the central density of the neutron stars with $\alpha>0$, beyond which the effective $f(T)$ fluid would have a steplike phase transition in density and pressure profiles, collapsing the numerical system. For negative $\alpha$, the $f(T)$ modification provides additional support for neutron stars to contain larger amount of matter. We obtain the mass-radius relations of the realistic models of neutron stars and subject them to the joint constraints from the observed massive pulsars PSR J0030+0451, PSR J0740+6620, and PSR J2215+5135, and gravitational wave events GW170817 and GW190814. For BSk19 equation of state, the neutron star model in $f(T)$ gravity can accommodate all the mentioned data when $\alpha\le 3.5 G^2M_\odot^2/c^4$. For BSk20, BSk21 and SLy equations of state, the observational data constrain the model parameter $\alpha$ to be negative. If one considers the unknown compact object in the event GW190814 not to be a neutron star and hence excludes this dataset, the constraints for BSk20 and BSk21 models can be loosened to $\alpha\le 0.4 G^2M_\odot^2/c^4$ and $\alpha\le 1.9 G^2M_\odot^2/c^4$, respectively.

We present converged ab initio calculations of structure factors for elastic spin-dependent WIMP scattering off all nuclei used in dark matter direct-detection searches: $^{19}$F, $^{23}$Na, $^{27}$Al, $^{29}$Si, $^{73}$Ge, $^{127}$I, and $^{129,131}$Xe. From a set of established two- and three-nucleon interactions derived within chiral effective field theory, we construct consistent WIMP-nucleon currents at the one-body level, including effects from axial-vector two-body currents. We then apply the in-medium similarity renormalization group to construct effective valence-space Hamiltonians and consistently transformed operators of nuclear responses. Combining the recent advances of natural orbitals with three-nucleon forces expressed in large spaces, we obtain basis-space converged structure factors even in heavy nuclei. Generally results are consistent with previous calculations but in certain cases can differ by as much as 80-90\% at low momentum transfer.

We review the equation of state (EoS) models covering a large range of temperatures, baryon number densities and electron fractions presently available on the \textsc{CompOSE} database. These models are intended to be directly usable within numerical simulations of core-collapse supernovae, binary neutron star mergers and proto-neutron star evolution. We discuss their compliance with existing constraints from astrophysical observations and nuclear data. For a selection of purely nucleonic models in reasonable agreement with the above constraints, after discussing the properties of cold matter, we review thermal properties for thermodynamic conditions relevant for core-collapse supernovae and binary neutron star mergers. We find that the latter are strongly influenced by the density dependence of the nucleon effective mass. The selected bunch of models is used to investigate the EoS dependence of hot star properties, where entropy per baryon and electron fraction profiles are inspired from proto-neutron star evolution. The $\Gamma$-law analytical thermal EoS used in many simulations is found not to describe well these thermal properties of the EoS. However, it may offer a fair description of the structure of hot stars whenever thermal effects on the baryonic part are small, as shown here for proto-neutron stars starting from several seconds after bounce.

Beyond the Newtonian approximation, gravitational fields in general relativity can be described using a formalism known as gravitoelectromagnetism. In this formalism a vector potential, the gravitomagnetic potential, arises as a result of moving masses, in strong analogy with the magnetic force due to moving charges in Maxwell's theory. Gravitomagnetism can affect orbits in the gravitational field of a massive, rotating body. This raises the possibility that gravitomagnetism may serve as the dominant physics behind the anomalous rotation curves of spiral galaxies, eliminating the need for dark matter. In this essay, we methodically work out the magnitude of the gravitomagnetic equivalent of the Lorentz force and apply the result to the Milky Way. We find that the resulting contribution is too small to produce an observable effect on these orbits. We also investigate the impact of cosmological boundary conditions on the result and find that these, too, are negligible.

The formation and disappearance of spiral arms are studied by focusing on jammed Keplerian gas in a coupled map lattice (CML) with a minimal set of procedures for simulating diverse patterns in astronomical objects. The CML shows that a spiral arm is a type of traffic jam, and its motion is governed by both a gas inflow into and outflow from the jam. In particular, a new mechanism for the disappearance of spiral arms is found. It is caused not by conventional differential rotation, but by the gas flow rate difference between the light inflow and heavy outflow, here called "light-in and heavy-out", leading to the disappearance of traffic jams. Furthermore, we propose a general approximate formula for the lifetime of spiral arms, which is simply derived from the mechanism of the "light-in and heavy-out". The proposed formula is successfully applied to the CML simulations, and moreover, to the observational data of the spiral galaxy M51.