Papers for Wednesday, Jun 05 2024

We investigate the emergence phase of young star clusters in the nearby spiral galaxy NGC 628. We use JWST NIRCam and MIRI observations to create spatially resolved maps of the Pa$\alpha$-1.87 $\mu$m and Br$\alpha$-4.05 $\mu$m hydrogen recombination lines, as well as the 3.3 $\mu$m and 7.7 $\mu$m emission from polycyclic aromatic hydrocarbons (PAHs). We extract 953 compact HII regions and analyze the PAH emission and morphology at ~10 pc scales in the associated photo-dissociation regions (PDRs). While HII regions remain compact, radial profiles help us to define three PAH morphological classes: compact (~ 42%), extended (~ 34%) and open (~24%). The majority of compact and extended PAH morphologies are associated with very young star clusters (<5 Myr), while open PAH morphologies are mainly associated with star clusters older than 3 Myr. We observe a general decrease in the 3.3 $\mu$m and 7.7 $\mu$m PAH band emission as a function of cluster age, while their ratio remains constant with age out to 10 Myr and morphological class. The recovered PAH$_{3.3 \mu{\rm m}}$/PAH$_{7.7 \mu{\rm m}}$ ratio is lower than values reported in the literature for reference models and analyses conducted at larger physical scales. The 3.3 $\mu$m and 7.7 $\mu$m bands are typically associated to neutral and ionised PAHs, respectively. While we expected neutral PAHs to be suppressed in proximity of the ionizing source, the constant PAH$_{3.3 \mu{\rm m}}$/PAH$_{7.7 \mu{\rm m}}$ ratio would indicate that both families of molecules disrupt at similar rates in proximity of the HII regions.

Stars grazing supermassive black holes (SMBHs) on bound orbits may produce periodic flares over many passages, known as repeating partial tidal disruption events (TDEs). Here we present 3D hydrodynamic simulations of sun-like stars over multiple tidal encounters. The star is significantly restructured and becomes less concentrated as a result of mass loss and tidal heating. The vulnerability to mass loss depends sensitively on the stellar density structure, and the strong correlation between the fractional mass loss $\Delta M/M_*$ and the ratio of the central and average density $\rho_{\mathrm{c}}/\bar\rho$, which was initially derived in disruption simulations of main-sequence stars, also applies for stars strongly reshaped by tides. Over multiple orbits, the star loses progressively more mass in each encounter and is doomed to a complete disruption. Throughout its lifetime, the star may produce numerous weak flares (depending on the initial impact parameter), followed by a couple of luminous flares whose brightness increases exponentially. Flux-limited surveys are heavily biased towards the brightest flares, which may appear similar to the flare produced by the same star undergoing a full disruption on its first tidal encounter. This places new challenges on constraining the intrinsic TDE rates, which needs to take repeating TDEs into account. Other types of stars with different initial density structure (e.g., evolved stars with massive cores) follow distinct evolution tracks, which might explain the diversity of the long-term luminosity evolution seen in recently uncovered repeaters.

We report on the final two days of a multiwavelength campaign of Sgr A* observing in the radio, submillimeter, infrared, and X-ray bands in July 2019. Sgr A* was remarkably active, showing multiple flaring events across the electromagnetic spectrum. We detect a transient $\sim35$-minute periodicity feature in Spitzer Space Telescope light curves on 21 July 2019. Time-delayed emission was detected in ALMA light curves, suggesting a hotspot within the accretion flow on a stable orbit. On the same night, we observe a decreased flux in the submillimeter light curve following an X-ray flare detected by the Chandra X-ray Observatory and model the feature with an adiabatically expanding synchrotron hotspot occulting the accretion flow. The event is produced by a plasma $0.55~R_{\text{S}}$ in radius with an electron spectrum $p=2.84$. It is threaded by a $\sim130$ Gauss magnetic field and expands at $0.6\%$ the speed of light. Finally, we reveal an unambiguous flare in the infrared, submillimeter, and radio, demonstrating that the variable emission is intrinsically linked. We jointly fit the radio and submillimeter light curves using an adiabatically expanding synchrotron hotspot and find it is produced by a plasma with an electron spectrum $p=0.59$, $187$ Gauss magnetic field, and radius $0.47~R_{\text{S}}$ that expands at $0.029c$. In both cases, the uncertainty in the appropriate lower and upper electron energy bounds may inflate the derived equipartition field strengths by a factor of 2 or more. Our results confirm that both synchrotron- and adiabatic-cooling processes are involved in the variable emission's evolution at submillimeter and infrared wavelengths.

Baryonification algorithms model the impact of galaxy formation and feedback on the matter field in gravity-only simulations by adopting physically motivated parametric prescriptions. In this paper, we extend these models to describe gas temperature and pressure, allowing for a self-consistent modelling of the thermal Sunyaev-Zeldovich effect, weak gravitational lensing, and their cross-correlation, down to small scales. We validate our approach by showing that it can simultaneously reproduce the electron pressure, gas, stellar, and dark matter power spectra as measured in all BAHAMAS hydrodynamical simulations. Specifically, with only two additional free parameters, we can fit the electron pressure auto- and cross-power spectra at 10% while reproducing the suppression in the matter power spectrum induced by baryons at the per cent level, for different AGN feedback strengths in BAHAMAS. Furthermore, we reproduce BAHAMAS convergence and thermal Sunyaev Zeldovich angular power spectra within 1% and 10% accuracy, respectively, down to l = 5000. When used jointly with cosmological rescaling algorithms, the baryonification presented here will allow a fast and accurate exploration of cosmological and astrophysical scenarios. Therefore, it can be employed to create mock catalogues, lightcones, and large training sets for emulators aimed at interpreting forthcoming multi-wavelength observations of the large-scale structure of the Universe.

Galaxy mergers play a critical role in galaxy evolution - altering the size, morphology, dynamics and composition of galaxies. So far, galaxy mergers have mostly been identified through visual inspection of their rest-frame optical and NIR emission. But, dust can obscure this emission, resulting in the misclassification of mergers as single galaxies, and the incorrect interpretation of their baryonic properties. Having serendipitously discovered a dust-obscured galaxy merger at z = 1.17, we aim to determine the baryonic properties of the two merging galaxies, including the star formation rate, and stellar, molecular gas, and dust masses. Using Band 3 and 6 observations from the Atacama Large Millimeter and submillimeter Array (ALMA), and ancillary data, we study the morphology of this previously misclassified merger. We deblend the emission, derive the gas masses from CO observations, and model the spectral energy distributions, to determine the properties of each galaxy. Using the rare combination of ALMA CO(2-1), CO(5-4) and dust-continuum (rest-frame 520um) observations, we provide insights into the gas and dust content and ISM properties of each merger component. We find that only one of the two galaxies is highly dust-obscured, whereas both are massive (> 10^10.5 Msun), highly star-forming (SFR = 60-900Msun/yr), have a moderate-to-low depletion time (tdepl < 0.7Gyr) and high gas fraction ( fgas >= 1). These properties can be interpreted as the positive impact of the merger. With this serendipitous discovery, we highlight the power of (sub)millimeter observations to identify and characterise the individual components of obscured galaxy mergers.

We report long baseline interferometric observations with the CHARA Array that resolve six previously known double-lined spectroscopic binary systems in the Hyades cluster, with orbital periods ranging from 3 to 358 days: HD 27483, HD 283882, HD 26874, HD 27149, HD 30676, and HD 28545. We combine those observations with new and existing radial-velocity measurements, to infer the dynamical masses for the components as well as the orbital parallaxes. For most stars the masses are determined to better than 1%. Our work significantly increases the number of systems with mass determinations in the cluster. We find that while current models of stellar evolution for the age and metallicity of the Hyades are able to reproduce the overall shape of the empirical mass-luminosity relation, they overestimate the $V$-band fluxes by about 0.1 mag between 0.5 and 1.4 $M_{\odot}$. The disagreement is smaller in $H$, and near zero in $K$, and depends somewhat on the model. We also make use of the TESS light curves to estimate rotation periods for our targets, and detect numerous flares in one of them (HD 283882), estimating an average flaring rate of 0.44 events per day.

The infall of the Large Magellanic Cloud (LMC) is predicted to displace the inner Milky Way (MW), imprinting an apparent 'reflex motion' on the observed velocities of distant halo stars. We construct the largest all-sky spectroscopic dataset of luminous red giant stars from $50-160$ kpc, including a new survey of the southern celestial hemisphere. We fit the full 6D kinematics of our data to measure the amplitude and direction of the inner MW's motion towards the outer halo. The observed velocity grows with distance such that, relative to halo stars at $100$ kpc, the inner MW is lurching at $\approx 40$ km s$^{-1}$ towards a recent location along the LMC's past orbit. Our measurements align with N-body simulations of the halo's response to a $1.8 \times 10^{11} M_\odot$ LMC on first infall, suggesting that the LMC is at least 15% as massive as the MW. Our findings highlight the dramatic disequilibrium of the MW outskirts, and will enable more accurate measurements of the total mass of our Galaxy.

We present new deep, wide-field Large Binocular Telescope (LBT) $g$ and $r$ imaging data from the Smallest Scale of Hierarchy Survey (SSH) revealing previously undetected tidal features and stellar streams in the outskirts of six dwarf irregular galaxies (NGC 5238, UGC 6456, UGC 6541, UGC 7605, UGC 8638, and UGC 8760) with stellar masses in the range $1.2 \times 10^7$ M$_{\odot}$ to $1.4 \times 10^8$ M$_{\odot}$. The six dwarfs are located 1-2 Mpc away from large galaxies, implying that the observed distortions are unlikely to be due to tidal effects from a nearby, massive companion. At the dwarfs' distances of $\sim$3-4 Mpc, the identified tidal features are all resolved into individual stars in the LBT images and appear to be made of a population older than 1-2 Gyr, excluding the possibility that they result from irregular and asymmetric star formation episodes that are common in gas-rich dwarf galaxies. The most plausible explanation is that we are witnessing the hierarchical merging assembling of these dwarfs with their satellite populations, a scenario also supported by the peculiar morphology and disturbed velocity field of their HI component. From the SSH sample we estimate a fraction of late type dwarfs showing signs of merging with satellites of $\sim$13\%, in agreement with other recent independent studies and theoretical predictions within the $\Lambda$CDM cosmological framework.

