Papers for Friday, Aug 02 2024

The Lobster Eye Imager for Astronomy (LEIA) is the pathfinder of the wide-field X-ray telescope used in the Einstein Probe mission. In this study, we present an image of the Virgo Cluster taken by LEIA in the 0.5-4.5 keV band with an exposure time of $\sim$17.3 ks in the central region. This extended emission is generally consistent with the results obtained by ROSAT. However, the field is affected by bright point sources due to the instrument's Point Spread Function (PSF) effect. Through fitting of the LEIA spectrum of the Virgo Cluster, we obtained a temperature of $2.1^{+0.3}_{-0.1}$ keV, which is consistent with the XMM-Newton results ($\sim$2.3 keV). Above 1.6 keV, the spectrum is dominated by the X-ray background. In summary, this study validates LEIA's extended source imaging and spectral resolution capabilities for the first time.

Galaxy evolution emerges from the balance between cosmic gas accretion, fueling star formation, and supernova (SN) feedback, regulating the metal enrichment. Hence, the stellar mass ($M_*$) - gas metallicity relation (MZR) is key to understand the physics of galaxies. High-quality JWST data enable accurate measurements of the MZR up to redshift z=10. Our aims are to understand the observed MZR, its connection with the star formation rate (SFR), the role played by SFR stochasticity, and how it is regulated by SN feedback. We compare the MZR from the JADES, CEERS, and UNCOVER surveys, which comprise about 180 galaxies at $z=3-10$ with $10^6<M_*/M_\odot<10^{10}$, with 200 galaxies from the SERRA cosmological simulations. To interpret the MZR, we develop a minimal model for galaxy evolution that includes: cosmic accretion modulated with an amplitude $A_{100}$ on 100 Myr; a time delay $t_d$ between SFR and SN; SN-driven outflows with a varying mass loading factor $\epsilon_{SN}$. Using our minimal model, we find the observed mean MZR is reproduced by weak outflows ($\epsilon_{SN}=1/4$), in line with findings from JADES. Matching the observed MZR dispersion requires $t_d=20$ Myr and a $A_{100}=1/3$ modulation of the accretion rate. Successful models have low stochasticity ($\sigma_{SFR}=0.2$), yielding a MZR dispersion of $\sigma_{Z}=0.2$. Such values are close but lower than SERRA predictions ($\sigma_{SFR}=0.24$, $\sigma_{Z}=0.3$), clarifying why SERRA show no clear MZR trend and some tension with the observations. As the MZR is very sensitive to SFR stochasticity, models predicting high r.m.s. values ($\sigma_{SFR}=0.5$) result in a ``chemical chaos'' (i.e. $\sigma_{Z}=1.4$), virtually destroying the MZR. As a consequence, invoking a highly stochastic SFR ($\sigma_{SFR}=0.8$) to explain the overabundance of bright, super-early galaxies leads to inconsistencies with the observed MZR.

Studies of galaxy protoclusters yield insights into galaxy cluster formation complementary to those obtained via `archaeological' studies of present-day galaxy clusters. Submillimetre-selected galaxies (SMGs) are one class of sources used to find high-redshift protoclusters. However, due to the rarity of protoclusters (and thus the large simulation volume required) and the complexity of modeling dust emission from galaxies, the relationship between SMGs and protoclusters has not been adequately addressed in the theoretical literature. In this work, we apply the L-GALAXIES semi-analytic model (SAM) to the Millennium N-body simulation. We assign submillimetre (submm) flux densities to the model galaxies using a scaling relation from previous work, in which dust radiative transfer was performed on high-resolution galaxy zoom simulations. We find that the fraction of model galaxies that are submm-bright is higher in protocluster cores than in both protocluster `outskirts' and the field; the fractions for the latter two are similar. This excess is not driven by an enhanced starburst frequency. Instead, the primary reason is that overdense environments have a relative overdensity of high-mass halos and thus `oversample' the high-mass end of the star formation main sequence relative to less-dense environments. The fraction of SMGs that are optically bright is dependent on stellar mass and redshift but independent of environment. The fraction of galaxies for which the majority of star formation is dust-obscured is higher in protocluster cores, primarily due to the dust-obscured fraction being correlated with stellar mass. Our results can be used to guide and interpret multi-wavelength studies of galaxy populations in protoclusters.

Submillimetre galaxies (SMGs) are some of the most extreme star-forming systems in the Universe, whose place in the framework of galaxy evolution is as yet uncertain. It has been hypothesised that SMGs are progenitors of local early-type galaxies, requiring that SMGs generally reside in galaxy cluster progenitors at high redshift. We test this hypothesis and explore SMG environments using a narrowband VLT/HAWK-I+GRAAL study of H$\alpha$ and [OIII] emitters around an unbiased sample of three ALMA-identified and spectroscopically-confirmed SMGs at $z \sim 2.3$ and $z \sim 3.3$, where these SMGs were selected solely on spectroscopic redshift. Comparing with blank-field observations at similar epochs, we find that one of the three SMGs lies in an overdensity of emission-line sources on the $\sim4$ Mpc scale of the HAWK-I field of view, with overdensity parameter $\delta_{g} = 2.6^{+1.4}_{-1.2}$. A second SMG is significantly overdense only on $\lesssim 1.6$ Mpc scales and the final SMG is consistent with residing in a blank field environment. The total masses of the two overdensities are estimated to be $\log(M_{h}/{\rm M}_{\odot}) =$12.1--14.4, leading to present-day masses of $\log(M_{h,z=0}/{\rm M}_{\odot}) =$12.9--15.9. These results imply that SMGs occupy a range of environments, from overdense protoclusters or protogroups to the blank field, suggesting that while some SMGs are strong candidates for the progenitors of massive elliptical galaxies in clusters, this may not be their only possible evolutionary pathway.

It is commonly believed that our own Milky Way is on a collision course with the neighbouring Andromeda galaxy. As a result of their merger, predicted in around five billion years, the two large spiral galaxies that define the present Local Group would form a new elliptical galaxy. Here we consider the latest and most accurate observations by the Gaia and Hubble space telescopes, along with recent consensus mass estimates to derive possible future scenarios and identify the major sources of uncertainty in the evolution of the Local Group over the next 10 billion years. We find that the next most massive Local Group member galaxies -- namely, M33 and the Large Magellanic Cloud -- distinctly and radically affect the Milky Way - Andromeda orbit. While including M33 increases the merger probability, the orbit of the Large Magellanic Cloud runs perpendicular to the Milky Way - Andromeda orbit and makes their merger less likely. In the full system, we find that uncertainties in the present positions, motions, and masses of all galaxies leave room for drastically different outcomes, and a probability of close to 50% that there is no Milky Way - Andromeda merger during the next 10 billion years.

The Hubble Tension, a >5 sigma discrepancy between direct and indirect measurements of the Hubble constant (H0), has persisted for a decade and motivated intense scrutiny of the paths used to infer H0. Comparing independently-derived distances for a set of galaxies with different standard candles, such as the tip of the red giant branch (TRGB) and Cepheid variables, can test for systematics in the middle rung of the distance ladder. The I band is the preferred filter for measuring the TRGB due to constancy with color, a result of low sensitivity to population differences in age and metallicity supported by stellar models. We use James Webb Space Telescope (JWST) observations with the maser host NGC 4258 as our geometric anchor to measure I-band (F090W vs F090W-F150W) TRGB distances to 7 hosts of 9 Type Ia supernovae (SNe Ia) within 27 Mpc: NGC 1448, NGC 1559, NGC 2525, NGC 3370, NGC 3447, NGC 5584, and NGC 5643. We compare these with Hubble Space Telescope (HST) Cepheid-based relative distance moduli for the same galaxies and anchor. We find no evidence of a difference between their weighted means, 0.01 +/- 0.04 (stat) +/- 0.04 (sys) mag. We produce fourteen variants of the TRGB analysis, altering the smoothing level and color range used to measure the tips to explore their impact. For some hosts, this changes the identification of the strongest peak, but this causes little change to the sample mean difference producing a full range of 0.01 to 0.03 mag, all consistent at 1 sigma with no difference. The result matches past comparisons of I-band TRGB and Cepheids when both use HST. SNe and anchor samples observed with JWST are too small to yield a measure of H0 that is competitive with the HST sample of 42 SNe Ia and 4 anchors; however, they already provide a vital systematic crosscheck to HST measurements of the distance ladder.

