Papers for Friday, Sep 06 2024

Vaporized metal, silicates, and ices on the verge of re-condensing into solid or liquid particles appear in many contexts: behind shocks, in impact ejecta, and within the atmospheres and outflows of stars, disks, planets, and minor bodies. We speculate that a condensing gas might fragment, forming overdensities within relative voids, from a radiation-condensation instability. Seeded with small thermal fluctuations, a condensible gas will exhibit spatial variations in the density of particle condensates. Regions of higher particle density may radiate more, cooling faster. Faster cooling leads to still more condensation, lowering the local pressure. Regions undergoing runaway condensation may collapse under the pressure of their less condensed surroundings. Particle condensates will compactify with collapsing regions, into overdense clumps or macroscopic solids (planetesimals). As a first step toward realizing this hypothetical instability, we calculate the evolution of a small volume of condensing silicate vapor -- a spherical test "bubble" embedded in a background medium whose pressure and radiation field are assumed fixed for simplicity. Such a bubble condenses and collapses upon radiating its latent heat to the background, assuming its energy loss is not stopped by background irradiation. Collapse speeds can range up to sonic, similar to cavitation in terrestrial settings. Adding a non-condensible gas like hydrogen to the bubble stalls the collapse. We discuss whether cavitation can provide a way for mm-sized chondrules and refractory solids to assemble into meteorite parent bodies, focusing on CB/CH chondrites whose constituents likely condensed from silicate/metal vapor released from the most energetic asteroid collisions.

Precise and accurate mass and radius measurements of evolved stars are crucial to calibrating stellar models. Stars in detached eclipsing binaries (EBs) are excellent potential calibrators because their stellar parameters can be measured with fractional uncertainties of a few percent, independent of stellar models. The All-Sky Automated Survey for Supernovae (ASAS-SN) has identified tens of thousands of EBs, >35,000 of which were included in the ASAS-SN eclipsing binaries catalog. Here, we select eight EBs from this sample that contain giants based on their Gaia colors and absolute magnitudes. We use LBT/PEPSI, APF, and CHIRON to obtain multi-epoch spectra of these binaries and measure their radial velocities using two-dimensional cross-correlation methods. We simultaneously fit the ASAS-SN light curves and the radial velocities with PHOEBE to derive accurate and precise masses and radii with fractional uncertainties of $\lesssim 3\%$. For four systems, we also include Transiting Exoplanet Survey Satellite (TESS) light curves in our PHOEBE models, which significantly improves the radius determinations. In seven of our systems, both components have evolved off of the main sequence, and one system has a giant star component with a main sequence, Sun-like companion. Finally, we compare our mass and radius measurements to single-star evolutionary tracks and distinguish between systems that are first ascent red giant branch stars and those that are likely core helium-burning stars.

Cold gas and dust reprocess the X-ray emission of active galactic nuclei (AGN), producing spectro-polarimetric features in the X-ray band. The recent launch of IXPE allows for observations of this X-ray polarisation signal, encoding unique information on the circumnuclear medium of AGN. However, the models for interpreting these polarimetric data are under-explored and do not reach the same level of sophistication as the corresponding spectral models. We aim at closing the gap between the spectral and spectro-polarimetric modelling of AGN circumnuclear media by providing the tools for simulating X-ray polarisation in complex 3D transfer media alongside X-ray spectra. We lay out the framework for X-ray polarisation in 3D radiative transfer simulations and provide an implementation to the 3D radiative transfer code SKIRT, focussing on (de)polarisation due to scattering and fluorescent re-emission. As a first application, we studied a 2D toroidal reprocessor of cold gas, modelling the AGN circumnuclear medium. For the 2D torus model, we find a complex behaviour of the polarisation angle with photon energy, which we interpret as a balance between the reprocessed photon flux originating from different sky regions, with a direct link to the torus geometry. We calculated a large grid of AGN torus models and demonstrated how spatially resolved polarisation maps could form a useful tool for interpreting the geometrical information that is encoded in IXPE observations. With this work, we release high-resolution AGN torus templates that simultaneously describe X-ray spectra and spectro-polarimetry, for observational data fitting with XSPEC. The SKIRT code can now model X-ray polarisation simultaneously with X-ray spectra and provide synthetic spectro-polarimetric observations for complex 3D circumnuclear media, with all features of the established SKIRT framework available.

We summarize all available 13CO and accompanying 12CO measurements of local galaxy centers in transitions of J=5-4 and higher, including new APEX 13CO(6-5) and 12CO(6-5) observations of 11 galaxies. The observed integrated temperature ratios of 12CO(6-5) to 12CO(1-0) range from 0.10 to 0.45. Multi-aperture data indicate that 13CO(6-5) is more centrally concentrated than 12CO(6-5). The emission of 12CO(6-5) and HCO+ but not HCN may be correlated. The new data are essential to constrain the physical properties of the galaxy center molecular gas even in a simple two-phase model approximating the more complex multi-phase structure. In all galaxies, except the Seyfert AGN galaxy NGC 1068, high J emission from the center is dominated by a dense (n = 100 000) and relatively cool (T = 20 to 60 K) high-pressure gas. In contrast, the low-J lines are dominated in most galaxies by low-pressure gas of a moderate density (n = 1000) and more elevated temperature (T = 60 to 150 K). The three exceptions with significant high-pressure gas contributions to the low J emission are all associated with active central star formation.

Anisotropies of the Sunyaev Zel'dovich (SZ) effect serve as a powerful probe of the thermal history of the universe. At high redshift, hot galactic outflows driven by supernovae (SNe) can inject a significant amount of thermal energy into the intergalactic medium, causing a strong $y$-type distortion of the CMB spectrum through inverse Compton scattering. The resulting anisotropies of the $y$-type distortion are sensitive to key physical properties of high-$z$ galaxies pertaining to the launch of energetic SNe-driven outflows, such as the efficiency and the spatio-temporal clustering of star formation. We develop a simple analytic framework to calculate anisotropies of $y$-type distortion associated with SNe-powered outflows of galaxies at $z>6$. We show that galactic outflows are likely the dominant source of thermal energy injection, compared to contributions from reionized bubbles and gravitational heating. We further show that next-generation CMB experiments such as LiteBIRD can detect the contribution to $y$ anisotropies from high-$z$ galactic outflows through the cross-correlation with surveys of Lyman-break galaxies by e.g. the Roman Space Telescope. Our analysis and forecasts demonstrate that thermal SZ anisotropies are a promising probe of SNe feedback in early star-forming galaxies.

Recent JWST observations at $z > 6$ may imply galactic ionizing photon production in excess of prior expectations. Under observationally motivated assumptions about escape fractions, these suggest a $z \sim 8-9$ end to reionization, in strong tension with the $z < 6$ end required by the Ly$\alpha$ forest. In this work, we use radiative transfer simulations to understand what different observations tell us about when reionization ended and when it started. We consider a model that ends too early (at $z \approx 8$) alongside two more realistic scenarios that end late at $z \approx 5$: one that starts late ($z \sim 9$) and another that starts early ($z \sim 13$). We find that the latter requires up to an order-of-magnitude evolution in galaxy ionizing properties at $6 < z < 12$, perhaps in tension with recent measurements of $\xi_{\rm ion}$ by JWST, which indicate little evolution. We also study how these models compare to recent measurements of the Ly$\alpha$ forest opacity, mean free path, IGM thermal history, visibility of $z > 8$ Ly$\alpha$ emitters, and the patchy kSZ signal from the CMB. We find that neither of the late-ending scenarios is conclusively disfavored by any single data set. However, a majority of these observables, spanning several distinct types of observations, prefer a late start. Not all probes agree with this conclusion, hinting at a possible lack of concordance between observables. Observations by multiple experiments (including JWST, Roman, and CMB-S4) in the coming years will either establish a concordance picture of reionization's early stages or reveal systematics in data and/or theoretical modeling.

K2-24 is a planetary system composed of two transiting low-density Neptunians locked in an almost perfect 2:1 resonance and showing large TTVs, i.e., an excellent laboratory to search for signatures of planetary migration. Previous studies performed with K2, Spitzer and RV data tentatively claimed a significant non-zero eccentricity for one or both planets, possibly high enough to challenge the scenario of pure disk migration through resonant capture. With 13 new CHEOPS light curves (seven of planet -b, six of planet -c), we carried out a global photometric and dynamical (RV+TTV) re-analysis by including all the available literature data as well. We got the most accurate set of planetary parameters to date for the K2-24 system, including radii and masses at 1% and 5% precision (now essentially limited by the uncertainty on stellar parameters) and non-zero eccentricities $e_b=0.0498_{-0.0018}^{+0.0011}$, $e_c=0.0282_{-0.0007}^{+0.0003}$ detected at very high significance for both planets. Such relatively large values imply the need for an additional physical mechanism of eccentricity excitation during or after the migration stage. Also, while the accuracy of the previous TTV model had drifted by up to 0.5 days at the current time, we constrained the orbital solution firmly enough to predict the forthcoming transits for the next ~15 years, thus enabling an efficient follow-up with top-level facilities such as JWST or ESPRESSO.

We conduct a study on the relationship between galaxy environments and their active galactic nuclei (AGN) activity at high redshifts ($2.0<z<4.0$). Specifically, we study the AGN fraction in galaxies residing in a range of environments at these redshifts, from field galaxies to highly overdense peaks in the GOODS-S extragalactic field. Utilizing the extensive photometric and spectroscopic observations in this field, we measure local- and global-overdensities over a large a range of environments, including in several massive (M$_{tot}\geq10^{14.8}$M$_\odot$) protostructures. We employ a multi-wavelength AGN catalog consisting of AGN in nine different categories. Our analysis shows a higher AGN fraction (10.9$^{+3.6}_{-2.3}$%) for galaxies in the highest local-overdensity regions compared to the AGN fraction (1.9$^{+0.4}_{-0.3}$%) of coeval field galaxies (a ~4$\sigma$ difference). This trend of increasing AGN fraction in denser environments relative to the field is present in all redshift bins. We also find this trend consistently present in all five AGN categories that have a sufficient number of AGN to make a meaningful comparison: mid-IR SED, mid-IR color, X-ray luminosity, X-ray-luminosity-to-radio-luminosity-ratio, and optical-spectroscopy. Our results also demonstrate a clear trend of higher (~4x) AGN fractions in denser environments for a given stellar mass. Additionally, we observe the same trend (though at a lower significance) with the global environment of galaxies, measured using a metric based on the projected distance of galaxies from their nearest massive ($M_{tot}>10^{12.8}M_\odot$) overdense ($\sigma_\delta>5.0$) peak, normalized with respect to the size of the peak. These findings indicate that the prevalence of AGN activity is highly dependent on the environment in which a host galaxy resides, even at early times in the formation history of the Universe.

