Papers for Wednesday, Jun 29 2022

Measuring the amount of gas and dust in protoplanetary disks is a key challenge in planet formation studies. Here we provide a new set of dust depletion factors and relative mass surface densities of gas and dust for the innermost regions of a sample of protoplanetary disks. We do this by combining stellar theory with observed refractory element abundances in both disk hosts and open cluster stars. Our results are independent of, and complementary to, those obtained from spatially resolved disk observations.

We report a Karl G. Jansky Very Large Array search for redshifted CO(1-0) emission from three HI-absorption-selected galaxies at $z \approx 2$, identified earlier in their CO(3-2) or CO(4-3) emission. We detect CO(1-0) emission from DLA B1228-113 at $z\approx2.1933$ and DLA J0918+1636 at $z\approx2.5848$; these are the first detections of CO(1-0) emission in high-$z$ HI-selected galaxies. We obtain high molecular gas masses, $\rm M_{mol}\approx10^{11}\times(\alpha_{\rm CO}/4.36)\ M_\odot$, for the two objects with CO(1-0) detections, which are a factor of $\approx1.5-2$ lower than earlier estimates. We determine the excitation of the mid$-J$ CO rotational levels relative to the $J=1$ level, r$_{ J1}$, in HI-selected galaxies for the first time, obtaining r$_{\rm 31}=1.00\pm0.20$ and r$_{\rm 41}=1.03\pm0.23$ for DLA J0918+1636, and r$_{\rm 31}=0.86\pm0.21$ for DLA B1228-113. These values are consistent with thermal excitation of the $J=3,4$ levels. The excitation of the $J=3$ level in the HI-selected galaxies is similar to that seen in massive main-sequence and sub-mm galaxies at $z\gtrsim2$, but higher than that in main-sequence galaxies at $z\approx1.5$; the higher excitation of the galaxies at $z\gtrsim2$ is likely to be due to their higher star-formation rate (SFR) surface density. We use Hubble Space Telescope Wide Field Camera 3 imaging to detect the rest-frame near-ultraviolet emission of DLA B1228-113, obtaining an NUV SFR of $4.44\pm0.47$ M$_{\odot}$ yr$^{-1}$, significantly lower than that obtained from the total infrared luminosity, indicating significant dust extinction in the $z\approx2.1933$ galaxy.

In anticipation of the new era of high-redshift exploration marked by the commissioning of the James Webb Space Telescope (JWST), we present two sets of galaxy catalogues that are designed to aid the planning and interpretation of observing programs. We provide a set of 40 wide-field lightcones with footprints spanning approximately ~ 1,000 sq. arcmin, containing galaxies up to z = 10, and a new set of 8 ultra-deep lightcones with 132 sq. arcmin footprints, containing galaxies up to z ~ 12 down to the magnitudes expected to be reached in the deepest JWST surveys. These mock lightcones are extracted from dissipationless N-body simulations and populated with galaxies using the well-established, computationally efficient Santa Cruz semi-analytic model for galaxy formation. We provide a wide range of predicted physical properties, and simulated photometry from NIRCam and many other instruments. We explore the predicted counts and luminosity functions and angular two-point correlation functions for galaxies in these simulated lightcones. We also explore the predicted field-to-field variance using multiple lightcone realizations. We find that these lightcones reproduce the available measurements of observed clustering from 0.2 < z < 7.5 very well. We provide predictions for galaxy clustering at high redshift that may be obtained from future JWST observations. All of the lightcones presented here are made available through a web-based, interactive data release portal.

We present the deepest constraints yet on the median rest-UV+optical SED of $z\sim10$ galaxies, prior to JWST science operations. We constructed stacks based on four robust $J_{125}$-dropouts, previously identified across the GOODS fields. We used archival HST/WFC3 data and the full depth Spitzer/IRAC mosaics from the GREATS program, the deepest coverage at $\sim3-5\mu$m to date. The most remarkable feature of the SED is a blue IRAC $[3.6]-[4.5]=-0.18\pm0.25$ mag color. We also find a nearly flat $H_{160}-[3.6]=0.07\pm0.22$ mag color, corresponding to a UV slope $\beta= -1.92\pm0.25$. This is consistent with previous studies, and indicative of minimal dust absorption. The observed blue IRAC color and SED fitting suggest that $z\sim10$ galaxies have very young (few $\times10$ Myr) stellar populations, with $80\%$ of stars being formed in the last $\lesssim 160$ Myr ($2\sigma$). While an exciting result, the uncertainties on the SED are too large to allow us to place strong constraints on the presence of a nebular continuum in $z\sim10$ galaxies (as might be suggested by the blue $[3.6]-[4.5] < 0$ mag color). The resulting sSFR is consistent with the specific accretion rate of dark matter halos, indicative of a star-formation efficiency showing quite limited evolution at such early epochs.

The non-linear relation between the X-ray and ultraviolet (UV) luminosity in quasars has been used to derive quasar distances and to build a Hubble diagram at redshifts up to $z\sim$ 7. This cosmological application is based on the assumption of independence of the relation on redshift and luminosity. We want to test the reliability of this hypothesis by studying the spectroscopic properties of high-redshift quasars in the X-ray and UV bands. We performed a one-by-one analysis of a sample of 130 quasars at $z>$ 2.5 with high-quality X-ray and UV spectroscopic observations. We found that not only the X-ray to UV correlation still holds at these redshifts, but its intrinsic dispersion is as low as 0.12 dex (previous works reached 0.20$-$0.22 dex). For a sample of quasars at $z\sim$ 3 with particularly high-quality observations the dispersion further drops to 0.09 dex, a value entirely accountable for by intrinsic variability and source geometry effects. The composite spectra of these quasars, in both the X-rays and the UV, do not show any difference with respect to the average spectra of quasars at lower redshifts. The absence of any spectral difference between high- and low-$z$ quasars and the tightness of the X-ray to UV relation suggests that no evolutionary effects are present in the relation. Therefore, it can be safely employed to derive quasar distances. Under this assumption, we obtain a measurement of the luminosity distance at $z\sim$ 3 with 15 % uncertainty, and in a 4$\sigma$ tension with the concordance model.

The multiscale entropy assesses the complexity of a signal across different timescales. It originates from the biomedical domain and was recently successfully used to characterize light curves as part of a supervised machine learning framework to classify stellar variability. We explore the behavior of the multiscale entropy in detail by studying its algorithmic properties in a stellar variability context and by linking it with traditional astronomical time series analysis methods. We subsequently use the multiscale entropy as the basis for an interpretable clustering framework that can distinguish hybrid pulsators with both p- and g-modes from stars with only p-mode pulsations, such as $\delta$ Sct stars, or from stars with only g-mode pulsations, such as $\gamma$ Dor stars. We find that the multiscale entropy is a powerful tool for capturing variability patterns in stellar light curves. The multiscale entropy provides insights into the pulsation structure of a star and reveals how short- and long-term variability interact with each other based on time-domain information only. We also show that the multiscale entropy is correlated to the frequency content of a stellar signal and in particular to the near-core rotation rates of g-mode pulsators. We find that our new clustering framework can successfully identify the hybrid pulsators with both p- and g-modes in sets of $\delta$ Sct and $\gamma$ Dor stars, respectively. The benefit of our clustering framework is that it is unsupervised. It therefore does not require previously labeled data and hence is not biased by previous knowledge.

Many galaxy properties are known to correlate with the environment in which the galaxies are embedded. Their cold, neutral gas content, usually assessed through 21cm HI observations, is related to many other galaxy properties as it is the underlying fuel for star formation. With its high sensitivity and broad sky coverage the blind Arecibo Galaxy Environment Survey (AGES) survey brings significant improvement to the census of HI properties of galaxies in a wide range of environments, from voids to the core of a massive cluster. Here we present an HI census over a volume of ~44000 Mpc$^{3}$ towards the merging cluster Abell 1367 and the large-scale structure (LSS) surrounding the cluster out to cz = 20000 km/s. The survey is sensitive down to a column density of N$_{HI}$ = 1.5 x 10$^{17}$ cm$^{-2}$ for emission filling the beam and a line width of 10 km/s. As an approximate mass sensitivity limit, a member of A1367 (at a distance of 92 Mpc), containing M$_{HI}$ = 2.7x10$^{8}$ M$_{\odot}$ distributed over a top-hat profile of 50 km/s width would be detected at 4$\sigma$. The results are analysed in combination with optical spectroscopy data from SDSS which we use to estimate the local galaxy density based on the Voronoi-Delaunay method. In total we detect 457 HI sources, 213 of which are detected for the first time by the AGES survey, 134 of which are presented in this article for the first time. 225 of the detections are in the cluster and 232 in the remaining volume surveyed. Here we present the full catalogue of HI detections and their basic properties, including optical ones. We concentrate on the difference between the cluster and the foreground and background LSS, revealing a continuous correlation of HI detected fraction (and HI deficiency) with local galaxy density, independent of global environment.

Correlated red noise recently reported from pulsar timing observations may be an indication of stochastic gravitational waves emitted by cosmic strings that formed during a primordial phase transition near the Grand Unification energy scale. Unfortunately, known probes of cosmic strings, namely the Cosmic Microwave Background anisotropies and string lensing of extragalactic galaxies, are not sensitive enough for low string tensions of $G\mu = 10^{-10}-10^{-7}$ that are needed to explain this putative signal. We show that strong gravitational lensing of Fast Radio Bursts (FRBs) by cosmic strings is a potentially unambiguous avenue to probe that range of string tension values. The image pair of string lensing are expected to have identical magnification factor and parity, and have a typical time delay of $\sim 10^2\,\,(G\,\mu/10^{-8})^2$ seconds. The unique spectral fingerprint of each FRB, as well as the possibility to detect correlations in the time series of the electric field of the radio waves, will enable verification of the string lensing interpretation. Very-Long-Baseline Interferometry (VLBI) observations can spatially resolve the image pair and provide a lower bound on the string tension based on the image separation. We calculate the FRB lensing rate as a function of FRB detection number for several different models of the FRB redshift distribution. We find that a survey detecting $\sim 10^5$ FRBs, in line with estimates for the detection rate of the forthcoming survey CHORD, can uncover a strong lensing event for a string tension of $G\mu \simeq 10^{-7}$. Larger FRB surveys, such as Phase 2 of the Square Kilometre Array (SKA), have the potential to significantly improve the sensitivity on the string tension to $G\mu \simeq 10^{-9}$.

We present the design and performance of the GROWTH-India telescope, a 0.7 m robotic telescope dedicated to time-domain astronomy. The telescope is equipped with a 4k back-illuminated camera giving a 0.82-degree field of view and sensitivity of m_g ~20.5 in 5-min exposures. Custom software handles observatory operations: attaining high on-sky observing efficiencies (>~ 80%) and allowing rapid response to targets of opportunity. The data processing pipelines are capable of performing PSF photometry as well as image subtraction for transient searches. We also present an overview of the GROWTH-India telescope's contributions to the studies of Gamma-ray Bursts, the electromagnetic counterparts to gravitational wave sources, supernovae, novae and solar system objects.

Over the last ten years there has been a large increase in the number of projects using sound to represent astronomical data and concepts. Motivation for these projects includes the potential to enhance scientific discovery within complex datasets, by utilising the inherent multi-dimensionality of sound and the ability of our hearing to filter signals from noise. Other motivations include creating engaging multi-sensory resources, for education and public engagement, and making astronomy more accessible to people who are blind or have low vision, promoting their participation in science and related careers. We describe potential benefits of sound within these contexts and provide an overview of the nearly 100 sound-based astronomy projects that we identified. We discuss current limitations and challenges of the approaches taken. Finally, we suggest future directions to help realise the full potential of sound-based techniques in general and to widen their application within the astronomy community.