Observations of low-mass protostellar systems show evidence of rich complex organic chemistry. Their low luminosity, however, makes determining abundance distributions of complex organic molecules (COMs) within the water snowline challenging. However, the excitation conditions sampled by differing molecular distributions may produce substantive changes in the resulting emission. Thus, molecular excitation may recover spatial information from spatially unresolved data. By analyzing spatially-unresolved NOrthern Extended Millimeter Array (NOEMA) observations of CH$_3$OH and CH$_3$CN, we aim to determine if CH$_3$OH and CH$_3$CN are distributed differently in the protostellar disk around HOPS-370, a highly-luminous intermediate mass protostar. Rotational diagram analysis of CH$_3$OH and CH$_3$CN yields rotational temperatures of $198 \pm 1.2$ K and $448 \pm 19$ K, respectively, suggesting the two molecules have different spatial distributions. Source-specific 3D LTE radiative transfer models are used to constrain the spatial distribution of CH$_3$OH and CH$_3$CN within the disk. A uniform distribution with an abundance of $4\times10^{-8}$ reproduces the CH$_3$OH observations. In contrast, the spatial distribution of CH$_3$CN needs to be either more compact (within $\sim120$ au versus $\sim240$ au for CH$_3$OH) or exhibiting a factor of $\gtrsim 15$ increase in abundance in the inner $\sim55$ au. A possible explanation for the difference in spatial abundance distributions of CH$_3$OH and CH$_3$CN is carbon-grain sublimation.

The location of the baryon acoustic oscillation (BAO) feature in the two-point correlation function (2PCF) of matter produces a standard ruler that is useful for the measurement of the expansion history of the Universe. Inspired by the possibility of reconstructing the positions of protohalos in the initial density field with a novel method rooted in optimal transport theory, we revisit the BAO signal in the protohalo correlation function. Our work examines the performance of a template 2PCF built on a tracer bias relation that includes scale dependence -- a term that can be motivated by peaks theory or a general bias expansion. Working in protohalos, halos, and the linear combination of the protohalo and matter fields that is motivated by the continuity equation, we demonstrate that this model accurately captures the shape of the BAO feature and improves the precision of the BAO scale measurement relative to a model that does not include scale-dependent bias by 47% in protohalos, 15% in halos, and 14% in the linear combination of the protohalo and matter fields. Allowing for scale dependence does not appear to introduce any shift in the BAO feature. The precision of the BAO distance scale estimate is highest with the linear combination of the protohalo and matter fields, which offers a factor of 3.5 improvement over Eulerian-space measurements and a factor of 4-8 improvement over the estimate made with protohalos alone.

Mass estimates of black holes (BHs) in the centers of Active Galactic Nuclei (AGN) often rely on the radius-luminosity relation. However, this relation, usually probed by reverberation mapping (RM), is poorly constrained in the high-luminosity and high-redshift ends due to the very long expected lag times. Multiply imaged AGN may thus offer a unique opportunity to explore the radius-luminosity relation at these ends. In addition to comprising several magnified images which enable a more efficient light-curve sampling, the time delay between multiple images of strongly lensed quasars can also aid in making such RM measurements feasible on reasonable timescales: If the time delay is, for example, of the order of the expected time lag, changes in the emission lines in the leading image can be observed around the same time as the changes in the continuum in the trailing image. In this work we probe the typical time-delay distribution in galaxy-cluster lenses and estimate the number of both high-mass ($\sim10^9-10^{10}$ M$_{\odot}$), and high-redshift ($z\gtrsim4-12$) quasars that are expected to be strongly lensed by clusters. We find that less than one, very massive (and luminous, L$_{UV}>10^{46.5}$ erg s$^{-1}$) multiply imaged quasar is expected across the sky down to 30 AB mag. Nonetheless, up to several dozen thousand M$_{BH}\sim10^{6}$-$10^{8}$ M$_{\odot}$ broad-line AGN at $z>4$ should be multiply imaged by galaxy clusters and detectable with JWST, hundreds with $\textit{Euclid}$ and several thousands with the $\textit{Roman}$ Space Telescope, across the whole sky. These could supply an important calibration for the BH mass scaling in the early Universe.

The existence of compact stellar remnants in the mass range $2-5\,M_{\odot}$ has long been debated. This so-called lower mass gap was initially suggested by the lack of low-mass X-ray binary observations with accretors in this mass range, but it has recently been called into question following newer observations, including a lower-mass-gap candidate with a millisecond pulsar companion in the dense globular cluster NGC~1851. Here we model NGC~1851 with a grid of similar dense star clusters utilizing the state-of-the-art Monte Carlo $N$-body code \texttt{CMC}, and we specifically study the formation of lower-mass-gap black holes. We demonstrate that both massive star evolution and dynamical interactions can contribute to forming lower-mass-gap black holes. In general, the collapse of massive remnants formed through mergers of neutron stars or massive white dwarfs produces the largest number of lower-mass-gap black holes among all formation channels. However, in more massive clusters, supernova core collapse can contribute comparable numbers. Our NGC~1851-like models can reproduce millisecond pulsar -- lower-mass-gap black hole binaries similar to the observed system. Additionally, the lower-mass-gap black holes can also become components of dynamically assembled binaries, and some will be in merging black hole--neutron star systems similar to the recently detected gravitational wave source GW230529. However, the corresponding merger rate is probably $\lesssim 1~{\rm Gpc^{-3}\,yr^{-1}}$.

Cataclysmic variables (CVs) are compact binary systems in which a white dwarf accretes matter from a Roche-lobe-filling companion star. In this study, we searched for new CVs in the Milky Way in the Chandra Source Catalog v2.0, cross-matched with Gaia Data Release 3 (DR3). We identified new CV candidates by combining X-ray and optical data in a color-color diagram called the ``X-ray Main Sequence". We used two different cuts in this diagram to compile pure and optically variable samples of CV candidates. We undertook optical spectroscopic follow-up observations with the Keck and Palomar Observatories to confirm the nature of these sources. We assembled a sample of 25,887 Galactic X-ray sources and found 14 new CV candidates. Seven objects show X-ray and/or optical variability. All sources show X-ray luminosity in the $\rm 10^{29}-10^{32}$ $\rm erg\ s^{-1}$ range, and their X-ray spectra can be approximated by a power-law model with photon indices in the $\rm \Gamma \sim 1-3$ range or an optically thin thermal emission model in the $\rm kT \sim 1-70$ keV range. We spectroscopically confirmed four CVs, discovering two new polars, one low accretion rate polar and a WZ~Sge-like low accretion rate CV. X-ray and optical properties of the other 9 objects suggest that they are also CVs (likely magnetic or dwarf novae), and one other object could be an eclipsing binary, but revealing their true nature requires further observations. These results show that a joint X-ray and optical analysis can be a powerful tool for finding new CVs in large X-ray and optical catalogs. X-ray observations such as those by Chandra are particularly efficient at discovering magnetic and low accretion rate CVs, which could be missed by purely optical surveys.