Close to the Roche radius of a white dwarf (WD), an asteroid on a circular orbit sheds material that then adopts a very similar orbit. Observations of the resulting debris show a periodic behavior and changes in flux on short timescales, implying ongoing dynamical activity. Additional encounters from other minor planets may then yield co-orbital rings of debris at different inclinations. The structure, dynamics, and lifetime of these debris discs remains highly uncertain, but is important for understanding WD planetary systems. We aim to identify and quantify the locations of co-orbitals in WD-asteroid-dust particle 3-body systems by exploring the influence of 1:1 resonant periodic orbits. We begin this exploration with co-planar and inclined orbits in the circular restricted 3-body problem (CRTBP) and model the dynamical evolution of these exosystems over observable timescales. The mass ratio parameter for this class of systems ($~2\times 10^{-11}$) is one of the lowest ever explored in this dynamical configuration. We computed the periodic orbits, deduced their linear stability, and suitably seeded the dynamical stability maps. We carried out a limited suite of N-body simulations to provide direct comparisons with the maps. We derive novel results for this extreme mass ratio in the CRTBP, including new unstable 3D families. We illustrate through the maps and N-body simulations where dust can exist in a stable configuration over observable timescales across a wide expanse of parameter space in the absence of strong external forces. Over a timescale of 10 yr, the maximum orbital period deviations of stable debris due to the co-orbital perturbations of the asteroid is about a few seconds. Unstable debris in a close encounter with the asteroid typically deviates from the co-orbital configuration by more than about 20 km and is on a near-circular orbit with an eccentricity lower than ~0.01.

