Papers for Wednesday, Dec 01 2021

The advent of space-based photometry missions such as CoRoT, Kepler and TESS has sparkled the development of asteroseismology and exoplanetology. The advent of PLATO will further strengthen such multi-disciplinary studies. Testing asteroseismic modelling and its importance for our understanding of planetary systems is crucial. We carried out a detailed modelling of Kepler-93, an exoplanet host star observed by Kepler. This star is particularly interesting as it is very similar to the PLATO benchmark target (G spectral type, ~ 6000K, ~ 1 Msun and ~ 1 Rsun) and provides a real-life testbed for potential procedures to be used for PLATO. We use global and local minimization techniques for the seismic modelling of Kepler-93, varying the ingredients of our stellar models. We compute seismic inversions of the mean density. We use these revised stellar parameters to provide new planetary parameters and simulate the orbital evolution of the system under the effects of tides and atmospheric evaporation. Our fundamental parameters for Kepler-93: mean density = 1.654 +/- 0.004 g/cm3, M = 0.907 +/- 0.023 Msun , R = 0.918 +/- 0.008 Rsun and Age = 6.78 +/- 0.32 Gyr. The uncertainties we report for this benchmark are within the requirements of PLATO. For the exoplanet Kepler-93b, we find Mp = 4.01 +/- 0.67 Mearth, Rp = 1.478 +/- 0.014 Rearth and semi-major axis a = 0.0533 +/- 0.0005 AU. According to our simulations, it seems unlikely that Kepler-93b formed with a mass large enough to be impacted by stellar tides. For the benchmark of PLATO, detailed asteroseismic modelling procedures will be able to provide fundamental stellar parameters within the requirements. We illustrate what synergies can be achieved regarding the orbital evolution and atmospheric evaporation of exoplanets. We note the importance of the high-quality radial velocity follow-up to constrain the formation scenarii of exoplanets.

While gravitational microlensing by planetary systems can provide unique vistas on the properties of exoplanets, observations of such 2-body microlensing events can often be explained with multiple and distinct physical configurations, so-called model degeneracies. An understanding of the intrinsic and exogenous origins of different classes of degeneracy provides a foundation for phenomenological interpretation. Here, leveraging a fast machine-learning based inference framework, we present the discovery of a new regime of degeneracy--the offset degeneracy--which unifies the previously known close-wide and inner-outer degeneracies, generalises to resonant caustics, and upon reanalysis, is ubiquitous in previously published planetary events with 2-fold degenerate solutions. Importantly, our discovery suggests that the commonly reported close-wide degeneracy essentially never arises in actual events and should, instead, be more suitably viewed as a transition point of the offset degeneracy. While previous studies of microlensing degeneracies are largely studies of degenerate caustics, our discovery demonstrates that degenerate caustics do not necessarily result in degenerate events, which for the latter it is more relevant to study magnifications at the location of the source. This discovery fundamentally changes the way in which degeneracies in planetary microlensing events should be interpreted, suggests a deeper symmetry in the mathematics of 2-body lenses than has previously been recognised, and will increasingly manifest itself in data from new generations of microlensing surveys.

Short-period double white dwarf (DWD) binaries will be the most prolific source of gravitational waves (GWs) for the Laser Interferometer Space Antenna (LISA). Not only will tens of thousands of DWDs be individually resolved, but DWDs with GW frequencies below ~1 mHz will be the dominant contributor to a stochastic foreground caused by confusion from overlapping GW signals, limiting the detectability of individual sources of all kinds. Population modeling of Galactic DWDs typically assumes a standard binary fraction of 50%. However, recent observations have shown that the binary fraction of close (P$_{\rm{orb}} \leq 10^4$ days) solar-type stars exhibits a strong anti-correlation with metallicity. In this study we perform the first simulation of the Galactic DWD population observable by LISA which incorporates an empirically-derived metallicity-dependent binary fraction. We simulate DWDs using the binary population synthesis suite COSMIC and incorporate a metallicity-dependent star formation history to create a Galactic population of short-period DWDs. We compare two models: one which assumes a metallicity-dependent binary fraction, and one with a binary fraction of 50%. We find that while metallicity impacts the evolution and intrinsic properties of our simulated DWD progenitor binaries, the LISA-resolvable populations of the two models remain roughly indistinguishable. However, the size of the total Galactic DWD population orbiting in the LISA frequency band is reduced by more than half when accounting for a metallicity-dependent binary fraction. This effect serves to lower the confusion foreground, effectively increasing the sensitivity for detecting all types of low-frequency LISA sources. We repeat our analysis for three different assumptions for Roche-lobe overflow interactions and find the population reduction to be robust when a metallicity-dependent binary fraction is assumed.

Cosmic gas makes up about 90% of baryonic matter in the Universe and H$_2$ is the closest molecule to star formation. In this work we study cold neutral gas and its H$_2$ component at different epochs, exploiting state-of-the-art hydrodynamic simulations that include time-dependent atomic and molecular non-equilibrium chemistry coupled to star formation, feedback effects, different UV backgrounds presented in the recent literature and a number of additional processes - such as gas self-shielding, H$_2$ dust grain catalysis, photoelectric and cosmic-ray heating - occurring during structure formation (ColdSIM). We find neutral-gas mass density parameters $ \Omega_{\rm neutral} \simeq $10$^{-3}$ and increasing from lower to higher redshift, in agreement with available HI data. Resulting H$_2$ fractions can be as high as $\sim $50% at $z\sim $4-8, in line with the latest high-$z$ measurements. Albeit dependent on the adopted UV background, derived $ \Omega_{\rm H_2} $ values agree with observations up to $z\sim$7 and both HI and H$_2$ trends are better reproduced by our non-equilibrium H$_2$-based star formation modelling. The predicted gas depletion timescales decrease towards lower $z$, with H$_2$ depletion times remaining below the Hubble time and comparable to the dynamical time at all considered redshifts. This implies that non-equilibrium molecular cooling is efficient at driving cold-gas collapse in a broad variety of environments and since the very early cosmic epochs. In appendix, we show detailed analyses of individual processes, as well as simple numerical parameterizations and fits to account for them. Our findings suggest that, in addition to HI, non-equilibrium H$_2$ observations are pivotal probes for assessing cold-gas abundances and the role of UV background radiation - Abridged

The power spectrum of the nonlinearly evolved large-scale mass distribution recovers only a minority of the information available on the mass fluctuation amplitude. We investigate the recovery of this information in 2D "slabs" of the mass distribution averaged over $\approx100$~$h^{-1}$Mpc along the line of sight, as might be obtained from photometric redshift surveys. We demonstrate a Hamiltonian Monte Carlo (HMC) method to reconstruct the non-Gaussian mass distribution in slabs, under the assumption that the projected field is a point-transformed Gaussian random field, Poisson-sampled by galaxies. When applied to the \textit{Quijote} $N$-body suite at $z=1$ and at a transverse resolution of 2~$h^{-1}$Mpc, the method recovers $\sim 30$ times more information than the 2D power spectrum in the well-sampled limit, recovering the Gaussian limit on information. At a more realistic galaxy sampling density of $0.01$~$h^3$Mpc$^{-3}$, shot noise reduces the information gain to a factor of five improvement over the power spectrum at resolutions of 4~$h^{-1}$Mpc or smaller.

Future searches for gravitational waves from space will be sensitive to double compact objects (DCOs) in our Milky Way. We present new simulations of the populations of double black holes (BHBHs), black hole neutron stars (BHNSs) and double neutron stars (NSNSs) that will be detectable by the planned space-based gravitational wave detector LISA. For our estimates, we use an empirically-informed model of the metallicity dependent star formation history of the Milky Way. We populate it using an extensive suite of binary population-synthesis predictions for varying assumptions relating to mass transfer, common-envelope, supernova kicks, remnant masses and wind mass loss physics. For a 4(10)-year LISA mission, we predict between 30-370(50-550) detections over these variations, out of which 6-154(9-238) are BHBHs, 2-198(3-289) are BHNSs and 3-35(4-57) are NSNSs. We discuss how the variations in the physics assumptions alter the distribution of properties of the detectable systems, even when the detection rates are unchanged. In particular we discuss the observable characteristics such as the chirp mass, eccentricity and sky localisation and how the BHBH, BHNS and NSNS populations can be distinguished, both from each other and from the more numerous double white dwarf population. We further discuss the possibility of multi-messenger observations of pulsar populations with the Square Kilometre Array (SKA) and assess the benefits of extending the LISA mission.

We report on the masses ($M_\mathrm{WD}$), effective temperatures ($T_\mathrm{eff}$) and secular mean accretion rates ($\langle \dot{M} \rangle$) of 43 cataclysmic variable (CV) white dwarfs, 42 of which were obtained from the combined analysis of their $\mathit{Hubble~Space~Telescope}$ ultraviolet data with the parallaxes provided by the Early Third Data Release of the $\mathit{Gaia}$ space mission, and one from the white dwarf gravitational redshift. Our results double the number of CV white dwarfs with an accurate mass measurement, bringing the total census to 89 systems. From the study of the mass distribution, we derive $\langle M_\mathrm{WD} \rangle = 0.81^{+0.16}_{-0.20}\,\mathrm{M_\odot}$, in perfect agreement with previous results, and find no evidence of any evolution of the mass with orbital period. Moreover, we identify five systems with $M_\mathrm{WD} < 0.5\mathrm{M_\odot}$, which are most likely representative of helium-core white dwarfs, showing that these CVs are present in the overall population. We reveal the presence of an anti-correlation between the average accretion rates and the white dwarf masses for the systems below the $2-3\,$h period gap. Since $\langle \dot{M} \rangle$ reflects the rate of system angular momentum loss, this correlation suggests the presence of an additional mechanism of angular momentum loss that is more efficient at low white dwarf masses. This is the fundamental concept of the recently proposed empirical prescription of consequential angular momentum loss (eCAML) and our results provide observational support for it, although we also highlight how its current recipe needs to be refined to better reproduce the observed scatter in $T_\mathrm{eff}$ and $\langle \dot{M} \rangle$, and the presence of helium-core white dwarfs.

We have applied our method for weighing the Galactic disk using phase-space spirals to the Gaia EDR3 proper motion sample. For stars in distant regions of the Galactic disk, the latitudinal proper motion has a close projection with vertical velocity, such that the phase-space spiral in the plane of vertical position and vertical velocity can be observed without requiring that all stars have available radial velocity information. We divided the Galactic plane into 360 separate data samples, each corresponding to an area cell in the Galactic plane in the distance range of 1.4-3.4 kpc, with an approximate cell length of 200-400 pc. Roughly half of our data samples were disqualified altogether due to severe selection effects, especially in the direction of the Galactic centre. In the remainder, we were able to infer the vertical gravitational potential by fitting an analytic model of the phase-space spiral to the data. This work is the first of its kind, in the sense that we are weighing distant regions of the Galactic disk with a high spatial resolution, without relying on the strong assumptions of axisymmetry. Post-inference, we fit a thin disk scale length of $2.2\pm 0.1$ kpc, although this value is sensitive to the considered spatial region. We see surface density variations as a function of azimuth on the order of 10-20 %, which is roughly the size of our estimated sum of potential systematic biases. With this work, we have demonstrated that our method can be used to weigh distant regions of the Galactic disk despite strong selection effects. We expect to reach even greater distances and improve our accuracy with future Gaia data releases and further improvements to our method.

Using the state-of-the-art suite of hydrodynamic simulations Simba, as well as its dark-matter-only counterpart, we study the impact of the presence of baryons and of different stellar/AGN feedback mechanisms on large-scale structure, halo density profiles, and on the abundance of different baryonic phases within halos and in the intergalactic medium (IGM). The unified picture that emerges from our analysis is that the main physical drivers shaping the distribution of matter at all scales are star formation-driven galactic outflows at $z>2$ for lower mass halos and AGN jets at $z<2$ in higher mass halos. Feedback suppresses the baryon mass function with time relative to the halo mass function, and it even impacts the halo mass function itself at the ~20% level, particularly evacuating the centres and enhancing dark matter just outside halos. At early epochs baryons pile up in the centres of halos, but by late epochs and particularly in massive systems gas has mostly been evacuated from within the inner halo. AGN jets are so efficient at such evacuation that at low redshifts the baryon fraction within $\sim 10^{12}-10^{13} \, \rm M_{\odot}$ halos is only 25% of the cosmic baryon fraction, mostly in stars. The baryon fraction enclosed in a sphere around such halos approaches the cosmic value $\Omega_{\rm b}/\Omega_{\rm m}$ only at 10-20 virial radii. As a result, 87% of the baryonic mass in the Universe lies in the IGM at $z=0$, with 67% being in the form of warm-hot IGM ($T>10^5 \, \rm K$).

We study bulge formation in MW/M31-like galaxies in a $\Lambda$-cold dark matter scenario, focusing on the origin of high- and low-Sersic index bulges. For this purpose we use TNG50, a simulation of the IllustrisTNG project that combines a resolution of $\sim 8 \times 10^4 M_{\odot}$ in stellar particles with a cosmological volume 52 cMpc in extent. We parametrize bulge surface brightness profiles by the S\'ersic index and the bulge-to-total (B/T) ratio obtained from two-component photometric decompositions. In our sample of 287 MW/M31-like simulated galaxies, $17.1\%$ of photometric bulges exhibit high-S\'ersic indices and $82.9\%$ show low-S\'ersic indices. We study the impact that the environment, mergers and bars have in shaping the surface brightness profiles. We explore two different definitions for local environment and find no correlation between bulge properties and the environment where they reside. Simulated galaxies with higher S\'ersic indices show, on average, a higher fraction of ex-situ stars in their kinematically selected bulges. For this bulge population the last significant merger (total mass ratio $m_{\rm sat}/m_{\rm host} > 0.1$) occurs, on average, at later times. However, a substantial fraction of low-S\'ersic index bulges also experience a late significant merger. We find that bars play an important role in the development of the different types of photometric bulges. We show that the fraction of simulated galaxies with strong bars is smaller for the high- than for the low-S\'ersic index population, reaching differences of $20\%$ at $z > 1$. Simulated galaxies with high fractions of ex-situ stars in the bulge do not develop strong bars. Conversely, simulated galaxies with long-lived strong bars have bulges with ex-situ fractions, $f_{\rm ex-situ} < 0.2$.