We investigate the `intrinsic' H$\alpha$ EW distributions of $z \sim 0.4 - 2.2$ narrowband-selected H$\alpha$ samples from HiZELS and DAWN using a forward modeling approach. We find an EW - stellar mass anti-correlation with steepening slopes $-0.18\pm0.03$ to $-0.24^{+0.06}_{-0.08}$ at $z \sim 0.4$ and $z\sim 2.2$, respectively. Typical EW increases as $(1+z)^{1.78^{+0.22}_{-0.23}}$ for a $10^{10}$ M$_\odot$ emitter from $15^{+2.4}_{-2.3}$Å ($z \sim 0.4$) to $67.7^{+10.4}_{-10.0}$Å ($z \sim 2.2$) and is steeper with decreasing stellar mass highlighting the high EW nature of low-mass high-$z$ systems. We model this redshift evolving anti-correlation, $W_0(M,z)$, and find it produces H$\alpha$ luminosity and SFR functions strongly consistent with observations validating the model and allowing us to use $W_0(M,z)$ to investigate the relative contribution of H$\alpha$ emitters towards cosmic SF. We find EW$_0 > 200$ Å emitters contribute significantly to cosmic SF activity at $z \sim 1.5 - 2$ making up $\sim 40$% of total SF consistent with sSFR $> 10^{-8.5}$ yr$^{-1}$ ($\sim 45 - 55$%). Overall, this highlights the importance of high EW systems at high-$z$. Our $W_0(M,z)$ model also reproduces the cosmic sSFR evolution found in simulations and observations and show that tension between the two can simply arise from selection effects in observations. Lastly, we forecast Roman and Euclid grism surveys using $W_0(M,z)$ including observational efficiency and limiting resolution effects where we predict $\sim 24000$ and $\sim 30000$ $0.5 < z < 1.9$ H$\alpha$ emitters per deg$^{-2}$, respectively, down to $>5\times10^{-17}$ erg s$^{-1}$ cm$^{-2}$ including $10^{7.2 - 8}$ M$_\odot$ galaxies at $z > 1$ with EW$_0 >1000$Å. Both Roman and Euclid will enable us to observe with unprecedented detail some of the most bursty/high EW, low-mass star-forming galaxies near cosmic noon.

This research introduces an innovative application of physics-informed neural networks (PINNs) to tackle the intricate challenges of radiative transfer (RT) modeling in exoplanetary atmospheres, with a special focus on efficiently handling scattering phenomena. Traditional RT models often simplify scattering as absorption, leading to inaccuracies. Our approach utilizes PINNs, noted for their ability to incorporate the governing differential equations of RT directly into their loss function, thus offering a more precise yet potentially fast modeling technique. The core of our method involves the development of a parameterized PINN tailored for a modified RT equation, enhancing its adaptability to various atmospheric scenarios. We focus on RT in transiting exoplanet atmospheres using a simplified 1D isothermal model with pressure-dependent coefficients for absorption and Rayleigh scattering. In scenarios of pure absorption, the PINN demonstrates its effectiveness in predicting transmission spectra for diverse absorption profiles. For Rayleigh scattering, the network successfully computes the RT equation, addressing both direct and diffuse stellar light components. While our preliminary results with simplified models are promising, indicating the potential of PINNs in improving RT calculations, we acknowledge the errors stemming from our approximations as well as the challenges in applying this technique to more complex atmospheric conditions. Specifically, extending our approach to atmospheres with intricate temperature-pressure profiles and varying scattering properties, such as those introduced by clouds and hazes, remains a significant area for future development.

The extremes of Active Galactic Nuclei (AGN) variability offer valuable new insights into the drivers and physics of AGN. We discuss some of the most extreme cases of AGN variability; the highest amplitudes, deep minima states, extreme spectral states, Seyfert-type changes, and semi-periodic signals, including new X-ray observations. The properties of changing-look (CL) AGN are briefly reviewed and a classification scheme is proposed which encompasses the variety of CL phenomena; distinguishing slow and fast events, repeat events, and frozen-look AGN which do not show any emission-line response. Long-term light curves that are densely covered over multiple years, along with follow-up spectroscopy, are utilized to gain insight into the underlying variability mechanisms including accretion disk and broad-line region physics. Remarkable differences are seen, for instance, in the optical spectral response to extreme outbursts, implying distinct intrinsic variability mechanisms. Furthermore, we discuss methods for distinguishing between CL AGN and CL look-alike events (tidal disruption events or supernovae in dense media). Finally, semi-periodic light curve variability is addressed and the latest multiwavelength (MWL) light curve of the binary supermassive black hole (SMBH) candidate OJ 287 from the MOMO project is presented. Recent results from that project have clearly established the need for new binary SMBH modelling matching the tight new constraints from observations, including the measurement of a low (primary) SMBH mass of ~10^8 Msun which also implies that OJ 287 is no longer in the regime of near-future pulsar timing arrays.

Outflows driven by active galactic nuclei (AGN) are seen in numerous compact sources; however, it has remained unclear how to distinguish between the driving mechanisms, such as winds and jets. Therefore, our study aims to offer observational insights from simulations to aid in this distinction. Specifically, in this paper, we investigate the evolution of wide-angled, moderately relativistic, magnetized winds and analyze their non-thermal radio emission and polarization properties. We find that the evolution of winds varies depending on factors such as power, density, and opening angle, which in turn influence their observable characteristics. Additionally, different viewing angles can lead to varying observations. Furthermore, we note distinctions in the evolution of winds compared to jets, resulting in disparities in their observable features. Jets typically exhibit a thin spine and hotspot(s). Winds manifest broader spines or an "hourglass-shaped" bright emission in the cocoon, which are capped by bright arcs. Both display high polarization coinciding with the bright spine and hotspots/arcs, although these regions are relatively compact and localized in jets when compared to winds. We emphasize the importance of high resolution, as we demonstrate that emission features from both jets and winds can become indistinguishable at lower resolutions. The distribution of polarization is largely unaffected by resolution, though lower polarization becomes more noticeable when the resolution is decreased.

Dynamical friction and stellar orbital motion in spiral galaxies with dark matter composed of ultralight bosons in the state of {rotating} Bose-Einstein condensate (BEC) are studied. It is found that the dynamical friction force is significantly affected by the topological charge of the vortex structure of the BEC core with the strongest effect at distances near the galactic center. It is also shown that the ultralight dark matter self-interaction plays an important role in studying the dynamical friction.

Stellar theory enables us to understand the properties of stars at different stages of their evolution, and contributes to other fields of astrophysics such as galactic and exoplanet studies. Assessing the accuracy of stellar theories necessitates high precision, model-independent measurements of the properties of real stars, such as those obtainable for the components of double lined eclipsing binaries (DLEBs), while asteroseismology offers probing power of the stellar interior if one or both components pulsate. KIC 4851217 is a DLEB containing two late A-type stars and exhibits pulsations of the $\delta$ Scuti type. By analysing high resolution HERMES and moderate resolution ISIS spectra, jointly with Kepler and TESS light curves, we measured the masses, radii and effective temperatures of the components to precisions of ~0.5, ~1.1 and ~1 per cent, respectively. We additionally report the discovery and characterisation of a tertiary M-dwarf companion. Models of the system's spectral energy distribution agree with an age of 0.82 Gyr, with the more massive and larger secondary component near the end of the main sequence lifetime. An examination of the pulsating component's pulsation frequencies reveals 39 pulsation multiplets that are split by the orbital frequency. For most of these, it is evident that the pulsation axes have been tilted into the orbital plane. This makes KIC 4851217 a tidally tilted pulsator (TTP). This precisely characterized $\delta$ Scuti DLEB is an ideal candidate for advancing intermediate-mass stellar theory, contributing to our understanding of hierarchichal systems as well as to the topic of TTPs.

Recent studies have identified star clusters with multiple components based on accurate spatial distributions and/or proper motions from Gaia DR3, utilising diverse diagnostics to improve our understanding of subgroup evolution. These findings motivated us to search for subgroups among the objects examined in our previous work, which employed fractal statistics. The present study considers seven open clusters that exhibit significant dispersion in age and/or proper motion distributions, suggesting that they would not be single clusters. For characterizing the stellar groups, we calculate the membership probability using Bayesian multi-dimensional analysis by fitting the observed proper motion distribution of the candidates. A probability distribution is also used to determine the distance of the cluster, which is obtained from the mean value of the distance modes. The photometry from Gaia DR3 is compared with evolutionary models to estimate the cluster age and total mass. In our sample, double components are found only for Markarian 38 and NGC 2659. The other five clusters are confirmed as being single. The structural parameters, such as Q, Lambda_MSR and Sigma_LDR are compared with results from N-body simulations to investigate how the morphology of the stellar clustering evolves. The new results, for a more complete sample of cluster members, provide a better definition of the distribution type (central concentration or substructured region) inferred from the m-s plot.

The degree of coupling between dust particles and their surrounding gas in protoplanetary disks is quantified by the dimensionless Stokes number. The Stokes number (St) governs particle size and spatial distributions, in turn establishing the dominant mode of planetary accretion in different disk regions. In this paper, we model the characteristic St of particles across time in disks evolving under both turbulent viscosity and magnetohydrodynamic (MHD) disk winds. In both turbulence- and wind-dominated disks, we find that collisional fragmentation is the limiting mechanism of particle growth, and the water-ice sublimation line constitutes a critical transition point between dust settling, drift, and size regimes. The St dichotomy across the ice-line translates to distinct planet formation pathways between the inner and outer disk. While pebble accretion proceeds slowly for rocky embryos within the ice-line (across most of parameter space), it does so rapidly for volatile-rich embryos beyond it, allowing for the growth of giant planet cores before disk dissipation. Through simulations of rocky planet growth, we evaluate the competition between pebble accretion and classical pairwise collisions between planetesimals. We conclude that the dominance of pebble accretion can only be realized in disks that are driven by MHD winds, slow-evolving, and devoid of pressure maxima that may concentrate solids and give rise of planetesimal rings. Such disks are extremely quiescent, with Shakura-Sunyaev turbulence parameters $\alpha_{\nu} \sim 10^{-4}$. We conclude that for most of parameter space corresponding to values of $\alpha_{\nu}$ reflected in observations of protoplanetary disks ($\gtrsim 10^{-4}$), pairwise collisions constitute the dominant pathway of rocky planet accretion. Our results are discussed in the context of super-Earth origins and Earth's accretion history.

We present the first results of the holographic beam mapping program for the Canadian Hydrogen Intensity Mapping Experiment (CHIME). We describe the implementation of the holographic technique as adapted for CHIME, and introduce the processing pipeline which prepares the raw holographic timestreams for analysis of beam features. We use data from six bright sources across the full 400-800\,MHz observing band of CHIME to provide measurements of the co-polar and cross-polar beam response of CHIME in both amplitude and phase for the 1024 dual-polarized feeds instrumented on CHIME. In addition, we present comparisons with independent probes of the CHIME beam which indicate the presence of polarized beam leakage in CHIME. Holographic measurements of the CHIME beam have already been applied in science with CHIME, e.g. in estimating detection significance of far sidelobe FRBs, and in validating the beam models used for CHIME's first detections of \tcm emission (in cross-correlation with measurements of large-scale structure from galaxy surveys and the Lyman-$\alpha$ forest). Measurements presented in this paper, and future holographic results, will provide a unique data set to characterize the CHIME beam and improve the experiment's prospects for a detection of BAO.

As part of the Bristol PROVE mission, a nano satellite in low Earth orbit will be required to track a ground based target during a 400 second flyover. This requires agile attitude control that will be achieved using a system of flywheels. To calculate the necessary torque from these flywheels, a controller was designed. Using newly derived equations of motion for the system, an expression to optimise the gains was produced. With this controller, simulations were run to evaluate the largest causes of error in target pointing. Disturbance torques were safely handled by the controller, but led to a 12% increase in wheel speeds, reaching 8325 rpm. This higher speed led to an increased gyroscopic torque, reaching 10^-7 Nm in the worst case. However since the flywheels can deliver 10^-5 Nm of torque, the controller could also correct for this. Hardware performance was then varied to assess the effect of each component on pointing accuracy. Attitude sensor noise was found to increase pointing error by 1.9 degrees in the worst case. Minimum performance requirements were then determined for each component in order to maintain an acceptable pointing accuracy.

X-ray observations provide essential and valuable insights into the acceleration and propagation of non-thermal electrons during solar flares. Improved X-ray spectral analysis requires a deeper understanding of the dynamics of energetic electrons. Previous studies have demonstrated that the dynamics of accelerated electrons of a few thermal speeds are more complex. To better describe the energetic electrons after injection, a model considering energy diffusion and thermalization effects in flare conditions (warm-target model) has recently been developed for Hard X-ray spectral analysis. This model has demonstrated how the low-energy cut-off, which can hardly be constrained in cold-target modeling, can be determined. However, the power-law form may not be the most suitable representation of injected electrons. The kappa distribution, which is proposed as a physical consequence of electron acceleration, has shown successful application in RHESSI spectral analysis. In this study, we employ the kappa-form injected electrons in the warm-target model to analyze two M-class flares, observed by RHESSI and STIX, respectively. The best-fit results show that the kappa-form energetic electron spectrum generates lower non-thermal energy when producing a similar photon spectrum in the fit range compared to the power-law form. We also demonstrated that the fit parameters associated with kappa-form electron spectrum can be well determined with small fit uncertainty. Further, the kappa distribution, which covers the entire electron energy range, enables the determination of key electron properties such as total electron number density and average energy in the flare site, providing valuable information on electron acceleration processes.

The fraction of low-mass galaxies hosting an intermediate-mass black hole (IMBH, with masses $M_{\rm BH} \approx 10^2-10^5$ M$_\odot$), is sensitive to how black hole seeds formed in the early Universe but is observationally still unconstrained. In this paper, we assemble a sample of dwarf galaxies within 10 Mpc hosting bright nuclear star clusters (NSCs) that could host IMBHs. For a subset of them, we use their observed surface brightness from {\it Hubble Space Telescope} (\hst) images, an assumed synthetic spectrum of their stellar population, Jeans Anisotropic Model (JAM) of the stellar dynamics, and the {\tt HSIM} simulator software to create mock observations with the High Angular Resolution Monolithic Optical and Near-infrared Integral (HARMONI) field spectrograph for the Extremely Large Telescope (ELT). We analyze the simulated data cube like real data, using JAM to infer the IMBH mass and its error in a Bayesian framework. Our simulations show that the ELT/HARMONI instrument can clearly detect the existence of IMBH demographics in NSCs down to a mass of about 0.5\% of the NSC.

As a key science project of the Square Kilometre Array (SKA), the discovery and timing observations of radio pulsars in the Galactic Center would provide high-precision measurements of the spacetime around the supermassive black hole, Sagittarius A* (Sgr A*), and initiate novel tests of general relativity. The spin of Sgr A* could be measured with a relative error of $\lesssim 1\%$ by timing one pulsar with timing precision that is achievable for the SKA. However, the real measurements depend on the discovery of a pulsar in a very compact orbit, $P_b\lesssim0.5\,{\rm yr}$. Here for the first time we propose and investigate the possibility of probing the spin of Sgr A* with two or more pulsars that are in orbits with larger orbital periods, $P_b\sim 2- 5\,{\rm yr}$, which represents a more realistic situation from population estimates. We develop a novel method for directly determining the spin of Sgr A* from the timing observables of two pulsars and it can be readily extended for combining more pulsars. With extensive mock data simulations, we show that combining a second pulsar improves the spin measurement by $2-3$ orders of magnitude in some situations, which is comparable to timing a pulsar in a very tight orbit.

Morphological evolution of expanding shells of fast-mode magnetohydrodynamic (MHD) waves through an inhomogeneous ISM is investigated in order to qualitatively understand the complicated morphology of shell-type supernova remnants (SNR). Interstellar clouds with high Alfvén velocity act as concave lenses to diverge the MHD waves, while those with slow Alfvén velocity act as convex lenses to converge the waves to the focal points. By combination of various types of clouds and fluctuations with different Alfvén velocities, sizes, or wavelengths, the MHD-wave shells attain various morphological structures, exhibiting filaments, arcs, loops, holes, and focal strings, mimicking old and deformed SNRs.

Exoplanetary science is increasingly prioritizing efforts toward direct imaging of planetary systems, with emphasis on those that may enable the detection and characterization of potentially habitable exoplanets. The recent 2020 Astronomy and Astrophysics decadal survey recommended the development of a space-based direct imaging mission that has subsequently been referred to as the Habitable Worlds Observatory (HWO). A fundamental challenge in the preparatory work for the HWO search for exo-Earths is the selection of suitable stellar targets. Much of the prior efforts regarding the HWO targets has occurred within the context of exoplanet surveys that have characterized the stellar properties for the nearest stars. The preliminary input catalog for HWO consists of 164 stars, of which 30 are known exoplanet hosts to 70 planets. Here, we provide a dynamical analysis for these 30 systems, injecting a terrestrial planet mass into the Habitable Zone (HZ) and determining the constraints on stable orbit locations due to the influence of the known planets. For each system, we calculate the percentage of the HZ that is dynamically viable for the potential presence of a terrestrial planet, providing an additional metric for inclusion of the stars within the HWO target list. Our analysis shows that, for 11 of the systems, less than 50% of the HZ is dynamically viable, primarily due to the presence of giant planets whose orbits pass near or through the HZ. These results demonstrate the impact that known system architectures can have on direct imaging target selection and overall system habitability.

Precise cosmological measurements are essential for understanding the evolution of the universe and the nature of dark energy. The Five-hundred-meter Aperture Spherical Telescope (FAST), the most sensitive single-dish radio telescope, has the potential to provide the precise cosmological measurements through neutral hydrogen 21-cm intensity mapping sky survey. This paper primarily explores the potential of technological upgrades for FAST in cosmology. The most crucial upgrade begins with equipping FAST with a wide-band receiver ($0 < z < 2.5$). This upgrade can enable FAST to achieve higher precision in cosmological parameter estimation than the Square Kilometre Array Phase-1 Mid frequency. On this basis, expanding to a FAST array (FASTA) consisting of six identical FAST would offer significant improvements in precision compared to FAST. Additionally, compared with the current results from the data combination of cosmic microwave background, baryon acoustic oscillations (optical galaxy surveys), and type Ia supernovae, FASTA can provide comparable constraints. Specifically, for the dark-energy equation-of-state parameters, FASTA can achieve $\sigma(w_0) = 0.09$ and $\sigma(w_a) = 0.33$.

In this paper, we intend to investigate the dynamics of the Circular Restricted Three-Body Problem. Here we assumed the primaries as the source of radiation and have variable mass. The gravitational perturbation from disk-like structure are also considered in this study. There exist five equilibrium points in this system. By considering the combined effect from disk-like structure and the mass transfer, we found that the classical collinear equilibrium points depart from x-axis. Meanwhile, this combined effect also breaks the symmetry of tringular equlibrium point positions. We noted that the quasi-equilibrium points are unstable whereas the triangular equilibrium points are stable if the mass ratio $\mu$ smaller than critical mass $\mu_c$. It shows that the stability of triangular equilibrium points depends on time.

Based on the magnetization, an accretion disk with large-scale magnetic field can be separated into either standard and normal evolution (SANE) or magnetically arrested disk (MAD), which are difficult to identify from observations. It is still unclear whether all the radio-loud active galactic nuclei (RLAGNs) with a thin disk and strong radio emissions contain a MAD. We investigate this issue by utilizing the 3CRR catalog. We compile a sample of 35 quasars and 14 high-excitation radio galaxies powered by a thin accretion disk. In order to consistently compare with the MAD sample given by Li et al. (2022), the optical-UV emissions of our sample are all detected by the Hubble Space Telescope (HST). It is found that the average X-ray luminosity ($L_{\rm X}$) of our sample is about 5.0 times higher than that of radio-quiet AGNs (RQAGNs) with matching optical-UV luminosity ($L_{\rm UV}$), in general accord with the factor of 4.5 times in MAD sample within the uncertainty. The relationship between radio (5~GHz) and X-ray (2 keV) luminosities in the 3CRR sources is also found to be consistent with the MAD sample. Furthermore, the jet efficiencies of 3CRR sources are consistent with those from the GRMHD simulations of MAD. Therefore, we suggest that probably all the quasars and at least a fraction of high-excitation radio galaxies in the 3CRR catalog, and perhaps all the RLAGNs with strong radio emissions contain a MAD.

The first observation of a gravitational wave (GW) and a short gamma-ray burst (sGRB) emitted by the same binary neutron star (BNS) merger officially opened the field of GW multi-messenger astronomy. In this paper, we define and address \textit{lagging sirens}, a new class of multi-messenger BNSs for which associated GWs and sGRBs are observed without the identification of their host galaxy. We propose a new methodology to use the observed time delay of these sources to constrain the speed of gravity that is, the propagation speed of gravitational waves, the Hubble constant and the prompt time delay distribution between GWs and sGRBs, even though a direct redshift estimation from the host galaxy is unavailable. Our method exploits the intrinsic relation between GWs and sGRBs observed and prompt time delays to obtain a statistical redshift measure for the cosmological sources. We show that this technique can be used to infer the Hubble constant at the $10\%$~level of precision with future-generation GW detectors such as the Einstein Telescope and only 100 observations of this kind. The novel procedure that we propose has systematics that differ completely from the ones of previous GW methods for cosmology. Additionally, we demonstrate for the first time that the speed of gravity and the distribution of the prompt time-delays between GWs and sGRBs can be inferred conjointly with less than 10 sources even with current GW detector sensitivities.

We present the identification of changing-look active galactic nuclei (CL-AGNs) from the Dark Energy Spectroscopic Instrument First Data Release and Sloan Digital Sky Survey Data Release 16 at z \leq 0.9. To confirm the CL-AGNs, we utilize spectral flux calibration assessment via an [O\,{\sc iii}]-based calibration, pseudo-photometry examination, and visual inspection. This rigorous selection process allows us to compile a statistical catalog of 561 CL-AGNs, encompassing 527 $\rm H\beta$, 149$\rm H\alpha$, and 129 Mg II CL behaviors. In this sample, we find 1) a 283:278 ratio of turn-on to turn-off CL-AGNs. 2) the critical value for CL events is confirmed around Eddington ratio \sim 0.01. 3) a strong correlation between the change in the luminosity of the broad emission lines (BEL) and variation in the continuum luminosity, with Mg II and $\rm H\beta$ displaying similar responses during CL phases. 4) the Baldwin-Phillips-Terlevich diagram for CL-AGNs shows no statistically difference from the general AGN catalog. 5) five CL-AGNs are associated with asymmetrical mid-infrared flares, possibly linked to tidal disruption events. Given the large CL-AGNs and the stochastic sampling of spectra, we propose that some CL events are inherently due to typical AGN variability during low accretion rates, particularly for CL events of the singular BEL. Finally, we introduce a Peculiar CL phase, characterized by a gradual decline over decades in the light curve and the complete disappearance of entire BEL in faint spectra, indicative of a real transition in the accretion disk.

Using \textit{Swift} Burst Alert Telescope (BAT) event-mode data during Gamma Ray Burst (GRB) occurrences, we conducted spectral analysis for the Crab system. From 38 good observations, which spans over a period of 18 years from 2006 to 2023, we found that the Crab's X-ray flux does not only flicker, but also significantly anti-correlates to its spectral power-law photon index. Since emission contribution of the Crab pulsar in this energy range is small, this anti-correlation is mainly about the emission of the Crab nebula. We suggest that this anti-correlation is an observational supporting evidence for the long-standing notion that the nebula emission is due to synchrotron radiation of shocked pulsar winds in the nebula.

Neutral hydrogen (HI) is the primary component of the cool interstellar medium (ISM) and is the reservoir of fuel for star formation. Owing to the sensitivity of existing radio telescopes, our understanding of the evolution of the ISM in galaxies remains limited, as it is based on only a few hundred galaxies detected in HI beyond the local Universe. With the high sensitivity of the Five-hundred-meter Aperture Spherical radio Telescope (FAST), we carried out a blind HI search, the FAST Ultra-Deep Survey (FUDS), which extends to redshifts up to 0.42 and a sensitivity of 50 $\rm \mu Jy \cdot beam^{-1}$. Here, we report the first discovery of six galaxies in HI at $z>0.38$. For these galaxies, the FAST angular resolution of $\sim\,4'$ corresponds to a mean linear size of $\sim1.3\,h_{70}^{-1}\,$Mpc. These galaxies are among the most distant HI emission detections known, with one having the most massive HI content ($10^{10.93 \pm 0.04}~h_{70}^{-2}\, \rm M_\odot$). Using recent data from the DESI survey, and new observations with the Hale, BTA, and Keck telescopes, optical counterparts are detected for all galaxies within the 3-$\sigma$ positional uncertainty ($0.5\,h_{70}^{-1}\,$Mpc) and $\rm 200\,km \cdot s^{-1}$ in recession velocity. Assuming that the dominant source of HI is the identified optical counterpart, we find an evidence of evolution in the HI content of galaxies over the last 4.2 Gyr. Our new high-redshift HI galaxy sample provides the opportunity to better investigate the evolution of cool gas in galaxies. A larger sample size in the future will allow us to refine our knowledge of the formation and evolution of galaxies.

The AGORA Cosmorun (arXiv:2106.09738) is a set of hydrodynamical cosmological zoom-in simulations carried out within the AGORA High-resolution Galaxy Simulations Comparison Project (arXiv:1308.2669,arXiv:1610.03066). These simulations show the formation and evolution of a Milky Way-sized galaxy using eight of the most widely used numerical codes in the community (Art-I, Enzo, Ramses, Changa, Gadget-3, Gear, Gizmo, and Arepo). In this short report, we describe the public release of the raw output data from all of these simulations at z = 8, 7, 6, 5, 4, 3, 2 (plus at z=1, 0 when available), and several metadata files containing the halo centers, virial quantities, and merger trees. The data from even thinner timesteps will be released as soon as the upcoming collaboration papers (VII-IX) are submitted and accepted.

Evidently, low-luminosity active galactic nuclei LLAGNs comprise of inner advective disk and outer geometrically thin disk. The wind is inevitable in LLAGNs, mainly interpreted indirect way, also the evidences are growing for the presence of wind in outer thin disk. We present a hydrodynamics HD model for wind from outer thin disk where the main driver is the inner disk irradiation (which is parameterized by a number $x$ in hydrostatic equilibrium equation with a slightly unbalancing role). The model works for a low intense irradiation or from a height $z_s$ in optically thin medium. We solve the model equations in cylindrical coordinate along the $z$-axis for a given radius $r$ with assuming a tiny vertical speed $v_z$ ($\ll c_s$ sound speed). The sonic point conditions assure an isobaric regime above the sonic height ($z^{max}$), in addition from the height $z_f (\ll z^{max}$) the radial pressure gradient also supports the fluid rotation, and both conjointly assure a wind ejection from the $z^{max}$ with fluid speed. The $z^{max}$ increases with $x$, and beyond a large $z^{max}$ (say $z^{max}_t$ corresponded to maximum $x$) there is no physical solution. We start the computation from outer radius $r_o^{thin}$ to inner $r_{in}^{thin}$ with Bondi mass accretion rate $\dot{M}_{Bondi}$, to explore the $r$ dependency of mass inflow rate $\dot{M}$ and wind properties. We constrain the model by fixing $\dot{M}$ at $r_{in}^{thin}$ from the observations of NGC 1097, and check the feasibility of model by comparing the energetics with observed bolometric luminosity. The wind is an equatorial with viewing angle $i>85$degree, and capable to generate red-/blue-shifted lines, it would be a general characteristics for LLAGNs.

The escape of Lyman-alpha photons at redshifts greater than two is an ongoing subject of study and an important quantity to further understanding of Lyman-alpha emitters (LAEs), the transmission of Lyman-alpha photons through the interstellar medium and intergalactic medium, and the impact these LAEs have on cosmic reionisation. This study aims to assess the Lyman-alpha escape fraction over the redshift range 2.9<z<6.7, focusing on VLT/MUSE-selected, gravitationally lensed, intrinsically faint LAEs. These galaxies are of particular interest as the potential drivers of cosmic reionisation. We assessed the Lyman-alpha escape fraction in two ways: through an individual study of 96 LAEs behind the A2744 lensing cluster, with JWST/NIRCam and HST data, and through a study of the global evolution of the escape fraction using the state-of-the-art luminosity functions for LAEs and the UV-selected `parent' population (dust-corrected). We compared these studies to those in the literature based on brighter samples. We find a negligible redshift evolution of the escape fraction for our individual galaxies; it is likely that it was washed out by significant intrinsic scatter. We observed a more significant evolution towards higher escape fractions with decreasing UV magnitude and fit this relation. When comparing the two luminosity functions to derive the escape fraction in a global sense, we saw agreement with previous literature when integrating the luminosity functions to a bright limit. However, when integrating using a faint limit equivalent to the observational limits of our samples, we observed enhanced escape fraction values, particularly around z~6, where the escape fraction becomes consistent with 100%. This indicates for the faint regimes we sampled that galaxies towards reionisation tend to allow very large fractions of Lyman-alpha photons to escape. (shortened)