We report evidence of galaxy assembly bias - the correlation between galaxy properties and biased secondary halo properties at fixed halo mass (M$_H$) - in the stellar-to-halo mass relation (SHMR) for red central galaxies from the Sloan Digital Sky Survey. In the M$_H = 10^{11.5}-10^{13.5} h^{-1}$ M$_{\odot}$ range, central galaxy stellar mass (M$_*$) is correlated with the number density of galaxies within $10 h^{-1}$ Mpc ($\delta_{10}$), a common proxy for halo formation time. This galaxy assembly bias signal is also present when M$_H$, M$_*$, and $\delta_{10}$ are substituted with group luminosity, galaxy luminosity, and metrics of the large-scale density field. To associate differences in $\delta_{10}$ with variations in halo formation time, we fitted a model that accounts for (1) errors in the M$_H$ measured by the Tinker 2021, 2022 group catalog and (2) the level of correlation between halo formation time and M$_*$ at fixed M$_H$. Fitting of this model yields that (1) errors in M$_H$ are 0.15 dex and (2) halo formation time and M$_*$ are strongly correlated (Spearman's rank correlation coefficient ~0.85). At fixed M$_H$, variations of ~0.4 dex in M$_*$ are associated with ~1-3 Gyr variations in halo formation time and in galaxy formation time (from stellar population fitting; Oyarzún et al. 2022). These results are indicative that halo properties other than M$_H$ can impact central galaxy assembly.

Much of modern cosmology relies on the Cosmological Principle, the assumption that the Universe is isotropic and homogeneous on sufficiently large scales, but it remains worthwhile to examine cosmological models that violate this principle slightly. We examine a class of such spacetimes that maintain homogeneity but break isotropy through their underlying local spatial geometries. These spacetimes are endowed with one of five anisotropic model geometries of Thurston's geometrization theorem, and their evolution is sourced with perfect fluid dust and cosmological constant. We show that the background evolution of these spacetimes induces fluctuations in the observed cosmic microwave background (CMB) temperature with amplitudes coupled to the curvature parameter $\Omega_K$. In order for these fluctuations to be compatible with the observed CMB angular power spectrum, we find $|\Omega_K| \lesssim 10^{-5}$ is required in all five geometries, with two geometries requiring $|\Omega_K| \lesssim 10^{-10}$. This strongly limits the cosmological consequences of these models.

The IW And phenomenon in cataclysmic variables presents a significant challenge to the accretion disk instability model. Using photometric data from the All-Sky Automated Survey for Supernovae, the Zwicky Transient Facility, and the Transiting Exoplanet Survey Satellite, we identify Karachurin 12 as a new IW And-type object with a cycle period of 35.69(3) days. We also report for the first time that Karachurin 12 is a negative superhump (NSH) system featuring a precessing tilted disk, with precession, orbital, and NSH periods of 4.9588(2) days, 0.3168895(13) days, and 0.2979861(8) days, respectively. Our analysis, using dips as indexes and NSHs as probes, reveals diverse cycle patterns in Karachurin 12, with NSH amplitude varying throughout the cycle. Analyses of NSH behaviour using multiple methods shows that the amplitude decreases with increasing outbursts and increases with weakening outbursts, which we suggest is related to changes in the radius of the accretion disk. The mass transfer outburst model fails to explain the observed NSH amplitude variations, whereas the tilted thermally unstable disk model accounts for most of the observed cycles and dip patterns, particularly the disappearance of the IW And phenomenon when NSHs become undetectable. However, no significant changes in the NSH period were observed. These findings suggest a potential link between the IW And phenomenon and a tilted disk. Further improvements in simulations are needed to better understand the detailed dynamics of these cycles.

We report the discovery of jets emanating from pre-main-sequence objects exclusively at mid-infrared wavelengths, enabled by the superb sensitivity of JWST's Mid-InfraRed Medium-Resolution Spectrometer (MIRI MRS) instrument. These jets are observed only in lines of [NiII], [FeII], [ArII], and [NeII]. The H$_2$ emission, imaged in eight distinct transitions, has a completely different morphology, exhibiting a wide-angled, biconical shape, symmetrically distributed about the jet axes. Synergistic high-resolution Atacama Large Millimeter/submillimeter Array (ALMA) observations resolve a pair of side-by-side edge-on accretion disks lying at the origin of the twin mid-infrared jets. Assuming coevality of the components of the young multiple system under investigation, the system age is at least (2 $-$ 2.5) $\times$ 10$^6$ yr, despite the discrepantly younger age inferred from the spectral energy distribution of the combined edge-on disk sources. The later system evolutionary stage is corroborated by ALMA observations of CO(2$-$1), $^{13}$CO(2$-1$), and C$^{18}$O(2$-$1), which show no traces of molecular outflows or remnant cavity walls. Consequently, the observed H$_2$ structures must have their origins in wide-angled disk winds, in the absence of any ambient, swept-up gas. In the context of recent studies of protostars, we propose an outflow evolutionary scenario in which the molecular gas component dominates in the youngest sources, whereas the fast, ionized jets dominate in the oldest sources, as is the case for the twin jets discovered in the WL 20 system.

GEO 600, Kagra, LIGO, and Virgo were built to observe gravitational waves at frequencies in the audio band, where the highest event rates combined with the largest signal to noise ratios had been predicted. Currently, hypothetical sources of cosmological origin that could have produced signals at higher frequencies are under discussion. What is not widely known is that current interferometric GW observatories have a frequency comb of high optical sensitivity that encompasses these high frequencies. Here we calculate the high-frequency noise spectral densities of operating GW observatories under the justified assumption that photon shot noise is the dominant noise source. We explain the underlying physics of why high sensitivity is achieved for all integer multiples of the free spectral ranges of the observatory's resonators when an interferometer arm is not orientated perpendicular to the propagation direction of the GW. Proposals for new concepts of high-frequency GW detectors must be compared with the high-frequency sensitivities presented here.

Basilisk is a novel Bayesian hierarchical method for inferring the galaxy-halo connection, including its scatter, using the kinematics of satellite galaxies extracted from a redshift survey. In this paper, we introduce crucial improvements, such as updated central and satellite selection, advanced modelling of impurities and interlopers, extending the kinematic modelling to fourth order by including the kurtosis of the line-of-sight velocity distribution, and utilizing satellite abundance as additional constraint. This drastically enhances Basilisk's performance, resulting in an unbiased recovery of the full conditional luminosity function (central and satellite) and with unprecedented precision. After validating Basilisk's performance using realistic mock data, we apply it to the SDSS-DR7 data. The resulting inferences on the galaxy-halo connection are consistent with, but significantly tighter than, previous constraints from galaxy group catalogues, galaxy clustering and galaxy-galaxy lensing. Using full projected phase-space information, Basilisk breaks the mass-anisotropy degeneracy, thus providing precise global constraint on the average orbital velocity anisotropy of satellite galaxies across a wide range of halo masses. Satellite orbits are found to be mildly radially anisotropic, in good agreement with the mean anisotropy for subhaloes in dark matter-only simulations. Thus, we establish Basilisk as a powerful tool that is not only more constraining than other methods on similar volumes of data, but crucially, is also insensitive to halo assembly bias which plagues the commonly used techniques like galaxy clustering and galaxy-galaxy lensing.

Due to the unprecedented signal strengths offered by the newest high-resolution spectrographs on 10-m class telescopes, exploring the 3D nature of exoplanets is possible with an unprecedented level of precision. In this paper, we present a new technique to probe the vertical structure of exoplanetary winds and dynamics using ensembles of planet absorption lines of varying opacity, and apply it to the well-studied ultra-hot Jupiter WASP-76b. We then compare these results to state-of-the-art global circulation models (GCMs) with varying magnetic drag prescriptions. We find that the known asymmetric velocity shift in Fe I absorption during transit persists at all altitudes, and observe tentative trends for stronger blueshifts and more narrow line profiles deeper in the atmosphere. By comparing three different model prescriptions (a hydrodynamical model with no drag, a magnetic drag model, and a uniform drag model) we are able to rule out the uniform drag model due to inconsistencies with observed trends in the data. We find that the magnetic model is slightly favored over the the hydrodynamic model, and note that this 3-Gauss kinematic magnetohydrodynamical GCM is also favored when compared to low-resolution data. Future generation high-resolution spectrographs on Extremely large telescopes (ELTs) will greatly increase signals and make methods like these possible with higher precision and for a wider range of objects.

In physics, accurately modeling large-scale phenomena such as core-collapse supernovae, (CCSN), and neutron star mergers are computationally challenging and require solving large sets of partial and ordinary differential equations. Traditional methods used widely in the scientific community are predominantly implicit, which are approximations that often require drastic simplifications and can be computationally inefficient. This thesis presents results on a new software suite titled "Fast Explicit Neutrino Networks" or "FENN", that introduces a suite of algebraically stabilized explicit methods known as explicit asymptotic for modeling Neutrino Electron Scattering, (NES), presenting a novel approach that combines the stability of traditional methods with enhanced computational efficiency. Initial results show that FENN can deliver accurate solutions for neutrino networks at improved computational speeds. This thesis further covers new results for scaled networks beyond the constraints of standard energy groupings, as well as the dynamics of neutrino interactions such as the scattering of various neutrino flavors -- electron neutrinos ($\nu_e$), electron anti-neutrinos ($\bar{\nu}_e$), and muon/tau neutrinos ($\nu_{\mu,\tau}$) as well as their anti-particles ($\bar{\nu}_{\mu,\tau}$) -- off electrons.

Globular clusters (GCs) are among the oldest objects in the Universe for which an age can be directly measured, thus playing an important cosmological role. This age, on the other hand, depends sensitively on the He abundance, which cannot be reliably measured from spectroscopy in GC stars. Detached eclipsing binaries (DEBs) near the turnoff (TO) point may play an important role in this regard. The aim of this study is to explore the possibility that, by working with differential measurements of stars that comprise a TO binary system, and assuming both stars have the same age and He abundance, one can achieve tighter, more robust, and less model-dependent constraints on the latter two quantities than otherwise possible by working with the absolute parameters of the stars. We compare both absolute and differential parameters of the stars in V69, a TO DEB pair in the GC 47 Tuc, with two different sets of stellar evolutionary tracks, making use of a Monte Carlo technique to estimate its He abundance and age, along with their uncertainties. We find that the relative approach can produce age and He abundance estimates that are in good agreement with those from the literature. We show that our estimates are also less model-dependent, less sensitive to [Fe/H], and more robust to inherent model systematics than those obtained with an absolute approach. On the other hand, the relative analysis finds larger statistical uncertainties than does its absolute counterpart, at least in the case of V69, where both stars have very similar properties. For binary pairs in which one of the components is less evolved than the other, the statistical uncertainty can be reduced. Our study suggests that the method proposed in this work may be useful to robustly constrain the He abundance and ages of GCs.

Void number counts (VNC) indicates the number of low-density regions in the large-scale structure (LSS) of the Universe, and we propose to use it as an effective cosmological probe. By generating the galaxy mock catalog based on Jiutian simulations and considering the spectroscopic survey strategy and instrumental design of the China Space Station Telescope (CSST), which can reach a magnitude limit $\sim$23 AB mag and spectral resolution $R\gtrsim200$ with a sky coverage 17,500 deg$^2$, we identify voids using the watershed algorithm without any assumption of void shape, and obtain the mock void catalog and data of the VNC in six redshift bins from $z=0.3$ to1.3. We use the Markov Chain Monte Carlo (MCMC) method to constrain the cosmological and VNC parameters. The void linear underdensity threshold $\delta_{\rm v}$ in the theoretical model is set to be a free parameter at a given redshift to fit the VNC data and explore its redshift evolution. We find that, the VNC can correctly derive the cosmological information, and the constraint strength on the cosmological parameters is comparable to that from the void size function (VSF) method, which can reach a few percentage levels in the CSST full spectroscopic survey. This is because that, since the VNC is not sensitive to void shape, the modified theoretical model can match the data better by integrating over void features, and more voids could be included in the VNC analysis by applying simpler selection criteria, which will improve the statistical significance. It indicates that the VNC can be an effective cosmological probe for exploring the LSS.

The north-south (NS) anisotropy of galactic cosmic rays (GCRs) is dominated by a diamagnetic drift flow of GCRs in the interplanetary magnetic field (IMF), allowing us to derive key parameters of cosmic-ray propagation, such as the density gradient and diffusion coefficient. We propose a new method to analyze the rigidity spectrum of GCR anisotropy and reveal a solar cycle variation of the NS anisotropy's spectrum using ground-based muon detectors in Nagoya, Japan, and Hobart, Australia. The physics-based correction method for the atmospheric temperature effect on muons is used to combine the different-site detectors free from local atmospheric effects. NS channel pairs in the multi-directional muon detectors are formed to enhance sensitivity to the NS anisotropy, and in this process, general graph matching in graph theory is introduced to survey optimized pairs. Moreover, Bayesian estimation with the Gaussian process allows us to unfold the rigidity spectrum without supposing any analytical function for the spectral shape. Thanks to these novel approaches, it has been discovered that the rigidity spectrum of the NS anisotropy is dynamically varying with solar activity every year. It is attributed to a rigidity-dependent variation of the radial density gradient of GCRs based on the nature of the diamagnetic drift in the IMF. The diffusion coefficient and mean-free-path length of GCRs as functions of the rigidity are also derived from the diffusion-convection flow balance. This analysis expands the estimation limit of the mean-free-path length into $\le200$ GV rigidity region from $<10$ GV region achieved by solar energetic particle observations.

Generating low-thrust transfer trajectories between Earth and the Near Rectilinear Halo Orbit (NRHO), that is selected for NASA's Gateway, can be challenging due to the low control authority available from the propulsion system and the important operational constraint that the duration of all eclipses has to be less than a prescribed 90-minute threshold. We present a method for generating eclipse-feasible, minimum-time solutions to the aforementioned trajectory design problem using a Lyapunov control law. Coasting is enforced during solar eclipses due to both the Earth and Moon. We used particle swarm optimization to optimize the NRHO insertion date, time of flight, and control law parameters according to a cost function that prioritizes 1) convergence to the target orbit, 2) satisfaction of eclipse-duration constraints, and 3) minimization of time of flight. Trajectories can serve as initial guesses for NASA's high-fidelity trajectory design tools such as Copernicus and GMAT.

We present NIR spectroscopy using MMT/MMIRS for a sample of twenty-nine massive galaxies ($\mathrm{log\ M_* / M_{\odot} \gtrsim10}$) at $\mathrm{z\sim0.7}$ with optical spectroscopy from the LEGA-C survey. Having both optical and NIR spectroscopy at this redshift allows us to measure the full suite of rest-optical strong emission lines, enabling the study of ionization sources and the rest-optical selection of active galactic nuclei (AGN), as well as the measurement of dust-corrected $\mathrm{H\alpha}$-based SFRs. We find that eleven out of twenty-nine galaxies host AGN. We infer the nonparametric star formation histories with the SED fitting code \texttt{Prospector} and classify galaxies as star-forming, green valley, or quiescent based on their most recent sSFRs. We explore the connection between AGN activity and suppressed star formation and find that $89\pm15\%$ of galaxies in the green valley or below host AGN, while only $15\%\pm8\%$ of galaxies above the green valley host AGN. We construct the star-forming main sequence (SFMS) and find that the AGN host galaxies are 0.37 dex below the SFMS while galaxies without detectable AGN are consistent with being on the SFMS. However, when compared to a bootstrapped mass-matched sample, the SFRs of our sample of AGN host galaxies are consistent with the full LEGA-C sample. Based on this mass-matched analysis, we cannot rule out that this suppression of star formation is driven by other processes associated with the higher mass of the AGN sample. We therefore cannot link the presence of AGN activity to the quenching of star formation.

Massive stars play a crucial role in the cosmic dynamics and chemical evolution of galaxies, however our understanding of their evolution and properties remains limited. An accurate determination of stellar parameters is essential for advancing our knowledge. Detached eclipsing binaries (DEBs) are particularly valuable for these determinations due to the minimal interaction between their stellar components, allowing for precise measurements. This study introduces the Young Massive Detached Binary (YMDB) catalog, designed to address the gap in the high-precision absolute parameter determination for young massive stars. By focusing on DEBs within the spectral range O9-B1, this catalog seeks to provide a reliable database for future studies and improve our understanding of massive star evolution. We conducted a photometric analysis of 87 young massive stars in DEBs using TESS light curves (LCs) that were processed through a custom pipeline. This analysis involved determining the amplitude of magnitude variations, orbital periods, times of minima, eccentricities, and the presence of apsidal motion and heartbeat phenomena. A thorough literature review was performed to obtain MK spectral classifications. We performed our own spectral classification of 19 systems to support the sample where a new classification was lacking or inconclusive. The analysis identified 20 previously unreported binary systems, with 13 newly recognized as variables. Among the 87 stars examined, 30 are confirmed as YMDB members, and 25 are candidates pending spectral classification. The remaining 32 stars belong to unsuitable spectral types or nondetached binary nature. Notable findings include the identification of new LC classifications, eccentricities in 13 systems, and heartbeat phenomena in several targets. The YMDB catalog offers a resource of high-quality LCs and reliable stellar classifications.

Observations of the 21 cm signal face significant challenges due to bright astrophysical foregrounds that are several orders of magnitude higher than the brightness of the hydrogen line, along with various systematics. Successful 21 cm experiments require accurate calibration and foreground mitigation. Errors introduced during the calibration process such as systematics, can disrupt the intrinsic frequency smoothness of the foregrounds, leading to power leakage into the Epoch of Reionisation (EoR) window. Therefore, it is essential to develop strategies to effectively address these challenges. In this work, we adopt a stringent approach to identify and address suspected systematics, including malfunctioning antennas, frequency channels corrupted by radio frequency interference (RFI), and other dominant effects. We implement a statistical framework that utilises various data products from the data processing pipeline to derive specific criteria and filters. These criteria and filters are applied at intermediate stages to mitigate systematic propagation from the early stages of data processing. Our analysis focuses on observations from the Murchison Widefield Array (MWA) Phase I configuration. Out of the observations processed by the pipeline, our approach selects 18%, totalling 58 hours, that exhibit fewer systematic effects. The successful selection of observations with reduced systematic dominance enhances our confidence in achieving 21 cm measurements.

1226 ATLASGAL clumps with HII regions were matched with radio sources in the CORNISH-North/South surveys, and 392 of them have corresponding radio sources. We determined the stellar luminosity according to the Lyman continuum flux. When the bolometric luminosity of HII-clumps is less than $\approx$ 10$^{3.7}$ L$_{\odot}$, corresponding to a clump mass $\approx$ 10$^{2.55}$ M$_{\odot}$, the stellar luminosities derived from the Lyman continuum flux overestimate the actual stellar luminosities, because the accretion onto the protostars contributes significantly to the radio emission. After subtracting the accretion luminosity, we obtained reasonable estimates of the stellar luminosity. Using the 0.5 Myr isochrone, we calculated the stellar masses according to the stellar luminosities, and found that they roughly follow the $m_{\rm max}-M_{\rm ecl}$ relation of embedded clusters, consistent with the ionizing sources representing the most massive stars in the embedded clusters of HII-clumps. We also studied the contribution of the possible flaring activity to the observed stellar luminosity and found that they can be neglected. We further studied the change of SFE with the clump mass. According to the derived mass of the most massive star in each HII-clump, using the theoretical $m_{\rm max}-M_{\rm ecl}$ relation, we calculated the mass of the corresponding embedded cluster and then the SFE of the clump. The SFE decreases with increasing clump mass, with a median value of $\approx$0.3. We also independently derived the SFE for each HII-clump based on the model developed in our previous work. The SFEs of HII-clumps derived from the observation and the model are in good agreement. Concerning the star formation histories of the ATLASGAL clumps, low-mass clumps may reach the peak of star formation earlier than high-mass clumps, consistent with the shorter free-fall time of low-mass clumps.

The Varying Speed of Light (VSL) model describes how the speed of light in a vacuum changes with cosmological redshift. Despite numerous models, there is little observational evidence for this variation. While the speed of light can be accurately measured by physical means, cosmological methods are rarely used. Previous studies quantified the speed of light at specific redshifts using Gaussian processes and reconstructed the redshift-dependent function $c(z)$. It is crucial to quantify the speed of light across varying redshifts. We use the latest data on angular diameter distances $D_A(z)$ and Hubble parameters $H(z)$ from baryon acoustic oscillation (BAO) and cosmic chronometer measurements in the redshift interval $z\in[0.07,1.965]$. The speed of light $c(z)$ is determined using Gaussian and deep Gaussian processes to reconstruct $H(z)$, $D_A(z)$, and $D^{\prime}_A(z)$. Furthermore, we conduct comparisons across three distinct models, encompassing two renowned VSL models. We get the result of the parameters constraints in the models (1) for the ``$c$-c" model, $c_0=29492.6 \pm^{6.2}_{5.3} \mathrm{~km} \mathrm{~s}^{-1}$. (2) For the ``$c$-cl" model, $c_0=29665.5 \pm^{11.2}_{11.4}\mathrm{~km} \mathrm{~s}^{-1}$ and $n=0.05535 \pm^{0.00008}_{0.00007}$. (3) For the ``$c$-CPL" model, $c_0=29555.7 \pm^{13.3}_{13.2} \mathrm{~km} \mathrm{~s}^{-1}$ and $n=-0.0607 \pm 0.0001$. Based on our findings, it may be inferred that Barrow's classical VSL model is not a suitable fit for our data. In contrast, the widely recognized Chevallier-Polarski-Linder (CPL) VSL model, under some circumstances, as well as the universal ``c is constant" model, demonstrate a satisfactory ability to account for our findings.

We propose different estimators to probe the epoch of reionization (EoR) intergalactic medium (IGM) using the dispersion measure (${\rm DM}$) of the FRBs. We consider three different reionization histories which we can distinguish with a total of $\lesssim 1000$ ${\rm DM}$ measurements during EoR if their redshifts are known. We note that the redshift derivatives of ${\rm DM}$ are also directly sensitive to the reionization history. The major point of this work is exploring the variance in the ${\rm DM}$ measurements and the information encoded in them. We find that the all-sky average $\overline{\rm DM}(z)$ gets biased from the LoS fluctuations in the ${\rm DM}$ measurements introduced by the ionization of IGM during EoR. We find that the ratio $\sigma_{\rm DM}/\overline{\rm DM}$ depends directly on the ionization bubble sizes as well as the reionization history. On the other hand, we also find that angular variance (coined as $structure$ $function$) of ${\rm DM}$ encodes the information about the duration of reionization and the typical bubble sizes as well. We establish the usefulness of variances in ${\rm DM}$ using toy models of reionization and later verify it with the realistic reionization simulations.

Previous studies have revealed that the Galactic warp is a long-lived, nonsteady, and asymmetric structure. There is a need for a model that accounts for the warp's long-term evolution. Given that this structure has persisted for over 5 Gyrs, its timeline may coincide with the completion of Gaia-Sausage-Enceladus (GSE) merger. Recent studies indicate that the GSE, the significant merger of our Galaxy, was likely a gas-rich merger and the large amount of gas introduced could have created a profound impact on the Galactic morphology. This study utilizes GIZMO simulation code to construct a gas-rich GSE merger. By reconstructing the observed characteristics of the GSE, we successfully reproduce the disk warp and capture nearly all of its documented features that aligns closely with observational data from both stellar and gas disks. This simulation demonstrates the possibility that the single major merger could generate the Galactic warp amplitude and precession. Furthermore, the analysis of the warp's long-term evolution may offer more clues into the formation history of the Milky Way.

The solar wind drag-based model is a widely used framework for predicting the propagation of Coronal Mass Ejections (CMEs) through interplanetary space. This model primarily considers the aerodynamic drag exerted by the solar wind on CMEs. However, factors like magnetic forces, pressure gradients, and the internal dynamics within CMEs justify the need of introducing an additional small-scale acceleration term in the game. Indeed, by accounting for this extra acceleration, the extended drag-based model is shown to offer improved accuracy in describing the evolution of CMEs through the heliosphere and, in turn, in forecasting CME trajectories and arrival times at Earth. This enhancement is crucial for better predicting Space Weather events and mitigating their potential impacts on space-based and terrestrial technologies.

EX Lupi is the prototype of EX Lup-type stars, meaning classical T Tauri stars (cTTSs) showing luminosity bursts and outbursts of 1 to 5 magnitudes lasting for a few months to a few years. These events are ascribed to an episodic accretion that can occur repeatedly but whose physical mechanism is still debated. In this work, we aim to investigate the magnetically-driven accretion of EX Lup in quiescence, including for the first time a study of the small and large-scale magnetic field. This allows us to provide a complete characterisation of the magnetospheric accretion process of the system. We use spectropolarimetric times series acquired in 2016 and 2019 with the Echelle SpectroPolarimetric Device for the Observation of Stars and in 2019 with the SpectroPolarimètre InfraRouge at the Canada-France-Hawaii telescope, during a quiescence phase of EX Lup. We were thus able to perform a variability analysis of the radial velocity, the emission lines and surface averaged longitudinal magnetic field along different epochs and wavelength domains. We also provide a small-scale magnetic field analysis using Zeeman intensification of photospheric lines and large-scale magnetic topology reconstruction using Zeeman-Doppler Imaging. Our study reveals a typical magnetospheric accretion ongoing on EX Lup, with a main accretion funnel flow connecting the inner disc to the star in a stable fashion and producing an accretion shock on the stellar surface close to the pole of the magnetic dipole component. We also measure one of the strongest fields ever observed on cTTSs. Such a strong field indicates that the disc is truncated by the magnetic field close but beyond the corotation radius, where the angular velocity of the disc equals the angular velocity of the star. Such a configuration is suitable for a magnetically-induced disc instability yielding episodic accretion onto the star.

Stars play a decisive role in our Universe, from its beginning throughout its complete evolution. For a thorough understanding of their properties, evolution, and physics of their outer envelopes, stellar spectra need to be analyzed by comparison with numerical models of their atmospheres. We discuss the foundations of how to calculate such models (in particular, density and temperature stratification, affected by convective energy transport in low-mass stars), which requires a parallel treatment of hydrodynamics, thermodynamics and radiative transfer. We stress the impact of emissivities, opacities, and particularly their ratio (source function), and summarize how these quantities are calculated, either adopting or relaxing the assumption of LTE (local thermodynamic equilibrium). Subsequently, we discuss the influence and physics of stellar winds (and their various driving mechanisms as a function of stellar type), rotation, magnetic fields, inhomogeneities, and multiplicity. Finally, we outline the basics of quantitative spectroscopy, namely how to analyze observed spectra in practice.

LSST will provide galaxy cluster catalogs up to z$\sim$1 that can be used to constrain cosmological models once their selection function is well-understood. We have applied the deep convolutional network YOLO for CLuster detection (YOLO-CL) to LSST simulations from the Dark Energy Science Collaboration Data Challenge 2 (DC2), and characterized the LSST YOLO-CL cluster selection function. We have trained and validated the network on images from a hybrid sample of (1) clusters observed in the Sloan Digital Sky Survey and detected with the red-sequence Matched-filter Probabilistic Percolation, and (2) simulated DC2 dark matter haloes with masses $M_{200c} > 10^{14} M_{\odot}$. We quantify the completeness and purity of the YOLO-CL cluster catalog with respect to DC2 haloes with $M_{200c} > 10^{14} M_{\odot}$. The YOLO-CL cluster catalog is 100% and 94% complete for halo mass $M_{200c} > 10^{14.6} M_{\odot}$ at $0.2<z<0.8$, and $M_{200c} > 10^{14} M_{\odot}$ and redshift $z \lesssim 1$, respectively, with only 6% false positive detections. All the false positive detections are dark matter haloes with $ 10^{13.4} M_{\odot} \lesssim M_{200c} \lesssim 10^{14} M_{\odot}$. The YOLO-CL selection function is almost flat with respect to the halo mass at $0.2 \lesssim z \lesssim 0.9$. The overall performance of YOLO-CL is comparable or better than other cluster detection methods used for current and future optical and infrared surveys. YOLO-CL shows better completeness for low mass clusters when compared to current detections in surveys using the Sunyaev Zel'dovich effect, and detects clusters at higher redshifts than X-ray-based catalogs. The strong advantage of YOLO-CL over traditional galaxy cluster detection techniques is that it works directly on images and does not require photometric and photometric redshift catalogs, nor does it need to mask stellar sources and artifacts.

Unlike the Solar System planets, thousands of smaller bodies beyond Neptune orbit the Sun on eccentric ($e >$ 0.1) and ($i>$ 3$^\circ$) orbits. While migration of the giant planets during the early stages of Solar System evolution can induce substantial scattering of trans-Neptunian objects (TNO), this process cannot account for the small number of distant TNOs ($r_p >$ 60 au) outside the planets' reach. The alternative scenario of the close flyby of another star can instead produce all these TNO features simultaneously, but the possible parameter space for such an encounter is vast. Here, we compare observed TNO properties with thousands of flyby simulations to determine the specific properties of a flyby that reproduces all the different dynamical TNO populations, their location and their relative abundance and find that a 0.8$^{+0.1}_{-0.1}$ $M_{\odot}$ star passing at a distance of $ r_p =$ 110 $\mathbf{\pm}$ 10 au, inclined by $i$ = 70$^\circ$ $^{+5}_{-10}$ gives a near-perfect match. This flyby also replicates the retrograde TNO population, which has proved difficult to explain. Such a flyby is reasonably frequent; at least 140 million solar-type stars in the Milky Way are likely to have experienced a similar one. In light of these results, we predict that the upcoming Vera Rubin telescope will reveal that distant and retrograde TNOs are relatively common.

Gaia/DR3 has provided a large sample of new g-mode pulsators, among which~11,600 are Gam Dor stars. This work present the spectroscopic parameters of these Gam Dor pulsators estimated by the GSP-spec module that analysed millions of Gaia spectra. The Galactic positions, kinematics, and orbital properties of these new Gaia pulsators were examined in order to define a sub-sample belonging to the Milky Way thin disc, in which these young stars should preferentially be found. The stellar luminosities, radii, and astrometric surface gravities were estimated without adopting any priors from uncertain stellar evolution models. These parameters, combined with the GSP-spec effective temperatures, spectroscopic gravities, and metallicities were then validated by comparison with recent literature studies. Most stars are found to belong to the Galactic thin disc, as expected. It is also found that the derived luminosities, radii, and astrometric surface gravities are of high quality and have values typical of genuine Gam Dor pulsators. Moreover, we have shown that Teff and [M/H] of pulsators with high enough SNR spectra or slow to moderate rotation rates are robust. This allowed to define a sub-sample of genuine slow-rotating Gaia Gam Dor pulsators. Their Teff are found between ~6500 and ~7800K, log(g) around 4.2 and luminosities and stellar radii peak at ~5Lsun and ~1.7Rsun. [M/H] is close to the Solar value, although 0.5dex more metal-poor and metal-rich Gam Dor are identified. The [alpha/Fe] content is fully consistent with the chemical properties of the Galactic disc. Gaia/DR3 spectroscopic properties of Gam Dor stars therefore confirm the nature of these pulsators and allow to chemo-physically parametrise a new large sample of such stars. Moreover, future Gaia DR should drastically increase the number of Gam Dor stars with good-precision spectroscopically derived parameters.

In this chapter, we will cover how stars form from the stellar nurseries that are giant molecular clouds. We will first review the physical processes that compete to regulate star formation. We then review star formation in turbulent, magnetized molecular clouds and the associated statistics giving rise to the star formation rate and the initial mass function of stars. We then present the protostellar stages in detail from an observational perspective. We will primarily discuss low-mass ($<1.5\msun$) stars. Finally, we examine how multiplicity complicates the single-star formation picture. This chapter will focus on star formation at redshift~0

With the expanding integration of infrared instruments in astronomical missions, accurate per-pixel flux estimation for near-infrared hybrid detectors has become critical to the success of these missions. Based on CPPM's involvement in both SVOM/Colibri and Euclid missions, this study introduces universally applicable methods and framework for characterizing IR hybrid detectors and decorrelating their intrinsic properties. The characterization framework, applied to the ALFA detector and \Euclid's H2RG, not only validates the proposed methods but also points out subtle behaviors inherent to each detector.

Neutron star (NS) has many extreme physical conditions, and one may obtain some important informations about NS via accreting neutron star binary (ANSB) systems. The upcoming Chinese Space Station Telescope (CSST) provides an opportunity to search for a large sample of ANSB candidates. Our goal is to check the completeness of the potential ANSB samples from CSST data. In this paper, we generate some ANSBs and normal binaries under CSST photometric system by binary evolution and binary population synthesis method and use a machine learning method to train a classification model. Although the Precision ($94.56~ \%$) of our machine learning model is as high as before study, the Recall is only about $63.29~ \%$. The Precision/Recall is mainly determined by the mass transfer rate between the NSs and their companions. In addition, we also find that the completeness of ANSB samples from CSST photometric data by the machine learning method also depends on the companion mass and the age of the system. ANSB candidates with low initial mass companion star ($0.1~ {\rm M}_\odot$ to $1~ {\rm M}_\odot$) have a relatively high Precision ($94.94~ \%$) and high Recall ($86.32~ \%$), whereas ANSB candidates with higher initial mass companion star ($1.1~ {\rm M}_\odot$ to $3~ {\rm M}_\odot$) have similar Precision ($93.88~ \%$) and quite low Recall ($42.67~ \%$). Our results indicate that although the machine learning method may obtain a relative pure sample of ANSBs, a completeness correction is necessary for one to obtain a complete sample.

We study a mid-order wavefront sensor (MOWFS) to address fine cophasing errors in exoplanet imaging with future large segmented aperture space telescopes. Observing Earth analogs around Sun-like stars requires contrasts down to $10^{-10}$ in visible light. One promising solution consists of producing a high-contrast dark zone in the image of an observed star. In a space observatory, this dark region will be altered by several effects, and among them, the small misalignments of the telescope mirror segments due to fine thermo-mechanical drifts. To correct for these errors in real time, we investigate a wavefront control loop based on a MOWFS with a Zernike sensor. Such a MOWFS was installed on the high-contrast imager for complex aperture telescopes (HiCAT) testbed in Baltimore in June 2023. The bench uses a 37-segment Iris-AO deformable mirror to mimic telescope segmentation and some wavefront control strategies to produce a dark zone with such an aperture. In this contribution, we first use the MOWFS to characterize the Iris-AO segment discretization steps. For the central segment, we find a minimal step of 125\,$\pm$31\,pm. This result will help us to assess the contribution of the Iris-AO DM on the contrast in HiCAT. We then determine the detection limits of the MOWFS, estimating wavefront error amplitudes of 119 and 102\,pm for 10\,s and 1\,min exposure time with a SNR of 3. These values inform us about the measurement capabilities of our wavefront sensor on the testbed. These preliminary results will be useful to provide insights on metrology and stability for exo-Earth observations with the Habitable Worlds Observatory.

This work focuses on the study of the AME, an important emission mechanism between 10 and 60 GHz whose polarization properties are not yet fully understood, and is therefore a potential contaminant for future CMB polarization observations. We use new QUIJOTE-MFI maps 11, 13, 17 and 19 GHz, together with other public ancillary data including WMAP and Planck, to study the polarization properties of the AME in three Galactic regions: rho-Ophiuchi, Perseus and W43. We have obtained the SEDs for those three regions over the frequency range 0.4-3000 GHz, both in intensity and polarization. The intensity SEDs are well described by a combination of free-free emission, thermal dust, AME and CMB anisotropies. In polarization, we extracted the flux densities using all available data between 11 and 353 GHz. We implemented an improved intensity-to-polarization leakage correction that has allowed for the first time to derive reliable polarization constraints well below the 1% level from Planck-LFI data. A frequency stacking of maps in the range 10-60 GHz has allowed us to reduce the statistical noise and to push the upper limits on the AME polarization level. We have obtained upper limits on the AME polarization fraction of order <1% (95% confidence level) for the three regions. In particular we get Pi_AME < 1.1% (at 28.4 GHz), Pi_AME < 1.1% (at 22.8 GHz) and Pi_AME < 0.28% (at 33 GHz) in rho-Ophiuchi, Perseus and W43 respectively. At the QUIJOTE 17 GHz frequency band, we get Pi_AME< 5.1% for rho-Ophiuchi, Pi_AME< 3.5% for Perseus and Pi_AME< 0.85% for W43. Our final upper limits derived using the stacking procedure are Pi_AME < 0.58% for rho-Ophiuchi, Pi_AME < 1.64% for Perseus and Pi_AME < 0.31% for W43. Altogether, these are the most stringent constraints to date on the AME polarization fraction of these three star-forming regions.

The aim of our study is to characterise the spin states of the members of the Athor and Zita collisional families and test whether these members have a spin distribution consistent with a common origin from the break up of their respective family parent asteroids. Our method is based on the asteroid family evolution, which indicates that there should be a statistical predominance of retrograde-rotating asteroids on the inward side of family's V-shape, and prograde-rotating asteroids on the outward side. We used photometric data from our campaign and the literature in order to reveal the spin states of the asteroids belonging to these families. We combined dense and sparse photometric data in order to construct lightcurves; we performed the lightcurve inversion method to estimate the sidereal period, spin axis and convex shape of several family members. We obtained 34 new asteroid models for Athor family members and 17 for Zita family members. Along with the literature and revised models, the Athor family contains 60% of retrograde asteroids on the inward side and, 76% of prograde asteroids on the outward side. We also found that the Zita family exhibits 80% of retrograde asteroids on the inward side and an equal amount of prograde and retrograde rotators on the outward side. However, when we applied Kernel density estimation, we also found a clear peak for prograde asteroids on the outward side, as expected from the theory. The spin states of these asteroids validate the existence of both families, with the Athor family exhibiting a stronger signature for the presence of retrograde-rotating and prograde-rotating asteroids on the inner and outer side of the family, respectively. Our work provides an independent confirmation and characterisation of these very old families, whose presence and characteristics offer constraints for theories and models of the Solar System's evolution.

Context. The complex dynamics of asymptotic giant branch (AGB) stars and the resulting stellar winds have a significant impact on the measurements of stellar parameters and amplify their uncertainties. Three-dimensional (3D) radiative hydrodynamic (RHD) simulations of convection suggest that convection-related structures at the surface of AGB star affect the photocentre displacement and the parallax uncertainty measured by Gaia. Aims. We explore the impact of the convection on the photocentre variability and aim to establish analytical laws between the photocentre displacement and stellar parameters to retrieve such parameters from the parallax uncertainty. Methods. We used a selection of 31 RHD simulations with CO5BOLD and the post-processing radiative transfer code Optim3D to compute intensity maps in the Gaia G band [320-1050 nm]. From these maps, we calculated the photocentre position and temporal fluctuations. We then compared the synthetic standard deviation to the parallax uncertainty of a sample of 53 Mira stars observed with Gaia. Results. The simulations show a displacement of the photocentre across the surface ranging from 4 to 13 % of the corresponding stellar radius, in agreement with previous studies. We provide an analytical law relating the pulsation period of the simulations and the photocentre displacement as well as the pulsation period and stellar parameters. By combining these laws, we retrieve the surface gravity, the effective temperature, and the radius for the stars in our sample. Conclusions. Our analysis highlights an original procedure to retrieve stellar parameters by using both state-of-the-art 3D numerical simulations of AGB stellar convection and parallax observations of AGB stars. This will help us refine our understanding of these giants.

Most domains of science are experiencing a paradigm shift due to the advent of a new generation of instruments and detectors which produce data and data streams at an unprecedented rate. The scientific exploitation of these data, namely Data Driven Discovery, requires interoperability, massive and optimal use of Artificial Intelligence methods in all steps of the data acquisition, processing and analysis, the access to large and distributed computing HPC facilities, the implementation and access to large simulations and interdisciplinary skills that usually are not provided by standard academic curricula. Furthermore, to cope with this data deluge, most communities have leveraged solutions and tools originally developed by large corporations for purposes other than scientific research and accepted compromises to adapt them to their specific needs. Through the presentation of several astrophysical use cases, we show how the Data Driven based solutions could represent the optimal playground to achieve the multi-disciplinary methodological approach.

In an attempt to further characterise the near-Earth asteroid (NEA) population we present 38 new light-curves acquired between September 2020 and November 2023 for NEAs (7335) 1989 JA, (7822) 1991 CS, (154244) 2002 KL6 and (159402) 1999 AP10, obtained from observations taken at the Teide Observatory (Tenerife, Spain). With these new observations along with archival data, we computed their first shape models and spin solutions by applying the light curve inversion method. The obtained rotation periods are in good agreement with those reported in previous works, with improved uncertainties. Additionally, besides the constant period models for (7335) 1989 JA, (7822) 1991 CS and (159402) 1999 AP10, our results for (154244) 2002 KL6 suggest that it could be affected by a Yarkovsky-O'Keefe-Radzievskii-Paddack acceleration with a value of $\upsilon \simeq -7 \times 10^{-9}$ rad d$^{-2}$. This would be one of the first detections of this effect slowing down an asteroid.

Context. Tailed radio galaxies are shaped by ram pressure owing to the high-velocity motion of their host through the intracluster medium (ICM). Recent works have reported on the increasing complexity of the phenomenology of tailed galaxies, with departures from theoretical ageing models and evidence of re-energising mechanisms, which are yet unclear. Aims. The nearby (z = 0.0894) galaxy cluster Abell 2142 hosts two tailed galaxies, namely T1 and T2, which exhibit peculiar morphological features. We aim to investigate the properties of T1 and T2 and constrain their spectral evolution, dynamics, and interactions with the ICM. Methods. We combined LOw Frequency Array (LOFAR), upgraded Giant Metrewave Radio Telescope (uGMRT), Very Large Array (VLA), and MeerKAT data (from 30 MHz to 6.5 GHz) to carry out a detailed spectral analysis of T1 and T2. We analysed surface brightness profiles, measured integrated and spatially-resolved spectral indices, and performed a comparison with single injection ageing models. Chandra X-ray data were used to search for discontinuities in the ICM properties in the direction of the targets. Results. The spectral properties of T1 at low frequencies are predicted by ageing models, and provide constraints on the 3D dynamics of the host by assuming a constant velocity. However, sharp transitions along sub-regions of the tail, local surface brightness enhancements, and a spectral shape at high frequencies that is not predicted by models suggest a more complex scenario, possibly involving hydrodynamical instabilities and particle mixing. T2 exhibits unusual morphological and surface brightness features, and its spectral behaviour is not predicted by standard models. Two AGN outburst events during the infall of T2 towards the cluster centre could explain its properties.

To prepare for the analyses of the future PLATO light curves, we develop a deep learning model, Panopticon, to detect transits in high precision photometric light curves. Since PLATO's main objective is the detection of temperate Earth-size planets around solar-type stars, the code is designed to detect individual transit events. The filtering step, required by conventional detection methods, can affect the transit, which could be an issue for long and shallow transits. To protect transit shape and depth, the code is also designed to work on unfiltered light curves. We trained the model on a set of simulated PLATO light curves in which we injected, at pixel level, either planetary, eclipsing binary, or background eclipsing binary signals. We also include a variety of noises in our data, such as granulation, stellar spots or cosmic rays. The approach is able to recover 90% of our test population, including more than 25% of the Earth-analogs, even in the unfiltered light curves. The model also recovers the transits irrespective of the orbital period, and is able to retrieve transits on a unique event basis. These figures are obtained when accepting a false alarm rate of 1%. When keeping the false alarm rate low (<0.01%), it is still able to recover more than 85% of the transit signals. Any transit deeper than 180ppm is essentially guaranteed to be recovered. This method is able to recover transits on a unique event basis, and does so with a low false alarm rate. Thanks to light curves being one-dimensional, model training is fast, on the order of a few hours per model. This speed in training and inference, coupled to the recovery effectiveness and precision of the model make it an ideal tool to complement, or be used ahead of, classical approaches.

Integrated optics are used to achieve astronomical interferometry inside robust and compact materials, improving the instrument's stability and sensitivity. In order to perform differential phase measurements at the H$\alpha$ line (656.3 nm) with the 600-800 nm spectro-interferometer FIRST, a photonic integrated circuit (PIC) is being developed in collaboration with TEEM Photonics. This PIC performs the interferometric combination of the beams coming from sub-apertures selected in the telescope pupil, thus implementing the pupil remapping technique to restore the diffraction limit of the telescope. In this work, we report on the latest developments carried out within the FIRST project to produce a high performance visible PIC. The PICs are manufactured by TEEM Photonics, using their technology based on $K_+:Na_+$ ion exchange in glass. The first part of the study consists in the experimental characterization of the fundamental properties of the waveguides, in order to build an accurate model, which is the basis for the design of more complex functions. In the second part, theoretical designs and their optimization for three types of combiner architectures are presented: symmetric directional coupler, asymmetric directional couplers and ABCD cells including achromatic phase shifters.

Coronal jets are the captivating eruptions which are often found in the solar atmosphere, and primarily formed due to magnetic reconnection. Despite their short-lived nature and lower energy compared to many other eruptive events, e.g. flares and coronal mass ejections, they play an important role in heating the corona and accelerating charged particles. However, their generation in the ambience of non-standard flare regime is not fully understood, and warrant a deeper investigation, in terms of their onset, growth, eruption processes, and thermodynamic evolution. Toward this goal, this paper reports the results of a data-constrained three-dimensional (3D) magnetohydrodynamics (MHD) simulation of an eruptive jet; initialized with a Non-Force-Free-Field (NFFF) extrapolation and carried out in the spirit of Implicit Large Eddy Simulation (ILES). The simulation focuses on the magnetic and dynamical properties of the jet during its onset, and eruption phases, that occurred on February 5, 2015 in an active region NOAA AR12280, associated with a seemingly three-ribbon structure. In order to correlate its thermal evolution with computed energetics, the simulation results are compared with differential emission measurement (DEM) analysis in the vicinity of the jet. Importantly, this combined approach provides an insight to the onset of reconnection in transients in terms of emission and the corresponding electric current profiles from MHD evolutions. The presented study captures the intricate topological dynamics, finds a close correspondence between the magnetic and thermal evolution in and around the jet location. Overall, it enriches the understanding of the thermal evolution due to MHD processes, which is one of the broader aspects to reveal the coronal heating problem.

In this study, we explored the fundamental contrast limit of NIRCam coronagraphy observations, representing the achievable performance with post-processing techniques. This limit is influenced by photon noise and readout noise, with complex noise propagation through post-processing methods like principal component analysis. We employed two approaches: developing a formula based on simplified scenarios and using Markov Chain Monte Carlo (MCMC) methods, assuming Gaussian noise properties and uncorrelated pixel noise. Tested on datasets HIP\,65426, AF\,Lep, and HD\,114174, the MCMC method provided accurate but computationally intensive estimates. The analytical approach offered quick, reliable estimates closely matching MCMC results in simpler scenarios. Our findings showed the fundamental contrast curve is significantly deeper than the current achievable contrast limit obtained with post-processing techniques at shorter separations, being 10 times deeper at $0.5''$ and 4 times deeper at $1''$. At greater separations, increased exposure time improves sensitivity, with the transition between photon and readout noise dominance occurring between $2''$ and $3''$. We conclude the analytical approach is a reliable estimate of the fundamental contrast limit, offering a faster alternative to MCMC. These results emphasize the potential for greater sensitivity at shorter separations, highlighting the need for improved or new post-processing techniques to enhance JWST NIRCam sensitivity or contrast curve.

The irregular moons orbit the giant planets on distant, inclined, and eccentric trajectories, in sharp contrast with the coplanar and quasicircular orbits of the regular moons. The origin of these irregular moons is still an open question, but these moons have a lot in common with the objects beyond Neptune (trans-Neptunian objects -- TNOs), suggestive of a common origin. Here, we show that the close flyby of a star may be the connecting element. A stellar flyby can simultaneously reproduce the complex TNO dynamics quantitatively while explaining the origin of the irregular moons and the colour distributions of both populations. This flyby would have catapulted 7.2% of the original TNO population into the region of the planets, many on retrograde orbits. Most injected TNOs would have been subsequently ejected from the solar system (85%). However, a considerable fraction would have had the potential to be captured by the planets. The exclusively distant origin of the injected TNOs may also explain the lack of very red irregular moons.

SN 2023ixf is one of the most well-observed core-collapse supernova in recent decades, yet there is inconsistency in the inferred zero-age-main-sequence (ZAMS) mass $M_{\rm ZAMS}$ of its progenitor. Direct observations of the pre-SN red supergiant (RSG) estimate $M_{\rm ZAMS}$ spanning widely from 11 to 18 $M_{\rm \odot}$. Additional constraints, including host environment and the pulsation of its progenitor RSG, suggest a massive progenitor with $M_{\rm ZAMS}$ > 17 $M_{\rm \odot}$. However, the analysis of the properties of supernova, from light curve modeling to late phase spectroscopy, favor a relatively low mass scenario ($M_{\rm ZAMS}$ < 15 $M_{\rm \odot}$). In this work, we conduct systematic analysis of SN 2023ixf, from the RSG progenitor, plateau phase light curve to late phase spectroscopy. Using MESA+STELLA to simulate the RSG progenitor and their explosions, we find that, despite the ZAMS mass of the RSG models being varied from 12.0 to 17.5 $M_{\rm \odot}$, they can produce light curves that well match with SN 2023ixf if the envelope mass and the explosion energy are allowed to vary. Using late phase spectroscopy as independent measurement, the oxygen emission line [O I] suggests the ZAMS is intermediate massive (~16.0 $M_{\rm \odot}$), and the relatively weak H$\alpha$ emission line indicates the hydrogen envelope has been partially removed before the explosion. By incorporating the velocity structure derived from the light curve modeling into an axisymmetric model, we successfully generated [O I] line profiles that are consistent with the [O I] line observed in late phase spectroscopy of SN 2023ixf. Bringing these analyses together, we conclude that SN 2023ixf is the aspherical explosion of an intermediate massive star ($M_{\rm ZAMS}$ = 15-16 $M_{\rm \odot}$) with the hydrogen envelope being partially stripped to 4-5 $M_{\rm \odot}$ prior to its explosion.

We use polarization data from SOFIA HAWC+ to investigate the interplay between magnetic fields and stellar feedback in altering gas dynamics within the high-mass star-forming region RCW 36, located in Vela C. This region is of particular interest as it has a bipolar HII region powered by a massive star cluster which may be impacting the surrounding magnetic field. To determine if this is the case, we apply the Histogram of Relative Orientations (HRO) method to quantify the relative alignment between the inferred magnetic field and elongated structures observed in several datasets such as dust emission, column density, temperature, and spectral line intensity maps. The HRO results indicate a bimodal alignment trend, where structures observed with dense gas tracers show a statistically significant preference for perpendicular alignment relative to the magnetic field, while structures probed by photo-dissociation region (PDR) tracers tend to align preferentially parallel relative to the magnetic field. Moreover, the dense gas and PDR associated structures are found to be kinematically distinct such that a bimodal alignment trend is also observed as a function of line-of-sight velocity. This suggests that the magnetic field may have been dynamically important and set a preferred direction of gas flow at the time that RCW 36 formed, resulting in a dense ridge developing perpendicular to the magnetic field. However on filament-scales near the PDR region, feedback may be energetically dominating the magnetic field, warping its geometry and the associated flux-frozen gas structures, causing the observed the preference for parallel relative alignment.

Gamma-ray burst (GRB) afterglows are powerful probes for studying the different properties of their host galaxies (e.g., the interstellar dust) at all redshifts. By fitting their spectral energy distribution (SED) over a large range of wavelengths, we can gain direct insights into the properties of the interstellar dust by studying the extinction curves. Unlike the dust extinction templates, such as the average Milky Way (MW) or the Small and Large Magellanic Cloud (SMC and LMC), the extinction curves of galaxies outside the Local Group exhibit deviation from these laws. Altogether, X-ray and gamma-ray satellites as well as ground-based telescopes, such as Neil Gehrels Swift Observatory (Swift) and Gamma-Ray Optical and Near-Infrared Detector (GROND), provide measurements of the afterglows from the X-ray to the NIR, which can be used to extract information on dust extinction curves along their lines of sight. The study presented in this paper undertakes such a photometric study, comprising a preparatory work for the SVOM mission and its ground-based follow-up telescope COLIBRI. We propose a simple parameterization of the dust extinction curve of GRB host galaxies. Our model is based on a power law form with the addition of a Loretzian-like Drude profile with two parameters: the extinction slope, $\gamma$, and the 2175 $Å$ bump amplitude, $E_\rm{b}$. Using the seven GROND filter bands, we tested our dust extinction model and explored the parameter space in extinction and redshift by fitting SEDs of simplified simulations of GRB afterglow spectra based on different extinction curve templates. From a final sample of 10 real Swift/GROND extinguished GRBs, we determined the quantities of the dust extinction in their host and measured their extinction curves. The measured average extinction curve is equivalent to a quasi-featureless in-between SMC-LMC template.

Magnetic field evolution of neutron stars is a long-standing debate. The rate of magnetic field decay for isolated, non-accreting neutron stars can be quantified by measuring the negative second derivative of the spin period. Alternatively, this rate can be estimated by observing an excess of thermal emission with respect to the standard cooling without additional heating mechanisms involved. One of the nearby cooling isolated neutron stars -- RX J0720.4-3125, -- offers a unique opportunity to probe the field decay as for this source there are independent measurements of the surface X-ray luminosity, the second spin period derivative, and magnetic field. We demonstrate that the evolution rate of the spin period derivative is in correspondence with the rate of dissipation of magnetic energy of the dipolar field if a significant part of the released energy is emitted in X-rays. The instantaneous time scale for the magnetic field decay is $\sim 10^4$ years.

At high frequencies beyond the acoustic cut-off, a peak-like structure is visible in the solar power spectrum. Known as the pseudo-modes, their frequencies have been shown to vary in anti-phase with solar magnetic activity. In this work, we determined temporal variations in these frequencies across the solar disc, with the aim of identifying any potential latitudinal dependence of pseudo-mode frequency shifts. We utilised nearly 22 years of spatially resolved GONG data for all azimuthal orders, $\textit{m}$, for harmonic degrees 0 $\leq$ $\textit{l}$ $\leq$ 200, and determined shifts using the resampled periodogram method. Periodogram realisations were created from overlapping, successive 216d-long segments in time, and cropped to 5600-6800$\mu$Hz. Cross-correlation functions were then repeatedly generated between these realisations to identify any variation in frequency and the uncertainty. We categorised each mode by its latitudinal sensitivity and used this categorisation to produce average frequency shifts for different latitude bands (15$^\circ$ and 5$^\circ$ in size) which were compared to magnetic proxies, the $F_{\mathrm{10.7}}$ index and GONG synoptic maps. Morphological differences in the pseudo-mode shifts between different latitudes were found, which were most pronounced during the rise to solar maximum where shifts reach their minimum values. At all latitudes, shift behaviour was strongly in anti-correlation with the activity proxy. Additionally, periodicities shorter than the 11-year cycle were observed. Wavelet analysis was used to identify a periodicity of four years at all latitudes.

Understanding the global rotational profile of the solar atmosphere and its variation is fundamental to uncovering a comprehensive understanding of the dynamics of the solar magnetic field and the extent of coupling between different layers of the Sun. In this study, we employ the method of image correlation to analyze the extensive dataset provided by the Atmospheric Imaging Assembly of the Solar Dynamic Observatory in different wavelength channels. We find a significant increase in the equatorial rotational rate ($A$) and a decrease in absolute latitudinal gradient ($|B|$) at all temperatures representative of the solar atmosphere, implying an equatorial rotation up to $4.18\%$ and $1.92\%$ faster and less differential when compared to the rotation rates for the underlying photosphere derived from Doppler measurement and sunspots respectively. In addition, we also find a significant increase in equatorial rotation rate ($A$) and a decrease in differential nature ($|B|$ decreases) at different layers of the solar atmosphere. We also explore a possible connection from the solar interior to the atmosphere and interestingly found that $A$ at $r=0.94\,\mathrm{R}_{\odot}, 0.965\,\mathrm{R}_{\odot}$ show an excellent match with 171 Angstrom, 304 Angstrom and 1600 Angstrom, respectively. Furthermore, we observe a positive correlation between the rotational parameters measured from 1600 Angstrom, 131 Angstrom, 193 Angstrom and 211 Angstrom with the yearly averaged sunspot number, suggesting a potential dependence of the solar rotation on the appearance of magnetic structures related to the solar cycle or the presence of cycle dependence of solar rotation in the solar atmosphere.

According to the well established hierarchical framework for galaxy evolution, galaxies grow through mergers with other galaxies and the LambdaCDM cosmological model predicts that the stellar halos of galaxies are rich in remnants from minor mergers. The Stellar Streams Legacy Survey has provided a first release of a catalogue with a statistically significant sample of stellar streams in the Local Universe that can be used to study minor mergers and test the cosmological models. The main objective is to compare the results of the observations of stellar tidal streams with the predictions of state-of-the-art cosmological simulations regarding the formation of stellar streams up to a redshift z < 0.02, according to the LambdaCDM model. We use the predictions of the cosmological simulations Copernicus Complexio, TNG50 of the IllustrisTNG project and Auriga to generate 225 mock-images of nearby halos at a distance of 70 Mpc, and search for stellar streams. We compare the obtained stream frequency and characteristics with those obtained from the Stellar Streams Legacy Survey. We find good agreement between the results of analysing real images from the Dark Energy Survey and mock-images from cosmological simulations. We obtained predictions for the detection rate of stellar streams to a surface brightness limit of 35 mag arcsec^-2. The cosmological simulations predict that for a surface brightness limit of 32 mag arcsec^-2 a frequency of almost 70% in the detection of streams around galaxies can be achieved.

PSR J1618-3921 is one of five known millisecond pulsars (MSPs) in eccentric orbits (eMPSs) located in the Galactic plane, whose formation is poorly understood. Earlier studies of these objects revealed significant discrepancies between observation and predictions from standard binary evolution scenarios of pulsar-Helium white dwarf binaries. We conducted observations with the L-band receiver of the MeerKAT radio telescope and the UWL receiver of the Parkes Murriyang radio telescope between 2019 and 2021. These data were added to archival observations. We perform an analysis of this joint 23-year-dataset. We use the recent observations to give a brief account of the emission properties of J1618-3921, including a Rotating Vector model fit of the linear polarisation position angle of the pulsar. The long timing baseline allowed for a highly significant measurement of the rate of advance of periastron of $\dot{\omega}$. We can only report a low significance detection of the orthometric Shapiro delay parameters $h_3$ and $\varsigma$, leading to mass estimates of the total and individual binary masses. We detect an unexpected change in the orbital period of, which is an order of magnitude larger and carries an opposite sign to what is expected from Galactic acceleration and the Shklovskii effect. We also detect a significant second derivative of the spin frequency. Furthermore, we report an unexpected, abrupt change of the mean pulse profile in June 2021 with unknown origin. We propose that the anomalous $\dot{P_b}$ and $\ddot{f}$ indicate an additional varying acceleration due to a nearby mass, i.e., the J1618-3921 binary system is likely part of a hierarchical triple. This finding suggests that at least some eMSPs might have formed in triple star systems. Although the uncertainties are large, the binary companion mass is consistent with the $P_b$ - $M_{WD}$ relation.

We reanalyse the second data release of the European Pulsar Timing Array (EPTA) using an observationally-driven model for ensemble properties of pulsar noise. We show that the revised gravitational wave background properties are in better agreement with theoretical expectations for the strain spectrum. Our improved model for ensemble pulsar noise properties reduces a systematic error at $1\sigma$ level and increases Bayesian odds of Hellings-Downs correlations by $\sim 10\%$.

We describe a novel splitting approach to numerical relativistic magnetohydrodynamics (RMHD) designed to expand its applicability to the domain of ultra-high magnetisation (high-$\sigma$). In this approach, the electromagnetic field is split into the force-free component and its perturbation due to the plasma inertia. Accordingly, the system of RMHD equations is extended to include the subsystem of force-free degenerate electrodynamics and the subsystem governing the plasma dynamics and the perturbation of the force-free field. The combined system of conservation laws is integrated simultaneously, to which aim various numerical techniques can be used, and the force-free field is recombined with its perturbation at the end of every timestep. To explore the potential of this splitting approach, we combined it with a 3rd-order WENO method, and carried out a variety of 1D and 2D test simulations. The simulations confirm the robustness of the splitting method in the high-$\sigma$ regime, and also show that it remains accurate in the low-$\sigma$ regime, all the way down to $\sigma$ = 0. Thus, the method can be used for simulating complex astrophysical flows involving a wide range of physical parameters. The numerical resistivity of the code obeys a simple ansatz and allows fast magnetic reconnection in the plasmoid-dominated regime. The results of simulations involving thin and long current sheets agree very well with the theory of resistive magnetic reconnection.

Since Parker Solar Probe's (Parker's) first perihelion pass at the Sun, large amplitude Alfvén waves grouped in patches have been observed near the Sun throughout the mission. Several formation processes for these magnetic switchback patches have been suggested with no definitive consensus. To provide insight to their formation, we examine the heavy ion properties of several adjacent magnetic switchback patches around Parker's 11th perihelion pass capitalizing on a spacecraft lineup with Solar Orbiter where each samples the same solar wind streams over a large range of longitudes. Heavy ion properties (Fe/O, C$^{6+}$/C$^{5+}$, O$^{7+}$/O$^{6+}$) related to the wind's coronal origin, measured with Solar Orbiter can be linked to switchback patch structures identified near the Sun with Parker. We find that switchback patches do not contain distinctive ion and elemental compositional signatures different than the surrounding non-switchback solar wind. Both the patches and ambient wind exhibit a range of fast and slow wind qualities, indicating coronal sources with open and closed field lines in close proximity. These observations and modeling indicate switchback patches form in coronal hole boundary wind and with a range of source region magnetic and thermal properties. Furthermore, the heavy ion signatures suggest interchange reconnection and/or shear driven processes may play a role in their creation.

