Papers for Friday, Apr 21 2023

Dark sector theories naturally lead to multi-component scenarios for dark matter where a sub-component can dissipate energy through self-interactions, allowing it to efficiently cool inside galaxies. We present the first cosmological hydrodynamical simulations of Milky Way analogues where the majority of dark matter is collisionless Cold Dark Matter (CDM), but a sub-component (6%) is strongly dissipative minimal Atomic Dark Matter (ADM). The simulations, implemented in GIZMO and utilizing FIRE-2 galaxy formation physics to model the standard baryonic sector, demonstrate that the addition of even a small fraction of dissipative dark matter can significantly impact galactic evolution despite being consistent with current cosmological constraints. We show that ADM gas with roughly Standard-Model-like masses and couplings can cool to form a rotating "dark disk" with angular momentum closely aligned with the visible stellar disk. The morphology of the disk depends sensitively on the parameters of the ADM model, which affect the cooling rates in the dark sector. The majority of the ADM gas gravitationally collapses into dark "clumps" (regions of black hole or mirror star formation), which form a prominent bulge and a rotating thick disk in the central galaxy. These clumps form early and quickly sink to the inner ~kpc of the galaxy, affecting the galaxy's star-formation history and present-day baryonic and CDM distributions.

The intracluster medium (ICM) in the centers of galaxy clusters is heavily influenced by the ``feedback'' from supermassive black holes (SMBHs). Feedback can drive turbulence in the ICM and turbulent dissipation can potentially be an important source of heating. Due to the limited spatial and spectral resolutions of X-ray telescopes, direct observations of turbulence in the hot ICM have been challenging. Recently, we developed a new method to measure turbulence in the ICM using multiphase filaments as tracers. These filaments are ubiquitous in cluster centers and can be observed at very high resolution using optical and radio telescopes. We study the kinematics of the filaments by measuring their velocity structure functions (VSFs) over a wide range of scales in the centers of $\sim 10$ galaxy clusters. We find features of the VSFs that correlate with the SMBHs activities, suggesting that SMBHs are the main driver of gas motions in the centers of galaxy clusters. In all systems, the VSF is steeper than the classical Kolmogorov expectation and the slopes vary from system to system. One theoretical explanation is that the VSFs we have measured so far mostly reflect the motion of the driver (jets and bubbles) rather than the cascade of turbulence. We show that in Abell 1795, the VSF of the outer filaments far from the SMBH flattens on small scales to a Kolmogorov slope, suggesting that the cascade is only detectable farther out with the current telescope resolution. The level of turbulent heating computed at small scales is typically an order of magnitude lower than that estimated at the driving scale. Even though SMBH feedback heavily influences the kinematics of the ICM in cluster centers, the level of turbulence it drives is rather low, and turbulent heating can only offset $\lesssim10\%$ of the cooling loss, consistent with the findings of numerical simulations.

We build a deep learning framework that connects the local formation process of dark matter halos to the halo bias. We train a convolutional neural network (CNN) to predict the final mass and concentration of dark matter halos from the initial conditions. The CNN is then used as a surrogate model to derive the response of the halos' mass and concentration to long-wavelength perturbations in the initial conditions, and consequently the halo bias parameters following the "response bias" definition. The CNN correctly predicts how the local properties of dark matter halos respond to changes in the large-scale environment, despite no explicit knowledge of halo bias being provided during training. We show that the CNN recovers the known trends for the linear and second-order density bias parameters $b_1$ and $b_2$, as well as for the local primordial non-Gaussianity linear bias parameter $b_\phi$. The expected secondary assembly bias dependence on halo concentration is also recovered by the CNN: at fixed mass, halo concentration has only a mild impact on $b_1$, but a strong impact on $b_\phi$. Our framework opens a new window for discovering which physical aspects of the halo's Lagrangian patch determine assembly bias, which in turn can inform physical models of halo formation and bias.

Motivated by observations of multiphase galaxy outflows, we explore the impact of isotropic and anisotropic electron thermal conduction on the evolution of radiatively-cooled, cold clouds embedded in hot, magnetized winds. Using the adaptive mesh refinement code AthenaPK, we conduct simulations of clouds impacted by supersonic and transonic flows with magnetic fields initially aligned parallel and perpendicular to the flow direction. In cases with isotropic thermal conduction, an evaporative wind forms, stabilizing against instabilities and leading to a mass loss rate that matches the hydrodynamic case. In anisotropic cases, the impact of conduction is more limited and strongly dependent on the field orientation. In runs with initially transverse fields, the field lines are folded back into the tail, strongly limiting conduction, but magnetic fields act to dampen instabilities and slow the stretching of the cloud in the flow direction. In the aligned case, anisotropic conduction aids cloud survival by forming a radiative wind near the front of the cloud, which suppresses instabilities with minimal mass loss. In all cases, anisotropic conduction has a minimal impact on the acceleration of the cloud.

About ten percent of Sun-like ($1$-$2 M_\odot$) stars will engulf a $1$-$10 M_{\rm J}$ planet as they expand during the red giant branch (RGB) or asymptotic giant branch (AGB) phase of their evolution. Once engulfed, these planets experience a strong drag force in the star's convective envelope and spiral inward, depositing energy and angular momentum. For these mass ratios, the inspiral takes $\sim 10$-$10^{2}$ years ($\sim 10^{2}$-$10^{3}$ orbits); the planet undergoes tidal disruption at a radius of $\sim R_\odot$. We use the Modules for Experiments in Stellar Astrophysics (MESA) software instrument to track the stellar response to the energy deposition while simultaneously evolving the planetary orbit. For RGB stars, as well as AGB stars with $M_{\rm p} \lesssim 5 M_{\rm J}$ planets, the star responds quasistatically but still brightens measurably on a timescale of years. In addition, asteroseismic indicators, such as the frequency spacing or rotational splitting, differ before and after engulfment. For AGB stars, engulfment of a $M_{\rm p} \gtrsim 5 M_{\rm J}$ planet drives supersonic expansion of the envelope, causing a bright, red, dusty eruption similar to a "luminous red nova." Based on the peak luminosity, color, duration, and expected rate of these events, we suggest that engulfment events on the AGB could be a significant fraction of low-luminosity red novae in the Galaxy. We do not find conditions where the envelope is ejected prior to the planet's tidal disruption, complicating the interpretation of short-period giant planets orbiting white dwarfs as survivors of common-envelope evolution.

Next-generation ground-based gravitational-wave detectors, such as Cosmic Explorer (CE), are expected to be sensitive to gravitational-wave signals with frequencies as low as 5 Hz, allowing signals to spend a significant amount of time in the detector frequency band. As a result, the effects caused by the rotation of the Earth become increasingly important for such signals. Additionally, the length of the arms of these detectors can be comparable to the wavelength of detectable gravitational waves, which introduces frequency-dependent effects that are not significant in current-generation detectors. These effects are expected to improve the ability to localize compact binary coalescences in the sky even when using only one detector. This study aims to understand how much these effects can help in localization. We present the first comprehensive Bayesian parameter estimation framework that accounts for all these effects using \textsc{Bilby}, a commonly used Bayesian parameter estimation tool. We focus on sky localization constraints for binary neutron star events with an optimal signal-to-noise ratio of 1000 with one detector at the projected CE sensitivity. We find that these effects help localize sources using one detector with sky areas as low as 10 square degrees. Moreover, we explore and discuss how ignoring these effects in the parameter estimation can lead to biases in the inference.

For the majority of short period exoplanets transiting massive stars with radiative envelopes, the spin angular momentum of the host star is greater than the planetary orbital angular momentum. In this case, the orbits of the planets will undergo nodal precession, which can significantly impact the probability that the planets transit their parent star. In particular, for some combinations of the spin-orbit angle $\psi$ and the inclination of the stellar spin $i_*$, all such planets will eventually transit at some point over the duration of their precession period. Thus, as the time over which the sky has been monitored for transiting planets increases, the frequency of planets with detectable transits will increase, potentially leading to biased estimates of exoplanet occurrence rates, especially orbiting more massive stars. Furthermore, due to the dependence of the precession period on orbital parameters such as spin-orbit misalignment, the observed distributions of such parameters may also be biased. We derive the transit probability of a given exoplanet in the presence of nodal precession induced by a rapidly spinning host star. We find that the effect of nodal precession has already started to become relevant for some short-period planets, i.e., Hot Jupiters, orbiting massive stars, by increasing transit probabilities by of order a few percent for such systems within the original $Kepler$ field. We additionally derive simple expressions to describe the time evolution of the impact parameter $b$ for applicable systems, which should aid in future investigations of exoplanet nodal precession and spin-orbit alignment.

