Papers for Thursday, Dec 07 2023

Models derived in 2009 to fit mid-infrared (8-24 micron) source counts from the IRAS, ISO and Spitzer missions, provide an excellent fit to deep counts with JWST, demonstrating that the evolution of dusty star-forming galaxies is well understood. The evolution of dust in galaxies at high redshifts is discussed and a simple prescription is proposed to model this. This allows more realistic models for source-counts at submillimetre wavelength. A reasonable fit to 250, 500, 850 and 1100 micron counts is obtained. This paper therefore draws together the IRAS, ISO, Spitzer, Akari, Herschel, submillimetre ground-based, and JWST surveys into a single picture.

The theory of how low mass stars form from the collapse of a dense molecular cloud core has been well-established for decades. Thanks to significant progress in computing and numerical modelling, more physical models have been developed and a wider parameter space explored to understand the early stages of star formation more fully. In this review, I describe the expected physical properties of the first and second core stages and how the inclusion of different physics affects those predicted characteristics. I provide an overview of chemical models and synthetic observations, looking towards the positive identification of the first core in nature, which remains elusive. However, there are a few likely candidate first cores, which are listed, and I briefly discuss the recent progress in characterising the youngest protostellar sources. Chemistry will be instrumental in the firm identification of the first core so we require robust theoretical predictions of the chemical evolution of protostellar cores, especially of the first and second core outflows. Looking ahead, simulations can shed light on how the protostellar collapse phase shapes the evolution of the protostellar disc. Simulations of dust evolution during protostellar core collapse show there is significant enhancement in grain size and abundance towards the centre of the core. Chemical models show that the warm, dense conditions of the first core drive chemical evolution. There is a wide scope for further study of the role that the first and second core stages play in determining the structure and composition of the protostellar disc and envelope and, of course, the eventual influence on the formation of planets.

We report on the X-ray spectral and spatial evolution of the Symbiotic star R Aqr. Through a multi-epoch observational campaign performed with Chandra between 2017 and 2022, we study the X-ray emission of this binary system, composed of an evolved red giant star and a white dwarf (WD). This analysis is particularly timely as the WD approached the periastron in late 2018/early 2019, thus mass transfer, jet emission and outburst phenomena are to be expected. Through detailed spectral analysis, we detect a significant rise in the soft X-ray (0.5-2 keV) emission of R Aqr, likely linked to jet emission, followed by a decay towards the previous quiescent state. The hard X-ray emission (5-8 keV), is not immediately affected by the periastron passage; the hard component, after maintaining the same flux level between 2017 and 2021, rapidly decays after 2022. Possible explanations for this are a change in the reflection properties of the medium surrounding the binary, obscuration of the central region by material ejected during the periastron passage, or even the partial/complete destruction of the inner regions of the accretion disc surrounding the WD. In addition to this activity in the central region, extended emission is also detected, likely linked to a hot spot in a pre-outburst-emitted jet, which can be observed moving away from the system's central region.

Pulsar-timing collaborations have recently reported evidence for the detection of an isotropic stochastic gravitational-wave background consistent with one sourced by a population of inspiralling supermassive black hole binaries. However, a certain degree of anisotropy and polarization may be present. Thus, the characterization of the energy density and polarization of the background at different angular scales is important. In this paper, we describe the signatures of linear polarization in the stochastic gravitational-wave background on the timing residuals obtained with pulsar-timing arrays. We expand the linear polarization map in terms of spin-weighted spherical harmonics and recast it into the $E$-mode (parity even) and $B$-mode (parity odd) basis. We provide expressions for the minimum-variance estimators for the coefficients of that expansion and evaluate the smallest detectable signal as a function of the signal-to-noise ratio with which the isotropic GW signal is detected and the number of pulsars in the survey. We evaluate the covariance between the estimators for the spherical-harmonic coefficients of the linear polarization $E$-modes and those for the intensity anisotropy. We also show that there is no covariance between the spherical-harmonic coefficients for the $B$-modes of the linear polarization and those for the circular polarization, even though both have the same parity. Our approach results in simple, elegant, and easily evaluated expressions for the overlap reduction functions for linear polarization.

We present a new physically-motivated model for estimating the molecular line emission in active galaxies. The model takes into account (i) the internal density structure of giant molecular clouds (GMCs), (ii) the heating associated both to stars and to the active galactic nuclei (AGN), respectively producing photodissociation regions (PDRs) and X-ray dominated regions (XDRs) within the GMCs, and (iii) the mass distribution of GMCs within the galaxy volume. The model needs, as input parameters, the radial profiles of molecular mass, far-UV flux and X-ray flux for a given galaxy, and it has two free parameters: the CO-to-H2 conversion factor $\alpha_{CO}$, and the X-ray attenuation column density $N_H$. We test this model on a sample of 24 local ($z \leq 0.06$) AGN-host galaxies, simulating their carbon monoxide spectral line energy distribution (CO SLED). We compare the results with the available observations and calculate, for each galaxy, the best ($\alpha_{CO}$, $N_H$) with a Markov chain Monte Carlo algorithm, finding values consistent with those present in the literature. We find a median $\alpha_{CO} = 4.8$ M$_{\odot}$ (K km s$^{-1}$ pc$^{2}$)$^{-1}$ for our sample. In all the modelled galaxies, we find the XDR component of the CO SLED to dominate the CO luminosity from $J_{\text{upp}} \geq 4$. We conclude that, once a detailed distribution of molecular gas density is taken into account, PDR emission at mid-/high-$J$ becomes negligible with respect to XDR.

We present the first star formation history (SFH) and age-metallicity relation (AMR) derived from resolved stellar populations imaged with the JWST NIRCam instrument. The target is the Local Group star-forming galaxy WLM at 970 kpc. The depth of the color-magnitude diagram (CMD) reaches below the oldest main sequence turn-off with a SNR=10 at M_F090W=+4.6 mag; this is the deepest CMD for any galaxy that is not a satellite of the Milky Way. We use Hubble Space Telescope (HST) optical imaging that overlaps with the NIRCam observations to directly evaluate the SFHs derived based on data from the two great observatories. The JWST and HST-based SFHs are in excellent agreement. We use the metallicity distribution function measured from stellar spectra to confirm the trends in the AMRs based on the JWST data. Together, these results confirm the efficacy of recovering a SFH and AMR with the NIRCam F090W-F150W filter combination and provide validation of the sensitivity and accuracy of stellar evolution libraries in the near-infrared relative to the optical for SFH recovery work. From the JWST data, WLM shows an early onset to star formation, followed by an extended pause post-reionization before star formation re-ignites, which is qualitatively similar to what has been observed in the isolated galaxies Leo~A and Aquarius. Quantitatively, 15% of the stellar mass formed in the first Gyr, while only 10% formed over the next ~5 Gyr; the stellar mass then rapidly doubled in ~2.5 Gyr, followed by constant star formation over the last ~5 Gyr.

Astronomical images often have regions with missing or unwanted information, such as bad pixels, bad columns, cosmic rays, masked objects, or residuals from imperfect model subtractions. In certain situations it can be essential, or preferable, to fill in these regions. Most existing methods use low order interpolations for this task. In this paper a method is described that uses the full information that is contained in the pixels just outside masked regions. These edge pixels are extrapolated inwards, using iterative median filtering. This leads to a smoothly varying spatial resolution within the filled-in regions, and ensures seamless transitions between masked pixels and good pixels. Gaps in continuous, narrow features can be reconstructed with high fidelity, even if they are large. The method is implemented in maskfill, an open-source MIT licensed Python script. Its performance is illustrated with several examples.

Recently, several non-interacting black hole-stellar binaries have been identified in Gaia data. For example, Gaia BH1, where a Sun-like star is in moderate eccentricity (e=0.45), 185-day orbit around a black hole. This orbit is difficult to explain through binary evolution. The present-day separation suggests the progenitor system would have undergone an episode of common envelope evolution, but a common envelope should shrink the period below the observed one. Since the majority of massive stars form in higher multiplicity systems, a triple evolution scenario is more likely for the progenitors of BH binaries. Here we show that such systems can indeed be more easily explained via evolution in hierarchical triple systems. von Zeipel--Lidov--Kozai oscillations or instabilities can delay the onset of the common envelope phase in the inner binary of the triple, so that the black hole progenitor and low-mass star are more widely separated when it begins, leading to the formation of wider binaries. There are also systems with similar periods but larger eccentricities, where the BH progenitor is a merger product of the inner binary in the triple. Such mergers lead to a more top-heavy black hole mass function.

We study the properties of line bisectors in the spectrum of the Sun-as-a-star, as observed using the Integrated Sunlight Spectrometer (ISS) of the SOLIS project. Our motivation is to determine whether changes in line shape, due to magnetic modulation of photospheric convection can be separated from the 9 cm/sec Doppler reflex of the Earth's orbit. Measuring bisectors of 21 lines over a full solar cycle, our results overwhelmingly indicate that solar magnetic activity modulates photospheric convection so as to reduce the asymmetries of line profiles in the spectrum of the Sun-as-a-star (having both C-shaped and reversed-C-shaped bisectors). However, some lines are constant or have variations in shape that are too small to measure. We inject a 9 cm/sec radial velocity signal with a 1-year period into the ISS spectra. Informed by a principal component analysis of the line bisectors, we fit the most significant components to the bisectors of each line by linear regression, including a zero-point offset in velocity that is intended to capture the injected radial velocity signal. Averaging over lines, we are able to recover that signal to solid statistical significance in the presence of much larger changes in the line shapes. Although our work has limitations (that we discuss), we establish that changes in absorption line shapes do not in themselves prevent the detection of an Earth-like planet orbiting a Sun-like star using precise radal velocity techniques.