The null hypothesis in Pulsar Timing Array (PTA) analyses includes assumptions about ensemble properties of pulsar time-correlated noise. These properties are encoded in prior probabilities for the amplitude and the spectral index of the power-law power spectral density of temporal correlations of the noise. In this work, we introduce a new procedure for numerical marginalisation over the uncertainties in pulsar noise priors. The procedure may be used in searches for nanohertz gravitational waves and other PTA analyses to resolve prior misspecification at negligible computational cost. Furthermore, we infer the distribution of amplitudes and spectral indices of the power spectral density of spin noise and dispersion measure variation noise based on the observation of 25 millisecond pulsars by the European Pulsar Timing Array (EPTA). Our results may be used for the simulation of realistic noise in PTAs.

There is an intriguing and growing population of Neptune-sized planets with stellar obliquities near $\sim90^{\circ}$. One previously proposed formation pathway is a disk-driven resonance, which can take place at the end stages of planet formation in a system containing an inner Neptune, outer cold Jupiter, and protoplanetary disk. This mechanism occurs within the first $\sim10$ Myr, but most of the polar Neptunes we see today are $\sim$Gyrs old. Up until now, there has not been an extensive analysis of whether the polar orbits are stable over $\sim$Gyr timescales. Tidal realignment mechanisms are known to operate in other systems, and if they are active here, this would cause theoretical tension with a primordial misalignment story. In this paper, we explore the effects of tidal evolution on the disk-driven resonance theory. We use both $N$-body and secular simulations to study tidal effects on both the initial resonant encounter and long-term evolution. We find that the polar orbits are remarkably stable on $\sim$Gyr timescales. Inclination damping does not occur for the polar cases, although we do identify sub-polar cases where it is important. We consider two case study polar Neptunes, WASP-107 b and HAT-P-11 b, and study them in the context of this theory, finding consistency with present-day properties if their tidal quality factors are $Q \gtrsim 10^4$ and $Q \gtrsim 10^5$, respectively.