A widely-held assumption is that each single white dwarf containing observable rocky debris requires the presence of at least one terrestrial or giant planet to have gravitationally perturbed the progenitor of the debris into the star. However, these planets could have been previously been engulfed by the star or escaped the system, leaving behind asteroids, boulders, cobbles, pebbles, sand and dust. These remaining small bodies could then persist throughout the host star's evolution into a white dwarf at ~2-100 au scales, and then be radiatively dragged into the white dwarf without the help of a planet. Here we identify the parameter space and cooling ages for which this one metal-pollution mechanism is feasible by, for the first time, coupling Poynting-Robertson drag, the Yarkovsky effect and the YORP effect solely from rapidly dimming white dwarf radiation. We find that this no-planet pollution scenario is efficient for remnant 10^{-5}-10^{-4} m dust up to about 80 au, 10^{-4}-10^{-3} m sand up to about 25 au and 10^{-3}-10^{-2} m small pebbles up to about 8 au, and perhaps 10^{-1}-10^{0} m small boulders up to tens of au. Further, young white dwarf radiation can spin up large strength-less boulders with radii 10^{2}-10^{3} m to destruction, breaking them down into smaller fragments which then can be dragged towards the white dwarf. Our work hence introduces a planet-less metal-pollution mechanism that may be active in some fraction of white dwarf planetary systems.

The tidal tails of globular clusters have been shown to be sensitive to the external tidal field. We investigate how Galactic globular clusters with observed tails are affected by satellite dwarf galaxies by simulating tails in galaxy models with and without dwarf galaxies. The simulations indicate that tidal tails can be subdivided into into three categories based on how they are affected by dwarf galaxies: 1) dwarf galaxies perturb the progenitor cluster's orbit (NGC 4590, Pal 1, Pal 5), 2) dwarf galaxies perturb the progenitor cluster's orbit and individual tail stars (NGC 362, NGC 1851, NGC 4147, NGC 5466, NGC 7492, Pal 14, Pal 15), and 3) dwarf galaxies negligibly affect tidal tails (NGC 288, NGC 5139, NGC 5904, Eridanus). Perturbations to a cluster's orbit occur when dwarf galaxies pass within its orbit, altering the size and shape of the orbital and tail path. Direct interactions between one or more dwarf galaxies and tail stars lead to kinks and spurs, however we find that features are more difficult to observe in projection. We further find that the tails of Pal 5 are shorter in the galaxy model with dwarf galaxies as it is closer to apocentre, which results in the tails being compressed. Additional simulations reveal that differences between tidal tails in the two galaxy models are primarily due to the Large Magellanic Cloud. Understanding how dwarf galaxies affect tidal tails allows for tails to be used to map the distribution of matter in dwarf galaxies and the Milky Way.

The empirical upper limit to Red Supergiant (RSG) luminosity, known as the Humphreys-Davidson (HD) limit, has been commonly explained as being caused by the stripping of stellar envelopes by metallicity-dependent, line-driven winds. As such, the theoretical expectation is that the HD limit should be higher at lower metallicity, where weaker mass-loss rates mean that higher initial masses are required for an envelope to be stripped. In this paper, we test this prediction by measuring the luminosity function of RSGs in M31 and comparing to those in the LMC and SMC. We find that $\log (L_{\rm max}/L_{\odot}) = 5.53 \pm 0.03$ in M31 (Z $\gtrsim$ Z$_{\odot}$), consistent with the limit found for both the LMC (Z $\sim$ 0.5 Z$_{\odot}$) and SMC (Z $\sim$ 0.25 Z$_{\odot}$), while the RSG luminosity distributions in these 3 galaxies are consistent to within 1$\sigma$. We therefore find no evidence for a metallicity dependence on both the HD limit and the RSG luminosity function, and conclude that line-driven winds on the main sequence are not the cause of the HD limit.

We explore a new approach for extracting reionization-era contributions to the kinetic Sunyaev-Zel'dovich (kSZ) effect. Our method utilizes the cross-power spectrum between filtered and squared maps of the cosmic microwave background (CMB) and photometric galaxy surveys during the Epoch of Reionization (EoR). This kSZ$^2$-galaxy cross-power spectrum statistic has been successfully detected at lower redshifts ($z \lesssim 1.5$). Here we extend this method to $z \gtrsim 6$ as a potential means to extract signatures of patchy reionization. We model the expected signal across multiple photometric redshift bins using semi-numeric simulations of the reionization process. In principle, the cross-correlation statistic robustly extracts reionization-era contributions to the kSZ signal, while its redshift evolution yields valuable information regarding the timing of reionization. Specifically, the model cross-correlation signal near $\ell \sim 1,000$ peaks during the early stages of the EoR, when about 20% of the volume of the universe is ionized. Detectible $\ell$-modes mainly reflect squeezed triangle configurations of the related bispectrum, quantifying correlations between the galaxy overdensity field on large scales and the smaller-scale kSZ power. We forecast the prospects for detecting this signal using future wide-field samples of Lyman-break galaxies from the Roman Space Telescope and next-generation CMB surveys including the Simons Observatory, CMB-S4, and CMB-HD. We find that a roughly 13$\sigma$ detection is possible for CMB-HD and Roman after summing over all $\ell$-modes. We discuss the possibilities for improving this approach and related statistics, with the aim of moving beyond simple detections to measure the scale and redshift dependence of the cross-correlation signals.

Virial black hole mass ($M_{BH}$) determination directly involves knowing the broad line region (BLR) clouds velocity distribution, their distance from the central supermassive black hole ($R_{BLR}$) and the virial factor ($f$). Understanding whether biases arise in $M_{BH}$ estimation with increasing obscuration is possible only by studying a large (N$>$100) statistical sample of obscuration unbiased (hard) X-ray selected active galactic nuclei (AGN) in the rest-frame near-infrared (0.8-2.5$\mu$m) since it penetrates deeper into the BLR than the optical. We present a detailed analysis of 65 local BAT-selected Seyfert galaxies observed with Magellan/FIRE. Adding these to the near-infrared BAT AGN spectroscopic survey (BASS) database, we study a total of 314 unique near-infrared spectra. While the FWHMs of H$\alpha$ and near-infrared broad lines (He\textsc{i}, Pa$\beta$, Pa$\alpha$) remain unbiased to either BLR extinction or X-ray obscuration, the H$\alpha$ broad line luminosity is suppressed when $N_H\gtrsim10^{21}$ cm$^{-2}$, systematically underestimating $M_{BH}$ by $0.23-0.46$ dex. Near-infrared line luminosities should be preferred to H$\alpha$ until $N_H<10^{22}$ cm$^{-2}$, while at higher obscuration a less biased $R_{BLR}$ proxy should be adopted. We estimate $f$ for Seyfert 1 and 2 using two obscuration-unbiased $M_{BH}$ measurements, i.e. the stellar velocity dispersion and a BH mass prescription based on near-infrared and X-ray, and find that the virial factors do not depend on redshift or obscuration, but for some broad lines show a mild anti-correlation with $M_{BH}$. Our results show the critical impact obscuration can have on BLR characterization and the importance of the near-infrared and X-rays for a less biased view of the BLR.

Hydrogen emission lines can provide extensive information about star-forming galaxies in both the local and high-redshift Universe. We present a detailed Lyman continuum (LyC), Lyman-alpha (Ly{\alpha}), and Balmer line (H{\alpha} and H\b{eta}) radiative transfer study of a high-resolution isolated Milky-Way simulation using the Arepo-RT radiation hydrodynamics code with the SMUGGLE galaxy formation model. The realistic framework includes stellar feedback, non-equilibrium thermochemistry, and dust grain evolution in the interstellar medium (ISM). We extend our Cosmic Ly{\alpha} Transfer (COLT) code with photoionization equilibrium Monte Carlo radiative transfer for self-consistent end-to-end (non-)resonant line predictions. Accurate LyC reprocessing to recombination emission requires modelling pre-absorption by dust (27.5%), helium ionization (8.7%), and anisotropic escape fractions (7.9%), as these reduce the available budget for hydrogen line emission (55.9%). We investigate the role of the multiphase dusty ISM, disc geometry, gas kinematics, and star formation activity in governing the physics of emission and escape, focusing on the time variability, gas phase structure, and spatial, spectral, and viewing angle dependence of the emergent photons. Isolated disc simulations are well-suited for comprehensive observational comparisons with local H{\alpha} surveys, but would require a proper cosmological circumgalactic medium (CGM) environment as well as less dust absorption and rotational broadening to serve as analogs for high-redshift Ly{\alpha} emitting galaxies. Future applications of our framework to next-generation cosmological simulations of galaxy formation including radiation-hydrodynamics that resolve <10 pc multiphase ISM and <1 kpc CGM structures will provide crucial insights and predictions for current and upcoming Ly{\alpha} observations.

Stellar binaries represent a substantial fraction of stellar systems, especially among young stellar objects. Accordingly, binaries play an important role in setting the architecture of a large number of protoplanetary discs. Binaries in coplanar and polar orientations with respect to the circumbinary disc are stable configurations and could induce non-axisymmetric structures in the dust and gas distributions. In this work, we suggest that the structures shown in the central region of the protoplanetary disc HD 169142 are produced by the presence of an inner stellar binary and a circumbinary (P-type) planet. We find that a companion with a mass-ratio of 0.1, semi-major axis of 9.9 au, eccentricity of 0.2, and inclination of 90{\deg}, together with a 2 Jupiter Mass coplanar planet on a circular orbit at 45 au reproduce the structures at the innermost ring observed at 1.3 mm and the shape of spiral features in scattered light observations. The model predicts changes in the disc's dust structure, and star's astrometric parameters, which would allow testing its veracity by monitoring this system over the next 20 years.

We investigated the hydrogen isotopic compositions and water contents of pyroxenes in two recent ordinary chondrite falls, namely, Chelyabinsk (2013 fall) and Benenitra (2018 fall), and compared them to three ordinary chondrite Antarctic finds, namely Graves Nunataks GRA 06179, Larkman Nunatak LAR 12241, and Dominion Range DOM 10035. The pyroxene minerals in Benenitra and Chelyabinsk are hydrated ($\sim$0.018-0.087 wt.$\%$ H$_2$O) and show D-poor isotopic signatures ($\delta$D$_{SMOW}$ from -444$\unicode{x2030}$ to -49$\unicode{x2030}$). On the contrary, the ordinary chondrite finds exhibit evidence of terrestrial contamination with elevated water contents ($\sim$0.039-0.174 wt.$\%$) and values (from -199$\unicode{x2030}$ to -14$\unicode{x2030}$). We evaluated several small parent body processes that are likely to alter the measured compositions in Benenitra and Chelyabinsk, and inferred that water-loss in S-type planetesimals is minimal during thermal metamorphism. Benenitra and Chelyabinsk hydrogen compositions reflect a mixed component of D-poor nebular hydrogen and water from the D-rich mesostases. 45-95$\%$ of water in the minerals characterized by low $\delta$D$_{SMOW}$ values was contributed by nebular hydrogen. S-type asteroids dominantly composed of nominally anhydrous minerals can hold 254-518 ppm of water. Addition of a nebular water component to nominally dry inner Solar System bodies during accretion suggests a reduced need of volatile delivery to the terrestrial planets during late accretion.

Astronomical data is full of holes. While there are many reasons for this missing data, the data can be randomly missing, caused by things like data corruptions or unfavourable observing conditions. We test some simple data imputation methods(Mean, Median, Minimum, Maximum and k-Nearest Neighbours (kNN)), as well as two more complex methods (Multivariate Imputation by using Chained Equation (MICE) and Generative Adversarial Imputation Network (GAIN)) against data where increasing amounts are randomly set to missing. We then use the imputed datasets to estimate the redshift of the galaxies, using the kNN and Random Forest ML techniques. We find that the MICE algorithm provides the lowest Root Mean Square Error and consequently the lowest prediction error, with the GAIN algorithm the next best.

Partially ionized plasmas, such as the solar chromosphere, require a generalized Ohm's law including the effects of ambipolar and Hall drift. While both describe transport processes that arise from the multifluid equations and are therefore of hyperbolic nature, they are often incorporated in models as a diffusive, i.e. parabolic process. While the formulation as such is easy to include in standard MHD models, the resulting diffusive time-step constraints do require often a computationally more expensive implicit treatment or super-time-stepping approaches. In this paper we discuss an implementation that retains the hyperbolic nature and allows for an explicit integration with small computational overhead. In the case of ambipolar drift, this formulation arises naturally by simply retaining a time derivative of the drift velocity that is typically omitted. This alone leads to time-step constraints that are comparable to the native MHD time-step constraint for a solar setup including the region from photosphere to lower solar corona. We discuss an accelerated treatment that can further reduce time-step constraints if necessary. In the case of Hall drift we propose a hyperbolic formulation that is numerically similar to that for the ambipolar drift and we show that the combination of both can be applied to simulations of the solar chromosphere at minimal computational expense.

The periodic pulsations of stars teach us about their underlying physical process. We present a convolutional autoencoder-based pipeline as an automatic approach to search for out-of-distribution anomalous periodic variables within The Zwicky Transient Facility Catalog of Periodic Variable Stars (ZTF CPVS). We use an isolation forest to rank each periodic variable by its anomaly score. Our overall most anomalous events have a unique physical origin: they are mostly highly variable and irregular evolved stars. Multiwavelength data suggest that they are most likely Red Giant or Asymptotic Giant Branch stars concentrated in the Milky Way galactic disk. Furthermore, we show how the learned latent features can be used for the classification of periodic variables through a hierarchical random forest. This novel semi-supervised approach allows astronomers to identify the most anomalous events within a given physical class, significantly increasing the potential for scientific discovery.