Elucidating the nature of Dark Matter (DM), which does not interact with light and which interacts with known matter primarily or only through gravity, is one of the principal quests in physics. Leading candidates for DM are weakly interacting massive particles (WIMPs) or ultralight bosons (axions), at opposite extremes in mass scales, that have been postulated by competing theories to solve deficiencies in the Standard Model of particle physics. Whereas DM WIMPs behave like discrete particles ($\varrho$DM), quantum interference between DM axions is manifested as waves ($\psi$DM). Here, we show that gravitational lensing leaves signatures in multiply-lensed images of background galaxies that reveal whether the foreground lensing galaxy inhabits a $\varrho$DM or $\psi$DM halo. Whereas $\varrho$DM lens models leave well documented anomalies between the predicted and observed brightnesses and positions of multiply-lensed images, $\psi$DM lens models correctly predict the level of anomalies left over by $\varrho$DM lens models. More challengingly, when subjected to a battery of tests for reproducing the quadruply-lensed triplet images in the system HS 0810+2554, $\psi$DM is able to reproduce all aspects of this system whereas $\varrho$DM often fails. The growing success of $\psi$DM in reproducing astrophysical observations tilt the balance toward new physics invoking axions.

During the first half of their main-sequence lifetimes, stars rapidly lose angular momentum to their magnetized winds, a process known as magnetic braking. Recent observations suggest a substantial decrease in the magnetic braking efficiency when stars reach a critical value of the Rossby number, the stellar rotation period normalized by the convective overturn timescale. Cooler stars have deeper convection zones with longer overturn times, reaching this critical Rossby number at slower rotation rates. The nature and timing of the transition to weakened magnetic braking has previously been constrained by several solar analogs and two slightly hotter stars. In this Letter, we derive the first direct constraints from stars cooler than the Sun. We present new spectropolarimetry of the old G8 dwarf $\tau$ Cet from the Large Binocular Telescope, and we reanalyze a published Zeeman Doppler image of the younger G8 star 61 UMa, yielding the large-scale magnetic field strengths and morphologies. We estimate mass-loss rates using archival X-ray observations and inferences from Ly$\alpha$ measurements, and we adopt other stellar properties from asteroseismology and spectral energy distribution fitting. The resulting calculations of the wind braking torque demonstrate that the rate of angular momentum loss drops by a factor of 300 between the ages of these two stars (1.4-9 Gyr), well above theoretical expectations. We summarize the available data to help constrain the value of the critical Rossby number, and we identify a new signature of the long-period detection edge in recent measurements from the Kepler mission.

Rapidly outflowing cold H-I gas is ubiquitously observed to be co-spatial with a hot phase in galactic winds, yet the ablation time of cold gas by the hot phase should be much shorter than the acceleration time. Previous work showed efficient radiative cooling enables clouds to survive in hot galactic winds, as can magnetic fields even in purely adiabatic simulations for sufficiently small density contrasts between the wind and cloud. In this work, we study the interplay between radiative cooling and magnetic draping via three dimensional radiative magnetohydrodynamic simulations. We find magnetic fields decrease the critical cloud radius for survival by two orders of magnitude (i.e., to sub-pc scales) in the strongly magnetized ($\beta_{\rm wind}=1$) case. Our results show magnetic fields (i) accelerate cloud entrainment through magnetic draping, (ii) can cause faster cloud destruction in cases of inefficient radiative cooling, (iii) do not significantly suppress mass growth for efficiently cooling clouds, and, crucially, in combination with radiative cooling (iv) reduce the average overdensity by providing non-thermal pressure support of the cold gas. This substantially reduces the acceleration time compared to the destruction time (more than due to draping alone), enhancing cloud survival. Our results may help to explain the cold, tiny, rapidly outflowing cold gas observed in galactic winds and the subsequent high covering fraction of cold material in galactic halos.

The heavy element content ("metallicity") of the Universe is a record of the total star formation history. Gas-phase metallicity in galaxies, as well as its evolution with time, is of particular interest as a tracer of accretion and outflow processes. However, metallicities from the widely-used electron temperature ($T_e$) method are typically ~2x lower than the values based on the recombination line method. This "abundance discrepancy factor" (ADF) is well known and is commonly ascribed to bias due to temperature fluctuations. We present a measurement of oxygen abundance in the nearby (3.4 Mpc) system, Mrk 71, using a combination of optical and far-IR emission lines to measure and correct for temperature fluctuation effects. Our far-IR result is inconsistent ($> 2 \sigma$ significance) with the metallicity from recombination lines and instead indicates little to no bias in the standard $T_e$ method, ruling out the long-standing hypothesis that the ADF is explained by temperature fluctuations for this object. Our results provide a framework to accurately measure metallicity across cosmic history, including with recent data reaching within the first billion years with JWST and the Atacama Large Millimeter Array (ALMA).

We present analysis of nearly 1,000 near-infrared, integrated light spectra from APOGEE in the inner $\sim$7 kpc of M31. We utilize full spectrum fitting with A-LIST simple stellar population spectral templates that represent a population of stars with the same age, [M/H], and [$\alpha$/M]. With this, we determine the mean kinematics, metallicities, $\alpha$ abundances, and ages of the stellar populations of M31's bar, bulge, and inner disk ($\sim$4-7 kpc). We find a non-axisymmetric velocity field in M31 resulting from the presence of a bar. The bulge of M31 is metal-poor relative to the disk ([M/H] = $-0.149^{+0.067}_{-0.081}$ dex), features minima in metallicity on either side of the bar ([M/H] $\sim$ -0.2), and is enhanced in $\alpha$ abundance ([$\alpha$/M] = $0.281^{+0.035}_{-0.038}$). The disk of M31 within $\sim$7 kpc is enhanced in both metallicity ([M/H] = $-0.023^{+0.050}_{-0.052}$) and $\alpha$ abundance ([$\alpha$/M] = $0.274^{+0.020}_{-0.025}$). Both of these structural components are uniformly old at $\simeq$ 12 Gyr. We find the metallicity increases with distance from the center of M31, with the steepest gradient along the disk major axis ($0.043\pm0.021$ dex/kpc). This gradient is the result of changing light contributions from the metal-poor bulge and metal-rich disk. The chemodynamics of stellar populations encodes information about a galaxy's chemical enrichment, star formation history, and merger history, allowing us to discuss new constraints on M31's formation. Our results provide a stepping stone between our understanding of the Milky Way and other external galaxies.

Context. Lightcurve variability is well-suited for characterising objects in surveys with high cadence and long baseline. This is especially relevant in view of the large datasets to be produced by the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST). Aims. We aim to determine variability parameters for objects in the Palomar Transient Factory (PTF) and explore differences between quasars (QSOs), stars and galaxies. We will relate variability and colour information in preparation for future surveys. Methods. We fit joint likelihoods to structure functions (SFs) of 71 million PTF lightcurves with a Markov Chain Monte Carlo method. For each object, we assume a power law SF and extract two parameters: the amplitude on timescales of one year, $A$, and a power law index, $\gamma$. With these parameters and colours in the optical (Pan-STARRS1) and mid infrared (WISE), we identify regions of parameter space dominated by different types of spectroscopically confirmed objects from SDSS. Candidate QSOs, stars and galaxies are selected to show their parameter distributions. Results. QSOs have high amplitude variations in the $R$ band, and the strongest timescale dependence of variability. Galaxies have a broader range of amplitudes and low timescale dependency. With variability and colours, we achieve a photometric selection purity of 99.3 % for QSOs. Even though hard cuts in monochromatic variability alone are not as effective as seven-band magnitude cuts, variability is useful in characterising object sub-classes. Through variability, we also find QSOs that were erroneously classified as stars in the SDSS. We discuss perspectives and computational solutions in view of the upcoming LSST.