We introduce a new set of zoom-in cosmological simulations with sub-pc resolution, intended to model extremely faint, highly magnified star-forming stellar clumps, detected at z=6.14 thanks to gravitational lensing. The simulations include feedback from individual massive stars (in both the pre-supernova and supernova phases), generated via stochastic, direct sampling of the stellar initial mass function. We adopt a modified 'delayed cooling' feedback scheme, specifically created to prevent artificial radiative loss of the energy injected by individual stars in very dense gas (n~10^3-10^5 cm^{-3}). The sites where star formation ignites are characterised by maximum densities of the order of 10^5 cm^{-3} and pressures P/k>10^7 K/cm^3, corresponding to the values of the local, turbulent regions where the densest stellar aggregates form. The total stellar mass at z=6.14 is 3.4x10^7 M_sun, in satisfactory agreement with the observed stellar mass of the observed systems. The most massive clumps have masses of ~10^6 M_sun and half-mass sizes of ~100 pc. These sizes are larger than the observed ones, including also other samples of lensed high-redshift clumps, and imply an average density one order of magnitude lower than the observed one. In the size-mass plane, our clumps populate a sequence that is intermediate between the ones of observed high-redshift clumps and local dSph galaxies.

A multi-disciplinary team recently came together online to discuss the application of sonification in astronomy, focussing on the effective use of sound for scientific discovery and for improving accessibility to astronomy research and education. Here we provide a meeting report.

We present a catalog of ~100,000 periodic variable stars in TESS FFI data among members of widely distributed moving groups identified with Gaia in the previous papers in the series. By combining the periods from our catalog attributable to rotation with previously derived rotation periods for benchmark open clusters, we develop an empirical gyrochronology model of angular momentum evolution that is valid for stars with ages 10-1000 Myr. Excluding stars rotating faster than 2 days, which we find are predominantly binaries, we achieve a typical age precision of ~0.2-0.3 dex and improving at older ages. Importantly, these empirical relations apply to not only FGK-type stars but also M-type stars, due to the angular momentum distribution being much smoother, simpler, continuous and monotonic as compared to the rotation period distribution. As a result, we are also able to begin tracing in fine detail the nature of angular momentum loss in low-mass stars as functions of mass and age. We characterize the stellar variability amplitudes of the cool stars as functions of mass and age, which may correlate with the starspot covering fractions. We also identify pulsating variables among the hotter stars in the catalog, including $\delta$ Scuti, $\gamma$ Dor and SPB-type variables. These data represent an important step forward in being able to estimate precise ages of FGK- and M-type stars in the field, starting as early as the pre-main-sequence phase of evolution.

We study the initiation of thermonuclear detonations in tidally disrupted white dwarf stars by intermediate-mass ($10^3 M_\odot$) black holes. The length scales required to resolve the initiation mechanism are not easily reached in three-dimensions, so instead we have devised two-dimensional proxy models which, together with a logarithmic gridding strategy, can adequately capture detonation wave fronts as material undergoes simultaneous compression and stretching from tidal forces. We consider 0.15 and 0.6 solar mass white dwarf stars parameterized by tidal strengths in the range $\beta=4~\text{to}~23$. High spatial resolution elucidates the manner and conditions leading to thermonuclear detonation, linking the initiation sequence to stellar composition and tidal strength. All of our models suffer sustained detonations triggered by a combination of adiabatic compression, mild thermonuclear preconditioning, and collisional heating, in degrees depending primarily on tidal strength. We find many diagnostics, such as temperature, total released energy, and iron group products, are fairly well-converged (better than 10%) at resolutions below 10 km along the scale height of the orbital plane. The exceptions are intermediate mass transients like calcium, which remain uncertain up to factors of two even at 1 km resolution.

We present high resolution Karl G. Jansky Very Large Array (VLA) observations of the protostar L1527 IRS at 7 mm, 1.3 cm, and 2 cm wavelengths. We detect the edge-on dust disk at all three wavelengths and find that it is asymmetric, with the southern side of the disk brighter than the northern side. We confirm this asymmetry through analytic modeling and also find that the disk is flared at 7 mm. We test the data against models including gap features in the intensity profile, and though we cannot rule such models out, they do not provide a statistically significant improvement in the quality of fit to the data. From these fits, we can however place constraints on allowed properties of any gaps that could be present in the true, underlying intensity profile. The physical nature of the asymmetry is difficult to associate with physical features due to the edge-on nature of the disk, but could be related to spiral arms or asymmetries seen in other imaging of more face-on disks.

Clusters of galaxies are excellent laboratories for studying recurring nuclear activity in galactic nuclei since their hot gaseous medium can vastly prolong the detectability of their radio lobes via better confinement. We report here a multi-band study of the sparsely studied galaxy cluster Abell 980, based on our analysis of {\it{Chandra}} X-ray and the GMRT (150 and 325 MHz) and EVLA (1.5 GHz) radio archival data, revealing an unusually rich phenomenology. It is shown to be a quasi-relaxed cluster with a cool core ($T\sim4.2$ keV) surrounded by a hot and extensive intracluster medium (ICM) at $T\sim6.8$ keV. The radio emission shows a rich diversity, having (i) two large diffuse sources of ultra-steep spectrum (USS), extending to opposite extremities of the ICM, each associated with an X-ray brightness discontinuity (cold front), (ii) a bright radio-double of size $\sim55$ kpc coinciding with the central BCG, and (iii) a diffuse radio source, likely a mini-halo of size $\sim110$~kpc around the BCG which possesses a huge ellipsoidal stellar halo of extent $\sim 80$~kpc. The association of cold fronts with two highly aged (~ 260 Myr) USS sources in a cool-core cluster, makes it a very rare system. These USS sources are probably radio lobes from a previous episode of jet activity in the BCG, driven buoyantly towards the outskirts of the X-ray halo, thereby creating the cold fronts. A deeper radio image of this cluster may provide a rare opportunity to verify the recently proposed alternative model which explains radio mini-haloes as the aggregate radio emission from Type I supernovae occurring in the giant stellar halo extended across the cluster core.

In a companion paper, we presented BayesCal, a mathematical formalism for mitigating sky-model incompleteness in interferometric calibration. In this paper, we demonstrate the use of BayesCal to calibrate the degenerate gain parameters of full-Stokes simulated observations with a HERA-like hexagonal close-packed redundant array, for three assumed levels of completeness of the a priori known component of the calibration sky model. We compare the BayesCal calibration solutions to those recovered by calibrating the degenerate gain parameters with only the a priori known component of the calibration sky model both with and without imposing physically motivated priors on the gain amplitude solutions and for two choices of baseline length range over which to calibrate. We find that BayesCal provides calibration solutions with up to four orders of magnitude lower power in spurious gain amplitude fluctuations than the calibration solutions derived for the same data set with the alternate approaches, and between $\sim10^7$ and $\sim10^{10}$ times smaller than in the mean degenerate gain amplitude on the full range of spectral scales accessible in the data. Additionally, we find that in the scenarios modelled only BayesCal has sufficiently high fidelity calibration solutions for unbiased recovery of the 21 cm power spectrum on large spectral scales ($k_\parallel \lesssim 0.15~h\mathrm{Mpc}^{-1}$). In all other cases, in the completeness regimes studied, those scales are contaminated.

High fidelity radio interferometric data calibration that minimises spurious spectral structure in the calibrated data is essential in astrophysical applications, such as 21 cm cosmology, which rely on knowledge of the relative spectral smoothness of distinct astrophysical emission components to extract the signal of interest. Existing approaches to radio interferometric calibration have been shown to impart spurious spectral structure to the calibrated data if the sky model used to calibrate the data is incomplete. In this paper, we introduce BayesCal: a novel solution to the sky-model incompleteness problem in interferometric calibration, designed to enable high fidelity data calibration. The BayesCal data model supplements the a priori known component of the forward model of the sky with a statistical model for the missing and uncertain flux contribution to the data, constrained by a prior on the power in the model. We demonstrate how the parameters of this model can be marginalised out analytically, reducing the dimensionality of the parameter space to be sampled from and allowing one to sample directly from the posterior probability distribution of the calibration parameters. Additionally, we show how physically motivated priors derived from theoretical and measurement-based constraints on the spectral smoothness of the instrumental gains can be used to constrain the calibration solutions. In a companion paper, we apply this algorithm to simulated observations with a HERA-like array and demonstrate that it enables up to four orders of magnitude suppression of power in spurious spectral fluctuations relative to standard calibration approaches.

The search for planets orbiting other stars has recently expanded to include stars from galaxies outside the Milky Way. With the TESS and Gaia surveys, photometric and kinematic information can be combined to identify transiting planet candidates of extragalactic origin. Here, 1,080 low-luminosity red giant branch stars observed by Gaia and TESS with kinematics suggesting a high likelihood of extragalactic origin were searched for planet transits. Transit injection-recovery tests were performed to measure the sensitivity of the TESS data and completeness of the transit search. Injected signals of planets larger than Jupiter with orbital periods of 10 days or less were recovered in $\approx$44% of cases. Although no planet transits were detected in this sample, we find an upper limit on planet occurrence of 0.52% for hot Jupiters, consistent with previous studies of planet occurrence around similar host stars. As stars in the halo tend to be lower metallicity, and short period giant planet occurrence tends to be strongly correlated with stellar metallicity, we predict that relative to the Galactic disk population, a smaller fraction of halo stars will host planets detectable by transit surveys. Thus, applying the known planet occurrence trends to potential planet detection around halo stars, we predict $\gtrsim$7,000 stars must be searched with similar cadence and precision as the stars studied here before a detection of a planet of extragalactic origin is likely. This may be possible with future data releases from the TESS and Gaia missions.

Small high-speed impact craters formed in rocks, ice, and other brittle materials consist of an outer, broad shallow concentric region formed by tensile fracture (spall), surrounding a smaller central "pit" crater of greater depth. On the Earth, that "spall crater" morphology ceases to exist for craters greater than a few meters in diameter. They are not commonly recognized for craters in the solar system but might be an issue for cratering on the small brittle asteroids. We consider the physics of the processes of shock-wave spall cratering and formulate the scaling laws to apply those processes to the bodies of the solar system. Our scaling is based upon analyses of shock-wave propagation and tensile fracture mechanisms, including the important feature of size-dependent tensile fracture, and the role of gravity in lofting spalled material to form the outer parts of the spall craters. We consider the existing scaling laws for cratering in the strength regime and derive the conditions for which spall features will be present or absent. The conditions giving rise to spall cratering are found to be a distinct subset of the 'strength' regime, forming a new sub-regime of cratering. We find that this regime may be very consequential for planetary cratering; in fact, it might dominate all cratering on small rocky asteroids. That has important implications in the interpretation of crater counts and the expected surface effects for rocky, 10-100 km objects.