In the canonical tachyonic resonance preheating scenario, only an order one fraction of energy density in the inflaton is transferred to radiation, due to backreaction effects. One possible way to improve the energy transfer efficiency is to allow for the perturbative decays of the resonantly produced daughter particles, which serve as the "spillway" to drain the direct decay products from inflaton and to reduce the backreaction. In this article, we study two observational consequences of spillway preheating. The first is on the inflationary observables: the scalar spectrum tilt $n_s$ and tensor-to-scalar ratio $r$. The spillway scenario modifies the evolution of the equation of state between the end of inflation and the thermal big bang. As a result, it affects the time elapsed from inflation to the Cosmic Microwave Background (CMB), as well as the fits of inflationary models and their corresponding prediction for $n_s$ and $r$. We map out the equation of state by systematically scanning the parameter space of the spillway scenario, and show that the most efficient spillway scenario predicts a bluer spectrum, compared to the tachyonic preheating scenario. Another consequence is the production of high-frequency gravitational waves (GWs). Comparing the simulation results with those of tachyonic preheating, we find that the existence of spillways leads to sharper-peaked GW spectra with a mildly damped amplitude.

Dusty star-forming galaxies emit most of their light at far-IR to mm wavelengths as their star formation is highly obscured. Far-IR and mm observations have revealed their dust, neutral and molecular gas properties. The sensitivity of JWST at rest-frame optical and near-infrared wavelengths now allows the study of the stellar and ionized gas content. We investigate the spatially resolved distribution and kinematics of the ionized gas in GN20, a dusty star forming galaxy at $z$=4.0548. We present deep MIRI/MRS integral field spectroscopy of the near-infrared rest-frame emission of GN20. We detect spatially resolved \paa, out to a radius of 6 kpc, distributed in a clumpy morphology. The star formation rate derived from \paa\ (144 $\pm$ 9 \msunperyear) is only 7.7 $\pm 0.5 $\% of the infrared star formation rate (1860 $\pm$ 90 \msunperyear). We attribute this to very high extinction (A$_V$ = 17.2 $\pm$ 0.4 mag, or A$_{V,mixed}$ = 44 $\pm$ 3 mag), especially in the nucleus of GN20, where only faint \paa\ is detected, suggesting a deeply buried starburst. We identify four, spatially unresolved, clumps in the \paa\ emission. Based on the double peaked \paa\ profile we find that each clump consist of at least two sub-clumps. We find mass upper limits consistent with them being formed in a gravitationally unstable gaseous disk. The UV bright region of GN20 does not have any detected \paa\ emission, suggesting an age of more than 10 Myrs for this region of the galaxy. From the rotation profile of \paa\ we conclude that the gas kinematics are rotationally dominated and the $v_{rot}/\sigma_{m} = 3.8 \pm 1.4$ is similar to low-redshift LIRGs. We speculate that the clumps seen in GN20 could contribute to building up the inner disk and bulge of GN20.

We present a new five-epoch Chandra X-ray Observatory monitoring survey of the nearby spiral galaxy M33 which probes X-ray variability with time sampling between two weeks and four months. We characterize the X-ray variability of 55 bright point sources outside of the nucleus, many of which are expected to be high-mass X-ray binaries (HMXBs). We detect eight new candidate transients not detected in previous X-ray catalogs of M33 and discuss their possible nature. The final catalog includes 26 known HMXB candidates identified in the literature. We extend the baseline of the X-ray light curves up to 21 years by including archival X-ray observations of these sources. We compare the detection and non-detection epochs of the sources to suites of simulated source duty cycles and infer that most of our detected sources have duty cycles > 30%. We find only four sources whose detection patterns are consistent with having duty cycles below 30%. This large fraction of sources with high duty cycles is unexpected for a population of HMXBs, thus more frequent X-ray monitoring will likely reveal many more low duty cycle HMXBs in M33.

We report new analyses of spectra of the $3.2-3.3~\mu$m absorption feature observed in the diffuse interstellar medium toward three Milky Way sources: 2MASS $J17470898-2829561$ (2M1747) and the Quintuplet Cluster, both located in the Galactic center, and Cygnus OB2-12. The $3.2-3.3~\mu$m interval coincides with the CH-stretching region for compact polycyclic aromatic hydrocarbons (PAHs). We focus on the 2M1747 spectrum. Its published optical depth spectrum contains residual telluric transmission features, which arise from the 0.06 difference in mean airmasses between the observations of the source and its telluric standard star. We corrected the published spectrum by adding the airmass residual optical depth spectrum. The corrected spectrum is well fit by a superposition of four Gaussians. The absorption spectra of the other two sources were also fit by four Gaussians, with similar central wavelengths, widths, and relative peak opacities. We associate the three longer wavelength Gaussians covering the $3.23-3.31~\mu$m interval with compact PAHs in positive, neutral, and negative charge states. We identify the shortest wavelength Gaussian, near 3.21 $\mu$m, with irregularly-shaped PAHs. Constraints imposed by spectral smoothness on the corrected 2M1747 spectrum, augmented by a PAH cluster formation model for post-asymptotic giant branch stars, suggests that $> 99$\%\ of the PAHs in the diffuse interstellar medium reside in small clusters. This study supports the PAH hypothesis, and suggests that a family of primarily compact PAHs with a C$_{66}$H$_{20}$ (circumvalene) parent is consistent with the observed mid-infrared and ultraviolet interstellar absorption spectrum.

The origin of repeating fast radio bursts (RFRBs) is still a mystery. We propose that short-lived RFRBs might be triggered from the tidal disruption of white dwarfs (WDs) by intermediate-mass black holes (BHs). In this model, we show that the remnant WD clusters after tidal collapse cuts the magnetic lines on the BH accretion discs, and during each fall of the clump, so that electrons are torn from the surface of the mass and instantly accelerated to the relativistic energy. The subsequent movement of these electrons along magnetic field lines will result in coherent curvature radiation. This short-lived radio transients might accompany with the accretion process. The luminosity and the timescale can be estimated to be $L_\mathrm{tot}\sim 1.96\times10^{40}~{\rm erg~s^{-1}}$ and $\Delta t\sim1.14~{\rm ms}$, respectively, which are consistent with the typical properties of RFRBs. Moreover, the total event rate of our model for generating RFRBs might be as high as $\sim 10~\rm {yr^{-1}~Gpc^{-3}}$.

The Blandford-Znajek (BZ) mechanism in stellar-mass black hole (BH) hyperaccretion systems is generally considered to power gamma-ray bursts (GRBs). Based on observational GRB data, we use the BZ mechanism driven by the BH hyperaccretion disc to investigate the evolution of the BH mass and spin after the jets break out from the progenitors. We find that the BH growths are almost independent of initial BH masses. Meanwhile, the BH growths will be more efficient with smaller initial spin parameters. We conclude that (i) the BZ mechanism is efficient for triggering BH growths for only 1 of 206 typical long-duration GRBs; (ii) the mean BH mass growths of ultra-long GRBs are marginal for all 7 samples collected; (iii) for the short-duration GRBs, the results that BHs show minimal growths is consistent with the mass supply limitation in the scenario of compact object mergers.

A tension between the expansion rate of the universe obtained from direct measurements and that inferred from fits to the acoustic power spectrum of the Cosmic Microwave Background (CMB) radiation has emerged as higher multipoles have been incorporated into CMB fits. This temporal evolution is suggestive of a systematic effect at fine angular scales, and may be related to the observation of unexpectedly massive galaxies at high redshift. Such objects could cause anomalous gravitational lensing of the surface of last scattering, which in turn may affect the best-fit value of $H_0$. If so, the Hubble tension may be the result of a systematic effect in the fitting of CMB data rather than a systematic error in local measurements.

We present the extended ALMA MaNGA QUEnching and STar formation survey, a combination of the original 46 ALMaQUEST galaxies plus new ALMA observations for a further 20 interacting galaxies. Three well-studied scaling relations are fit to the 19,999 star-forming spaxels in the extended sample, namely the resolved Schmidt-Kennicutt (rSK) relation, the resolved star forming main sequence (rSFMS) and the resolved molecular gas main sequence (rMGMS). We additionally investigate the relationship between the dynamical equilibrium pressure (PDE) and star formation rate surface density (Sigma_SFR), which we refer to as the resolved PDE (rPDE) relation. Contrary to previous studies that have focussed on normal star-forming galaxies and found an approximately linear rPDE relation, the presence of more vigourously star-forming galaxies in the extended ALMaQUEST sample reveals a marked turnover in the relation at high pressures. Although the scatter around the linear fit to the rPDE relation is similar to the other three relations, a random forest analysis, which can extract non-linear dependences, finds that PDE is unambiguously more important than either Sigma_H2 or Sigma_star for predicting Sigma_SFR. We compare the observed rPDE relation to the prediction of the pressure-regulated feedback-modulated (PRFM) model of star formation, finding that galaxies residing on the global SFMS do indeed closely follow the rPDE relation predicted by the PRFM theory. However, galaxies above and below the global SFMS show significant deviations from the model. Galaxies with high SFR are instead consistent with models that include other contributions to turbulence in addition to the local star formation feedback.