Active dwarf galaxies are important because they contribute to the evolution of dwarf galaxies and can reveal their hosted massive black holes. However, the sample size of such sources beyond the local universe is still highly limited. In this work, we search for active dwarf galaxies in the recently completed XMM-Spitzer Extragalactic Representative Volume Survey (XMM-SERVS). XMM-SERVS is currently the largest medium-depth X-ray survey covering 13 $\mathrm{deg}^2$ in three extragalactic fields, which all have well-characterized multi-wavelength information. After considering several factors that may lead to misidentifications, we identify 73 active dwarf galaxies at $z<1$, which constitutes the currently largest X-ray-selected sample beyond the local universe. Our sources are generally less obscured than predictions based on the massive-AGN (active galactic nucleus) X-ray luminosity function and have a low radio-excess fraction. We find that our sources reside in similar environments to inactive dwarf galaxies. We further quantify the accretion distribution of the dwarf-galaxy population after considering various selection effects and find that it decreases with X-ray luminosity, but redshift evolution cannot be statistically confirmed. Depending upon how we define an AGN, the active fraction may or may not show a strong dependence on stellar mass. Their Eddington ratios and X-ray bolometric corrections significantly deviate from the expected relation, which is likely caused by several large underlying systematic biases when estimating the relevant parameters for dwarf galaxies. Throughout this work, we also highlight problems in reliably measuring photometric redshifts and overcoming strong selection effects for distant active dwarf galaxies.

Context. Large, high-dimensional astronomical surveys require efficient data analysis. Automatic fitting of lightcurve variability and machine learning may assist in identification of sources including candidate quasars. Aims. We aim to classify sources from the Palomar Transient Factory (PTF) as quasars, stars or galaxies, and to examine model performance using variability and colours. We determine the added value of variability information as well as quantifying the performance when colours are not available. Methods. We use supervised learning in the form of a histogram-based gradient boosting classifier to predict spectroscopic SDSS classes using photometry. For comparison, we create models with structure function variability parameters only, magnitudes only and using all parameters. Results. We achieve highly accurate predictions for 71 million sources with lightcurves in PTF. The full model correctly identifies 92.49 % of spectroscopically confirmed quasars from the SDSS with a purity of 95.64 %. With only variability, the completeness is 34.97 % and the purity is 58.71 % for quasars. The predictions and probabilities of PTF objects belonging to each class are made available in a catalogue, VILLAIN-Cat, including magnitudes and variability parameters. Conclusions. We have developed a method for automatic and effective classification of PTF sources using magnitudes and variability. For similar supervised models, we recommend using at least 100,000 labeled objects, and we show how performance scales with data volume.

Lateral distribution functions of particles in extensive air showers with the energy $E_0 \simeq 10^{19}$ eV recorded by ground-based and underground scintillation detectors with a threshold of $E_{\mu} \simeq 1.0 \times \sec\theta$ GeV at the Yakutsk array during the continuous observations from 1986 to 2016 have been analyzed using events with zenith angles $\theta \le 60^{\circ}$ functions have been compared to the predictions obtained with the QGSJet01 hadron interaction model by applying the CORSIKA code. The entire dataset indicates that cosmic rays consist predominantly of protons.

We report the discovery of millihertz quasi-periodic oscillations (mHz QPOs) from the neutron star (NS) low-mass X-ray binary 4U 1730--22 using the Neutron Star Interior Composition Explorer (NICER). After being inactive for almost 50 years, 4U 1730--22 went into outburst twice between June and August 2021, and between February and July 2022. We analyse all the NICER observations of this source, and detect mHz QPOs with a significance > $4\sigma$ in 35 observations. The QPO frequency of the full data set ranged between ~4.5 and ~8.1 mHz with an average fractional rms amplitude of the order of ~2%. The X-ray colour analysis strongly suggests that 4U 1730--22 was in a soft spectral state during the QPO detections. Our findings are consistent with those reported for other sources where the mHz QPOs have been interpreted as the result of a special mode of He burning on the NS surface called marginally stable nuclear burning (MSNB). We conclude that the mHz QPOs reported in this work are also associated with the MSNB, making 4U 1730--22 the eighth source that shows this phenomenology. We discuss our findings in the context of the heat flux from the NS crust.

The dipper is a novel class of young stellar object associated with large drops in flux on the order of 10 to 50 per cent lasting for hours to days. Too significant to arise from intrinsic stellar variability, these flux drops are currently attributed to disk warps, accretion streams, and/or transiting circumstellar dust. Dippers have been previously studied in young star forming regions including the Orion Complex. Using Next Generation Transit Survey (NGTS) data, we identified variable stars from their lightcurves. We then applied a machine learning random forest classifier for the identification of new dipper stars in Orion using previous variable classifications as a training set. We discover 120 new dippers, of which 83 are known members of the Complex. We also investigated the occurrence rate of disks in our targets, again using a machine learning approach. We find that all dippers have disks, and most of these are full disks. We use dipper periodicity and model-derived stellar masses to identify the orbital distance to the inner disk edge for dipper objects, confirming that dipper stars exhibit strongly extended sublimation radii, adding weight to arguments that the inner disk edge is further out than predicted by simple models. Finally, we determine a dipper fraction (the fraction of stars with disks which are dippers) for known members of 27.8 plus minus 2.9 per cent. Our findings represent the largest population of dippers identified in a single cluster to date.

Accurate detection of the cosmological 21-cm global signal requires galactic foreground models that can fit spectra down to $\sim 20$ mK or less, representing a removal of power over nearly six orders of magnitude. Rarely are such models tested to this level, let alone their dependence upon model inputs like sky temperature maps. We therefore test the ability of seven commonly employed foreground models -- including nonlinear and linear forward-models, polynomials, and maximally-smooth polynomials -- to fit realistic simulated mock spectra, as well as their dependence upon model inputs. The mock spectra are synthesized from intrinsic foregrounds with realistic spatial and spectral structure, chromatic beams, horizon profiles, and discrete time-sampling. For a single LST bin spectrum, the nonlinear-forward model with 4 parameters is preferred using a KS-test of the noise-normalized residuals, while the linear forward-model fits well with 6-7 parameters. The polynomials and maximally-smooth polynomials, like those employed by the EDGES and SARAS3 experiments, cannot produce good fits with 5 parameters. However, we find that polynomials with 6 parameters pass the KS-test, although a 9 parameter fit produces the highest p-value. When fitting multiple LST bins simultaneously to decrease overlap with global signal models, we find that the linear forward-model outperforms the nonlinear for 2, 5 and 10 LST bins. In addition, the nonlinear forward-model fails to produce good fits to spectra with 10 LST bins, in contrast to the linear. Importantly, the KS-test consistently identifies best-fit \textit{and} preferred models as opposed to the $\chi^2_{red}$ and Bayesian evidence, especially in cases involving nonlinear models.

Using spatially-resolved fossil record analysis on a large sample of 'red and dead' elliptical galaxies (classical ellipticals, CLEs) from the MaNGA/SDSS-IV DR15 survey, we reconstruct the archaeological evolution of their radial gradients in mass-to-luminosity ratio ($M/L$), $g-r$ color, and specific star formation (SF) rate. We also calculate other metrics that quantify the inside-out SF quenching and external mass growth processes. The $M/L$ gradients, $\nabla\Upsilon_{\star}$, are approximately flat at high look-back times ($t_{\rm lb}$), but then they become negative and steeper until an epoch, when this trend reverses. These trends are shifted to later epochs the less massive the galaxies are. Color gradients follow qualitatively similar trends. We find that these trends are mainly driven by strong inside-out quenching, without significant outer growth or structural changes overall. Our results suggest a scenario where the main progenitors of local CLE galaxies evolved quasi-passively after an early dissipative phase, but underwent radial photometric changes due to the inside-out quenching that led to the systematic decrease of $\nabla\Upsilon_{\star}$ and to an increase of the light-weighted radius. The late reversing of $\nabla\Upsilon_{\star}$, $t_{\rm lb}\approx2-4$ Gyr, roughly coincides with the global quenching of the CLE galaxies. We have pushed archaeological inferences to the limit, but thanks to the large number of objects and an understanding of how the caveats and assumptions affect our results, we conclude that they offer an average description of evolutionary behaviors of CLE progenitors that is valid at least qualitatively.