The recent JWST detections of carbon-bearing molecules in a habitable-zone sub-Neptune have opened a new era in the study of low-mass exoplanets. The sub-Neptune regime spans a wide diversity of planetary interiors and atmospheres not witnessed in the solar system, including mini-Neptunes, super-Earths, and water worlds. Recent works have investigated the possibility of gas dwarfs, with rocky interiors and thick H$_2$-rich atmospheres, to explain aspects of the sub-Neptune population, including the radius valley. Interactions between the H$_2$-rich envelope and a potential magma ocean may lead to observable atmospheric signatures. We report a coupled interior-atmosphere modelling framework for gas dwarfs to investigate the plausibility of magma oceans on such planets and their observable diagnostics. We find that the surface-atmosphere interactions and atmospheric composition are sensitive to a wide range of parameters, including the atmospheric and internal structure, mineral composition, volatile solubility and atmospheric chemistry. While magma oceans are typically associated with high-temperature rocky planets, we assess if such conditions may be admissible and observable for temperate sub-Neptunes. We find that a holistic modelling approach is required for this purpose and to avoid unphysical model solutions. We find using our model framework and considering the habitable-zone sub-Neptune K2-18 b as a case study that its observed atmospheric composition is incompatible with a magma ocean scenario. We identify key atmospheric molecular and elemental diagnostics, including the abundances of CO$_2$, CO, NH$_3$ and, potentially, S-bearing species. Our study also underscores the need for fundamental material properties for accurate modelling of such planets.