Theta Mus is a remarkable spectroscopic binary (SB) consisting of a carbon-type Wolf-Rayet star and OV companion (WC6+O6-7V) in a 19-day orbit. In addition an O-supergiant is visually detected at a small offset of 46 mas and if gravitationally bound to the SB system would have an orbital period of many decades. Theta Mus is X-ray bright and a nonthermal radio source as commonly observed in massive colliding wind (CW) binaries. We present new Chandra X-ray observations of Theta Mus which complement previous XMM-Newton observations. The X-ray emission consists of a cool nearly steady weakly-absorbed plasma component with broad redshifted emission lines located in an extended region far from the SB system. Hotter plasma is also present traced by Fe XXV emission. The observed flux in the 2-5 keV range dropped significantly on a timescale of less than 5 years. The flux decrease can be attributed to an increase in absorption toward the hotter plasma which is likely located in the confined wind interaction region of the short-period SB system. The X-ray emission of Theta Mus is remarkably similar to the WC+O binary gamma^2 Vel including carbon recombination spectral lines but both systems show unusual line centroid properties that challenge CW models.

This is a model-independent analysis that investigates the statistical isotropy in the Local Universe using the ALFALFA survey data ($0 < z < 0.06$). We investigate the angular distribution of HI extra-galactic sources from the ALFALFA catalogue and study whether they are compatible with the statistical isotropy hypothesis using the two-point angular correlation function (2PACF). Aware that the Local Universe is plenty of clustered structures and large voids, we compute the 2PACF with the Landy-Szalay estimator performing directional analyses to inspect 10 sky regions. We investigate these 2PACF using power-law best-fit analyses, and determine the statistical significance of the best-fit parameters for the 10 ALFALFA regions by comparison with the ones obtained through the same procedure applied to a set of mock catalogues produced under the homogeneity and isotropy hypotheses. Our conclusion is that the Local Universe, as mapped by the HI sources of the ALFALFA survey, is in agreement with the hypothesis of statistical isotropy within $2\,\sigma$ confidence level, for small and large angle analyses, with the only exception of one region --located near the Dipole Repeller-- which appears slightly outlier ($2.4\,\sigma$). Interestingly, regarding the large angular distribution of the HI sources, we found 3 regions where the presence of cosmic voids reported in the literature left their signature in our 2PACF, suggesting projected large underdensities there, with number-density contrast $\delta \simeq -0.7$. According to the current literature these regions correspond, partially, to the sky position of the void structures known as Local Cosmic Void and Dipole Repeller.

Over this last year we have published four independent refereed studies confirming the presence of a gravitational anomaly, from studying the relative velocities and separations on the plane of the sky, $v_{2D}$ and $s_{2D}$ respectively, of wide binary stars observed by the {\it Gaia} satellite. These studies show results which are fully consistent with Newtonian dynamics in a high acceleration $s_{2D}<2000$ au regime, but which conclusively identify MOND phenomenology for the low acceleration $s_{2D}>2000$ au regime. These four studies span a range of sample selection strategies and cover also a range of statistical techniques, in all cases, results are consistent in identifying a change in the effective value of the gravitational constant in the low acceleration $s_{2D}>2000$ au regime of $G \to \gamma G$ with $\gamma=1.5 \pm \sigma_{\gamma}$, with $0.06<\sigma_{\gamma}<0.2$, depending on the sample selection strategy and the statistical modelling implemented. Recently, a contradictory study appeared, Banik et al. (2024) (originally arXiv:2311.03436 in 2023), claiming a 19$\sigma$ statistical preference of a purely Newtonian model over a MOND alternative, looking also at distributions of $v_{2D}$ and $s_{2D}$ for wide binaries from the {\it Gaia} satellite, although only in the $s_{2D}>2000$ au regime. Inspection of said study readily shows the use of a statistical treatment inconsistent with the error structure of the data used. Related to the above, the best fit posterior physical parameters found in Banik et al. (2024), show a lack of correlation with their inferred gravity index. In this brief comment we expand upon the above to show that the results of Banik et al. (2024) are not due to the physics of the problem being treated, but to the methodological problems of that study.

We present the analysis of 200-ks NuSTAR observation of the Vela pulsar and the pulsar wind nebula (PWN). The phase-resolved spectra corresponding to the two main peaks in the folded pulse profile differ significantly. The spectrum of Peak 1 is significantly harder than that of Peak 2 in qualitative agreement with the earlier RXTE results. However, for both spectra, the best-fit values of photon indices for the power-law (PL) fit are noticeably larger than the previously reported values. The hardest (Peak 1) spectrum has a photon index of $1.10\pm0.15$ which is close to those measured for the bright inner jets of the PWN. We used the off-pulse interval to isolate pulsar emission and measure the compact pulsar wind nebula (PWN) spectrum in hard X-rays. We also measured the spectrum from the south-western (SW) region of the PWN which is resolved by NuSTAR from the compact PWN. For both regions, we fit the NuSTAR spectra by themselves and together with the Chandra X-ray Observatory spectra. We found that the compact PWN spectrum requires a more complex model than a simple PL with evidence for exponential cutoff in 50-80 keV. We do not find such evidence for the spectrum extracted from the SW PWN region, located farther from the pulsar, implying the energies over 600 TeV for the pulsar wind particles. This may indicate in-situ particle acceleration in that region.

The origin of the solar system abundances of several proton-rich isotopes, especially $^{92,94}$Mo and $^{96,98}$Ru, has been an enduring mystery in nuclear astrophysics. An attractive proposal to solve this problem is the $\nu p$-process, which can operate in neutrino-driven outflows in a core-collapse supernova after the shock is launched. Years of detailed studies, however, have cast doubt over the ability of this process to generate sufficiently high absolute and relative amounts of various $p$-nuclei. The $\nu p$-process is also thought to be excluded by arguments based on the long-lived radionuclide $^{92}$Nb. Here, we present explicit calculations, in which both the abundance ratios and the absolute yields of the $p$-nuclei up to $A\lesssim105$ are successfully reproduced, even when using the modern (medium enhanced) triple-$\alpha$ reaction rates. The process is also shown to produce the necessary amounts of $^{92}$Nb. The models are characterized by subsonic outflows and by the protoneutron star masses in the {$\gtrsim1.7 M_\odot$ range}. This suggests that the Mo and Ru $p$-nuclides observed in the Solar System were made in CCSN explosions characterized by an extended accretion stage.

Stars approaching supermassive black holes can be tidally disrupted. Despite being expected to emit X-rays, TDEs have been largely observed in optical bands, which is poorly understood. In this Letter, we simulate the tidal disruption of a $1~M_\odot$ main sequence star on an eccentric ($e=0.95$) orbit with a periapsis distance one or five times smaller than the tidal radius ($\beta = 1$ or $5$) using general relativistic smoothed particle hydrodynamics. We follow the simulation for up to a year post-disruption. We show that accretion disks in eccentric TDEs are masked by unbound material outflowing at $\sim10,000~$km/s. Assuming electron scattering opacity, this material would be visible as a $\sim100~$au photosphere at $\sim10^4~$K, in line with observations of candidate TDEs.

HD 148937 is a peculiar massive star (Of?p) with a strong magnetic field (1kG). The hourglass-shaped emission nebula NGC 6164/5 surrounds this star. This nebula is presumed to originate from episodic mass-loss events of the central O-type star, but the detailed formation mechanism is not yet well understood. Grasping its three-dimensional structure is essential to uncover the origin of this nebula. Here we report the high-resolution multi-object spectroscopic observation of NGC 6164/5 using the GIRAFFE on the 8.2m Very Large Telescope. Integrated intensity maps constructed from several spectral lines delineate well the overall shape of this nebula, such as the two bright lobes and the inner gas region. The position-velocity diagrams show that the two bright lobes are found to be redshifted and blueshifted, respectively, while the inner region has multiple layers. We consider a geometric model composed of a bilateral outflow harboring nitrogen-enriched knots and expanding inner shells. Its spectral features are then simulated by using a Monte-Carlo radiative transfer technique for different sets of velocities. Some position-velocity diagrams from simulations are very similar to the observed ones. According to the model that best reproduces the observational data, the two bright lobes and the nitrogen-enriched knots are moving away from HD 148937 at about 120 km s$^{-1}$. Their minimum kinematic age is estimated to be about 7,500 years. We discuss possible formation mechanisms of this nebula in the context of binary interaction.