We explore how the galaxy stellar spins acquire a peculiar tendency of being aligned with the major principal axes of the local tidal fields, in contrast to their DM counterparts which tend to be perpendicular to them, regardless of their masses. Analyzing the halo and subhalo catalogs from the IllustrisTNG 300 hydrodynamic simulations at $z\le 1$, we determine the cosines of the alignment angles, $\cos\alpha$, between the galaxy stellar and DM spins. Creating four $\cos\alpha$-selected samples of the galaxies and then controlling them to share the same density and mass distributions, we determine the average strengths of the alignments between the galaxy stellar spins and the tidal major axes over each sample. It is clearly shown that at $z\le 0.5$ the more severely the galaxy stellar spin directions deviate from the DM counterparts, the stronger the peculiar tidal alignments become. Taking the ensemble averages of such galaxy properties as central blackhole to stellar mass ratio, specific star formation rate, formation epoch, stellar-to-total mass ratio, velocity dispersions, average metallicity, and degree of the cosmic web anisotropy over each sample, we also find that all of these properties exhibit either strong correlations or anti-correlations with $\cos\alpha$. Our results imply that the peculiar tidal alignments of the galaxy stellar spins may be caused by anisotropic occurrence of some baryonic process responsible for discharging stellar materials from the galaxies along the tidal major directions at $z<1$.

Proper elements are quasi-invariants of a Hamiltonian system, obtained through a normalization procedure. Proper elements have been successfully used to identify families of asteroids, sharing the same dynamical properties. We show that proper elements can also be used within space debris dynamics to identify groups of fragments associated to the same break-up event. The proposed method allows to reconstruct the evolutionary history and possibly to associate the fragments to a parent body. The procedure relies on different steps: (i) the development of a model for an approximate, though accurate, description of the dynamics of the space debris; (ii) the construction of a normalization procedure to determine the proper elements; (iii) the production of fragments through a simulated break-up event. We consider a model that includes the Keplerian part, an approximation of the geopotential, and the gravitational influence of Sun and Moon. We also evaluate the contribution of Solar radiation pressure and the effect of noise on the orbital elements. We implement a Lie series normalization procedure to compute the proper elements associated to semi-major axis, eccentricity and inclination. Based upon a wide range of samples, we conclude that the distribution of the proper elements in simulated break-up events (either collisions and explosions) shows an impressive connection with the dynamics observed immediately after the catastrophic event. The results are corroborated by a statistical data analysis based on the check of the Kolmogorov-Smirnov test and the computation of the Pearson correlation coefficient.

We study the main cosmological properties of the Generalized Chaplygin Gas (GCG) dark energy model at the background and perturbation levels. By using the latest cosmological data in both the background and perturbation levels, we implement a joint likelihood analysis to constrain the cosmological parameters of the model. Using the available expansion and growth rate data, we place constraints on the free parameters of the GCG model based on the statistical Markov chain Monte Carlo method. Then, the best-fit values of cosmological parameters and those of confidence regions are found. We obtain the best-fit value of the current expansion rate of the universe in the GCG model and show that it is in good agreement with the $\Lambda$CDM model. Moreover, the growth rate of matter perturbations is investigated in the context of a unified GCG model. It is shown that in this model, the dark energy component, like the $\Lambda$ sector in the $\Lambda$CDM model, can suppress the amplitude of matter perturbations. We show that the growth rate of perturbations in GCG parametrization is consistent with cluster-scale observations similar to the case of the concordance $\Lambda$CDM model. Our results show that the tension on $\sigma_{8}$ appeared in concordance model can be alleviated in GCG cosmology.