We decomposed the molecular gas in the spiral galaxy NGC 628 (M74) into multi-scale hub-filament structures using the CO (2-1) line by the dendrogram algorithm. All leaf structures as potential hubs were classified into three categories, i.e. leaf-HFs-A, leaf-HFs-B and leaf-HFs-C. leaf-HFs-A exhibit the best hub-filament morphology, which also have the highest density contrast, the largest mass and the lowest virial ratio. We employed the FILFINDER algorithm to identify and characterize filaments within 185 leaf-HFs-A structures, and fitted the velocity gradients around the intensity peaks. Measurements of velocity gradients provide evidence for gas inflow within these structures. The numbers of the associated 21 $\mu$m and H$_{\alpha}$ structures and the peak intensities of 7.7 $\mu$m, 21 $\mu$m and H$_{\alpha}$ emissions decrease from leaf-HFs-A to leaf-HFs-C. The spatial separations between the intensity peaks of CO and 21 $\mu$m structures of leaf-HFs-A are larger than those of leaf-HFs-C. These evidence indicate that leaf-HFs-A are more evolved than leaf-HFs-C. There may be an evolutionary sequence from leaf-HFs-C to leaf-HFs-A. Currently, leaf-HFs-C lack a distinct gravitational collapse process that would result in a significant density contrast. The density contrast can effectively measure the extent of the gravitational collapse and the depth of the gravitational potential of the structure which, in turn, shapes the hub-filament morphology. Combined with the kinematic analysis presented in previous studies, a picture emerges that molecular gas in spiral galaxies is organized into network structures through the gravitational coupling of multi-scale hub-filament structures. Molecular clouds, acting as knots within these networks, serve as hubs, which are local gravitational centers and the main sites of star formation.