We investigate the stochastic gravitational wave background (SGWB) from neutrino-dominated accretion flows (NDAFs) based on the results of our fallback core-collapse supernova (CCSN) simulations. We find that the predicted SGWB is mainly determined by the typical CCSN initial explosion energy and progenitor metallicity. For the optimistic cases in which the typical initial explosion energy is low, the SGWB from NDAFs without disk outflows might be detected by next-generation space-based interferometers such as Decihertz Interferometer Gravitational wave Observatory (DECIGO) and Big Bang Observer (BBO). In the low-frequency regime $\sim10^{-3}-10^{-1}$ Hz, this background is comparable to that expected from standard inflationary models. Therefore, the SGWB from NDAFs may become a foreground for searches of the SGWB generated in the inflationary epoch. Combining the diffuse NDAF neutrino background and SGWB from NDAFs, one may constrain the properties of the CCSNe and NDAFs.

Present cosmic microwave background (CMB) observations have significantly advanced our understanding of the universe's origin, especially with primordial gravitational waves (PGWs). Currently, ground-based CMB telescopes are mainly located in the southern hemisphere, leaving an untapped potential for observations in the northern hemisphere. In this work, we investigate the perspective of a northern hemisphere CMB polarization telescope (NHT) to detect PGWs and present mock data for such a project. We forecast the detection sensitivity on the tensor-to-scalar ratio r of NHT and compare it with the existed ground-based experiments, also search for optimal experimental configurations that can achieve the best sensitivity of r. Our results indicate that, considering realistic experimental conditions, the first year of NHT observations combined with Planck can achieve a precision of \sigma (r)= 0.015, reaching the level of BICEP2/Keck, with significant potential for improvement with subsequent instrumentation parameter enhancements.

