Papers for Thursday, Dec 16 2021

We use artificial intelligence (AI) to learn and infer the physics of higher order gravitational wave modes of quasi-circular, spinning, non precessing binary black hole mergers. We trained AI models using 14 million waveforms, produced with the surrogate model NRHybSur3dq8, that include modes up to $\ell \leq 4$ and $(5,5)$, except for $(4,0)$ and $(4,1)$, that describe binaries with mass-ratios $q\leq8$ and individual spins $s^z_{\{1,2\}}\in[-0.8, 0.8]$. We use our AI models to obtain deterministic and probabilistic estimates of the mass-ratio, individual spins, effective spin, and inclination angle of numerical relativity waveforms that describe such signal manifold. Our studies indicate that AI provides informative estimates for these physical parameters. This work marks the first time AI is capable of characterizing this high-dimensional signal manifold. Our AI models were trained within 3.4 hours using distributed training on 256 nodes (1,536 NVIDIA V100 GPUs) in the Summit supercomputer.

In this paper we investigate how massive central galaxies cease their star formation by comparing theoretical predictions from cosmological simulations: EAGLE, Illustris and IllustrisTNG with observations of the local Universe from the Sloan Digital Sky Survey (SDSS). Our machine learning (ML) classification reveals supermassive black hole mass ($M_{\rm BH}$) as the most predictive parameter in determining whether a galaxy is star forming or quenched at redshift $z=0$ in all three simulations. This predicted consequence of active galactic nucleus (AGN) quenching is reflected in the observations, where it is true for a range of indirect estimates of $M_{\rm BH}$ via proxies as well as its dynamical measurements. Our partial correlation analysis shows that other galactic parameters lose their strong association with quiescence, once their correlations with $M_{\rm BH}$ are accounted for. In simulations we demonstrate that it is the integrated power output of the AGN, rather than its instantaneous activity, which causes galaxies to quench. Finally, we analyse the change in molecular gas content of galaxies from star forming to passive populations. We find that both gas fractions ($f_{\rm gas}$) and star formation efficiencies (SFEs) decrease upon transition to quiescence in the observations but SFE is more predictive than $f_{\rm gas}$ in the ML passive/star-forming classification. These trends in the SDSS are most closely recovered in IllustrisTNG and are in direct contrast with the predictions made by Illustris. We conclude that a viable AGN feedback prescription can be achieved by a combination of preventative feedback and turbulence injection which together quench star formation in central galaxies.

Before the end of the epoch of reionization, the Hydrogen in the Universe was predominantly neutral. This leads to a strong attenuation of Ly$\alpha$ lines of $z\gtrsim6$ galaxies in the intergalactic medium. Nevertheless, Ly$\alpha$ has been detected up to very high redshifts ($z\sim9$) for several especially UV luminous galaxies. Here, we test to what extent the galaxy's local environment might impact the Ly$\alpha$ transmission of such sources. We present an analysis of dedicated Hubble Space Telescope (HST) imaging in the CANDELS/EGS field to search for fainter neighbours around three of the most UV luminous and most distant spectroscopically confirmed Ly$\alpha$ emitters: EGS-zs8-1, EGS-zs8-2 and EGSY-z8p7 at $z_\mathrm{spec}=7.73$, 7.48, and 8.68, respectively. We combine the multi-wavelength HST imaging with Spitzer data to reliably select $z\sim7-9$ galaxies around the central, UV-luminous sources. In all cases, we find a clear enhancement of neighbouring galaxies compared to the expected number in a blank field (by a factor $\sim 3-9\times$). Our analysis thus reveals ubiquitous overdensities around luminous Ly$\alpha$ emitting sources in the heart of the cosmic reionization epoch. We show that our results are in excellent agreement with expectations from the Dragons simulation, confirming the theoretical prediction that the first ionized bubbles preferentially formed in overdense regions. JWST follow-up observations of the neighbouring galaxies identified here will be needed to confirm their physical association and to map out the ionized regions produced by these sources.

We report the detection and characterization of a new magnetospheric star, HD 135348, based on photometric and spectropolarimetric observations. The TESS light curve of this star exhibited variations consistent with stars known to possess rigidly rotating magnetospheres (RRMs), so we obtained spectropolarimetric observations using the Robert Stobie Spectrograph (RSS) on the Southern African Large Telescope (SALT) at four different rotational phases. From these observations, we calculated the longitudinal magnetic field of the star $\langle B_z \rangle$, as well as the Alfv\'en and Kepler radii, and deduced that this star contains a centrifugal magnetosphere. However, an archival spectrum does not exhibit the characteristic "double-horned" emission profile for H$\alpha$ and the Brackett series that has been observed in many other RRM stars. This could be due to the insufficient rotational phase coverage of the available set of observations, as the spectra of these stars significantly vary with the star's rotation. Our analysis underscores the use of TESS in photometrically identifying magnetic star candidates for spectropolarimetric follow-up using ground-based instruments. We are evaluating the implementation of a machine learning classifier to search for more examples of RRM stars in TESS data.

In this work we demonstrate that the Perseus Arm is not a continuous structure of molecular gas in the second quadrant. We first show that the observed, distanced-resolved velocity structure of the Galaxy in the outer disk is capable of creating illusory spiral arms, as was first proposed by Burton (1971). Second, we measure the distances to a collection of CO clouds at velocities consistent with the Perseus arm with $135^\circ < l < 160^\circ$. We find these distances using 3D dust maps from Green et al. (2019). We determine that these molecular cloud do not preferentially lie at the distance of a purported Perseus arm, but rather extend over 3 kpc in distance, with some evidence for a closer, high pitch angle structure between 1 and 1.5 kpc away. Finally, we demonstrate that velocity perturbations of the amplitude found near the Perseus arm can wreak havoc on our interpretation of the longitude-velocity diagram for more than half of the Milky Way disk.

Observations have shown that the star-formation activity and the morphology of galaxies are closely related, but the underlying physical connection is not well understood. Using the TNG50 simulation, we explore the quenching and the morphological evolution of the 102 massive quiescent galaxies in the mass range of $10.5<\log(M_{\rm stellar}/M_{\odot})<11.5$ selected at $z=0$. The morphology of galaxies is quantified based on their kinematics, and we measure the quenching timescale of individual galaxies directly from star formation history. We show that galaxies tend to be quenched more rapidly if they: (i) are satellites in massive halos, (ii) have lower star-forming gas fractions, or (iii) inject a larger amount of black hole kinetic feedback energy. By following the global evolutionary pathways, we conclude that quiescent discs are mainly disc galaxies that are recently and slowly quenched. Approximately half of the quiescent ellipticals at $z=0$ are rapidly quenched at higher redshifts while still disc-like. While being quiescent, they gradually become more elliptical mostly by disc heating, yet these ellipticals still retain some degree of rotation. The other half of quiescent ellipticals with the most random motion-dominated kinematics build up large spheroidal components before quenching primarily by mergers, or in some cases, misaligned gas accretion. However, the mergers that contribute to morphological transformation do not immediately quench galaxies in many cases. In summary, we find that quenching and morphological transformation are decoupled. We conclude that the TNG black hole feedback -- in combination with the stochastic merger history of galaxies -- leads to a large diversity of quenching timescales and a rich morphological landscape.

We explore synergies between the Nancy Grace Roman Space Telescope and CMB lensing data to constrain dark energy and modified gravity scenarios. A simulated likelihood analysis of the galaxy clustering and weak lensing data from the Roman Space Telescope High Latitude Survey combined with CMB lensing data from the Simons Observatory is undertaken, marginalizing over important astrophysical effects and calibration uncertainties. Included in the modeling are the effects of baryons on small-scale clustering, scale-dependent growth suppression by neutrinos, as well as uncertainties in the galaxy clustering biases, in the intrinsic alignment contributions to the lensing signal, in the redshift distributions, and in the galaxy shape calibration. The addition of CMB lensing roughly doubles the dark energy figure-of-merit from Roman photometric survey data alone, varying from a factor of 1.7 to 2.4 improvement depending on the particular Roman survey configuration. Alternatively, the inclusion of CMB lensing information can compensate for uncertainties in the Roman galaxy shape calibration if it falls below the design goals. Furthermore, we report the first forecast of Roman constraints on a model-independent structure growth, parameterized by $\sigma_8 (z)$, and on the Hu-Sawicki f(R) gravity as well as an improved forecast of the phenomenological $(\Sigma_0,\mu_0)$ model. We find that CMB lensing plays a crucial role in constraining $\sigma_8(z)$ at z>2, with percent-level constraints forecasted out to z=4. CMB lensing information does not improve constraints on the f(R) models substantially. It does, however, increase the $(\Sigma_0,\mu_0)$ figure-of-merit by a factor of about 1.5.

Microlensing is a powerful technique to study the Galactic population of "dark" objects such as exoplanets both bound and unbound, brown dwarfs, low-luminosity stars, old white dwarfs, neutron stars, and almost the only way to study isolated stellar-mass black holes. The majority of previous efforts to search for gravitational microlensing events have concentrated towards high-density fields such as the Galactic bulge. Microlensing events in the Galactic plane have the advantage of closer proximity and better constrained relative proper motions, leading to better constrained lens mass estimates at the expense of a lower optical depth compared to events towards the Galactic bulge. We use the Zwicky Transient Facility (ZTF) Data Release 5 (DR5) compiled from 2018--2021 to survey the Galactic plane in the region of $|b| < 20^\circ$. We find a total of 60 candidate microlensing events including three that show a strong microlensing parallax effect. The rate of events traces Galactic structure, decreasing exponentially as a function Galactic longitude with scale length $\ell_0 \sim 37^\circ$. On average, we find Einstein timescales of our microlensing events to be about three times as long ($\sim60$ days) compared to those towards the Galactic bulge ($\sim20$ days). This pilot project demonstrates that microlensing towards the Galactic plane shows strong promise for characterization of dark objects within the Galatic disk.

Context: R Coronae Borealis (RCB) variables and their non-variable counterparts, the dustless Hydrogen-Deficient Carbon (dLHdC) stars have been known to exhibit enhanced s-processed material on their surfaces, especially Sr, Y, and Ba. No comprehensive work has been done to explore the s-process in these types of stars, however one particular RCB star, U Aqr, has been under scrutiny for its extraordinary Sr enhancement. Aims: We aim to identify RCB and dLHdC stars that have significantly enhanced Sr abundances, such as U Aqr, and use stellar evolution models to begin to estimate the type of neutron exposure that occurs in a typical HdC star. Methods: We compare the strength of the Sr II 4077 $\AA$ spectral line to Ca II H to identify the new subclass of Sr-rich HdCs. We additionally use the structural and abundance information from existing RCB MESA models to calculate the neutron exposure parameter, $\tau$ Results: We identify six stars in the Sr-rich class. Two are RCBs, and four are dLHdCs. We additionally find that the preferred RCB MESA model has a neutron exposure $\tau$ ~ 0.1 mb$^{-1}$, which is lower than the estimated $\tau$ between 0.15 and 0.6 mb$^{-1}$ for the Sr-rich star U Aqr found in the literature. We find trends in the neutron exposure corresponding to He-burning shell temperature, metallicity, and assumed s-processing site. Conclusions: We have found a sub-class of 6 HdCs known as the Sr-rich class, which tend to lie in the halo, outside the typical distribution of RCBs and dLHdCs. We find that dLHdC stars are more likely to be Sr-rich than RCBs, with an occurrence rate of ~13\% for dLHdCs and ~2\% for RCBs. This is one of the first potential spectroscopic differences between RCBs and dLHdCs, along with dLHdCs having stronger surface abundances of $^{18}$O.