Triple AGN systems are expected to be the result of the hierarchical model of galaxy formation. Since there are very few of them confirmed as such, we present the results of a new study of the triple-AGN candidate SDSS J102700.40+174900.8 (center nucleus) through observations with $\it{GTC}$-$\it{MEGARA}$ Integral Field Unit. 1D and 2D analysis of the line ratios of the three nuclei allow us to locate them in the EW(H$\alpha$) vs. [Nii] /H$\alpha$ diagram. The central nucleus is found to be a retired galaxy (or fake AGN). The neighbors are found to be a strong AGN (southeastern nucleus, J102700.55+174900.2) compatible with a Sy2 galaxy, and a weak AGN (northern nucleus, J102700.38+174902.6) compatible with a LINER2. We find evidence that the neighbors constitute a dual AGN system (Sy2-LINER2) with a projected separation of 3.98 kpc in the optical bands. The H$\alpha$ velocity map shows that the northern nucleus has an H$\alpha$ emission with a velocity offset of $\sim$-500 km s$^{-1}$, whereas the southeastern nucleus has a rotating disk and H$\alpha$ extended emission at kpc scales. Chandra archival data confirm that the neighbors have X-ray (0.5-2) keV and (2-7) keV emission, whereas the center nucleus shows no X-ray emission. A collisional ring with knots is observed in the HST images of the southeastern nucleus. These knots coincide with star formation regions that along with the ring are predicted in a head-on collision. In this case, the morphology changes are probably due to a minor merger that was produced by the passing of the northern through the southeastern nucleus.

The classification of galaxy morphology is a hot issue in astronomical research. Although significant progress has been made in the last decade in classifying galaxy morphology using deep learning technology, there are still some deficiencies in spatial feature representation and classification accuracy. In this study, we present a multi-scale convolutional capsule network (MSCCN) model for the classification of galaxy morphology. First, this model improves the convolutional layers through using a multi-branch structure to extract multi-scale hidden features of galaxy images. In order to further explore the hidden information in the features, the multi-scale features are encapsulated and fed into the capsule layer. Second, we use a sigmoid function to replace the softmax function in dynamic routing, which can enhance the robustness of MSCCN. Finally, the classification model achieving 97% accuracy, 96% precision, 98% recall, and 97% F1-score under macroscopic averaging. In addition, a more comprehensive model evaluation were accomplished in this study. We visualized the morphological features for the part of sample set, which using the t-distributed stochastic neighbor embedding (t-SNE) algorithm. The results shows that the model has the better generalization ability and robustness, it can be effectively used in the galaxy morphological classification.

X-ray binaries are known to show state transitions related to accretion rate changes which are often accompanied with dramatic changes in the jet emission. However, it is not clear whether this characteristics of stellar-mass black hole systems can be scaled up to the accretion disk of active galactic nuclei. The Seyfert 1 galaxy, KUG 1141+371 has been showing a steadily increasing X-ray flux since 2007, and exhibited variability behaviour similar to the state transitions observed in X-ray binaries. It was hypothesised to undergo a rapid boost of mass accretion. If the X-ray binary analogy holds then the appearance of jet emission can also be expected in KUG 1141+371. While the source was not detected in the Faint Images of the Radio Sky at Twenty-centimeters in 1994, it appears in the VLA Sky Survey in 2019 and at 22 GHz in a VLA observation in 2018 at mJy flux density level. Our VLBI observations revealed a compact, flat-spectrum radio feature. Its high brightness temperature indicates the radio emission originates from an AGN.

While collisionless cold dark matter models have been largely successful in explaining a wide range of observational data, some tensions still exist, and it remains possible that dark matter possesses a non-negligible level of self interactions. In this paper, we investigate a possible observable consequence of self-interacting dark matter: offsets between the central galaxy and the center of mass of its parent halo. We examine 23 relaxed galaxy clusters in a redshift range of 0.1 to 0.3 drawn from clusters in the Dark Energy Survey and the Sloan Digital Sky Survey which have archival Chandra X-ray data of sufficient depth for center and relaxation determination. We find that most clusters in our sample show non-zero offsets between the X-ray center, taken to be the centroid within the cluster core, and the central galaxy position. All of the measured offsets are larger, typically by an order of magnitude, than the uncertainty in the X-ray position due to Poisson noise. In all but six clusters, the measured offsets are also larger than the estimated, combined astrometric uncertainties in the X-ray and optical positions. A more conservative cut on concentration to select relaxed clusters marginally reduces but does not eliminate the observed offset. With our more conservative sample, we find an estimated mean X-ray to central galaxy offset of $\mu = 5.5 \pm 1.0$ kpc. Comparing to recent simulations, this distribution of offsets is consistent with some level of dark matter self interaction, though further simulation work is needed to place constraints.

An open question in SN Ia research is where the boundary lies between 'normal' Type Ia supernovae (SNe Ia) that are used in cosmological measurements and those that sit off the Phillips relation. We present the spectroscopic modelling of one such '86G-like' transitional SN Ia, SN 2021rhu, that has recently been employed as a local Hubble Constant calibrator using a tip of the red-giant branch measurement. We detail its modelling from -12 d until maximum brightness using the radiative-transfer spectral-synthesis code tardis. We base our modelling on literature delayed-detonation and deflagration models of Chandrasekhar mass white dwarfs, as well as the double-detonation models of sub-Chandrasekhar mass white dwarfs. We present a new method for 'projecting' abundance profiles to different density profiles for ease of computation. Due to the small velocity extent and low outer densities of the W7 profile, we find it inadequate to reproduce the evolution of SN 2021rhu as it fails to match the high-velocity calcium components. The host extinction of SN 2021rhu is uncertain but we use modelling with and without an extinction correction to set lower and upper limits on the abundances of individual species. Comparing these limits to literature models we conclude that the spectral evolution of SN 2021rhu is also incompatible with double-detonation scenarios, lying more in line with those resulting from the delayed detonation mechanism (although there are some discrepancies, in particular a larger titanium abundance in SN 2021rhu compared to the literature). This suggests that SN 2021rhu is likely a lower luminosity, and hence lower temperature, version of a normal SN Ia.

We present the dynamical evolution of ten open clusters which were part of our previous studies. These clusters include both young and intermediate-age open clusters with ages ranging from 25$\pm$19 Myr to 1.78$\pm$0.20 Gyr. The total mass of these clusters ranges from 356.18$\pm$142.90 to 1811.75$\pm$901.03 M$_{\odot}$. The Galactocentric distances to the clusters are in the range of 8.91$\pm$0.02 to 11.74$\pm$0.18 kpc. The study is based on the ground-based UBVRI data supplemented by the astrometric data from the Gaia archive. We studied the minimum spanning tree of the member stars for these clusters. The mass segregation in these clusters was quantified by mass segregation ratios calculated from the mean edge length obtained through the minimum spanning tree. The clusters NGC 2360, NGC 1960, IC 1442, King 21, and SAI 35 have ${\Gamma}_{MSR}$ to be 1.65$\pm$0.18, 1.94$\pm$0.22, 2.21$\pm$0.20, 1.84$\pm$0.23, and 1.96$\pm$0.25, respectively which indicate moderate mass segregation in these clusters. The remaining five clusters are found to exhibit weak or no mass segregation. We used the ratio of half mass radius to the tidal radius i.e. R$_{h}$/R$_{t}$ to investigate the effect of the tidal interactions on the cluster structure and dynamics. The ratios of half mass radii to tidal radii are found to be positively correlated with the Galactocentric distances with a linear slope of 0.06$\pm$0.01 having linear regression coefficient r-square = 0.93 for the clusters.

We study the behaviour of dust in galaxies at cosmic dawn, z=6-8, by coupling the FirstLight simulations with the radiative transfer code POLARIS. The starburst nature of these galaxies and their complex distribution of dust lead to a large diversity of attenuation curves. These follow the Calzetti model only for relatively massive galaxies, Mstars=10^9 Msun. Galaxies with lower masses have steeper curves, consistent with the model for the Small Magellanic Cloud (SMC). The ultraviolet and optical slopes of the attenuation curves are closer to the modified Calzetti model, with a slight preference for the power-law model for galaxies with the highest values of attenuation. We have also examined the relation between the slope in the far-ultraviolet, beta_UV , and the infrared excess, IRX. At z=6, it follows the Calzetti model with a shift to slightly lower beta_UV values due to lower metallicities at lower attenuation. The same relation at z=8 shows a shift to higher IRX values due to a stronger CMB radiation at high-z.