We explore the multi-scatter capturing of the massive dark matter (DM) particle inside the neutron star via a momentum-dependent dark matter-nucleon scattering cross-section. We find that the capturing enhanced for the positive velocity and momentum transfer dependent DM-nucleon scattering in comparison with the constant cross-section case. Further, a large capture of the DM particles can be thermalized and lead to black hole formation and, therefore, destroy the neutron star. Using the observation of the old neutron star in the DM-dominated region, we obtain strong constraints on massive DM parameters.

We homogeneously measured the elliptical shapes of 163 globular clusters (GCs) using the on-sky distribution of their cluster members and the third data release of the ESA mission Gaia (DR3). The astrometry enables the differentiation of stars within clusters from those in the field. This feature is particularly valuable for clusters located in densely populated areas of the sky, where conventional methods for measuring the geometry of the GCs are not applicable. The median axial ratio of our full sample is $\langle b/a \rangle = 0.935^{+0.033}_{-0.090}$ and $0.986^{+0.009}_{-0.004}$ for the subset of 11 GCs previously studied based on Hubble Space Telescope imaging. We investigated whether the minor axis of the ellipses can be interpreted as a pseudo-rotation axis by comparing it to measurements of cluster rotation. Using the radial velocities from Gaia, we detected rotation for three clusters, NGC 5139, NGC 104, and NGC 6341, and observed an alignment between the pseudo-rotation axis and the 2D projection of the real rotation axis. To expand the set of clusters for which rotation has been detected, we analyzed multiple literature references. Depending on the reference used for comparison, we observed an alignment in between 76% to 100% of the clusters. The lack of an alignment observed in some clusters may be linked to different scales analyzed in various studies. Several studies have demonstrated that the orientation of rotation varies with the distance from the center. We estimate that the next Gaia release will increase the number of stars with radial velocities in GCs from $\sim 10,000$ in Gaia DR3 to $\sim 55,000$ in Gaia DR4. This will enable the measurement of rotation and ellipticities at identical angular scales for additional clusters, which will help us to clarify whether the previously mentioned alignment occurs in all clusters.