We present a new hydrodynamic scheme named Godunov Density-Independent Smoothed Particle Hydrodynamics (GDISPH), that can accurately handle shock waves and contact discontinuities without any manually tuned parameters. This is in contrast to the standard formulation of smoothed particle hydrodynamics (SSPH), which requires the parameters for an artificial viscosity term to handle the shocks and struggles to accurately handle the contact discontinuities due to unphysical repulsive forces, resulting in surface tension that disrupts pressure equilibrium and suppresses fluid instabilities. While Godunov SPH (GSPH) can handle the shocks without the parameters by using solutions from a Riemann solver, it still cannot fully handle the contact discontinuities. Density-Independent Smoothed Particle Hydrodynamics (DISPH), one of several schemes proposed to handle contact discontinuities more effectively than SSPH, demonstrates superior performance in our tests involving strong shocks and contact discontinuities. However, DISPH still requires the artificial viscosity term. We integrate the Riemann solver into DISPH in several ways, yielding some patterns of GDISPH. The results of standard tests such as the one-dimensional Riemann problem, pressure equilibrium, Sedov-Taylor, and Kelvin-Helmholtz tests are favourable to GDISPH Case 1 and GDISPH Case 2, as well as DISPH. We conclude that GDISPH Case 1 has an advantage over GDISPH Case 2, effectively handling shocks and contact discontinuities without the need for specific parameters or kernels and without introducing any additional numerical diffusion.

Astrochemistry has been widely developed as a power tool to probe physical properties of the interstellar medium (ISM) in various conditions of the Milky Way (MW) Galaxy, and in near and distant galaxies. Most current studies conventionally apply linear scaling to all elemental abundances based on the gas-phase metallicity. However, these elements, including carbon and oxygen, are enriched differentially by stellar nucleosynthesis and the overall galactic chemical evolution, evident from $\alpha$-enhancement in multiple galactic observations such as starbursts, high-redshift star-forming galaxies, and low-metallicity dwarfs. We perform astrochemical modeling to simulate the impact of an $\alpha$-enhanced ISM gas cloud on the abundances of the three phases of carbon (C$^+$, C, CO) dubbed as `the carbon cycle'. The ISM environmental parameters considered include two cosmic-ray ionization rates ($\zeta_{\rm CR}=10^{-17}$ and $10^{-15}\,{\rm s}^{-1}$), two isotropic FUV radiation field strengths ($\chi/\chi_0=1$ and $10^2$), and (sub-)linear dust-to-gas relations against metallicity, mimicking the ISM conditions of different galaxy types. In galaxies with [C/O] $<$ 0, CO, C and C$^+$ all decrease in both abundances and emission, though with differential biases. The low-$J$ CO emission is found to be the most stable tracer for the molecular gas, while C and C$^+$ trace H$_2$ gas only under limited conditions, in line with recent discoveries of [CI]-dark galaxies. We call for caution when using [CII]~$158\mu$m and [CI](1-0) as alternative H$_2$-gas tracers for both diffuse and dense gas with non-zero [C/O] ratios.

Large-scale structure (LSS) surveys will increasingly provide stringent constraints on our cosmological models. Recently, the density-marked correlation function (MCF) has been introduced, offering an easily computable density-correlation statistic. Simulations have demonstrated that MCFs offer additional, independent constraints on cosmological models beyond the standard two-point correlation (2PCF). In this study, we apply MCFs for the first time to SDSS CMASS data, aiming to investigate the statistical information regarding clustering and anisotropy properties in the Universe and assess the performance of various weighting schemes in MCFs. Upon analyzing the CMASS data, we observe that, by combining different weights ($\alpha = [-0.2, 0, 0.2, 0.6]$), the MCFs provide a tight and independent constraint on the cosmological parameter $\Omega_m$, yielding $\Omega_m = 0.293 \pm0.006$ at the $1\sigma$ level, which represents a significant reduction in the statistical error by a factor of 3.4 compared to that from 2PCF. Our constraint is consistent with recent findings from the small-scale clustering of BOSS galaxies \cite{arXiv:2203.08999v2} within the 1$\sigma$ level. However, we also find that our estimate is lower than the Planck measurements by about 2.6$\sigma$, indicating the potential presence of new physics beyond the standard cosmological model if all the systematics are fully corrected. The method outlined in this study can be extended to other surveys and datasets, allowing for the constraint of other cosmological parameters. Additionally, it serves as a valuable tool for forthcoming emulator analysis on the Chinese Space Station Telescope (CSST).

We present archive Atacama Large Millimeter/Submillimeter Array (ALMA) Band 6 observations of the $^{13}$CO (J=2--1) and $^{12}$CO (J=2--1) molecular line emission of the protostellar system associated with HH 30. The $^{13}$CO molecular line shows the accretion disk while the molecular outflow is traced by the emission of the $^{12}$CO molecular line. We estimated a dynamical mass for the central object of $0.45\pm0.14\,\msun$, and a mass for the molecular outflow of $1.83\pm0.19\times10^{-4}\,\msun$. The molecular outflow presents an internal cavity as well as multiple outflowing shell structures. We distinguish three different shells with constant expansion ($\sim4-6\,\kms$) and possible rotation signatures ($\leq0.5\,\kms$). We find that the shells can be explained by magnetocentrifugal disk winds with launching radii $R_\mathrm{launch}\lesssim4\,\au$ and a small magnetic lever arm $\lambda\sim1.6-1.9$. The multiple shell structure may be the result of episodic ejections of the material from the accretion disk associated with three different epochs with dynamical ages of $497\pm15$ yr, $310\pm9$ yr, and $262\pm11$ yr for the first, second, and third shells, respectively. The outermost shell was ejected $187\pm17$ yr before the medium shell, while the medium shell was launched $48\pm14$ yr before the innermost shell. Our estimations of the linear and angular momentum rates of the outflow as well as the accretion luminosity are consistent with the expected values if the outflow of HH 30 is produced by a wide-angle disk wind.

Observations of gravitational waves from binary black hole (BBH) mergers have measured the redshift evolution of the BBH merger rate. The number density of galaxies in the Universe evolves differently with redshift based on their physical properties, such as their stellar masses and star formation rates. In this work we show that the measured population-level redshift distribution of BBHs sheds light on the properties of their probable host galaxies. We first assume that the hosts of BBHs can be described by a mixture model of galaxies weighted by stellar mass or star formation rate, and find that we can place upper limits on the fraction of mergers coming from a stellar mass weighted sample of galaxies. We then constrain parameters of a physically motivated power-law delay-time distribution using GWTC-3 data, and self-consistently track galaxies in the \textsc{UniverseMachine} simulations with this delay time model to infer the probable host galaxies of BBHs over a range of redshifts. We find that the inferred host galaxy distribution at redshift $z=0.21$ has a median star formation rate $\sim 0.9\,M_\odot\mathrm{yr}^{-1}$ and a median stellar mass of $\sim 1.9 \times 10^{10}\,M_\odot$. We also provide distributions for the mean stellar age, halo mass, halo radius, peculiar velocity, and large scale bias associated with the host galaxies, as well as their absolute magnitudes in the B- and ${ \rm K_s}$-bands. Our results can be used to design optimal electromagnetic follow-up strategies for BBHs, and also to aid the measurement of cosmological parameters using the statistical dark siren method.

We present the results of our search for nearby planetary companions of transiting hot Jupiters in 12 planetary systems: HAT-P-24, HAT-P-39, HAT-P-42, HAT-P-50, KELT-2, KELT-15, KELT-17, WASP-23, WASP-63, WASP-76, WASP-79, and WASP-161. Our analysis was based on multi-sector time-series photometry from the Transiting Exoplanet Survey Satellite and precise transit timing data sets. We detected no additional transiting planets down to the 2-4 Earth radii regime. For ten hot Jupiters, no departure from linear transit ephemerides was observed. Whilst we refute long-term variations of the orbital period for WASP-161 b, which were claimed in the literature, we notice a tentative hint of the orbital period shortening for WASP-79 b. In addition, we spot a short-period transit timing variation for KELT-2A b with the characteristics typical of the so-called exomoon corridor. We conclude, however, that further observations are required to confirm these findings.

During the first four all-sky surveys eRASS:4 carried out from December 2019 to 2021, the extended Roentgen Survey with an Imaging Telescope Array (eROSITA) on board Spektrum-Roentgen-Gamma (Spektr-RG, SRG) observed the Galactic HII region Carina nebula. We analysed the eRASS:4 data to study the distribution and the spectral properties of the hot interstellar plasma and the bright stellar sources in the Carina nebula. Spectral extraction regions of the diffuse emission were defined based on X-ray spectral morphology and multi-wavelength data. The spectra were fit with a combination of thermal and non-thermal emission models. X-ray bright point sources in the Carina nebula are the colliding wind binary $\eta$ Car, several O stars, and Wolf-Rayet (WR) stars. We extracted the spectrum of the brightest stellar sources, which can be well fit with a multi-component thermal plasma model. The spectra of the diffuse emission in the brighter parts of the Carina nebula is well reproduced by two thermal models, a lower-temperature component ($\sim$0.2 keV) and a higher-temperature component (0.6 - 0.8 keV). An additional non-thermal component dominates the emission above $\sim$1 keV in the central region around $\eta$ Car and the other massive stars. Significant orbital variation of the X-ray flux was measured for $\eta$ Car, WR22 and WR25. $\eta$ Car requires an additional time-variable thermal component in the spectral model, which is associated to the wind-wind-collision zone. Properties like temperature, pressure, and luminosity of the X-ray emitting plasma in the Carina nebula derived from the eROSITA data are consistent with theoretical calculations of emission from superbubbles. It confirms that the X-ray emission is caused by the hot plasma inside the Carina nebula which has been shocked-heated by the stellar winds of the massive stars, in particular, of $\eta$ Car.

Alignments of disc galaxies were thought to result from tidal torquing, where tidal field of the cosmic large-scale structure exert torquing moments onto dark matter haloes, determining their angular momentum and ultimately the orientation of galactic discs. In this model, resulting intrinsic ellipticity correlations are typically present on small scales, but neither observations nor simulations have found empirical evidence; instead, simulations point at the possibility that alignments of disc galaxies follow a similar alignment model as elliptical galaxies, but with a weaker alignment amplitude. In our article we make the case for the theory of linear alignments resulting from tidal distortions of the galactic disc, investigate the physical properties of this model and derive the resulting angular ellipticity spectra, as they would appear as a contribution to weak gravitational lensing in surveys such as Euclid's. We discuss in detail on the statistical and physical properties of tidally induced alignments in disc galaxies, as they are relevant for mitigation of alignment contamination in weak lensing data, comment on the consistency between the alignment amplitude in spiral and elliptical galaxies and finally, estimate their observability with a Euclid-like survey.

We report the discovery of a geocoronal solar wind charge exchange (SWCX) event corresponding to the well-known 2006 December 13th coronal mass ejection (CME) event. Strong evidence for the charge exchange origin of this transient diffuse emission is provided by prominent non-thermal emission lines at energies of $\rm O^{7+}$, $\rm Ne^{9+}$, $\rm Mg^{11+}$, $\rm Si^{12+}$, $\rm Si^{13+}$. Especially, a 0.53 keV emission line that most likely arises from the $\rm N^{5+}$ $1s^1 5p^1 \to 1s^2$ transition is detected. Previously, the forecastability of SWCX occurrence with proton flares has been disputed. In this particular event, we found that the SWCX signal coincided with the arrival of the magnetic cloud inside CME, triggered with a time delay after the proton flux fluctuation as the CME shock front passed through the Earth. Moreover, a spacecraft orbital modulation in SWCX light curve suggests that the emission arises close to the Earth. The line of sight was found to always pass through the northern magnetospheric cusp. The SWCX intensity was high when the line of sight passed the dusk side of the cusp, suggesting an azimuthal anisotropy in the flow of solar-wind ions inside the cusp. An axisymmetric SWCX emission model is found to underestimate the observed peak intensity by a factor of about 50. We suggest this discrepancy is related to the azimuthal anisotropy of the solar-wind flow in the cusp.

Aiming at improving the survey efficiency of the Wide Field Survey Telescope, we have developed a basic scheduling strategy that takes into account the telescope characteristics, observing conditions, and weather conditions at the Lenghu site. The sky area is divided into rectangular regions, referred to as `tiles', with a size of 2.577 deg * 2.634 deg slightly smaller than the focal area of the mosaic CCDs. These tiles are continuously filled in annulars parallel to the equator. The brightness of the sky background, which varies with the moon phase and distance from the moon, plays a significant role in determining the accessible survey fields. Approximately 50 connected tiles are grouped into one block for observation. To optimize the survey schedule, we perform simulations by taking into account the length of exposures, data readout, telescope slewing, and all relevant observing conditions. We utilize the Greedy Algorithm for scheduling optimization. Additionally, we propose a dedicated dithering pattern to cover the gaps between CCDs and the four corners of the mosaic CCD array, which are located outside of the 3 deg field of view. This dithering pattern helps to achieve relatively uniform exposure maps for the final survey outputs.