Afterglow light curves of gamma-ray bursts (GRBs) exhibit very complex temporal and spectral features, such as a sudden intensity jump about one hour after the prompt emission in the optical band. We model this feature through the late collision of two relativistic shells and investigate the corresponding high-energy neutrino emission within a multi-messenger framework, while contrasting our findings with the ones from the classic fireball model. For a constant density circumburst medium, the total number of emitted neutrinos can increase by about an order of magnitude within a dynamical time when an optical jump occurs with respect to the self-similar afterglow scenario. By exploring the detection prospects with the IceCube Neutrino Observatory and future radio arrays such as IceCube-Gen2 radio, RNO-G and GRAND200k, as well as the POEMMA spacecraft, we conclude that the detection of neutrinos with IceCube-Gen2 radio could enable us to constrain the fraction of GRB afterglows with a jump as well as the properties of the circumburst medium. We also investigate the neutrino signal expected for the afterglows of GRB 100621A and a GRB 130427A-like burst with an optical jump. The detection of neutrinos from GRB afterglows could be crucial to explore the yet-to-be unveiled mechanism powering the optical jumps.

We present statistics of HeII Lya transmission spikes and large-scale absorption troughs using archival high-resolution ($R=\lambda /\Delta \lambda \simeq 12,500$-$18,000$) far-UV spectra of eight HeII-transparent quasars obtained with the Cosmic Origins Spectrograph on the Hubble Space Telescope. The sample covers the redshift range 2.5<z<3.8, thereby probing the rapidly evolving HeII absorption at the end of the HeII reionization epoch. The measured lengths of the troughs decrease dramatically from L>100cMpc at z>3 to L~30cMpc at z~2.7, signaling a significant progression of HeII reionization at these redshifts. Furthermore, unexpectedly long L~65cMpc troughs detected at z~2.9 suggest that the UV background fluctuates at larger scales than predicted by current models. By comparing the measured incidence of transmission spikes to predictions from forward-modeled mock spectra created from the outputs of a (146cMpc)^3 optically thin Nyx hydrodynamical simulation employing different UV background models, we infer the redshift evolution of the HeII photoionization rate $\Gamma_\mathrm{He\,II}(z)$. The photoionization rate decreases with increasing redshift from $\simeq 4.6\times 10^{-15}\mathrm{\,s^{-1}}$ at z~2.6 to $\simeq 1.2 \times 10^{-15}\mathrm{\,s^{-1}}$ at z~3.2, in agreement with previous inferences from the HeII effective optical depth, and following expected trends of current models of a fluctuating HeII-ionizing background.

R Coronae Borealis (RCB) and dustless Hydrogen-deficient Carbon (dLHdC) stars are believed to be remnants of low mass white dwarf mergers. These supergiant stars have peculiar hydrogen-deficient carbon-rich chemistries and stark overabundances of $^{18}$O. RCB stars undergo dust formation episodes resulting in large-amplitude photometric variations that are not seen in dLHdC stars. Recently, the sample of known dLHdC stars in the Milky Way has more than quintupled with the discovery of 27 new dLHdC stars. We present medium resolution (R$\approx3000$) near-infrared spectra of 20 newly discovered dLHdC stars. We confirm that unlike RCB stars, dLHdC stars do not show strong blueshifted ($>200$ km s$^{-1}$) He I 1.0833 $\mu$m absorption, suggesting the absence of strong, dust-driven winds around them. We also present medium resolution (R$\approx3000-8000$) $K$-band spectra for 47 RCB stars. We measure the $^{16}$O/$^{18}$O ratios of 7 dLHdC and 31 RCB stars that show $^{12}$C$^{16}$O and $^{12}$C$^{18}$O absorption bands, and present the largest sample of values of $^{16}$O/$^{18}$O for dLHdC and RCB stars to date. We find that six of the seven dLHdC stars have $^{16}$O/$^{18}$O $<0.5$, while 26 of the 31 RCB stars have $^{16}$O/$^{18}$O $>1$. We conclude that most dLHdC stars have lower $^{16}$O/$^{18}$O than most RCB stars. This confirms one of the first spectroscopic differences between RCB and dLHdC stars. Our results rule out the existing picture that RCB stars represent an evolved stage of dLHdC stars. Instead, we suggest that whether the white dwarf merger remnant is a dLHdC or RCB star depends on the mass ratios, masses and compositions of the merging white dwarfs.

Decades after their discovery, only four hydrogen-deficient carbon (HdC) stars were known to have no circumstellar dust shell. This is in complete contrast to the $\sim$130 known Galactic HdC stars that are notorious for being heavy dust producers, i.e. the R Coronae Borealis (RCB) stars. Together they form a rare class of supergiant stars that are thought to originate from the merger of CO/He white dwarf binary systems, otherwise known as the double-degenerate scenario. We searched for new dustless HdC (dLHdC) stars. We used primarily the 2MASS and GAIA eDR3 catalogues to select candidates that were followed-up spectroscopically. We discovered 27 new dLHdC stars, one new RCB star and two new EHe stars. Surprisingly, 20 of the new dLHdC stars share a characteristic of the known dLHdC star HD 148839, having lower atmospheric hydrogen deficiencies. The uncovered population of dLHdC stars exhibit a Bulge-like distribution, like the RCB stars, but show multiple differences from those that indicate they are a different population of HdC stars following its own evolutionary sequence with a fainter luminosity and also a narrow range of effective temperature. We found indication of a current low dust production activity for four of the new dLHdC stars which could be typical RCB stars passing through a transition time. We have evidence for the first time of a large range of absolute magnitudes in the overall population of HdC stars, spanning over 3 mag. In the favoured formation framework, this is explained by a large range in the initial total WD binary mass which leads to a series of evolutionary sequences with distinct maximum brightness and initial temperature. The cold Galactic RCB stars are also noticeably fainter than the Magellanic ones, possibly due to a difference in metallicity resulting in different WD mass ratio. In our Galaxy, there could be as many dLHdC stars as RCB stars.

The Cold Spot is a puzzling large-scale feature in the Cosmic Microwave Background temperature maps and its origin has been subject to active debate. As an important foreground structure at low redshift, the Eridanus supervoid was recently detected, but it was subsequently determined that, assuming the standard $\Lambda$CDM model, only about 10-20$\%$ of the observed temperature depression can be accounted for via its Integrated Sachs-Wolfe imprint. However, $R\gtrsim100~h^{-1}\mathrm{Mpc}$ supervoids elsewhere in the sky have shown ISW imprints $A_{\mathrm{ISW}}\approx5.2\pm1.6$ times stronger than expected from $\Lambda$CDM ($A_{\mathrm{ISW}}=1$), which warrants further inspection. Using the Year-3 redMaGiC catalogue of luminous red galaxies from the Dark Energy Survey, here we confirm the detection of the Eridanus supervoid as a significant under-density in the Cold Spot's direction at $z<0.2$. We also show, with $\mathrm{S/N}\gtrsim5$ significance, that the Eridanus supervoid appears as the most prominent large-scale under-density in the dark matter mass maps that we reconstructed from DES Year-3 gravitational lensing data. While we report no significant anomalies, an interesting aspect is that the amplitude of the lensing signal from the Eridanus supervoid at the Cold Spot centre is about $30\%$ lower than expected from similar peaks found in N-body simulations based on the standard $\Lambda$CDM model with parameters $\Omega_{\rm m} = 0.279$ and $\sigma_8 = 0.82$. Overall, our results confirm the causal relation between these individually rare structures in the cosmic web and in the CMB, motivating more detailed future surveys in the Cold Spot region.

The local position invariance (LPI) is one of the three major pillars of Einstein equivalence principle, ensuring the space-time independence on the outcomes of local experiments. The LPI has been tested by measuring the gravitational redshift effect in various depths of gravitational potentials. We propose a new cosmological test of the LPI by observing the asymmetry in the cross-correlation function between different types of galaxies, which predominantly arises from the gravitational redshift effect induced by the gravitational potential of halos at which the galaxies reside. We show that the ongoing and upcoming galaxy surveys can give a fruitful constraint on the LPI-violating parameter, $\alpha$, at distant universes (redshift $z\sim0.1-1.8$) over the cosmological scales (separation $s\sim5-10\, {\rm Mpc}/h$) that have not yet been explored, finding that the expected upper limit on $\alpha$ can reach $0.03$.

Since its discovery in 1995, V2400 Ophiuchi (V2400 Oph) has stood apart from most known intermediate polar cataclysmic variables due to its proposed magnetic field strength (9-27 MG) and disk-less accretion. To date, the exact accretion mechanism of the system is still unknown, and standard accretion models fail to accurately predict the peculiar behavior of its lightcurve. We present the K2 Campaign~11 light curve of V2400 Oph recording 74.19 days of photometric data cadenced at 1 minute. The light curve is dominated by aperiodic flickering and quasi-periodic oscillations, which make the beat and spin signals inconspicuous on short timescales. Notably, a log-log full power spectrum shows a break frequency at $\sim10^2$ cycles~d$^{-1}$ similar to some disk-fed systems. Through power spectral analysis, the beat and spin periods are measured as $1003.4\pm0.2$ seconds and $ 927.7\pm 0.1$ seconds respectively. A power spectrum of the entire K2 observation demonstrates beat period dominance. However, time-resolved power spectra reveals a strong dependence between observation length and the dominant frequency of the light curve. For short observations (2-12 hrs) the beat, spin, or first beat harmonic can be observed as the dominant periodic signal. Such incoherence and variability indicate a dynamical accretion system more complex than current intermediate polar theories can explain. We propose that a diamagnetic blob accretion model may serve as a plausible explanation for the accretion mechanism.

We present the first multi-wavelength simultaneous detection of QPP in a superflare (more than a thousand times stronger than known solar flares) on a cool star, in soft X-rays (SXR, with XMM-Newton) and white light (WL, with Kepler). It allowed for the first-ever analysis of oscillatory processes in a stellar flare simultaneously in thermal and non-thermal emissions, conventionally considered to come from the corona and chromosphere of the star, respectively. The observed QPP have periods $1.5 \pm 0.15$ hours (SXR) and $3 \pm 0.6$ hours (WL), and correlate well with each other. The unique relationship between the observed parameters of QPP in SXR and WL allowed us to link them with oscillations of the electric current in the flare loop, which directly affect the dynamics of non-thermal electrons and indirectly (via Ohmic heating) the thermal plasma. These findings could be considered in favour of the equivalent LCR-contour model of a flare loop, at least in the extreme conditions of a stellar superflare.