We report on the discovery of polarized X-ray emission from an accreting millisecond pulsar. During a 10-day-long coverage of the February 2024 outburst of SRGA J144459.2-604207, the Imaging X-ray Polarimetry Explorer (IXPE) detected an average polarization degree of the 2-8 keV emission of 2.3% +/- 0.4% at an angle of 59° +/- 6° (East of North; uncertainties quoted at the 1${\sigma}$ confidence level). The polarized signal shows a significant energy dependence with a degree of 4.0% +/- 0.5% between 3 and 6 keV and < 2% (90% c.l.) in the 2-3 keV range. We used NICER, XMM-Newton, and NuSTAR observations to obtain an accurate pulse timing solution and perform a phase-resolved polarimetric analysis of IXPE data. We did not detect any significant variability of the Stokes parameters Q and U with the spin and the orbital phases. We used the relativistic rotating vector model to show that a moderately fan-beam emission from two point-like spots at a small magnetic obliquity ($\simeq$ 10°) is compatible with the observed pulse profile and polarization properties. IXPE also detected 52 type-I X-ray bursts, with a recurrence time $\Delta t_{rec}$ increasing from 2 to 8 h as a function of the observed count rate $C$ as as $\Delta t_{rec} \simeq C^{-0.8}$ We stacked the emission observed during all the bursts and obtained an upper limit on the polarization degree of 8.5% (90% c.l.).