This thesis summarises my scientific works in the field of thermo-chemical modelling of planet-forming discs since 2009, in particular the development of the Protoplanetary Disc Model (ProDiMo). By combining chemical rate networks with continuum & line radiative transfer, and the calculation of all relevant dust and gas heating & cooling rates in an axisymmetric disc structure, these models make detailed predictions about the molecular composition of the disc, the discs' internal gas and dust temperature structure, and the composition of ice layers which form on the refractory dust grain surfaces. The development of ProDiMo was a process that took about 15 years. I have started the ProDiMo project and am the main developer, but without the involvement of an international team of scientists, in particular Inga Kamp, Wing-Fai Thi and Christian Rab, ProDiMo would not have its capabilities and would not be at the level of international recognition that is has achieved to date. Using formal solutions of the line & continuum radiative transfer, we can predict the spectral appearance of these discs from optical to millimetre wavelengths, for example the continuum and line fluxes, monochromatic images, radial intensity profiles, high-resolution line profiles that probe the disc dynamics, visibilities and channel maps. A large part of this thesis describes the publications that compared these predictions to disc observations that have been obtained by various space-borne and ground-based astronomical instruments, in particular Herschel/PACS, Spitzer/IRS, VLT/CRIRES, JWST/MIRI and ALMA. These observations probe the gas and the dust in different radial disc regions and in different layers above the midplane. This thesis summarises the conclusions drawn from the ProDiMo thermo-chemical disc models about the chemical and physical state of protoplanetary discs as the birth places of exoplanets.

The dusty winds of cool evolved stars are a major contributor of the newly synthesised material enriching the Galaxy and future generations of stars. However, the details of the physics and chemistry behind dust formation and wind launching have yet to be pinpointed. Recent spatially resolved observations show the importance of gaining a more comprehensive view of the circumstellar chemistry, but a comparative study of the intricate interplay between chemistry and physics is still difficult because observational details such as frequencies and angular resolutions are rarely comparable. Aiming to overcome these deficiencies, ATOMIUM is an ALMA Large Programme to study the physics and chemistry of the circumstellar envelopes of a diverse set of oxygen-rich evolved stars under homogeneous observing conditions at three angular resolutions between ~0.02"-1.4". Here we summarize the molecular inventory of these sources, and the correlations between stellar parameters and molecular content. Seventeen oxygen-rich or S-type asymptotic giant branch (AGB) and red supergiant (RSG) stars have been observed in several tunings with ALMA Band 6, targeting a range of molecules to probe the circumstellar envelope and especially the chemistry of dust formation close to the star. We systematically assigned the molecular carriers of the spectral lines and measured their spectroscopic parameters and the angular extent of the emission of each line from integrated intensity maps. Across the ATOMIUM sample, we detect 291 transitions of 24 different molecules and their isotopologues. This includes several first detections in oxygen-rich AGB/RSG stars: PO v=1, SO2 v1=1 and v2=2, and several high energy H2O transitions. We also find several first detections in S-type AGB stars: vibrationally excited HCN v2=2,3 and SiS v=4,5,6, as well as first detections of the molecules SiC, AlCl, and AlF in W Aql...

We present a scheme based on artificial neural networks (ANN) to estimate the line-of-sight velocities of individual galaxies from an observed redshift-space galaxy distribution. By training the network with environmental characteristics surrounding each galaxy in redshift space, our ANN model can accurately predict the line-of-sight velocity of each individual galaxy. When this velocity is used to eliminate the RSD effect, the two-point correlation function (TPCF) in real space can be recovered with an accuracy better than 1% at $s$ > 8 $h^{-1}\mathrm{Mpc}$, and 4% on all scales compared to ground truth. The real-space power spectrum can be recovered within 3% on $k$< 0.5 $\mathrm{Mpc}^{-1}h$, and less than 5% for all $k$ modes. The quadrupole moment of the TPCF or power spectrum is almost zero down to $s$ = 10 $h^{-1}\mathrm{Mpc}$ or all $k$ modes, indicating an effective correction of the spatial anisotropy caused by the RSD effect. We demonstrate that on large scales, without additional training with new data, our network is adaptable to different cosmological models, and mock galaxy samples at high redshifts and high biases, achieving less than 10% error for scales greater than 15 $h^{-1}\mathrm{Mpc}$. As it is sensitive to large-scale densities, it does not manage to remove Fingers of God in large clusters, but works remarkably well at recovering real-space galaxy positions elsewhere. Our scheme provides a novel way to predict the peculiar velocity of individual galaxies, to eliminate the RSD effect directly in future large galaxy surveys, and to reconstruct the 3-D cosmic velocity field accurately.

We show that the spectral energy distribution (SED) of the tightly focused radiation generated by the superluminally moving current sheet in the magnetosphere of a non-aligned neutron star fits the gamma-ray spectra of the Crab, Vela and Geminga pulsars over the entire range of photon energies so far detected by Fermi-LAT, MAGIC and H.E.S.S. from them: over $10^2$ MeV to $20$ TeV. While emblematic of any emission that entails caustics, the SED introduced here radically differs from those of the disparate emission mechanisms currently invoked in the literature to fit the data in different sections of these spectra. We specify, moreover, the connection between the values of the fit parameters for the analysed spectra and the physical characteristics of the central neutron stars of the Crab, Vela and Geminga pulsars and their magnetospheres.

We present a novel method for measuring the lags of (weak) variability components in neutron-star and black-hole low-mass X-ray binaries (LMXBs). For this we assume that the power and cross spectra of these sources consists of a number of components that are coherent in different energy bands, but are incoherent with one another. The technique is based on fitting simultaneously the power spectrum (PS) and the Real and Imaginary parts of the cross spectrum (CS) with a combination of Lorentzian functions. We show that, because the PS of LMXBs is insensitive to signals with a large Imaginary part and a small Real part in the CS, this approach allows us to uncover new variability components that are only detected in the CS. We also demonstrate that, contrary to earlier claims, the frequency of the type-C quasi-periodic oscillation (QPO) in the black-hole binary GRS 1915+105 does not depend on energy. Rather, the apparent energy dependence of the QPO frequency can be explained by the presence of a separate QPO component with a slightly higher frequency than that of the QPO, whose rms amplitude increases faster with energy than the rms amplitude of the QPO. From all the above we conclude that, as in the case of the PS, the CS of black-hole and neutron-star binaries can be fitted by a combination of Lorentzian components. Our findings provide evidence that the frequency-dependent part of the transfer function of these systems can be described by a combination of responses, each of them acting over relatively well-defined time scales. This conclusion challenges models that assume that the main contribution to the lags comes from a global, broadband, transfer function of the accreting system.