Assuming that gas and dust separate in the interstellar medium (ISM) so that high-density regions, where stars can form, are almost devoid of dust, the amount of metals being removed form the ISM can be significantly reduced (minimized astration). Here, it is shown by simple analytical models that this may increase the total metal budget of a galaxy considerably. It is suggested that these extra metals may increase the mass of dust such that the "dust budget crisis", i.e., the fact that there seem to be more dust at high redshifts than can be accounted for, can be ameliorated. Reducing the amount of astration, the metal budget can be more than doubled, in particular for systems that evolve under continuous gas accretion.

We discuss the design, construction, and operation of a new intensity interferometer, based on the campus of Southern Connecticut State University in New Haven, Connecticut. While this paper will focus on observations taken with an original two-telescope configuration, the current instrumentation consists of three portable 0.6-m Dobsonian telescopes with single-photon avalanche diode (SPAD) detectors located at the Newtonian focus of each telescope. Photons detected at each station are time-stamped and read out with timing correlators that can give cross-correlations in timing to a precision of 48 ps. We detail our observations to date with the system, which has now been successfully used at our university in 16 nights of observing. Components of the instrument were also deployed on one occasion at Lowell Observatory, where the Perkins and Hall telescopes were made to function as an intensity interferometer. We characterize the performance of the instrument in detail. In total, the observations indicate the detection of a correlation peak at the level of 6.76-sigma when observing unresolved stars, and consistency with partial or no detection when observing at a baseline sufficient to resolve the star. Using these measurements we conclude that the angular diameter of Arcturus is larger than 15 mas, and that of Vega is between 0.8 and 17 mas. While the uncertainties are large at this point, both results are consistent with measures from amplitude-based long baseline optical interferometers.

Our major aim is a height-time model $r(t)$ of the propagation of {\sl Coronal Mass Ejections (CMEs)}, where the lower corona is self-consistently connected to the heliospheric path. We accomplish this task by using the Neupert effect to derive the peak time, duration, and rate of the CME acceleration phase, as obtained from the time derivative of the {\sl soft X-ray (SXR)} light curve. This novel approach offers the advantage to obtain the kinematics of the CME height-time profile $r(t)$, the CME velocity profile $v(t)=dr(t)/dt$, and the CME acceleration profile $a(t)=dv(t)/dt$ from {\sl Geostationary Orbiting Earth Satellite (GOES)} and white-light data, without the need of {\sl hard X-ray (HXR)} data. We apply this technique to a data set of 576 (GOES X and M-class) flare events observed with GOES and the {\sl Large Angle Solar Coronagraph (LASCO)}. Our analysis yields acceleration rates in the range of $a_A = 0.1-13$ km s$^{-2}$, acceleration durations of $\tau_A = 1.2-45$ min, and acceleration distances in the range of $d_A = 3-1063$ Mm, with a median of $d_A=39$ Mm, which corresponds to the hydrostatic scale height of a corona with a temperature of $T_e \approx 0.8$ MK. The results are consistent with standard flare/CME models that predict magnetic reconnection and synchronized (primary) acceleration of CMEs in the low corona (at a height of ~0.1 R_sun), while secondary (weaker) acceleration may occur further out at heliospheric distances.

A stellar occultation by Pluto was observed on 6 June 2020 with the 1.3-m and 3.6-m telescopes located at Devasthal, Nainital, India, using imaging systems in the I and H bands, respectively. From this event, we derive a surface pressure for Pluto's atmosphere of $p_{\rm surf}= 12.23^{+0.65}_{-0.38} $~$\mu$bar. This shows that Pluto's atmosphere is in a plateau phase since mid-2015, a result which is in excellent agreement with the Pluto volatile transport model of Meza et al. (2019}. This value does not support the pressure decrease reported by independent teams, based on occultations observed in 2018 and 2019, see Young et al. (2021} and Arimatu et al. (2020), respectively.

Aims. We explore a strategy for how the Schwarzschild and mass precessions can be separated from each other despite their secular interference, by pinpointing their signatures within a single orbit. From these insights, we then seek to assess the prospects for improving the dark mass constraints in the coming years. Methods. We analysed the dependence of the osculating orbital elements and of the observables on true anomaly, and we compared these functions for models with and without extended mass. We then translated the maximum astrometric impacts within one orbit to detection thresholds given hypothetical data of different accuracies. These theoretical investigations were then supported and complemented by an extensive mock-data fitting analysis. Results. We have four main results. 1. While the mass precession almost exclusively impacts the orbit in the apocentre half, the Schwarzschild precession almost exclusively impacts it in the pericentre half, allowing for a clear separation of the effects. 2. Data that are limited to the pericentre half are not sensitive to a dark mass, while data limited to the apocentre half are, but only to a limited extent. 3. A full orbit of data is required to substantially constrain a dark mass. 4. For a full orbit of astrometric and spectroscopic data, the astrometric component in the pericentre half plays the stronger role in constraining the dark mass than the astrometric data in the apocentre half. Furthermore, we determine the 1{\sigma} dark mass detection thresholds given different datasets on one full orbit. In particular, with a full orbit of data of 50 microarcseconds (VLTI/GRAVITY) and 10 km/s (VLT/SINFONI) precision, the 1{\sigma} bound would improve to about 1000 solar masses, for example.

Broad line regions (BLRs) in high-redshift quasars provide crucial information of chemical enrichment in the early universe. Here we present a study of BLR metallicities in 33 quasars at redshift $5.7<z<6.4$. Using the near-IR spectra of the quasars obtained from the Gemini telescope, we measure their rest-frame UV emission line flux and calculate flux ratios. We then estimate BLR metallicities with empirical calibrations based on photoionization models. The inferred median metallicity of our sample is a few times the solar value, indicating that the BLR gas had been highly metal-enriched at $z\sim6$. We compare our sample with a low-redshift quasar sample with similar luminosities and find no evidence of redshift evolution in quasar BLR metallicities. This is consistent with previous studies. The Fe II$/$Mg II flux ratio, a proxy for the Fe$/\alpha$ element abundance ratio, shows no redshift evolution as well, further supporting rapid nuclear star formation at $z\sim6$. We also find that the black hole mass-BLR metallicity relation at $z\sim6$ is consistent with the relation measured at $2<z<5$, suggesting that our results are not biased by a selection effect due to this relation.

We present a sample of Ly{\alpha} emitters (LAEs) at $z\approx6.6$ from our spectroscopic survey of high-redshift galaxies using the multi-object spectrograph M2FS on the Magellan Clay telescope. The sample consists of 36 LAEs selected by the narrow-band (NB921) technique over nearly 2 deg$^2$ in the sky. These galaxies generally have high Ly{\alpha} luminosities spanning a range of ${\sim}3\times10^{42}{-}7\times10^{43}$ erg s$^{-1}$, and include some of the most Ly{\alpha}-luminous galaxies known at this redshift. They show a positive correlation between the Ly{\alpha} line width and Ly{\alpha} luminosity, similar to the relation previously found in $z\approx5.7$ LAEs. Based on the spectroscopic sample, we calculate a sophisticated sample completeness correction and derive the Ly{\alpha} luminosity function (LF) at $z\approx6.6$. We detect a density bump at the bright end of the Ly{\alpha} LF that is significantly above the best-fit Schechter function, suggesting that very luminous galaxies tend to reside in overdense regions that have formed large ionized bubbles around them. By comparing with the $z\approx5.7$ Ly{\alpha} LF, we confirm that there is a rapid LF evolution at the faint end, but a lack of evolution at the bright end. The fraction of the neutral hydrogen in the intergalactic medium at $z\approx6.6$ estimated from such a rapid evolution is about $\sim0.3\pm0.1$, supporting a rapid and rather late process of cosmic reionization.

We present a well-designed sample of more than 1000 type 1 quasars at $3.5<z<5$ and derive UV quasar luminosity functions (QLFs) in this redshift range. These quasars were selected using the Sloan Digital Sky Survey (SDSS) imaging data in SDSS Stripe 82 and overlap regions with repeat imaging observations. They are about one magnitude fainter than those found using the SDSS single-epoch data. The spectroscopic observations were conducted by the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS) as one of the BOSS ancillary programs. This quasar sample reaches $i\sim21.5$ mag and bridges previous samples from brighter surveys and deeper surveys. We use a $1/V_\mathrm{a}$ method to derive binned QLFs at $3.6<z<4.0$, $4.0<z<4.5$, and $4.5<z<4.9$, and use a double-power law model to parameterize the QLFs. We also combine our data with those in the literature to better constrain the QLFs in the context of a much wider luminosity baseline. We find that the faint-end and bright-end slopes of the QLFs in this redshift range are around $-1.7$ and $-3.7$, respectively, with uncertainties from 0.2$-$0.3 to $>0.5$. The evolution of the QLFs from $z\sim5$ to $3.5$ can be described by a pure density evolution model ($\propto10^{kz}$) and the parameter $k$ is similar to that at $5<z<7$, suggesting a nearly uniform evolution of the quasar density at $z=3.5-7$.

In Jiang et al., we detected a bright flash (hereafter GN-z11-flash) that appeared as compact continuum emission during our Keck MOSFIRE observations of the galaxy GN-z11 at z ~ 11. We performed a comprehensive analysis of the origin of the flash using all available information and our current understanding of known man-made objects or moving objects in the solar system. We found that GN-z11-flash was likely a rest-frame UV flash associated with a long gamma-ray burst (GRB) from GN-z11. Recently, Steinhardt et al., Michalowski et al., and Nir et al. reported that GN-z11-flash was more likely from a satellite. While one cannot completely rule out the possibility of unknown satellites (or debris), we find that either the chance probabilities of being a satellite estimated by these authors have been largely overestimated or their identified satellites have been ruled out in our original analysis. Our new calculations show that the probability of GN-z11-flash being a satellite is still lower than that of it being a signal originated from GN-z11.

We explore the possible advantages of extending the standard $\Lambda$CDM model by more realistic backgrounds compared to its spatially flat Robertson-Walker spacetime assumption, while preserving the underpinning physics; in particular, by simultaneously allowing non-zero spatial curvature and anisotropic expansion on top of $\Lambda$CDM, viz., the An-$o\Lambda$CDM model. This is to test whether the latest observational data still support spatial flatness and/or isotropic expansion in this case, and, if not, to explore the roles of spatial curvature and expansion anisotropy (due to its stiff fluid-like behavior) in addressing some of the current cosmological tensions associated with $\Lambda$CDM. We first present the theoretical background and explicit mathematical construction of An-$o\Lambda$CDM. Then we constrain the parameters of this model and its particular cases, namely, An-$\Lambda$CDM, $o\Lambda$CDM, and $\Lambda$CDM, by using the data sets from different observational probes, viz., Planck CMB(+Lens), BAO, SnIa Pantheon, and CC data, and discuss the results in detail. Ultimately, we conclude that (i) the observational data confirm the spatial flatness and isotropic expansion assumptions of $\Lambda$CDM, though a very small amount of expansion anisotropy cannot be excluded, e.g., $\Omega_{\sigma0}\lesssim10^{-18}$ (95% C.L.) for An-$\Lambda$CDM from CMB+Lens data, (ii) the introduction of spatial curvature or anisotropic expansion, or both, on top $\Lambda$CDM does not offer a possible relaxation to the $H_0$ tension, and (iii) the introduction of anisotropic expansion neither affects the closed space prediction from the CMB(+Lens) data nor does it improve the drastically reduced value of $H_0$ led by the closed space.