Short period, massive planets, known as hot Jupiters (HJs), have been discovered around $\sim 1$ percent of local field stars. The inward migration necessary to produce HJs may be `low eccentricity', due to torques in the primordial disc, or `high eccentricity' (HEM). The latter involves exciting high orbital eccentricity, allowing sufficiently close passages with the host star to raise circularising tides in the planet. We present an analytic framework for quantifying the role of dynamical encounters in high density environments during HEM. We show that encounters can enhance or suppress HEM, depending on the local stellar density and the initial semi-major axis $a_0$. For moderate densities, external perturbations can excite large eccentricities that allow a planet to circularise over the stellar lifetime. At extremely high densities, these perturbations can instead result in tidal disruption of the planet, thus yielding no HJ. This may explain the apparent excess of HJs in M67 compared with their local field star abundance versus their apparent deficit in 47 Tuc. Applying our analytic framework, we demonstrate that for an initial massive planet population similar to the field, the expected HJ occurrence rate in 47 Tuc is $f_\mathrm{HJ}=2.2\times 10^{-3}$, which remains consistent with present constraints. Future large (sample sizes $\gtrsim 10^5$) or sensitive transit surveys of stars in globular clusters are required to refute the hypothesis that the initial planet population is similar to the solar neighbourhood average. Non-detection in such surveys would have broad consequences for planet formation theory, implying planet formation rates in globular clusters must be suppressed across a wide range of $a_0$.

We searched for and found no higher mass (>3Msun) unbound binary stellar companions to the progenitor of pulsar J1124-5916. There are lower mass candidates, but they all have high probabilities of being false positives. There are no candidates for it now being a fully unbound triple system. Even if one of the lower mass candidates is an unbound companion, it seems unlikely that it could have contributed to stripping the progenitor prior to the supernova. The stars are too low mass to be significant mass gainers, and they are too slowly moving to be the survivors of a compact, post-common envelope binary. The addition of one more system slightly improves the statistical constraints on the binary and triple status of supernova progenitors just before and after death.

Nuclear reaction rates are a fundamental yet uncertain ingredient in stellar evolution models. The astrophysical S-factor pertaining to the initial reaction in the proton-proton chain is uncertain at the 1% level, which contributes a systematic but generally unpropagated error of similar order in the theoretical ages of stars. In this work, we study the prospect of improving the measurement of this and other reaction rates in the pp chain and CNO cycle using helioseismology and solar neutrinos. We show that when other aspects of the solar model are improved, then it shall be possible using current solar data to improve the precision of this measurement by nearly an order of magnitude, and hence the corresponding uncertainty on the ages of low-mass stars by a similar amount.

Aims. In this work, we aim to make a differential comparison of the neutron-capture and p-process element molybdenum (Mo) in the stellar populations in the local disk(s) and the bulge, focusing on minimising possible systematic effects in the analysis. Methods. The stellar sample consists of 45 bulge and 291 local disk K-giants, observed with high-resolution optical spectra. The abundances are determined by fitting synthetic spectra using the SME-code. The disk sample is separated into thin- and thick-disk components using a combination of abundances and kinematics. The cosmic origin of Mo is investigated and discussed by comparing with previous published abundances of Mo and the neutron-capture elements cerium (Ce) and europium (Eu). Results. We determine reliable Mo abundances for 35 bulge and 282 disk giants with a typical uncertainty of [Mo/Fe]~0.2 and ~0.1 dex for the bulge and disk, respectively. Conclusions. We find that the bulge possibly is enhanced in [Mo/Fe] compared to the thick disk, which we do not observe in either [Ce/Fe] nor [Eu/Fe]. This might suggest a higher past star-formation rate in the bulge, however, since we do not observe the bulge to be enhanced in [Eu/Fe], the origin of the molybdenum enhancement is yet to be constrained. Although, the scatter is large, we may be observing evidence of the p-process contributing to the heavy element production in the chemical evolution of the bulge.

The prediction of solar flares is still a significant challenge in space weather research, with no techniques currently capable of producing reliable forecasts performing significantly above climatology. In this paper, we present a flare forecasting technique using data assimilation coupled with computationally inexpensive cellular automata called sandpile models. Our data assimilation algorithm uses the simulated annealing method to find an optimal initial condition that reproduces well an energy-release time series. We present and empirically analyze the predictive capabilities of three sandpile models, namely the Lu and Hamilton model (LH) and two deterministically-driven models (D). Despite their stochastic elements, we show that deterministically-driven models display temporal correlations between simulated events, a needed condition for data assimilation. We present our new data assimilation algorithm and demonstrate its success in assimilating synthetic observations produced by the avalanche models themselves. We then apply our method to GOES X-Ray time series for 11 active regions having generated multiple X-class flares in the course of their lifetime. We demonstrate that for such large flares, our data assimilation scheme substantially increases the success of ``All-Clear'' forecasts, as compared to model climatology.

The temporal and spectral evolution of gamma-ray burst (GRB) afterglows can be used to infer the density and density profile of the medium through which the shock is propagating. In long-duration (core-collapse) GRBs, the circumstellar medium (CSM) is expected to resemble a wind-blown bubble, with a termination shock separating the stellar wind and the interstellar medium (ISM). A long standing problem is that flat density profiles, indicative of the ISM, are often found at lower radii than expected for a massive star progenitor. Furthermore, the presence of both wind-like environments at high radii and ISM-like environments at low radii remains a mystery. In this paper, we perform a 'CSM population synthesis' with long GRB progenitor stellar evolution models. Analytic results for the evolution of wind blown bubbles are adjusted through comparison with a grid of 2D hydrodynamical simulations. Predictions for the emission radii, ratio of ISM to wind-like environments, wind and ISM densities are compared with the largest sample of afterglow-derived parameters yet compiled, which we make available for the community. We find that high ISM densities around 1000/cm3 best reproduce observations. If long GRBs instead occur in typical ISM densities of approximately 1/cm3, then the discrepancy between theory and observations is shown to persist at a population level. We discuss possible explanations for the origin of variety in long GRB afterglows, and for the overall trend of CSM modelling to over-predict the termination shock radius.

Low mass X-ray binaries (LMXBs) show strong variability over a broad range of time scales. The analysis of this variability, in particular of the quasi-periodic oscillations (QPO), is key to understanding the properties of the innermost regions of the accretion flow in these systems. We present a time-dependent Comptonisation model that fits the energy-dependent rms-amplitude and phase-lag spectra of low-frequency QPOs in black-hole (BH) LMXBs. We model the accretion disc as a multi-temperature blackbody source emitting soft photons which are then Compton up-scattered in a spherical corona, including feedback of Comptonised photons that return to the disc. We compare our results with those obtained with a model in which the seed-photons source is a spherical blackbody: at low energies the time-averaged, rms and phase-lag spectra are smoother for the disk-blackbody than for a blackbody, while at high energies both models give similar spectra. In general, we find that the rms increases with energy, the slope of the phase-lag spectrum depends strongly on the feedback, while the minimum-lag energy is correlated with the disc temperature. We fit the model to a 4.45-Hz type-B QPO in the BH LMXB MAXI J1438-630 and find statistically-better fits and more compatible parameters with the steady-state spectrum than those obtained with a blackbody seed-photons source. Furthermore, we successfully apply the model to the type-C QPO in the BH LMXB GRS 1915+105, and thus conclude that this variable-Comptonisation model reproduces the rms and phase-lags of both type B and C low-frequency QPOs in BH LMXBs.

Several hundred gamma-ray burst (GRB) redshifts have been determined to date. One of the other important properties-besides the distance-of the GRBs is the duration of the burst. In this paper, we analyse these two important quantities of the phenomena. In this paper, we map the two-dimensional distribution and explore some suspicious areas. As it is well known that the short GRBs are closer than the others, we search for parts in the Universe where the GRB duration is different from the others. We also analyse whether there are any ranges in the duration where the redshifts differ. We find some suspicious areas, however, no other significant region was found than the short GRB region.

Photometric redshifts are essential in studies of both galaxy evolution and cosmology, as they enable analyses of objects too numerous or faint for spectroscopy. The Rubin Observatory, Euclid, and Roman Space Telescope will soon provide a new generation of imaging surveys with unprecedented area coverage, wavelength range, and depth. To take full advantage of these datasets, further progress in photometric redshift methods is needed. In this review, we focus on the greatest common challenges and prospects for improvement in applications of photo-$z$'s to the next generation of surveys: - Gains in $performance$ -- i.e., the precision of redshift estimates for individual galaxies -- could greatly enhance studies of galaxy evolution and some probes of cosmology. - Improvements in $characterization$ -- i.e., the accurate recovery of redshift $distributions$ of galaxies in the presence of uncertainty on individual redshifts -- are urgently needed for cosmological measurements with next-generation surveys. - To achieve both of these goals, improvements in the scope and treatment of the samples of spectroscopic redshifts which make high-fidelity photo-$z$'s possible will also be needed. For the full potential of the next generation of surveys to be reached, the characterization of redshift distributions will need to improve by roughly an order of magnitude compared to the current state of the art, requiring progress on a wide variety of fronts. We conclude by presenting a speculative evaluation of how photometric redshift methods and the collection of the necessary spectroscopic samples may improve by the time near-future surveys are completed.

The rapid increase in serendipitous X-rays source detections requires the development of novel approaches to efficiently explore the nature of X-ray sources. If even a fraction of these sources could be reliably classified, it would enable population studies for various astrophysical source types on a much larger scale than currently possible. Classification of large numbers of sources from multiple classes characterized by multiple properties (features) must be done automatically and supervised machine learning (ML) seems to provide the only feasible approach. We perform classification of Chandra Source Catalog version 2 (CSCv2) sources to explore the potential of the ML approach and identify various biases, limitations, and bottlenecks that present themselves in these kinds of studies. We establish the framework and present a flexible and expandable Python pipeline, which can be used and improved by others. We also release the training dataset of 2,941 X-ray sources with confidently established classes. In addition to providing probabilistic classifications of 66,369 CSCv2 sources (21% of the entire CSCv2 catalog), we perform several narrower-focused case studies (high-mass X-ray binary candidates and X-ray sources within the extent of the H.E.S.S. TeV sources) to demonstrate some possible applications of our ML approach. We also discuss future possible modifications of the presented pipeline, which are expected to lead to substantial improvements in classification confidences.

We have carried out Chandra, HST, and VLA observations of four MOJAVE blazars that have been previously classified as 'hybrid' (FR I/II) blazars in terms of radio morphology but not total radio power. The motivation of this study is to determine the X-ray emission mechanism in jets, these being different in FR I and FR II jets. We detected X-ray jet emission with sufficient SNR in two blazars viz. PKS 0215+015 and TXS 0730+504. We carried out spectral energy distribution (SED) modeling of the broad-band emission from the jet regions in these sources and found that a single synchrotron emission model is ruled out due to the deep upper limits obtained from HST optical and IR data. The IC- CMB model can reproduce the X-ray jet emission in both sources although the model requires extreme jet parameters. Both our sources possess FR II like radio powers and our results are consistent with previous studies suggesting that radio power is more important than FR morphology in determining the emission mechanism of X-ray jets.

Motivated by the newly discovery of a millisecond pulsar (MSP) J1835$-$3259B in the globular cluster (GC) NGC~6652, we analyze the gamma-ray data obtained with the Large Area Telescope (LAT) onboard {\it Fermi Gamma-ray Space Telescope (Fermi)} for the GC. From timing analysis of the data, a pulse profile that is similar to the radio one is obtained. Also the weighted H-test value from the analysis is $\sim$41, corresponding to a chance probability of $7.5\times 10^{-8}$ ($\simeq 5.4\sigma$). We thus consider that we have detected the gamma-ray emission of the MSP, and discuss the implications. Based on different studies of the sources in the GC, the observed gamma-ray emission from the GC could mainly arise from this MSP, like the previous two cases the GCs NGC~6624 and NGC~6626. If this is the case, the pulsar would have a gamma-ray luminosity of $\sim 5\times 10^{34}$\,erg\,s$^{-1}$ and a gamma-ray efficiency of $\sim 0.12$.