We use parameterized post-Friedmann (PPF) description for dark energy and apply ellipsoidal nested sampling to perform the Bayesian model selection method on different time-dependent dark energy models using a combination of $Planck$ and data based on distance measurements, namely baryon acoustic oscillations and supernovae luminosity distance. Models with two and three free parameters described in terms of linear scale factor $a$, or scaled in units of e-folding $\ln a$ are considered. Our results show that parameterizing dark energy in terms of $\ln a$ provides better constraints on the free parameters than polynomial expressions. In general, two free-parameter models are adequate to describe the dynamics of the dark energy compared to their three free-parameter generalizations. According to the Bayesian evidence, determining the strength of support for cosmological constant $\Lambda$ over polynomial dark energy models remains inconclusive. Furthermore, considering the $R$ statistic as the tension metric shows that one of the polynomial models gives rise to a tension between $Planck$ and distance measurements data sets. The preference for the logarithmic equation of state over $\Lambda$ is inconclusive, and the strength of support for $\rm \Lambda$CDM over the oscillating model is moderate.

Many asymptotic giant branch (AGB) and supergiant stars exhibit extended detached shells in the far-infrared, resembling rings or arcs. These structures have long been interpreted as the bow shock formed in the interface between the stellar wind and the interstellar medium, the astrosphere. To date, only a few AGB stars have been observed showing an extended shell in the ultraviolet: the cometary tail drifting away from $o$ Ceti, and a bubble around IRC+10216, CIT6, and U Hya. This paper describes a search of UV extended shells around AGB stars using archival GALEX far-UV images. After inspecting visually 282 GALEX images, we identified the fourth discovery of a UV bubble around the AGB star R Dor. The bubble is seen as a 26'x29' ring, corresponding to an actual diameter of 0.41x0.46 parsec$^2$. The mass of the thin UV bubble is estimated to be $\simeq$0.003 $M_{\odot}$. The morphological asymmetry (less than $\sim 20$\%) and brightness variations of this shell are uncorrelated with the stellar proper motion and thus they can rather be ascribed to inhomogeneities in the ISM. Archival \emph{IRAS} 60 and 100$\mu$m images reveal that the bubble is filled with cold (i.e. < 32 K) dust. All UV bubbles known to date are limited to be within a distance < 350 pc and at high Galactic latitudes (|b| > 35 degree), which suggests that their detection is hampered in most cases by the strong UV interstellar extinction.