We present results of using a basic binary classification neural network model to identify likely catastrophic outlier photometric redshift estimates of individual galaxies, based only on the galaxies' measured photometric band magnitude values. We find that a simple implementation of this classification can identify a significant fraction of galaxies with catastrophic outlier photometric redshift estimates while falsely categorizing only a much smaller fraction of non-outliers. These methods have the potential to reduce the errors introduced into science analyses by catastrophic outlier photometric redshift estimates.

The optical light curves of quiescent black hole low-mass X-ray binaries often exhibit significant non-ellipsoidal variabilities, showing the photospheric radiation of the companion star is veiled by other source of optical emission. Assessing this "veiling" effect is critical to the black hole mass measurement. Here in this work, we carry out a strictly simultaneous spectroscopic and photometric campaign on the prototype of black hole low-mass X-ray binary A0620-00. We find that for each observation epoch, the extra optical flux beyond a pure ellipsoidal modulation is positively correlated with the fraction of veiling emission, indicating the accretion disk contributes most of the non-ellipsoidal variations. Meanwhile, we also obtain a K2V spectral classification of the companion, as well as the measurements of the companion's rotational velocity $v \sin i = 83.8\pm1.9$ km s$^{-1}$ and the mass ratio between the companion and the black hole $q=0.063\pm0.004$.

With the third observing run the number of gravitational wave emitting events has increased significantly. The data of the recorded events is inspected to search for overall properties on the population. The properties of subpopulations are determined and compared to predictions from simulations. It appears that the most outstanding systems follow linear relations in the parameter space of the total binary and the chirp mass. Those relations are too tight to have a stochastic origin and are supported by at least five independent events each. The origin of the correlations is still open to be confirmed, while possible sources, ranging from instrumental artefacts to unknown physics, are discussed and partly excluded. Depending on the relation's source the two events (GW190814 and GW200210_092254) having a smaller mass component between 2.5 and 2.9 Msun may reveal in a very different light.

Magnetic reconnection is a fundamental physical process converting magnetic energy into not only plasma energy but also particle energy in various astrophysical phenomena. In this letter, we show a unique dataset of a solar flare where various plasmoids were formed by a continually stretched current sheet. EUV images captured reconnection inflows, outflows, and particularly the recurring plasma blobs (plasmoids). X-ray images reveal nonthermal emission sources at the lower end of the current sheet, presumably as large plasmoids with a sufficiently amount of energetic electrons trapped in. In the radio domain, an upward slowly drifting pulsation structure, followed by a rare pair of oppositely drifting structures, was observed. These structures are supposed to map the evolution of the primary and the secondary plasmoids formed in the current sheet. Our results on plasmoids at different locations and scales shed important light on the dynamics, plasma heating, particle acceleration, and transport processes in the turbulent current sheet and provide observational evidence for the cascading magnetic reconnection process.

We use an updated version of the halo-based galaxy group catalog of Yang et al., and take surface brightness of galaxy group ($\mu_{\rm lim}$) based on projected positions and luminosities of galaxy members as a compactness proxy to divide groups into sub-systems with different compactness. By comparing various properties, including galaxy conditional luminosity function, stellar population, AGN activity, and X-ray luminosity of intra-cluster medium of carefully controlled high (HC) and low compactness (LC) group samples, we find that the group compactness plays an essential role in characterizing the detailed physical properties of group themselves and their group members, especially for low mass groups with $M_h \lesssim 10^{13.5}h^{-1}M_{\odot}$. We find that the low mass HC groups have systematically lower magnitude gap $\Delta m_{12}$ and X-ray luminosity than their LC counterparts, indicating that the high compactness groups are probably in the early stage of group merging. On the other hand, a higher fraction of passive galaxies is found in HC groups, which however is a result of systematically smaller halo-centric distance distribution of their satellite population. After controlling of both $M_h$ and halo-centric distance, we do not find any differences on both the quenching faction and AGN activity of the member galaxies between the HC and LC groups. Therefore, we conclude that the halo quenching effect, which makes the halo-centric dependence of galaxy population, is a faster process comparing to the dynamical relaxing time-scale of galaxy groups.

We present a simple analytic description of atmospheric mass loss by aerial bursts and demonstrate that mass loss from aerial bursts becomes significant when the maximum impactor size that leads to an aerial burst rather than a ground explosion, $r_o$, is larger than the minimum impactor size needed to achieve atmospheric loss, $r_{min}$. For vertical trajectories, which give the most stringent limit, this condition is approximately satisfied when $\rho_o/\rho_i \gtrsim 0.4 v_e/v_\infty$, which implies atmospheric densities need to be comparable to impactor densities for impactor velocities that are a few times the escape velocity of the planet. The range of impactor radii resulting in aerial burst-induced mass loss, $r_o-r_{min}$, increases with the ratio of the atmosphere to the impactor density and with the trajectory angle of the impactor. The range of impactor radii that result in aerial burst-induced mass loss and the atmospheric mass lost is larger in adiabatic atmospheres than isothermal atmospheres of equivalent total mass, scale height, and atmospheric surface density. Our results imply that aerial bursts are not expected to significantly contribute to the atmospheric mass-loss history of Earth, but are expected to play an important role for planets and exoplanets similar to Neptune with significant atmospheres. For Neptune-like atmospheres, the atmospheric mass ejected per impactor mass by aerial bursts is comparable to that lost by ground explosions, which implies that, for impactors following a Dohnanyi size distribution, overall loss by aerial busts is expected to exceed that by ground explosions by a factor of $(r_{ground}/r_{aerial})^{0.5}$.

We performed a detailed modelling of the background counts observed in a phoswich scintillator X-ray detector at balloon altitude, used for astronomical observations, on board small scientific balloon. We used Monte Carlo simulation technique in Geant4 simulation environment, to estimate the detector background from various plausible sources. High energy particles and radiation generated from the interaction of Galactic Cosmic Rays with the atmospheric nuclei is a major source of background counts (under normal solar condition) for such detectors. However, cosmogenic or induced radioactivity in the detector materials due to the interaction of high energy particles and natural radioactive contamination present in the detector can also contribute substantially to the detector background. We considered detailed 3D modelling of the earth's atmosphere and magnetosphere to calculate the radiation environment at the balloon altitude and deployed a proper mass model of the detector to calculate the background counts in it. The calculation satisfactorily explains the observed background in the detector at 30 km altitude (atmospheric depth: 11.5 $g/cm^{2}$) during the balloon flight experiment from a location near 14.5$^{\circ}$N geomagnetic latitude.

Abridged. Fast rotating and self-gravitating astrophysical objects suffer strong deformations from centrifugal forces. If moreover they are magnetized, they generate an electromagnetic wave that is perturbed accordingly. When stellar objects are also surrounded by an ideal plasma, a magnetosphere is formed. We study the electromagnetic configuration of a force-free magnetosphere encompassing an ideal spheroidal rotating conductor as an inner boundary. We put special emphasize to millisecond period neutron star magnetospheres, those showing a significant oblate shape. Force-free solutions are computed by numerical integration of the time-dependent Maxwell equations in spheroidal coordinates. Relevant quantities such as the magnetic field structure, the spin down luminosity, the polar cap rims and the current density are shown. We find that the force-free magnetic field produced by spheroidal stars remains very similar to their spherical counterpart. However, the spin down luminosity slightly decreases with increasing oblateness or prolateness. Moreover the polar cap area increases and always mostly encompasses the equivalent spherical star polar cap rims. The polar cap current density is also drastically affected.

With the possible spacial association to the Fermi/LAT source 3FGL J1330.0-3818, TOL 1326-379 may be the first one that is identified as a gamma-ray emitting Fanaroff-Riley type 0 radio galaxy (FR0 RG). We analyze the ~12 yr Fermi/LAT observation data of this gamma-ray source and examine its association to TOL 1326-379. We show that the gamma-ray source (named as J1331.0-3818) is tentatively detected with a TS value of 28.7, 3FGL J1330.0-3818 is out of the 95% containment of J1331.0-3818, and their positions are spatially separated ~0.2 degr. 4FGL J1331.3-3818 falls into the 68% containment of J1331.0-3818, suggesting that our result agrees with that reported in the Fourth Fermi LAT Source Catalog. TOL 1326-379 is out of the 95% containment of J1331.0-3818, and their positions are spatially separated ~0.4 degr, indicating that the association between J1331.0-3818 and TOL 1326-379 is quite ambiguous. However, we do not find other possible potential radio or X-ray counterpart within the circle centered at J1331.0-3818 with a radius of 0.4 degr. The spectral energy distribution (SED) of TOL 1326-379 shows a bimodal feature as seen in the gamma-ray emitting RGs. We fit the SED with the one-zone leptonic model and find that the average energy spectrum of J1331.0-3818 agrees with the model prediction. Assuming that J1331.0-3818 is an unidentified gamma-ray source, we derive the upper-limit of the gamma-ray flux for TOL 1326-379. It is not tight enough to exclude this possibility with the SED modeling. Based on these results, we cautiously argue that the gamma-ray source J1331.0-3818 is associated with TOL 1326-379 and its jet radiation physic is similar to those gamma-ray emitting RGs.

We present constraints on the tensor-to-scalar ratio r using a combination of BICEP/Keck 2018 and Planck PR4 data allowing us to fit for r consistently with the six parameters of the $\Lambda$CDM model without fixing any of them. In particular, we are able to derive a constraint on the reionization optical depth $\tau$ and thus propagate its uncertainty onto the posterior distribution for r. While Planck sensitivity to r is no longer comparable with ground-based measurements, combining Planck with BK18 and BAO gives results consistent with r = 0 and tightens the constraint to r < 0.032.

We address the issue of the postmaximum bump observed in the light curve of some superluminous supernovae. We rule out the popular mechanism of a circumstellar interaction suggested for the bump explanation. Instead we propose that the postmaximum bump is caused by the magnetar dipole field enhancement several months after the explosion. The modeling of SN 2019stc light curve based on the thin shell approximation implies that at the age of $\sim 90$ days the initial dipole magnetic field should be amplified by a factor of 2.8 to account for the postmaximum bump. The specific mechanism for the field amplification of the newborn magnetar on the timescale of several months has yet to be identified.

TAP, the Table Access Protocol, is a widely used Virtual Observatory specification allowing client software to interact with remote database services in a standardised way. This paper presents taplint, a tool for assessing the compliance of deployed TAP services with the the dozen or so formal specifications that form the TAP protocol stack. We provide an overview of its capabilities and operation, and the context within which it is used to improve robustness of data services.