We report the discovery of 40 new satellite dwarf galaxy candidates in the sphere of influence of the Sombrero galaxy (M104) the most luminous galaxy in the Local Volume. Using the Subaru Hyper Suprime-Cam, we surveyed 14.4 square degrees of its surroundings, extending to the virial radius. Visual inspection of the deep images and GALFIT modelling yielded a galaxy sample highly complete down to $M_{g}\sim-9$ ($L_{g}\sim3\times 10^{5}\,L_\odot$) and spanning magnitudes $-16.4 < M_g < -8$ and half-light radii $50\,pc\ <\ r_e\ <\ 1600\,pc$ assuming the distance of M104. These 40 new, out of which 27 are group members with high confidence, double the number of potential satellites of M104 within the virial radius, placing it among the richest hosts in the Local Volume. Using a Principle Component Analysis (PCA), we find that the entire sample of candidates consistent with an almost circular on-sky distribution, more circular than any comparable environment found in the Illustris TNG100-1 simulation. However the distribution of the high probability sample is more oblate and consistent with the simulation. The cumulative satellite luminosity function is broadly consistent with analogues from the simulation, albeit it contains no bright satellite with $M_{g}<-16.4$ ($L_{g}\sim3 \times 10^{8}\,L_\odot$), a $2.3\,\sigma$ occurrence. Follow-up spectroscopy to confirm group membership will begin to demonstrate how these systems can act as probes of the structure and formation history of the halo of M104.

Giant Star-forming Clumps (GSFCs) are areas of intensive star-formation that are commonly observed in high-redshift (z>1) galaxies but their formation and role in galaxy evolution remain unclear. High-resolution observations of low-redshift clumpy galaxy analogues are rare and restricted to a limited set of galaxies but the increasing availability of wide-field galaxy survey data makes the detection of large clumpy galaxy samples increasingly feasible. Deep Learning, and in particular CNNs, have been successfully applied to image classification tasks in astrophysical data analysis. However, one application of DL that remains relatively unexplored is that of automatically identifying and localising specific objects or features in astrophysical imaging data. In this paper we demonstrate the feasibility of using Deep learning-based object detection models to localise GSFCs in astrophysical imaging data. We apply the Faster R-CNN object detection framework (FRCNN) to identify GSFCs in low redshift (z<0.3) galaxies. Unlike other studies, we train different FRCNN models not on simulated images with known labels but on real observational data that was collected by the Sloan Digital Sky Survey Legacy Survey and labelled by volunteers from the citizen science project `Galaxy Zoo: Clump Scout'. The FRCNN model relies on a CNN component as a `backbone' feature extractor. We show that CNNs, that have been pre-trained for image classification using astrophysical images, outperform those that have been pre-trained on terrestrial images. In particular, we compare a domain-specific CNN -`Zoobot' - with a generic classification backbone and find that Zoobot achieves higher detection performance and also requires smaller training data sets to do so. Our final model is capable of producing GSFC detections with a completeness and purity of >=0.8 while only being trained on ~5,000 galaxy images.

The Event Horizon Telescope (EHT) observations carried out in 2018 April at 1.3 mm wavelengths included 9 stations in the array, comprising 7 single-dish telescopes and 2 phased arrays. The metadata package for the 2018 EHT observing campaign contains calibration tables required for the a-priori amplitude calibration of the 2018 April visibility data. This memo is the official documentation accompanying the release of the 2018 EHT metadata package, providing an overview of the contents of the package. We describe how telescope sensitivities, gain curves and other relevant parameters for each station in the EHT array were collected, processed, and validated to produce the calibration tables.

This paper presents the first HI results extracted from the SARAO MeerKAT Galactic Plane Survey (SMGPS) $-$ a narrow strip ($b \sim 3^\circ$) along the southern Milky Way. The primary goal consisted in tracing the Great Attractor (GA) Wall across the innermost Zone of Avoidance. We reduced a segment spanning the longitude range $302^\circ \leq \ell \leq 332^\circ$ for the redshift range $z \leq 0.08$. The superb SMGPS sensitivity (rms = 0.3-0.5 mJy beam$^{-1}$ per 44 kms$^{-1}$ channel) and angular resolution ($\sim$ 31" $\times$ 26") lead to a detection limit of log$(M_{\rm HI}/$M$_\odot) \geq$ 8.5 at the GA distance ($V_{\rm hel} \sim 3500 - 6500$ kms$^{-1}$). A total of 477 galaxy candidates were identified over the full redshift range. A comparison of the few HI detections with counterparts in the literature (mostly HIZOA) found the HI fluxes and other HI parameters to be highly consistent. The continuation of the GA Wall is confirmed through a prominent overdensity of $N = 214$ detections in the GA distance range. At higher latitudes, the wall moves to higher redshifts, supportive of a possible link with the Ophiuchus cluster located behind the Galactic Bulge. This deep interferometric HI survey demonstrates the power of the SMGPS in improving our insight of large-scale structures at these extremely low latitudes, despite the high obscuration and continuum background.

We perform a systematic, multiwavelength spectral energy distribution (SED) analysis of X-ray detected Active Galactic Nuclei (AGNs) at $z=0.2-0.8$ with SDSS counterparts in the Stripe 82 region, consisting of 60 type-1 and 137 type-2 AGNs covering a 2--10 keV luminosity range of $41.6 < {\rm log}\ L_{\rm x} < 44.7$. The latest CIGALE code, where dusty polar components are included, is employed. To obtain reliable host and AGN parameters in type-1 AGNs, we utilize the image decomposed optical SED of host galaxies by Li et al. (2021) based on the Subaru Hyper Suprime-Cam (HSC) images. The mean ratio of black hole masses ($M_{\rm BH}$) and stellar masses ($M_{\rm stellar}$) of our X-ray detected type-1 AGN sample, $\log (M_{\rm BH}/M_{\rm stellar}) = -2.7\pm0.5$, is close to the local relation between black hole and stellar masses, as reported by Li et al. (2021) for SDSS quasars. This ratio is slightly lower than that found for more luminous ($\log L_{\rm bol} > 45$) type-1 AGNs at $z\sim1.5$. This can be explained by the AGN-luminosity dependence of $\log (M_{\rm BH}/M_{\rm stellar})$, which little evolves with redshift. We confirm the trend that the UV-to-X-ray slope ($\alpha_{\rm OX}$) or X-ray-to-bolometric correction factor ($\kappa_{2-10}$) increases with AGN luminosity or Eddington ratio. We find that type-1 and type-2 AGNs with the same luminosity ranges share similar host stellar-mass distributions, whereas type-2s tend to show smaller AGN luminosities than type-1s. This supports the luminosity (or Eddington ratio) dependent unified scheme.

Aims. Recurring jets are observed in the solar atmosphere, which can erupt intermittently. By the observation of intermittent jets, we want to understand the causes of periodic eruption characteristics. Methods. We report intermittent jets observed by the GST. The analysis was aided by 1400 {\AA} and 2796 {\AA} data from IRIS.These observational instruments allowed us to analyze the temporal characteristics of jet events. Results. The jet continued for up to 4 hours. The time distance diagram shows that the peak of the jet has obviously periodic eruption characteristics (5 minutes) during 18:00 UT-18:50 UT. We also found periodic brightening phenomenon (5 minutes) during jets bursts in the observed bands in the Transition Region (1400 {\AA} and 2796 {\AA}), which may be a response to intermittent jets in the upper solar atmosphere.The time lag is 3 minutes. Evolutionary images in the TiO band revealed the horizontal movement of granulation at the location of jet. Compared to the quiet region of the Sun, we found that the footpoint of the jet is enhanced at the center of the H {\alpha} spectral line profile, with no significant changes in the line wings. This suggests the presence of prolonged heating at the footpoint of the jet. In the mixed-polarity magnetic field region of the jet, we observed magnetic flux emergence, cancellation, and shear indicating possible intermittent magnetic reconnection. That is confirmed by the NLFFF model reconstructed using the magneto-friction method. Conclusions. The multi-wavelength analysis indicates that the events we studied were triggered by magnetic reconnection caused by mixed-polarity magnetic fields. We suggest that the horizontal motion of the granulation in the photosphere drives the magnetic reconnection, which is modulated by p-mode oscillations.