We put forward a pressure-parametric model to study the tiny deviation from cosmological constant(CC) behavior of the dark sector accelerating the expansion of the Universe. Data from cosmic microwave background (CMB) anisotropies, baryonic acoustic oscillations (BAO), Type Ia supernovae (SN Ia) observation are applied to constrict the model parameters. The constraint results show that such model suffers with $H_0$ tension as well. To realize this model more physically, we reconstruct it with the quintessence and phantom scalar fields, and find out that although the model predicts a quintessence-induced acceleration of the Universe at past and present, at some moment of the future, dark energy's density have a disposition to increase.

Galaxy photometric redshift (photo-$z$) is crucial in cosmological studies, such as weak gravitational lensing and galaxy angular clustering measurements. In this work, we try to extract photo-$z$ information and construct its probabilistic distribution function (PDF) using the Bayesian neural networks (BNN) from both galaxy flux and image data expected to be obtained by the China Space Station Telescope (CSST). The mock galaxy images are generated from the Advanced Camera for Surveys of Hubble Space Telescope (HST-ACS) and COSMOS catalog, in which the CSST instrumental effects are carefully considered. And the galaxy flux data are measured from galaxy images using aperture photometry. We construct Bayesian multilayer perceptron (B-MLP) and Bayesian convolutional neural network (B-CNN) to predict photo-$z$ along with the PDFs from fluxes and images, respectively. We combine the B-MLP and B-CNN together, and construct a hybrid network and employ the transfer learning techniques to investigate the improvement of including both flux and image data. We find that the accuracy and outlier fraction of photo-$z$ can achieve $\sigma_{NMAD}$ = 0.022 and $\eta$ = 2.83% for the B-MLP using the flux data only, and $\sigma_{NMAD}$ = 0.025 and $\eta$ = 2.32% for the B-CNN using the image data only. The Bayesian hybrid transfer network can improve the result to $\sigma_{NMAD}$ = 0.021 and $\eta$ = 1.62%, and can provide the most confident predictions with the lowest average uncertainty.

Transient sources such as supernovae (SNe) and tidal disruption events are candidates for high energy neutrino sources. However, SNe commonly occur in the Universe and a chance coincidence of their detection to a neutrino signal cannot be avoided, which may lead to a challenge of claiming their association with neutrino emissions. In order to overcome this difficulty, we propose a search for $\sim10-100$ TeV neutrino multiple events within a timescale of $\sim 30$ days coming from the same direction, called neutrino multiplets. We show that demanding multiplet detection by a $\sim 1$ km$^3$ neutrino telescope limits distances of detectable neutrino sources, which enables us to identify source counterparts by multiwavelength observations owing to the substantially reduced rate of the chance coincidence detection of transients. We apply our results to construct a feasible strategy for optical followup observations and demonstrate that wide-field optical telescopes with a $\gtrsim4$ m dish should be capable of identifying a transient associated with a neutrino multiplet. We also present the resultant sensitivity of multiplet neutrino detection as a function of the released energy of neutrinos and burst rate density. A model of neutrino transient sources with emission energy greater than ${\rm a~few}\times 10^{51}$ erg and burst rate rarer than ${\rm a~few}\times 10^{-8}\ {\rm Mpc}^{-3}\ {\rm yr}^{-1}$ is constrained by the null detection of multiplets by a $\sim 1$ km$^3$ scale neutrino telescope. This already disfavors the canonical high-luminosity gamma ray bursts and jetted tidal disruption events as major sources in the TeV-energy neutrino sky.

Presolar grains are microscopic dust grains that formed in the stellar winds or explosions of ancient stars that died before the formation of the solar system. The majority (~90% in number) of presolar silicon carbide (SiC) grains, including types mainstream (MS), Y, and Z, came from low-mass C-rich asymptotic giant branch (AGB) stars, which is supported by the ubiquitous presence of SiC dust observed in the circumstellar envelope of AGB stars and the signatures of slow neutron-capture process preserved in these grains. Here, we review the status of isotope studies of presolar AGB SiC grains with an emphasis on heavy-element isotopes and highlight the importance of presolar grain studies for nuclear astrophysics. We discuss the sensitives of different types of nuclei to varying AGB stellar parameters and how their abundances in presolar AGB SiC grains can be used to provide independent, detailed constraints on stellar parameters, including 13C formation, stellar temperature, and nuclear reaction rates.

Of all the light elements, the evolution of lithium (Li) in the Milky Way is perhaps the most difficult to explain. Li is difficult to synthesize and is easily destroyed, making most stellar sites unsuitable for producing Li in sufficient quantities to account for the proto-solar abundance. For decades, novae have been proposed as a potential explanation to this 'Galactic Li problem', and the recent detection of 7Be in the ejecta of multiple nova eruptions has breathed new life into this theory. In this work, we assess the viability of novae as dominant producers of Li in the Milky Way. We present the most comprehensive treatment of novae in a galactic chemical evolution code to date, testing theoretical- and observationally-derived nova Li yields by integrating metallicity-dependent nova ejecta profiles computed using the binary population synthesis code binary c with the galactic chemical evolution code OMEGA+. We find that our galactic chemical evolution models which use observationally-derived Li yields account for the proto-solar Li abundance very well, while models relying on theoretical nova yields cannot reproduce the proto-solar observation. A brief exploration of physical uncertainties including single-stellar yields, the metallicity resolution of our nova treatment, common-envelope physics, and nova accretion efficiencies indicates that this result is robust to physical assumptions. Scatter within the observationally-derived Li yields in novae is identified as the primary source of uncertainty, motivating further observations of 7Be in nova ejecta.

High luminosity accretion onto a strongly magnetized neutron star results in a radiation pressure dominated, magnetically confined accretion column. We investigate the dynamics of these columns using two-dimensional radiative relativistic magnetohydrodynamic simulations, restricting consideration to modest accretion rates where the height of the column is low enough that Cartesian geometry can be employed. The column structure is dynamically maintained through high-frequency oscillations of the accretion shock at $\simeq 10-25$~kHz. These oscillations arise because it is necessary to redistribute the power released at the accretion shock through bulk vertical motions, both to balance the cooling and to provide vertical pressure support against gravity. Sideways cooling always dominates the loss of internal energy. In addition to the vertical oscillations, photon bubbles form in our simulations and add additional spatial complexity to the column structure. They are not themselves responsible for the oscillations, and they do not appear to affect the oscillation period. However, they enhance the vertical transport of radiation and increase the oscillation amplitude in luminosity. The time-averaged column structure in our simulations resembles the trends in standard 1D stationary models, the main difference being that the time-averaged height of the shock front is lower because of the higher cooling efficiency of the 2D column shape.

Observations of the rest-frame UV emission of high-redshift galaxies suggest that the early stages of galaxy formation involve disturbed structures. Imaging the cold interstellar medium can provide a unique view of the kinematics associated with the assembly of galaxies. In this paper, we analyzed the spatial distribution and kinematics of the cold ionized gas of the normal star-forming galaxy COS-2987030247 at z = 6.8076, based on new high-resolution observations of the [C II] 158um line emission obtained with the Atacama Large Millimeter/submillimeter Array. These observations allowed us to compare the spatial distribution and extension of the [C II] and rest-frame UV emission, model the [C II] line data-cube using 3DBarolo, and measure the [C II] luminosity and star formation rate (SFR) surface densities in the galaxy subregions. The system is found to be composed of a main central source, a fainter north extension, and candidate [C II] companions located 10-kpc away. We find similar rest-frame UV and [C II] spatial distributions, suggesting that the [C II] emission emerges from the star-forming regions. The agreement between the UV and [C II] surface brightness radial profiles rules out diffuse, extended [C II] emission in the main galaxy component. The [C II] velocity map reveals a velocity gradient in the north-south direction suggesting ordered motion, as commonly found in rotating-disk galaxies. But higher-resolution observations would be needed to rule out a compact merger scenario. Our model indicates a low average velocity dispersion, $\sigma$ < 30 km s$^{-1}$. This result implies a dispersion lower than the expected value from observations and semi-analytic models of high redshift galaxies. We argue that COS-2987030247 is a candidate rotating disk experiencing a short period of stability which will be possibly perturbed at later times by accreting sources.

Magnetic fields and turbulence are important components of the interstellar medium (ISM) of star-forming galaxies. It is challenging to measure the properties of the small-scale ISM magnetic fields (magnetic fields at scales smaller than the turbulence driving scale). Using numerical simulations, we demonstrate how the second-order rotation measure (RM, which depends on thermal electron density, $n_{\rm e}$, and magnetic field, $b$) structure function can probe the properties of small-scale $b$. We then apply our results to observations of the Small and Large Magellanic Clouds (SMC and LMC). First, using Gaussian random $b$, we show that the characteristic scale where the RM structure function flattens is approximately equal to the correlation length of $b$. We also show that computing the RM structure function with a higher-order stencil (more than the commonly-used two-point stencil) is necessary to accurately estimate the slope of the structure function. Then, using Gaussian random $b$ and lognormal $n_{\rm e}$ with known power spectra, we derive an empirical relationship between the slope of the power spectrum of $b$, $n_{\rm e}$, and RM. We apply these results to the SMC and LMC and estimate the following properties of small-scale $b$: correlation length ($160~\pm 21~{\rm pc}$ for the SMC and $87~\pm~17~{\rm pc}$ for the LMC), strength ($14~\pm 2~\mu{\rm G}$ for the SMC and $15~\pm 3~\mu{\rm G}$ for the LMC), and slope of the magnetic power spectrum ($-1.3~\pm~0.4$ for the SMC and $-1.6~\pm~0.1$ for the LMC). We also find that $n_{\rm e}$ is practically constant over the estimated $b$ correlation scales.

We spoke with four researchers to understand the accessibility challenges in astronomy research, education and outreach for blind and visually impaired (BVI) persons, as well as solutions to these challenges and how it innovates data analysis methods for all astronomers. Those interviewed: Nicolas Bonne (University of Portsmouth); Cheryl Fogle-Hatch (Museum Senses); Garry Foran (Swinburne University of Technology) and Enrique Perez Montero (Instituto de Astrof\'isica de Andaluc\'ia).

The life and death of a star are controlled by its internal rotation dynamics through subtle transport and mixing mechanisms, which so far remain poorly understood. While magnetic fields must play a crucial role in transporting angular momentum and chemical species, the very origin of magnetism in radiative stellar layers and its influence on spinning dynamics are yet to be unraveled. Using global numerical modeling, we report the existence of a dynamo sharing many characteristics with the (never observed) Tayler-Spruit model, which can generate strong magnetic fields and significantly enhance transport in radiative zones. The resulting, deep toroidal fields are screened by the outer medium, allowing for the existence of intense magnetism in radiative stars where no magnetic fields could be directly observed so far.

Combining Supernovae, Baryon Acoustic Oscillations and Redshift-Space Distortions data from the next generation of (Stage-IV) cosmological surveys, we aim to reconstruct the expansion history up to large redshifts using forward-modeling of $f_{\mathrm DE}(z) = \rho_\mathrm{DE}(z)/\rho_\mathrm{DE,0}$ with Gaussian processes (GP). In order to reconstruct cosmological quantities at high redshifts where few or no data are available, we adopt a new approach to GP which enforces the following minimal assumptions: a) Our cosmology corresponds to a flat Friedman-Lema\^itre-Robertson-Walker (FLRW) universe; b) An Einstein de Sitter (EdS) universe is obtained on large redshifts. This allows us to reconstruct the perturbations growth history from the reconstructed background expansion history. Assuming various DE models, we show the ability of our reconstruction method to differentiate them from $\Lambda$CDM at $\gtrsim2\sigma$.