Mini-EUSO is a detector observing the Earth in the ultraviolet band from the International Space Station through a nadir-facing window, transparent to the UV radiation, in the Russian Zvezda module. Mini-EUSO main detector consists in an optical system with two Fresnel lenses and a focal surface composed of an array of 36 Hamamatsu Multi-Anode Photo-Multiplier tubes, for a total of 2304 pixels, with single photon counting sensitivity. The telescope also contains two ancillary cameras, in the near infrared and visible ranges, to complement measurements in these bandwidths. The instrument has a field of view of 44 degrees, a spatial resolution of about 6.3 km on the Earth surface and of about 4.7 km on the ionosphere. The telescope detects UV emissions of cosmic, atmospheric and terrestrial origin on different time scales, from a few micoseconds upwards. On the fastest timescale of 2.5 microseconds, Mini-EUSO is able to observe atmospheric phenomena as Transient Luminous Events and in particular the ELVES, which take place when an electromagnetic wave generated by intra-cloud lightning interacts with the ionosphere, ionizing it and producing apparently superluminal expanding rings of several 100 km and lasting about 100 microseconds. These highly energetic fast events have been observed to be produced in conjunction also with Terrestrial Gamma-Ray Flashes and therefore a detailed study of their characteristics (speed, radius, energy...) is of crucial importance for the understanding of these phenomena. In this paper we present the observational capabilities of ELVE detection by Mini-EUSO and specifically the reconstruction and study of ELVE characteristics.

We propose a mixture model of open clusters (OCs) in the color-magnitude diagram (CMD) to measure the OC properties, including isochrone parameters (age, distance, metallicity, and dust extinction), stellar mass function, and binary parameters (binary fraction and mass ratio distribution), with unprecedented precision and reliability. The model treats an OC in the CMD as a mixture of single and binary member stars, and field stars in the same region. The cluster members are modeled based on the theoretical stellar model, mass function, and binary properties. The field component is modeled nonparametrically using a separate field star sample in the vicinity of the cluster. Unlike conventional methods that rely on stringent member selection, ours allows us to use a sample of more cluster members and attendant field stars. The larger star sample reduces the statistical error and minimizes the potential bias by retaining more stars, which is crucial for age estimation and mass function measurement. After testing with 1000 mock clusters, we measured the parameters of ten real OCs using Gaia EDR3 data. The best-fit isochrones are consistent with previous measurements in general, and provide more reasonable estimates of age for several OCs. Our method shows significantly better precision, 0.01 dex for logarithmic age and 0.01 mag for distance modulus. The inferred MF slope is -2.7 to -1.6 for clusters younger than 2 Gyr, while older clusters appear to have significantly flatter MFs. The binary fraction is 30 to 50%. The photometric and astrometric distances accord well.

We report adaptive-optics (AO) follow-up imaging of OGLE-2014-BLG-1050, which is the second binary microlensing event with space-based parallax measurements. The degeneracy in microlens parallax pi_E led to two sets of solutions, either a ~(0.9, 0.35) M_Sun binary at ~3.5 kpc, or a ~(0.2, 0.07) M_Sun binary at ~1.1 kpc. We measure the flux blended with the microlensed source by conducting Magellan AO observations, and find that the blending is consistent with the predicted lens flux from the higher-mass solution. From the combination of the AO flux measurement together with previous lensing constraints, it is estimated that} the lens system consists of a $1.05^{+0.08}_{-0.07}$ M_Sun primary and a $0.38^{+0.07}_{-0.06}$ M_Sun secondary at $3.43^{+0.19}_{-0.21}$ kpc.

We constrain the host-star flux of the microlensing planet OGLE-2014-BLG-0676Lb using adaptive optics (AO) images taken by the Magellan and Keck telescopes. We measure the flux of the light blended with the microlensed source to be K = 16.79 +/- 0.04 mag and J = 17.76 +/- 0.03 mag. Assuming that the blend is the lens star, we find that the host is a $0.73_{-0.29}^{+0.14}$ M_Sun star at a distance of $2.67_{-1.41}^{+0.77}$ kpc, where the relatively large uncertainty in angular Einstein radius measurement is the major source of uncertainty. With mass of $M_p = 3.68_{-1.44}^{+0.69}$ M_J, the planet is likely a "super Jupiter" at a projected separation of $r_{\perp} = 4.53_{-2.50}^{+1.49}$ AU, and a degenerate model yields a similar $M_p = 3.73_{-1.47}^{+0.73}$ M_J at a closer separation of $r_{\perp} = 2.56_{-1.41}^{+0.84}$ AU. Our estimates are consistent with the previous Bayesian analysis based on a Galactic model. OGLE-2014-BLG-0676Lb belongs to a sample of planets discovered in a "second-generation" planetary microlensing survey, and we attempt to systematically constrain host properties of this sample with high-resolution imaging to study the distribution of planets.

We study emission line profiles of 21 nearby low-mass ($M_*=10^4-10^7~M_\odot$) galaxies in deep medium-high resolution spectra taken with Magellan/MagE. These low-mass galaxies are actively star-forming systems with high specific star-formation rates of $\mathrm{sSFR}\sim100-1000~\mathrm{Gyr}^{-1}$ that are well above the star-formation main sequence and its extrapolation. We identify broad-line components of H$\alpha$ and [OIII]$\lambda 5007$ emission in 14 out of the 21 galaxies that cannot be explained by the MagE instrumental profile or the natural broadening of line emission. We conduct double Gaussian profile fitting to the emission of the 14 galaxies, and find that the broad-line components have line widths significantly larger than those of the narrow-line components, indicative of galactic outflows. The board-line components have moderately large line widths of $\sim 100$ km s$^{-1}$. We estimate the maximum outflow velocities $v_\mathrm{max}$ and obtain values of $\simeq 60-200$ km s$^{-1}$, which are found to be comparable to or slightly larger than the escape velocities. Positive correlations of $v_\mathrm{max}$ with star-formation rates, stellar masses, and circular velocities, extend down into this low-mass regime. Broad- to narrow-line flux ratios BNRs are generally found to be smaller than those of massive galaxies. The small $v_\mathrm{max}$ and BNRs suggest that the mass loading factors $\eta$ can be as small as 0.1 - 1 or below, in contrast to the large $\eta$ of energy-driven outflows predicted by numerical simulations.

We employ a set of high resolution N-body simulations to study the merger rate of dark matter halos. For halos of a certain mass, we define a specific merger rate by normalizing the average merger rate per halo with the logarithmic mass growth rate of the hosts at the time of accretion. Based on the simulation results, we find that this specific merger rate, $\rm{d}N_{\rm{merge}}(\xi|M,z)/\rm{d}\xi/\rm{d}\rm{lg}M(z)$, has a universal form, which is only a function of the merger mass ratio, $\xi$, and does not depend on the host halo mass, $M$, or redshift, $z$, over a wide range of masses ($10^{12}\lesssim M \lesssim10^{14}\,M_\odot/h$) and merger ratios ($\xi\ge 1e-2$). We further test with simulations of different $\Omega_m$ and $\sigma_8$, and get the same specific merger rate. The universality of the specific merger rate shows that halos in the universe are built up self-similarly, with a universal composition in the mass contributions and an absolute merger rate that grows in proportion to the halo mass growth. As a result, the absolute merger rate relates with redshift and cosmology only through the halo mass variable, whose evolution can be readily obtained from the universal mass accretion history (MAH) model of \cite{2009ApJ...707..354Z}. Lastly, we show that this universal specific merger rate immediately predicts an universal un-evolved subhalo mass function that is independent on the MAH, redshift or the final halo mass.

Gyrochronology is one of the methods currently used to estimate the age of stellar open clusters. Hundreds of new clusters, associations, and moving groups unveiled by Gaia and complemented by accurate rotation period measurements provided by recent space missions such as Kepler and TESS are allowing us to significantly improve the reliability of this method. We use gyrochronology, that is, the calibrated age-mass-rotation relation valid for low-mass stars, to measure the age of the recently discovered moving group Group X. We extracted the light curves of all candidate members from the TESS full frame images and measured their rotation periods using different period search methods. We measured the rotation period of 168 of a total of 218 stars and compared their period-colour distribution with those of two age-benchmark clusters, the Pleiades (125 Myr) and Praesepe (625 Myr), as well as with the recently characterised open cluster NGC3532 (300 Myr). As result of our analysis, we derived a gyro age of 300$\pm$60 Myr. We also applied as independent methods the fitting of the entire isochrone and of the three brightest candidate members individually with the most precise stellar parameters, deriving comparable values of 250 Myr and 290 Myr, respectively. Our dating of Group X allows us to definitively rule out the previously proposed connection with the nearby but much older Coma Berenices cluster.

Pulsar timing uses the highly stable pulsar spin period to investigate many astrophysical topics. In particular, pulsar timing arrays make use of a set of extremely well-timed pulsars and their time correlations as a challenging detector of gravitational waves. It turns out that pulsar timing arrays are particularly sensitive to ultra-low-frequency gravitational waves, which makes them complementary to other gravitational-wave detectors. Here, we summarize the basics, focusing especially on supermassive black-hole binaries and cosmic strings, which have the potential to form a stochastic gravitational-wave background in the pulsar timing array detection band, and the scientific goals on this challenging topic. We also briefly outline the recent interesting results of the main pulsar timing array collaborations, which have found strong evidence of a common-spectrum process compatible with a stochastic gravitational-wave background and mention some new perspectives that are particularly interesting in view of the forthcoming radio observatories such as the Five hundred-meter Aperture Spherical Telescope, the MeerKAT telescope, and the Square Kilometer Array.

Hot corinos are compact regions around solar-mass protostellar objects that are very rich in interstellar complex organic molecules (iCOMs). They are believed to represent the very early phases of our Solar System's birth, which was very likely also characterized by rich organic chemistry. While most of the studied hot corinos are either isolated or born in a loose protocluster, our Sun was born in a densely packed star cluster, near massive stars whose ultraviolet radiation must have contributed to shaping the evolution of the surrounding environment. In addition, internal irradiation from energetic particles ($>$10 Mev), whose imprint is seen today in the products of short-lived radionuclides in meteoritic material, is also known to have occurred during the Solar System formation. How did all these conditions affect the chemistry of the proto-Sun and its surroundings is still an open question. To answer this question, we studied HOPS-108, the hot corino located in the protosolar analogue OMC-2 FIR4. The study was carried out with ALMA at 1.3mm with an angular resolution of $\sim$100 AU. We detected 11 iCOMs such as CH$_{3}$OH HCOOCH$_{3}$ and CH$_{3}$OCH$_{3}$. Our results can be summarized as follows: (1) an enhancement of HCOOCH3 with respect to other hot corinos, (2) a [CH$_{3}$OCH$_{3}$]/[HCOOCH$_{3}$] abundance ratio of $\sim$0.2 marginally deviating from the usual trend seen in other sources ([CH$_{3}$OCH$_{3}$]/[HCOOCH$_{3}$] $\sim$1), (3) a [CH$_{2}$DOH]/[CH$_{3}$OH] ratio of 2.5\% which is lower than what is seen in Perseus and Ophiuchus hot corinos ($\sim$7\%-9\%) and similar to that seen in HH212 another source located in Orion. This might result from different physical conditions in the Orion molecular complex compared to other regions.