The observed hard TeV gamma-ray spectrum of the nearby gamma-ray burst (GRB) 190829A may challenge the conventional leptonic GRB afterglow model. It has been proposed that an ultra-high-energy (UHE; $\varepsilon^{'}_{\rm p}\sim 10^{20}$ eV) proton population can be pre-accelerated by internal shocks in GRB jets. We study possible signatures of the UHE protons embedded in the TeV afterglows when they escape the afterglow fireball. We show that the leptonic model can represent the observed multiwavelength lightcurves and spectral energy distributions of GRB 190829A by considering the uncertainties of the model parameters. Attributing the TeV gamma-ray afterglows to the emission of both the electron self-Compton scattering process and the UHE proton synchrotron radiations in the afterglow fireball, we obtain tentative upper limits of $\log_{10} \varepsilon_{\rm p}^{\prime}/{\rm eV}\sim 20.46$ and $\log_{10}E_{\rm p, total}/{\rm erg}\leq 50.75$, where $E_{\rm p, total}$ is the total energy of the proton population. The synchrotron radiations of the UHE protons should dominate the early TeV gamma-ray afterglows, implying that early observations are critical for revealing the UHE proton population.

We study the formation and evolution of solitons supported by repulsive self-interactions inside extended halos, for scalar-field dark matter scenarios. We focus on the semiclassical regime where the quantum pressure is typically much smaller than the self-interactions. We present numerical simulations, with initial conditions where the halo is described by the WKB approximation for its eigenfunction coefficients. We find that when the size of the system is of the order of the Jeans length associated with the self-interactions, a central soliton quickly forms and makes about 50% of the total mass. However, if the halo is ten times greater than this self-interaction scale, a soliton only quickly forms in cuspy halos where the central density is large enough to trigger the self-interactions. If the halo has a flat core, it takes a longer time for a soliton to appear, after small random fluctuations on the de Broglie wavelength size build up to reach a large enough density. In some cases, we observe the co-existence of several narrow density spikes inside the larger self-interaction-supported soliton. All solitons appear robust and slowly grow, unless they already make up 40% of the total mass. We develop a kinetic theory, valid for an inhomogeneous background, to estimate the soliton growth rate for low masses. It explains the fast falloff of the growth rate as resonances between the ground state and halo excited states disappear. Our results suggest that cosmological halos would show a large scatter for their soliton mass, depending on their assembly history.

ESA's Jupiter Icy Moons Explorer (JUICE) will provide a detailed investigation of the Jovian system in the 2030s, combining a suite of state-of-the-art instruments with an orbital tour tailored to maximise observing opportunities. We review the Jupiter science enabled by the JUICE mission, building on the legacy of discoveries from the Galileo, Cassini, and Juno missions, alongside ground- and space-based observatories. We focus on remote sensing of the climate, meteorology, and chemistry of the atmosphere and auroras from the cloud-forming weather layer, through the upper troposphere, into the stratosphere and ionosphere. The Jupiter orbital tour provides a wealth of opportunities for atmospheric and auroral science: global perspectives with its near-equatorial and inclined phases, sampling all phase angles from dayside to nightside, and investigating phenomena evolving on timescales from minutes to months. The remote sensing payload spans far-UV spectroscopy (50-210 nm), visible imaging (340-1080 nm), visible/near-infrared spectroscopy (0.49-5.56 $\mu$m), and sub-millimetre sounding (near 530-625\,GHz and 1067-1275\,GHz). This is coupled to radio, stellar, and solar occultation opportunities to explore the atmosphere at high vertical resolution; and radio and plasma wave measurements of electric discharges in the Jovian atmosphere and auroras. Cross-disciplinary scientific investigations enable JUICE to explore coupling processes in giant planet atmospheres, to show how the atmosphere is connected to (i) the deep circulation and composition of the hydrogen-dominated interior; and (ii) to the currents and charged particle environments of the external magnetosphere. JUICE will provide a comprehensive characterisation of the atmosphere and auroras of this archetypal giant planet.

Since a black hole does not emit light from its interior, nor does it have a surface on which light from nearby sources can be reflected, observational study of black hole physics requires observing the gravitational impact of the black hole on its surroundings. A massive black hole leaves a dynamical imprint on stars and gas close by. Gas in the immediate vicinity of an accreting massive black hole can, due to the presence of the black hole, shine so brightly that it outshines the light of the billions of stars in its host galaxy and be detected across the Universe. By observing the emission from stars and gas and determining their kinematics scientists can extract vital information not only on the fundamental properties of the black holes themselves but also the impact they have on their surroundings. As it turns out, supermassive black holes appear to play a vital role in shaping the Universe as we know it, as they can profoundly impact the star formation history in galaxies. As a consequence, these black holes indirectly impact the cosmic build up of chemical elements heavier than Helium and thus affect when and where life can form. For these reasons alone, observations of massive black holes constitute a very active research area of modern astrophysics. In this chapter we aim to provide a general overview -- fit for a non-expert -- of what scientists have learned, and hope to learn, from analyzing observations of massive black holes and the material around them.

If dark matter annihilates into standard model particles with a cross-section which is velocity dependent, then Local Group dwarf galaxies will not be the best place to search for the resulting gamma ray emission. A greater flux would be produced by more distant and massive halos, with larger velocity dispersions. We construct full-sky predictions for the gamma-ray emission from galaxy- and cluster-mass halos within $\sim 200 \, {\mathrm{Mpc}}$ using a suite of constrained $N$-body simulations (CSiBORG) based on the Bayesian Origin Reconstruction from Galaxies algorithm. Comparing to observations from the Fermi Large Area Telescope and marginalising over reconstruction uncertainties and other astrophysical contributions to the flux, we obtain constraints on the cross-section which are two (seven) orders of magnitude tighter than those obtained from dwarf spheroidals for $p$-wave ($d$-wave) annihilation. We find no evidence for either type of annihilation from dark matter particles with masses in the range $m_\chi = 2-500 \, {\mathrm{GeV}}/c^2$, for any channel. As an example, for annihilations producing bottom quarks with $m_\chi = 10 \, {\mathrm{GeV}}/c^2$, we find $a_{1} < 2.4 \times 10^{-21} \, {\mathrm{cm^3 s^{-1}}}$ and $a_{2} < 3.0 \times 10^{-18} \, {\mathrm{cm^3 s^{-1}}}$ at 95% confidence, where the product of the cross-section, $\sigma$, and relative particle velocity, $v$, is given by $\sigma v = a_\ell (v/c)^{2\ell}$ and $\ell=1, 2$ for $p$-, $d$-wave annihilation, respectively. Our bounds, although failing to exclude the thermal relic cross-section for velocity-dependent annihilation channels, are among the tightest to date.

We propose a new probe of inflationary gravitational waves (IGWs): the cross-correlation of the lensing of inflationary $B$-mode polarization with a large-scale structure (LSS) tracer, which can also be a CMB lensing map. This is equivalent to measuring a three-point function of two CMB $B$-modes and an LSS tracer. We forecast expected $1\,\sigma$ constraints on the tensor-to-scalar ratio, $r$, albeit with a simplistic foreground treatment, and find constraints of $\sigma_r \simeq 7 \times 10^{-3}$ from the correlation of CMB-S4-Deep $B$-mode lensing and LSST galaxies, $\sigma_r \simeq 5 \times 10^{-3}$ from the correlation of CMB-S4-Deep $B$-mode lensing and CMB-S4-Deep CMB lensing, and $\sigma_r \simeq 10^{-2}$ from the correlation of LiteBIRD $B$-mode lensing and CMB-S4-Wide lensing. Because this probe is inherently non-Gaussian, simple Gaussian foregrounds will not produce any biases to the measurement of $r$. While a detailed investigation of non-Gaussian foreground contamination for different cross-correlations will be essential, this observable has the potential to be a powerful probe of IGWs, complementary to standard methods for constraining $r$.

We present APO and Gemini time-series photometry of WD J004917.14$-$252556.81, an ultramassive DA white dwarf with $T_{\rm eff} = 13020$ K and $\log{g} = 9.34$. We detect variability at two significant frequencies, making J0049$-$2525 the most massive pulsating white dwarf currently known with $M_\star=1.31~M_{\odot}$ (for a CO core) or $1.26~M_{\odot}$ (for an ONe core). J0049$-$2525 does not display any of the signatures of binary mergers, there is no evidence of magnetism, large tangential velocity, or rapid rotation. Hence, it likely formed through single star evolution and is likely to have an ONe core. Evolutionary models indicate that its interior is $\gtrsim99$% crystallized. Asteroseismology offers an unprecedented opportunity to probe its interior structure. However, the relatively few pulsation modes detected limit our ability to obtain robust seismic solutions. Instead, we provide several representative solutions that could explain the observed properties of this star. Extensive follow-up time-series photometry of this unique target has the potential to discover a significant number of additional pulsation modes that would help overcome the degeneracies in the asteroseismic fits, and enable us to probe the interior of an $\approx1.3~M_{\odot}$ crystallized white dwarf.

Eccentricity is a smoking gun for the formation channel of stellar-mass binary black holes (sBBHs). Space-based gravitational wave observatories can determine binary eccentricity to $e_0\gtrsim\mathcal{O}(10^{-4}) $, but the detection of these systems can be very challenging. A targeted search of archival data triggered by ground-based detectors shrinks the search range thus making the task tractable. Previous studies ignored the effect of eccentricity. For the first time, we constructed a template bank for space-borne gravitational wave detectors that includes the impact of eccentricity. We find that even for a mild upper limit of $0.1$, the inclusion of eccentricity can still boost the template bank size by five orders of magnitudes. Our work marked a solid step towards the detection of a realistic sBBH, and it demonstrated that with the appropriate extension, the template bank method can still identify the early inspiral of sBBHs.

We used Bands 6 and 7 of the Atacama Large Millimeter/submillimeter Array (ALMA) in Cycles 7 and 8 to search for $\mathrm{CO}\,(3-2)$ emission from a sample of five fast radio burst (FRB) host galaxies discovered by the Commensal Real-time ASKAP Fast Transients (CRAFT) survey and the Fast and Fortunate for FRB Follow-up (F$^4$) team. These galaxies have redshifts $z \approx 0.16-0.48$, masses log$(M_{\rm star}/M_{\odot})\approx 9.30-10.4$ characteristic of field galaxies, and emission lines indicative of ongoing star formation. We detected three of the five galaxies with luminosities $L'(3-2)\approx0.2-4\times10^8\,\rm K\,km \, s^{-1}\,pc^2$ and set upper limits for the other two. Adopting standard metallicity-dependent CO-to-H$_2$ conversion factors, we estimate molecular gas masses $M_{\rm gas}\approx 0.2-3\times 10^9 \, M_{\odot}$. As a population, FRB host galaxies track the main $M_{\rm star}-M_{\rm gas}$ locus of star-forming galaxies in the present-day universe, with gas fractions of $\mu_{\rm gas}\approx0.1$ and gas depletion times $t_{\rm dep} \gtrapprox 1\,$Gyr. We employ the Kaplan-Meier estimator to compare the redshift-corrected $\mu_{\rm gas}$ and $t_{\rm dep}$ for all known FRB hosts with measurements or upper limits with those from the xCOLD GASS survey and find statistically different gas fractions. The difference is not statistically significant when we consider only the five hosts studied here with consistently determined properties, suggesting more FRB hosts with measured molecular gas masses are needed to robustly study the population. Lastly, we present a multi-wavelength analysis of one host (HG20180924B) combining high-spatial resolution imaging and integral field spectroscopy to demonstrate that future high-resolution observations will allow us to study the host galaxy environments local to the FRBs.

The quenched fraction of satellite galaxies is aligned with the orientation of the halo's central galaxy, such that on average, satellites form stars at a lower rate along the major axis of the central. This effect, called anisotropic satellite galaxy quenching (ASGQ), has been found in observational data and cosmological simulations. Analyzing the IllustrisTNG simulation, Mart\'in-Navarro et al. (2021) recently argued that ASGQ is caused by anisotropic energetic feedback and constitutes "compelling observational evidence for the role of black holes in regulating galaxy evolution." In this letter, we study the causes of ASGQ in state-of-the-art galaxy formation simulations to evaluate this claim. We show that cosmological simulations predict that on average, satellite galaxies along the major axis of the dark matter halo tend to have been accreted at earlier cosmic times and are hosted by subhalos of larger peak halo masses. As a result, a modulation of the quenched fraction with respect to the major axis of the central galaxy is a natural prediction of hierarchical structure formation. We show that ASGQ is predicted by the UniverseMachine galaxy formation model, a model without anisotropic feedback. Furthermore, we demonstrate that even in the IllustrisTNG simulation, anisotropic satellite accretion properties are the main cause of ASGQ. Ultimately, we argue that ASGQ is not a reliable indicator of supermassive black hole feedback in galaxy formation simulations and, thus, should not be interpreted as such in observational data.

The role of AGN in quenching galaxies and driving the evolution from star-forming to quiescent remains a key question in galaxy evolution. We present evidence from the Mapping Nearby Galaxies at APO (MaNGA) survey for fading AGN activity in 6/93 post-starburst galaxies. These six galaxies show extended emission line regions (EELRs) consistent with ionization from past AGN activity, analogous to "Hanny's voorwerp" and other systems where the OIII5007 emission is bright enough to be visible in broadband imaging. Using the infrared luminosities from IRAS to estimate the current AGN luminosity, we find that 5/6 of the post-starburst galaxies have current AGN which have faded from the peak luminosity required to have ionized the EELRs. Given the rate at which we observe EELRs, the typical EELR visibility timescale, and an estimate of how often EELRs would be visible, we estimate the duty cycle of AGN activity during the post-starburst phase. The timescale for the galaxy to cycle between peaks in AGN luminosity is $t_{\rm EELR}\sim1.1-2.3\times10^5$ yr. Given the rate at which we observe current AGN activity during this phase, we estimate that the AGN spends only 5.3% of this time (or $t_{\rm ON} = 0.6-1.3\times10^4$ yr) in its luminous phase, with the rest of the time spent "off" or in a low-luminosity phase. The length of this duty cycle may explain why so few luminous AGN have been observed during the post-starburst phase, despite evidence for AGN feedback at work.

Deep Cosmic Microwave Background polarization experiments allow in principle very precise internal reconstruction of the gravitational lensing signal. To this aim, likelihood-based or Bayesian methods are typically necessary, where performing a sometimes very large number of lensing and delensing remappings on the sphere is required before satisfactory convergence. We discuss here in some detail an optimized piece of numerical code able to perform both the lensing operation and its adjoint (closely related to delensing) efficiently, and to arbitrary accuracy, using non-uniform Fast Fourier Transform technology. Where applicable, we find the code outperforms by massive amounts current widespread software, being able to produce high-resolution maps accurate enough for next-generation CMB experiments on the timescale of seconds on a modern laptop. The adjoint operation performs similarly, and removes the need for computation of inverse deflection fields. This publicly available code enables de facto efficient spherical harmonic transforms on completely arbitrary grids, and could possibly find applications also in other areas.

The measurement of the Hubble constant $H_0$ plays an important role in the study of cosmology. In this letter, we propose a new method to constrain the Hubble constant using the strongly lensed gravitational wave (GW) signals. By reparameterizing the waveform, we find that the lensed waveform is sensitive to the $H_0$. Assuming the scenario that no electromagnetic counterpart of the GW source can be identified, our method can still give meaningful constraints on the $H_0$ with the information of the lens redshift. We then apply Fisher information matrix and Markov Chain Monte Carlo to evaluate the potential of this method. For the space-based GW detector, TianQin, the $H_0$ can be constrained within a relative error of $\sim$ 0.3-2\%, using a single strongly lensed GW event. Precision varies according to different levels of electromagnetic information.

The most luminous quasars at $z > 6$ are suspected to be both highly clustered and reside in the most massive dark matter halos in the early Universe, making them prime targets to search for galaxy overdensities and/or protoclusters. We search for Lyman-break dropout-selected galaxies using HST WFC3/ACS broadband imaging in the fields of three $6 < z < 7$ quasars, as well as their simultaneously observed coordinated-parallel fields, and constrain their photometric redshifts using EAZY. One field, J0305-3150, shows a volume density 10$\times$ higher than the blank-field UV luminosity function (UVLF) at M$_{UV} < -20$, with tentative evidence of a 3$\sigma$ overdensity in its parallel field located 15 cMpc away. Another field, J2054-0005, shows an angular overdensity within 500 ckpc from the quasar but still consistent with UVLF predictions within 3$\sigma$, while the last field, J2348-3054, shows no enhancement. We discuss methods for reducing uncertainty in overdensity measurements when using photometric selection and show that we can robustly select LBGs consistent with being physically associated with the quasar, corroborated by existing JWST/NIRCam WFSS data in the J0305 field. Even accounting for incompleteness, the overdensities in J0305 and J2054 are higher for brighter galaxies at short angular separations, suggesting preferential enhancement of more massive galaxies in the immediate vicinity of the quasar. Finally, we compare the LBG population with previously-identified [CII] and mm-continuum companions; the LBG overdensities are not accompanied by an enhanced number of dusty galaxies, suggesting that the overdense quasar fields are not in the bursty star-forming phase sometimes seen in high-redshift protoclusters.

We present 3D hydrodynamical simulations of core convection with a stably stratified envelope of a 25 $\mathrm{M}_\odot$ star in the early phase of the main-sequence. We use the explicit gas-dynamics code $\texttt{PPMstar}$ which tracks two fluids and includes radiation pressure and radiative diffusion. Multiple series of simulations with different luminosities and radiative thermal conductivities are presented. The entrainment rate at the convective boundary, internal gravity waves in and above the boundary region, and the approach to dynamical equilibrium shortly after a few convective turnovers are investigated. From the results of these simulations we extrapolate to find the entrainment rate at the nominal heating rate and thermal diffusion given by the $\texttt{MESA}$ stellar evolution model on which the 3D stratification is based. Further, to study the effect of radiative diffusion on the thermal timescale, we perform very long simulations accelerated by 10000 times their nominal luminosities. In these simulations the growing penetrative convection reduces the initially unrealistically large entrainment. This reduction is enabled by a spatial separation that develops between the entropy gradient and the composition gradient. The convective boundary moves outward much more slowly at the end of these simulations. Finally, we present a method to predict the extent and character of penetrative convection beyond the Schwarzschild boundary. This method is intended to be ultimately deployed in 1D stellar evolution calculations and is based on the properties of penetrative convection in our simulations carried forward through the local thermal timescale.

We present high spatial resolution ($\approx24$ pc) Atacama Large Millimeter/sub-millimeter Array $^{12}$CO(2-1) observations of the central region of the nearby barred spiral galaxy NGC 5806. NGC 5806 has a highly structured molecular gas distribution with a clear nucleus, a nuclear ring and offset dust lanes. We identify $170$ spatially- and spectrally-resolved giant molecular clouds (GMCs). These clouds have comparable sizes ($R_{\mathrm{c}}$) and larger gas masses, observed linewidths ($\sigma_{\mathrm{obs,los}}$) and gas mass surface densities than those of clouds in the Milky Way disc. The size -- linewidth relation of the clouds is one of the steepest reported so far ($\sigma_{\mathrm{obs,los}}\propto R_{\mathrm{c}}^{1.20}$), the clouds are on average only marginally bound (with a mean virial parameter $\langle\alpha_{\mathrm{vir}}\rangle\approx2$), and high velocity dispersions are observed in the nuclear ring. These behaviours are likely due to bar-driven gas shocks and inflows along the offset dust lanes, and we infer an inflow velocity of $\approx120$ kms$^{-1}$ and a total molecular gas mass inflow rate of $\approx5$ M$_\odot$ yr$^{-1}$ into the nuclear ring. The observed internal velocity gradients of the clouds are consistent with internal turbulence. The number of clouds in the nuclear ring decreases with azimuthal angle downstream from the dust lanes without clear variation of cloud properties. This is likely due to the estimated short lifetime of the clouds ($\approx6$ Myr), which appears to be mainly regulated by cloud-cloud collision and/or shear processes. Overall, it thus seems that the presence of the large-scale bar and gas inflows to the centre of NGC 5806 affect cloud properties.

We use the ESPRESSO spectrograph at the Very Large Telescope to measure velocity shifts and gravitational redshifts of eight bona fide Hyades white dwarfs, with an accuracy better than 1.5 percent. By comparing the gravitational redshift measurements of the mass-to-radius ratio with the same ratios derived by fitting the \textit{Gaia} photometry with theoretical models, we find an agreement to better than one per cent. It is possible to reproduce the observed white dwarf cooling sequence and the trend of the mass-to-radius ratios as a function of colour using isochrones with ages between 725 and 800 Myr, tuned for the Hyades. One star, EGGR\,29, consistently stands out in all diagrams, indicating that it is possibly the remnant of a blue straggler. We also computed mass-to-radius ratios from published gravities and masses, determined from spectroscopy. The comparison between photometric and spectroscopic stellar parameters reveals that spectroscopic effective temperature and gravity are systematically larger than the photometric values. Spectroscopic mass-to-radius ratios disagree with those measured from gravitational redshift, indicating the presence of systematics affecting the white dwarf parameters derived from the spectroscopic analysis.

The sensitivity and field of view of the Canadian Hydrogen Intensity Mapping Experiment (CHIME) has enabled its fast radio burst (FRB) backend to detect thousands of FRBs. However, the low angular resolution of CHIME prevents it from localizing most FRBs to their host galaxies. Very long baseline interferometry (VLBI) can readily provide the subarcsecond resolution needed to localize many FRBs to their hosts. Thus we developed TONE: an interferometric array of eight $6~\mathrm{m}$ dishes to serve as a pathfinder for the CHIME/FRB Outriggers project, which will use wide field of view cylinders to determine the sky positions for a large sample of FRBs, revealing their positions within their host galaxies to subarcsecond precision. In the meantime, TONE's $\sim3333~\mathrm{km}$ baseline with CHIME proves to be an excellent testbed for the development and characterization of single-pulse VLBI techniques at the time of discovery. This work describes the TONE instrument, its sensitivity, and its astrometric precision in single-pulse VLBI. We believe that our astrometric errors are dominated by uncertainties in the clock measurements which build up between successive Crab pulsar calibrations which happen every $\approx 24~\mathrm{h}$; the wider fields of view and higher sensitivity of the Outriggers will provide opportunities for higher-cadence calibration. At present, CHIME-TONE localizations of the Crab pulsar yield systematic localization errors of ${0.1}-{0.2}~\mathrm{arcsec}$ - comparable to the resolution afforded by state-of-the-art optical instruments ($\sim 0.05 ~\mathrm{arcsec}$).

Compact Steep Spectrum (CSS) radio sources are active galactic nuclei that have radio jets propagating only on galactic scales, defined as having projected linear sizes (LS) of up to $20\,$kpc. CSS sources are generally hosted by massive early-type galaxies with little on-going star formation, however a small fraction are known to have enhanced star formation. Using archival data from the Faint Images of the Radio Sky at Twenty cm survey, the Very Large Array Sky Survey and the Sloan Digital Sky Survey we identify a volume-limited sample of $166$ CSS sources at $z<0.2$ with $L_{1.4\,\text{GHz}}>10^{24}\,\text{W}\,\text{Hz}^{-1}$. Comparing the star formation rates and linear sizes of these CSS sources, we find that the $\approx14\,\%$ of CSS sources with specific star formation rates above $0.01\,\text{Gyr}^{-1}$ all have $\text{LS}<10\,$kpc. We discuss the possible mechanisms driving this result, concluding that it is likely the excess star formation in these sources occurred in multiple bursts and ceased prior to the AGN jet being triggered.

Despite achieving sensitivities capable of detecting the extremely small amplitude of gravitational waves (GWs), LIGO and Virgo detector data contain frequent bursts of non-Gaussian transient noise, commonly known as 'glitches'. Glitches come in various time-frequency morphologies, and they are particularly challenging when they mimic the form of real GWs. Given the higher expected event rate in the next observing run (O4), LIGO-Virgo GW event candidate validation will require increased levels of automation. Gravity Spy, a machine learning tool that successfully classified common types of LIGO and Virgo glitches in previous observing runs, has the potential to be restructured as a signal-vs-glitch classifier to accurately distinguish between glitches and GW signals. A signal-vs-glitch classifier used for automation must be robust and compatible with a broad array of background noise, new sources of glitches, and the likely occurrence of overlapping glitches and GWs. We present GSpyNetTree, the Gravity Spy Convolutional Neural Network Decision Tree: a multi-CNN classifier using CNNs in a decision tree sorted via total GW candidate mass tested under these realistic O4-era scenarios.

The Kerr-Newman metric is the unique vacuum solution of the General Relativistic field equations, in which any singularities or spacetime pathologies are hidden behind horizons. They are believed to describe the spacetimes of massive astrophysical objects with no surfaces, which we call black holes. This spacetime, which is defined entirely by the mass, spin, and charge of the black hole, gives rise to a variety of phenomena in the motion of particles and photons outside the horizons that have no Newtonian counterparts. Moreover, the Kerr-Newman spacetime remains remarkably resilient to many attempts in modifying the underlying theory of gravity. The monitoring of stellar orbits around supermassive black holes, the detection of gravitational waves from the coalescence of stellar-mass black holes, and the observation of black-hole shadows in images with horizon-scale resolution, all of which have become possible during the last decade, are offering valuable tools in testing quantitatively the predictions of this remarkable solution to Einstein's equations.

The following annotated bibliography contains a reasonably complete survey of contemporary work in the philosophy of astrophysics. Spanning approximately forty years from the early 1980s to the present day, the bibliography should help researchers entering the field to acquaint themselves with its major texts, while providing an opportunity for philosophers already working on astrophysics to expand their knowledge base and engage with unfamiliar material.

In this work, we assess the sensitivity reach of pulsar timing array (PTA) measurements to probe pointlike primordial black holes (PBHs), with an extended mass distribution, which originate from collapsed Fermi balls that are formed through the aggregation of asymmetric U(1) dark fermions trapped within false vacuum bubbles during a dark first order phase transition (FOPT). The PBH formation scenario is mainly characterized by the dark asymmetry, strength of the FOPT, rate of FOPT, and the percolation temperature. Meanwhile, for PBH masses of interest lying within $10^{-10} M_\odot - 10^{-3}M_\odot$, the relevant signal for PTA measurements is the Doppler phase shift in the timing signal, due to the velocity change induced by transiting PBHs on pulsars. Taking the dark asymmetry parameter to be $10^{-4}$ and $10^{-5}$, we find that percolation temperatures within the 0.1-10 keV range, FOPT rates above $10^3$ times the Hubble parameter at percolation, and FOPT strengths within $10^{-6}-0.1$ can give rise to PBHs that can be probed by an SKA-like PTA observation. On the other hand, the accompanying gravitational wave (GW) signal from the FOPT can be used as a complementary probe, assuming that the peak frequency lies within the $\mathcal{O}(10^{-9})-\mathcal{O}(10^{-6})$ Hz range, and the peak GW abundance is above the peak-integrated sensitivity curves associated with pulsar timing observations that search for stochastic GWs. At the fundamental level, a quartic effective potential for a dark scalar field can trigger the FOPT. By performing a parameter scan, we obtained the class of effective potentials that lead to FOPT scenarios that can be probed by SKA through pulsar timing and GW observations.

When we impose the discrete symmetry in the standard model we have Dai-Freed global anomalies. However, interestingly if we introduce three right-handed neutrinos we can have an anomaly-free discrete $Z_4$ gauge symmetry. This $Z_4$ symmetry should be spontaneously broken down to the $Z_2$ symmetry to generate the heavy Majorana masses for the right-handed neutrinos. We show that this symmetry breaking naturally generates topological inflation, which is consistent with the CMB observations at present and predicts a significant tensor mode with scalar-tensor ratio $r > 0.03$. The right-handed neutrinos play an important role in reheating processes. The reheating temperature is as high as $\sim 10^8$GeV, and non-thermal leptogenesis successfully takes place.

GW191219\_163120 is a gravitational wave signal that is believed to have originated from a neutron star-black hole (NSBH) coalescence with an extreme mass ratio. In this work, we use data of GW191219\_163120 from LIGO and Virgo to test the first law of black hole mechanics by considering the neutron star as a perturbation to the black hole before the merger, and the remnant black hole as a stationary black hole after the merger. Our results demonstrate consistency with the first law of black hole mechanics, with an error level of about 6\% at 68\% credibility and 10\% at 95\% credibility. We also find that the higher the mass ratio of the gravitational wave source, the more consistent our results are with the first law of black hole mechanics. Overall, our study sheds light on the nature of NSBH coalescences and their implications for black hole mechanics.

Methane is a powerful greenhouse gas, and a primary target for mitigating climate change in the short-term future due to its relatively short atmospheric lifetime and greater ability to trap heat in Earth's atmosphere compared to carbon dioxide. Top-down observations of atmospheric methane are possible via drone and aircraft surveys as well as satellites such as the TROPOspheric Monitoring Instrument (TROPOMI). Recent work has begun to apply the divergence method to produce regional methane emission rate estimates. Here we show that spatially incomplete observations of methane can produce negatively biased time-averaged regional emission rate estimates via the divergence method, but that this effect can be counteracted by adopting a procedure in which daily advective fluxes of methane are time-averaged before the divergence method is applied. Using such a procedure with TROPOMI methane observations, we calculate yearly Permian emission rates of 3.1, 2.4 and 2.7 million tonnes per year for the years 2019 through 2021. We also show that highly-resolved plumes of methane can have negatively biased estimated emission rates by the divergence method due to the presence of turbulent diffusion in the plume, but this is unlikely to affect regional methane emission budgets constructed from TROPOMI observations of methane. The results from this work are expected to provide useful guidance for future implementations of the divergence method for emission rate estimation from satellite data - be it for methane or other gaseous species in the atmosphere.