$\omega$ Centauri, the most massive globular cluster in the Milky Way, has long been suspected to be the stripped nucleus of a dwarf galaxy that fell into the Galaxy a long time ago. There is considerable evidence for this scenario including a large spread in metallicity and an unusually large number of distinct sub-populations seen in photometric studies. In this work, we use new MUSE spectroscopic and HST photometric catalogs to investigate the underlying metallicity distributions as well as the spatial variations of the populations within the cluster up to its half-light radius. Based on 11,050 member stars, the [M/H] distribution has a median of $ (-1.614 \pm 0.003)$ dex and a large spread of $\sim$ 1.37 dex reaching from $ -0.67$ dex to $ -2.04$ dex for 99.7 % of the stars. In addition, we show the chromosome map of the cluster, which separates the red giant branch stars into different sub-populations, and analyze the sub-populations of the metal-poorest component. Finally, we do not find any metallicity gradient within the half-light radius, and the different sub-populations are well mixed.

We report the X-ray polarization properties of the high-synchrotron-peaked (HSP) blazar PKS 2155$-$304 based on observations with the Imaging X-ray Polarimetry Explorer (IXPE). We observed the source between Oct 27 and Nov 7, 2023. We also conducted an extensive contemporaneous multiwavelength (MW) campaign. We find that during the first half ($T_1$) of the IXPE pointing, the source exhibited the highest X-ray polarization degree detected for an HSP blazar thus far, (30.7$\pm$2.0)%, which dropped to (15.3$\pm$2.1)% during the second half ($T_2$). The X-ray polarization angle remained stable during the IXPE pointing at 129.4$^\circ$$\pm$1.8$^\circ$ and 125.4$^\circ$$\pm$3.9$^\circ$ during $T_1$ and $T_2$, respectively. Meanwhile, the optical polarization degree remained stable during the IXPE pointing, with average host-galaxy-corrected values of (4.3$\pm$0.7)% and (3.8$\pm$0.9)% during the $T_1$ and $T_2$, respectively. During the IXPE pointing, the optical polarization angle changed achromatically from $\sim$140$^\circ$ to $\sim$90$^\circ$ and back to $\sim$130$^\circ$. Despite several attempts, we only detected (99.7% conf.) the radio polarization once (during $T_2$, at 225.5 GHz): with degree (1.7$\pm$0.4)% and angle 112.5$^\circ$$\pm$5.5$^\circ$. The direction of the broad pc-scale jet is rather ambiguous and has been found to point to the east and south at different epochs; however, on larger scales (> 1.5 pc) the jet points toward the southeast ($\sim$135$^\circ$), similar to all of the MW polarization angles. Moreover, the X-ray to optical polarization degree ratios of $\sim$7 and $\sim$4 during $T_1$ and $T_2$, respectively, are similar to previous IXPE results for several HSP blazars. These findings, combined with the lack of correlation of temporal variability between the MW polarization properties, agree with an energy-stratified shock-acceleration scenario in HSP blazars.