We have developed the 2.5D version of the basic physically motivated 1D model of Czerny & Hryniewicz (2011), i.e. Failed Radiatively Accelerated Dusty Outflow (FRADO) model. This model is based on the idea that radiation pressure acting on dust is responsible for the formation of the low ionized part of the Broad Line Region (BLR). Such radiation pressure is strong enough to form a fast outflow from the disk surface in the inner part of low ionized BLR. The outflow properties depend on the basic physical parameters, like black hole mass, Eddington ratio and gas metallicity. We here aim at estimating the disk mass loss rate due to this process, and comparing the results with outflows detected in Broad Absorption Line (BAL) quasars.

We study X-ray and soft gamma-ray spectra from the hard state of the accreting black-hole binary MAXI J1820+070. We perform analysis of joint spectra from HXMT, NuSTAR and INTEGRAL. We find an overall agreement between the spectra from all three satellites. Satisfactory fits to the data require substantial spectral complexity, with our models including two Comptonization regions and their associated disk reflection, a disk blackbody and a narrow Fe K$\alpha$ line. Our fits confirm the presence of the truncation of the reflecting optically-thick disk at least at $>$10 gravitational radii. However, we find that the HXMT data alone cannot significantly constrain the disk inner radii.

We present the results of a study aiming at retrieving the fundamental parameters of M dwarfs from spectra secured with SPIRou, the near-infrared high-resolution spectropolarimeter installed at the Canada-France-Hawaii Telescope (CFHT), in the framework of the SPIRou Legacy Survey (SLS). Our study relies on comparing observed spectra with two grids of synthetic spectra, respectively computed from PHOENIX and MARCS model atmospheres, with the ultimate goal of optimizing the precision at which fundamental parameters can be determined. In this first step, we applied our technique to 12 inactive M dwarfs with effective temperatures ($T_{\rm eff}$) ranging from 3000 to 4000 K. We implemented a benchmark to carry out a comparison of the two models used in this study. We report that the choice of model has a significant impact on the results and may lead to discrepancies in the derived parameters of 30 K in $T_{\rm eff}$ and 0.05 dex to 0.10 dex in surface gravity ($\log{g}$) and metallicity ([M/H]), as well as systematic shifts of up to 50 K in $T_{\rm eff}$ and 0.4 dex $\log{g}$ and [M/H]. The analysis is performed on high signal-to-noise ratio template SPIRou spectra, averaged over multiple observations corrected from telluric absorption features and sky lines, using both a synthetic telluric transmission model and principal component analysis. With both models, we retrieve $T_{\rm eff}$ , $\log{g}$ and [M/H] estimates in good agreement with reference literature studies, with internal error bars of about 30 K, 0.05 dex and 0.1 dex, respectively.

We report on recent improvements of the Atacama Pathfinder Experiment Control System (APECS) to cope with the ever increasing data rates and volumes. Also the very wide bandwidths of current instruments required switching to vectorized atmospheric opacity corrections using parallelization to speed these computations up for the quasi-realtime online pipeline. We look ahead at the coming years of continued APEX operations.

A sample of 22 narrow-line Seyfert 1 (NLS1) and 47 broad-line Seyfert 1 (BLS1) galaxies observed with Suzaku is used to examine the Fe K band properties of each group. Three different models are used to examine the presence of: narrow neutral Fe Ka line at 6.4 keV and ionised Fe xxv and Fe xxvi emission lines (model A); a broad emission feature at around 6 - 7 keV (model B); and an absorption edge at 7.1 keV (model C). In all three models, the neutral Fe Ka line is weaker (lower luminosity and equivalent width) in NLS1s than in BLS1s. Model (B) also finds a more significant broad component (larger equivalent width) in NLS1s than in BLS1s. The feature does not appear to be an artefact of steeper spectra in NLS1s, but rather an intrinsic property of these sources. From model (C), the optical depth of the absorption edge appears comparable between the two samples. When comparing the absorption with the emission line properties, NLS1s seem to exhibit a lower ratio of emission-to-absorption of iron than BLS1s, and than expected based on the fluorescence yield. The observed differences may arise from different torus geometries (e.g. larger opening angle in NLS1s), and/or additional sources of Fe K emission and absorption in NLS1s beyond pure fluorescence (e.g. originating in the disc and broad line region).

The presented study gives a comprehensive overview of the theory and the evidence for a systematically varying stellar initial mass function (IMF). Then we focus on the impact of this paradigm change, that is, from the universal invariant IMF to a variable IMF, on galaxy chemical evolution (GCE) studies. For this aim, we developed the first GCE code, GalIMF, that is able to incorporate the empirically calibrated environment-dependent IMF variation theory, the integrated galactic initial mass function (IGIMF) theory. In this theory, the galaxy-wide IMF is calculated by summing all the IMFs in all embedded star clusters which formed throughout the galaxy in 10 Myr time epochs. The GalIMF code recalculates the galaxy-wide IMF at each time step because the integrated galaxy-wide IMF depends on the galactic star formation rate and metallicity. The resulting galaxy-wide IMF and metal abundance evolve with time. Using this code, we examine the chemical evolution of early-type galaxies (ETGs) from dwarf to the most massive. We find that the introduction of the non-canonical IMF affects the best estimation of the galaxy properties such as their mass, star formation history, and star formation efficiency. Moreover, we are able to provide an independent estimation on the stellar formation timescale of galaxies, the type Ia supernova production efficiency, and constrain the IMF variation law of the low-mass stars. This work provides to the community the publicly available GalIMF code with improved constraints on the IMF variation and has, for the first time, resolved the discrepancy between the galaxy formation timescales obtained from stellar population synthesis and chemical enrichment studies.

The ALMA-IMF Large Program imaged a total noncontiguous area of 53pc2, covering 15 extreme, nearby protoclusters of the Milky Way. They were selected to span relevant early protocluster evolutionary stages. Our 1.3mm and 3mm observations provide continuum images that are homogeneously sensitive to point-like cores with masses of 0.2 and 0.6Msun, respectively, with a matched spatial resolution of 2000au. We also detect lines that probe the protocluster structure, kinematics, chemistry, and feedback over scales from clouds to filaments to cores. We classify ALMA-IMF protoclusters as Young, Intermediate, or Evolved based on the amount of dense gas in the cloud that has potentially been impacted by HII regions. The ALMA-IMF catalog contains 700 cores that span a mass range of 0.15-250Msun at a typical size of 2100au. We show that this core sample has no significant distance bias and can be used to build core mass functions at similar physical scales. Significant gas motions, which we highlight here in the G353.41 region, are traced down to core scales and can be used to look for inflowing gas streamers and to quantify the impact of the possible associated core mass growth on the shape of the CMF with time. Our first analysis does not reveal any significant evolution of the matter concentration from clouds to cores or from the youngest to more evolved protoclusters, indicating that cloud dynamical evolution and stellar feedback have for the moment only had a slight effect on the structure of high-density gas in our sample. Furthermore, the first-look analysis of the line richness toward bright cores indicates that the survey encompasses several tens of hot cores, of which we highlight the most massive in the G351.77 cloud. Their homogeneous characterization can be used to constrain the emerging molecular complexity in protostars of high to intermediate masses.

We present the first data release of the ALMA-IMF Large Program, which covers the 12m-array continuum calibration and imaging. The ALMA-IMF Large Program is a survey of fifteen dense molecular cloud regions spanning a range of evolutionary stages that aims to measure the core mass function (CMF). We describe the data acquisition and calibration done by the Atacama Large Millimeter/submillimeter Array (ALMA) observatory and the subsequent calibration and imaging we performed. The image products are combinations of multiple 12m array configurations created from a selection of the observed bandwidth using multi-term, multi-frequency synthesis imaging and deconvolution. The data products are self-calibrated and exhibit substantial noise improvements over the images produced from the delivered data. We compare different choices of continuum selection, calibration parameters, and image weighting parameters, demonstrating the utility and necessity of our additional processing work. Two variants of continuum selection are used and will be distributed: the ``best-sensitivity'' data, which include the full bandwidth, including bright emission lines that contaminate the continuum, and ``cleanest'', which select portions of the spectrum that are unaffected by line emission. We present a preliminary analysis of the spectral indices of the continuum data, showing that the ALMA products are able to clearly distinguish free-free emission from dust emission, and that in some cases we are able to identify optically thick emission sources. The data products are made public with this release.

Axionlike particles (ALPs) could mix with photons in the presence of astrophysical magnetic fields. Searching for this effect in gamma-ray observations of blazars has provided some of the strongest constraints on ALP parameter space so far. Previously, photon-photon dispersion of gamma-rays off of the CMB has been shown to be important for these calculations, and is universally included in ALP-photon mixing models. Here, we assess the effects of dispersion off of other photon fields within the blazar (produced by the accretion disk, the broad line region, the dust torus, starlight, and the synchrotron field) by modelling the jet and fields of the flat spectrum radio quasar 3C454.3 and propagating ALPs through the model both with and without the full dispersion calculation. We find that the full dispersion calculation can strongly affect the mixing, particularly at energies above 100 GeV -- often reducing the ALP-photon conversion probability. This could have implications for future searches planned with, e.g., the Cherenkov Telescope Array, particularly those looking for a reduced opacity of the universe at the highest energies.

During the last decades, our understanding of the universe has reached a remarkable level, being able to test cosmological predictions with an astonishing precision. Nonetheless, the nature, composition, mass and interactions of the Dark Matter still remain unknown, presenting one of the most intriguing conundrums in current cosmology. In this doctoral thesis, signatures of Dark Matter candidates which can leave an impact on the process of formation of structures and on the evolution of the Intergalactic Medium are studied. This thesis is organized in three parts. Part I is devoted to a broad introduction to the fundamentals, describing the state of the art of the topics considered. The basics of the $\Lambda$CDM are presented in Chapter 1. Chapter 2 overviews the historical progress of evidences of Dark Matter, followed by a discussion of the status and small-scale issues of the Cold Dark Matter paradigm, examining two alternative non-standard scenarios: Warm Dark Matter and Interacting Dark Matter. Chapter 3 considers Primordial Black Holes as another Dark Matter candidate, discussing the effects of accretion of surrounding matter and the enhancement of small-scale fluctuations due to the Poisson shot noise, both of which could leave an observational impact in the Intergalactic Medium. The fundamentals of the 21 cm cosmological signal are reviewed in Chapter 4, summarizing the main processes which drive the brightness temperature, and discussing its spatial fluctuations via the power spectrum. Finally, Chapter 5 is dedicated to the ionization and thermal evolution of the Intergalactic Medium during the Cosmic Dawn and the Reionization epochs. Part II includes seven original scientific articles published during the development of the PhD, which constitute the main work of this thesis. Finally, Part III contains a summary of the main results in Spanish.

Telescope Array collaboration has reported an evidence for existence of a source of ultra-high-energy cosmic ray (UHECR) events in Perseus-Pisces supercluster. We show that the mere existence of such a source imposes an upper bound on the strength of intergalactic magnetic field (IGMF) in the Taurus void lying between the Perseus-Pisces supercluster and the Milky Way galaxy. This limit is at the level of 10^{-10} G for a field with correlation length larger than the distance of the super-cluster ~ 70 Mpc. This bound is an order-of-magnitude stronger that the previously known bound on IGMF from radio Faraday rotation measurements and it is the first upper bound on magnetic field in the voids of the Large Scale Structure.