We use K-band multi-object near-infrared spectroscopy with Keck/MOSFIRE to search for environmental imprints on the gas properties of 27 narrow-band selected H$\alpha$ emitters (HAEs) across the three major clumps of the assembling USS1558--003 protocluster at $z=2.53$. We target the H$\alpha$ and [NII]$\lambda$6584 emission lines to obtain star-formation rates (SFR) and gas-phase oxygen abundances for our sources, confirming the membership of 23 objects. HAEs belonging to this protocluster display enhanced SFRs with respect to the main sequence of star formation at the same cosmic epoch. This effect is more prominent for low-mass galaxies ($\mathrm{\log M_*/M_\odot<10.0}$), which may be experiencing a vigorous phase of mass assembly shortly after they were formed. We compute the individual and stacked gas-phase metallicities for our sources finding a metallicity deficit for low-mass objects when compared against the field mass-metallicity relation and the massive Spiderweb protocluster at $z=2.16$. These results suggest that HAEs within USS1558--003 may be less evolved than those in the Spiderweb protocluster. Finally, we explore the gas metallicity - gas fraction relation for a small sample of five galaxies with CO(3-2) molecular gas information. Assuming our objects are in equilibrium, we obtain a relatively wide range of mass loading factors ($\mathrm{\lambda=0.5-2}$) matching field samples at the cosmic noon but in contrast with our previous results in the Spiderweb protocluster. We speculate that these discrepancies between protoclusters may be (partly) driven by differences in their current dynamical and mass assembly stages, hinting at the co-evolution of protoclusters and their galaxy populations at $2<z<3$.

We present and discuss the results of the Giant Metrewave Radio Telescope (GMRT) HI 21-cm line mapping for five isolated low-mass (M_bary ~(2--8)*10^7 Mo) eXtremely Metal Poor (XMP) dwarfs [12+log(O/H)=7.13-7.28], selected from the Nearby Void Galaxy (NVG) sample. All the studied void dwarfs show the disturbed morphology in the HI maps with the angular resolutions of ~11" to ~40". We examine the HI morphology and velocity field and the relative orientation of their stellar and gas body spins. We discuss the overall non-equilibrium state of their gas and the possible origin and evolution of the studied void dwarfs. The most straightforward interpretation of the ubiquitous phenomenon of the gas component non-equilibrium state in these and similar void dwarfs is the cold accretion from the void filaments and/or minor mergers. The cold gas accretion in voids could be linked to the presence of small filaments that constitute the substructure of voids.

Using data from ALFALFA, xGASS, HI-MaNGA and SDSS, we identify a sample of 49 "red but HI-rich"(RR) galaxies with $NUV-r > 5$ and unusually high HI-to-stellar mass ratios. We compare the optical properties and local environments between the RR galaxies and a control sample of "red and HI-normal"(RN) galaxies that are matched in stellar mass and color. The two samples are similar in the optical properties typical of massive red (quenched) galaxies in the local Universe. The RR sample tends to be associated with lower density environments and has lower clustering amplitudes and smaller neighbor counts at scales from several $\times$100kpc to a few Mpc. The results are consistent with that the RR galaxies preferentially locate at the center of low-mass halos, with a median halo mass $\sim 10^{12}h^{-1}M_{\odot}$ compared to $\sim 10^{12.5}h^{-1}M_{\odot}$ for the RN sample. This result is confirmed by the SDSS group catalog which reveals a central fraction of 90% for the RR sample, compared to $\sim 60\%$ for the RN sample. If assumed to follow the HI size-mass relation of normal galaxies, the RR galaxies have an average HI-to-optical radius ratio of $R_{HI}/R_{90}\sim 4$, four times the average ratio for the RN sample. We compare our RR sample with similar samples in previous studies, and quantify the population of RR galaxies using the SDSS complete sample. We conclude that the RR galaxies form a unique but rare population, accounting for only a small fraction of the massive quiescent galaxy population. We discuss the formation scenarios of the RR galaxies.

Using integral field spectroscopy from MaNGA, we study the resolved microstructures in a shocked region in Criss Cross Nebula (CCN), with an unprecedentedly high resolution of $\lesssim$1000 AU. We measure surface brightness maps for 34 emission lines, which can be broadly divided into three categories: (1) the [OIII] $\lambda$5007-like group including seven high-ionization lines and two [OII] auroral lines which uniformly present a remarkable lane structure, (2) the H$\alpha$ $\lambda$6563-like group including 23 low-ionization or recombination lines which present a clump-like structure, and (3) [OII] $\lambda$3726 and [OII] $\lambda$3729 showing high densities at both the [OIII] $\lambda$5007 lane and the H$\alpha$ clump. We use these measurements to constrain resolved shock models implemented in MAPPINGS V. We find our data can be reasonably well-fitted by a model which includes a plane-parallel shock with a velocity of $133\pm5$ km/s, plus an isotropic two-dimensional Gaussian component which is likely another clump of gas ionized by photons from the shocked region, and a constant background. We compare the electron density and temperature profiles as predicted by our model with those calculated using observed emission line ratios. We find different line ratios to provide inconsistent temperature maps, and the discrepancies can be attributed to observational effects caused by limited spatial resolution and projection of the shock geometry, as well as contamination of the additional Gaussian component. Implications on shock properties and perspectives on future IFS-based studies of CCN are discussed.

We study the size-mass relation (SMR) and recent star formation history (SFH) of post-starburst (PSB) galaxies in the local Universe, using spatially resolved spectroscopy from the final data release of MaNGA. Our sample includes 489 PSB galaxies: 94 cPSB galaxies with central PSB regions, 85 rPSB galaxies with ring-like PSB regions and 310 iPSB galaxies with irregular PSB regions. When compared to control galaxies of similar SFR, redshift and mass, a similar SMR is found for all types of PSB samples except the cPSB galaxies which have smaller sizes at intermediate masses ($9.5\lesssim \log_{10}(\rm M_\ast/M_\odot)\lesssim 10.5$). The iPSB galaxies in the star-forming sequence (iPSB-SF) show no/weak gradients in $\textrm{D}_{n}(4000)$, $\textrm{EW}(\textrm{H}\delta_{A})$ and $\textrm{EW}(\textrm{H}\alpha)$, consistent with the global star-forming status of this type of galaxies, while the quiescent iPSB (iPSB-Q) sample shows negative gradients in $\textrm{D}_{n}(4000)$ and positive gradients in $\textrm{EW}(\textrm{H}\delta_{A})$, indicating older stellar populations in the inner regions. Both cPSB and rPSB samples show positive gradients in $\textrm{D}_{n}(4000)$ and negative gradients in $\textrm{EW}(\textrm{H}\delta_{A})$, indicating younger stellar populations in the inner regions. These results imply that the four types of PSB galaxies can be broadly divided into two distinct categories in terms of evolutionary pathway: (1) iPSB-SF and iPSB-Q which have SMRs and SFHs similar to control galaxies, preferring an inside-out quenching process, (2) rPSB and cPSB which appear to be different stages of the same event, likely to follow the outside-in quenching process driven by disruption events such as mergers that result in a more compact structure as quenching proceeds.

Quasi-periodic oscillation (QPOs) analysis is important for understanding the dynamical behavior of many astrophysical objects during transient events such as gamma-ray bursts, solar flares, magnetar flares, and fast radio bursts. In this paper, we analyze QPO data in low-mass X-ray binary (LMXB) systems, using the Lense-Thirring, Kerr, and approximate Zipoy-Voorhees metrics. We demonstrate that the inclusion of spin and quadrupole parameters modifies the well-established results for the fundamental frequencies in the Schwarzschild spacetime. We interpret the QPO data within the framework of the standard relativistic precession model, allowing us to infer the values of the mass, spin, and quadrupole parameters of neutron stars in LMXBs. We explore recent QPO data sets from eight distinct LMXBs, assess their optimal parameters, and compare our findings with results in the existing literature. Finally, we discuss the astrophysical implications of our findings.

A complete accounting of nearby objects -- from the highest-mass white dwarf progenitors down to low-mass brown dwarfs -- is now possible, thanks to an almost complete set of trigonometric parallax determinations from Gaia, ground-based surveys, and Spitzer follow-up. We create a census of objects within a Sun-centered sphere of 20-pc radius and check published literature to decompose each binary or higher-order system into its separate components. The result is a volume-limited census of $\sim$3,600 individual star formation products useful in measuring the initial mass function across the stellar ($<8 M_\odot$) and substellar ($\gtrsim 5 M_{Jup}$) regimes. Comparing our resulting initial mass function to previous measurements shows good agreement above 0.8$M_\odot$ and a divergence at lower masses. Our 20-pc space densities are best fit with a quadripartite power law, $\xi(M) = dN/dM \propto M^{-\alpha}$ with long-established values of $\alpha = 2.3$ at high masses ($0.55 < M < 8.00 M_\odot$) and $\alpha = 1.3$ at intermediate masses ($0.22 < M < 0.55 M_\odot$), but at lower masses we find $\alpha = 0.25$ for $0.05 < M <0.22 M_\odot$ and $\alpha = 0.6$ for $0.01 < M < 0.05 M_\odot$. This implies that the rate of production as a function of decreasing mass diminishes in the low-mass star/high-mass brown dwarf regime before increasing again in the low-mass brown dwarf regime. Correcting for completeness, we find a star to brown dwarf number ratio of, currently, 4:1, and an average mass per object of 0.41 $M_\odot$.