This paper presents a multiwavelength investigation of the Galactic HII IRAS 17149$-$3916. Using the Giant Meterwave Radio Telescope, India, first low-frequency radio continuum observations at 610 and 1280 MHz for this region are presented. The ionized gas emission displays an interesting cometary morphology which is likely powered by the early type source, E4 (IRS-1). The origin of the cometary morphology is discussed under the framework of the widely accepted bow shock, champagne flow, and clumpy cloud mechanisms. The mid- and far-infrared data from Spitzer-GLIMPSE and Herschel-Hi-GAL reveal a complex network of pillars, clumps, bubble, filaments, and arcs suggesting the profound influence of massive stars on the surrounding medium. Triggered star formation at the tip of an observed pillar structure is reported. High-resolution ALMA continuum data show a string of cores detected within the identified clumps. The core masses are well explained by thermal Jeans fragmentation and support the hierarchical fragmentation scenario. Four `super-Jeans' cores are identified which, at the resolution of the present data set, are suitable candidates to form high-mass stars.

Long gamma-ray bursts (GRBs) are considered to be originated from core collapse of massive stars. It is believed that the afterglow property is determined by the density of the material in the surrounding interstellar medium. Therefore, the circumburst density can be used to distinguish between an interstellar wind $n(R) \propto R^{\rm -k}$, and a constant density medium (ISM), $n(R) = const$. Previous studies with different afterglow samples, show that the circumburst medium of GRBs is neither simply supported by an interstellar wind, nor completely favored by an ISM. In this work, our new sample is consisted of 39 GRBs with smoothly onset bump-like features in early X-ray afterglows, in which 20 GRBs have the redshift measurements. By using a smooth broken power law function to fit the bumps of X-ray light curves, we derive the full width at half-maximum (FWHM) as the feature width ($\omega$), as well as the rise and decay time scales of the bumps ($T_{\rm r}$ and $T_{\rm d}$). The correlations between the timescales of X-ray bumps are similar to those found previously in the optical afterglows. Based on the fireball forward shock (FS) model of the thin shell case, we obtain the distribution of the electron spectral index $p$, and further constrain the medium density distribution index $k$. The new inferred $k$ is found to be concentrated at 1.0, with a range from 0.2 to 1.8. This finding is consistent with previous studies. The conclusion of our detailed investigation for X-ray afterglows suggests that the ambient medium of the selected GRBs is not homogeneous, i.e., neither ISM nor the typical interstellar wind.

The theory of General Relativity has successfully passed a large number of observational tests. The theory has been extensively tested in the weak-field regime with experiments in the Solar System and observations of binary pulsars. Recently, there have seen significant advancements in the study of the strong-field regime, which can now be tested with gravitational waves, X-ray data, and mm Very Long Baseline Interferometry observations. Here we summarize the state-of-the-art of the tests of General Relativity with black hole X-ray data, discussing its recent progress and future developments.

I have used mid-infrared (IR) photometry from the Wide-field Infrared Survey Explorer (WISE) to perform a census of circumstellar disks among ~10,000 candidate members of the Sco-Cen complex that were recently identified with data from the Gaia mission. IR excesses are detected for more than 1200 of the WISE counterparts that are within the commonly adopted boundary for Sco-Cen, ~400 of which are newly reported in this work. The richest population in Sco-Cen, UCL/LCC, contains the largest available sample of disks (>500) for any population near its age (~20 Myr). UCL/LCC also provides the tightest statistical constraints to date on the disk fractions of low-mass stars for any single age beyond that of Upper Sco (~11 Myr). For Upper Sco and UCL/LCC, I have measured the disk fractions as a function of spectral type. The disk fraction in Upper Sco is higher at later spectral types, which is consistent with the results for previous samples of candidate members. In UCL/LCC, that trend has become more pronounced; the disk fractions in UCL/LCC are lower than those in Upper Sco by factors of ~10, 5.7, and 2.5 at B7-K5.5, K6-M3.5, and M3.75-M6, respectively. The data in UCL/LCC also demonstrate that the disk fraction for low-mass stars remains non-negligible at an age of 20 Myr (0.09+/-0.01). Finally, I find no significant differences in the ages of disk-bearing and diskless low-mass stars in Upper Sco and UCL/LCC based on their positions in color-magnitude diagrams.

We present spectroscopy of 285 previously identified candidate members of populations in the Sco-Cen complex, primarily Ophiuchus, Upper Sco, and Lupus. The spectra are used to measure spectral types and diagnostics of youth. We find that 269 candidates exhibit signatures of youth in our spectra or previous data, which is consistent with their membership in Sco-Cen. We have constructed compilations of candidate members of Ophiuchus, Upper Sco, and Lupus that have spectral classifications and evidence of youth, which contain a total of 2274 objects. In addition, we have used spectra from previous studies to classify three sources in Ophiuchus that have been proposed to be protostellar brown dwarfs: ISO Oph 70, 200, and 203. We measure spectral types of early M from those data, which are earlier than expected for young brown dwarfs based on evolutionary models (>=M6.5) and instead are indicative of stellar masses (~0.6 Msun).

I have used high-precision photometry and astrometry from the early installment of the third data release of Gaia (EDR3) to perform a survey for members of the stellar populations within the Sco-Cen complex, which consist of Upper Sco, UCL/LCC, the V1062 Sco group, Ophiuchus, and Lupus. Among Gaia sources with sigma(pi)<1 mas, I have identified 10,509 candidate members of those populations. I have compiled previous measurements of spectral types, Li equivalent widths, and radial velocities for the candidates, which are available for 3169, 1420, and 1740 objects, respectively. In a subset of candidates selected to minimize field star contamination, I estimate that the contamination is <1% and the completeness is ~90% at spectral types of <=M6-M7 for the populations with low extinction (Upper Sco, V1062 Sco, UCL/LCC). I have used that cleaner sample to characterize the stellar populations in Sco-Cen in terms of their initial mass functions, ages, and space velocities. For instance, all of the populations in Sco-Cen have histograms of spectral types that peak near M4--M5, which indicates that they share similar characteristic masses for their initial mass functions (~0.15-0.2 Msun). After accounting for incompleteness, I estimate that the Sco-Cen complex contains nearly 10,000 members with masses above ~0.01 Msun. Finally, I also present new estimates for the intrinsic colors of young stars and brown dwarfs (<=20 Myr) in bands from Gaia EDR3, the Two Micron All Sky Survey, the Wide-field Infrared Survey Explorer, and the Spitzer Space Telescope.

Vamana is the mixture model framework that infers the chirp mass, mass ratio, and aligned spin distributions of the binary black hole population. We extend the mixing components to also model the redshift evolution of merger rate and report all the major one and two-dimensional features in the binary black hole population using the 69 gravitational wave signals detected with a false alarm rate $<1\mathrm{yr}^{-1}$ in the third Gravitational Wave Transient Catalog. Endorsing our previous report and corroborating recent report from LIGO Scientific, Virgo, and KAGRA Collaborations, we observe the chirp mass distribution has multiple peaks and a lack of mergers with chirp masses $10 \mathrm{-} 12M_\odot$. In addition, we observe aligned spins show mass dependence with heavier binaries exhibiting larger spins, mass ratio does not show a notable dependence on either the chirp mass or the aligned spin, and the redshift evolution of the merger rate for the peaks in the mass distribution is disparate. These features possibly reflect the astrophysics associated with the binary black formation channels. However, additional observations are needed to improve our limited confidence in them.

Unraveling the geological processes ongoing at Io's numerous sites of active volcanism requires high spatial resolution to, for example, measure the areal coverage of lava flows or identify the presence of multiple emitting regions within a single volcanic center. In de Kleer et al. (2017) we described observations with the Large Binocular Telescope (LBT) during an occultation of Io by Europa at ~6:17 UT on 2015 March 08, and presented a map of the temperature distribution within Loki Patera derived from these data. Here we present emission maps of three other volcanic centers derived from the same observation: Pillan Patera, Kurdalagon Patera, and the vicinity of Ulgen Patera/PV59/N Lerna Regio. The emission is localized by the light curves and resolved into multiple distinct emitting regions in two of the cases. Both Pillan and Kurdalagon Paterae had undergone eruptions in the months prior to our observations, and the location and intensity of the emission is interpreted in the context of the temporal evolution of these eruptions observed from other facilities. The emission from Kurdalagon Patera is resolved into two distinct emitting regions separated by only a few degrees in latitude that were unresolved by Keck observations from the same month.

CTA 102 is a blazar implying that its relativistic jet points towards Earth and emits synchrotron radiation produced by energetic particles gyrating in the magnetic field. This study aims to figure out the physical origins of radio flares in the jet, including the connection between the magnetic field and the radio flares. The dataset in the range 2.6-343.5 GHz was collected over a period of 5.5 years (2012 November 20-2018 September 23). During the data collection period, seven flares at 15 GHz with a range of the variability time-scale of roughly 26-171 days were detected. The quasi-simultaneous radio data were used to investigate the synchrotron spectrum of the source. We found that the synchrotron radiation is self-absorbed. The turnover frequency and the peak flux density of the synchrotron self-absorption (SSA) spectra are in the ranges of 42-167 GHz and 0.9-10.2 Jy, respectively. From the SSA spectra, we derived the SSA magnetic field strengths to be 9.20 mG, 12.28 mG, and 50.97 mG on 2013 December 24, 2014 February 28, and 2018 January 13, respectively. We also derived the equipartition magnetic field strengths to be in the range 24-109 mG. The equipartition magnetic field strengths are larger than the SSA magnetic field strengths in most cases, which indicates that particle energy mainly dominates in the jet. Our results suggest that the flares in the jet of CTA 102 originated due to particle acceleration. We propose the possible mechanisms of particle acceleration.

One objective of the Polstar spectropolarimetry mission is to characterize the degree of nonconservative mass transfer that occurs at various stages of binary evolution, from the initial mass reversal to the late Algol phase. The proposed instrument combines spectroscopic and polarimetric capabilities, where the spectroscopy can resolve Doppler shifts in UV resonance lines with 10 km/s precision, and polarimetry can resolve linear polarization with 1e-3 precision or better. The spectroscopy will identify absorption by mass streams seen in projection against the stellar disk as a function of orbital phase, as well as scattering from extended splash structures and a circumbinary disk of material that fails to be transferred conservatively. The polarimetry affects more the light coming from material not seen against the stellar disk, allowing the geometry of the scattering to be tracked, resolving ambiguities left by the spectroscopy and light-curve information. For example, nonconservative mass streams ejected in the polar direction will produce polarization of the opposite sign from conservative transfer accreting in the orbital plane. Also, time domain coverage over a range of phases of the binary orbit are well supported by the Polstar observing strategy. Combining these elements will significantly improve our understanding of the mass transfer process and the amount of mass that can escape from the system, an important channel for changing the final mass, and ultimate supernova, of the large number of massive stars found in binaries at close enough separation to undergo interaction.

Even though Pulsar Timing Arrays already have the potential to detect the gravitational wave background by finding a quadrupole correlation in the timing residuals, this goal has not yet been achieved. Motivated by some theoretical arguments, we analyzed some advantages of including the millisecond pulsars within globular clusters, especially those in their cores, in current and future Pulsar Timing Array projects for detecting the gravitational waves emitted by an ensemble of supermassive black holes.

Shocks, modelled over a broad range of parameters, are used to construct a new tool to deduce the mechanical energy and physical conditions from observed atomic or molecular emission lines. We compute magnetised, molecular shock models with velocities $V_s=5$-$80$ km s$^{-1}$, preshock proton densities $n_{\rm H}=10^2$-$10^6$ cm$^{-3}$, weak or moderate magnetic field strengths, and in the absence or presence of an external UV radiation field. We develop a simple emission model of an ensemble of shocks for connecting any observed emission lines to the mechanical energy and physical conditions of the system. For this range of parameters we find the full diversity (C-, C$^*$-, CJ-, and J-type) of magnetohydrodynamic shocks. H$_2$ and H are dominant coolants, with up to 30% of the shock kinetic flux escaping in Ly$\alpha$ photons. The reformation of molecules in the cooling tail means H$_2$ is even a good tracer of dissociative shocks and shocks that were initially fully atomic. For each shock model we provide integrated intensities of rovibrational lines of H$_2$, CO, and CH$^+$, atomic H lines, and atomic fine-structure and metastable lines. We demonstrate how to use these shock models to deduce the mechanical energy and physical conditions of extragalactic environments. As a template example, we interpret the CH$^+$(1-0) emission from the Eyelash starburst galaxy. A mechanical energy injection rate of at least $10^{11}$ $L_\odot$ into molecular shocks is required to reproduce the observed line. The low-velocity, externally irradiated shocks are at least an order magnitude more efficient than the most efficient shocks with no external irradiation, in terms of the total mechanical energy required. We predict differences of more than 2 orders of magnitude in intensities of the pure rotational lines of CO, Ly$\alpha$, metastable lines of O, S$^+$, and N, between representative models.

We calculate the relativistic corrections to hydrostatic X-ray masses for galaxy clusters in Kottler spacetime, which is the spherically symmetric solution to Einstein's equations in General relativity endowed with a cosmological constant. The hydrostatic masses for clusters (calculated assuming Newtonian gravity) have been found to be underestimated compared to lensing masses, and this discrepancy is known as hydrostatic mass bias. Since the relativistic hydrostatic X-ray masses are automatically lower than lensing masses, under the edifice of Kottler metric, we check if the hydrostatic mass bias problem gets alleviated using this {\it ansatz}. We consider a sample of 18 galaxy clusters for this pilot test. We find that the ratio of X-ray to lensing mass is close to unity even in Kottler spacetime. Therefore, the effect of relativistic corrections to hydrostatic X-ray masses for galaxy clusters is negligible.

The huge variety of planetary systems discovered in recent decades likely depends on the early history of their formation. In this contribution we introduce the FAUST Large Program, which focuses specifically on the early history of Solar-like protostars and their chemical diversity at scales of $\sim$ 50 au, where planets are expected to form. In particular, the goal of the project is to reveal and quantify the variety of chemical composition of the envelope/disk system at scales of 50 au in a sample of Class 0 and I protostars representative of the chemical diversity observed at larger scales. For each source, we propose a set of molecules able to: (1) disentangle the components of the 50-2000 au envelope/disk system; (2) characterise the organic complexity in each of them; (3) probe their ionization structure; (4) measure their molecular deuteration. The output will be a homogeneous database of thousands of images from different lines and species, i.e., an unprecedented source-survey of the chemical diversity of Solar-like protostars. FAUST will provide the community with a legacy dataset that will be a milestone for astrochemistry and star formation studies.

Observations show that the surface brightness of disc galaxies can be well-described by a single exponential (TI), up-bending (TIII) or down-bending (TII) profiles in the outskirts. Here we characterize the mass surface densities of simulated late-type galaxies from the EAGLE project according to their distribution of mono-age stellar populations, the star formation activity and angular momentum content. We find that the inner scale-lengths of TII galaxies correlate with their stellar spin parameter $\lambda$, while those with TI and TIII profiles show a correlation only for $\lambda>0.35$. The outer scale-lengths of TII and TIII discs show a positive trend with $\lambda$, albeit weaker for the latter. TII discs prefer fast rotator galaxies. With regards to the stellar age distribution, negative and U-shape age profiles are the most common for all disc types. Positive age profiles are determined by a more significant contributions of young stars in the central regions, which decrease rapidly in the outer parts. TII discs prefer relative higher contributions of old stars compared to other mono-age populations across the discs whereas TIII discs become progressively more dominated by intermediate age (2-6 Gyrs) stars for increasing radius. The change in slope of the age profiles is located after the break of the mass surface density. We find evidence of larger flaring for the old stellar populations in TI and TIII systems compared to TII. Overall, the relative distributions of mono-age stellar populations and the dependence of star formation on radius is found to modulate the different disc types and age profiles.

Clouds are an integral part of planetary atmospheres, with most planets hosting clouds. The understanding of not only the formation, but also the composition of clouds is crucial to the understanding of future observations. As observations of the planet's surface will remain very difficult, it is essential to link the observable high atmosphere gas and cloud composition to the surface conditions. We present a fast and simple chemical equilibrium (eq.) model for the troposphere of rocky exoplanets, which is in chemical and phase eq. with the crust. The hydrostatic eq. atmosphere is built from bottom to top. In each atmospheric layer chemical eq. is solved and all thermally stable condensates are removed, depleting the atmosphere above in the effected elements. These removed condensates build an upper limit for cloud formation and can be separated into high and low temperature condensates. The most important cloud condensates for 1000K>T>400K are KCl, NaCl, FeS, FeS2, FeO, Fe2O3, Fe3O4. For T<400K H2O, C, NH3, NH4Cl, NH4SH are thermally stable. For even lower temperatures of T<150K CO2, CH4, NH3, H2S become stable. The inclusion of clouds with trace abundances results in the thermal stability of a total of 72 condensates for atmospheres with the different surface conditions (300K<T<1000K and p=1bar,100bar). The different cloud condensates are not independent of each other, but follow sequences of condensation, which are robust against changes in crust composition, surface pressure, and surface temperature. Independent of the existence of water as a crust condensate H2O is a thermally stable cloud condensate for all investigated elemental abundances. However, the water cloud base depends on the hydration level of the crust. Therefore, the detection of water condensates alone does not necessarily imply stable water on the surface, even if the temperature could allow water condensation.

Solar-type stars generate spherical winds, which are pressure driven flows, that start subsonic, reach the sound speed at the sonic point and transition to supersonic flows. The sonic point, mathematically corresponds to a singularity of the system of differential equations describing the flow. In the problem of an isothermal wind, the Parker solution provides an exact analytical expression tuned appropriately so that the singularity does not affect the solution. However, if the wind is polytropic it is not possible to find an analytical solution and a numerical approach needs to be followed. We study solutions of spherical winds that are driven by pressure within a gravitational field. The solutions pass smoothly from the critical point and allow us to study the impact of the changes of the polytropic index to these winds. We explore the properties of these solutions as a function of the polytropic index and the boundary conditions used. We apply the Complex Plane Strategy (CPS) and we obtain numerically solutions of polytropic winds. This allows us to avoid the singularity appearing in the equations through the introduction complex variables and integration on the complex plane. Applying this method, we obtain solutions with physical behaviour at the stellar surface, the sonic point and at large distances from the star. We further explore the role of the polytropic index in the flow and the effect of mass--loss rate and temperature on the solution. We find that the increase of polytropic index as well as the decrease of flow parameter both yield to a smoother velocity profile and lower velocities and shifts the transition point from subsonic to supersonic behaviour further from the star. Finally, we verify that the increasing of coronal temperature yields higher wind velocities and a weaker dependence on polytropic index.

We propose a backscattering dominated prompt emission model for gamma-ray bursts (GRB) prompt phase in which the photons generated through pair annihilation at the centre of the burst are backscattered through Compton scattering by an outflowing stellar cork. We show that the obtained spectra are capable of explaining the low and high energy slopes as well as the distribution of spectral peak energies in their observed prompt spectra.

The identification and analysis of different variable sources is a hot issue in astrophysical research. The Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) spectroscopic survey has accumulated massive spectral data but contains no information about variable sources. Although a few related studies present variable source catalogs for the LAMOST, the studies still have a few deficiencies regarding the type and number of variable sources identified. In this study, we presented a statistical modeling approach to identify variable source candidates. We first crossed the Kepler, Sloan Digital Sky Survey (SDSS), and Zwicky Transient Facility (ZTF) catalogs to obtain light curves data of variable and non-variable sources. The data are then modeled statistically using commonly used variability parameters, respectively. And then, an optimal variable source identification model is determined using the Receiver Operating Characteristic (ROC) curve and four credible evaluation indices such as precision, accuracy, recall, and F1score. Based on this identification model, a catalog of LAMOST variable sources (including 631,769 variable source candidates with a probability greater than 95% and so on) is obtained. To validate the correctness of the catalog, we performed a two-by-two cross-comparison with the GAIA catalog and other published variable source catalogs. We achieved the correct rate ranging from 50% to 100%. Among the 123,756 sources cross-matched, our variable source catalog identifies 85,669 with a correct rate of 69%, which indicates that the variable source catalog presented in this study is credible.

We consider a dynamical model for dark energy based on an ultralight mass scalar field with very large-scale inhomogeneities. This model may cause observable impacts on the anisotropic properties of the cosmic microwave background (CMB) intensity and luminosity distance. We formulate the model as the cosmological perturbations of the superhorizon scales, focusing on the local region of our universe. Moreover, we investigated the characteristic properties of the late-time evolution of inhomogeneous dark energy. Our numerical solutions show that the model can mimic the standard $\Lambda$CDM cosmology while including spatially dependent dark energy with flexible ranges of the model parameters. We put a constraint on the amplitude of these inhomogeneities of the dark energy on very large scales with the observations of the CMB anisotropies. We also discuss their influence on the estimation of the luminosity distance.

In this paper we present a class of models in order to explain the production of Primordial Black Holes (PBHs) and Gravitational Waves (GWs) in the Universe. These models are based on no-scale theory. By breaking the SU(2,1)/SU(2)$\times$U(1) symmetry we fix one of the two chiral fields and we derive effective scalar potentials which are capable of generating PBHs and GWs. As it is known in the literature there is an important unification of the no-scale models, which leads to the Starobinsky effective scalar potential based on the coset SU(2,1)/SU(2)$\times$U(1). We use this unification in order to additionally explain the generation of PBHs and GWs. Concretely, we modify well-known superpotentials, which reduce to Starobinsky-like effective scalar potentials. Thus, we derive scalar potentials which, on the one hand, explain the production of PBHs and GWs and, on the other hand, they conserve the transformation laws, which yield from the parametrization of the coset SU(2,1)/SU(2)$\times$U(1) as well as the unification between the models which are yielded this coset. We numerically evaluate the scalar power spectra with the effective scalar potential based on this theory. Furthermore, we evaluate the fractional abundances of PBHs by comparing two methods the Press-Schechter approach and the peak theory, while focusing on explaining the dark matter in the Universe. By using the resulting scalar power spectrum we evaluate the amount of GWs. All models are in complete consistence with Planck constraints.

Calculations of stellar evolution up to the early white dwarf stage were carried out for stars with mass on the main sequence $M_0=0.82M_\odot$, $0.85M_\odot$, $0.9M_\odot$ and with initial abundances of helium and heavier elements $Y=0.25$ and $Z=10^{-3}$, respectively. For each value of $M_0$ the AGB and post--AGB evolutionary phases were computed with three values of the mass loss parameter in the Bl\"ocker formula: $\eta_B=0.02$, 0.05 and 0.1. The variable star RU~Cam with pulsation period $\Pi\approx 22$ day is shown to be in the post--AGB stage and the pulsation amplitude decrease in years 1962--1963 is due to movement of the star across the HR diagram beyond the pulsation instability region. Theoretical estimates of the mass and the luminosity of RU~Cam are $0.524M_\odot\le M\le 0.532M_\odot$ and $2.20\times 10^3L_\odot\le L\le 2.33\times 10^3L_\odot$, respectively. Hydrodynamic calculations of nonlinear stellar pulsations show that while the star approaches the instability boundary a significant reduction ($\approx 90\%$) in the pulsation amplitude occurs for nearly two years with subsequent slow decay of low--amplitude oscillations. Solution of the equations of hydrodynamics with time--dependent inner boundary conditions describing evolutionary changes in the radius and the luminosity at the bottom of the pulsating envelope allows us to conclude that decay of radial oscillations in RU~Cam is accompanied by the effect of oscillation hysteresis. In particular, the stage of large--amplitude limit cycle oscillations extends by $\approx 12$ years and the subsequent stage of decaying small--amplitude oscilations spreads beyond the formal boundary of pulsation instability.

We present a new publicly available code, $\texttt{class_rot}$, which modifies $\texttt{class}$ to enable fast non-perturbative calculation of cosmic microwave background polarization power spectra due to both isotropic and anisotropic polarization rotation from cosmic birefringence. Cosmic birefringence can arise from new parity-violating physics such as axion dark matter with a Chern-Simons coupling to photons or Faraday rotation due to a primordial magnetic field. Constraints on these effects can be obtained by comparing measurements to precise numerical calculations of the polarization power spectra. We describe the implementation of $\texttt{class_rot}$ in terms of both mathematical formalism and coding architecture. We also provide usage examples and demonstrate the accuracy of the code by comparing with simulations.

The aim of the present study is to show the varying asymmetries during the decay of sunspot groups. The source of input data is the SoHO/MDI-Debrecen Database (SDD) sunspot catalog that contains the magnetic polarity data for time interval 1996-2010. Several types of asymmetries were examined on the selected sample of 142 sunspot groups. The leading-following asymmetry increases in three phases during the decay and exhibits anticorrelation with size. It is also related to a hemispheric asymmetry, during the decay the area asymmetry index has higher values in the southern hemisphere which may be due to the higher activity level in the southern hemisphere in cycle 23. The total umbral area is inversely proportional to the umbra/penumbra ratio but it is directly proportional to the umbral decay rate. During the decay the umbra/penumbra (U/P) ratio decreases unambiguously in the trailing parts but in most cases in the leading parts as well. The U/P variation is a consequence of the different depths of umbral and penumbral fields.

CO2 ice is present on the trailing hemisphere of Ariel but is mostly absent from its leading hemisphere. The leading/trailing hemispherical asymmetry in the distribution of CO2 ice is consistent with radiolytic production of CO2, formed by charged particle bombardment of H2O ice and carbonaceous material in Ariel's regolith. This longitudinal distribution of CO2 on Ariel was previously characterized using 13 near-infrared reflectance spectra collected at 'low' sub-observer latitudes between 30S to 30N. Here, we investigated the distribution of CO2 ice on Ariel using 18 new spectra: two collected over low sub-observer latitudes, five collected at 'mid' sub-observer latitudes (31 - 44N), and eleven collected over 'high' sub-observer latitudes (45 - 51N). Analysis of these data indicates that CO2 ice is primarily concentrated on Ariel's trailing hemisphere. However, CO2 ice band strengths are diminished in the spectra collected over mid and high sub-observer latitudes. This sub-observer latitudinal trend may result from radiolytic production of CO2 molecules at high latitudes and subsequent migration of this constituent to low latitude cold traps. We detected a subtle feature near 2.13 microns in two spectra collected over high sub-observer latitudes, which might result from a 'forbidden' transition mode of CO2 ice that is substantially stronger in well mixed substrates composed of CO2 and H2O ice, consistent with regolith-mixed CO2 ice grains formed by radiolysis. Additionally, we detected a 2.35-micron feature in some low sub-observer latitude spectra, which might result from CO formed as part of a CO2 radiolytic production cycle.

We describe how to implement the spectral kurtosis method of interference removal (zapping) on a digitized signal of averaged power values. Spectral kurtosis is a hypothesis test, analogous to the t-test, with a null hypothesis that the amplitudes from which power is formed belong to a `good' distribution -- typically Gaussian with zero mean -- where power values are zapped if the hypothesis is rejected at a specified confidence level. We derive signal-to-noise ratios (SNRs) as a function of amount of zapping for folded radio pulsar observations consisting of a sum of signals from multiple telescopes in independent radio-frequency interference (RFI) environments, comparing four methods to compensate for lost data with coherent (tied-array) and incoherent summation. For coherently summed amplitudes, scaling amplitudes from non-zapped telescopes achieves a higher SNR than replacing zapped amplitudes with artificial noise. For incoherently summed power values, the highest SNR is given by scaling power from non-zapped telescopes to maintain a constant mean. We use spectral kurtosis to clean a tied-array radio pulsar observation by the Large European Array for Pulsars (LEAP): the signal from one telescope is zapped with time and frequency resolutions of 6.25 $\mu$s and 0.16 MHz, removing interference along with 0.27 per cent of `good' data, giving an uncertainty of 0.25 $\mu$s in pulse time of arrival (TOA) for PSR J1022+1001. We use a single-telescope observation to demonstrate recovery of the pulse profile shape, with 0.6 per cent of data zapped and a reduction from 1.22 to 0.70 $\mu$s in TOA uncertainty.