We provide a detailed map of the ages and metallicities of turnoff stars in the Milky Way disc based on data from GALAH DR3 and Gaia EDR3. From this map, we identify previously undetected features in the age-metallicity distribution of disc stars and interpret these results as indicating a three-phase formation history of the Milky Way. In the first phase, inner disc stars form along a single age-metallicity sequence and are today kinematically hot. The end of this phase is marked by a local minimum in the inner disc age distribution 10 Gyr ago. At this time, we find the stellar populations to transition from high to low alpha-element abundances and from high to low vertical velocity dispersion. In the second phase, stars form across the disc with outwardly decreasing metallicity. In this phase, inner disc stars form at super-solar metallicites in a continuation of the early age-metallicity relation while outer disc stars begin forming at metallicities at least 0.5 dex lower. Finally, the third phase is associated with a recent burst of star formation across the local disc marked by a local minimum in the age-metallicity distribution 4 to 6 Gyr ago. Future quantitative comparisons between the observed age-metallicity distribution and those of simulated galaxies could help constrain the processes driving each of the star formation phases.

We show that a small, but \textit{measurable} shift in the eclipse mid-point time of eclipsing binary (EBs) stars of $\sim$ 0.1 seconds over a decade baseline can be used to directly measure the Galactic acceleration of stars in the Milky Way at $\sim$ kpc distances from the Sun. We consider contributions to the period drift rate from dynamical mechanisms other than the Galaxy's gravitational field, and show that the Galactic acceleration can be reliably measured using a sample of $\textit{Kepler}$ EBs with orbital and stellar parameters from the literature. Given the uncertainties on the formulation of tidal decay, our approach here is necessarily approximate, and the contribution from tidal decay is an upper limit assuming the stars are not tidally synchronized. We also use simple analytic relations to search for well-timed sources in the \textit{Kepler} field, and find $\sim$ 70 additional detached EBs with low eccentricities that have estimated timing precision better than 1 second. We illustrate the method with a prototypical, precisely timed EB using an archival \textit{Kepler} light curve and a modern synthetic \textit{HST} light curve (which provides a decade baseline). This novel method establishes a realistic possibility for obtaining fundamental Galactic parameters using eclipse timing to measure Galactic accelerations, along with other emerging new methods, including pulsar timing and extreme precision radial velocity observations. This acceleration signal grows quadratically with time. Therefore, given baselines established in the near-future for distant EBs, we can expect to measure the period drift in the future with space missions like \textit{JWST} and the \textit{Roman Space Telescope}.

HD 93521 is a massive, rapidly rotating star that is located about 1 kpc above the Galactic disk, and the evolutionary age for its estimated mass is much less than the time-of-flight if it was ejected from the disk. Here we present a re-assessment of both the evolutionary and kinematical timescales for HD 93521. We calculate a time-of-flight of 39 +/- 3 Myr based upon the distance and proper motions from Gaia EDR3 and a summary of radial velocity measurements. We then determine the stellar luminosity using a rotational model combined with the observed spectral energy distribution and distance. A comparison with evolutionary tracks for rotating stars from Brott et al. yields an evolutionary age of about 5 +/- 2 Myr. We propose that the solution to the timescale discrepancy is that HD 93521 is a stellar merger product. It was probably ejected from the Galactic disk as a close binary system of lower mass stars that eventually merged to create the rapidly rotating and single massive star we observe today.

Feedback processes are expected to shape galaxy evolution by ejecting gas from galaxies and their associated dark matter haloes, and also by preventing diffuse gas from ever being accreted. We present predictions from the EAGLE simulation project for the mass budgets associated with "ejected" and "prevented" gas, as well as for ejected metals. We find that most of the baryons that are associated with haloes of mass $10^{11} < M_{200} \, /\mathrm{M_\odot} < 10^{13}$ at $z=0$ have been ejected beyond the virial radius after having been accreted. When the gas ejected from satellites (and their progenitors) is accounted for, the combined ejected mass represents half of the total baryon budget even in the most massive simulated galaxy clusters ($M_{200} \approx 10^{14.5} \, \mathrm{M_\odot}$), with the consequence that the total baryon budget exceeds the cosmic average if ejected gas is included. We find that gas is only prevented from being accreted onto haloes for $M_{200} < 10^{12} \, \mathrm{M_\odot}$, and that this component accounts for about half the total baryon budget for $M_{200} < 10^{11} \, \mathrm{M_\odot}$, with ejected gas making up most of the remaining half. For metals, most of the mass that is not locked into stars has been ejected beyond the virial radius, at least for $M_{200} < 10^{13} \, \mathrm{M_\odot}$. Finally, within the virial radius we find that most of the mass in the circum-galactic medium (CGM) has not passed through the ISM of a progenitor galaxy, for all halo masses and redshifts. About half of the CGM within half the virial radius has passed through the ISM in the past, however.

Aims: We examine the movements of mass elements within dense fibrils using passive tracer particles (corks) in order to understand fibril creation and destruction processes. Methods: Simulated fibrils were selected at times when they were visible in an H$\alpha$ image proxy. The corks were selected within fibril H$\alpha$ formation regions. From this set, a cork was selected, and the field line passing through it was constructed. Other fibrilar corks close to this fieldline were also selected. Pathlines were constructed, revealing the locations of the mass elements forward and backward in time. The forces acting on these mass elements were analysed. Results: The main process of fibrilar loading in the simulation is different to the mass loading scenario in which waves steepen into shocks and push material upwards along the fieldlines from near their footpoints. Twisted low lying fieldlines were destabilised and then they untwisted, lifting the material trapped above their apexes via the Lorentz force. Subsequently, the majority of the mass drained down the fieldlines towards one or both footpoints under gravity. Material with large horizontal velocities could also elevated in rising fieldlines, creating somewhat parabolic motions, but material was not generally moving upward along a stationary magnetic fieldline during loading. Conclusions: The processes observed in the simulation are plausible additional scenarios. Criteria for observing such events are described. It is desirable that our simulations can also form more densely-packed fibrils from material fed from the base of field footpoints. Experimental parameters required to achieve this are discussed.

Recently, microquasar jets have aroused the interest of many researchers focusing on the astrophysical plasma outflows and various jet ejections. In this work, we concentrate on the investigation of electromagnetic radiation and particle emissions from the jets of stellar black hole binary systems characterized by the hadronic content in their jets. Such emissions are reliably described within the context of the relativistic magneto-hydrodynamics. Our model calculations are based on the Fermi acceleration mechanism through which the primary particles (mainly protons) of the jet are accelerated. As a result, a small portion of thermal protons of the jet acquire relativistic energies, through shock-waves generated into the jet plasma. From the inelastic collisions of fast (non-thermal) protons with the thermal (cold) ones, secondary charged and neutral particles (pions, kaons, muons, $\eta$-particles, etc.) are created as well as electromagnetic radiation from the radio wavelength band, to X-rays and even to very high energy gamma-rays. One of our main goals is, through the appropriate solution of the transport equation and taking into account the various mechanisms that cause energy losses to the particles, to study the secondary particle concentrations within hadronic astrophysical jets. After testing our method on the Galactic MQs SS 433 and Cyg X-1, as a concrete extragalactic binary system, we examine the LMC X-1 located in the Large Magellanic Cloud, a satellite galaxy of our Milky Way Galaxy. It is worth mentioning that, for the companion O star (and its extended nebula structure) of the LMC X-1 system, new observations using spectroscopic data from VLT/UVES have been published few years ago.

There may be a chance of small-scale ephemeral liquid water formation on present day Mars, even though the current climate does not support the existence of larger bodies of water. Through a process called deliquescence, hygroscopic salts can enter solution by absorbing water vapor directly from the atmosphere. Due to the absence of in-situ deliquescence experiments so far, the most reliable way to forecast deliquescence is through atmospherical modeling, however, the locations and times when salty liquid water could emerge are not yet well known. In this paper we present our results of likely brine formation on Mars, their proposed locations and seasons, as well as the possible limiting factors. For our calculations we used the data of Laboratoire de M\'et\'eorologie Dynamique Mars General Circulation Model version 5. The results show that from L$_s$ 35$^\circ$ - L$_s$ 160$^\circ$, between 9 PM and 11 PM there is a good chance for calcium perchlorate deliquescence above 30$^\circ$ N, while in this zone the ideal regions are concentrated mostly to Acidalia Planitia and Utopia Planitia between 1 AM and 3 AM. We found that in the Southern Hemisphere, between L$_s$ 195$^\circ$ and L$_s$ 320$^\circ$, there is a noticeable, but weaker band in the vicinity of 60$^\circ$ S, and both Argyre Planitia and Hellas Planitia show some chance for brine formation. According to our statistics the key limiting factor of deliquescence could be relative humidity in most cases. Our results suggest that during summer -- early fall seasons, there could be deliquescence in both hemispheres in specific areas from the late evening until the early morning hours. There are only few studies detailing the geological and temporal distribution of brine formation through deliquescence, thus this work could be used as a good guide for future landing site analysis or in choosing a specific location for further research.

The detailed characterization of scaling laws relating the observables of cluster of galaxies to their mass is crucial for obtaining accurate cosmological constraints with clusters. In this paper, we present a comparison between the hydrostatic and lensing mass profiles of the cluster \psz\ at $z=0.59$. The hydrostatic mass profile is obtained from the combination of high resolution NIKA2 thermal Sunyaev-Zel'dovich (tSZ) and \xmm\ X-ray observations of the cluster. Instead, the lensing mass profile is obtained from an analysis of the CLASH lensing data based on the lensing convergence map. We find significant variation on the cluster mass estimate depending on the observable, the modelling of the data and the knowledge of the cluster dynamical state. This {\bf might} lead to significant systematic effects on cluster cosmological analyses for which only a single observable is generally used. From this pilot study, we conclude that the combination of high resolution SZ, X-ray and lensing data could allow us to identify and correct for these systematic effects. This would constitute a very interesting extension of the NIKA2 SZ Large Program.

We use a numerical Galactic chemical evolution model and find that the common envelope jets supernova (CEJSN) r-process scenario can account for both the very early average ratio of europium to iron and its evolution at later times in the Milky-Way (MW) Galaxy. In the CEJSN scenario a neutron star (NS) spirals-in inside a red supergiant (RSG) star all the way to the core and destroys it. According to this scenario r-process isotopes are nucleosynthesized inside neutron-rich jets that the accretion disk around the NS launches inside the core. The merger of a NS with an RSG core already takes place in the very young Galaxy. The wide variety of these mergers properties can explain the scatter in the ratio of europium to iron around the mean at low metalicites. We conclude that CEJSNe can be a major contributor to r-process nucleosynthesis.