Here we report the detection of dust depletion in a misaligned inner disk around UX Tau A using JWST MIRI spectra. Mid-infrared (MIR) continuum "seesaw" variability was detected in this disk by Spitzer and attributed to variable shadows cast on the outer disk by the inner disk. The JWST MIRI spectrum of UX Tau A also shows seesaw variability, but with a significant decrease of emission shortwards of 10 micron to nearly photospheric levels. We argue that UX Tau A's MIR continuum variability is due to depletion of dust in a misaligned inner disk. We speculate that this dust depletion occurs because the inner disk is misaligned from the outer disk, which can disrupt the replenishment of the inner disk from the outer disk. Using contemporaneous measurements of the mass accretion rate of UX Tau A and estimating the amount of dust necessary to produce the MIR excess in the Spitzer observations, we estimate a minimum dust depletion timescale of ~0.1 yr. These observations show that we can indirectly detect the signatures of misaligned inner disks through MIR continuum variability and that in some cases the inner disk may be significantly depleted of dust and become optically thin.

The sub-Jovian desert is a region in the mass-period and radius-period parameter space, typically encompassing short-period ranges between super-Earths and hot Jupiters, that exhibits an intrinsic dearth of planets. This scarcity is likely shaped by photoevaporation caused by the stellar irradiation received by giant planets that have migrated inward. We report the detection and characterization of TOI-3568 b, a transiting super-Neptune with a mass of $26.4\pm1.0$ M$_\oplus$, a radius of $5.30\pm0.27$ R$_\oplus$, a bulk density of $0.98\pm0.15$ g cm$^{-3}$, and an orbital period of 4.417965(5) d situated in the vicinity of the sub-Jovian desert. This planet orbiting a K dwarf star with solar metallicity, was identified photometrically by TESS. It was characterized as a planet by our high-precision radial velocity monitoring program using MAROON-X at Gemini North, supplemented by additional observations from the SPICE large program with SPIRou at CFHT. We performed a Bayesian MCMC joint analysis of the TESS and ground-based photometry, MAROON-X and SPIRou radial velocities, to measure the orbit, radius, and mass of the planet, as well as a detailed analysis of the high-resolution flux and polarimetric spectra to determine the physical parameters and elemental abundances of the host star. Our results reveal TOI-3568 b as a hot super-Neptune, rich in hydrogen and helium with a core of heavier elements with a mass between 10 and 25 M$_\oplus$. We analyzed the photoevaporation status of TOI-3568 b and found that it experiences one of the highest EUV luminosities among planets with a mass M$_{\rm p}$ $<2$ M$_{\rm Nep}$, yet it has an evaporation lifetime exceeding 5 Gyr. Positioned in the transition between two significant populations of exoplanets on the mass-period and energy diagrams, this planet presents an opportunity to test theories concerning the origin of the sub-Jovian desert.