Galactic fountains driven by star formation result in a variety of kinematic structures such as ionised winds and thick gas disks, both of which manifest as complex emission line profiles that can be parametrised by multiple Gaussian components. We use integral field spectroscopy (IFS) from the SAMI Galaxy Survey to spectrally resolve these features, traced by broad H$\alpha$ components, and distinguish them from the star-forming thin disk, traced by narrow components, in 3068 galaxies in the local Universe. Using a matched sample analysis technique, we demonstrate that the presence of complex emission line profiles in star-forming galaxies is most strongly correlated with the global star formation rate (SFR) surface density of the host galaxy measured within $1R_{\rm e}$ ($\Sigma_{{\rm SFR},R_{\rm e}}$), even when controlling for both observational biases, including inclination, amplitude-to-noise and angular scale, and sample biases in parameters such as stellar mass and SFR. Leveraging the spatially resolved nature of the dataset, we determine that the presence of complex emission line profiles within individual spaxels is driven not only by the local $\Sigma_{\rm SFR}$, but by the $\Sigma_{{\rm SFR},R_{\rm e}}$ of the host galaxy. We also parametrise the clumpiness of the SFR within individual galaxies, and find that $\Sigma_{{\rm SFR},R_{\rm e}}$ is a stronger predictor of the presence of complex emission line profiles than clumpiness. We conclude that, with a careful treatment of observational effects, it is possible to identify structures traced by complex emission line profiles, including winds and thick ionised gas disks, at the spatial and spectral resolution of SAMI using the Gaussian decomposition technique.

We propose to employ emission line luminosities obtained via optical spectroscopy to estimate the content of neutral hydrogen (HI) in galaxies. We use the optical spectroscopy data from the Mapping of Nearby Galaxies at APO (MaNGA) survey released in the frames of public DR17 of the Sloan Digital Sky Survey (SDSS). We compare the HI mass estimated by us for a large sample of SDSS/MaNGA galaxies with direct HI measurements from the ALFALFA survey and find a tight correlation between the masses with the correlation coefficient (CC) of 0.91 and the r.m.s scatter of 0.15 dex for the logarithmic mass. The obtained relationship is verified via another sample of MaNGA galaxies with HI masses measured with the Green Bank Telescope. Despite the coarser angular resolution of the radio data, the relation between the estimated and measured directly HI mass is tight as well - in this case CC=0.74 and the r.nm.s. is 0.29 dex. The established relations allow us to estimate the total mass of neutral hydrogen as well as the spatial distribution of HI surface density in galaxies from optical spectroscopic observations only in a simple and efficient way.

We present a novel machine-learning approach for detecting faint point sources in high-contrast adaptive optics imaging datasets. The most widely used algorithms for primary subtraction aim to decouple bright stellar speckle noise from planetary signatures by subtracting an approximation of the temporally evolving stellar noise from each frame in an imaging sequence. Our approach aims to improve the stellar noise approximation and increase the planet detection sensitivity by leveraging deep learning in a novel direct imaging post-processing algorithm. We show that a convolutional autoencoder neural network, trained on an extensive reference library of real imaging sequences, accurately reconstructs the stellar speckle noise at the location of a potential planet signal. This tool is used in a post-processing algorithm we call Direct Exoplanet Detection with Convolutional Image Reconstruction, or ConStruct. The reliability and sensitivity of ConStruct are assessed using real Keck/NIRC2 angular differential imaging datasets. Of the 30 unique point sources we examine, ConStruct yields a higher S/N than traditional PCA-based processing for 67$\%$ of the cases and improves the relative contrast by up to a factor of 2.6. This work demonstrates the value and potential of deep learning to take advantage of a diverse reference library of point spread function realizations to improve direct imaging post-processing. ConStruct and its future improvements may be particularly useful as tools for post-processing high-contrast images from the James Webb Space Telescope and extreme adaptive optics instruments, both for the current generation and those being designed for the upcoming 30 meter-class telescopes.

Observations with radio arrays that target the 21-cm signal originating from the early Universe suffer from a variety of systematic effects. An important class of these are reflections and spurious couplings between antennas. We apply a Hamiltonian Monte Carlo sampler to the modelling and mitigation of these systematics in simulated Hydrogen Epoch of Reionisation Array (HERA) data. This method allows us to form statistical uncertainty estimates for both our models and the recovered visibilities, which is an important ingredient in establishing robust upper limits on the Epoch of Reionisation (EoR) power spectrum. In cases where the noise is large compared to the EoR signal, this approach can constrain the systematics well enough to mitigate them down to the noise level for both systematics studied. Where the noise is smaller than the EoR, our modelling can mitigate the majority of the reflections with there being only a minor level of residual systematics, while cross-coupling sees essentially complete mitigation. Our approach performs similarly to existing filtering/fitting techniques used in the HERA pipeline, but with the added benefit of rigorously propagating uncertainties. In all cases it does not significantly attenuate the underlying signal.

We analyze the impact of resonant conversions mediated by non-vanishing magnetic moments between active neutrinos and a heavy sterile neutrino on the supernova neutrino flux. We present the level-crossing scheme for such a scenario and derive the neutrino fluxes after conversion, paying special attention to the order in which the resonances occur. We then compute the expected event rates from the neutronization burst of a future supernova at DUNE and Hyper-Kamiokande to derive new constraints on the neutrino magnetic moment. With this, we find a sensitivity down to a few $10^{-15} \mu_B$ for a sterile neutrino in the $O(\rm{eV})$ mass range.

The excess radio background detected by ARCADE 2 represents a puzzle within the standard cosmological model. There is no clear viable astrophysical solution and, therefore, it might indicate the presence of new physics. Radiative decays of a relic neutrino $\nu_i$ into a sterile neutrino $\nu_{\rm s}$, assumed to be quasi-degenerate, provide a solution that currently evades all constraints posed by different cosmological observations and reproduces very well the ARCADE 2 data. We find a very good fit to the ARCADE 2 data with best fit values $\tau_i = 1.46 \times 10^{21}\,{\rm s}$ and $\Delta m_i = 4.0 \times 10^{-5}\,{\rm eV}$, where $\tau_i$ is the lifetime and $\Delta m_i$ is the mass difference between the decaying active neutrino and the sterile neutrino. On the other hand, if relic neutrino decays do not explain ARCADE 2 data, then these place a stringent constraint $\D m_i^{3/2} \tau_i \gtrsim 2 \times 10^{14}\,{\rm eV}^{3/2}\,{\rm s}$ in the range $1.4 \times 10^{-5} \, {\rm eV} < \D m_i < 2.5 \times 10^{-4}\,{\rm eV}$. The solution also predicts a stronger 21 cm absorption global signal than the predicted one from the $\L$CDM model, with a contrast brightness temperature $T_{21} = -238^{+21}_{-20}\,{\rm mK}$ ($99\%$ C.L.) at redshift $z\simeq 17$. This is in mild tension with the even stronger signal found by the EDGES collaboration, $T_{21} = - 500^{+200}_{-500}\,{\rm mK} $, suggesting that this might have been overestimated, maybe receiving a contribution from some unidentified foreground source.

We compute the quasinormal modes of static and spherically symmetric black holes (BHs) with electric and magnetic charges. For the electrically charged case, the dynamics of perturbations separates into the odd- and even-parity sectors with two coupled differential equations in each sector. In the presence of both electric and magnetic charges, the differential equations of four dynamical degrees of freedom are coupled with each other between odd- and even-parity perturbations. Despite this notable modification, we show that, for a given total charge and mass, a BH with mixed electric and magnetic charges gives rise to the same quasinormal frequencies for fundamental modes. This includes the case in which two BHs have equal electric and magnetic charges for each of them. Thus, the gravitational-wave observations of quasinormal modes during the ringdown phase alone do not distinguish between electrically and magnetically charged BHs.

The recent compelling observation of the nanohertz stochastic gravitational wave background has brought to light a new galactic arena to test gravity. In this paper, we derive a formula for the most general expression of the stochastic gravitational wave background correlation that could be tested with pulsar timing and future square kilometer arrays. Our expressions extends the harmonic space analysis, also often referred to as the power spectrum approach, to predict the correlation signatures of an anisotropic polarized stochastic gravitational wave background with subluminal tensor, vector, and scalar gravitational degrees of freedom. We present the first few nontrivial anisotropy and polarization signatures in the correlation and discuss their dependence on the gravitational wave speed and pulsar distances. Our results set up tests that could potentially be used to rigorously examine the isotropy of the stochastic gravitational wave background and strengthen the existing constraints on possible non-Einsteinian polarizations in the nanohertz gravitational wave regime.

The physics of primordial black holes can be affected by the non-Gaussian statistics of the density fluctuations that generate them. Therefore, it is important to have good theoretical control of the higher-order correlation functions for primordial curvature perturbations. By working at leading order in a $1/|\eta|$ expansion, we analytically determine the bispectrum of curvature fluctuations for single field inflationary scenarios producing primordial black holes. The bispectrum has a rich scale and shape dependence, and its features depend on the dynamics of the would-be decaying mode. We apply our analytical results to study gravitational waves induced at second order by enhanced curvature fluctuations. Their statistical properties are derived in terms of convolution integrals over wide momentum ranges, and they are sensitive on the scale and shape dependence of the curvature bispectrum we analytically computed.

We present a framework for constraining the automatic sequential generation of equations to obey the rules of dimensional analysis by construction. Combining this approach with reinforcement learning, we built $\Phi$-SO, a Physical Symbolic Optimization method for recovering analytical functions from physical data leveraging units constraints. Our symbolic regression algorithm achieves state-of-the-art results in contexts in which variables and constants have known physical units, outperforming all other methods on SRBench's Feynman benchmark in the presence of noise (exceeding 0.1%) and showing resilience even in the presence of significant (10%) levels of noise.