We calculate the arrival direction distribution of ultra-high-energy cosmic rays (UHECRs) with a new suite of models of the Galactic magnetic field (GMF), assuming sources follow the large-scale structure of the Universe. We find a significantly reduced dipole amplitude compared to previous GMF models, and trace the change to the accidental position of the peak of the extragalactic UHECR flux, which falls at the boundary of the strong flux de-magnification due to the GMF toward the central region of the Galaxy. This serendipitous sensitivity of UHECR anisotropies to the GMF model will be a powerful probe of the source distribution as well as Galactic and extragalactic magnetic fields. Demagnification by the GMF also impacts visibility of some popular source candidates.

We estimate the loss of nitrogen from Pluto over its lifetime, including the giant planet instability period, which we term the "Wild Years." We analyze the orbital migration of 53 simulated Plutinos, which are Kuiper Belt Objects (KBOs) captured into 3:2 mean-motion resonance with Neptune during the instability. This orbital migration brought the Plutinos from 20 to 30 au to their present-day orbits near 40 au along a nonlinear path that includes orbits with semimajor axes from 10 to 100 au. We model the thermal history that results from this migration and estimate the volatile loss rates due to the ever-changing thermal environment. Due to the early Sun's enhanced ultraviolet radiation, the photochemical destruction rate during the Wild Years was a factor of 100 higher than the present-day rate, but this only results in a loss of ~10 m global equivalent layer (GEL). The enhanced Jeans escape rate varies wildly with time, and a net loss of ~100 cm GEL is predicted. Additionally, we model the impact history during the migration and find that impacts are a net source, not loss, of N2, contributing ~100 cm GEL. The 100 cm GEL is 0.1% of the amount of N2 in Sputnik Planitia. We therefore conclude that Pluto did not lose an excessive amount of volatiles during the Wild Years, and its primordial volatile inventory can be approximated as its present-day inventory. However, significant fractions of this small total loss of N2 occurred during the Wild Years, so estimates made using present-day rates will be underestimates.

AB Dor is a young solar-type star with a surface large-scale magnetic field $10^2$ to $10^3$ times stronger than the that of the Sun. Although strong magnetic fields are thought to inhibit coronal mass ejections (CMEs), dimming signatures typically associated with an eruptive CME were recently observed in AB Dor. The uninterrupted, long-duration dimming signal suggests that a CME took place at a high latitude, where it remained in view as the star rotates. A high-latitude CME is also consistent with observations that indicate that AB Dor hosts polar active regions. To investigate magnetic confinement in AB Dor, we conduct a parametric modelling study of twenty-one CMEs at latitudes $\sim 60^\circ$, varying the location, mass and magnetic field strength of an injected flux rope. Twelve models had the flux rope located in an open magnetic field region, while the remaining nine were in a closed region. Results show that CMEs in open-field regions are in general more likely to erupt. The four eruptive CMEs from closed regions had high free magnetic energies $\gtrsim 3\times 10^{35}$ erg, and ten CMEs, predominantly from the closed-field regions (8/10) were confined. CMEs in closed-field regions exhibited lower kinetic energies, since part of the CME energy was expended to overcome magnetic tension and break open the overlying field. In conclusion our work suggests that eruptive CMEs in AB Dor may occur in high-latitude regions of open magnetic field, as the magnetic tension in such regions does not significantly inhibit the eruption.

Detecting active black holes in dwarf galaxies has proven to be a challenge due to their small size and weak electromagnetic signatures. Mid-infrared variability has emerged as a promising tool that can be used to detect active low-mass black holes in dwarf galaxies. We analyzed 10.4 years of photometry from the ALL$WISE$/NEO$WISE$ multi-epoch catalogs, identifying 25 objects with AGN-like variability. Independent confirmation of AGN activity was found in 68% of these objects using optical and near-infrared diagnostics. Notably, we discovered a near-infrared coronal line [S IX] $\lambda$ 1.252 $\mu$m in J1205, the galaxy with the lowest stellar mass (log M$_{*}$ = 7.5 M$_{\odot}$) and low metallicity (12 + log(O/H) = 7.46) in our sample. Additionally, we found broad Pa$\alpha$ potentially from the BLR in two targets, and their implied black hole masses are consistent with black hole-stellar mass relations. Comparing non-variable galaxies with similar stellar masses and $WISE$ $W1-W2$ colors, we found no clear trends between variability and large-scale galaxy properties. However, we found that AGN activity likely causes redder $W1-W2$ colors in variable targets, while for the non-variable galaxies, the contribution stems from strong star formation activity. A high incidence of optical broad lines was also observed in variable targets. Our results suggest that mid-infrared variability is an effective method for detecting AGN activity in low-mass galaxies and can help uncover a larger sample of active low-mass ($<$ 10$^{6}$ M$_{\odot}$) black holes in the universe.