The fine magnetic structure is vitally important to understanding the formation, stabilization and eruption of solar filaments, but so far, it is still an open question yet to be resolved. Using stereoscopic observations taken by the Solar Dynamics Observatory and Solar TErrestrial RElations Obsevatory, we studied the generation mechanism of a two-sided-loop jet (TJ) and the ejection process of the jet plasma into the overlying filament-cavity system. We find that the generation of the two-sided-loop jet was due to the magnetic reconnection between an emerging flux loop and the overlying filament. The jet's two arms ejected along the filament axis during the initial stage. Then, the north arm bifurcated into two parts at about 50 Mm from the reconnection site. After the bifurcation, the two bifurcated parts were along the filament axis and the cavity which hosted the filament, respectively. By tracing the ejecting plasma flows of the TJ inside the filament, we not only measured that the magnetic twist stored in the filament was at least 5$\pi$ but also found that the fine magnetic structure of the filament-cavity flux rope system is in well agreement with the theoretical results of Magnetic flux rope models.

Aims. The Chemical Evolution of R-process Elements in Stars (CERES) project aims to provide a homogeneous analysis of a sample of metal-poor stars ([Fe/H]<-1.5). We present the stellar parameters and the chemical abundances of elements up to Zr for a sample of 52 giant stars.Methods. We relied on a sample of high signal-to-noise UVES spectra. We determined stellar parameters from Gaia photometry and parallaxes. Chemical abundances were derived using spectrum synthesis and model atmospheres.Results. We determined chemical abundances of 26 species of 18 elements: Na, Mg, Al, Si, Ca, Sc, Ti, V, Cr, Mn, Fe, Co, Ni, Cu, Zn, Sr, Y, and Zr. For several stars, we were able to measure both neutral and ionised species, including Si, Sc, Mn, and Zr. We have roughly doubled the number of measurements of Cu for stars at [Fe/H] <= -2.5. The homogeneity of the sample made it possible to highlight the presence of two Zn-rich stars ([Zn/Fe]~+0.7), one r-rich and the other r-poor. We report the existence of two branches in the [Zn/Fe] versus [Ni/Fe] plane and suggest that the high [Zn/Fe] branch is the result of hypernova nucleosynthesis. We discovered two stars with peculiar light neutron-capture abundance patterns: CES1237+1922 (also known as BS 16085-0050), which is ~1 dex underabundant in Sr, Y, and Zr with respect to the other stars in the sample, and CES2250-4057 (also known as HE 2247-4113), which shows a ~1 dex overabundance of Sr with respect to Y and Zr.Conclusions. The high quality of our dataset allowed us to measure hardly detectable ions. This can provide guidance in the development of line formation computations that take deviations from local thermodynamic equilibrium and hydrodynamical effects into account.

Current wisdom accounts to the diversity of neutron star observational manifestations to their birth scenarios, influencing their thermal and magnetic field evolution. Among the kind of observed neutron stars, radio pulsars represent by far the largest population of neutron stars. In this paper, we aim at constraining the observed population of the canonical neutron star period, magnetic field and spatial distribution at birth in order to understand the radio and high-energy emission processes in a pulsar magnetosphere. For this purpose we design a population synthesis method self-consistently taking into account the secular evolution of a force-free magnetosphere and the magnetic field decay. We generate a population of pulsars and evolve them from their birth to the present time, working in the force-free approximation. We assume a given initial distribution for the spin period, surface magnetic field and spatial galactic location. Radio emission properties are accounted by the polar cap geometry, whereas the gamma-ray emission is assumed to be produced within the striped wind model. We found that a decaying magnetic field gave better agreement with observations compared to a constant magnetic field model. Starting from an initial mean magnetic field strength of $B=2.5\times 10^8$~T with a characteristic decay timescale of $4.6 \times 10^5~$yr, a neutron star birth rate of $1/70~$yr and a mean initial spin period of 60~ms, we found that the force-free model satisfactorily reproduces the distribution of pulsars in the $P-\dot{P}$ diagram with simulated populations of radio-loud, radio-only and radio quiet gamma-ray pulsars similar to the observed populations.

With the remarkable success of the LVK consortium in detecting binary black hole mergers, it has become possible to use the population properties to constrain our understanding of the progenitor stars' evolution. The most striking features of the observed primary black hole mass distributions are the extended tail up to 100M$_\odot$ and an excess of masses at 35M$_\odot$. Currently, isolated binary population synthesis have difficulty explaining these features. Using the well-tested BPASS detailed stellar binary evolution models to determine mass transfer stability, accretion rates, and remnant masses, we postulate that stable mass transfer with super-Eddington accretion is responsible for the extended tail. Furthermore, that the excess is not due to pulsation-pair instability, as previously thought, but due to stable mass transfer. These systems are able to merge within the Hubble time due to more stable mass transfer with extreme mass ratios that allows the orbits to shrink sufficiently to allow for a merger. These finding are at odds with those from other population synthesis codes but in agreement with other recent studies using detailed binary evolution models.