We show with Gaia XP spectroscopy that extremely metal-rich stars in the Milky Way (EMR; $[M/H]_{XP} > 0.5$) - but only those - are largely confined to a tight "knot" at the center of the Galaxy. This EMR knot is round in projection, has a fairly abrupt edge near $\sim 1.5$kpc, and is a dynamically hot system. This central knot also contains very metal-rich (VMR; $+0.2\le [M/H]_{XP} \le +0.4$) stars. However, in contrast to EMR stars, the bulk of VMR stars form an extended, highly flattened distribution in the inner Galaxy ($R_{\mathrm{GC}}\lesssim 5$ kpc). We draw on TNG50 simulations of Milky Way analogs for context and find that compact, metal-rich knots confined to $<1.5$kpc are a universal feature. In typical simulated analogs, the top 5-10% most metal-rich stars are confined to a central knot; however, in our Milky Way data this fraction is only 0.1%. Dust-penetrating wide-area near-infrared spectroscopy, such as SDSS-V, will be needed for a rigorous estimate of the fraction of stars in the Galactic EMR knot. Why in our Milky Way only EMR giants are confined to such a central knot remains to be explained. Remarkably, the central few kiloparsecs of the Milky Way harbor both the highest concentration of metal-poor stars (the `poor old heart') and almost all EMR stars. This highlights the stellar population diversity at the bottom of galactic potential wells.

Although supermassive black holes (SMBHs) reside in the heart of virtually every massive galaxy, it remains debated whether dwarf galaxies also commonly host SMBHs. Because low-mass galaxies may retain a memory of the assembly history of their black holes, probing the black hole occupation fraction of local dwarf galaxies might offer insights into the growth and seeding mechanisms of the first black holes. In this work, we exploit the Western half of the eROSITA all-sky survey (covering $20,000~\rm{deg^2}$) and compile a catalog of accreting SMBHs in local ($D<200$~Mpc) dwarf galaxies. After cleaning our sample from cosmic X-ray background sources, X-ray binaries, and ultraluminous X-ray sources, we identify 74 AGN-dwarf galaxy pairs. Using this large and uniform sample, we derive a luminosity function of dwarf galaxy AGN, fitting it with a power law function and obtaining ${\rm d}N/{\rm d}L_{\rm X} = (15.9\pm2.2)\times L_{\rm X}^{-1.63\pm0.05}$. Measuring the offset between the centroid of dwarf galaxies and the X-ray sources, we find that about $50\%$ of the AGN are likely off-nuclear, in agreement with theoretical predictions. We also compare the black hole-to-stellar mass relation of the AGN in our sample with the local and high-redshift relations, finding that our sources better adhere to the former. This suggests that local AGN across different mass scales underwent a similar growth history. Finally, we compare our sources with semi-analytical models: while our sample is too shallow to distinguish between different seeding models, it favors a growth mechanism linked to the star-formation rate of the host galaxy.

CHEOPS is a space telescope specifically designed to monitor transiting exoplanets orbiting bright stars. In September 2023, CHEOPS completed its nominal mission and remains in excellent operational conditions. The mission has been extended until the end of 2026. Scientific and instrumental data have been collected throughout in-orbit commissioning and nominal operations, enabling a comprehensive analysis of the mission's performance. In this article, we present the results of this analysis with a twofold goal. First, we aim to inform the scientific community about the present status of the mission and what can be expected as the instrument ages. Secondly, we intend for this publication to serve as a legacy document for future missions, providing insights and lessons learned from the successful operation of CHEOPS. To evaluate the instrument performance in flight, we developed a comprehensive monitoring and characterisation programme. It consists of dedicated observations that allow us to characterise the instrument's response. In addition to the standard collection of nominal science and housekeeping data, these observations provide input for detecting, modelling, and correcting instrument systematics, discovering and addressing anomalies, and comparing the instrument's actual performance with expectations. The precision of the CHEOPS measurements has enabled the mission objectives to be met and exceeded. Careful modelling of the instrumental systematics allows the data quality to be significantly improved during the light curve analysis phase, resulting in more precise scientific measurements. CHEOPS is compliant with the driving scientific requirements of the mission. Although visible, the ageing of the instrument has not affected the mission's performance.

Accurate photometric redshift (photo-$z$) estimation requires support from multi-band observational data. However, in the actual process of astronomical observations and data processing, some sources may have missing observational data in certain bands for various reasons. This could greatly affect the accuracy and reliability of photo-$z$ estimation for these sources, and even render some estimation methods unusable. The same situation may exist for the upcoming Chinese Space Station Telescope (CSST). In this study, we employ a deep learning method called Generative Adversarial Imputation Networks (GAIN) to impute the missing photometric data in CSST, aiming to reduce the impact of data missing on photo-$z$ estimation and improve estimation accuracy. Our results demonstrate that using the GAIN technique can effectively fill in the missing photometric data in CSST. Particularly, when the data missing rate is below 30\%, the imputation of photometric data exhibits high accuracy, with higher accuracy in the $g$, $r$, $i$, $z$, and $y$ bands compared to the $NUV$ and $u$ bands. After filling in the missing values, the quality of photo-$z$ estimation obtained by the widely used Easy and Accurate Zphot from Yale (EAZY) software is notably enhanced. Evaluation metrics for assessing the quality of photo-$z$ estimation, including the catastrophic outlier fraction ($f_{out}$), the normalized median absolute deviation ($\rm {\sigma_{NMAD}}$), and the bias of photometric redshift ($bias$), all show some degree of improvement. Our research will help maximize the utilization of observational data and provide a new method for handling sample missing values for applications that require complete photometry data to produce results.

Palomar Gattini-IR is a wide-field, synoptic infrared time domain survey covering $\approx 15000$ sq. deg. of the sky at $\approx 1-3$ night cadence to a depth of $J\approx 13.0$ and $\approx 14.9$ Vega mag in and outside the Galactic plane, respectively. Here, we present the first data release of $J$-band light curves of 2MASS sources within the survey footprint covering approximately the first four years of operations. We describe the construction of the source catalog based on the 2MASS point sources, followed by exposure filtering criteria and forced PSF photometry. The catalog contains light curves of $\approx 286$ million unique sources with 2MASS magnitudes of $J < 15.5$ mag, with a total of $\approx 50$ billion photometric measurements and $\approx 20$ billion individual source detections at signal-to-noise-ratio $> 3$. We demonstrate the photometric fidelity of the catalog by i) quantifying the magnitude-dependent accuracy and uncertainty of the photometry with respect to 2MASS and ii) comparing against forced PGIR aperture photometry for known variable sources. We present simple filtering criteria for selecting reliable photometric measurements as well as example Python notebooks for users. This catalog is the largest compilation of nightly cadence, synoptic infrared light curves to date, comparable to those in the largest optical surveys, providing a stepping stone to upcoming infrared surveys in the coming decade.

We examined archival Far Ultraviolet Spectroscopic Explorer data to search for far-ultraviolet emission lines in the starburst galaxy M82. The observations were made in an outflow region that extends beyond the galactic disk. We found the O VI $\lambda\lambda$ 1032, 1038 emission lines from the galaxy's southern outflow region. The O VI lines suggest that the outflowing warm-hot gas is undergoing radiative cooling. We measured a radial velocity of $\sim$420 km s$^{-1}$ from the O VI lines, which is faster than the velocity seen in H$\alpha$ observations. The O VI $\lambda$1038 emission line seems to be blended with the C II $\lambda$1037 emission line, which has a radial velocity of $\sim$300 km s$^{-1}$, similar to what is observed in H$\alpha$ observations. The outflow medium of M82 appears to be composed of gas in multiple phases with varying temperatures and kinematics. Future spectroscopic observations in high energy regimes covering a wider spatial area are necessary to understand better the properties of the warm-hot gas medium in the outflow.

Scaling relations of galaxy clusters are a powerful probe of cosmic isotropy in the late Universe. Owing to their strong cosmological dependence, galaxy cluster scaling relations can obtain tight constraints on the spatial variation of cosmological parameters, such as the Hubble constant ($H_0$), and detect large-scale bulk flow motions. Such tests are crucial to scrutinise the validity of $\Lambda$CDM in the local Universe and determine at what cosmic scales (if any) extra-galactic objects converge to isotropy within the Cosmic Microwave Background rest frame. This review describes the methodology for conducting cosmic isotropy tests with cluster scaling relations and examines possible systematic biases. We also discuss the results of past studies that reported statistically significant observed anisotropies in the local Universe. Finally, we explore the future potential of cluster scaling relations as a cosmic isotropy probe given the large amount of multi-wavelength cluster data expected in the near future.

We present an analysis of gravitational wave (GW) predictions from five two-dimensional Core Collapse Supernova (CCSN) simulations that varied only in the Equation of State (EOS) implemented. The GW signals from these simulations are used to produce spectrograms in the absence of noise, and the emergent high-frequency feature (HFF) is found to differ quantitatively between simulations. Below 1 kHz, the HFF is well approximated by a first-order polynomial in time. The resulting slope was found to vary between 10-50% across all models. Further, using real interferometric noise we investigated the current capabilities of GW detectors to resolve these differences in HFF slope for a Galactic CCSN. We find that for distances up to 1 kpc, current detectors can resolve HFF slopes that vary by at least 30%. For further Galactic distances, current detectors are capable of distinguishing the upper and lower bounds of the HFF slope for groupings of our models that varied in EOS. With the higher sensitivity of future GW detectors, and with improved analysis of the HFF, our ability to resolve properties of the HFF will improve for all Galactic distances. This study shows the potential of using the HFF of CCSN produced GWs to provide insight into the physical processes occurring deep within CCSN during collapse, and in particular its potential to further constrain the EOS through GW detection.

Context. The carbamoyl radical (H2NCO) is believed to play a central role in the ice-grain chemistry of crucial interstellar complex organic molecules as formamide and acetamide. Yet, little is known about this radical that remains elusive in laboratory gas-phase experiments. Aims. In order to enable interstellar searches of H2NCO, we have undertaken a mandatory laboratory characterisation of its pure rotational spectrum. Methods. We report the gas-phase laboratory detection of H2NCO, produced by H-atom abstraction from formamide, using pure rotational spectroscopy at millimetre and submillimetre wavelengths. Millimetre-wave data were acquired using chirped-pulse Fourier-transform spectroscopy while submillimetre-wave ones were obtained using Zeeman-modulated spectroscopy. Experimental measurements were guided by quantum-chemical calculations at the $\omega$B97X-D/cc-pVQZ level of theory. Interstellar searches for the radical have been undertaken on the Protostellar Interferometric Line Survey (PILS) towards the solar-type protostar IRAS 16293-2422. Results. From the assignment and fit of experimental transitions up to 660 GHz, reliable spectroscopic parameters for H2NCO in its ground vibrational state have been derived, enabling accurate spectral predictions. No transitions of the radical were detected on the PILS survey. The inferred upper limit shows that H2NCO abundance is at least 60 times below that of formamide and 160 times below that of HNCO in this source; a value that is in agreement with predictions from a physico-chemical model of this young protostar.

We measure the clustering of Lyman Alpha Emitting galaxies (LAEs) selected from the One-hundred-square-degree DECam Imaging in Narrowbands (ODIN) survey, with spectroscopic follow-up from Dark Energy Spectroscopic Instrument (DESI). We use DESI spectroscopy to optimize our selection and to constrain the interloper fraction and redshift distribution of our narrow-band selected sources. We select samples at z=2.45 and 3.1 in the COSMOS field with median Ly-alpha fluxes of 10^{-16}erg/s/cm2. Covariances and cosmological inferences are obtained from a series of mock catalogs built upon high-resolution N-body simulations that match the footprint, number density, redshift distribution and observed clustering of the sample. We find that both samples have a correlation length of r_0=3.0+/-0.2 Mpc/h. Within our fiducial cosmology these correspond to 3D number densities of 10^{-3}h3/Mpc3 and, from our mock catalogs, biases of 1.7 and 2.0 at z=2.45 and 3.1, respectively. We discuss the implications of these measurements for the use of LAEs as large-scale structure tracers for high-redshift cosmology.

We investigate the evolution of magnetic fields generated by the crystallization-driven dynamo in carbon-oxygen white dwarfs (WDs) with masses $\lesssim1.05\ M_{\odot}$. We use scalings for the dynamo to demonstrate that the initial magnetic field strength ($B_{0}$) has an upper limit that depends on the initial convection zone size ($R_{\mathrm{out},0}$) and the WD mass. We solve the induction equation to follow the magnetic field evolution after the dynamo phase ends. We show that the predicted surface magnetic field strength ($B_{\mathrm{surf}}$) differs from $B_{0}$ by at least a factor of $\sim$0.3. This reduction depends on $R_{\mathrm{out},0}$, where values smaller than half of the star radius give $B_{\mathrm{surf}}\lesssim0.01\ B_{0}$. We implement electrical conductivities that account for the solid phase effect on the Ohmic diffusion. We observe that the conductivity increases as the solid core grows, freezing in the magnetic field at a certain point of the evolution and slowing its outwards transport. We study the effect of turbulent magnetic diffusivity induced by the convection and find that for a small $R_{\mathrm{out},0}$, $B_{\mathrm{surf}}$ is stronger than the non-turbulent diffusion cases because of the more rapid transport, but still orders of magnitude smaller than $B_{0}$. Given these limitations, the crystallization-driven dynamo theory could explain only magnetic C/O WDs with field strengths less than a few MG for the mass range 0.45-1.05 $M_{\odot}$. Our results also suggest that a buried fossil field must be at least 100 times stronger than observed surface fields if crystallization-driven convection is responsible for its transport to the surface.

The main challenge of exoplanet high-contrast imaging (HCI) is to separate the signal of exoplanets from their host stars, which are many orders of magnitude brighter. HCI for ground-based observations is further exacerbated by speckle noise originating from perturbations in the Earth's atmosphere and imperfections in the telescope optics. Various data post-processing techniques are used to remove this speckle noise and reveal the faint planet signal. Often, however, a significant part of the planet signal is accidentally subtracted together with the noise. In the present work, we use explainable machine learning to investigate the reason for the loss of the planet signal for one of the most used post-processing methods: Principal Component Analysis (PCA). We find that PCA learns the shape of the telescope point spread function for high numbers of PCA components. This representation of the noise captures not only the speckle noise, but also the characteristic shape of the planet signal. Building upon these insights, we develop a new post-processing method (4S) that constrains the noise model to minimize this signal loss. We apply our model to 11 archival HCI datasets from the VLT-NACO instrument in the L'-band and find that our model consistently outperforms PCA. The improvement is largest at close separations to the star ($\leq 4 \lambda /D$) providing up to 1.5 magnitudes deeper contrast. This enhancement enables us to detect the exoplanet AF Lep b in data from 2011, 11 years before its subsequent discovery. We present updated orbital parameters for this object.

Observations of transient phenomena, such as GRBs, FRBs, novae/supernovae explosions, coupled with the detection of cosmic messengers like high-energy neutrinos and gravitational waves, have transformed astrophysics. Maximizing the discovery potential necessitates tools for swiftly acquiring an overview of the most relevant information for each new detection. Introducing Astro-COLIBRI, a comprehensive platform designed to meet this challenge. Astro-COLIBRI features a public API, real-time databases and alert systems, a discussion forum, and a website and iOS/Android apps as user clients. In real time, it evaluates incoming astronomical observation messages from all available alert streams, filters them based on user-defined criteria, and contextualizes them in the multi-wavelength (MWL) and multi-messenger (MM) context. User clients offer a graphical representation, providing a succinct summary for quick identification of interesting phenomena and assessing observing conditions globally.

Many accretion-powered pulsars rotate in magnetocentrifugal disequilibrium, spinning up or down secularly over multi-year intervals. The magnetic dipole moment $\mu$ of such systems cannot be inferred uniquely from the time-averaged aperiodic X-ray flux $\langle L(t) \rangle$ and pulse period $\langle P(t) \rangle$, because the radiative efficiency of the accretion is unknown and degenerate with the mass accretion rate. Here we circumvent the degeneracy by tracking the fluctuations in the unaveraged time series $L(t)$ and $P(t)$ using an unscented Kalman filter, whereupon $\mu$ can be estimated uniquely, up to the uncertainties in the mass, radius and distance of the star. The analysis is performed on Rossi X-ray Timing Explorer observations for $24$ X-ray transients in the Small Magellanic Cloud, which have been monitored regularly for $\sim 16$ years. As well as independent estimates of $\mu$, the analysis yields time-resolved histories of the mass accretion rate and the Maxwell stress at the disk-magnetosphere boundary for each star, and hence auto- and cross-correlations involving the latter two state variables. The inferred fluctuation statistics convey important information about the complex accretion physics at the disk-magnetosphere boundary.

New JWST near-infrared imaging of the nearby galaxy NGC 628 from the Cycle 1 program JWST-FEAST is combined with archival JWST mid-infrared imaging to calibrate the 21 $\mu$m emission as a star formation rate indicator (SFR) at $\sim$120 pc scales. The Pa$\alpha$ ($\lambda$1.8756 $\mu$m) hydrogen recombination emission line targeted by FEAST provides a reference SFR indicator that is relatively insensitive to dust attenuation, as demonstrated by combining this tracer with the HST H$\alpha$ imaging. Our analysis is restricted to regions that appear compact in nebular line emission and are sufficiently bright to mitigate effects of both age and stochastic sampling of the stellar initial mass function. We find that the 21 $\mu$m emission closely correlates with the nebular line emission, with a power-law with exponent=1.07$\pm$0.01, in agreement with past results. We calibrate a hybrid SFR indicator using a combination of H$\alpha$ and 24 $\mu$m (extrapolated from 21 $\mu$m) tracers and derive the proportionality constant between the two tracers $b=0.095\pm0.007$, which is $\sim$ 3-5 times larger than previous derivations using large regions/entire galaxies. We model these discrepancies as an increasing contribution to the dust heating by progressively older stellar populations for increasing spatial scales, in agreement with earlier findings that star formation is hierarchically distributed in galaxies. Thus, use of hybrid SFR indicators requires prior knowledge of the mean age of the stellar populations dominating the dust heating, which makes their application uncertain. Conversely, non-linear calibrations of SFRs from L(24) alone are more robust, with a factor $\lesssim$2.5 variation across the entire range of L(24) luminosities from HII regions to galaxies.

The Simons Observatory is a new ground-based cosmic microwave background experiment, which is currently being commissioned in Chile's Atacama Desert. During its survey, the observatory's small aperture telescopes will map 10% of the sky in bands centered at frequencies ranging from 27 to 280 GHz to constrain cosmic inflation models, and its large aperture telescope will map 40% of the sky in the same bands to constrain cosmological parameters and use weak lensing to study large-scale structure. To achieve these science goals, the Simons Observatory is deploying these telescopes' receivers with 60,000 state-of-the-art superconducting transition-edge sensor bolometers for its first five year survey. Reading out this unprecedented number of cryogenic sensors, however, required the development of a novel readout system. The SMuRF electronics were developed to enable high-density readout of superconducting sensors using cryogenic microwave SQUID multiplexing technology. The commissioning of the SMuRF systems at the Simons Observatory is the largest deployment to date of microwave multiplexing technology for transition-edge sensors. In this paper, we show that a significant fraction of the systems deployed so far to the Simons Observatory's large aperture telescope meet baseline specifications for detector yield and readout noise in this early phase of commissioning.

We outline the development of the readout software for the Prime-Cam and Mod-Cam instruments on the CCAT Fred Young Submillimeter Telescope (FYST), primecam_readout. The instruments feature lumped-element kinetic inductance detector (LEKID) arrays driven by Xilinx ZCU111 RFSoC boards. In the current configuration, each board can drive up to 4000 KIDs, and Prime-Cam is implementing approximately 25 boards. The software runs on a centralized control computer connected to the boards via dedicated ethernet, and facilitates such tasks as frequency-multiplexed tone comb driving, comb calibration and optimization, and detector timestream establishment. The control computer utilizes dynamically generated control channels for each board, allowing for simultaneous parallel control over all, while uniquely tracking diagnostics for each. This work demonstrates a scalable RFSoC readout architecture where computational demands increase linearly with the number of detectors, enabling control of tens-of-thousands of KIDs with modest hardware, and opening the door to the next generation of KID arrays housing millions of detectors.

We examine the sky distribution of radio galaxies in the NRAO VLA Sky Survey (NVSS) and the Rapid ASKAP Continuum Survey (RACS). Analyses of these samples have reported tension between their inferred dipoles and the kinematic dipole of the Cosmic Microwave Background (CMB). This represents a challenge to the traditional assumption that the Universe is homogeneous and isotropic on large scales: the cosmological principle. We find that NVSS and RACS contain local radio sources which give a non-negligible contribution to the overall dipole signal. These need to be adequately accounted for since the aim is to probe the composition of the Universe at large scales. By appropriately considering these sources, the inferred dipole amplitude in either sample is reduced. Nonetheless, we find support for a dipole aligning with that of the CMB but larger in amplitude, especially in the joint analysis. However, the 'clustering dipole' - the contribution of local sources to the net inferred dipole - appears to align with the direction of the CMB dipole, and its magnitude increases as deeper nearby sources are considered up to a comoving distance of $\approx 130$ Mpc ($h=0.7$). The significance of this observation in the context of the cosmological principle is unclear and prompts further inquiry.

We present an overview of the LIGHTS (LBT Imaging of Galactic Halos and Tidal Structures) survey, which currently includes 25 nearby galaxies that are on average $\sim$ 1 mag fainter than the Milky Way, and a catalog of 54 low central surface brightness (24 $< \mu_{0,g}$/mag arcsec$^{-2} < 28$) satellite galaxy candidates, most of which were previously uncatalogued. The depth of the imaging exceeds the full 10-year depth of the Rubin Observatory's Legacy Survey of Space and Time (LSST). We find, after applying completeness corrections, rising numbers of candidate satellites as we approach the limiting luminosity (M$_r \sim -8$ mag) and central surface brightness ($\mu_{0,g} \sim 28$ mag arcsec$^{-2}$). Over the parameter range we explore, each host galaxy (excluding those that are in overdense regions, apparently groups) has nearly 4 such candidate satellites to a projected radius of $\sim$ 100 kpc. These objects are mostly just at or beyond the reach of spectroscopy unless they are H I rich or have ongoing star formation. We identify 3, possibly 4, ultra-diffuse satellite galaxies (UDGs; effective radius $ > 1.5$ kpc). This incidence rate falls within expectations of the extrapolation of the published relationship between the number of UDGs and host halo mass. Lastly, we visually identify 12 candidate satellites that host a nuclear star cluster (NSC). The NSC occupation fraction for the sample (12/54) matches that published for satellites of early-type galaxies, suggesting that the parent's morphological type plays at most a limited role in determining the NSC occupation fraction.

Understanding the envelope composition of sub-Neptune-type exoplanets is challenging due to the inherent degeneracy in their interior composition scenarios. Particularly, the H2O/H2 ratio, or can be expressed as the O/H ratio, in the planetary envelope provides crucial insights into the origin of these exoplanets relative to the ice line during formation. Using self-consistent radiative transfer modeling and a rate-based automatic chemical network generator combined with 1D photochemical kinetic-transport atmospheric modeling, we investigate atmospheres of temperate sub-Neptunes, ranging from H2-dominated to H2O-dominated scenarios with Teq = 250-400 K, using K2-18 b (Teq = 255 K), LP 791-18 c (Teq = 324 K), and TOI-270 d (Teq = 354 K) as examples. Our models indicate that using the atmospheric CO2/CH4 ratio to infer the deep-interior H2O/H2 ratio. Applying to recent JWST observations, our findings suggest K2-18 b likely has an interior highly enriched in water (approximately 50% H2O), exceeding the amount of water in a 100x solar metallicity scenario and suggesting a formation history that involved substantial accretion of ices. In contrast, TOI-270 d has an interior composition of approximately 25% H2O, which is comparable to the conventional metallicity framework with a metallicity higher than 100x solar metallicity. Furthermore, our models identify carbonyl sulfide (OCS) and sulfur dioxide (SO2) as strong indicators of at least a 10% water-rich envelope in temperate sub-Neptunes. These results provide a method to delineate the internal composition and formation mechanisms of sub-Neptunes with Teq< ~500 K via atmospheric characterization through transmission spectroscopy.

T CrB is a symbiotic star that experiences nova outbursts every $\sim$ 80~yr. The next, long-anticipated nova outburst should occur during the 2024-2026 period. Here, we present results of high-resolution optical spectroscopy of T CrB in the period 2016 - 2023. In these spectra, we measured the equivalent widths of the H$\alpha$, H$\beta$, HeI and HeII emission lines. The maximum equivalent width (EW) was recorded on May 2021, when the EW of $H\alpha$ reached -44.6 Åand H$\beta$ = -21.5 Å. At the other extreme, the minimum of EW($H\alpha$)= -2.9 Åwas recorded in October 2023. After October 2023, the B-band emission brightened, suggesting a re-appearance of the orbital modulation. In addition to the optical data, we study the X-ray behaviour in the same period. We find a strong correlation between $EW(H\alpha)$ and X-ray flux with a correlation coefficient -0.78 and a significance of 2.6$\times 10^{-5}$.

The Dark Energy Spectroscopic Instrumnet (DESI) collaboration recently released the first year data of baryon acoustic oscillations (BAOs). Basing on the five different tracers, the cosmological constraint shows a hint of deviation from the standard $\Lambda$CDM model. In this letter, We combine the DESI BAOs with other cosmic probes to constrain the evolution of Hubble constant as a function of redshift in flat $\Lambda$CDM model. The non-parametric method is used to estimate the value of Hubble constant at different redshift bins. The correlation among different bins are removed by diagonalizing the covariance matrix. The joint data sample demonstrate a decreasing trend of Hubble constant with a significance of $8.6 \sigma$, which can naturally resolve the Hubble tension. It may be due to dynamical dark energy or modified gravity.

Low-mass X-ray binaries (LMXBs) with neutron stars show quite different features which depend on the rate of mass transfer from the donor star. With a high transfer rate the Z sources are in a persistent soft spectral state, with a moderate rate the transient Atoll sources have outburst cycles like the black hole X-ray binaries. The observations document very long outburst recurrence times for quite a number of sources. We follow with our computations the evolution of the accretion disc until the onset of the ionization instability. For sources with a low mass transfer rate the accumulation of matter in the disc is essentially reduced due to the continuous evaporation of matter from the disc to the coronal flow. Different mass transfer rates result in nearly the same amount of matter accumulated for the outburst which means the outburst properties are similar for sources with short and sources with long outburst cycles, contrary to some expectations. Then of systems with long recurrence time less sources will be detected and the total population of LMXBs could be larger than it appears. This would relieve the apparent problem that the observed number of LMXBs as progenitors of millisecond pulsars (MSP) is too small compared to the number of MSP. Concerning the few quasi-persistent sources with year-long soft states we argue that these states are not outbursts, but quasi-stationary hot states as in Z sources.

LS I + 61$^\circ$303 is a high-mass X-ray binary system comprising a massive Be star and a rapidly rotating neutron star. Its spectral energy distribution across multi-wavelengths categorizes it as a $\gamma$-ray binary system. In our analysis of LS I + 61$^\circ$303 using Fermi-LAT observations, we not only confirmed the three previously discussed periodicities of orbital, superorbital, and orbital-superorbital beat periods observed in multi-wavelength observations, but also identified an additional periodic signal. This newly discovered signal exhibits a period of $\sim$26.3 day at a $\sim7\sigma$ confidence level. Moreover, the power spectrum peak of the new signal gradually decreases as the energy increases across the energy ranges of 0.1-0.3, 0.3-1.0, and 1.0-500.0 GeV. Interestingly, a potential signal with a similar period was found in data obtained from the Owens Valley Radio Observatory 40 m telescope. We suggest that the newly discovered periodic signal may originate from a coupling between the orbital period and the retrograde stellar precession period.

As Type Ia supernova cosmology transitions from a statistics dominated to a systematics dominated era, it is crucial to understand leftover unexplained uncertainties affecting their luminosity, such as the ones stemming from astrophysical biases. Indeed, SNe Ia are standardisable candles, whose absolute magnitude reach a 0.15~mag scatter once empirical correlations with their lightcurve stretch and colour and with their environment are accounted for. In this paper, we investigate how the standardisation process of SNe Ia depends on environment, to ultimately reduce their scatter in magnitude, focusing on colour standardisation. We use the volume-limited ZTF SN Ia DR2 sample, which offers unprecedented statistics for the low redshift ($z<0.06$) range. We first study the colour distribution, focusing on the effects of dust, to then select a dustless subsample of objects from low stellar mass environments and from the outskirts of their host galaxies. We then look at the colour-residuals relation and its associated parameter $\beta$. Finally, we investigate the colour dependency of the environment-dependent magnitude offsets (steps), to try to disentangle intrinsic and extrinsic colour origin. Our sample probes well the red tail of the colour distribution, up to $c=0.8$. The dustless sample exhibits a significantly lower red tail ($4.6\sigma$) in comparison to the whole sample. This suggests that reddening above $c\geq0.2$ is dominated by host interstellar dust absorption. Looking at the colour-residuals relation, we find it to be linear with lightcurve colour. We show hints of a potential evolution of $\beta$ with host stellar mass at a $2.5\sigma$ level. Finally, unlike recent claims from the literature, we see no evolution of steps as a function of lightcurve colour, suggesting that dust may not be the dominating mechanism responsible for the environmental dependency of SNe Ia magnitude.

Type Ia supernova (SN Ia) cosmology relies on the estimation of lightcurve parameters to derive precision distances that leads to the estimation of cosmological parameters. The empirical SALT2 lightcurve modeling that relies on only two parameters, a stretch x1, and a color c, has been used by the community for almost two decades. In this paper we study the ability of the SALT2 model to fit the nearly 3000 cosmology-grade SN Ia lightcurves from the second release of the Zwicky Transient Facility (ZTF) cosmology science working group. While the ZTF data was not used to train SALT2, the algorithm is modeling the ZTF SN Ia optical lightcurves remarkably well, except for lightcurve points prior to -10 d from maximum, where the training critically lacks statistics. We find that the lightcurve fitting is robust against the considered choice of phase-range, but we show the [-10; +40] d range to be optimal in terms of statistics and accuracy. We do not detect any significant features in the lightcurve fit residuals that could be connected to the host environment. Potential systematic population differences related to the SN Ia host properties might thus not be accountable for by the addition of extra lightcurve parameters. However, a small but significant inconsistency between residuals of blue- and red-SN Ia strongly suggests the existence of a phase-dependent color term, with potential implications for the use of SNe Ia in precision cosmology. We thus encourage modellers to explore this avenue and we emphasize the importance that SN Ia cosmology must include a SALT2 retraining to accurately model the lightcurves and avoid biasing the derivation of cosmological parameters.

We present spectropolarimetric observations of an active region recorded simultaneously in the H$\alpha$ Ca II 8662 Å lines. The sunspot exhibits multiple structures, including a lightbridge and a region where Ca II 8662 Å line core is in emission. Correspondingly, the H$\alpha$ line core image displays brightening in the emission region, with the spectral profiles showing elevated line cores. The stratification of the line-of-sight magnetic field is inferred through non-LTE multiline inversions of the Ca II 8662 Å line and the weak field approximation over the H$\alpha$ line. The field strength inferred from the H$\alpha$ line core is consistently smaller than that inferred from inversions at $\log \tau_{500}$ = $-$4.5. However, the study finds no correlation between the WFA over the core of the H$\alpha$ line and that inferred from inversions at $\log \tau_{500}$ = $-$4.5. In regions exhibiting emission features, the morphology of the magnetic field at $\log \tau_{500}$ = $-$4.5 resembles that at $\log \tau_{500}$ = $-$1, with slightly higher or comparable field strengths. The magnetic field morphology inferred from the core of the H$\alpha$ line is also similar to that inferred from the full spectral range of the H$\alpha$ line in the emission region. The field strength inferred in the lightbridge at $\log \tau_{500}$ = $-$1 is smaller than the surrounding umbral regions and comparable at $\log \tau_{500}$ = $-$4.5. Similarly, the field strength inferred in the lightbridge from the WFA over the H$\alpha$ line appears lower compared to the surrounding umbral regions.

The formation of optical fluorescent lines in moving media has not yet been studied in detail, so this work represents a first step in investigating the fluorescence process in different types of macroscopic velocity fields: (a) accelerated outflows, (b) accelerated infalls, and (c) non-monotonic velocity fields (such as an accelerating outflow followed by a deceleration region or an accretion shock front). We solve the radiative transfer equations for the lines involved in the fluorescent process, assuming spherical symmetry and a simplified atomic model. We use the framework of the generalized Sobolev theory for computing the interacting, non-local source functions. The emergent line fluxes are then integrated exactly. Because of Doppler shifts in the moving gaseous envelope, photons of the three lines involved in TTS FeI fluorescence CaII H, FeI 3969, and H_epsilon interact with each other in a complex way, so that fluorescent amplification of the line flux occurs not only for FeI 3969, but also for the other two lines, in all velocity fields that we investigated. With the assumption of LTE populations, the line source functions of moderately optically thick lines are not strongly affected by line interactions, while they are depressed in the inner envelope for optically thick lines because of stellar photon absorption in the interaction regions. Fluorescent amplification takes place mainly in the observer's reference frame during the flux integration. Further comparison with observations will require solving the rate equations for the atomic populations involved, along with the radiation field computed with the method presented here. The main product of this research is the open-source computer code SLIM2 (Spectral Line Interactions in Moving Media), written in Python/Numpy, which numerically solves the fluorescence problem for arbitrary 2D velocities.

The suppression of small-scale matter power spectrum is a distinct feature of Warm Dark Matter (WDM), which permits a constraint on the WDM mass from galaxy surveys. In the thermal relic WDM scenario, quantum statistical effects are not manifest. In a unified framework, we investigate the quantum statistical effects for a fermion case with a degenerate pressure and a boson case with a Bose-Einstein condensation (BEC). Compared to the thermal relic case, the degenerate fermion case only slightly lowers the mass bound while the boson case with a high initial BEC fraction ($\gtrsim90\%$) significantly lowers it. On the other hand, the BEC fraction drops during the relativistic-to-nonrelativistic transition and completely disappears if the initial fraction is below $\sim64\%$. Given the rising interest in resolving the late-time galaxy-scale problems with boson condensation, a question is posed on how a high initial BEC fraction can be dynamically created so that a DM condensed component remains today.

$Context.$ On October 31, 2009, the Photodetector Array Camera and Spectrometer (PACS) on board the Herschel Space Observatory observed far-infrared spectra of Jupiter between 50 and 220$\,\mu$m as part of the program "Water and Related Chemistry in the Solar System". $Aims.$ We investigate the disk-averaged chemical composition of Jupiter's atmosphere as a function of height using these observations. $Methods.$ We used the Planetary Spectrum Generator (PSG) and the least-squares fitting technique to infer the abundances of trace constituents. $Results.$ The PACS data include numerous spectral lines attributable to ammonia (NH$_3$), methane (CH$_4$), phosphine (PH$_3$), water (H$_2$O), and deuterated hydrogen (HD) in the Jovian atmosphere. We infer an ammonia abundance profile that decreases from a mole fraction of $(1.7\pm 0.8)\times 10^{-4}$ at $p\sim 900\,$mbar to $(1.7\pm 0.9)\times 10^{-8}$ at $p\sim 275\,$mbar, following a fractional scale height of about 0.114. For phosphine, we find a mole fraction of $(7.2\pm 1.2)\times 10^{-7}$ at pressures higher than $(550\pm 100)\,$mbar and a decrease of its abundance at lower pressures following a fractional scale height of $(0.09\pm 0.02)$. Our analysis delivers a methane mole fraction of $(1.49\pm 0.09)\times 10^{-3}$. Analyzing the HD $R(0)$ line at $112.1\,\mu$m yields a new measurement of Jupiter's D/H ratio, $\text{D/H}=(1.5\pm 0.6)\times 10^{-5}$. Finally, the PACS data allow us to put the most stringent $3\sigma$ upper limits yet on the mole fractions of hydrogen halides in the Jovian troposphere. These new upper limits are $<1.1\times 10^{-11}$ for hydrogen fluoride (HF), $<6.0\times 10^{-11}$ for hydrogen chloride (HCl), $<2.3\times 10^{-10}$ for hydrogen bromide (HBr) and $<1.2\times 10^{-9}$ for hydrogen iodide (HI) and support the proposed condensation of hydrogen halides into ammonium halide salts in the Jovian troposphere.

We use the `Middle Ages Galaxy Properties with Integral field spectroscopy' (MAGPI) survey to investigate whether galaxies have evolved in the distribution of their stellar angular momentum in the past 3-4 Gyr, as probed by the observational proxy for spin, $\lambda_{R}$. We use 2D stellar kinematics to measure $\lambda_{R}$ along with detailed photometric models to estimate galaxy ellipticity. The combination of these measurements quantifies the kinematic classes of `fast rotators' and the rarer `slow rotators', which show no regular rotation in their line-of-sight velocity fields. We compare 51 MAGPI galaxies with $\log_{10} (M_{\star}/\mathrm{M}_\odot) > 10$ to carefully drawn samples of MaNGA galaxies in the local Universe, selected to represent possible descendants of the MAGPI progenitors. The EAGLE simulations are used to identify possible evolutionary pathways between the two samples, explicitly accounting for progenitor bias in our results and the varied evolutionary pathways a galaxy might take between the two epochs. We find that the occurrence of slow rotating galaxies is unchanged between the MAGPI ($z \sim 0.3$) and MaNGA ($z \sim 0$) samples, suggesting the massive slow rotator population was already in place $\sim 4$ Gyr ago and has not accumulated since. There is a hint of the MAGPI sample having an excess of high $\lambda_{R}$ galaxies compared to the MaNGA sample, corresponding to more ordered rotation, but statistically the samples are not significantly different. The large-scale stellar kinematics, as quantified through the $\lambda_{R}$ parameter, of galaxies at $z \sim 0.3$ have already evolved into the diversity of structures seen today in the local Universe.

We present an optical polarimetric study of a nearby star-forming region, Lambda-Orionis, to map plane-of-the-sky magnetic field geometry to understand the magnetized evolution of the HII region and associated small molecular clouds. We made multi-wavelength polarization observations of 34 bright stars distributed across the region. R-band polarization measurements focused on small molecular clouds BRC 17 and BRC 18 located at the periphery of the HII region are also presented. The magnetic field lines exhibit a large-scale ordered orientation consistent with the Planck sub-mm polarization measurements. The magnetic field lines in both the BRCs are found to be roughly in north-south directions; however, a larger dispersion is noticed in the orientation for BRC 17 compared to BRC 18. Using structure-function analysis, the strength of the plane-of-the-sky component of the magnetic field is estimated as $\sim$28 $\mu$G for BRC 17 and $\sim$40 $\mu$G for BRC 18. The average dust grain size and the mean value of the total-to-selective extinction ratio (R$_{V}$) in the HII region are found to be $\sim$0.51 $\pm$ 0.05 $\mu$m and $\sim$2.9 $\pm$ 0.3, respectively. The distance of the whole HII region is estimated as $\sim$392 $\pm$ 8 pc by combining astrometry information from GAIA EDR3 for YSOs associated with BRCs and confirmed members of central cluster Collinder 69.

In order to examine where, how and why the quenching of star formation begins in the outskirts of galaxy clusters, we investigate the de-projected radial distribution of a large sample of quenched and star-forming galaxies (SFGs) out to $30R_{500}$ around clusters. We identify the SFG sample using radio continuum emission from the Low-Frequency Array Two-metre Sky Survey. We find that the SFG fraction starts to decrease from the field fraction as far out as $10R_{500}$, well outside the virial radius of the clusters. We investigate how the SFG fraction depends on both large-scale and local environments, using radial distance from a cluster to characterise the former, and distance from 5th nearest neighbour for the latter. The fraction of SFGs in high-density local environments is consistently lower than that found in low-density local environments, indicating that galaxies' immediate surroundings have a significant impact on star formation. However, for high-mass galaxies -- and low mass galaxies to a lesser extent -- high-density local environments appear to act as a protective barrier for those SFGs that survived this pre-processing, shielding them from the external quenching mechanisms of the cluster outskirts. For those galaxies that are not in a dense local environment, the global environment causes the fraction of SFGs to decrease toward the cluster centre in a manner that is independent of galaxy mass. Thus, the fraction of SFGs depends on quite a complex interplay between the galaxies' mass, their local environment, and their more global cluster-centric distance.

With the exceptional sensitivity of the Five-hundred-meter Aperture Spherical radio Telescope (FAST), we conducted observations of the neutral hydrogen (HI) in the Circular Galactical Medium (CGM) of Andromeda's (M31) satellite galaxies, specifically Andromeda II, NGC 205, and NGC 185. Initially, three drift scans were executed for these satellites, with a detection limit of $4\times10^{18}$ cm$^{-2}$ ( approximately $1.88\times10^3 M_{\odot}$ of HI mass), followed by a more in-depth scan of a specific region. We discovered a C-shaped HI arc structure sharing a position and line-of-sight velocity similar to a stellar ring structure around Andromeda II, hinting at a potential connection with Andromeda II. In the context of NGC 205, we identified two mass concentrations in the northeast direction, which could be indicative of tidal streams resulting from the interaction between this galaxy and M31. These new lumps discovered could be very helpful in solving the missing interstellar medium (ISM) problem for NGC 205. Observations regarding NGC 185 are consistent with previous studies, and we did not detect any additional HI material around this galaxy. These observational results enhance our understanding of the evolution of these satellite galaxies and provide insight into their historical interactions with the M31 galaxy.

We present a novel approach based on Bayesian field-level inference capable of resolving individual galaxies within the Local Group (LG), enabling detailed studies of its structure and formation via posterior simulations. We extend the Bayesian Origin Reconstruction from Galaxies (BORG) algorithm with a multi-resolution approach, allowing us to reach smaller mass scales and apply observational constraints based on LG galaxies. Our updated data model simultaneously accounts for observations of mass tracers within the dark haloes of the Milky Way (MW) and M31, their observed separation and relative velocity, and the quiet surrounding Hubble flow represented through the positions and velocities of galaxies at distances from one to four Mpc. Our approach delivers representative posterior samples of $\Lambda$CDM realisations that are statistically and simultaneously consistent with all these observations, leading to significantly tighter mass constraints than found if the individual datasets are considered separately. In particular, we estimate the virial masses of the MW and M31 to be $\log_{10}(M_{200c}/M_\odot) = 12.07\pm0.08$ and $12.33\pm0.10$, respectively, their sum to be $\log_{10}(\Sigma M_{200c}/M_\odot)= 12.52\pm0.07$, and the enclosed mass within spheres of radius $R$ to be $\log_{10}(M(R)/M_\odot)= 12.71\pm0.06$ and $12.96\pm0.08$ for $R=1$ Mpc and 3 Mpc, respectively. The M31-MW orbit is nearly radial for most of our $\Lambda$CDM LG's, and most lie in a dark matter sheet that aligns approximately with the Supergalactic Plane, even though the surrounding density field was not used explicitly as a constraint. The approximate simulations employed in our inference are accurately reproduced by high-fidelity structure formation simulations, demonstrating the potential for future high-resolution, full-physics $\Lambda$CDM posterior simulations of LG look-alikes.

General circulation models are a useful tool in understanding the three dimensional structure of hot Jupiter and sub-Neptune atmospheres; however, understanding the validity of the results from these simulations requires an understanding the artificial dissipation required for numerical stability. In this paper, we investigate the impact of the longitudinal filter and vertical ``sponge'' used in the Met Office's {\sc Unified Model} when simulating gaseous exoplanets. We demonstrate that excessive dissipation can result in counter-rotating jets and a catastrophic failure to conserve angular momentum. Once the dissipation is reduced to a level where a super-rotating jet forms, however, the jet and thermal structure are relatively insensitive to the dissipation, except in the nightside gyres where temperatures can vary by $\sim 100\,\mathrm{K}$. We do find, however, that flattening the latitudinal profile of the longitudinal filtering alters the results more than a reduction in the strength of the filtering itself. We also show that even in situations where the temperatures are relatively insensitive to the dissipation, the vertical velocities can still vary with the dissipation, potentially impacting physical processes that depend on the local vertical transport.

Aims. Our aim is to test how different binning strategies previously studied in one-dimensional models perform in three-dimensional radiative hydrodynamic simulations of stellar atmospheres. Methods. Realistic box-in-a-star simulations of the near-surface convection and photosphere of three spectral types (G2V, K0V, and M2V) were run with the MANCHA code with grey opacity. After reaching the stationary state, one snapshot of each of the three stellar simulations was used to compute the radiative energy exchange rate with grey opacity, opacity binned in four $\tau$-bins, and opacity binned in 18 $\{ \tau, \lambda \}$-bins. These rates were compared with the ones computed with opacity distribution functions. Then, stellar simulations were run with grey, four-bin, and 18-bin opacities to see the impact of the opacity setup on the mean stratification of the temperature and its gradient after time evolution. Results. The simulations of main sequence cool stars with the MANCHA code are consistent with those in the literature. For the three stars, the radiative energy exchange rates computed with 18 bins are remarkably close to the ones computed with the opacity distribution functions. The rates computed with four bins are similar to the rates computed with 18 bins, and present a significant improvement with respect to the rates computed with the Rosseland opacity, especially above the stellar surface. The Rosseland mean can reproduce the proper rates in sub-surface layers, but produces large errors for the atmospheric layers of the G2V and K0V stars. In the case of the M2V star, the Rosseland mean fails even in sub-surface layers, owing to the importance of the contribution from molecular lines in the opacity, underestimated by the harmonic mean. Similar conclusions are reached studying the mean stratification of the temperature and its gradient after time evolution.

Chromospheric evaporation (CE) and coronal rain (CR) represent two crucial phenomena encompassing the circulation of mass and energy during solar flares. While CE marks the start of the hot inflow into the flaring loop, CR marks the end, indicating the outflow in the form of cool and dense condensations. With \textit{IRIS} and \textit{AIA/SDO}, we examine and compare the evolution, dynamics, morphology, and energetics of the CR and CE during a C2.1 flare. The CE is directly observed in imaging and spectra in the \ion{Fe}{XXI} line with \textit{IRIS} and in the \ion{Fe}{XVIII} line of AIA, with upward average total speeds of $138\pm[35]~$km~s$^{-1}$ and a temperature of $[9.03\pm3.28]\times10^{6}$~K. An explosive to gentle CE transition is observed, with an apparent reduction in turbulence. From quiescent to gradual flare phase, the amount and density of CR increases by a factor of $\approx4.4$ and 6, respectively. The rain's velocity increases by a 1.4, in agreement with gas pressure drag. In contrast, the clump widths variation is negligible. The location and morphology of CE match closely those of the rain showers, with similar CE sub-structure to the rain strands, reflecting fundamental scales of mass and energy transport. We obtain a CR outflow mass three times larger than the CE inflow mass, suggesting the presence of unresolved CE, perhaps at higher temperatures. The CR energy corresponds to half that of the CE. These results suggest an essential role of coronal rain in the mass-energy cycle of a flare.

The signal from a transiting planet can be diluted by astrophysical contamination. In the case of circumstellar debris disks, this contamination could start in the mid-infrared and vary as a function of wavelength, which would then change the observed transmission spectrum for any planet in the system. The MIRI/LRS WASP-39b transmission spectrum shows an unexplained dip starting at $\sim$10 $\mu$m that could be caused by astrophysical contamination. The spectral energy distribution displays excess flux at similar levels to that which are needed to create the dip in the transmission spectrum. In this article, we show that this dip is consistent with the presence of a bright circumstellar debris disk, at a distance of $>$2 au. We discuss how a circumstellar debris disk like that could affect the atmosphere of WASP-39b. We also show that even faint debris disks can be a source of contamination in MIRI exoplanet spectra.

With the recent announcement by NASA's Planetary Science and Astrobiology Decadal Survey 2023-2032, a priority flagship mission to the planet Uranus is anticipated. Here, we explore the prospects of using the mission's radio Doppler tracking equipment to detect gravitational waves (GWs) and other analogous signals related to dark matter (DM) over the duration of its interplanetary cruise. By employing a methodology to stack tracking data in combination with Monte-Carlo Markov-Chain parameter recovery tests, we show that the mission will be sensitive to GWs over the wide frequency range of $3\times 10^{-9}$ Hz to $10^{-1}$ Hz, provided that tracking data is taken consistently over a large fraction of the cruise duration. Thus, the mission has the potential to fill the gap between pulsar timing and space-based-interferometry GW observatories. Within this assumption, we forecast the detection of $\mathcal{\mathcal{O}}(1 - 100)$ individual massive black hole binaries using two independent population models. Additionally, we determine the mission's sensitivity to both astrophysical and primordial stochastic gravitational wave backgrounds, as well as its capacity to test, or even confirm via detection, ultralight DM models. In all these cases, the tracking of the spacecraft over its interplanetary cruise would enable coverage of unexplored regions of parameter space, where signals from new phenomena in our Universe may be lurking.

The problem of Astrophysical Jet formation from relativistic accretion disks through the establishment of relativistic disk-powerful jet equilibrium structure is studied applying the Beltrami-Bernoulli equilibrium approach of Shatashvili & Yoshida 2011; Arshilava et al. 2019. Accretion disk is weakly magnetized consisting of fully ionized relativistic electron-ion plasma and photon gas strongly coupled to electrons due to Thompson Scattering. %hence, making the behavior of photon gas similar to that of "a charged fluid". Analysis is based on the generalized Shakura-Sunyaev $\alpha $-turbulent dissipation model for local viscosity (being the main source of accretion), in which the contributions from both the photon and ion gases are taken into account. Ignoring the self-gravitation in the disk we constructed the analytical self-similar solutions for the equilibrium relativistic disk-jet structure characteristic parameters in the field of gravitating central compact object for the force-free condition. It is shown, that the magnetic field energy in the Jet is several orders greater compared to that of accretion disk, while jet-outflow is locally Super-Alfvénic with local {\it Plasma-beta} $< 1$ near the jet-axis. The derived solutions can be used to analyze the astrophysical jets observed in binary systems during the star formation process linking the jet properties with the parameters of relativistic accretion disks of electron-ion-photon gas.

Time-domain astrophysics continues to grow rapidly, with the inception of new surveys drastically increasing data volumes. Democratised, distributed approaches to training sets for machine learning classifiers are crucial to make the most of this torrent of discovery -- with citizen science approaches proving effective at meeting these requirements. In this paper, we describe the creation of and the initial results from the $\textit{Kilonova Seekers}$ citizen science project, built to find transient phenomena from the GOTO telescopes in near real-time. $\textit{Kilonova Seekers}$ launched in July 2023 and received over 600,000 classifications from approximately 2,000 volunteers over the course of the LIGO-Virgo-KAGRA O4a observing run. During this time, the project has yielded 20 discoveries, generated a `gold-standard' training set of 17,682 detections for augmenting deep-learned classifiers, and measured the performance and biases of Zooniverse volunteers on real-bogus classification. This project will continue throughout the lifetime of GOTO, pushing candidates at ever-greater cadence, and directly facilitate the next-generation classification algorithms currently in development.

In previous works, we have investigated the star-forming region G29.96$-$0.02 where the massive young stellar object (MYSO) G29.862$-$0.0044 (hereafter G29) is embedded in a hot molecular core. In one of them, of multiwavelength nature, using data from the Atacama Submillimeter Telescope Experiment (ASTE), data from the Atacama Large Millimeter Array (ALMA), and photometric data from NIRI-Gemini, G29 was investigated at different spatial scales. However, the intriguing morphology of G29 in the near-infrared, together with the distribution of the associated molecular gas, reveals that the star-formation scenario is far from being understood. This work incorporates the analysis of the emission of several molecular lines acquired with ALMA that were not previously examined (eg.,~CH$_{3}$OH, HC$_{3}$N, H$_{ 2}$CO, C$^{34}$S, H$_{2}$CS) as well as a new determination of the temperature of the region. Additionally, we present the progress of results obtained through new observations in the near-infrared, in this case spectroscopic, using NIFS-Gemini, and in radio continnum obtained with the Karl G. Jansky Very Large Array (JVLA). This research allows us to carry out a detailed chemical study of the region, which will contribute to the understanding of the physical processes involved in the high-mass star formation.

Third generation ground-based gravitational wave (GW) detectors, such as Einstein Telescope and Cosmic Explorer, will operate in the $(\text{few}-10^4)$ Hz frequency band, with a boost in sensitivity providing an unprecedented reach into primordial cosmology. Working concurrently with pulsar timing arrays in the nHz band, and LISA in the mHz band, these 3G detectors will be powerful probes of beyond the standard model particle physics on scales $T\gtrsim 10^{7}$GeV. Here we focus on their ability to probe phase transitions (PTs) in the early universe. We first overview the landscape of detectors across frequencies, discuss the relevance of astrophysical foregrounds, and provide convenient and up-to-date power-law integrated sensitivity curves for these detectors. We then present the constraints expected from GW observations on first order PTs and on topological defects (strings and domain walls), which may be formed when a symmetry is broken irrespective of the order of the phase transition. These constraints can then be applied to specific models leading to first order PTs and/or topological defects. In particular we discuss the implications for axion models, which solve the strong CP problem by introducing a spontaneously broken Peccei-Quinn (PQ) symmetry. For post-inflationary breaking, the PQ scale must lie in the $10^{8}-10^{11}$ GeV range, and so the signal from a first order PQ PT falls within reach of ground based 3G detectors. A scan in parameter space of signal-to-noise ratio in a representative model reveals their large potential to probe the nature of the PQ transition. Additionally, in heavy axion type models domain walls form, which can lead to a detectable GW background. We discuss their spectrum and summarise the expected constraints on these models from 3G detectors, together with SKA and LISA.

The Wide-Field Spectroscopic Telescope (WST) is a concept for a 12-m class seeing-limited telescope providing two concentric fields of view for simultaneous Multi-Object Spectroscopy and Integral Field Spectroscopy. The specified wavelength range is 0.35-1.6 microns. The baseline optical design relies on a corrected Cassegrain solution feeding Multi-Object spectrographs through fibres, while the central area of the field is propagated down to a gravity-stable Integral Field Station housing 144 spectrographs. The Cassegrain corrector also provides for atmospheric dispersion compensation. All optical components are within commercially available dimensions. With a view to minimizing risks and costs, to the maximum possible extent the telescope relies on proven subsystem solutions. An exception is the tip-tilt secondary mirror, which would likely have to provide some rejection of wind shake. An iteration of the optical design is ongoing, with a view to mitigating the weaknesses of the first baseline design. The telescope would be wavefront-controlled on-sky at the common-path MOS focus. Controls in the IFS path will need to compensate for the effect of subsequent differentials - wavefront and line of sight. There is no shortage of degrees of freedom and metrology solution to do so. The size of the dome is driven by the Nasmyth footprint and the height of the pier, which houses the IFS station. The baseline assumption is that a VLT-like enclosure would provide suitable shielding and ventilation.

Astrochemical models are important tools to interpret observations of molecular and atomic species in different environments. However, these models are time-consuming, precluding a thorough exploration of the parameter space, leading to uncertainties and biased results. Using neural networks to simulate the behavior of astrochemical models is a way to circumvent this problem, providing fast calculations that are based on real astrochemical models. In this paper, we present a fast neural emulator of the astrochemical code Nautilus based on conditional neural fields. The resulting model produces the abundance of 192 species for arbitrary times between 1 and 10$^7$ years. Uncertainties well below 0.2 dex are found for all species, while the computing time is of the order of 10$^4$ smaller than Nautilus. This will open up the possibility of performing much more complex forward models to better understand the physical properties of the interstellar medium. As an example of the power of these models, we ran a feature importance analysis on the electron abundance predicted by Nautilus. We found that the electron density is coupled to the initial sulphur abundance in a low density gas. Increasing the initial sulphur abundance from a depleted scenario to the cosmic abundance leads to an enhancement of an order of magnitude of the electron density. This enhancement can potentially influence the dynamics of the gas in star formation sites.

We have developed a novel method of determining 2D radial density profiles for astronomical systems of discrete objects using Voronoi tessellations. This Voronoi-based method was tested against the standard annulus-based method on 5 simulated systems of objects, following known Hubble density profiles of varying parameters and sizes. It was found that the Voronoi-based method returned radial density fits with lower uncertainties on the fitting parameters across all 5 systems compared to the annulus-based method. The Voronoi-based method also consistently returned more accurate estimates of the total number of objects in each system than the annulus-based method, and this accuracy increased with increasing system size. Finally, the Voronoi-based method was applied to two observed globular cluster systems around brightest cluster galaxies ESO 444-G046 and 2MASX J13272961-3123237 and the results were compared to previous results for these galaxies obtained with the annulus-based method. Again, it was found that the Voronoi-based method returned fits with lower uncertainties on the fitting parameters, and the total number of globular clusters returned are within errors of the annulus-based method estimates, however also with lower uncertainties.

Asteroseismology gives us the opportunity to look inside stars and determine their internal properties. Based on these observations, estimations can be made for the amount of the convective boundary mixing and envelope mixing of such stars, and the shape of the mixing profile in the envelope. However, these results are not typically included in stellar evolution models. We aim to investigate the impact of varying convective boundary mixing and envelope mixing in a range based on asteroseismic modelling in stellar models, both for the stellar structure and for the nucleosynthetic yields. In this first study, we focus on the pre-explosive evolution of a 20Msun star and evolve the models to the final phases of carbon burning. We vary the convective boundary mixing, implemented as step-overshoot, with the overshoot parameter in the range 0.05-0.4 and the amount of envelope mixing in the range 1-10$^{6}$ with a mixing profile based on internal gravity waves. We use a large nuclear network of 212 isotopes to study the nucleosynthesis. We find that enhanced mixing according to asteroseismology of main-sequence stars, both at the convective core boundary and in the envelope, has significant effects on the nucleosynthetic wind yields. Our evolutionary models beyond the main sequence diverge in yields from models based on rotational mixing, having longer helium burning lifetimes and lighter helium-depleted cores. We find that the asteroseismic ranges of internal mixing calibrated from core hydrogen burning stars lead to similar wind yields as those resulting from the theory of rotational mixing. Adopting the seismic mixing levels beyond the main sequence, we find earlier transitions to radiative carbon burning compared to models based on rotational mixing. This influences the compactness and the occurrence of shell-mergers, which may affect the supernova properties and explosive nucleosynthesis.

The 21cm line of neutral hydrogen is a powerful probe of the high-redshift universe (Cosmic Dawn and the Epoch of Reionization), with an unprecedented potential to inform us about key processes of early galaxy formation, the first stars and even cosmology and structure formation, via intensity mapping. It is the subject of a number of current and upcoming low-frequency radio experiments. This paper presents 21cmSense v2.0, which is a Python package that provides a modular framework for calculating the sensitivity of these experiments, in order to enhance the process of their design and forecasting their power for parameter inference. Version 2.0 of 21cmSense has been re-written from the ground up to be more modular and extensible than its venerable predecessor (Pober et al., 2013, 2014), and to provide a more user-friendly interface. The package is freely available both to use and contribute towards at this https URL.

We show, both analytically and numerically, that non-Gaussian tails in the probability density function of curvature perturbations arise in ultra-slow-roll inflation from the $\delta N$ formalism, without invoking stochastic inflation. Previously reported discrepancies between both approaches are a consequence of not correctly accounting for momentum perturbations. Once they are taken into account, both approaches agree to an excellent degree. The shape of the tail depends strongly on the phase space of inflation.

There is a growing number of peculiar events that cannot be assigned to any of the main supernova (SN) classes. SN 1987A and a handful of similar objects, thought to be explosive outcomes of blue supergiant stars, belong to them: while their spectra closely resemble those of H-rich (IIP) SNe, their light-curve (LC) evolution is very different. Here we present the detailed photometric and spectroscopic analysis of SN 2021aatd, a peculiar Type II explosion: while its early-time evolution resembles that of the slowly evolving, double-peaked SN 2020faa (however, at a lower luminosity scale), after $\sim$40 days, its LC shape becomes similar to that of SN 1987A-like explosions. Beyond comparing LCs, color curves, and spectra of SN 2021aatd to that of SNe 2020faa, 1987A, and of other objects, we compare the observed spectra with our own SYN++ models and with the outputs of published radiative transfer models. We also modeled the pseudo-bolometric LCs of SNe 2021aatd and 1987A assuming a two-component (core+shell) ejecta, and involving the rotational energy of a newborn magnetar in addition to radioactive decay. We find that both the photometric and spectroscopic evolution of SN 2021aatd can be well described with the explosion of a $\sim$15 $M_\odot$ blue supergiant star. Nevertheless, SN 2021aatd shows higher temperatures and weaker Na ID and Ba II 6142 A lines than SN 1987A, which is reminiscent of rather to IIP-like atmospheres. With the applied two-component ejecta model (counting with both decay and magnetar energy), we can successfully describe the bolometric LC of SN 2021aatd, including the first $\sim$40-day long phase showing an excess compared to 87A-like SNe but being strikingly similar to that of the long-lived SN 2020faa. Nevertheless, finding a unified model that also explains the LCs of more luminous events (like SN 2020faa) is still a matter of concern.

Modified gravity (MG) theories have emerged as a promising alternative to explain the late-time acceleration of the Universe. However, the detection of MG in observations of the large-scale structure remains challenging due to the screening mechanisms that obscure any deviations from General Relativity (GR) in high-density regions. The marked two-point correlation function offers a promising approach to potentially detect MG signals. This work investigates novel marks based on large-scale environment estimates but also that exploit the anti-correlation between objects in low- and high-density regions. This is the first time discreteness effects in density-dependent marked correlation functions are investigated in depth. We assess the performance of various marks to distinguish GR from MG by using the ELEPHANT simulations, comprised of realisations of GR as well as $f(R)$ and nDGP gravity. In addition, discreteness effects are studied using the high-density Covmos catalogues. We establish a robust method to correct for shot-noise effects that allows the recovery of the true signal with an accuracy below $5\%$ over a wide range of scales. We find such correction to be crucial to measure the amplitude of the marked correlation function in an unbiased manner. Furthermore, we demonstrate that marks, anti-correlating objects in low- and high-density regions, are among the most effective in distinguishing between MG and GR. We report differences in the marked correlation function between $f(R)$ with $|f_{R0}|=10^{-6}$ and GR simulations of the order of 3-5$\sigma$ in real space up to scales of about $80\, h^{-1} \, {\rm Mpc}$. The redshift-space monopole exhibits similar features and performances. The combination of the proposed $\tanh$-mark with shot-noise correction paves the way towards an optimal approach for the detection of MG in current and future galaxy spectroscopic surveys.

We construct a hybrid-inflation model where the inflaton potential is generated radiatively, as gauge symmetries guarantee it to be accidentally flat at tree level. The model can be regarded as a small-field version of Natural Inflation, with inflation ending when the mass of a second scalar, the waterfall field, turns tachyonic. This provides a minimal, robust realisation of hybrid inflation, which predicts specific correlations among CMB observables. Tachyonic preheating leads to the production of gravitational waves which, for a low inflationary scale, might be detected by upcoming experiments. Simple variations of the model can give rise to topological defects, such as unstable domain walls. Their dynamics produces a stochastic gravitational-wave background, which can be compatible with the recent detection by pulsar timing arrays.

Over two hundred protoplanetary disk systems have been resolved by ALMA, and the vast majority suggest the presence of planets. The dust gaps in transition disks are considered evidence of giant planets sculpting gas and dust under appropriate disk viscosity. However, the unusually high accretion rates in many T Tauri stars hosting transition disks challenge this theory. As the only disk currently observed with high turbulence, the high accretion rate ($\sim10^{-8.3}M_{\odot}/yr$) observed in DM Tau indicates the presence of strong turbulence may within the system. Considering the recent theoretical advancements in magnetized disk winds is challenging the traditional gap-opening theories and viscosity-driven accretion models, our study presents a pioneering simulation incorporating a simplified magnetized disk wind model to explain the observed features in DM Tau. Employing multi-fluid simulations with an embedded medium mass planet, we successfully replicate the gap formation and asymmetric structures evident in ALMA Band 6 and the recently JVLA 7 mm observations. Our results suggest that when magnetized disk wind dominate the accretion mode of the system, it's entirely possible for a planet with a medium mass to exist within the gap inside 20 au of DM Tau. This means that DM Tau may not be as turbulence as imagined. However, viscosity within the disk should also contribute a few turbulence to maintain disk stability.