The advent of all-sky facilities, such as the Neil Gehrels Swift observatory, the All Sky Automated Search for Supernovae (ASASSN), eROSITA and Gaia has led to a new appreciation of the importance of transient sources in solving outstanding astrophysical questions. Identification and catalogue cross-matching of transients has been eased over the last two decades by the Virtual Observatory but we still lack a client capable of providing a seamless, self-consistent, analysis of all observations made of a particular object by current and historical facilities. HILIGT is a web-based interface which polls individual servers written for XMM-Newton, INTEGRAL and other missions, to find the fluxes, or upper limits, from all observations made of a given target. These measurements are displayed as a table or a time series plot, which may be downloaded in a variety of formats. HILIGT currently works with data from X-ray and Gamma-ray observatories.

Big-bang nucleosynthesis (BBN), a pillar of modern cosmology, begins with the trailblazing 1948 paper of Alpher, Bethe and Gamow \cite{abc} which proposed nucleosynthesis during an early ($t \sim 1 - 1000\,$sec), radiation-dominated phase of the Universe to explain the abundances of the chemical elements. While their model was flawed, they called attention to the importance of a hot beginning and the radiation now known as the cosmic microwave background (CMB). In subsequent papers they made estimates for the temperature of the CMB, from 5\,K to 50\,K, all based upon wrong physics. To illustrate the prediction {\it that could have been made} and to elucidate the key physics underpinning BBN, I show that in the absence of a detailed model of BBN and the desired abundances, at best one could have estimated an {\it upper} limit to the CMB temperature of between $10\,$K and $60\,$K, predicated upon the assumption of some nucleosynthesis.

The chemical element nickel is of particular interest in stellar physics. In the layers in which the Fe-peak elements dominate the mean opacity (the so-called Z-bump), Ni is the second contributor to the Rosseland opacity after iron, according to the Opacity Project data. Reliable nickel cross sections are therefore mandatory for building realistic stellar models, especially for main-sequence pulsators such as $\beta$ Cep and slowly pulsating B stars, whose oscillations are triggered by the $\kappa$-mechanism of the Fe-peak elements. Unfortunately, the Opacity Project data for Ni were extrapolated from those of Fe, and previous studies have shown that they were underestimated in comparison to detailed calculations. We investigate the impact of newly computed monochromatic cross sections on the Rosseland mean opacity of Ni and on the structure of main-sequence massive pulsators. We compare our results with the widely used Opacity Project and OPAL data. Monochromatic cross sections for Ni were obtained with the SCO-RCG code. The Toulouse-Geneva evolution code was used to build the stellar models. With the new data, the Rosseland opacities of Ni are roughly the same as those of the Opacity Project or OPAL at high temperatures ($\log\ T>6$). At lower temperatures, significant departures are observed; the ratios are up to six times higher with SCO-RCG. These discrepancies span a wider temperature range in the comparison with OPAL than in comparison with the Opacity Project. For massive star models, the results of the comparison with a structure computed with Opacity Project data show that the Rosseland mean of the global stellar mixture is only marginally altered in the Z-bump. The maximum opacity is shifted towards slightly more superficial layers. A new maximum appears in the temperature derivative of the mean opacity, and the driving of the pulsations should be affected.

The Astrophysics Source Code Library (ASCL ascl.net), started in 1999, is a free open registry of software used in refereed astronomy research. Over the past few years, it has spearheaded an effort to form a consortium of scientific software registries and repositories. In 2019 and 2020, ASCL contacted editors and maintainers of discipline and institutional software registries and repositories in math, biology, neuroscience, geophysics, remote sensing, and other fields to develop a list of best practices for these research software resources. At the completion of that project, performed as a Task Force for a FORCE11 working group, members decided to form SciCodes as an ongoing consortium. This presentation covered the consortium's work so far, what it is currently working on, what it hopes to achieve for making scientific research software more discoverable across disciplines, and how the consortium can benefit astronomers.

The Sun's magnetic field varies in multiple time scales. Observations show that the minimum between cycles 24 and 25 was the second consecutive minimum which was deeper and wider than several earlier minima. Since the active regions observed at the Sun's surface are manifestations of the magnetic field generated in the interior, it is crucial to investigate/understand the dynamics below the surface. In this context, we report, by probing the solar interior with helioseismic techniques applied to long-term oscillations data from the Global Oscillation Network Group (GONG), that the seismic minima in deeper layers have been occurring about a year earlier than that at the surface for the last two consecutive solar cycles. Our findings also demonstrate a decrease in strong magnetic fields at the base of the convection zone, the primary driver of the surface magnetic activity. We conclude that the magnetic fields located in the core and near-surface shear layers, in addition to the tachocline fields, play an important role in modifying the oscillation frequencies. This provides evidence that further strengthens the existence of a relic magnetic field in the Sun's core since its formation.

Low-frequency quasi-periodic oscillations (LFQPOs) have been routinely observed in black hole X-ray binaries (BHXRBs). These LFQPOs can be explained by axisymmetric shock oscillation in accretion flow around a rotating black hole. We address the physical origin of Type-C LFQPOs in BHXRBs observed by the Rossi X-ray Timing Explorer satellite considering a minimum number of free parameters, namely, specific energy and specific angular momentum of the infalling matter for a given set of BH mass and spin parameter. We apply the solution for a large number of BH candidates to further strengthen the scenario of an anti-correlation between the QPO frequency and the location of the shock. Our study also confirms that Compton cooling can be sufficient to explain the observed QPOs.

Here we present a low frequency (170-200MHz) search for coherent radio emission associated with nine short GRBs detected by the Swift and/or Fermi satellites using the Murchison Widefield Array (MWA) rapid-response observing mode. The MWA began observing these events within 30 to 60s of their high-energy detection, enabling us to capture any dispersion delayed signals emitted by short GRBs for a typical range of redshifts. We conducted transient searches at the GRB positions on timescales of 5s, 30s and 2min, resulting in the most constraining flux density limits on any associated transient of 0.42, 0.29, and 0.084Jy, respectively. We also searched for dispersed signals at a temporal and spectral resolution of 0.5s and 1.28MHz but none were detected. However, the fluence limit of 80-100Jy ms derived for GRB 190627A is the most stringent to date for a short GRB. We compared the fluence and persistent emission limits to short GRB coherent emission models, placing constraints on key parameters including the radio emission efficiency of the nearly merged neutron stars ($\lesssim10^{-4}$), the fraction of magnetic energy in the GRB jet ($\lesssim2\times10^{-4}$), and the radio emission efficiency of the magnetar remnant ($\lesssim10^{-3}$). Comparing the limits derived for our full GRB sample to the same emission models, we demonstrate that our 30-min flux density limits were sensitive enough to theoretically detect the persistent radio emission from magnetar remnants up to a redshift of $z\sim0.6$. Our non-detection of this emission could imply that some GRBs in the sample were not genuinely short or did not result from a binary neutron star merger, the GRBs were at high redshifts, these mergers formed atypical magnetars, the radiation beams of the magnetar remnants were pointing away from Earth, or the majority did not form magnetars but rather collapse directly into black holes.

Over the past decades, many explosion scenarios for Type Ia supernovae have been proposed and investigated including various combinations of deflagrations and detonations in white dwarfs of different masses up to the Chandrasekhar mass. One of these is the gravitationally confined detonation model. In this case a weak deflagration burns to the surface, wraps around the bound core, and collides at the antipode. A subsequent detonation is then initiated in the collision area. Since the parameter space for this scenario, that is, varying central densities and ignition geometries, has not been studied in detail, we used pure deflagration models of a previous parameter study dedicated to Type Iax supernovae as initial models to investigate the gravitationally confined detonation scenario. We aim to judge whether this channel can account for one of the many subgroups of Type Ia supernovae, or even normal events. To this end, we employed a comprehensive pipeline for three-dimensional Type Ia supernova modeling that consists of hydrodynamic explosion simulations, nuclear network calculations, and radiative transfer. The observables extracted from the radiative transfer are then compared to observed light curves and spectra. The study produces a wide range in masses of synthesized 56 Ni ranging from 0.257 to 1.057 $M_\odot$ , and, thus, can potentially account for subluminous as well as overluminous Type Ia supernovae in terms of brightness. However, a rough agreement with observed light curves and spectra can only be found for 91T-like objects. Although several discrepancies remain, we conclude that the gravitationally confined detonation model cannot be ruled out as a mechanism to produce 91T-like objects. However, the models do not provide a good explanation for either normal Type Ia supernovae or Type Iax supernovae.