Direct mid-infrared signatures of silicate clouds in substellar atmospheres were first detected in Spitzer observations of brown dwarfs, although their existence was previously inferred from near-infrared spectra. With the Mid-Infrared Instrument (MIRI) instrument on JWST, we can now more deeply probe silicate features from 8 to 10 microns, exploring specific particle composition, size, and structure. Recent characterization efforts have led to the identification in particular of silica (silicon dioxide, SiO$_2$) cloud features in brown dwarfs and giant exoplanets. Previous modeling, motivated by chemical equilibrium considerations, has primarily focused on magnesium silicates (forsterite, enstatite), crystalline quartz, and amorphous silica to match observations. Here, we explore the previously neglected possibility that other crystalline structures of silica, i.e. polymorphs, may be more likely to form at the pressure and temperature conditions of substellar upper atmospheres. We show how these polymorphs may be distinguished from each other with current JWST observations. We also explore how such particles could form and be dynamically lofted and sedimented throughout the atmosphere, and implications for the underlying chemical and dynamical processes governing these objects. We ultimately propose that accounting for the distinct opacities arising from the possible crystalline structure of cloud materials may act as a powerful, observable diagnostic tracer of atmospheric conditions, where particle crystallinity records the history of the atmospheric regions through which clouds formed and evolved. Finally, we highlight that high fidelity, accurate laboratory measurements of silica polymorphs are critically needed to draw meaningful conclusions about the identities and structures of clouds in substellar atmospheres.

We present the validation of two TESS super-Earth candidates transiting the mid-M dwarfs TOI-6002 and TOI-5713 every 10.90 and 10.44 days, respectively. The first star (TOI-6002) is located $32.038\pm0.019$ pc away, with a radius of $0.2409^{+0.0066}_{-0.0065} R_\odot$, a mass of $0.2105^{+0.0049}_{-0.0048} M_\odot$ and an effective temperature of $3229^{+77}_{-57}$ K. The second star (TOI-5713) is located $40.946\pm0.032$ pc away, with a radius of $0.2985^{+0.0073}_{-0.0072} R_\odot$, a mass of $0.2653\pm0.0061 M_\odot$ and an effective temperature of $3225^{+41}_{-40}$ K. We validated the planets using TESS data, ground-based multi-wavelength photometry from many ground-based facilities, as well as high-resolution AO observations from Keck/NIRC2. TOI-6002 b has a radius of $1.65^{+0.22}_{-0.19} R_\oplus$ and receives $1.77^{+0.16}_{-0.11} S_\oplus$. TOI-5713 b has a radius of $1.77_{-0.11}^{+0.13} R_\oplus$ and receives $2.42\pm{0.11} S_\oplus$. Both planets are located near the radius valley and near the inner edge of the habitable zone of their host stars, which makes them intriguing targets for future studies to understand the formation and evolution of small planets around M-dwarf stars.

Brown dwarfs occupy a middle ground in mass space between gaseous giant planets and ultra-cool dwarf stars, and the characterisation of their orbital orientations may shed light on how these neighbouring objects form. We present an analysis of the Rossiter-McLaughlin (RM) effect across the transit of TOI-2533 $b$, a brown dwarf on a moderately eccentric ($e_b = 0.2476\pm0.0090$) and wide-separation ($a_b/R_\star = 13.34\pm0.30$) orbit around an F8-type star, using data from the NEID/WIYN spectrograph in combination with archival photometry and radial velocity observations. Spin-orbit analyses of brown dwarfs are relatively rare, and TOI-2533 stands out as the fifth brown dwarf system with a measured spin-orbit constraint. We derive a sky-projected stellar obliquity of $\lambda = -7\pm14^{\circ}$ for TOI-2533 $b$, finding that the brown dwarf is consistent with spin-orbit alignment. Our joint model also indicates that TOI-2533 $b$ falls near the lower bound of the hydrogen-burning minimum mass range (M$_b$ = $74.9\pm5.3$ M$_{\rm \tiny Jup}$). Ultimately, we find that TOI-2533 $b$ is consistent with formation from disc fragmentation in a primordially spin-orbit aligned orientation, although we cannot rule out the possibility that the system has been tidally realigned during its lifetime.

In standard cosmology, the late Universe is assumed to be statistically homogeneous and isotropic. However, a recent study based on galaxy clusters by Migkas et al. (2021, arXiv:2103.13904) found an apparent spatial variation of approximately $9\%$ in the Hubble constant, $H_0$, across the sky. The authors utilised galaxy cluster scaling relations between various cosmology-dependent cluster properties and a cosmology-independent property, i.e., the temperature of the intracluster gas $(T)$. A position-dependent systematic bias of $T$ measurements can, in principle, result in an overestimation of apparent $H_0$ variations. In this study, we search for directional $T$ measurement biases by examining the scaling relation between the member galaxy velocity dispersion and the gas temperature $(\sigma_\mathrm{v}-T)$. Additionally, we search for apparent $H_0$ angular variations independently of $T$ by analysing the relations between the X-ray luminosity and Sunyaev-Zeldovich signal with the velocity dispersion, $L_\mathrm{X}-\sigma_\mathrm{v}$ and $Y_\mathrm{SZ}-\sigma_\mathrm{v}$. We utilise Monte Carlo simulations of isotropic cluster samples to quantify the statistical significance of any observed anisotropies. We find no significant directional $T$ measurement biases, and the probability that a directional $T$ bias causes the previously observed $H_0$ anisotropy is only $0.002\%$. On the other hand, from the joint analysis of the $L_\mathrm{X}-\sigma_\mathrm{v}$ and $Y_\mathrm{SZ}-\sigma_\mathrm{v}$ relations, the maximum variation of $H_0$ is found in the direction of $(295^\circ\pm71^\circ, -30^\circ\pm71^\circ)$ with a statistical significance of $3.64\sigma$, fully consistent with arXiv:2103.13904. Our findings strongly corroborate the previously detected spatial anisotropy of galaxy cluster scaling relations using a new independent cluster property, $\sigma_\mathrm{v}$.

Recent studies have revealed that certain nuclear parameters are more dominant than others in governing global neutron star properties, such as its structure or oscillation mode characteristics. Although neutron stars can in general assumed to be cold, in astrophysical scenarios such as newly born neutron stars or remnants of binary neutron star mergers, finite temperature effects play a non-negligible role. In this work, we perform a consistent and systematic investigation of the role of nuclear parameters and thermal effects on neutron star properties and fluid oscillation modes within a full general relativistic scheme. We impose constraints on the parameter space of the relativistic mean field model using state-of-the-art information from terrestrial experiments and multi-messenger astrophysical data. We find effective nucleon mass to be the most important nuclear parameter controlling astrophysical observables of hot neutron stars, similar to the cold beta equilibrated matter. However, we conclude that the interplay among saturation properties and astrophysical observables depends not only on the thermal configurations considered but also on the constraints imposed. We also investigated the role of nuclear saturation parameters on some universal relations for hot NSs which are important in gravitational wave asteroseismology. Our investigation confirmed that these relations are mostly insensitive to nuclear saturation properties and mainly affected by variation of charge fraction in the star.

We argue that the North Polar Spur (NPS) and many less prominent structures are formed by gaseous metal-rich plumes associated with star-forming regions (SFRs). The SFRs located at the tangent to the 3-5 kpc rings might be particularly relevant to NPS. A multi-temperature mixture of gaseous components and cosmic rays rises above the Galactic disk under the action of their initial momentum and buoyancy. Eventually, the plume velocity becomes equal to that of the ambient gas, which rotates with different angular speed than the stars in the disk. As a result, the plumes acquire characteristic bent shapes. An ad hoc model of plumes' trajectories shows an interesting resemblance to the morphology of structures seen in the radio continuum and X-rays.

Binary systems with comparable masses and a surrounding accretion disk can accrete gas through spiral accretion streams penetrating the central cavity formed by tidal interactions. Using three-dimensional Newtonian magnetohydrodynamics simulations, we investigate the possibility of a magnetically arrested accretion flow through the cavity. Rather than solely continuously feeding the binary through spiral accretion streams, the accretion is regulated by the strong magnetic field inside the cavity. Transport of mass and angular momentum onto the binary then proceeds largely periodically in magnetic flux eruption episodes. The ejected flux tubes carry angular momentum outwards and away from the binary, inject hot plasma into the disk and can launch flares. This likely intermittent scenario could have potential implications for the emission signatures of supermassive black hole binaries, and potentially shed light onto the role magnetic fields play in the binary's orbital evolution.