We report on the observations of four well-localized binary black hole (BBH) mergers by the High Energy Stereoscopic System (H.E.S.S.) during the second and third observing runs of Advanced LIGO and Advanced Virgo, O2 and O3. H.E.S.S. can observe $\mathrm{20\,deg^2}$ of the sky at a time and follows up gravitational-wave (GW) events by ``tiling'' localization regions to maximize the covered localization probability. During O2 and O3, H.E.S.S. observed large portions of the localization regions, between 35\% and 75\%, for four BBH mergers (GW170814, GW190512\_180714, GW190728\_064510, and S200224ca). For these four GW events, we find no significant signal from a pointlike source in any of the observations, and set upper limits on the very high energy ($>$100 GeV) $\gamma$-ray emission. The 1-10 TeV isotropic luminosity of these GW events is below $10^{45}$ erg s$^{-1}$ at the times of the H.E.S.S. observations, around the level of the low-luminosity GRB 190829A. Assuming no changes are made to how follow-up observations are conducted, H.E.S.S. can expect to observe over 60 GW events per year in the fourth GW observing run, O4, of which eight would be observable with minimal latency.

A model that combines celestial geometry and atmospheric physics is used to calculate the dimming of artificial satellites as they enter and exit the Earth's shadow. Refraction of sunlight by the terrestrial atmosphere can illuminate a satellite while it is inside the eclipse region determined from geometry alone. Meanwhile, refraction combines with atmospheric absorption to dim the satellites for tens of km outside of that region. Spacecraft brightness is reduced more in blue light than in red because absorption of sunlight is stronger at shorter wavelengths. Observations from the MMT-9 robotic observatory are consistent with the model predictions. Tables of satellite brightness as functions of their location in the eclipse region are provided.

Soft cosmology is an extension of standard cosmology allowing for a scale-dependent equation-of-state (EoS) parameter in the dark sectors, which is one of the properties of soft materials in condensed-matter physics, that may arise either intrinsically or effectively. We use data from Supernovae Type Ia (SNIa), Baryonic Acoustic Oscillations (BAO), and Cosmic Microwave Background (CMB) probes, in order to impose observational constraints on soft dark energy and soft dark matter. We examine three simple models, corresponding to the minimum extensions of $\Lambda$CDM scenario, namely we consider that at large scales the dark sectors have the EoS's of $\Lambda$CDM model (dust dark matter and cosmological constant respectively), while at intermediate scales either dark energy or dark matter or both, may have a different EoS according to constant ``softness'' parameters $s_{de}$ and $s_{dm}$. The observational confrontation shows that for almost all datasets the softness parameters deviate from their $\Lambda$CDM values (at 2$\sigma$ confidence level for one of the models, and at $1\sigma$ for the other two), and thus the data favor soft cosmology. Finally, performing a Bayesian evidence analysis we find that the examined models are preferred over $\Lambda$CDM cosmology.

Terrestrial planets are easier to detect around M dwarfs than other types of stars, making them promising for next-generation atmospheric characterization studies. The $TESS$ mission has greatly increased the number of known M dwarf planets that we can use to perform population studies, allowing us to explore how the rocky planet occurrence rate varies with host radius, following in the footsteps of past work with $Kepler$ data. In this paper, we use simulations to assess $TESS$'s yield of small ($R_p < 2 R_\oplus$) planet candidates around nearby ($d < 30$ pc) M dwarfs. We highlight the underappreciated fact that while $TESS$ was indeed expected to find a large number of planets around M dwarfs overall, it was not expected to have a high planetary yield for the latest M dwarfs. Furthermore, we find that $TESS$ has detected significantly fewer planets around stars with $R_\star<0.3 R_\odot$ than even was expected (11 observed vs. $23\,\pm\,5$ expected). We find evidence that the photometric noise of stars in the $TESS$ bandpass increases with decreasing radius for M dwarfs. However, this trend cannot explain the observed distribution of planets. Our main conclusions are: (1) the planetary occurrence rate likely doesn't increase, and may decrease for the latest M dwarfs; and (2) there are at least 17, and potentially three times that number, transiting planets around nearby late M dwarfs that will still not be detected by the end of $TESS$'s 4th year.

As long-lived quasi-solitons from the fragmentation of a scalar condensate, oscillons may dominate the preheating era after inflation. During this period, stochastic gravitational waves can also be generated. We quantify the gravitational wave production in this period with simulations accounting for full general relativity to capture all possible non-perturbative effects. We compute the gravitational wave spectra across a range of choices of the oscillon preheating models and compare our results to a conventional perturbative approach on an FLRW background. We clarify the gauge ambiguities in computing induced gravitational waves from scenarios where dense non-perturbative objects such as oscillons are being formed. In particular, we find that the synchronous gauge tends to contain large artificial enhancements in the gravitational wave spectrum due to gauge modes if gravity plays an important role in the formation of the oscillons, while other gauge choices, such as the radiation gauge or a suitably chosen "1+log" gauge, can efficiently reduce the contributions of gauge modes. The full general relativistic simulations indicate that gravitational wave spectra obtained from the perturbative approach on the FLRW background are fairly accurate, except when oscillon formation induces strong gravitational effects, for which case there can be an order unity enhancement.

During a cosmological first-order phase transition, particles of the plasma crossing the bubble walls can radiate a gauge boson. The resulting pressure cannot be computed perturbatively for large coupling constant and/or large supercooling. We resum the real and virtual emissions at all leading-log orders, both analytically and numerically using a Monte-Carlo simulation. We find that radiated bosons are dominantly soft and that the resulting retarding pressure on relativistic bubble walls is linear both in the Lorentz boost and in the order parameter, up to a log. We further quantitatively discuss IR cut-offs, wall thickness effects, the impact of various approximations entering the calculation, and comment on the fate of radiated bosons that are reflected.

A relativistic mean field hadronic model with a dark matter (DM) particle coupled to nucleons including short-range correlations (SRC) is applied to study neutron stars (NS). The lightest neutralino is chosen as the dark particle candidate, which interacts with nucleons by the exchange of Higgs bosons. A detailed thermodynamical analysis shows that the contribution of the DM fermions to the energy density of the matter composed by these particles and nucleons is completely dominated by the DM kinetic terms. The model reproduces satisfactorily the constraints on the mass-radius diagram obtained from the analysis of the combined data from the NICER mission, LIGO collaboration, and mass measurements from radio observations. We show that the SRC balance the reduction of the neutron star mass due to the DM component, and because of that the model is able to present more massive NS. We also present a study of the effect, in the NS mass-radius profiles, of the uncertainties in some bulk parameters related to the hadronic sector. We find that it is possible to generate parametrizations, with DM content, compatible with the recent astrophysical constraints and with the uncertainty in the symmetry energy slope obtained from the results reported by the updated Lead Radius EXperiment (PREX-2).

We here present a method of performing integrals of products of spherical Bessel functions (SBFs) weighted by a power-law. Our method, which begins with double-SBF integrals, exploits a differential operator $\hat{D}$ defined via Bessel's differential equation. Application of this operator raises the power-law in steps of two. We also here display a suitable base integral expression to which this operator can be applied for both even and odd cases. We test our method by showing that it reproduces previously-known solutions. Importantly, it also goes beyond them, offering solutions in terms of singular distributions, Heaviside functions, and Gauss's hypergeometric,$\;_2{\rm F}_1$ for $all$ double-SBF integrals with positive semi-definite integer power-law weight. We then show how our method for double-SBF integrals enables evaluating $arbitrary$ triple-SBF overlap integrals, going beyond the cases currently in the literature. This in turn enables reduction of arbitrary quadruple, quintuple, and sextuple-SBF integrals and beyond into tractable forms.

We construct and simulate the dynamics of gauged vortons - circular loops of cosmic string supported by the angular momentum of trapped charge and current and provide additional details on the fully stable vorton that we have previously presented. We find that their existence and dynamical properties can be accurately predicted by an analysis based on infinite, straight superconducting strings if an additional constraint on their phase frequency is satisfied. We show a good quantitative agreement with the thin string approximation (TSA) and provide evidence that curvature corrections are inversely proportional to the vorton radius. This is verified with an energy minimisation algorithm that produces vorton solutions and subsequent axial and full three dimensional evolution codes. We find that we can predict the frequencies of each mode of oscillation, determine which modes are unstable and calculate the growth rate of the unstable modes to a high degree of accuracy.

What if Big Bang was hot from its very inception? This is possible in a bimetric theory where the source of fluctuations is thermal, requiring the model to live on a critical boundary in the space of parameters and can be realized when an anti-DBI brane moves within an $EAdS_2 \times E_3$ geometry. This setup renders the model unique, with sharp predictions for the scalar spectral index and its running. We investigate the non-Gaussian signatures of this thermal bimetric model, or ``bi-thermal'' for short. We adapt the standard calculation of non-Gaussianities for $P(X,\phi)$ models to the thermal nature of the model, emphasising how the bi-thermal peculiarities affect the calculation and alter results. This leads to precise predictions for the shape and amplitude of the three-point function of the bi-thermal model (at tree-level): $f^{\rm local} _{\rm NL} = -3/2$ and $f^{\rm equil} _{\rm NL} = -2 + 4 \sqrt{3}\pi/9 \simeq 0.4$. We also discover a new shape of flattened non-gaussianity $\propto (k_1+k_2-k_3)^{-3/2} +$ permutations, which is expected due to the excited thermal initial conditions. These results, along with our earlier predictions for the scalar power spectrum, provide sharp targets for the future generation of cosmological surveys.

Dark matter captured by interaction with electrons inside the Sun may annihilate via long-lived mediator to produce observable gamma ray signals. We utilize the Fermi Large Area Telescope data to put bounds on the dark matter electron scattering cross-section which are three to four orders of magnitude stronger than the existing limits for dark matter masses ranging between GeV to TeV scale.

It is a well known fact that, in the absence of Dark Matter, the observation of the rotation curves of galaxies cannot be explained in terms of Newtonian gravity. Rotation curves become flat in the outer regions, in contrast to what is expected according to Keplerian motion. The gravitational field far from the galactic center is weak enough to use Newtonian gravity, however even in the weak-field approximation, there are general relativistic effects without a Newtonian counterpart, such as the gravitomagnetic effects originating from mass currents. Using the gravitoelectromagnetic approach to the solution of Einstein equations in the weak-field and slow-motion approximation, we discuss some simple arguments that suggest the surprising result that gravitomagnetic effects may have an important role in understanding galactic dynamics. In addition, treating matter as a fluid of dust, we study the impact of post-Newtonian effects on the fluid vorticity.

We argue that black holes admit vortex structure. This is based both on a graviton-condensate description of a black hole as well as on a correspondence between black holes and generic objects with maximal entropy compatible with unitarity, so-called saturons. We show that due to vorticity, a $Q$-ball-type saturon of a calculable renormalisable theory obeys the same extremality bound on the spin as the black hole. Correspondingly, a black hole with extremal spin emerges as a graviton condensate with vorticity. Next, we show that in the presence of mobile charges, the global vortex traps a magnetic flux of the gauge field. This can have macroscopically-observable consequences. For instance, the most powerful jets observed in active galactic nuclei can potentially be accounted for. As a signature, such emissions can occur even without a magnetized accretion disk surrounding the black hole. The flux entrapment can provide an observational window to various hidden sectors, such as millicharged dark matter.