We assess the dominant low-redshift anisotropic signatures in the distance-redshift relation and redshift drift signals. We adopt general-relativistic irrotational dust models allowing for gravitational radiation -- the `quiet universe models' -- which are extensions of the silent universe models. Using cosmological simulations evolved with numerical relativity, we confirm that the quiet universe model is a good description on scales larger than those of collapsing structures. With this result, we reduce the number of degrees of freedom in the fully general luminosity distance and redshift drift cosmographies by a factor of $\sim 2$ and $\sim 2.5$, respectively, for the most simplified case. We predict a dominant dipolar signature in the distance-redshift relation for low-redshift data, with direction along the gradient of the large-scale density field. Further, we predict a dominant quadrupole in the anisotropy of the redshift drift signal, which is sourced by the electric Weyl curvature tensor. The signals we predict in this work should be tested with present and near-future cosmological surveys.

In this study, we model a pulsar as a general relativistic oblique rotator, where the oblique rotator is a rotationally deformed neutron star whose rotation and magnetic axis are inclined at an angle. The oblique rotator spins down, losing rotational energy through the magnetic poles. The magnetic field is assumed to be dipolar; however, the star has a non-zero azimuthal component due to the misalignment. The magnetic field induces an electric field for a force-free condition. The magnetic field decreases as the misalignment increases and is minimum along the equatorial plane of the star. In contrast, the electric field remains almost constant initially but decreases rapidly at a high misalignment angle. The charge separation at the star surface is qualitatively similar to that of Newtonian calculation. We find that the power loss for a general relativistic rotator is minimum for either an aligned or an orthogonal rotator, which contrasts with Newtonian calculation, where the power loss increases with an increase in the misalignment angle.

We are presenting a novel, Deep Learning based approach to estimate the normalized broadband spectral energy distribution (SED) of different stellar populations in synthetic galaxies. In contrast to the non-parametric multiband source separation algorithm, SCARLET - where the SED and morphology are simultaneously fitted - in our study we provide a morphology-independent, statistical determination of the SEDs, where we only use the color distribution of the galaxy. We developed a neural network (sedNN) that accurately predicts the SEDs of the old, red and young, blue stellar populations of realistic synthetic galaxies from the color distribution of the galaxy-related pixels in simulated broadband images. We trained and tested the network on a subset of the recently published CosmoDC2 simulated galaxy catalog containing about 3,600 galaxies. The model performance was compared to the results of SCARLET, where we found that sedNN can predict the SEDs with 4-5% accuracy on average, which is about two times better than applying SCARLET. We also investigated the effect of this improvement on the flux determination accuracy of the bulge and disk. We found that using more accurate SEDs decreases the error in the flux determination of the components by approximately 30%.

FAIR principles have the intent to act as a guideline for those wishing to enhance the reusability of their data holdings and put specific emphasis on enhancing the ability of machines to automatically find and use the data, in addition to supporting its reuse by individuals. Interoperability, one core of these principles, especially when dealing with automated systems' ability to interface with each other, requires open standards to avoid restrictions that negatively impact the user's experience. Open-ness of standards is best supported when the governance itself is open and includes a wide range of community participation. In this contribution we report our experience with the FAIR principles, interoperable systems and open governance in astrophysics. We report on activities that have matured within the ESCAPE project with a focus on interfacing the EOSC architecture and Interoperability Framework.

Recent studies have shown that both 21cm HI absorption and soft X-ray absorption serve as excellent tracers of the dense and dusty gas near the active nucleus of young radio galaxies, offering new insight into the physical nature of the circumnuclear medium. To date, a correlation between the column densities derived using these absorption processes has been observed within Compact Steep Spectrum (CSS) and Gigahertz-Peaked Spectrum (GPS) radio sources. While it is possible that this correlation exists within the broader radio population, many samples of radio galaxies make this difficult to test due to selection effects. This paper explores the possibility of a correlation in the broader radio population by analysing a historic sample of 168 radio sources compiled from the literature in such a way so as to minimise selection bias. From this historic sample we conclude that there is some evidence for a correlation between HI and soft X-ray absorption outside of peaked spectrum sources, but that the selection bias towards these sources makes further analysis difficult using current samples. We discuss this in the context of the SEAFOG project and how forthcoming data will change the landscape of future absorption studies.

The final architecture of planetary systems depends on the extraction of angular momentum and mass-loss processes of the discs in which they form. Theoretical studies proposed that magnetized winds launched from the discs (MHD disc winds) could govern accretion and disc dispersal. In this work, we revisit the observed disc demographics in the framework of MHD disc winds, combining analytical solutions of disc evolution and a disc population synthesis approach. We show that MHD disc winds alone can account for both disc dispersal and accretion properties. The decline of disc fraction over time is reproduced by assuming that the initial accretion timescale (a generalization of the viscous timescale) varies from disc to disc and that the decline of the magnetic field strength is slower than that of the gas. The correlation between accretion rate and disc mass, and the dispersion of the data around the mean trend as observed in Lupus is then naturally reproduced. The model also accounts for the rapidity of the disc dispersal. This paves the way for planet formation models in the paradigm of wind-driven accretion.

The Cherenkov Telescope Array (CTA) will be the world's largest and most sensitive ground-based gamma-ray observatory in the energy range from a few tens of GeV to tens of TeV. The LST-1 prototype, currently in its commissioning phase, is the first of the four largest CTA telescopes, that will be built in the northern site of CTA in La Palma, Canary Islands, Spain. In this contribution, we present a full-image reconstruction method using a modified InceptionV3 deep convolutional neural network applied on non-parametrized shower images. We evaluate the performance of optimized networks on Monte Carlo simulations of LST-1 shower images, and compare the results with the performance of the standard reconstruction method. We also show how both methods work on real-data reconstruction.

The presence of HI gas in galaxies is inextricably linked to their morphology and evolution. This paper aims to understand the HI content of the already identified 2210 dwarfs located in the low-to-moderate density environments of the MATLAS deep imaging survey. We combine the HI observations from the ATLAS$^{3D}$ survey, with the extragalactic HI sources from the ALFALFA survey, to extract the HI line width, velocity and mass of the MATLAS dwarfs. From the 1773 dwarfs in our sample with available HI observations, 8% (145) have an HI line detection. The majority of the dwarfs show irregular morphology, while 29% (42) are ellipticals, the largest sample of HI-bearing dwarf ellipticals (dEs) to date. Of the HI dwarf sample, 2% (3) are ultra-diffuse galaxies (UDGs), 12% have a transition-type morphology, 5% are tidal dwarf candidates, and 10% appear to be disrupted objects. In our optically selected sample, 9.5% of the dEs, 7% of the UDGs and 10% of the classical dwarfs are HI-bearing. The HI-bearing dwarfs have on average bluer colors than the dwarfs without detected HI. We find relations between the stellar and HI masses, gas fraction, color and absolute magnitude consistent with previous studies of dwarfs probing similar masses and environments. For 79% of the dwarfs identified as satellites of massive early-type galaxies, we find that the HI mass increases with the projected distance to the host. Using the HI line width, we estimate dynamical masses and find that 5% (7) of the dwarfs are dark matter deficient.

The CaFe Project involves the study of the properties of the low ionization emission lines (LILs) pertaining to the broad-line region (BLR) in active galaxies. These emission lines, especially the singly-ionized iron (Fe II) in the optical and the corresponding singly-ionized calcium (Ca II) in the near-infrared (NIR) are found to show a strong correlation in their emission strengths, i.e., with respect to the broad H$\beta$ emission line, the latter also belonging to the same category of LILs. We outline the progress made in the past years that has developed our understanding of the location and the efficient production of these emission lines. We have yet to realize the full potential of Ca II emission and its connection to the black hole and the BLR parameters which can be useful in - (1) the classification of Type-1 active galactic nuclei (AGNs) in the context of the main sequence of quasars, (2) to realize an updated radius-luminosity relation wherein the inclusion of the strength of this emission line with respect to H$\beta$ can be an effective tracer of the accretion rate of the AGN, and, (3) the close connection of Ca II to Fe II can allow us to use the ratio of the two species to quantify the chemical evolution in these active galaxies across cosmic time. In this paper, we use our current sample and utilize a non-linear dimensionality reduction technique - t-distributed Stochastic Neighbour Embedding (tSNE), to understand the clustering in our dataset based on direct observables.

IXPE scientific payload comprises of three telescopes, each composed of a mirror and a photoelectric polarimeter based on the Gas Pixel Detector design. The three focal plane detectors, together with the unit which interfaces them to the spacecraft, are named IXPE Instrument and they will be built and calibrated in Italy; in this proceeding, we will present how IXPE Instrument will be calibrated, both on-ground and in-flight. The Instrument Calibration Equipment is being finalized at INAF-IAPS in Rome (Italy) to produce both polarized and unpolarized radiation, with a precise knowledge of direction, position, energy and polarization state of the incident beam. In flight, a set of four calibration sources based on radioactive material and mounted on a filter and calibration wheel will allow for the periodic calibration of all of the three IXPE focal plane detectors independently. A highly polarized source and an unpolarized one will be used to monitor the response to polarization; the remaining two will be used to calibrate the gain through the entire lifetime of the mission.

Using bonafide black hole (BH) mass estimates from reverberation mapping and the line ratio [Si VI]1.963$\mu$m/Br$\gamma_{\rm broad}$ as tracer of the Active Galactic Nuclei (AGN) ionizing continuum, we find a novel BH-mass scaling relation of the form log($M_{\rm BH}) = (6.40\pm 0.17) - (1.99\pm 0.37) \times$ log ([Si VI]/Br$\gamma_{\rm broad})$, over the BH mass interval, $10^6 - 10^8$ M$_{\odot}$. The current sample consists of 21 Type-1 AGNs with the overall dispersion in our scaling relation at 0.47 dex, one that emulates the well-established M-$\sigma$ relation which shows a dispersion $\sim$0.44 dex. The new scaling offers an economic, and physically motivated alternative for BH estimate using single epoch spectra, avoiding large telescope time (reverberation mapping) or absolute flux calibration (the continuum luminosity method). With the advent of big surveys in the infrared in the near future, we aim to reduce the scatter in the relation by supplementing with more sources.

Broad-band spectra of active galaxies contain a wide range of information that help reveal the nature and activity of the central continuum source and their immediate surroundings. Understanding the evolution of metals in the spectra of Active Galactic Nuclei (AGN) and linking them with the various fundamental black hole (BH) parameters, for example, BH mass, the bolometric luminosity of the source, its accreting power, can help address the connection between the growth of the BH across cosmic time. We investigate the role of selected metallicity indicators utilizing the rich spectroscopic database of emission lines covering a wide range in redshift in the recent spectroscopic data release of Sloan Digital Sky Survey (SDSS). We make careful filtering of the parent sample to prepare a pair of high-quality, redshift-dependent sub-samples and present the first results of the analysis here. To validate our findings from the simple correlations, we execute and evaluate the performance of a linear dimensionality reduction technique - principal component analysis (PCA), over our sub-samples and present the projection maps highlighting the primary drivers of the observed correlations. The projection maps also allow us to isolate peculiar sources of potential interest.

Astrophysical shocks are very common and are interesting as they are responsible for particle acceleration in supernovas, blazers, and neutron stars. In this work, we study general relativistic shocks from the frame of the front. We derive the jump conditions and the Taub adiabat equation for both the space-like and time-like shocks. We solve these equations in a neutron star system where the shock is followed by a combustion front which is deconfining hadronic matter to quark matter. The maximum mass of the daughter quark star (generated from the combustion of the parent neutron star) is consistent with the maximum mass limit for the EoS sequence. We find that matter velocities for GR shocks under suitable conditions can break the speed of light limit indicating a very fast combustion process. Also, the matter velocities imply that for space-like shocks the combustion process is most probably a deflagration and for time-like shocks, it is a detonation and can even proceed with velocities that are super-luminous.

We investigate the orbital dynamics of circumbinary planetary systems with two planets around a circular or eccentric orbit binary. The orbits of the two planet are initially circular and coplanar to each other, but misaligned with respect to the binary orbital plane. The binary-planet and planet-planet interactions result in complex planet tilt oscillations. We use analytic models and numerical simulations to explore the effects of various values of the planet semi-major axes, binary eccentricity, and initial inclination. Around a circular orbit binary, secular tilt oscillations are driven by planet-planet interactions and are periodic. In that case, planets undergo mutual libration if close together and circulation if far apart with an abrupt transition at a critical separation. Around an eccentric orbit binary, secular tilt oscillations are driven by both planet-planet interactions and binary-planet interactions. These oscillations generally display more than one frequency and are generally not periodic. The transition from mutual planet libration to circulation is not sharp and there is a range of separations for which the planets are on orbits that are sometimes mutually librating and sometimes circulating. In addition, at certain separations, there are resonances for which tilt oscillations are complicated but periodic. For planets that are highly misaligned with respect to an eccentric orbit binary, there are stationary (non-oscillating) tilt configurations that are generalisations of polar configurations for the single planet case. Tilt oscillations of highly inclined planets occur for initial tilts that depart from the stationary configuration.

Dust polarization observations are a powerful, practical tool to probe the geometry (and to some extent, the strength) of magnetic fields in star-forming regions. In particular, Planck polarization data have revealed the importance of magnetic fields on large scales in molecular clouds. However, due to insufficient resolution, Planck observations are unable to constrain the B-field geometry on prestellar and protostellar scales. The high angular resolution of 11.7 arcsec provided by NIKA2-Pol 1.15 mm polarimetric imaging, corresponding to $\sim$ 0.02 pc at the distance of the Orion molecular cloud (OMC), makes it possible to advance our understanding of the B-field morphology in star-forming filaments and dense cores (IRAM 30m large program B-FUN). The commissioning of the NIKA2-Pol instrument has led to several challenging issues, in particular, the instrumental polarization or intensity-to-polarization (leakage) effect. In the present paper, we illustrate how this effect can be corrected for, leading to reliable exploitable data in a structured, extended source such as OMC-1. We present a statistical comparison between NIKA2-Pol and SCUBA2-Pol2 results in the OMC-1 region. We also present tentative evidence of local pinching of the B-field lines near Orion-KL, in the form of a new small-scale hourglass pattern, in addition to the larger-scale hourglass already seen by other instruments such as Pol2.

We use 20 years of magnetic field measurements from the Oersted, CHAMP and Swarm satellite missions, supplemented by calibrated platform magnetometer data from the CryoSat2 satellite, to study time variations of the Earth's core field at satellite altitude and at the core-mantle boundary (CMB). From the satellite data we derive composite time series of the core field secular variation (SV) with 4month cadence, at 300 globally distributed Geomagnetic Virtual Observatories (GVO). GVO radial SV series display regional fluctuations with 5-10 years duration and amplitudes reaching 20 nT/yr, most notably at low latitudes over Indonesia (2014), over South America and the South Atlantic (2007, 2011 and 2014), and over the central Pacific (2017). Applying the Subtractive Optimally Localized Averages (SOLA) method, we map the SV at the CMB as a collection of locally averaged SV estimates. We demonstrate that using 2-year windows of CryoSat2 data, it is possible to reliably estimate the SV and its time derivative, the secular acceleration (SA), at the CMB, with a spatial resolution, corresponding to spherical harmonic degree 10. Along the CMB geographic equator, we find strong SA features under Indonesia from 2011-2014, under central America from 2015 to 2019, and sequences of SA with alternating sign under the Atlantic during 2004-2019. We find that data from CryoSat2 make a valuable contribution to the emerging picture of sub-decadal core field variations. Using 1 year windows of data from the Swarm satellites, it is possible to study SA changes at low latitudes on timescales down to 1 year, with spatial resolution corresponding to spherical harmonic degree 10. We find strong positive and negative SA features appearing side-by-side in the Pacific in 2017, and thereafter drift westward.

The coalescences of massive black hole binaries (MBHBs) are one of the main targets of space-based gravitational wave observatories. Such gravitational wave sources are expected to be accompanied by electromagnetic emission. Low latency time of gravitational wave searches and accurate sky localization are keys in triggering successful follow-up observations on the electromagnetic counterparts. Here we present a deep learning method for the first time to rapidly search for MBHB signals in the strain data. Our model is capable to process 1-year of data in just several seconds, identifying all MBHB coalescences with no false alarms. We test the performance of our model on the simulated data from the LISA data challenge. We demonstrate that the model shows a robust resistance for a wide range of generalization for MBHB signals. This method is supposed to be an effective approach, which combined the advances of artificial intelligence to open a new pathway for space-based gravitational wave observations.

Astronomy is entering an era of data-driven discovery, due in part to modern machine learning (ML) techniques enabling powerful new ways to interpret observations. This shift in our scientific approach requires us to consider whether we can trust the black box. Here, we overview methods for an often-overlooked step in the development of ML models: building community trust in the algorithms. Trust is an essential ingredient not just for creating more robust data analysis techniques, but also for building confidence within the astronomy community to embrace machine learning methods and results.

We analyze the evolution of current driven kink instabilities of a highly magnetized relativistic plasma column, focusing in particular on its dissipation properties. The instability evolution leads to the formation of thin current sheets where the magnetic energy is dissipated. We find that the total amount of dissipated magnetic energy is independent of the dissipation properties. Dissipation occurs in two stages: a peak when the instability saturates, which is characterized by the formation of a helicoidal current sheet at the boundary of the deformed plasma column, followed by a weaker almost flat phase, in which turbulence develops. The detailed properties of these two phases depend on the equilibrium configuration and other parameters, in particular on the steepness of the pitch radial profile, on the presence of an external axial magnetic field and on the amount of magnetization. These results are relevant for high energy astrophysical sources, since current sheets can be the sites of magnetic reconnection where particles can be accelerated to relativistic energies and give rise to the observed radiation.

We explore the linearly quantised primordial power spectra associated with palindromic universes. Extending the results of Lasenby et al. [1] and Bartlett et al. [2], we improve the modelling of recombination and include higher orders in the photonic Boltzmann hierarchy. In so doing, we find that the predicted power spectra become largely consistent with observational data. The improved recombination modelling involves developing further techniques for dealing with the future conformal boundary, by integrating the associated perturbation equations both forwards and backwards in conformal time. The resulting wavevector quantisation gives a lowest allowed wavenumber ${k_0 = 9.93 \times 10^{-5} \textrm{Mpc}^{-1}}$ and linear spacing ${\Delta k = 1.63 \times 10^{-4} \textrm{Mpc}^{-1}}$, providing fits consistent with observational data equivalent in quality to the $\Lambda$CDM model.

The sources of ultra-high-energy cosmic rays are still unknown, but assuming standard physics, they are expected to lie within a few hundred megaparsecs from us. Indeed, over cosmological distances cosmic rays lose energy to interactions with background photons, at a rate depending on their mass number and energy and properties of photonuclear interactions and photon backgrounds. The universe is not homogeneous at such scales, hence the distribution of the arrival directions of cosmic rays is expected to reflect the inhomogeneities in the distribution of galaxies; the shorter the energy loss lengths, the stronger the expected anisotropies. Galactic and intergalactic magnetic fields can blur and distort the picture, but the magnitudes of the largest-scale anisotropies, namely the dipole and quadrupole moments, are the most robust to their effects. Measuring them with no bias regardless of any higher-order multipoles is not possible except with full-sky coverage. In this work, we achieve this in three energy ranges (approximately 8--16 EeV, 16--32 EeV, and 32--$\infty$ EeV) by combining surface-detector data collected at the Pierre Auger Observatory until 2020 and at the Telescope Array (TA) until 2019, before the completion of the upgrades of the arrays with new scintillator detectors. We find that the full-sky coverage achieved by combining Auger and TA data reduces the uncertainties on the north-south components of the dipole and quadrupole in half compared to Auger-only results.

Two transiting planet candidates with super-Earth radii around the nearby K7--M0 dwarf star TOI-1238 were announced by TESS. We aim to validate their planetary nature using precise radial velocities (RV) taken with the CARMENES spectrograph. We obtained 55 CARMENES RV data that span 11 months. For a better characterization of the parent star's activity, we also collected contemporaneous optical photometric observations and retrieved archival photometry from the literature. We performed a combined TESS+CARMENES photometric and spectroscopic analysis by including Gaussian processes and Keplerian orbits to account for the stellar activity and planetary signals simultaneously. We estimate that TOI-1238 has a rotation period of 40 $\pm$ 5 d based on photometric and spectroscopic data. The combined analysis confirms the discovery of two transiting planets, TOI-1238 b and c, with orbital periods of $0.764597^{+0.000013}_{-0.000011}$ d and $3.294736^{+0.000034}_{-0.000036}$ d, masses of 3.76$^{+1.15}_{-1.07}$ M$_{\oplus}$ and 8.32$^{+1.90}_{-1.88}$ M$_{\oplus}$, and radii of $1.21^{+0.11}_{-0.10}$ R$_{\oplus}$ and $2.11^{+0.14}_{-0.14}$ R$_{\oplus}$. They orbit their parent star at semimajor axes of 0.0137$\pm$0.0004 au and 0.036$\pm$0.001 au, respectively. The two planets are placed on opposite sides of the radius valley for M dwarfs and lie between the star and the inner border of TOI-1238's habitable zone. The inner super-Earth TOI-1238 b is one of the densest ultra-short-period planets ever discovered ($\rho=11.7^{+4.2}_{-3.4}$ g $\rm cm^{-3}$). The CARMENES data also reveal the presence of an outer, non-transiting, more massive companion with an orbital period and radial velocity amplitude of $\geq$600 d and $\geq$70 m s$^{-1}$, which implies a likely mass of $M \geq 2 \sqrt{1-e^2}$ M$_{\rm Jup}$ and a separation $\geq$1.1 au from its parent star.

We quantify the relative importance of environmental quenching versus pre-processing in $z\sim1$ clusters by analysing the infalling galaxy population in the outskirts of 15 galaxy clusters at $0.8<z<1.4$ drawn from the GOGREEN and GCLASS surveys. We find significant differences between the infalling galaxies and a control sample; in particular, an excess of massive quiescent galaxies in the infalling region. These massive infalling galaxies likely reside in larger dark matter haloes than similar-mass control galaxies because they have twice as many satellite galaxies. Furthermore, these satellite galaxies are distributed in an NFW profile with a larger scale radius compared to the satellites of the control galaxies. Based on these findings, we conclude that it may not be appropriate to use 'field' galaxies as a substitute for infalling pre-cluster galaxies when calculating the efficiency and mass dependency of environmental quenching in high redshift clusters. By comparing the quiescent fraction of infalling galaxies at $1<R/R_{200}<3$ to the cluster sample ($R/R_{200}<1$) we find that almost all quiescent galaxies with masses $>10^{11}M_{\odot}$ were quenched prior to infall, whilst up to half of lower mass galaxies were environmentally quenched after passing the virial radius. This means most of the massive quiescent galaxies in $z\sim1$ clusters were self-quenched or pre-processed prior to infall.

We report new photometric and spectroscopic observations of the K2-99 planetary system. Asteroseismic analysis of the short-cadence light curve from K2's Campaign 17 allows us to refine the stellar properties. We find K2-99 to be significantly smaller than previously thought, with $R_{\star} = 2.55\pm0.02$ $\mathrm{R_\odot}$. The new light curve also contains four transits of K2-99b, which we use to improve our knowledge of the planetary properties. We find the planet to be a non-inflated warm Jupiter, with $R_\mathrm{b} = 1.06 \pm 0.01$ $\mathrm{R_{Jup}}$. Sixty new radial velocity measurements from HARPS, HARPS-N, and HIRES enable the determination of the orbital parameters of K2-99c, which were previously poorly constrained. We find that this outer planet has a minimum mass $M_\mathrm{c} \sin i_\mathrm{c} = 8.4\pm0.2$ $\mathrm{M_{Jup}}$, and an eccentric orbit ($e_\mathrm{c} = 0.210 \pm 0.009$) with a period of $522.2\pm1.4$ d. Upcoming TESS observations in 2022 have a good chance of detecting the transit of this planet, if the mutual inclination between the two planetary orbits is small.

We show that a scalar excited state with large occupation numbers during inflation leads to an enhancement of tensor modes and a characteristic pattern of order-one oscillations in the associated stochastic gravitational wave background (SGWB) sourced during inflation. An effective excited state, i.e. a departure from the Bunch-Davies vacuum, can emerge dynamically as the result of a transient non-adiabatic evolution, e.g. a sharp feature along the inflationary history. We provide an explicit example in a multifield context where the sharp feature triggering the excited state is identified with a strong turn in the inflationary trajectory. En passant, we derive a universal expression for the tensor power spectrum sourced at second order by an arbitrary number of scalar degrees of freedom during inflation, crucially taking into account the nontrivial structure of the Hilbert space in multifield setups. The SGWB sourced during inflation can overcome the standard scalar-induced SGWB sourced at horizon re-entry of the fluctuations after inflation, while being less constrained by perturbativity and backreaction bounds. In addition, one may entertain the possibility of detecting both since they peak at different frequencies exhibiting oscillations with distinct periods.

In this paper we use the model of extragalactic background light to investigate the factors that have influence on the confusion noise. It was shown that (1) Large-Scale Structure of the Universe is an important factor; (2) gravitational lensing does not have a significant effect on the confusion noise; (3) lower redshift limit of objects that contribute to the confusion noise does not depend on the wavelength and is about $z_{min}\sim 0.5-0.6$, while upper redshift limit gradually changes from $\sim4$ to $\sim3$ with the increase of wavelength from 70$\mu m$ up to 2000$\mu m$; (4) at rather short wavelengths ($\simeq70\mu m$) galaxies with luminosities in the range $10^7L_\odot$ -- $10^9L_\odot$ give the most contribution to the confusion noise, while at larger wavelengths (650-2000$\mu m$) their luminosities are greater than $L\geq10^{10}L_\odot$; (5) contribution from objects with different color characteristics is considered; (6) the variability of the extragalactic background on the timescale from 1 day to 1 year is noticeable at short wavelengths (70--350$\mu m$) and manifests at fluxes ${}^<_\sim$ 1~mJy.

Using Hubble Space Telescope imaging of the resolved stellar population of KK~242 = NGC6503-d1 = PGC~4689184, we measure the distance to the galaxy to be $6.46\pm0.32$ Mpc and find that KK~242 is a satellite of the low-mass spiral galaxy NGC~6503 located on the edge of the Local Void. Observations with the Karl G. Jansky Very Large Array show signs of a very faint HI-signal at the position of KK~242 within a velocity range of $V_{hel} = -80\pm10$ km\,s$^{-1}$. This velocity range is severely contaminated by HI emission from the Milky Way and from NGC6503. The dwarf galaxy is classified as the transition type, dIrr/dSph, with a total HI-mass of $< 10^6 M_{\odot}$ and a star formation rate SFR(H$\alpha$) = --4.82 dex ($M{\odot}$/yr). Being at a projected separation of 31 kpc with a radial velocity difference of -- 105 km\,s$^{-1}$ relative to NGC~6503, KK~242 gives an estimate of the halo mass of the spiral galaxy to be $\log(M/M_{\odot}$) = 11.6. Besides NGC~6503, there are 8 more detached low-luminosity spiral galaxies in the Local Volume: M~33, NGC~2403, NGC~7793, NGC~1313, NGC~4236, NGC~5068, NGC~4656 and NGC~7640, from whose small satellites we have estimated the average total mass of the host galaxies and their average total mass-to-K-band-luminosity $\langle M_T/M_{\odot}\rangle = (3.46\pm0.84)\times 10^{11}$ and $(58\pm19) M_{\odot}/L_{\odot}$, respectively.

A reduced-speed-of-light (RSOL) approximation is usually adopted in magnetohydrodynamic (MHD)-particle-in-cell (PIC) simulations, in which relativistic cosmic ray (CR) particles moving at nearly the speed of light $c$ and background (non-relativistic) plasma are evolved concurrently. With a RSOL, some "in-code" speed-of-light $\tilde{c}$ is set to much lower values than the true $c$, allowing simulations to take larger timesteps (which are restricted by the Courant condition given the large CR speeds). However, due to the absence of a well-formulated RSOL implementation in the literature for MHD-PIC simulations, the CR properties in simulations (e.g.\ CR energy or momentum density, gyro radius) vary artificially with respect to each other and with respect to the converged ($\tilde{c} \rightarrow c$) solutions with different implementations of a RSOL, breaking the correspondence between simulation results and physical reality. Here, we derive a new formulation of the MHD-PIC equations with a RSOL, and show that (1) it guarantees all steady-state properties of the CR distribution function and background plasma/MHD quantities are independent of the RSOL $\tilde{c}$ even for $\tilde{c} \ll c$, (2) ensures that the simulation can simultaneously represent the real physical values of CR number, mass, momentum, and energy density, (3) retains the correct physical meaning of various terms like the electric field, and (4) ensures the numerical timestep for CRs can always be safely increased by a factor $\sim c/\tilde{c}$. This should enable greater self-consistency and reduced CPU cost in simulations of CR-MHD interactions.

The different elemental abundances of the photosphere and the corona are striking features of not only the Sun, but other stars as well. This phenomenon is known as the FIP effect (FIP stands for first ionization potential), and its strength can be characterized by the FIP bias, the logarithmic abundance difference between low- and high-FIP elements in the corona, compared to the photosphere. The FIP bias was shown to depend on the surface temperature of the star. We compiled FIP bias and other parameters for 59 stars for which coronal composition is available, now including evolved stars. Using principal component analysis and linear discriminant analysis, we searched for correlations with other astrophysical parameters within the sample which may influence the stellar FIP bias. Adding stars to the $T_{\rm eff}-$FIP bias diagram unveiled new features in its structure. In addition to the previously known relationship, there appears to be a second branch, a parallel sequence about 0.5 dex above it. While the $T_{\rm eff}$ remains the main determinant of the FIP bias, other parameters such as stellar activity indicators also have influence. We find three clusters in the FIP bias determinant parameter space. One distinct group is formed by the evolved stars. Two groups contain main sequence stars in continuation separated roughly by the sign change of the FIP-bias value. The new branch of the $T_{\rm eff}-$FIP bias diagram contains stars with higher activity level, in terms of X-ray flux and rotational velocity. The two main sequence clusters run from the earliest spectral types of A-F with shallow convection zones through G-K-early M stars with gradually deeper convection zones, and end up with the fully convective M dwarf stars, depicting the change of the dynamo type with the internal differences of the main sequence stars in connection with the FIP-bias values.

The BICEP/Keck series of experiments target the Cosmic Microwave Background at degree-scale resolution from the South Pole. Over the next few years, the "Stage-3" BICEP Array (BA) telescope will improve the program's frequency coverage and sensitivity to primordial B-mode polarization by an order of magnitude. The first receiver in the array, BA1, began observing at 30/40 GHz in early 2020. The next two receivers, BA2 and BA3, are currently being assembled and will map the southern sky at frequencies ranging from 95 GHz to 150 GHz. Common to all BA receivers is a refractive, on-axis, cryogenic optical design that focuses microwave radiation onto a focal plane populated with antenna-coupled bolometers. High-performance antireflective coatings up to 760 mm in aperture are needed for each element in the optical chain, and must withstand repeated thermal cycles down to 4 K. Here we present the design and fabrication of the 30/40 GHz anti-reflection coatings for the recently deployed BA1 receiver, then discuss laboratory measurements of their reflectance. We review the lamination method for these single- and dual-layer plastic coatings with indices matched to various polyethylene, nylon and alumina optics. We also describe ongoing efforts to optimize coatings for the next BA cryostats, which may inform technological choices for future Small-Aperture Telescopes of the CMB "Stage 4" experiment.

The BICEP/Keck Collaboration is currently leading the quest to the highest sensitivity measurements of the polarized CMB anisotropies on degree scale with a series of cryogenic telescopes, of which BICEP Array is the latest Stage-3 upgrade with a total of $\sim32,000$ detectors. The instrument comprises 4 receivers spanning 30 to 270 GHz, with the low-frequency 30/40 GHz deployed to the South Pole Station in late 2019. The full complement of receivers is forecast to set the most stringent constraints on the tensor to scalar ratio $r$. Building on these advances, the overarching small-aperture telescope concept is already being used as the reference for further Stage-4 experiment design. In this paper I will present the development of the BICEP Array 150 GHz detector module and its fabrication requirements, with highlights on the high-density time division multiplexing (TDM) design of the cryogenic circuit boards. The low-impedance wiring required between the detectors and the first-stage SQUID amplifiers is crucial to maintain a stiff voltage bias on the detectors. A novel multi-layer FR4 Printed Circuit Board (PCB) with superconducting traces, capable of reading out up to 648 detectors, is presented along with its validation tests. I will also describe an ultra-high density TDM detector module we developed for a CMB-S4-like experiment that allows up to 1,920 detectors to be read out. TDM has been chosen as the detector readout technology for the Cosmic Microwave Background Stage-4 (CMB-S4) experiment based on its proven low-noise performance, predictable costs and overall maturity of the architecture. The heritage for TDM is rooted in mm- and submm-wave experiments dating back 20 years and has since evolved to support a multiplexing factor of 64x in Stage-3 experiments.

What the progenitors of Type Ia supernovae (SNe Ia) are, whether they are near-Chandrasekhar mass or sub-Chandrasekhar mass white dwarfs, has been the matter of debate for decades. Various observational hints are supporting both models as the main progenitors. In this paper, we review the explosion physics and the chemical abundance patterns of SNe Ia from these two classes of progenitors. We will discuss how the observational data of SNe Ia, their remnants, the Milky Way Galaxy, and galactic clusters can help us to determine the essential features where numerical models of SNe Ia need to match.

The H.E.S.S. Galactic Plane Survey has detected very-high-energy (VHE) gamma-ray emission from 78 sources in the Milky Way. These sources belong to different object classes (pulsar wind nebulae, supernova remnants or binary systems) and some of these sources remain unidentified. The gamma-ray emission of these objects may be of leptonic or hadronic origin and gamma-ray observations alone cannot distinguish between these two scenarios. The detection of neutrino emission would provide evidence for a hadronic scenario in these objects. Based on the observed gamma-ray spectra we predict the neutrino emission of these sources under the hypothesis that the emission is solely of hadronic origin. This prediction relies entirely on observation and is independent of the source class, the distance or the ambient target material. We use these predictions to create an empirical model for the neutrino emission of the Milky Way. This model can be used to search for neutrino emission from individual gamma-ray sources as well as testing for neutrino emission from potential source populations in the Milky Way.

Dark matter (DM) freeze-in through a light mediator is an appealing model with excellent detection prospects at current and future experiments. Light mediator freeze-in is UV-insensitive insofar as most DM is produced at late times, and thus the DM abundance does not depend on the unknown early evolution of our universe. However the final DM yield retains a dependence on the initial DM population, which is usually assumed to be exactly zero. We point out that in models with light mediators, the final DM yield will also depend on the initial conditions assumed for the light mediator population. We describe a class of scenarios we call "glaciation" where DM freezing in from the SM encounters a pre-existing thermal bath of mediators, and study the dependence of the final DM yield on the initial temperature of this dark radiation bath. To compute DM scattering rates in this cosmology, we derive for the first time an exact integral expression for the Boltzmann collision term describing interactions between two species at different temperatures. We quantify the dependence of the DM yield on the initial dark temperature and find that it can be sizeable in regions near the traditional (zero initial abundance) freeze-in curve. We generalize the freeze-in curve to a glaciation band, which can extend as much as an order of magnitude below the traditional freeze-in direct detection target, and point out that the DM phase space distribution as well as the yield can be strongly dependent on initial conditions.

The successful transition from core-collapse supernova simulations using classical neutrino transport to simulations using quantum neutrino transport will require the development of methods for calculating neutrino flavor transformations that mitigate the computational expense. One potential approach is the use of angular moments of the neutrino field, which has the added appeal that there already exist simulation codes which make use of moments for classical neutrino transport. Evolution equations for quantum moments based on the quantum kinetic equations can be straightforwardly generalized from the evolution of classical moments based on the Boltzmann equation. We present an efficient implementation of neutrino transformation using quantum angular moments in the free streaming, spherically symmetric bulb model. We compare the results against analytic solutions and the results from more exact multi-angle neutrino flavor evolution calculations. We find that our moment-based methods employing scalar closures predict, with good accuracy, the onset of collective flavor transformations seen in the multi-angle results. However, they over-estimate the coherence between neutrinos traveling along different trajectories in tests that neglect neutrino-neutrino interactions. More sophisticated quantum closures may improve the agreement between the inexpensive moment-based methods and the multi-angle approach.

Using a quantum hadrodynamics (QHD) and MIT based models we construct hybrid stars within the Maxwell criteria of hadron-quark phase transition. We are able to produce a hybrid star with maximum mass of 2.15$M_\odot$. Furthermore, a 2.03$M_\odot$ star with a quark core corresponding to more than $80\%$ of both, its total mass and radius, is also possible.

Gravitational waves from binary black hole merger leave permanent imprint on space-time, called gravitational wave memory. Amounts of memory events can form a ``stochastic gravitational wave memory background'' (SGWMB) which has not been investigated before. Here we find that SGWMB can be described as a Brownian motion. The power spectral density (PSD) of SGWMB corresponding to a set of binary black hole coalescence is proportional to $1/f^2$ where $f$ is the frequency of the background. And the strength of the SGWMB PSD depends on the event rate of the binary black hole coalescence. Not like the stochastic non-memory gravitational wave background of a set of binary black hole coalescence whose PSD has a upper frequency bound due to the merger frequency, SGWMB can expand to much higher frequency because of the Brownian motion behavior. Consequently SGWMB dominates the stochastic gravitational wave background in the relatively higher frequency band. If not counting SGWMB properly ones may confuse it with the stochastic non-memory gravitational wave background. Investigating SGWMB properly may let LIGO to detect the population of supermassive binary black hole systems. Also SGWMB provides a brand new mean to detect gravitational wave memory and it opens a new window to explore the gravity theory.

The time evolution of the field perturbations in the wormhole and black bounce backgrounds are investigated. We find that the asymmetry of spacetime results in the asymmetry of the effective potential of the perturbed equation. The quasinormal modes are strongly dependent on the shapes of the effective potentials. Specifically, the signals of echoes arise in some wormhole cases and reflect the asymmetric properties of wormholes. We examine the features of echoes within different circumstances. Besides, the negative values of effective potentials usually imply the instability of the system. By analyzing some specific metrics, we find that the negative regions of effective potentials are enclosed by the black hole horizons in these cases. But this statement could be broken in asymmetric cases.

We measured the relative positions between two pairs of compact extragalactic sources (CESs), J1925-2219 \& J1923-2104 (C1--C2) and J1925-2219 \& J1928-2035 (C1--C3) on 2020 October 23--25 and 2021 February 5 (totaling four epochs), respectively, using the Very Long Baseline Array (VLBA) at 15 GHz. Accounting for the deflection angle dominated by Jupiter, as well as the contributions from the Sun, planets other than Earth, the Moon and Ganymede (the most massive of the solar system's moons), our theoretical calculations predict that the dynamical ranges of the relative positions across four epochs in R.A. of the C1--C2 pair and C1--C3 pair are 841.2 and 1127.9 $\mu$as, respectively. The formal accuracy in R.A. is about 20 $\mu$as, but the error in Decl. is poor. The measured standard deviations of the relative positions across the four epochs are 51.0 and 29.7 $\mu$as in R.A. for C1--C2 and C1--C3, respectively. These values indicate that the accuracy of the post-Newtonian relativistic parameter, $\gamma$, is $\sim 0.061$ for C1--C2 and $\sim 0.026$ for C1--C3. Combining the two CES pairs, the measured value of $\gamma$ is $0.984 \pm 0.037$, which is comparable to the latest published results for Jupiter as a gravitational lens reported by Fomalont \& Kopeikin, i.e., $1.01 \pm 0.03$.

In the recent work [Dailey et al., Nature Astronomy 5, 150 (2021)], it was claimed that networks of quantum sensors can be used as sensitive multi-messenger probes of astrophysical phenomena that produce intense bursts of relativistic bosonic waves which interact non-gravitationally with ordinary matter. The most promising possibility considered in [Ibid.] involved clock-based searches for quadratic scalar-type interactions, with greatly diminished reach in the case of magnetometer-based searches for derivative-pseudoscalar-type interactions and clock-based searches for linear scalar-type interactions. In this note, we point out that the aforementioned work overlooked the ''back action'' of ordinary matter on scalar waves with quadratic interactions and that accounting for back-action effects can drastically affect the detection prospects of clock networks. In particular, back action can cause strong screening of scalar waves near Earth's surface and by the apparatus itself, rendering clock experiments insensitive to extraterrestrial sources of relativistic scalar waves. Additionally, back-action effects can retard the propagation of scalar waves through the interstellar and intergalactic media, significantly delaying the arrival of scalar waves at Earth compared to their gravitational-wave counterparts and thereby preventing multi-messenger astronomy on human timescales.

Sensor fusion is a technique used to combine sensors with different noise characteristics into a super sensor that has superior noise performance. To achieve sensor fusion, complementary filters are used in current gravitational-wave detectors to combine relative displacement sensors and inertial sensors for active seismic isolation. Complementary filters are a set of digital filters, which have transfer functions that are summed to unity. Currently, complementary filters are shaped and tuned manually rather than optimized, which can be suboptimal and hard to reproduce for future detectors. In this paper, an optimization-based method called $\mathcal{H}_\infty$ synthesis is proposed for synthesizing optimal complementary filters according to the sensor noises themselves. The complementary filter design problem is converted into an optimization problem that seeks minimization of an objective function equivalent to the maximum difference between the super sensor noise and the lower bound in logarithmic scale. The method is exemplified by synthesizing complementary filters for sensor fusion of 1) a relative displacement sensor and an inertial sensor, 2) a relative displacement sensor coupled with seismic noise and an inertial sensor, and 3) hypothetical displacement sensor and inertial sensor, which have slightly different noise characteristics compared to the typical ones. In all cases, the method produces complementary filters that suppress the super sensor noise equally close to the lower bound at all frequencies in logarithmic scale. The synthesized filters contain features that better suppress the sensor noises compared to the pre-designed complementary filters. Overall, the proposed method allows the synthesis of optimal complementary filters according to the sensor noises themselves and is a better and versatile method for solving sensor fusion problems.

We establish the nonlocal generalization of the Israel-Stewart model for the relativistic causal thermodynamics of the cosmic fluid, which evolves in the homogeneous isotropic Universe. Based on the second law of thermodynamics we derive the integro-differential master equation for the nonequilibrium pressure scalar and reduce it to the differential equation of the second order in time derivatives. We show that this master equation can be considered as the relativistic analog of the Burgers equation, which is known in the classical theory of viscoelasticity. We obtain the nonlinear key equation, which contains the energy density scalar only, and analyze two exact solutions of the model. The effective temperature is considered to be associated with the barotropic equation of state of the cosmic fluid.

Phenomenological success of inflation models with axion and SU(2) gauge fields relies crucially on control of backreaction from particle production. Most of the previous study only demanded the backreaction terms in equations of motion for axion and gauge fields be small on the basis of order-of-magnitude estimation. In this paper, we solve the equations of motion with backreaction for a wide range of parameters of the spectator axion-SU(2) model. First, we find a new slow-roll solution of the axion-SU(2) system in the absence of backreaction. Next, we obtain accurate conditions for stable slow-roll solutions in the presence of backreaction. Finally, we show that the amplitude of primordial gravitational waves sourced by the gauge fields can exceed that of quantum vacuum fluctuations in spacetime by a large factor, without backreaction spoiling slow-roll dynamics. Imposing additional constraints on the power spectra of scalar and tensor modes measured at CMB scales, we find that the sourced contribution can be more than ten times the vacuum one. Imposing further a constraint of scalar modes non-linearly sourced by tensor modes, the two contributions can still be comparable.

The search for neutrino events in correlation with several of the most intense fast radio bursts (FRBs) has been performed using the Borexino data. We have searched for signals with visible energies above $250$~keV within a time window of $\pm$1000~s corresponding to the detection time of a particular FRB. We also applied an alternative approach based on searching for specific shapes of neutrino-electron scattering spectra in the full exposure spectrum of the Borexino detector. In particular, two incoming neutrino spectra were considered: the monoenergetic line and the spectrum expected from supernovae. The same spectra were considered for electron antineutrinos detected through the inverse beta-decay reaction. No statistically significant excess over the background was observed. As a result, the strongest upper limits on FRB-associated neutrino fluences of all flavors have been obtained in the $0.5 - 50$~MeV neutrino energy range.