As part of our comprehensive, ongoing characterisation of the low-mass end of the main sequence in the Solar neighbourhood, we used the OSIRIS instrument at the 10.4 m Gran Telescopio Canarias to acquire low- and mid-resolution (R$\approx$300 and R$\approx$2500) optical spectroscopy of 53 late-M and L ultracool dwarfs. Most of these objects are known but poorly investigated and lacking complete kinematics. We measured spectral indices, determined spectral types (six of which are new) and inferred effective temperature and surface gravity from BT-Settl synthetic spectra fits for all objects. We were able to measure radial velocities via line centre fitting and cross correlation for 46 objects, 29 of which lacked previous radial velocity measurements. Using these radial velocities in combination with the latest Gaia DR3 data, we also calculated Galactocentric space velocities. From their kinematics, we identified two candidates outside of the thin disc and four in young stellar kinematic groups. Two further ultracool dwarfs are apparently young field objects: 2MASSW J1246467$+$402715 (L4$\beta$), which has a potential, weak lithium absorption line, and G 196$-$3B (L3$\beta$), which was already known as young due to its well-studied primary companion.

In the standard cosmological framework, neutrinos begin to free-stream after the weak interaction phase ends in the early universe, at a temperature of approximately $T \sim 1 \, {\rm MeV}$. However, the onset of neutrino free-streaming can be delayed if additional interactions occur in the early universe, leaving imprints on both the cosmic microwave background (CMB) angular power spectra and the large-scale structure (LSS) matter power spectra. We present a thorough analysis of early universe neutrino interactions with a fairly generalized parameterization of the interaction rates as a power law in neutrino temperature. In this (6+2) parameter scenario, we constrain the cosmological parameters along with the free-streaming redshift and the sum of the neutrino mass in presence of such interactions, with the help of full shape (FS) galaxy power spectra from BOSS Data Release 12. Our analysis reveals that a combined dataset of FS along with CMB and BAO offers improved constraints on the free-streaming redshift from present data, comparable to the forecast results from future CMB-S4 data. Additionally, we investigate the prospects of future galaxy surveys by forecasting on Euclid mission in combination with Planck and CMB-S4, and find significant improvement on both the free-streaming redshift and the sum of the neutrino mass than the existing constraints as well as than CMB-S4 alone.

Ultralight scalar fields can experience runaway `superradiant' amplification near spinning black holes, resulting in a macroscopic `axion cloud' which slowly dissipates via continuous monochromatic gravitational waves. For a particular range of boson masses, $\mathcal{O}(10^{-11}$ -- $10^{-10})$ eV, an axion cloud will radiate in the $10$ -- $100$ kHz band of the Levitated Sensor Detector (LSD). Using fiducial models of the mass, spin, and age distributions of stellar-origin black holes, we simulate the present-day Milky Way population of these hypothetical objects. As a first step towards assessing the LSD's sensitivity to the resultant ensemble of GW signals, we compute the corresponding signal-to-noise ratios which build up over a nominal integration time of $10^{7}$ s, assuming the projected sensitivity of the $1$-m LSD prototype currently under construction, as well as for future $10$-m and $100$-m concepts. For a $100$-m cryogenic instrument, hundreds of resolvable signals could be expected if the boson mass $\mu$ is around $3\times10^{-11}$ eV, and this number diminishes with increasing $\mu$ up to $\approx 5.5\times10^{-11}$ eV. The much larger population of unresolved sources will produce a confusion foreground which could be detectable by a $10$-m instrument if $\mu \in (3-4.5)\times10^{-11}$ eV, or by a $100$-m instrument if $\mu \in (3-6)\times10^{-11}$ eV.

Barred galaxies exhibit boxy/peanut or X-shapes (BP/X) protruding from their disks in edge-on views. Two types of BP/X morphologies exist depending on whether the X-wings meet at the center (CX) or are off-centered (OX). Orbital studies indicate that various orbital types can generate X-shaped structures. Here, we provide a classification approach that identifies the specific orbit families responsible for generating OX and CX-shaped structures. Applying this approach to three different N-body bar models, we show that both OX and CX structures are associated with the x1 orbit family, but OX-supporting orbits possess higher angular momentum (closer to x1 orbits) than orbits in CX structures. Consequently, as the bar slows down the contribution of higher angular momentum OX-supporting orbits decreases and that of lower angular momentum orbits increases resulting in an evolution of the morphology from OX to CX. If the bar does not slow down, the shape of the BP/X structure and the fractions of OX/CX supporting orbits remain substantially unchanged. Bars that do not undergo buckling but that do slow down initially show the OX structure and are dominated by high angular momentum orbits, transitioning to a CX morphology. Bars that buckle exhibit a combination of both OX and CX supporting orbits immediately after the buckling, but become more CX dominated as their pattern speed decreases. This study demonstrates that the evolution of BP/X morphology and orbit populations strongly depends on the evolution of the bar angular momentum.