Cold super-Earths which retain their primordial, H-He dominated atmosphere could have surfaces that are warm enough to host liquid water. This would be due to the collision induced absorption (CIA) of infra-red light by hydrogen, which increases with pressure. However, the long-term potential for habitability of such planet has not been explored yet. Here we investigate the duration of this potential exotic habitability by simulating planets of different core masses, envelope masses and semi-major axes. We find that terrestrial and super-Earth planets with masses of $\sim$1 - 10$M_{\oplus}$ can maintain temperate surface conditions up to 5 - 8 Gyr at radial distances larger than $\sim $2 AU. The required envelope masses are $\sim 10^{-4} \, M_{\oplus}$ (which is 2 orders of magnitude more massive than Earth's), but can be an order of magnitude smaller (when close-in) or larger (when far out). This result suggests that the concept of planetary habitability should be revisited and made more inclusive with respect to the classical definition.

We designed and built a novel model of a deployed space telescope which can reliably align its segments to achieve the finest possible resolution. An asymmetric design of both the segment shapes and their pupil locations were tested in simulation and experiment. We optimised the sparse aperture for better spatial frequency coverage and for smoother images with less artifacts. The unique segment shapes allow for an easier identification and alignment, and the feedback is based only upon the focal image. The autonomous alignment and fine tuning are governed by mechanical simplicity and reliability.

Axisymmetric dust rings containing tens to hundreds of Earth masses of solids have been observed in protoplanetary discs with (sub-)millimetre imaging. Here, we investigate the growth of a planetary embryo in a massive (150M$_\oplus$) axisymmetric dust trap through dust and gas hydrodynamics simulations. When accounting for the accretion luminosity of the planetary embryo from pebble accretion, the thermal feedback on the surrounding gas leads to the formation of an anticyclonic vortex. Since the vortex forms at the location of the planet, this has significant consequences for the planet's growth: as dust drifts towards the pressure maximum at the centre of the vortex, which is initially co-located with the planet, a rapid accretion rate is achieved, in a distinct phase of ``vortex-assisted'' pebble accretion. Once the vortex separates from the planet due to interactions with the disc, it accumulates dust, shutting off accretion onto the planet. We find that this rapid accretion, mediated by the vortex, results in a planet containing $\approx$ 100M$_\oplus$ of solids. We follow the evolution of the vortex, as well as the efficiency with which dust grains accumulate at its pressure maximum as a function of their size, and investigate the consequences this has for the growth of the planet as well as the morphology of the protoplanetary disc. We speculate that this extreme formation scenario may be the origin of giant planets which are identified to be significantly enhanced in heavy elements.

In this paper, we present a sample of 21 repeating fast radio bursts (FRBs) detected by different radio instruments before September 2021. Using the Anderson--Darling test, we compared the distributions of extra-Galactic dispersion measure ($DM_{\rm E}$) of non-repeating FRBs, repeating FRBs and all FRBs. It was found that the $ DM_{\rm E}$ values of three sub-samples are log-normally distributed. The $DM_{\rm E}$ of repeaters and non-repeaters were drawn from a different distribution on basis of the Mann--Whitney--Wilcoxon test. In addition, assuming that the non-repeating FRBs identified currently may be potentially repeators, i.e., the repeating FRBs to be universal and representative, one can utilize the averaged fluence of repeating FRBs as an indication from which to derive an apparent intensity distribution function (IDF) with a power-law index of $a_1=$ $1.10\pm 0.14$ ($a_2=$ $1.01\pm 0.16$, the observed fluence as a statistical variant), which is in good agreement with the previous IDF of 16 non-repeating FRBs found by Li et al. Based on the above statistics of repeating and non-repeating FRBs, we propose that both types of FRBs may have different cosmological origins, spatial distributions and circum-burst environments. Interestingly, the differential luminosity distributions of repeating and non-repeating FRBs can also be well described by a broken power-law function with the same power-law index of $-$1.4.

The research of infall motion is a common means to study molecular cloud dynamics and the early process of star formation. Many works had been done in-depth research on infall. We searched the literature related to infall study of molecular cloud since 1994, summarized the infall sources identified by the authors. A total of 456 infall sources are catalogued. We classify them into high-mass and low-mass sources, in which the high-mass sources are divided into three evolutionary stages: prestellar, protostellar and {H\small \hspace{0.1em}II} region. We divide the sources into clumps and cores according to their sizes. The H$_2$ column density values range from 1.21$\times$ 10$^{21}$ to 9.75 $\times$ 10$^{24}$ cm$^{-2}$, with a median value of 4.17$\times$ 10$^{22}$ cm$^{-2}$. The H$_2$ column densities of high-mass and low-mass sources are significantly separated. The median value of infall velocity for high-mass clumps is 1.12 km s$^{-1}$, and the infall velocities of low-mass cores are virtually all less than 0.5 km s$^{-1}$. There is no obvious difference between different stages of evolution. The mass infall rates of low-mass cores are between 10$^{-7}$ and 10$^{-4}$ M$_{\odot} \text{yr}^{-1}$, and those of high-mass clumps are between 10$^{-4}$ and 10$^{-1}$ M$_{\odot} \text{yr}^{-1}$ with only one exception. We do not find that the mass infall rates vary with evolutionary stages.

The design of mission scenarios for the flyby investigation of nearby star systems by probes launched using directed energy is addressed. Multiple probes are launched with a fixed launch infrastructure, and download of scientific data occurs following target encounter and data collection. Assuming the primary goal is to reliably recover a larger volume of collected scientific data with a smaller data latency (elapsed time from launch to complete recovery of the data), it is shown that there is an efficient frontier where volume cannot be increased for a given latency and latency cannot be reduced for a given volume. For each probe launch, increasing the volume along this frontier is achieved by increasing the probe mass, which results in a reduced probe speed. Thus choosing the highest feasible probe speed generally does not achieve an efficient tradeoff of volume and latency. Along this frontier the total distance traveled to the completion of data download does not vary significantly, implying that the download time duration is approximately a fixed fraction of the launch-to-target transit time. Due to longer propulsion duration when probe mass is increased, increasing data volume incurs a cost in the total launch energy expended, but with favorable economies of scale. An important characteristic of any probe technology is the scaling law that relates probe mass to transmit data rate, as this affects details of the efficient frontier.

In this work, we aimed at investigating the star formation rate of active galactic nuclei host galaxies in the early Universe. To this end, we constructed a sample of 149 luminous ($\rm L_{2-10keV} > 10^{44}\,erg\,s^{-1}$) X-ray AGNs at $\rm z \geq3.5$ selected in three fields with different depths and observed areas (Chandra COSMOS Legacy survey, XMM-XXL North and eFEDS). We built their spectral energy distributions (SED) using available multi-wavelength photometry from X-rays up to far-IR. Then, we estimated the stellar mass, M$_{*}$, and the SFR of the AGNs using the X-CIGALE SED fitting algorithm. After applying several quality criteria, we ended up with 89 high-z sources. More than half (55\%) of the X-ray sample have spectroscopic redshifts. Based on our analysis, our high-z X-ray AGNs live in galaxies with median $\rm M_{*}=5.6 \times10^{10}~M_\odot$ and $\rm SFR_{*}\approx240\,M_\odot yr^{-1}$. The majority of the high-z sources ($\sim89$\%) were found inside or above the main sequence (MS) of star-forming galaxies. Estimation of the normalised SFR, $\rm SFR_{NORM}$, defined as the ratio of the SFR of AGNs to the SFR of MS galaxies, showed that the SFR of AGNs is enhanched by a factor of $\sim 1.8$ compared to non-AGN star-forming systems. Combining our results with previous studies at lower redshifts, we confirmed that $\rm SFR_{NORM}$ does not evolve with redshift. Using the specific BHAR (i.e., $\rm L_X$ divided by $\rm M_{*}$), $\rm \lambda _{BHAR}$, that can be used as a tracer of the Eddington ratio, we found that the bulk of AGNs that lie inside or above the MS have higher specific accretion rates compared to sources below the MS. Finally, we found indications that the SFR of the most massive AGN host galaxies ($\rm log\,(M_{*}/ M_\odot) >10^{11.5-12}$) remains roughly constant as a function of M$_*$, in agreement with the SFR of MS star-forming galaxies.

Future far-infrared space missions require highly-sensitive spectroscopy as a primary diagnostics tool. However, these systems are sensitive to straylight, due to the ultra-sensitive few-mode detectors used, which affects the measurement and calibration of the spectrum, as revealed by the Herschel mission. To ensure that the science goals of future missions are met, the complex modal behaviour has to be understood, and appropriate verification and calibration strategies must be developed. We propose a modal framework to addresses these issues, using Herschel-SPIRE as a case study, and demonstrate how the technique can be used for the design and verification of spectrometers in future far-infrared missions.

We use panoramic optical spectroscopy obtained with MUSE@VLT to investigate the nature of five candidate extremely isolated low-mass star forming regions (Blue Candidates, BCs hereafter) toward the Virgo cluster of galaxies. Four of the five (BC1, BC3, BC4, BC5) are found to host several HII regions and to have radial velocities fully compatible with being part of the Virgo cluster. All the confirmed candidates have mean metallicity significantly in excess of that expected from their stellar mass, indicating that they originated from gas stripped from larger galaxies. In summary, these four candidates share the properties of the prototype system SECCO 1, suggesting the possible emergence of a new class of stellar systems, intimately linked to the complex duty cycle of gas within clusters of galaxies. A thorough discussion on the nature and evolution of these objects is presented in a companion paper, where the results obtained here from MUSE data are complemented with Hubble Space Telescope (optical) and Very Large Array (HI) observations.

We present Hubble Space Telescope WFC3/UVIS (F275W, F336W) and ACS/WFC optical (F435W, F555W, and F814W) observations of the nearby grand-design spiral galaxy M101 as part of the Legacy Extragalactic UV Survey (LEGUS). Compact sources detected in at least four bands were classified by both human experts and the convolutional neural network StarcNet. Human experts classified the 2,351 brightest sources, retrieving $N_{c} = 965$ star clusters. StarcNet, trained on LEGUS data not including M101, classified all 4,725 sources detected in four bands, retrieving $N_{c} = 2,270$ star clusters. The combined catalog represents the most complete census to date of compact star clusters in M101. We find that for the 2,351 sources with both a visual- and ML-classification StarcNet is able to reproduce the human classifications at high levels of accuracy ($\sim 80-90\%$), which is equivalent to the level of agreement between human classifiers in LEGUS. The derived cluster age distribution implies a disruption rate of $dN/d\tau \propto \tau^{-0.45 \pm 0.14}$ over $10^{7} < \tau < 10^{8.5}$yr for cluster masses $\geq 10^{3.55} M_{\odot}$ for the central region of M101 and $dN/d\tau \propto \tau^{-0.02 \pm 0.15}$ for cluster masses $\geq 10^{3.38} M_{\odot}$ in the northwest region of the galaxy. The trends we recover are weaker than those of other nearby spirals (e.g. M51) and starbursts, consistent with the M101 environment having a lower-density interstellar medium, and providing evidence in favor of environmentally-dependent cluster disruption in the central, southeast, and northwest regions of M101.

Multi-messenger astronomy is a vast and expanding field as electromagnetic observations (EM) are no longer the only way of exploring the Universe. Due to the new messengers, astrophysical events with both gravitational waves (GWs) and EM emission are no longer a dream of the astronomical community. A breakthrough for GW multi-messenger astronomy came when the LIGO-Virgo network detected a GW signal of two low-mass compact objects consistent with a binary neutron star (BNS, GW170817) an event that generated a short gamma-ray burst (GRBs) and a kilonova. While GW170817 represents the testimony to BNS mergers being the progenitor of at least some GRBs, a wide range of highly energetic astrophysical phenomena is expected to be accompanied by the emission of GWs and photons. Here we present an unmodelled method to search for GWs having gamma and radio counterparts, using the LIGO-Virgo data and observations of partner telescopes. We also discuss the most recent results of the unmodelled coherent searches targeting astrophysical events during the first part of the LIGO-Virgo third observing run (O3a): 105 GRBs detected by the Fermi and Swift satellites.

PAUCam is an innovative optical narrow-band imager mounted at the William Herschel Telescope built for the Physics of the Accelerating Universe Survey (PAUS). Its set of 40 filters results in images that are complex to calibrate, with specific instrumental signatures that cannot be processed with traditional data reduction techniques. In this paper we present two pipelines developed by the PAUS data management team with the objective of producing science-ready catalogues from the uncalibrated raw images. The Nightly pipeline takes care of all image processing, with bespoke algorithms for photometric calibration and scatter-light correction. The Multi-Epoch and Multi-Band Analysis (MEMBA) pipeline performs forced photometry over a reference catalogue to optimize the photometric redshift performance. We verify against spectroscopic observations that the current approach delivers an inter-band photometric calibration of 0.8% across the 40 narrow-band set. The large volume of data produced every night and the rapid survey strategy feedback constraints require operating both pipelines in the Port d'Informaci\'o Cientifica data centre with intense parallelization. While alternative algorithms for further improvements in photo-z performance are under investigation, the image calibration and photometry presented in this work already enable state-of-the-art photometric redshifts down to i_{AB}=22.5 .

Dedicated surveys searching for Fast Radio Bursts (FRBs) are subject to selection effects which bias the observed population of events. Software injection systems are one method of correcting for these biases by injecting a mock population of synthetic FRBs directly into the realtime search pipeline. The injected population may then be used to map intrinsic burst properties onto an expected signal-to-noise ratio (SNR), so long as telescope characteristics such as the beam model and calibration factors are properly accounted for. This paper presents an injection system developed for the Canadian Hydrogen Intensity Mapping Experiment Fast Radio Burst project (CHIME/FRB). The system was tested to ensure high detection efficiency, and the pulse calibration method was verified. Using an injection population of ~85,000 synthetic FRBs, we found that the correlation between fluence and SNR for injected FRBs was consistent with that of CHIME/FRB detections in the first CHIME/FRB catalog. We also noted that the sensitivity of the telescope varied strongly as a function of the broadened burst width, but not as a function of the dispersion measure. We conclude that some of the machine-learning based Radio Frequency Interference (RFI) mitigation methods used by CHIME/FRB can be re-trained using injection data to increase sensitivity to wide events, and that planned upgrades to the presented injection system will allow for determining a more accurate CHIME/FRB selection function in the near future.

Our recent re-analysis of the solar photospheric spectra with non-local thermodynamic equilibrium (non-LTE) models resulted in higher metal abundances compared to previous works. When applying the new chemical abundances to Standard Solar Model calculations, the new composition resolves the long-standing discrepancies with independent constraints on the solar structure from helioseismology. Critical to the determination of chemical abundances is the accuracy of the atomic data, specially the $f$-values, used in the radiative transfer models. Here we describe in detail the calculations of $f$-values for neutral oxygen and nitrogen used in our non-LTE models. Our calculations of $f$-values are based on a multi-method, multi-code approach and are the most detailed and extensive of its kind for the spectral lines of interest. We also report in this paper the details of extensive R-matrix calculation of photo-ionization cross sections for oxygen. Our calculation resulted in reliable $f$-values with well constrained uncertainties. We compare our results with previous theoretical and experimental determinations {of atomic data. We also quantify the influence of adopted photo-ionisation cross-sections on the spectroscopic estimate of the solar O abundance, using the data from different sources. We confirm that our 3D non-LTE value is robust and unaffected by the choice of photo-ionisation data, contrary to the recent claim made by Nahar.

The archetypical flare star UV Cet was observed by MeerKAT on 5-6 October 2021. A large radio outburst with a duration of $\sim\!2$ hr was observed between 886-1682 MHz with a time resolution of 8s and a frequency resolution of 0.84 MHz, enabling sensitive dynamic spectra to be formed. The emission is characterized by three peaks containing a multitude of broadband arcs or partial arcs in the time-frequency domain. In general, the arcs are highly right-hand circularly polarized. During end of the third peak, brief bursts occur that are significantly elliptically polarized. We present a simple model that appears to be broadly consistent with the characteristics of the radio emission from UV Cet. Briefly, the stellar magnetic field is modeled as a dipole aligned with the rotational axis of the star. The radio emission mechanism is assumed to be due to the cyclotron maser instability where x-mode radiation near the electron gyrofrequency is amplified. While the elliptically polarized bursts may be intrinsic to the source, rather stringent limits are imposed on the plasma density in the source and along the propagation path. We suggest that the elliptically polarized radiation may instead be the result of reflection on an over-dense plasma structure at some distance from the source. Radio emission from UV~Cet shares both stellar and planetary attributes.

The velocity profile of galaxies around voids is a key ingredient for redshift space distortion (RSD) measurements made using the void-galaxy correlation function. In this paper we use simulations to test whether the velocity profile of the tracers used to find the voids matches the velocity profile of the dark matter around these voids. A mismatch is expected and found in the inner part of voids, where tracers are very sparse. We discuss how this difference is caused by a selection effect where the void centre positions are correlated to the particular realization of the sparse tracers and their spatial distribution. In turn, this then affects the RSD void-galaxy correlation analysis. We show this by evaluating the Jacobian of the real to redshift space mapping using the tracer or matter velocity profile. Differences of the order of 20\% in the velocity profile translate into differences of the order of few percent in the Jacobian. This small discrepancy propagates to the monopole and quadrupole of the void-tracer correlation function, producing modifications of comparable magnitude to those from changes in $f\sigma_8$ at the level of the statistical uncertainties from current analyses.

M dwarfs are low-mass main-sequence stars, the most numerous type of stars in the solar neighbourhood, which are known to have significant magnetic activity. The aim of this work is to explore the dynamo solutions and magnetic fields of fully convective M dwarfs with varying magnetic Prandtl numbers ${\rm Pr_M}$, and a rotation period (${\rm P_{rot}}$) of 43 days. ${\rm Pr_M}$ is known to play an important role in the dynamo action; dynamos for low-$\rm Pr_M$ and large-$\rm Pr_M$ have very different properties. We performed three-dimensional magnetohydrodynamical (MHD) numerical simulations with the ``star-in-a-box'' model using stellar parameters for an M5 dwarf with $0.21 M_{\odot}$. We found that the dynamo solutions are sensitive to ${\rm Pr_M}$. The simulations at this rotation period present periodic cycles of the large-scale magnetic field up to ${\rm Pr_M} \leq 2$; for higher values the cycles disappear and irregular solutions to arise. Our results are consistent with previous studies and suggest that the dynamos operating in fully convective stars behave similarly as those in partially convective stars.

Off-axis collisions between galaxy clusters may induce the phenomenon of sloshing, causing dense gas to be dragged from the cool core of a cluster, resulting in a spiral of enhanced X-ray emission. Abell 2199 displays signatures of sloshing in its core and it is possible that the orbital plane of the collision is seen nearly edge-on. We aim to evaluate whether the features of Abell 2199 can be explained by a sloshing spiral seen under a large inclination angle. To address this, we perform tailored hydrodynamical $N$-body simulations of a non-frontal collision with a galaxy group of $M_{200}=1.6\times10^{13}\,{\rm M_{\odot}}$. We obtain a suitable scenario in which the group passed by the main cluster core 0.8 Gyr ago, with a pericentric separation of 292 kpc. Good agreement is obtained from the temperature maps as well as the residuals from a $\beta$-model fit to the simulated X-ray emission. We find that under an inclination of $i=70^{\circ}$ the simulation results remain consistent with the observations.

We present the first estimate of the Galactic nova rate based on optical transient surveys covering the entire sky. Using data from the All-Sky Automated Survey for Supernovae (ASAS-SN) and \textit{Gaia} -- the only two all-sky surveys to report classical nova candidates -- we find 39 confirmed Galactic novae and 7 additional unconfirmed candidates discovered from 2019--2021, yielding a nova discovery rate of $\approx 14$ yr$^{-1}$. Using accurate Galactic stellar mass models, three-dimensional dust maps, and incorporating realistic nova light curves, we have built a sophisticated Galactic nova model that allows an estimate of the recovery fraction of Galactic novae from these surveys over this time period. The observing capabilities of each survey are distinct: the high cadence of ASAS-SN makes it sensitive to fast novae, while the broad observing filter and high spatial resolution of \textit{Gaia} make it more sensitive to highly reddened novae across the entire Galactic plane and bulge. Despite these differences, we find that ASAS-SN and \textit{Gaia} give consistent Galactic nova rates, with a final joint nova rate of $26 \pm 5$ yr$^{-1}$. This inferred nova rate is substantially lower than found by many other recent studies. Critically assessing the systematic uncertainties in the Galactic nova rate, we argue that the role of faint fast-fading novae has likely been overestimated, but that subtle details in the operation of transient alert pipelines can have large, sometimes unappreciated effects on transient recovery efficiency. Our predicted nova rate can be directly tested with forthcoming red/near-infrared transient surveys in the southern hemisphere.

Combined action of helical motions of plasma (the $\alpha$ effect) and non-uniform (differential) rotation is a key dynamo mechanism of solar and galactic large-scale magnetic fields. Dynamics of magnetic helicity of small-scale fields is a crucial mechanism in a nonlinear dynamo saturation where turbulent magnetic helicity fluxes allow to avoid catastrophic quenching of the $\alpha$ effect. The convective zone of the Sun and solar-like stars as well as galactic discs are the source for production of turbulent magnetic helicity fluxes. In the framework of the mean-field approach and the spectral $\tau$ approximation, we derive turbulent magnetic helicity fluxes using the Coulomb gauge in a density-stratified turbulence. The turbulent magnetic helicity fluxes include non-gradient and gradient contributions. The non-gradient magnetic helicity flux is proportional to a nonlinear effective velocity (which vanishes in the absence of the density stratification) multiplied by small-scale magnetic helicity, while the gradient contributions describe turbulent magnetic diffusion of the small-scale magnetic helicity. In addition, the turbulent magnetic helicity fluxes contain source terms proportional to the kinetic $\alpha$ effect or its gradients, and also contributions caused by the large-scale shear (solar differential rotation). We have demonstrated that the turbulent magnetic helicity fluxes due to the kinetic $\alpha$ effect and its radial derivative in combination with the nonlinear magnetic diffusion of the small-scale magnetic helicity are dominant in the solar convective zone.

We present the results of our new observations of the young binary ZZ Tau with a circumbinary disc. The system was found to consist of two coeval (age $<2$ Myr) classical T Tauri stars with the total mass $0.86 \pm 0.09$ M$_\odot$, orbital period $46.8 \pm 0.8$ yr, semimajor axis $88.2 \pm 2.1$ mas, eccentricity $0.58 \pm 0.02$ and the orbital inclination $123.^{\rm o} 8 \pm 1.^{\rm o} 0.$ The accretion rate of ZZ Tau A and ZZ Tau B are approximately $7\times 10^{-10}$ and $2\times 10^{-10}$ M$_\odot$ yr$^{-1},$ respectively. No correlation was found between the long-term photometric variability of ZZ Tau and orbital position of its components. The periodic light variations with $P=4.171 \pm 0.002$ days was observed in the $BVRI$ bands presumably connected with an accretion (hot) spot on the surface of the primary (ZZ Tau A). At the same time no periodicity was observed in the $U$ band nor in the emission line profile variations probably due to the significant contribution of ZZ Tau B's emission, which dominates shortward of $\lambda \approx 0.4\,\mu$m. We argue that the extinction in the direction to the primary is noticeably larger than that to the secondary. It appeared that the rotation axis of the primary is inclined to the line of sight by $\approx 31^{\rm o} \pm 4^{\rm o}.$ We concluded also that ZZ Tau is the source of an CO molecular outflow, however, ZZ Tau IRS rather than ZZ Tau is the source of the Herbig-Haro object HH393.

It is well known that non-trivial squeezed tensor bispectra can lead to anisotropies in the inflationary stochastic gravitational wave (GW) background, providing us with an alternative and complementary window to primordial non-Gaussianities (NGs) with respect to the CMB. Previous works have highlighted the detection prospects of parity-even tensor NGs via the GW $I$-mode anisotropies. In this work we extend this by analysing for the first time the additional information carried by GW $V$-mode anisotropies due to squeezed NGs. We show that GW $V$ modes allow us to probe parity-odd squeezed $\langle \rm tts \rangle$ and $\langle \rm ttt \rangle$ bispectra. These bispectra break parity at the non-linear level and can be introduced by allowing alternative symmetry breaking patterns during inflation, like those comprised in solid inflation. Considering a BBO-like experiment, we find that a non-zero detection of squeezed $\langle \rm tts \rangle$ parity-odd bispectra in the $V$ modes dipole is possible without requiring any short-scale enhancement of the GW power spectrum amplitude over the constraints set by the CMB. We also briefly discuss the role of $V$-CMB cross-correlations. Our work can be extended in several directions and motivates a systematic search for polarized GW anisotropies in the next generations of GW experiments.

Cosmological models beyond $\Lambda$CDM, like those featuring massive neutrinos or modifications of gravity, often display a characteristic change (scale-dependent suppression or enhancement) in the matter power spectrum when compared to a $\Lambda$CDM baseline. It is therefore a widely held view that constraints on those models can be obtained by searching for such features in the clustering statistics of large-scale structure. However, when using biased tracers of matter in the analysis, the situation is complicated by the fact that the bias also depends on cosmology. Here we investigate how the selection of tracers affects the observed signatures for two examples of beyond-$\Lambda$CDM cosmologies: massive neutrinos and clustering dark energy ($k$-essence). We study the signatures in the monopole, quadrupole, and hexadecapole of the redshift-space power spectra for halo catalogues from large $N$-body simulations and argue that a fixed selection criterion based on local attributes like tracer mass leads to a near loss of signal in most cases. Instead, the full signal is recovered only if the selection of tracers is done at fixed bias. This emphasises the need to model or measure the bias parameters accurately in order to get meaningful constraints on the cosmological model.

We present LIMFAST, a semi-numerical code for computing the progress of reionization and line intensity mapping signals self-consistently, over large cosmological volumes and in short computational times. LIMFAST builds upon and extends the 21cmFAST code by implementing modern galaxy formation and evolution models. Furthermore, LIMFAST makes use of precomputed stellar population synthesis and photoionization results to obtain ensemble ionizing and line emission fields on large scales that vary with redshift, following the evolution of galaxy properties. We show LIMFAST calculations for the redshift evolution of the cosmic star formation rate, hydrogen neutral fraction, and metallicity in galaxies during reionization, which agree with current observational constraints. We also display the average signal with redshift, as well as the auto-power spectra at various redshifts, for the 21 cm line, the Ly$\alpha$ intergalactic and background emission, and the Ly$\alpha$, H$\alpha$, H$\beta$, [OII] $3727$\r{A}, and [OIII] $5007$\r{A} line emission from star formation. Overall, the LIMFAST results agree with calculations from other intensity mapping models, especially with those that account for the contribution of small halos during reionization. We further discuss the impact of considering redshift-space distortions, the use of local luminosity and star formation relations, and the dependence of line emission on the ionization parameter value. LIMFAST aims at being a resourceful tool for a broad range of intensity mapping studies, enabling the exploration of a variety of galaxy evolution and reionization scenarios and frequencies over large volumes in a short time scale.

It is widely believed that magnetic flux ropes are the key structure of solar eruptions; however, their observable counterparts are not clear yet. We study a flare associated with flux rope eruption in a comprehensive radiative magnetohydrodynamic simulation of flare-productive active regions, especially focusing on the thermodynamic properties of the plasma involved in the eruption and their relation to the magnetic flux rope. The pre-existing flux rope, which carries cold and dense plasma, rises quasi-statically before the eruption onsets. During this stage, the flux rope does not show obvious signatures in extreme ultraviolet (EUV) emission. After the flare onset, a thin `current shell' is generated around the erupting flux rope. Moreover, a current sheet is formed under the flux rope, where two groups of magnetic arcades reconnect and create a group of post-flare loops. The plasma within the `current shell', current sheet, and post-flare loops are heated to more than 10 MK. The post-flare loops give rise to abundant soft X-ray emission. Meanwhile a majority of the plasma hosted in the flux rope is heated to around 1 MK, and the main body of the flux rope is manifested as a bright arch in cooler EUV passbands such as AIA 171 \AA~channel.

We present 3D radiation-hydrodynamical (RHD) simulations of star cluster formation and evolution in massive, self-gravitating clouds, whose dust columns are optically thick to infrared (IR) photons. We use \texttt{VETTAM} -- a recently developed, novel RHD algorithm, which uses the Variable Eddington Tensor (VET) closure -- to model the IR radiation transport through the cloud. We also use realistic temperature ($T$) dependent IR opacities ($\kappa$) in our simulations, improving upon earlier works in this area, which used either constant IR opacities or simplified power laws ($\kappa \propto T^2$). We investigate the impact of the radiation pressure of these IR photons on the star formation efficiency (SFE) of the cloud, and its potential to drive dusty winds. We find that IR radiation pressure is unable to regulate star formation or prevent accretion onto the star clusters, even for very high gas surface densities ($\Sigma > 10^5 M_{\odot} \, \mathrm{pc}^{-2}$), contrary to recent semi-analytic predictions and simulation results using simplified treatments of the dust opacity. We find that the commonly adopted simplifications of $\kappa \propto T^2$ or constant $\kappa$ for the IR dust opacities leads to this discrepancy, as those approximations overestimate the radiation force. By contrast, with realistic opacities that take into account the micro-physics of the dust, we find that the impact of IR radiation pressure on star formation is very mild, even at significantly high dust-to-gas ratios ($\sim 3$ times solar), suggesting that it is unlikely to be an important feedback mechanism in controlling star formation in the ISM, and is also likely ineffective at launching galactic winds.

We investigate the cosmological constraints on the variable Chaplygin gas model from the latest observational data: SCP Union 2.1 compilation dataset of Type Ia supernovae (SNe Ia), Pantheon sample of SNe Ia and GWTC-3 of gravitational wave merger events. The variable Chaplygin gas is a model of interacting dark matter and dark energy which interpolates from dust-dominated era to quintessence dominated era. The variable Chaplygin gas model is shown to be compatible with Type Ia Supernovae and gravitational merger data. We have obtained tighter constraints on cosmological parameters $\Omega_m$ and $n$, using the Pantheon sample. By using the Markov chain Monte Carlo (MCMC) method on the Pantheon sample, we obtain $\Omega_m$=0.104 $\pm$ 0.027, n=0.435 $\pm$ 0.232 and $H_0$=70.453 $\pm$ 0.331 and on GWTC-3, we obtain $\Omega_m$=0.130 $\pm$ 0.091, n=0.815 $\pm$ 0.725 and $H_0$=69.359 $\pm$ 1.753.

This paper is devoted to a simple nonlocal de Sitter gravity model and its exact vacuum cosmological solutions. In the Einstein-Hilbert action with $\Lambda$ term, we introduce nonlocality by the following way: $R - 2 \Lambda = \sqrt{R-2\Lambda}\ \sqrt{R-2\Lambda} \to \sqrt{R-2\Lambda}\ F(\Box)\ \sqrt{R-2\Lambda} ,$ where ${F} (\Box) = 1 + \sum_{n= 1}^{+\infty} \big( f_n \Box^n + f_{-n} \Box^{-n} \big) $ is an analytic function of the d'Alembert-Beltrami operator $\Box$ and its inverse $\Box^{-1}$. By this way, $R$ and $\Lambda$ enter with the same form into nonlocal version as they are in the local one, and nonlocal operator $F(\Box)$ is dimensionless. The corresponding equations of motion for gravitational field $g_{\mu\nu}$ are presented. The first step in finding some exact cosmological solutions is solving the equation $\Box \sqrt{R-2\Lambda} = q \sqrt{R-2\Lambda} , $ where $ q =\zeta \Lambda \quad (\zeta \in \mathbb{R})$ is an eigenvalue and $\sqrt{R-2\Lambda}$ is an eigenfunction of the operator $\Box .$ We presented and discussed several exact cosmological solutions for homogeneous and isotropic universe. One of these solutions mimics effects that are usually assigned to dark matter and dark energy. Some other solutions are examples of the nonsingular bounce ones in flat, closed and open universe. There are also singular and cyclic solutions. All these cosmological solutions are a result of nonlocality and do not exist in the local de Sitter case.

Recent work has suggested that an additional $\lesssim 6.9 \rm{\, eV}$ per baryon of heating in the intergalactic medium is needed to reconcile hydrodynamical simulations with Lyman-$\alpha$ forest absorption line widths at redshift $z\simeq 0.1$. Resonant conversion of dark photon dark matter into low frequency photons is a viable source of such heating. We perform the first hydrodynamical simulations including dark photon heating and show that dark photons with mass $m_{A'}\sim 8\times 10^{-14}\rm\,eV\,c^{-2}$ and kinetic mixing $\epsilon \sim 5\times 10^{-15}$ can alleviate the heating excess. A prediction of this model is a non-standard thermal history for underdense gas at $z \gtrsim 3$.

We describe and study an instantaneous definition of eccentricity to be applied at the initial moment of full numerical simulations of binary black holes. The method consists of evaluating the eccentricity at the moment of maximum separation of the binary. We estimate it using up to third post-Newtonian (3PN) order, and compare these results with those of evolving (conservative) 3PN equations of motion for a full orbit and compute the eccentricity $e_r$ from the radial turning points, finding excellent agreement. We next include terms with spins up to 3.5PN, and then compare this method with the corresponding estimates of the eccentricity $e_r^{NR}$ during full numerical evolutions of spinning binary black holes, characterized invariantly by a fractional factor $0\leq f\leq1$ of the initial tangential momenta. It is found that our initial instantaneous definition is a very useful tool to predict and characterize even highly eccentric full numerical simulations.

The oscillating light axion field is known as wave dark matter. We propose an LC-resonance enhanced detection of the narrow band electric signals induced by the axion dark matter using a solenoid magnet facility. We provide full 3D electromagnetic simulation results for the signal electric field. The electric signal is enhanced by the high $Q$-factor of a resonant LC circuit and then amplified and detected by the state-of-the-art cryogenic electrical transport measurement technique. The amplifier noise is the leading noise in the setup. We demonstrate that the setup has promising sensitivity for axionic dark matter with mass $m_a$ below $10^{-6}$ eV. The projected sensitivities increase with the size of the magnetic field, and the electric signal measurement can be potentially sensitive to the QCD axion with $g_{a\gamma} \sim 10^{-16}$ GeV$^{-1}$ with a multi-meter scale magnetized region.

The existence of prolonged radiation domination prior to the Big Bang Nucleosynthesis (BBN), starting just after the inflationary epoch, is not yet established unanimously. If instead, the universe undergoes a non-standard cosmological phase, it will alter the Hubble expansion rate significantly and may also generate substantial entropy through non-adiabatic evolution. This leads to a thumping impact on the properties of relic species decoupled from the thermal bath before the revival of the standard radiation domination in the vicinity of the BBN. In this work, considering the Dirac nature of neutrinos, we have studied decoupling of ultra-relativistic right-handed neutrinos ($\nu_R$s) in presence of two possible non-standard cosmological phases. While in both cases we have modified Hubble parameters causing faster expansions in the early universe, one of the situations predicts a non-adiabatic evolution and thereby a slower redshift of the photon temperature due to the expansion. Considering the most general form of the collision term with Fermi-Dirac distribution and Pauli blocking factors, we have solved the Boltzmann equation numerically to obtain $\Delta{\rm N}_{\rm eff}$ for the three right-handed neutrinos. We have found that for a large portion of parameter space, the combined effect of early decoupling of $\nu_R$ as well as the slower redshift of photon bath can easily hide the signature of right-handed neutrinos, in spite of precise measurement of $\Delta{\rm N}_{\rm eff}$, at the next generation CMB experiments like CMB-S4, SPT-3G etc. This however will not be applicable for the scenarios with only fast expansion.

In this present work, we study the observational appearance of Kerr-Melvin black hole (KMBH) illuminated by an accretion disk. The accretion disk is assumed to be located on the equatorial plane and be thin both geometrically and optically. Considering the fact that outside the innermost stable circular orbit (ISCO) the accretion flow moves in prograde or retrograde circular orbit and falls towards the horizon along plunging orbit inside the ISCO, we develop the numerical backward ray-tracing method and obtain the images of KMBH accompanying with the accretion disk for various black hole spins, strengths of magnetic fields and inclination angles of observers. We present the intensity distribution horizontally and longitudinally and show the profiles of the red-shift for the direct and lensed images. Our study suggests that the inner shadow and critical curves can be used to estimate the magnetic field around a black hole without degeneration.

Gravitational waves from binary black hole mergers have allowed us to directly observe stellar-mass black hole binaries for the first time, and therefore explore their formation channels. One of the ways to infer how a binary system is assembled is by measuring the system's orbital eccentricity. Current methods of parameter estimation do not include all physical effects of eccentric systems such as spin-induced precession, higher-order modes and the initial argument of periapsis: an angle describing the orientation of the orbital ellipse. We explore how varying the argument of periapsis changes gravitational waveforms and study its effect on the inference of astrophysical parameters. We use the eccentric spin-aligned waveforms \texttt{TEOBResumS} and \texttt{SEOBNRE} to measure the change in the waveforms as the argument of periapsis is changed. We find that the argument of periapsis is likely to be resolvable in the foreseeable future only for the loudest events observed by LIGO--Virgo--KAGRA. The systematic error in previous analyses that have not taken into account the argument of periapsis is likely to be small.

We compute the red and blue shifts for astrophysical and cosmological sources. In particular, we consider low, intermediate and high gravitational energy domains. Thereby, we handle the binary system Earth - Mars as low energy landscape whereas white dwarfs and neutron stars as higher energy sources. To this end, we take into account a spherical Schwarzschild - de Sitter spacetime and an axially symmetric Zipoy - Voorhees metric to model all the aforementioned systems. Feasible outcomes come from modelling neutron stars and white dwarfs with the Zipoy - Voorhees metric, where quadrupole effects are relevant, and framing solar system objects using a Schwarzschild - de Sitter spacetime. In the first case, large $\delta$ parameters seem to be favorite, leading to acceptable bounds mainly for neutron stars. In the second case, we demonstrate incompatible red and blue shifts with respect to lunar and satellite laser ranging expectations, once the cosmological constant is taken to Planck satellite's best fit. To heal this issue, we suggest coarse-grained experimental setups and propose Phobos for working out satellite laser ranging in order to get more suitable red and blue shift intervals, possibly more compatible than current experimental bounds. Implications to cosmological tensions are also debated.

We construct a model of quintessential inflation in Palatini $R^2$ gravity employing a scalar field with a simple exponential potential and coupled to gravity with a running non-minimal coupling. At early times, the field acts as the inflaton, while later on it becomes the current dark energy. Combining the scalar sector with an ideal fluid, we study the cosmological evolution of the model from inflation all the way to dark energy domination. We interpret the results in the Einstein frame, where a coupling emerges between the fluid and the field, feeding energy from the former to the latter during the matter-dominated era. We perform a numerical scan over the parameter space and find points that align with observations for both the inflationary CMB data and the late-time behaviour. The final dark energy density emerges from an interplay between the model parameters, without requiring the extreme fine-tuning of the cosmological constant in $\Lambda$CDM.

We explore a mechanism of geometric cancellation of vacuum energy zero-point fluctuations, resolving \emph{de facto} the cosmological constant problem. To do so, we assume at primordial times that vacuum energy fuels an inflationary quadratic hilltop potential nonminimally coupled to gravity through a standard Yukawa-like interacting term, whose background lies on a perturbed Friedmann-Robertson-Walker metric. We demonstrate how vacuum energy release transforms into geometric particles, adopting a quasi-de Sitter phase where we compute the expected particle density and mass ranges. Hence, we propose the most suitable dark matter candidates, showing under which circumstances we can interpret dark matter as constituted by geometric quasiparticles. We confront our predictions with quantum particle production and constraints made using a Higgs portal. In addition, the role of the bare cosmological constant is reinterpreted to speed up the universe today. Thus, consequences on the standard $\Lambda$CDM paradigm are critically highlighted, showing how both coincidence and fine-tuning issues can be healed requiring the Israel-Darmois matching conditions between our involved inhomogeneous and homogeneous phases.

We study the Jesuit mission churches in America, which for almost two centuries were the most representative constructions in the process of Christian evangelisation on the continent until the expulsion of the Order in 1767. The main objective is to discern possible orientation patterns in the structures studied and to evaluate whether these orientations are related to the rising of the sun or other celestial bodies on the local horizon, as suggested by the texts of early Christian writers, which could provide important and novel information about their history and construction. Archaeoastronomical fieldwork has been carried out in the past in two large regions of South America: the historic provinces of Paraquaria and Chiquitan\'ia (eastern present-day Bolivia). These data must now be interpreted within a broader cultural and geographical context. Thus, in order to obtain a more complete picture of religious architecture on the continent, it is necessary to analyse and compare the results with the orientation of Jesuit churches built in the 17th and 18th centuries in North America (viceroyalty of New Spain). We approached this project through the historical and cultural analysis of the former mission sites and by means of satellite imagery. In this paper we present some preliminary results we have reached.