Papers for Wednesday, Oct 30 2024

For modeling the spectra of exoplanets one must know their atmospheric composition. This is necessary because the abundance of molecules, atoms, ions and condensates is needed to construct the total cross-section for the interaction between electro-magnetic radiation and matter. In addition, when solving for the temperature structure of an atmosphere the so-called adiabatic temperature gradient must be known, which describes the pressure-temperature dependence in convectively unstable regions well. Depending on the planetary properties, the composition and adiabatic gradients may be well described by equilibrium chemistry, which means that chemical reactions occur faster than any other relevant processes in the atmosphere, such as mixing. What is more, the equilibrium assumption often serves as a useful starting point for non-equilibrium calculations. Efficient and easy-to-use codes for determining equilibrium abundances are therefore needed. Here we report on our easyCHEM Python package that calculates atmospheric compositions and adiabatic temperature gradients in chemical equilibrium for any user-specified elemental composition.

The stellar halo of the Milky Way comprises an abundance of chemical signatures from accretion events and \textit{in-situ} evolution, that form an interweaving tapestry in kinematic space. To untangle this, we consider the mixtures of chemical information, in a given region of integral of motion space, as a variant of the blind source separation problem and utilise non-negative matrix factorisation (NMF). Specifically, we examine the variation in [Fe/H], [Mg/Fe], and [Al/Fe] distributions of APOGEE DR17 stars across the $(E,L_z)$ plane of the halo. When 2 components are prescribed, the NMF algorithm splits stellar halo into low- and high-energy components in the $(E,L_z)$ plane which approximately correspond to the accreted and \textit{in-situ} halo respectively. We use these two components to define a new boundary between the \textit{in-situ} and the accreted stellar halo. Moreover, we calculate the components fractional contribution to the stellar halo as a function of energy, galactocentric spherical radius, height, and galactocentric cylindrical radius. Using a stellar halo defined by kinematic cuts, we find that the halo transitions from \textit{in-situ} dominated to accretion dominated at $E \approx -1.67 \times 10^5$ (km/s)$^2$ (using the potential in McMillan 2017), and at $(r,z,R) \approx (8.7, 3.0, 8.1)$ kpc. The low-energy component is found to span a range of [Al/Fe] that falls beyond the typically accepted \textit{in-situ} floor of [Al/Fe] $=0$. Upon prescribing more components to the NMF model, we find hints of the existence of overlapping chemical evolution sequences that other techniques struggle to find. We also examine features within these components that resemble known substructures in the halo, such as \textit{Eos} and \textit{Aurora}. This work provides insight into their origin and the part they play in the Milky Way's formation.

Active galactic nuclei (AGN) feedback from supermassive black holes (SMBHs) at the centers of galaxy clusters plays a key role in regulating star formation and shaping the intracluster medium (ICM), often manifesting through prominent X-ray cavities embedded in the cluster's hot atmosphere. Here we show that X-ray cavities arise naturally due to AGN feedback in TNG-Cluster. This is a new suite of magnetohydrodynamic cosmological simulations of galaxy formation and evolution, and hence of galaxy clusters, whereby cold dark matter, baryon dynamics, galactic astrophysics, and magnetic fields are evolved together consistently. We construct mock Chandra X-ray observations of the central regions of the 352 simulated clusters at $z=0$ and find that $\sim$39 per cent contain X-ray cavities. Identified X-ray cavities vary in configuration (single, pairs, or multiples), with some still attached to SMBHs, while others have buoyantly risen. Their size ranges from a few to several tens of kpc. In terms of gas physical properties, TNG-Cluster X-ray cavities are underdense compared to the surrounding halo and filled with hot gas ($\sim$10$^8$K); 25 per cent of them are surrounded by an X-ray bright and compressed rim associated with a weak shock (Mach number $\sim 1.5$). Clusters exhibiting X-ray cavities are preferentially strong or weak cool-cores, are dynamically relaxed, and host SMBHs accreting at low Eddington rates. We show that TNG-Cluster X-ray cavities originate from episodic, wind-like energy injections from central AGN. Our results illustrate the existence and diversity of X-ray cavities simulated in state-of-the-art models within realistic cosmological environments and show that these can form without necessarily invoking bipolar, collimated, or relativistic jets.

Descriptions of the Galactic Center using Fermi gamma-ray data have so far modeled the Galactic Center Excess (GCE) as a template with fixed spatial morphology or as a linear combination of such templates. Although these templates are informed by various physical expectations, the morphology of the excess is a priori unknown. For the first time, we describe the GCE using a flexible, non-parametric machine learning model -- the Gaussian process (GP). We assess our model's performance on synthetic data, demonstrating that the model can recover the templates used to generate the data. We then fit the \Fermi data with our model in a single energy bin from 2-20 GeV (leaving a spectral GP analysis of the GCE for future work) using a variety of template models of diffuse gamma-ray emission to quantify our fits' systematic uncertainties associated with diffuse emission modeling. We interpret our best-fit GP in terms of GCE templates consisting of an NFW squared template and a bulge component to determine which bulge models can best describe the fitted GP and to what extent the best-fit GP is described better by an NFW squared template versus a bulge template. The best-fit GP contains morphological features that are typically not associated with traditional GCE studies. These include a localized bright source at around $(\ell,b) = (20^{\circ}, 0^{\circ})$ and a diagonal arm extending Northwest from the Galactic Center. In spite of these novel features, the fitted GP is explained best by a template-based model consisting of the bulge presented in Coleman et al. (2020) and a squared NFW component. Our results suggest that the physical interpretation of the GCE in terms of stellar bulge and NFW-like components is highly sensitive to the assumed morphologies, background models, and the region of the sky used for inference.

Despite the ubiquity of clumpy star-forming galaxies at high-redshift, the origin of clumps are still largely unconstrained due to the limited observations that can validate the mechanisms for clump formation. We postulate that if clumps form due to the accretion of metal-poor gas that leads to violent disk instability, clumpy galaxies should have lower gas-phase metallicities compared to non-clumpy galaxies. In this work, we obtain the near-infrared spectrum for 42 clumpy and non-clumpy star-forming galaxies of similar masses, SFRs, and colors at $z\approx0.7$ using the Gemini Near-Infrared Spectrograph (GNIRS) and infer their gas-phase metallicity from the {\nii} and {\halpha} line ratio. We find that clumpy galaxies have lower metallicities compared to non-clumpy galaxies, with an offset in the weighted average metallicity of $0.07\pm0.02$ dex. We also find an offset of $0.06\pm0.02$ dex between clumpy and non-clumpy galaxies in a comparable sample of 23 star-forming galaxies at $z\approx1.5$ using existing data from the FMOS-COSMOS survey. Similarly, lower {\nii}/{\halpha} ratio are typically found in galaxies that have more of their $\mathrm{UV_{rest}}$ luminosity originating from clumps, suggesting that \enquote{clumpier} galaxies are more metal poor. We also derive the intrinsic velocity dispersion and line-of-sight rotational velocity for galaxies from the GNIRS sample. The majority of galaxies have $\sigma_0/v_c \approx 0.2$, with no significant difference between clumpy and non-clumpy galaxies. Our result indicates that clump formation may be related to the inflow of metal-poor gas; however, the process that forms them does not necessarily require significant, long-term kinematic instability in the disk.

We consider approximating the linearly evolved 2-point correlation function (2pcf) of dark matter $\xi_{\rm lin}(r;\boldsymbol{\theta})$ in a cosmological model with parameters $\boldsymbol{\theta}$ as the linear combination $\xi_{\rm lin}(r;\boldsymbol{\theta})\approx\sum_i\,b_i(r)\,w_i(\boldsymbol{\theta})$, where the functions $\mathcal{B}=\{b_i(r)\}$ form a $\textit{model-agnostic basis}$ for the linear 2pcf. This decomposition is important for model-agnostic analyses of the baryon acoustic oscillation (BAO) feature in the nonlinear 2pcf of galaxies that fix $\mathcal{B}$ and leave the coefficients $\{w_i\}$ free. To date, such analyses have made simple but sub-optimal choices for $\mathcal{B}$, such as monomials. We develop a machine learning framework for systematically discovering a $\textit{minimal}$ basis $\mathcal{B}$ that describes $\xi_{\rm lin}(r)$ near the BAO feature in a wide class of cosmological models. We use a custom architecture, denoted $\texttt{BiSequential}$, for a neural network (NN) that explicitly realizes the separation between $r$ and $\boldsymbol{\theta}$ above. The optimal NN trained on data in which only $\{\Omega_{\rm m},h\}$ are varied in a $\textit{flat}$ $\Lambda$CDM model produces a basis $\mathcal{B}$ comprising $9$ functions capable of describing $\xi_{\rm lin}(r)$ to $\sim0.6\%$ accuracy in $\textit{curved}$ $w$CDM models varying 7 parameters within $\sim5\%$ of their fiducial, flat $\Lambda$CDM values. Scales such as the peak, linear point and zero-crossing of $\xi_{\rm lin}(r)$ are also recovered with very high accuracy. We compare our approach to other compression schemes in the literature, and speculate that $\mathcal{B}$ may also encompass $\xi_{\rm lin}(r)$ in modified gravity models near our fiducial $\Lambda$CDM model. Using our basis functions in model-agnostic BAO analyses can potentially lead to significant statistical gains.

We present new stellar population models, $\alpha$-MC, self-consistently taking into account non-solar $\rm [\alpha/Fe]$ abundances for both isochrones and stellar spectra. The $\alpha$-MC models are based on $\alpha$-enhanced MIST isochrones and C3K spectral libraries, which are publicly available in FSPS. Our new models cover a wide range of ages ($\rm \log (age/yr) = 5.0 - 10.3$), metallicities ($\rm [Fe/H]=[-2.5,+0.5]$ in steps of 0.25, $\rm [\alpha/Fe]=-0.2,+0.0,+0.2,+0.4,+0.6$), and wavelengths ($0.1-2.5\,\rm \mu m$). We investigate the separate and combined effects of $\alpha$-enhanced isochrones and stellar spectral libraries on simple stellar populations (SSPs), including their broadband colors, spectral indices, and full spectra. We find that the primary effect of $\alpha$-enhancement in isochrones is to lower the overall continuum levels and redden the continuum shapes, while $\alpha$-enhancement in stellar spectra mainly affects individual spectral lines. At constant $\rm [Fe/H]$, $\alpha$-enhancement has significant impacts on the broadband colors by $\rm \sim 0.1-0.4\,mag$ across all ages ($\rm 0.01 - 10\,Gyr$). The effects of $\alpha$-enhancement on colors at fixed $\rm [Z/H]$ are smaller, by $\rm \sim 0.1-0.2\,mag$. The spectral indices involving $\alpha$-elements, Ca4227 and Mg b, increase with $\rm [\alpha/Fe]$ (both at fixed $\rm [Fe/H]$ and fixed $\rm [Z/H]$) due to enhanced $\alpha$-abundances. At constant $\rm [Fe/H]$, $\alpha$-enhancement weakens most Fe-sensitive and Hydrogen Balmer lines. Our new self-consistent $\alpha$-enhanced models will be essential in deriving accurate physical properties of high-redshift galaxies, where $\alpha$-enhancement is expected to be common.

Within the Milky Way (MW), younger stellar populations exhibit steeper (more negative) metallicity radial gradients; the origin of this trend remains debated. The FIRE-2 cosmological simulations of MW-mass galaxies show the same trend as the MW, which in FIRE-2 arises because the metallicity gradient of the interstellar medium (ISM), and thus of stars at birth, became steeper over time. We seek to understand this evolution in the context of inside-out radial growth of galaxies. Most FIRE-2 galaxies grew radially inside-out in both gas and stars; specifically, their surface density profiles, $\Sigma(R)$, became shallower over time. Combined with a realized superlinear (Kennicutt-Schmidt-like) relation between star formation and total gas density, the profile of the ratio $\Sigma_{\rm star}(R)/\Sigma_{\rm gas}(R)$ became shallower (flatter) over time. Thus, if metals stayed where they were injected into the ISM from stars, the metallicity gradient would become shallower over time, as some models predict. However, metallicity gradients in FIRE-2 became steeper over time, because of the additional effects of (radial) mixing of metals in the ISM. Specifically, the velocity dispersion and net radial advection of gas declined over time, as ISM turbulence decreased and the disk settled, leading to upside-down vertical growth. In FIRE-2, this evolution in metal mixing of gas associated with upside-down growth dominates over inside-out radial growth, causing the metallicity radial gradient of the ISM and of stars at birth to become steeper over time. We argue that this reflects the ISM history of the MW and of typical MW-mass galaxies.

We first briefly review the adventure of scale invariance in physics, from Galileo Galilei, Weyl, Einstein, and Feynman to the revival by Dirac (1973) and Canuto et al. (1977). We then gather concrete observational evidence that scale-invariant effects are present and measurable in astronomical objects spanning a vast range of masses (0.5 M$_{\odot} <$ M $< 10^{14}$ M$_{\odot}$) and an equally impressive range of spatial scales (0.01 pc $<$ r $<$ 1 Gpc). Scale invariance accounts for the observed excess in velocity in galaxy clusters with respect to the visible mass, the relatively flat/small slope of rotation curves in local galaxies, the observed steep rotation curves of high-redshift galaxies, and the excess of velocity in wide binary stars with separations above 3000 kau found in Gaia DR3. Last but not least, we investigate the effect of scale invariance on gravitational lensing. We show that scale invariance does not affect the geodesics of light rays as they pass in the vicinity of a massive galaxy. However, scale-invariant effects do change the inferred mass-to-light ratio of lens galaxies as compared to GR. As a result, the discrepancies seen in GR between the total lensing mass of galaxies and their stellar mass from photometry may be accounted for. This holds true both for lenses at high redshift like JWST-ER1 and at low redshift like in the SLACS sample. Of note is that none of the above observational tests require dark matter or any adjustable parameter to tweak the theory at any given mass or spatial scale.

We present optical and IR observations from maximum light until around 600 d of SN 2022jli, a peculiar SE SN showing two maxima, each one with a peak luminosity of about 3 x 10^{42} erg/s and separated by 50 d. The second maximum is followed by periodic undulations with a period of P ~ 12.5 days. The spectra and the photometric evolution of the first maximum are consistent with the behaviour of a standard SE SN with an ejecta mass of 1.5 +/- 0.4 Msun, and a nickel mass of 0.12 +/- 0.01 Msun. The optical spectra after 400 d correspond to a standard SN Ic event, and at late times SN 2022jli exhibits a significant drop in the optical luminosity implying that the physical phenomena that produced the secondary maximum has ceased to power the SN light curve. One possibility is that the second maximum is powered by a magnetar with an initial spin period of P=48.5 ms and a magnetic field of B = 8.5x10^{14} G, while the light curve periodic undulations could be produced by accretion of material from a companion star onto the neutron star in a binary system. The near-IR spectra shows clear 1st CO overtone emission from about 190 d after the first maximum, and it becomes undetected at 400 d. A significant near-IR excess from hot dust emission is detected at 238 d produced by either newly formed dust in the SN ejecta or due to a strong near-IR dust echo. Depending on the assumptions of the dust composition, the estimated dust mass is 2-16 x 10^{-4} Msun. The magnetar power of the second maximum can fit in a more general picture where magnetars are the power source of super-luminous SNe, could produce their frequent bumps and undulations, and where pulsars could produce the late time excess observed in some SE SNe. The detection of CO and the potential detection of dust formed in the ejecta of SN2022jli are important to understand the formation molecules and dust in the ejecta of SE SNe.

Revealing the internal composition and structure of giant planets is fundamental for understanding planetary formation. However, the bulk composition can only be inferred through interior models. As a result, advancements in modelling aspects are essential to better characterise the interiors of giant planets. We investigate the effects of model assumptions such as the interior structure and the hydrogen-helium (H-He) equation of state (EOS) on the inferred interiors of giant exoplanets. We first assess these effects on a few test cases and compare H-He EOSs. We then calculate evolution models and infer the planetary bulk metallicity of 45 warm exoplanets, ranging from 0.1 to 10~$M_{\rm J}$. Planets with masses between about 0.2 and 0.6~$M_{\rm J}$ are most sensitive to the H-He EOS. Updating the H-He EOS reduces the inferred heavy-element mass, with an absolute difference in bulk metallicity of up to 13\%. Concentrating heavy elements in a core, rather than distributing them uniformly (and scaling opacities with metallicity), reduces the inferred metallicity (up to 17\%). The assumed internal structure, along with its effect on the envelope opacity, has the greatest effect on the inferred composition of massive planets ($M_{\rm p}>4~M_{\rm J}$). For $M_{\rm p}>0.6~M_{\rm J}$, the observational uncertainties on radii and ages lead to uncertainties in the inferred metallicity (up to 31\%) which are larger than the ones associated with the used H-He EOS and the assumed interior structure. However, for planets with $0.2<M_{\rm p}<0.6~M_{\rm J}$, the theoretical uncertainties are larger. Advancements in equations of state and our understanding of giant planet interior structures combined with accurate measurements of the planetary radius and age are crucial for characterising giant exoplanets.

The goal of this project is to construct an estimator for the masses of supermassive black holes in active galactic nuclei (AGNs) based on the broad Halpha emission line. We make use of published reverberation mapping data. We remeasure all Halpha time lags from the original data as we find that often the reverberation measurements are improved by detrending the light curves. We produce mass estimators that require only the Halpha luminosity and the width of the Halpha emission line as characterized by either the FWHM or the line dispersion. It is possible, on the basis of a single spectrum covering the Halpha emission line, to estimate the mass of the central supermassive black hole in AGNs, taking into account all three parameters believed to affect mass measurement: luminosity, line width, and Eddington ratio. The typical formal accuracy in such estimates is of order 0.2-0.3 dex.

The interstellar medium of the Milky Way contains turbulent plasma with structures driven by energetic processes that fuel star formation and shape the evolution of our Galaxy. Radio waves from pulsars are scattered off the small (au-scale and below) structures, resulting in frequency-dependent interference patterns that are modulated in time because of the relative motions of the pulsar, Earth, and plasma. Power spectral analyses of these patterns show parabolic arcs with curvatures that encode the locations and kinematics of individual structures. Here we report the discovery of at least 25 distinct plasma structures in the direction of the brilliant millisecond pulsar, PSR J0437$-$4715, in observations obtained with the MeerKAT radio telescope. Four arcs reveal structures within 5000 au of the pulsar, from a series of shocks induced as the pulsar and its wind interact with the ambient insterstellar medium. The measured radial distance and velocity of the main shock allows us to solve the shock geometry and space velocity of the pulsar in three dimensions, while the velocity of another structure unexpectedly indicates a back flow from the direction of the shock or pulsar-wind tail. The remaining 21 arcs represent a surprising abundance of structures sustained by turbulence within the Local Bubble -- a region of the interstellar medium thought to be depleted of gas by a series of supernova explosions about 14 Myr ago.

Melnick 34 (Mk 34) is one of the most massive binary systems known and is one of the brightest X-ray point sources in the 30 Doradus region. We investigated the impact of this massive system on the surrounding interstellar medium (ISM) using the optical spectroscopic capabilities of the narrow-field mode (NFM) of the Multi-Unit Spectroscopic Explorer (MUSE). MUSE-NFM spatially resolved the ISM in the vicinity of Mk 34 with a resolution comparable to that of the HST. The analysis of the [NII]$\lambda$6583 and [SII]$\lambda$6717 emission lines reveals a cone-like structure apparently originating from Mk 34 and extending southeast. Electron density maps and radial velocity measurements of the ISM lines further support an outflow scenario traced by these emissions. While no clear northwestern counterpart to this outflow was observed, we note increased extinction in that direction, towards the R136 cluster. The ISM material along the projected diagonal of the outflow on both sides of Mk 34 shows similar properties in terms of the emission line ratios seen in the Baldwin-Phillips-Terlevich diagram. These results are consistent across two observational epochs. Additionally, we examined the residual maps within a 0.5" radius of Mk 34 after modeling and subtracting the point spread function. The observed variations in the residuals could potentially be linked to Mk 34's known periodic behavior. However, further observations with appropriate cadence are needed to fully monitor the 155 day periodicity of Mk 34's X-ray emissions.

Among massive stars, binary interaction is the rule rather than the exception. The closest binaries, those with periods of less than about 10 days, undergo mass transfer during core-hydrogen burning, with many of them experiencing a nuclear-timescale contact phase. Current binary population synthesis models predict the mass-ratio distribution of contact binaries to be heavily skewed toward a mass ratio of unity, which is inconsistent with observations. It has been shown that effects of tidal deformation due to the Roche potential, as well as energy transfer in the common layers of a contact binary, alter the internal structure of close binary components. However, previous population studies neglected these effects. We model a population of massive binary stars that undergo mass transfer during core-hydrogen burning, while consistently considering the effects of tidal deformation and energy transfer in contact phases. We use the MESA binary-evolution code to compute large grids of models with primary star masses of $8\,M_\odot$ to $70\,M_\odot$ at Solar metallicity. We then perform a population synthesis study to predict distribution functions of the observational properties of close binary systems, focusing in particular on the mass and luminosity ratio distribution. We find that the effects of tidal deformation and energy transfer have a limited effect on the predicted mass-ratio distribution of massive contact binaries. Only a small fraction of the population has their mass ratio significantly shifted toward a more unequal configuration. However, we suggest that orbital hardening could affect the evolution of contact binaries and their progenitors, and we advocate for a homogeneous set of observed contact binary parameters.

Prime-Cam, one of the primary instruments for the Fred Young Submillimeter Telescope (FYST) developed by the CCAT Collaboration, will house up to seven instrument modules, with the first operating at 280 GHz. Each module will include three arrays of superconducting microwave kinetic inductance detectors (KIDs). The first KID array fabricated for the 280 GHz module uses titanium-nitride (TiN) as the superconducting material and has 3,456 individual detectors, while the other two arrays use aluminum. This paper presents the design and laboratory characterization of the 280 GHz TiN array, which is cooled below its critical temperature to ~0.1 K and read out over six RF feedlines. LED mapping, a technique for matching the measured resonant frequency of a detector to its physical position, was performed on the array so that the results can be used to lithographically trim the KID capacitors and increase the yield of the array by reducing frequency collisions. We present the methods and results of LED mapping the 280 GHz TiN KID array before deployment on FYST.

Over the past decade, several millimeter interferometer programs have mapped the nearby star-forming galaxy M51 at a spatial resolution of ${\le}170$ pc. This study combines observations from three major programs: the PdBI Arcsecond Whirlpool Survey (PAWS), the SMA M51 large program (SMA-PAWS), and the Surveying the Whirlpool at Arcseconds with NOEMA (SWAN). The dataset includes the (1-0) and (2-1) rotational transitions of $^{12}$CO, $^{13}$CO, and C$^{18}$O isotopologues. The observations cover the $r{<}\rm 3\,kpc$ region including center and part of the disk, thereby ensuring strong detections of the weaker $^{13}$CO and C$^{18}$O lines. All observations are convolved in this analysis to an angular resolution of 4$''$, corresponding to a physical scale of ${\sim}$170 pc. We investigate empirical line ratio relations and quantitatively evaluate molecular gas conditions such as temperature, density, and the CO-to-H$_2$ conversion factor ($\alpha_{\rm CO}$). We employ two approaches to study the molecular gas conditions: (i) assuming local thermal equilibrium (LTE) to analytically determine the CO column density and $\alpha_{\rm CO}$, and (ii) using non-LTE modeling with RADEX to fit physical conditions to observed CO isotopologue intensities. We find that the $\alpha_{\rm CO}$ values {in the center and along the inner spiral arm} are $\sim$0.5 dex (LTE) and ${\sim}$0.1 dex (non-LTE) below the Milky Way inner disk value. The average non-LTE $\alpha_{\rm CO}$ is $2.4{\pm}0.5$ M$_\odot$ pc$^{-2}$ (K km s$^{-1}$)$^{-1}$. While both methods show dispersion due to underlying assumptions, the scatter is larger for LTE-derived values. This study underscores the necessity for robust CO line modeling to accurately constrain the molecular ISM's physical and chemical conditions in nearby galaxies.

Correctly modelling the absorptive properties of dust and haze particles is of great importance for determining the abundance of solid matter within protoplanetary disks and planetary atmospheres. Rigorous analyses such as the discrete dipole approximation (DDA) can be used to obtain accurate absorption cross-sections, but these require significant computing time and are often impractical to use in models. A simple analytical equation exists for spherical particles in the long-wavelength limit (where the wavelength is much larger than the size of the dust particle), but we demonstrate that this can significantly underestimate the absorption. This effect is found to depend strongly on refractive index, with values of m = 1 + 11i corresponding to an underestimate in absorption by a factor of 1,000. Here we present MANTA-Ray (Modified Absorption of Non-spherical Tiny Aggregates in the RAYleigh regime): a simple model that can calculate absorption efficiencies within 10-20% of the values predicted by DDA, but 10^(13) times faster. MANTA-Ray is very versatile and works for any wavelength and particle size in the long wavelength regime. It is also very flexible with regards to particle shape, and can correctly model structures ranging from long linear chains to tight compact clusters, composed of any material with refractive index 1+0.01i < m < 11+11i. The packaged model is provided as publicly-available code for use by the astrophysical community.

As the number of gravitational-wave detections grows, the merger rate of binary black holes (BBHs) can help us to constrain their formation, the properties of their progenitors, and their birth environment. Here, we aim to address the impact of the metal-dependent star formation rate (SFR) on the BBH merger rate. To this end, we have developed a fully data-driven approach to model the metal-dependent SFR and coupled it to BBH evolution. We have adopted the most up-to-date scaling relations, based on recent observational results, and we have studied how the BBH merger rate density varies over a wide grid of galaxy and binary evolution parameters. Our results show that including a realistic metal-dependent SFR evolution yields a value of the merger rate density which is too high compared to the one inferred from GW data. Moreover, variations of the SFR in low-mass galaxies ($M_\ast \lesssim 10^8 \mathrm{M}_{\odot}$) do not contribute more than a factor $\sim 2$ to the overall merger rate density at redshift $z=0$. These results suggest that the discrepancy between the BBH merger rate density inferred from data and theoretical models is not caused by approximations in the treatment of the metal-dependent SFR, but rather stems from stellar evolution models and/or BBH formation channels.

The observed prevalence of galaxies exhibiting bursty star formation histories (SFHs) at $z\gtrsim6$ has created new challenges and opportunities for understanding their formation pathways. The degenerate effects of the efficiency and burstiness of star formation on the observed UV luminosity function are separable by galaxy clustering. However, quantifying the timescales of burstiness requires more than just the continuum UV measurements. Here we develop a flexible semi-analytic framework for modeling both the amplitude of star formation rate (SFR) variations and their temporal correlation, from which the luminosity function and clustering can be derived for SFR indicators tracing different characteristic timescales (e.g., UV continuum and H$\alpha$ luminosities). Based on this framework, we study the prospect of using galaxy summary statistics to distinguish models where SFR fluctuations are prescribed by different power spectral density (PSD) forms. Using the Fisher matrix approach, we forecast the constraints on parameters in our PSD-based model that can be extracted from mock JWST observations of the UV and H$\alpha$ luminosity functions and clustering bias factors at $z\sim6$. These constraints demonstrate the feasibility of constraining the burstiness of high-$z$ galaxies solely from their one-point and two-point statistics and underscore the importance of combining tracers of both long-term and short-term SFR variations. Our flexible framework can be readily extended to characterize the SFH of high-redshift galaxies with a wider range of observational diagnostics.

The timeline of cosmic reionization remains uncertain despite sustained efforts to study how the ionizing output of early galaxies shaped the intergalactic medium (IGM). Using the semi-numerical code LIMFAST, we investigate the prospects for timing the reionization process by cross-correlating the 21 cm signal with the cosmic near-infrared background (NIRB) contributed by galaxies at $z>5$. Tracing opposite phases of the IGM on large scales during reionization, the two signals together serve as a powerful probe for the reionization history. However, because long-wavelength, line-of-sight Fourier modes -- the only modes probed by NIRB fluctuations -- are contaminated by 21 cm foregrounds and thus inevitably lost to foreground cleaning or avoidance, a direct cross-correlation of the two signals vanishes. We show that this problem can be circumvented by squaring the foreground-filtered 21 cm signal and cross-correlating the squared field with the NIRB. This statistic is related to the 21 cm--21 cm--NIRB cross-bispectrum and encodes valuable information regarding the reionization timeline. Particularly, the 21 cm$^2$ and NIRB signals are positively correlated during the early phases of reionization and negatively correlated at later stages. We demonstrate that this behavior is generic across several different reionization models and compare our simulated results with perturbative calculations. We show that this cross-correlation can be detected at high significance by forthcoming 21 cm and NIRB surveys such as SKA and SPHEREx. Our methodology is more broadly applicable to cross-correlations between line intensity mapping data and 2D tracers of the large-scale structure, including photometric galaxy surveys and CMB lensing mass maps, among others.

The presence of dips in the gravito-inertial modes period-spacing pattern of gamma-Dor stars is now well established by recent asteroseismic studies. Such Lorentzian-shaped inertial dips arise from the interaction of gravito-inertial modes propagating in the radiative envelope of intermediate-mass main sequence stars with pure inertial modes that propagate in their convective core. We aim to investigate the signature of a differential rotation between the convective core and the near-core region inside gamma-Dor stars from the inertial dip properties. We first describe the bi-layer rotation profile we use and the approximations we adopt to maintain the analyticity of our study. We then describe our results on the inertial dip formation, location, and shape. We derive a modified Lorentzian profile and we compare it to the previously obtained results in the solid-body rotation case. This work highlights the inertial dips' probing power of the convective core rotation, an important observable in the context of the understanding of the angular momentum transport and chemicals mixing inside stars.

Planetary systems exhibiting mean-motion resonances (MMRs) offer unique opportunities to study the imprint of disk-induced migration on the orbital architectures of planetary systems. The HD 45364 system, discovered via the radial velocity (RV) method to host two giant planets in a 3:2 MMR, has been the subject of several studies attempting to reconstruct the system's orbital migration history based on its present-day resonant configuration. Recently, Li et al. (2022) called into question the system's residence in the 3:2 MMR based on a revised orbital solution derived from an expanded set of RV observations that extend the time baseline of the original discovery data by over a decade. However, we show that inferences about the planets' dynamical state with respect to the 3:2 MMR are sensitive to the particular prior assumptions adopted in the orbital modeling. Using $N$-body dynamical models, we show that orbital solutions constrained to reside deep in the 3:2 MMR fit the RV data with a similar quality to unconstrained orbital solutions. We conclude that the RV observations of HD 45364 are consistent with orbital configurations produced by smooth migration and resonance capture. We further show that past convergent orbital migration can reproduce the system's present-day orbital configuration provided that the ratio of migration to eccentricity damping timescales, $K$, was in the range $10\lesssim K \lesssim 175$. We also find that dynamical interactions in the system can break the usual mass-inclination degeneracy inherent to Keplerian models of RV observations and constrain the planets' absolute masses to within a factor of $\sim2$.

We study the importance of several processes that influence the evolution of dust and its grain size distribution on spatially resolved scales in nearby galaxies. Here, we compiled several multi-wavelength observations for the nearby galaxies NGC628(M74), NGC5457(M101), NGC598(M33), and NGC300. We applied spatially resolved spectral energy distribution fitting to the latest iteration of infrared data to get constraints on the galaxy dust masses and the small-to-large grain abundance ratio. For comparison, we took the radial profiles of the stellar mass and gas mass surface density for NGC628 combined with its metallicity gradient in the literature to calibrate a single-galaxy simulation using the GADGET4-OSAKA code. The simulations include a parametrization to separate the dense and diffuse phases of the ISM where different dust-evolution mechanisms are in action. We find that our simulation can reproduce the radial profile of dust mass surface density but overestimates the SLR in NGC628. Changing the dust-accretion timescale has little impact on the dust mass or SLR, as most of the available metals are accreted onto dust grains at early times (< 3Gyr), except in the outer regions of the galaxy. This suggests we can only constrain the accretion timescale of galaxies at extremely low metallicities where accretion still competes with other mechanisms controlling the dust budget. The overestimation of the SLR likely results from (i) overly efficient shattering processes in the diffuse interstellar medium, which were calibrated to reproduce Milky Way-type galaxies and/or (ii) our use of a diffuse and dense gas density subgrid model that does not entirely capture the intricacies of the small-scale structure present in NGC628.

The Seyfert 1 AGN Fairall 9 was targeted by NICER, Swift, and ground-based observatories for a $\sim$1000-day long reverberation mapping campaign. The following analysis of NICER spectra taken at a two-day cadence provides new insights into the structure and heating mechanisms of the central black hole environment. Observations of Fairall 9 with NICER and Swift revealed a strong relationship between the flux of the UV continuum and the X-ray soft excess, indicating the presence of a "warm" Comptonized corona which likely lies in the upper layers of the innermost accretion flow, serving as a second reprocessor between the "hot" X-ray corona and the accretion disk. The X-ray emission from the hot corona lacks sufficient energy and variability to power slow changes in the UV light curve on timescales of 30 days or longer, suggesting an intrinsic disk-driven variability process in the UV and soft X-rays. Fast variability in the UV on timescales shorter than 30 days can be explained through X-ray reprocessing, and the observed weak X-ray/UV correlation suggests that the corona changes dynamically throughout the campaign.

We present conditions for which X-ray spectra can be ``unfolded'' to present accurate representation of the true source spectra. The method we use to unfold the data is implemented in the \textit{Interactive Spectral Interpretation Software} \citep{Houck2000} and distinguishes itself as being model-independent. We find that this method of unfolding makes accurate representations of the true source spectra (1) The detector is high-resolution and (2) The spectrum is not steeply sloped. These criteria are not simple conditions that give concrete determinations; each detector and spectrum must be judged individually. We find that both grating and imaging detectors can be unfolded with minimal distortions as compared to both continuum and local spectral features; the latter CCD detectors being much more energy dependent. We also provide example use cases for unfolding in the context of current generation X-ray observatories and important caveats.

Strong stellar winds are an important feature in wind-accreting high-mass X-ray binary (HMXB) systems, providing insights into stellar evolution and their impact on surrounding environments. However, the long-term evolution and temporal variability of these winds are not fully understood. This work probes the archetypal wind-accreting HMXB Vela X-1 using MAXI observations over 14 years, focusing on orbit-to-orbit absorption variability in the 2-10 keV band. Additionally, the relation between hardness ratio trends in binary orbits and neutron star spin states is investigated. We calculate hardness ratios to track absorption variability, comparing flux changes across energy bands, as the effect of absorption on the flux is energy-dependent. Variability is analyzed by comparing hardness ratio trends across binary orbits to the MAXI long-term averaged evolution. The long-term averaged hardness ratio evolution displays a stable pattern. Yet, individual binary orbits reveal different hardness ratio evolutions between consecutive orbits with no evident periodicity. Less than half of the binary orbits align with the long-term evolution. Moreover, neutron star spin-up episodes exhibit harder-than-average hardness trends compared to spin-down episodes, although their distributions overlap considerably. The long-term averaged hardness ratio dispersion is consistent with absorption column densities reported in literature from shorter observations, suggesting that heterogeneous wind structures, including accretion wakes and wind clumps, drive observed variations. The orbit-to-orbit variability indicates that pointed X-ray observations provide limited insight into wind structure. The link between neutron star spin states and hardness trends underscores the influence of accretion on absorption, with variability tied to stellar wind density fluctuations.

We investigate the effect of primordial non-Gaussianities on halo number counts using N-body simulations with different values of $f_{\rm NL}^{\rm loc}$. We show how current theoretical models fail to adequately describe the non-Gaussian mass function of halos identified with different overdensity thresholds, $\Delta_{\rm b}$. We explain how these discrepancies are related to a variation in the density profile of dark matter halos, finding that the internal steepness (i.e. the compactness) of halos depends on the value of $f_{\rm NL}^{\rm loc}$. We then parametrize these deviations in halo number counts with a factor $\kappa(\Delta_{\rm b})$ that modifies the linear density threshold for collapse according to the halo identification threshold used, defined with respect to the Universe background density. We rely on a second-degree polynomial to describe $\kappa$ and employ a Bayesian analysis to determine the coefficients of this polynomial. In addition, we verify the independence of the latter on the sign and absolute value of $f_{\rm NL}^{\rm loc}$. Finally, we show how this re-parametrization prevents the extraction of biased constraints on $f_{\rm NL}^{\rm loc}$, correcting for large systematic errors especially in the case of halos identified with high density thresholds. This improvement is crucial in the perspective of deriving cosmological constraints with the non-Gaussian mass function from real data, as different mass definitions can be employed depending on the properties of the survey.

Early dark energy, an additional component of dark energy active in the decade of redshift before recombination, has emerged as one of the most effective models at reducing the $H_0$ tension between direct measurement of the Hubble parameter $H_0$ in the late-universe and the $\Lambda$CDM prediction when calibrated on Planck. However, it requires a slight increase in the dark matter density $\omega_{\rm cdm}$ and primordial tilt $n_s$ that worsens the $S_8$ tension between measurements of weak gravitational lensing at low redshifts and the Planck/$\Lambda$CDM prediction. Using a phenomenological fluid model, we investigate whether the inclusion of a drag term between dark matter and early dark energy can compensate for the effect of the increase in power at small-scales, such that both $H_0$ and $S_8$ tensions are simultaneously alleviated. We find that this works if the drag term is dynamically relevant in the post-recombination universe. However, a drag term active before or just around the time at which the early dark energy contribution to the energy density is maximum is significantly constrained due to its impact on the matter perturbations before recombination, and the subsequent modifications to the cosmic microwave background power spectra.

We propose a novel dust battery mechanism for generating seed magnetic fields in the early universe, in which charged dust grains are radiatively accelerated, inducing strong electric currents that subsequently generate magnetic fields. Our analysis demonstrates that this process is effective even at very low metallicities (approximately $ \sim 10^{-5} Z_\odot$), and capable of producing seed fields with significant amplitudes of $B \sim \rm \mu G$ around luminous sources over timescales of years to Myr and across spatial scales ranging from au to kpc. Crucially, we find that this mechanism is generically $\sim10^8$ times more effective than the radiatively-driven electron battery or Biermann battery in relatively cool gas ($\ll 10^{5}\,$K), including both neutral and ionized gas. Furthermore, our results suggest that, to first order, dissipation effects do not appear to significantly impede this process, and that it can feasibly generate coherent seed fields on macroscopically large ISM scales (much larger than turbulent dissipation scales or electron mean-free-paths in the ISM). These seed fields could then be amplified by subsequent dynamo actions to the observed magnetic fields in galaxies. Additionally, we propose a sub-grid model for integration into cosmological simulations, and the required electric-field expressions for magnetohydrodynamic-particle-in-a-cell (MHD-PIC) simulations that explicitly model dust dynamics. Finally, we explore the broad applicability of this mechanism across different scales and conditions, emphasizing its robustness compared to other known battery mechanisms.

Inferring atmospheric properties of exoplanets from observed spectra is key to understanding their formation, evolution, and habitability. Since traditional Bayesian approaches to atmospheric retrieval (e.g., nested sampling) are computationally expensive, a growing number of machine learning (ML) methods such as neural posterior estimation (NPE) have been proposed. We seek to make ML-based atmospheric retrieval (1) more reliable and accurate with verified results, and (2) more flexible with respect to the underlying neural networks and the choice of the assumed noise models. First, we adopt flow matching posterior estimation (FMPE) as a new ML approach to atmospheric retrieval. FMPE maintains many advantages of NPE, but provides greater architectural flexibility and scalability. Second, we use importance sampling (IS) to verify and correct ML results, and to compute an estimate of the Bayesian evidence. Third, we condition our ML models on the assumed noise level of a spectrum (i.e., error bars), thus making them adaptable to different noise models. Both our noise level-conditional FMPE and NPE models perform on par with nested sampling across a range of noise levels when tested on simulated data. FMPE trains about 3 times faster than NPE and yields higher IS efficiencies. IS successfully corrects inaccurate ML results, identifies model failures via low efficiencies, and provides accurate estimates of the Bayesian evidence. FMPE is a powerful alternative to NPE for fast, amortized, and parallelizable atmospheric retrieval. IS can verify results, thus helping to build confidence in ML-based approaches, while also facilitating model comparison via the evidence ratio. Noise level conditioning allows design studies for future instruments to be scaled up, for example, in terms of the range of signal-to-noise ratios.

WASP-12 b is an ultra-hot Jupiter (UHJ) of special interest for atmospheric studies since it is on an inspiraling orbit in an extreme environment of intense radiation and circumstellar gas. Previously claimed detections of active mass loss from this planet are controversial across the literature. To address this controversy, we obtain two new transit observations of WASP-12 b with the optical high-resolution PEPSI spectrograph on the Large Binocular Telescope. Contrary to previous work, we do not observe planetary H$\alpha$ absorption and rule out the amplitude of previously reported detections. Our non-detection may be limited by the sensitivity of our data or could indicate weaker mass loss than suggested by previous studies. We conduct injection-recovery experiments to place constraints on the radial extent of WASP-12 b's escaping atmosphere as probed by Balmer lines, but find that our data do not have the sensitivity to probe down to the planet's Roche Lobe. Using physically motivated models of atmospheric escape, we explore upper limit constraints on the planet's mass-loss rate and deem the data quality in the wavelength regime of Balmer lines insufficient to determine a physically meaningful constraint. We also conduct a spectral survey of other optical absorbers to trace atmospheric circulation but detect no additional absorption. We conclude that previous claims of H$\alpha$ absorption from the atmosphere of WASP-12 b should be reevaluated. Given the anticipated line strength of Balmer/optical features, observing the atmosphere of this faint target will require stacking more observations even with the largest telescope facilities available.

Testing possible variations in fundamental constants of nature is a crucial endeavor in observational cosmology. This paper investigates potential cosmological variations in the fine structure constant ($\alpha$) through a non-parametric approach, using galaxy cluster observations as the primary cosmological probe. We employ two methodologies based on galaxy cluster gas mass fraction measurements derived from X-ray and Sunyaev-Zeldovich observations, along with luminosity distances from type Ia supernovae. We also explore how different values of the Hubble constant ($H_0$) impact the variation of $\alpha$ across cosmic history. When using the Planck satellite's $H_0$ observations, a constant $\alpha$ is ruled out at approximately the 3$\sigma$ confidence level for $z \lesssim 0.5$. Conversely, employing local estimates of $H_0$ restores agreement with a constant $\alpha$.

We present a new analytical model for the attenuation to Epoch of Reionization (EoR) galaxies by proximate neutral hydrogen gas. Many galaxy spectra in the EoR taken by JWST have shown a flux deficit at wavelengths just redward of the Lyman break, and this has been regarded as resulting from Ly$\alpha$ damping wing absorption by the increasing amount of neutral hydrogen in the line-of-sight. However, previous attenuation models for the intergalactic medium (IGM) commonly used in photometric redshift template-fitting codes assume that the Lyman break is rather sharp, which leads to systematic overestimation of photometric redshifts at $z>7$. In this letter, we build and empirically calibrate a new attenuation model that takes the increased Ly$\alpha$ damping wing absorption into account. Our model consists of the canonical IGM absorption and an additional absorption component due to dense neutral hydrogen gas clouds proximate to the galaxy, and we derive the redshift evolution of HI column density of the proximate clouds by calibrating the model using CANUCS JWST observations. The resulting total transmission curve resolves the photometric redshift bias at $z>7$, an improvement that is robust to choice of template-fitting code, template set, and photometric catalog used. The new attenuation model can be easily implemented in existing template-fitting codes, and significantly improves the photometric redshift performance in the EoR.

Regular remote sensing of the magnetic field embedded within the million-degree solar corona is severely lacking. This reality impedes fundamental investigations of the nature of coronal heating, the generation of solar and stellar winds, and the impulsive release of energy into the solar system via flares and other eruptive phenomena. Resulting from advancements in large aperture solar coronagraphy, we report unprecedented maps of polarized spectra emitted at 1074 nm by Fe+12 atoms in the active corona. We detect clear signatures of the Zeeman effect that are produced by the coronal magnetic field along the optically thin path length of its formation. Our comparisons with global magnetohydrodynamic models highlight the valuable constraints that these measurements provide for coronal modeling efforts, which are anticipated to yield subsequent benefits for space weather research and forecasting.

Context. With the increasing number of discoveries of globular clusters in the inner Milky Way, the need for spectroscopic confirmation and further investigation of their stellar populations and chemodynamical properties has become crucial. Aims. Gran 5 is a newly reported low-mass globular cluster located close to the Galactic center, and it is thought to be an accreted object associated with the Gaia-Enceladus structure. This study aims to investigate the stellar populations of Gran 5 and their detailed chemical properties. Methods. We performed high-resolution near-infrared spectroscopy on seven stars in the field of Gran 5 using IGRINS on the Gemini-South telescope. Results. We identified six stars as cluster members and reveal that they are divided into two stellar populations with different metallicities, with mean [Fe/H] values of -0.76 dex and -0.55 dex, respectively. In addition, the chemodynamical properties of Gran 5 agree with those of in situ globular clusters. Conclusions. Our findings represent the first detection of two stellar populations with different metallicities in a low-mass globular cluster. This suggests that the metallicity variation in Gran 5 may have arisen from processes different from those in other globular clusters with metallicity variation, or that it may have lost a substantial amount of its initial mass during its evolution.

We examined the double red clump (RC) observed in the Galactic bulge, interpreted as a difference in distance ("X-shaped bulge scenario") or in chemical composition ("multiple population scenario"). To verify chemical differences between the RC groups, we performed low-resolution spectroscopy for RC and red giant branch (RGB) stars using Gemini-South/GMOS in three fields of the bulge, and collected diverse data from literature. We divided our sample stars not only into bright and faint RC groups, but also into bluer ([Fe/H] < -0.1) and redder ([Fe/H] > -0.1) groups following the recent u-band photometric studies. For the metal-poor stars, no statistically significant difference in CN index was detected between the bright and faint RC groups for all observed fields. However, we found, from cross-matching with high-resolution spectroscopic data, a sign of Na enhancement in the "metal-poor and bright" RC group compared to the "metal-poor and faint" group at (l,b)=(-1 deg,-8.5 deg). When the contributions of the RGB stars on the RC regimes are taken into account, the Na abundance difference between genuine RCs would correspond to approximately 0.23 dex, similar to globular cluster (GCs) with multiple populations. In contrast, the metal-rich stars do not show chemical differences between the bright and faint RCs. It implies that the double RC observed in the metal-poor component of the bulge might be linked to the multiple populations originated from GC-like subsystem, whereas that of the metal-rich component would have produced by the X-shaped structure. Our results support the previous studies suggesting composite nature of the Milky Way bulge.

The Milky Way hosts at least two modes in its present day distribution of Fe and alpha-elements. The exact cause of this bimodality is disputed, but one class of explanations involves the merger between the Milky Way and a relatively massive satellite (Gaia-Sausage-Enceladus) at z~2. However, reproducing this bimodality in simulations is not straightforward, with conflicting results on the prevalance, morphology, and mechanism behind multimodality. We present a case study of a galaxy in the Illustris TNG50 simulation which undergoes sequential phases of starburst, brief quiescence, and then rejuvenation. This scenario results in a pronounced abundance bimodality after a post-processing adjustment of the [alpha/Fe] of old stars designed to mimic a higher star formation efficiency in dense gas. The high- and low-alpha sequences are separated in time by the brief quiescent period, which is not associated with a merger but by the formation of a bar followed by AGN activity. This galaxy indicates a novel scenario in which the alpha-bimodality in the Milky Way is caused by the formation of the bar via AGN-induced quenching. In addition to a stellar age gap in the Milky Way, we predict that abundance bimodalities should be more common in barred as opposed to unbarred galaxies.

We present a novel python-based 1D sub-Neptune evolution model that emphasizes the thermal evolution and potential solidification of the rock/iron core and the structure of the radiative atmosphere. This model explores planetary structure from the molten center to nbar pressure levels. Treating the radiative atmosphere is crucial for sub-Neptunes, due to the large scale height and low gravity, which contributes up to 40\% of their observed radius, especially for low-mass, highly irradiated planets. Consequently, we generically find that lower H/He mass fractions are needed to match a given planetary radius, compared to previous work. While the presence of metal-enrichment in the H/He layers (here modeled as 50$\times$ solar) does not substantially influence the size of the convective envelope, it notably reduces the transit radius by shrinking the radiative atmospheric scale height. Sub-Neptunes cool differently from terrestrial planets, with the rock/iron core's cooling rate limited by the envelope, leading to longer solidification timescales. Complete solidification of the silicate mantle by 10 Gyr is found only for planets with very low masses ($\leq 1M_\oplus$) and small H/He envelopes ($\leq$ 0.1\%). Dynamo action in sub-Neptune iron cores persists as long as the mantle surface remains molten, often exceeding 10 Gyr, and becomes sensitive to core thermal conductivity after solidification. We examine aspects of ''boil-off,'' which sets the maximum allowed H/He mass and planetary radius for subsequent evolution. The rock/iron's cooling energy moderately decreases the post-boil-off H/He mass fraction in planets with large atmospheric scale heights only.

Context. Orbital similarity measures, such as the D-values, have been extensively used in meteor science to identify meteoroid streams and associate meteorite falls with near-Earth objects (NEOs). However, the chaotic nature of near-Earth space challenges the long-term reliability of these measures for stream identification, and the increasing size of our fireball, meteorite fall, and NEO databases make random associations more common. Despite this, many researchers erroneously continue to use orbital similarity beyond its inherent limits. Aims. We aim to assess the statistical significance of using orbital similarity measures for identifying streams of meteoroids or asteroids and explore the implications of chaotic dynamics on the long-term coherence of these streams. Conclusions. The rapid decoherence of meteoroid streams and the chaotic dynamics of near-Earth orbits suggest that no reported stream or NEO associations of meteorites or fireballs are statistically significant according to orbital discriminates. Many are likely coincidental rather than indicative of a true physical link. However, several statistically significant clusters found within the NEO population are consistent with a tidal disruption formation. This contrast and lack of statistically significant associations amongst the impact datasets is likely due to the fireball databases being 2 orders of magnitude smaller than the NEO database and the higher intrinsic uncertainties of fireball observation derived orbits.

We present ALMA Band 7 continuum (340 GHz) and CO J=3-2 observations for an extended sample of disks in the Upper Scorpius OB Association (Upper Sco, age ~ 10 Myr). The targets were selected from previous studies that identified new members of Upper Sco using photometry and astrometry from the Gaia mission, and the presence of a disk has been inferred from mid-infrared excess emission. The new ALMA observations are combined with previous ALMA data to define a sample of 202 Upper Sco members with disks that have spectral types between G0 and M5.5. Among these sources, 120 (59%) have been detected in the continuum with a signal-to-noise ratio >= 3, and 83 (41%) have been detected in CO J=3-2. Both the continuum and CO J=3-2 fluxes show a strong correlation with the spectral type of the central star and the type of disk inferred from the shape of the infrared spectral energy distribution, where disks around earlier type stars and full disks are more luminous than disks around later type stars and evolved and debris disks. The median dust continuum luminosity is lower for disks in Upper Sco than in younger regions, as found in previous studies, where the differences are more pronounced in later spectral types (M4-M5) than in earlier spectral types.

We present the results of our multi-wavelength (X-ray to radio) follow-up campaign of the Einstein Probe transient EP240408a. The initial 10 s trigger displayed bright soft X-ray (0.5-4 keV) radiation with peak luminosity $L_\textrm{X} \gtrsim 10^{49}$ ($10^{50}$) erg s$^{-1}$ for an assumed redshift z>0.5 (2.0). The Neil Gehrels Swift Observatory and Neutron star Interior Composition ExploreR discovered a fading X-ray counterpart lasting for $\sim$5 d (observer frame), which showed a long-lived (~4 d) plateau-like emission ($t^{-0.5}$) before a sharp powerlaw decline ($t^{-7}$). The plateau emission was in excess of $L_\textrm{X} \gtrsim 10^{46}$ ($10^{47}$) erg s$^{-1}$ at z>0.5 (2.0). Deep optical and radio observations resulted in non-detections of the transient. Our observations with Gemini South revealed a faint potential host galaxy ($r \approx 24$ AB mag) near the edge of the X-ray localization. The faint candidate host, and lack of other potential hosts ($r \gtrsim 26$ AB mag; $J \gtrsim 23$ AB mag), implies a higher redshift origin (z>0.5), which produces extreme X-ray properties that are inconsistent with many known extragalactic transient classes. In particular, the lack of a bright gamma-ray counterpart, with the isotropic-equivalent energy ($10 - 10,000$ keV) constrained by GECam and Konus-Wind to $E_{\gamma,\textrm{iso}} \lesssim 4\times10^{51}$ ($6\times10^{52}$) erg at z>0.5 (2.0), conflicts with known gamma-ray bursts (GRBs) of similar X-ray luminosities. We therefore favor a jetted tidal disruption event (TDE) as the progenitor of EP240408a at z>1.0, possibly caused by the disruption of a white dwarf by an intermediate mass black hole. The alternative is that EP240408a may represent a new, previously unknown class of transient.

We present a 34-year timing solution of the redback pulsar system Terzan 5A (Ter5A). Ter5A, also known as B1744$-$24A or J1748$-$2446A, has a 11.56 ms pulse period, a $\sim$0.1 solar mass dwarf companion star, and an orbital period of 1.82 hours. Ter5A displays highly variable eclipses and orbital perturbations. Using new timing techniques, we have determined a phase-connected timing solution for this system over 34 years. This is the longest ever published for a redback pulsar. We find that the pulsar's spin variability is much larger than most globular cluster pulsars. In fact, of the nine redback pulsars with published or in preparation long-term timing solutions, Ter5A is by far the noisiest. We see no evidence of strong correlations between orbital and spin variability of the pulsar. We also find that long-term astrometric timing measurements are likely too contaminated by this variability to be usable, and therefore require careful short-term timing to determine reasonable positions. Finally, we measure an orbital period contraction of $-2.5(3) \times 10^{-13}$, which is likely dominated by the general relativistic orbital decay of the system. The effects of the orbital variability due to the redback nature of the pulsar are not needed to explain the observed orbital period derivative, but they are constrained to less than $\sim$30% of the observed value.

High precision ephemerides not only support space missions, but can also be used to study the origin and future of celestial bodies. In this paper, a coupled orbit rotation dynamics model that fully takes into account the rotation of the Martian moons is developed. Phobos and Deimos rotation are firstly described by Eulerian rotational equations, and integrated simultaneously with the orbital motion equations. Orbital and orientational parameters of Mars satellites were simultaneously obtained by numerical integration for the first time. In order to compare the differences between our newly developed model and the one now used in the ephemerides, we first reproduced and simulated the current model using our own parameters, and then fit it to the IMCCE ephemerides using least square procedures. The adjustment test simulations show Phobos and Deimos orbital differences between the refined model and the current model is no more than 300 meters and 125 meters, respectively. The orientation parameters are confirmed and the results are in good agreement with the IAU results. Moreover, we simulated two perturbations (main asteroids and mutual torques) which were not included in our refined model, and find that their effects on the orbits are completely negligible. As for the effect on rotation, we propose to take care of the role of mutual attraction in future models.

The evolution equation of two-point correlation function $\xi$ of galaxies can analytically describe the large scale structure of the galaxy distribution, and the solution depends also upon the initial condition. The primeval spectrum of the baryon acoustic oscillations (BAO) contains multi peaks that survived the Silk damping, and, as a relevant portion, two peaks of the primeval BAO spectrum fall into the range of current galaxy surveys. Incorporating this portion, we use a twin-peak initial power spectrum of the galaxies, and obtain the evolution solution from a redshift $z=8$ to $z=0$ in the Gaussian approximation. The outcome $\xi(r)$ at $z=0.6$ still exhibits the 100 Mpc periodic bumps as observed by the WiggleZ survey, a feature largely determined by the Jeans length $\lambda_J$ in the equation. In particular, due to the superposition of the twin peaks in the initial condition, $\xi(r)$ shows a shallow trough at $\sim 70 h^{-1}$Mpc and a deep trough at $\sim 140 h^{-1}$Mpc, agreeing with the observational data, much better than our previous work that used a simple one-peak initial spectrum.

Based on the data from the Kodaikanal and Mount Wilson observatories, we investigate the relationships of the tilt angle of sunspot group (SG), including the mean tilt angle and the tilt-angle scatter, during the declining phase with the parameters of the next solar cycle (SC). The main findings are summarized in the following three points. (1) During the declining phase, the correlation between the mean tilt angle and the tilt-angle scatter is statistically insignificant. (2) Six quantities measured during the declining phase show significant anti-correlations with the strength and amplitude of the next SC, and positive correlations with the duration of the ascending phase of the next SC: the standard deviation of tilt angles, the root-mean-square tilt angle, the mean absolute value of tilt angles, the area-weighted absolute value of tilt angles, the latitude-weighted absolute value of tilt angles, and the area- and latitude-weighted absolute value of tilt angles. (3) The correlations of the mean tilt angle, the area-weighted tilt angle, the latitude-weighted tilt angle, and the area- and latitude- weighted tilt angle during the declining phase with the strength, amplitude, and duration of the ascending phase of the next SC are statistically insignificant. These findings demonstrate that the modulation of parameters of the next SC by the tilt-angle scatter during the declining phase plays a vital role in regulating SC variability.

In this study, we investigate the impact of covariance within uncertainties on the inference of cosmological and astrophysical parameters, specifically focusing on galaxy stellar mass functions derived from the CAMELS simulation suite. Utilizing both Fisher analysis and Implicit Likelihood Inference (ILI), we explore how different covariance structures, including simple toy models and physics-motivated uncertainties, affect posterior distributions and parameter variances. Our methodology utilizes forward modeling via emulators that are trained on CAMELS simulations to produce stellar mass functions based on input parameters, subsequently incorporating Gaussian noise as defined by covariance matrices. We examine both toy model covariance matrices and physically motivated covariance matrices derived from observational factors like the stellar Initial Mass Function (IMF) and photometric aperture size. Our results demonstrate that covariance terms significantly influence parameter inference, often leading to tighter constraints or revealing complex, multimodal posterior distributions. These findings underscore the necessity of accounting for covariance when interpreting astrophysical observations, especially in fields where accurate parameter estimation is critical for model validation and hypothesis testing.

The traditional approach to analyzing mean motion resonances is through canonical perturbation theory. While this is a powerful method, its generality leads to complicated combinations of variables that are challenging to interpret and require looking up numerical coefficients particular to every different resonance. In this paper we develop simpler scaling relations in the limit where orbits are closely spaced (period ratios $\lesssim 2$) and interplanetary interactions can be approximated by only considering the close-approaches each time the inner planet overtakes the outer at conjunction. We develop geometric arguments for several powerful results: (i) that $p$:$p-q$ MMRs of the same order $q$ are all rescaled versions of one another (ii) that the general case of two massive planets on closely spaced, eccentric, co-planar orbits can be approximately mapped onto the much simpler case of an eccentric test particle perturbed by a massive planet on a co-planar circular orbit and (iii) that while the effects of consecutive conjunctions add up coherently for first-order ($p$:$p-1$) MMRs, they partially cancel for $p$:$p-q$ MMRs with order $q>1$, providing a physical explanation for why these higher order MMRs are weaker and can often be ignored. Finally, we provide simple expressions for the widths of MMRs and their associated oscillation frequencies that are universal to all closely spaced MMRs of a given order $q$, in the pendulum approximation.

We have developed a new method of multi-wavelength data combination for the search of late-type radio dwarfs, and have put it into practice using GLEAM-X DR1 data. The initial sample is selected by cross-matching the Gaia/DR3 objects with the probability of being a star no less than 99$\%$, and removing the extragalactic objects assigned by the SIMBAD database. The late-type dwarf stars are judged according to their location in the $(BP-RP)_0/M_{\rm G}$ color-magnitude diagram and in the $(J-H)_0/(K-W1)_0$ near-infrared color-color diagram. Furthermore, stellar activity is searched by ultraviolet excess in the GALEX/NUV band and the Rossby number in the TESS light curves. In total, 12 stars are found to be late-type dwarf stars associated with radio source, which is consisted of five stars with the UV excess and seven stars with the Rossby number less than 0.13. Three of these 12 stars are previously studied to be associated with radio objects. All these 12 stars are considered to be reliable counterparts of radio sources.

Blind spectroscopy of massive lensing galaxy clusters with MUSE has revealed large numbers of gravitationally-lensed Lyman-$ \alpha $ emitters exhibiting asymmetric profiles at $ 2.9 \leq z \leq 6.7 $, suggesting abundant outflows from low-mass star-forming galaxies in the early universe. Are these primaeval galaxies experiencing their first bursts of star formation, or established galaxies experiencing rejuvenation? With JWST rest-frame optical/NIR continuum imaging now available for many of these objects, we can search for older stellar populations. Here, we search for spectroscopic confirmation of outflows from these galaxies, finding a few high-signal-to-noise cases in which blueshifted interstellar absorption lines are detected. Next, we analyse the star formation histories with combined HST + JWST photometry. We find most them to be well characterised by very young, low metallicity stellar populations. However, despite the rest-frame optical/NIR coverage of JWST, we cannot place strict upper bounds on the mass in old stars (age $ > 100\,\text{Myr} $).

Calvera (1RXS J141256.0+792204) is a pulsar of characteristic age 285 kyr at a high Galactic latitude of b=+37°, detected only in soft thermal X-rays. We measure a new and precise proper motion for Calvera using Chandra HRC-I observations obtained 10 years apart. We also derive a new phase-connected ephemeris using 6 years of NICER data, including the astrometric position and proper motion as fixed parameters in the timing solution. Calvera is located near the center of a faint, circular radio ring that was recently discovered by LOFAR and confirmed as a supernova remnant (SNR) by the detection of gamma-ray emission with Fermi/LAT. The proper motion of $78.5 \pm 2.9$ mas/yr at position angle $241°.3 \pm 2°.2$ (in Galactic coordinates) points away from the center of the ring, a result which differs markedly from a previous low-significance measurement, and greatly simplifies the interpretation of the SNR/pulsar association. It argues that the supernova indeed birthed Calvera <10 kyr ago, with an initial spin period close to its present value of 59 ms. The tangential velocity of the pulsar depends on its uncertain distance, $v_t=(372 \pm 14) d_{1 kpc}$ km/s, but is probably dominated by the supernova kick, while its progenitor could have been a runaway O or B star from the Galactic disk.

Recent JWST observations have revealed a large population of intermediate/low-luminosity AGN at early times with peculiar properties, different from local AGN or luminous quasars. To better understand the physical conditions in the BLRs of these early AGN, we used the optical FeII (4434--4684 Å) and the broad $\rm H \beta$ emission, and the ratio between their equivalent widths $R_{Fe}$, as a probe on a purposefully assembled sample. Specifically, we gathered a sample of 26 high redshift ($\langle z \rangle$=6.4) AGN, observed by JWST, with broad $\rm H\beta$ detection both in the high and low luminosity regimes (respectively 14 faint AGN and 12 quasars), to investigate their optical FeII emission properties. In addition, we carefully selected control samples at lower $z$. We found that the population of faint AGN ($\rm \log(L_{H \beta} / (erg \, s^{-1}))\lesssim 44$) exhibits a significantly lower FeII emission than their local counterparts ($R_{Fe}<$0.24 versus $R_{Fe}\simeq$0.85 in the control sample), while the quasars at the epoch of reionisation observed by JWST present a FeII emission profile that closely resembles that observed at $z<3$. We argue that the weakness of the FeII bump in the faint JWST AGN might be due to the reduced metallicity of their broad line region ($\lesssim 0.5~Z_{\odot}$), while luminous quasars have already reached chemical maturity ($\sim Z{_\odot}$ or higher). Lastly, we highlight an intriguing similarity between the spectral properties of the high redshift population of faint AGN with those harboured in local metal poor dwarf galaxies.

Causality is an important factor determining the maximum mass of a neutron star. Previous works studies causality for smooth equation of state. The density at the core of neutron stars can be few times nuclear saturation density, where the occurrence of first-order phase transition has not been ruled out. The causality condition for first-order phase transition is characteristically different from the causality condition of smooth equation of state which can be evident in the mass-radius expression. In this letter we examine the causality condition for first-order phase transition and its observational significance in the mass-radius sequence. We find that equation of state having first-order phase transition, the causality line deviates considerably from the that of smooth equation of state. Depending on the strength of the discontinuity and the onset density of phase transition there is a narrow band in the mass-radius plot which is available only to the stars having smooth equation of state. This can have significant consequence in the sense that if some pulsars are to lie in this band one can rule out equation of state having first-order phase transition occurring at that particular onset density.

High precision differential astrometry assesses the positions, distances, and motions of celestial objects in relation to the stars. The focal plane of such space telescope must be calibrated with a precision down to the level of 1e-5 pixel in order to be able to detect Earth-like planets in the close vicinity of the Sun. The presented characterization bench is designed to improve the technology readiness level for the following key points: calibration of new detectors with a high number of pixels and correcting the field distortion using stars in the field of view. The first aim of the project concentrates on the characterization of a 46 megapixels sensor from PYXALIS, to assess its typical parameters using an integrating sphere. The next objective intends to map the intra and extra pixel quantum yield of the detector with a precision of 1e-5 pixels and investigate the evolution of the pixel geometry in response to environment fluctuations. To conduct these tests, an optical bench is designed with an LCD screen and a doublet, used as a source that allows directing light to specific groups of pixels. Interferometric calibration of the detector pixel centroid position will be achieved using fibers that illuminate the detector with Young's fringes. To characterize the distortion of the detector, a diaphragm will produce adjustable optical aberrations to be corrected and therefore change the source sensor positional relationship. The final step involves the simulation of a star's field, which will be imaged on the detector to assess optical quality.

We discuss the neutrino mass and Hubble tension solutions and examine their effects on the Redshift-Space Distortion (RSD) observations. An analysis with RSD data indicates smaller amplitude of perturbation. Including RSD data results in a slightly weaker upper limit on the neutrino mass than that derived for data without RSD, which is common in other extended models too. We have evaluated the impacts of RSD observations on some extended models, including the varying electron mass model, a time-dependent dark energy model with two parameter equations of state (EOS), and a model where the number of neutrino species is free. When we estimate the cosmological parameters for data including RSD, we found that the EOS parameter for dark energy is larger than that of the cosmological constant, and the effective number of neutrino species is smaller than the standard value, which infers a smaller present Hubble parameter $H_0$. From the viewpoint of cosmological tensions, the varying electron mass model with non-zero neutrino mass option looks promising to relax the Hubble tension and the $S_8$ tension simultaneously.

Jupiter's equatorial eastward zonal flows reach wind velocities of ~100 m/s, while on Saturn they are three times as strong and extend about twice as wide in latitude, despite the two planets being overall dynamically similar. Recent gravity measurements obtained by the Juno and Cassini spacecraft uncovered that the depth of zonal flows on Saturn is about three times greater than on Jupiter. Here we show, using 3D deep convection simulations, that the atmospheric depth is the determining factor controlling both the strength and latitudinal extent of the equatorial zonal flows, consistent with the measurements for both planets. We show that the atmospheric depth is proportional to the convectively driven eddy momentum flux, which controls the strength of the zonal flows. These results provide a mechanistic explanation for the observed differences in the equatorial regions of Jupiter and Saturn, and offer new understandings about the dynamics of gas giants beyond the Solar System.

The rapid cooling of the neutron star in Cassiopeia A is speculated to arise from an enhanced neutrino emission caused by the onset of $^3P_2$-wave neutron superfluidity in the core. However, the neutrino emissivity due to Cooper-pair breaking and formation is too small to yield the observed cooling rate. Here, we show that such a rapid cooling can be explained once the non-Fermi liquid behavior of the non-superfluid neutron liquid induced by superfluid quantum criticality is included into the theoretical description of neutron star cooling, without assuming the existence of additional energy loss processes. Our results indicate that the neutron star in Cassiopeia A remains in the thermal relaxation stage, which is greatly prolonged by the non-Fermi liquid behavior. The good agreement between our theoretical results and recent observational cooling data points to the pivotal role played by superfluid quantum criticality in neutron stars.

We present capivara, a fast and scalable multi-decomposition package designed to study astrophysical properties within distinct structural components of galaxies. Our spectro-decomposition code for analyzing integral field unit (IFU) data enables a more holistic approach, moving beyond conventional radial gradients and the bulge-plus-disk dichotomy. It facilitates comprehensive comparisons of integrated stellar ages and metallicities across various galactic structures. Our classification method naturally identifies outliers and organizes the different pixels based on their dominant spectral features. The algorithm leverages the scalability and GPU acceleration of Torch, outputting both a one-dimensional spectrum and a full data cube for each galaxy component, without relying on Voronoi binning. We demonstrate the capabilities of our approach using a sample of galaxies from the MaNGA survey, processing the resulting data cubes with the starlight spectral fitting code to derive both stellar population and ionized gas properties of the galaxy components. Our method effectively groups regions with similar spectral properties in both the continuum and emission lines. By aggregating the spectra of these regions, we enhance the signal-to-noise ratio of our analysis while significantly speeding up computations by reducing the number of spectra processed simultaneously. capivara will be freely available on GitHub.

In this PhD thesis we investigate stellar evolution and nucleosynthesis in the low- and extremely-low metallicity regime - including models of stars with a pure Big Bang composition (i.e. $\rm{Z} = 0$). The metallicity range of the extremely metal-poor (EMP) models calculated is $-6.5 < \rm{[Fe/H]} < -3.0$, with a mass range $0.85 < \rm{M} < 3.0~\rm{M}_{\odot}$. We have also calculated a series of models with a metallicity of $\rm{[Fe/H]} = -1.4$, to compare with observations of abundance patterns in Galactic globular cluster stars. Many of the extremely metal-poor (EMP) and $\rm{Z} = 0$ models experience violent evolutionary episodes not seen at higher metallicities. We refer to these events as `Dual Flashes' (DF) since they are characterised by peaks in the hydrogen and helium burning luminosities occurring at the same time. Some of the material processed by these events is later dredged up by the convective envelope, causing very significant surface pollution. We have calculated the entire evolution of the $\rm{Z} = 0$ and EMP models, including detailed nucleosynthesis and yields. Although subject to many uncertainties these are, as far as we are aware, the only yields available in this mass and metallicity range. We find that our models predict an increased number of carbon-rich stars at the lowest metallicities. This is mainly due to the extra pollution provided by the DF events - which do not occur in higher metallicity models. This concurs well with the observations that show the proportion of carbon-enhanced metal-poor (CEMP) stars in the Galactic Halo to be higher at lower metallicities. We also compare the chemical pollution arising from our models with the detailed abundance patterns available for some of the most metal-poor CEMP stars, and find mixed results. Fluid dynamics calculations are likely needed to model the violent DF episodes. [Abridged]

Blazars are the most common sources of $\gamma$-ray photons in the extra-galactic sky. Their $\gamma$-ray light-curves are characterized by bright flaring episodes, similarly to what is observed at longer wavelengths. These Gamma-Ray Bursts from Blazars (GRBBLs) have been extensively studied individually, but never in terms of a population. The goal of this work is to provide a global characterization of GRBBLs, to investigate the parameter space of the population, and ultimately to classify GRBBLs. Their global properties could give insights on the physical mechanisms responsible for the $\gamma$-ray radiation and on the origin of the observed variability. I analyze a sample of publicly available Fermi-LAT light-curves, utilizing only blazars with certain redshift measurements. The redshift-corrected light-curves are then automatically scanned to identify GRBBLs. A simple flare profile, with exponential rise and decay, is then fitted to all events. The fit parameters, together with global properties from the LAT catalog, are then used as input for unsupervised machine learning classification. The analysis shows that the GRBBL population is remarkably homogeneous. When using only the properties of the integral light-curves, the classifier converges into a single population. When adding information on the evolution of the photon index, the classifier splits the population into achromatic (the large majority) and chromatic (the outliers) GRBBLs. As by product of this study, I identify a correlation between the rising/decay time-scales of the GRBBLs and their peak luminosity.

Kilonovae, possible electromagnetic counterparts to neutron star mergers, provide important information about high-energy transient phenomena and, in principle, also allow us to obtain information about the source properties responsible for powering the kilonova. Unfortunately, numerous uncertainties exist in kilonova modeling that, at the current stage, hinder accurate predictions. Hence, one has to account for possible systematic modeling uncertainties when interpreting the observed transients. In this work, we provide a data-driven approach to account for time-dependent and frequency-dependent uncertainties in kilonova models. Through a suite of tests, we find that the most reliable recovery of the source parameters and description of the observational data can be obtained through a combination of kilonova models with time- and frequency-dependent systematic uncertainties. We apply our new method to analyze AT2017gfo. While recovering a total ejecta mass consistent with previous studies, our approach gives insights into the temporal and spectral evolution of the systematic uncertainties of this kilonova. We consistently find a systematic error below $1$ mag between $1$ to $5$ days after the merger. Our work addresses the need for early follow-up of kilonovae at earlier times, and improved modeling of the kilonova at later times, to reduce the uncertainties outside of this time window.

Robust modeling of non-linear scales is critical for accurate cosmological inference in Stage IV surveys. For weak lensing analyses in particular, a key challenge arises from the incomplete understanding of how non-gravitational processes, such as supernovae and active galactic nuclei - collectively known as baryonic feedback - affect the matter distribution. Several existing methods for modeling baryonic feedback treat it independently from the underlying cosmology, an assumption which has been found to be inaccurate by hydrodynamical simulations. In this work, we examine the impact of this coupling between baryonic feedback and cosmology on parameter inference at LSST Y1 precision. We build mock 3$\times$2pt data vectors using the Magneticum suite of hydrodynamical simulations, which span a wide range of cosmologies while keeping subgrid parameters fixed. We perform simulated likelihood analyses for two baryon mitigation techniques: (i) the Principal Component Analysis (PCA) method which identifies eigenmodes for capturing the effect baryonic feedback on the data vector and (ii) HMCode2020 (Mead et al. 2021) which analytically models the modification in the matter distribution using a halo model approach. Our results show that the PCA method is robust to the coupling between cosmology and baryonic feedback, whereas, when using HMCode2020 there can be up to $0.5\sigma$ bias in $\Omega_\text{m}$-$S_8$. For HMCode2020, the bias also correlates with the input cosmology while for PCA we find no such correlation.

We present an exploration of the Milky Way's structural parameters using an all-sky sample of RC giants to map the stellar density from the inner to the outer parts of the Galactic disc. These evolved giants are considered to be standard candles due to their low intrinsic variance in their absolute luminosities, allowing us to estimate their distances with reasonable confidence. We exploit all-sky photometry from the AllWISE mid-infrared survey and the Gaia survey, along with astrometry from Gaia Data Release 3 and recent 3D extinction maps, to develop a probabilistic scheme in order to select with high confidence \rc{}-like stars. Our curated catalogue contains about 10 million sources, for which we estimate photometric distances based on the WISE $W1$ photometry. We then derive the selection function for our sample, which is the combined selection function of sources with both \gaia{} and \allwise{} photometry. Using the distances and accounting for the full selection function of our observables, we are able to fit a two-disc, multi-parameter model to constrain the scale height (\hz{}), scale-length (\rd{}), flaring, and the relative mass ratios of the two disc components. We illustrate and verify our methodology using mock catalogues of \rc{} stars. We find that the \rc{} population is best described by a flared thin disc with scale length \rd{}=$3.56\pm0.32$ kpc and scale height at the Sun of \hzsun{}=$0.17\pm0.01$ kpc, and a shorter and thicker disc with \rd{}=$2.59\pm0.11$ kpc, \hzsun{}=$0.45\pm0.11$ kpc, with no flare. The thicker disc constitutes 64\% of the \rc{} stellar mass beyond 3 kpc, while the thin disk shows evidence of being warped beyond 9 kpc from the Galactic center. The residuals between the predicted number density of RC stars from our axisymmetric model and the measured counts show possible evidence of a two-armed spiral perturbation in the disc of the Milky Way.

Infrared interferometry has seen a revolution over the last few years. The advent of GRAVITY+ is about to enable high-contrast observations, all-sky coverage and faint science up to K=21, with the implementation on 8m-class telescope of extreme adaptive optics, wide-field observations, and soon laser guide stars, following a long-term vision of technological and infrastructure development at VLTI. This major progress in sensitivity lift a fundamental limitation of infrared interferometry, namely the brightness temperature achievable with this technique down to milli-arcsecond resolution imaging. This change of paradigm is a crucial element for the expansion of current arrays to a facility up to one to ten kilometer baselines. Micro-arcsecond scales imaging in the infrared on thermal objects, reaching the highest angular resolution possible even compared to VLBI, could offer a unique window in observational astronomy for the next generation instrument.

We present the first LOFAR image of the centre of M31 at a frequency of 150 MHz. We clearly detect three supernova remnants, which, along with archival VLA data at 3 GHz and other published radio and X-ray data allows us to characterize them in detail. Our observations also allow us to obtain upper limits the historical SN 1885A which is undetected even at a low frequency of 150 MHz. From analytical modelling we find that SN 1885A will stay in its free-expansion phase for at least another couple of centuries. We find an upper limit of $n_{\rm H}~\lesssim 0.04$ cm$^{-3}$ for the interstellar medium of SN 1885A, and that the SN ejecta density is not shallower than $\propto r^{-9}$ (on average). From the $2.6\sigma$ tentative detection in X-ray, our analysis shows that non-thermal emission is expected to dominate the SN 1885A emission. Comparing our results with those on G1.9+0.3, we find that it is likely that the asymmetries in G1.9+0.3 make it a more efficient radio and X-ray emitter than SN 1885A. For Braun 80, 95 and 101, the other remnants in this region, we estimate ages of 5200, 8100, and 13,100 years, and shock speeds of 1150, 880, and 660 km s$^{-1}$}, respectively. Based on this, the supernova rate in the central 0.5 kpc $\times$ 0.6 kpc of M31 is at least one per $\sim 3000~{\rm yr}$. We estimate radio spectral indices of $-0.66\pm0.05$, $-0.37\pm0.03$ and $-0.50\pm0.03$ for the remnants, respectively, which match fairly well with previous studies.

In single-lined spectroscopic binaries (SB1s) where flux variations due to tidal deformation of the primary star (ellipsoidal variations, EVs) are detected, the binary mass can be determined by combining EVs with the primary's radial velocity (RV) variations from orbital motion and information about the primary's radius. This method has been used for mass estimation in close binaries including X-ray systems, but it has been pointed out that contaminating light from sources other than the primary star could introduce systematic errors in the mass and inclination estimates. Here, we focus on the apparent RV variations caused by asymmetric distortion of the absorption lines of the tidally deformed primary star (tidal RV). Because this signal contains information equivalent to that from photometric EVs, it enables mass estimation of the binary system using only the primary star's absorption lines from high-resolution spectroscopic data, providing a potentially more robust approach against contaminating light. We apply the method to the binary system V723 Monocerotis, where both photometric EV and tidal RV signals are detected, and successfully determine the component masses using only the primary star's RVs and projected rotational velocity, without relying on absolute flux measurements or on stellar evolutionary models. The masses derived from the tidal RV model show a reasonable agreement with those obtained from EVs after carefully modeling the flux contamination from the secondary. This result demonstrates that tidal RVs provide a useful alternative means for mass estimation in SB1s.

We investigate the production of dark radiation (DR) from axions and axion-like particles (ALPs) as potential origins of dark matter. Focusing on the dark matter misalignment mechanism, we examine non-thermal, pre-inflationary scenarios that could lead to the generation of DR. A key part of our analysis involves a Bayesian approach to confront ALP parameter space with current cosmological data. Additionally, we explore how DR production could offer solutions to persistent cosmological anomalies, particularly those related to the Hubble constant $H_0$ tension and the $S_8$ parameter discrepancy. Our findings aim to shed new light on the role of ALPs in addressing these open issues in cosmology.

The Local Universe (D < 120 Mpc) has been intensely studied for decades, with highly complete galaxy redshift surveys now publicly available. These data have driven density reconstructions of the underlying matter density field, as well as constrained simulations that aim to reproduce the observed structures. In this paper, we introduce a dispersion measure (DM) model that makes use of this detailed knowledge of our Local Universe within D < 120 Mpc. The model comprises three key components: (i) the DM from the Milky Way halo and the intra-group medium (up to 3.4 Mpc), derived from the HESTIA simulations, a series of constrained hydrodynamic simulations designed to reproduce our Local Group; (ii) the DM contribution from the large-scale intergalactic medium beyond the Local Group (3.4 Mpc < D < 120 Mpc), calculated using the HAMLET reconstructed matter density field; and (iii) the individual DM contributions from Local Universe galaxy halos and clusters based on data from the 2MASS Galaxy Group Catalog and the NASA/IPAC Extragalactic Database. This comprehensive model will be made available as a Python package. As the most realistic model to date for DM in the local volume, it promises to improve the constraints of DM contributions from the Intergalactic Medium and Circumgalactic Medium of FRBs, thereby enhancing the accuracy of cosmic baryon distribution calculations based on DM analysis of FRBs.

Many nucleosynthetic channels create the elements, but two-parameter models characterized by $\alpha$ and Fe nonetheless predict stellar abundances in the Galactic disk to accuracies of 0.02 to 0.05 dex for most measured elements, near the level of current abundance uncertainties. It is difficult to make individual measurements more precise than this to investigate lower-amplitude nucleosynthetic effects, but population studies of mean abundance patterns can reveal more subtle abundance differences. Here we look at the detailed abundances for 67315 stars from APOGEE DR17, but in the form of abundance residuals away from a best-fit two-parameter, data-driven nucleosynthetic model. We find that these residuals show complex structures with respect to age, guiding radius, and vertical action that are not random and are also not strongly correlated with sources of systematic error such as surface gravity, effective temperature, and radial velocity. The residual patterns, especially in Na, C+N, Ni, Mn, and Ce, trace kinematic structures in the Milky Way, such as the inner disk, thick disk, and flared outer disk. A principal component analysis suggests that most of the observed structure is low-dimensional and can be explained by a few eigenvectors. We find that some, but not all, of the effects in the low-$\alpha$ disk can be explained by dilution with fresh gas, so that abundance ratios resemble those of stars with higher metallicity. The patterns and maps we provide could be combined with accurate forward models of nucleosynthesis, star formation, and gas infall to provide a more detailed picture of star and element formation in different Milky Way components.

Early JWST observations have revealed the ubiquitous presence in the early Universe, up to z about16, of extreme baryon concentrations, namely forming globular clusters, extremely dense galaxies that may or may not be UV bright, and supermassive black holes in relatively low-mass galaxies. This paper is trying to pinpoint which physical conditions may have favored the formation of such concentrations, that appear to be very common at high redshifts while their formation being progressively more and more rare at lower redshifts. Building on local globular cluster evidence, it is argued that such conditions can consist in a combination of a 10 Myr extended feedback free time, coupled to low angular-momentum densities in deep local minima of the ISM vorticity field, where baryon concentrations are more likely to form. It is argued that the former condition would follow from more massive stars failing to explode as supernovae, and the latter one from low vorticity prevailing in the early Universe, in contrast to later times with their secular increase of the angular momentum density due to the cumulative effect of tidal interactions.

We present the photometric performance of SPIRIT, a ground-based near-infrared InGaAs CMOS-based instrument (1280 by 1024 pixels, 12 micron pitch), using on-sky results from the SPECULOOS-Southern Observatory during 2022 - 2023. SPIRIT was specifically designed to optimise time-series photometric precision for observing late M and L type stars. To achieve this, a custom wide-pass filter (0.81 - 1.33 microns, zYJ ) was used, which was also designed to minimise the effects of atmospheric precipitable water vapour (PWV) variability on differential photometry. Additionally, SPIRIT was designed to be maintenance-free by eliminating the need for liquid nitrogen for cooling. We compared SPIRIT's performance with a deeply-depleted (2048 by 2048 pixels, 13.5 micron pitch) CCD-based instrument (using an I+z' filter, 0.7 - 1.1 microns) through simultaneous observations. For L type stars and cooler, SPIRIT exhibited better photometric noise performance compared to the CCD-based instrument. The custom filter also significantly minimised red noise in the observed light curves typically introduced by atmospheric PWV variability. In SPIRIT observations, the detector's read noise was the dominant limitation, although in some cases, we were limited by the lack of comparison stars.

Sudden phase transitions during inflation can give rise to strongly enhanced primordial density perturbations on scales much smaller than those directly probed by cosmic microwave background anisotropies. In this paper, we study the effect of the incoming quantum state on the steepest growth found in the primordial power spectrum using a simple model of an instantaneous transition during single-field inflation. We consider the case of a general de Sitter-invariant initial state for the inflaton field (the $\alpha$-vacuum), and also an incoming state perturbed by a preceding transition. For the $\alpha$-vacua we find that $k^6$ growth is possible for $\alpha>0$, while $k^4$ growth is seen for $\alpha\leq0$, including the standard case of an initial Bunch-Davies vacuum state. The features of an enhanced primordial power spectrum on small scales are thus sensitive to the initial quantum state during inflation. We calculate the scalar-induced gravitational wave power spectrum for each case.

In the present-day universe, the most massive galaxies are ellipticals located in the cores of galaxy clusters, harboring the heaviest super-massive black holes (SMBHs). However the mechanisms that drive the early growth phase and subsequent transformation of these morphology and kinematics of galaxies remain elusive. Here we report (sub)kiloparsec scale observations of stars, gas, and dust in ADF22.A1, a bright dusty starburst galaxy at z=3.1, hosting a heavily obscured active galactic nucleus and residing in a proto-cluster core. ADF22.A1 is a giant spiral galaxy with the kinematics of a rotating disk with rotation velocity Vrot=530+/-10km/s and diameter larger than 30 kpc. The high specific stellar angular momentum of this system, j*=3400+/-600 kpc km/s, requires a mechanism to effectively spin-up ADF22.A1, indicating the importance of accretion from the cosmic web to supply both gas and angular momentum to galaxies in their early gas-rich starburst phase. In its inner region, gas flows along dust lanes in a bar connected with the bright dusty core and the estimated mass ratio of a bulge to SMBH matches the local relation, suggesting that bars are a key mechanism to shape the early co-evolution of these components. Comparison with cosmological simulations shows that ADF22.A1 will likely evolve into a massive elliptical at the present day, experiencing a significant reduction in angular momentum associated with subsequent galaxy mergers.

Statistical studies show that stars of GK spectral types, with masses below 1.1 Sun mass, are depleted in hot Jupiters. This finding is evidence of tidal orbital decay during the main-sequence lifetime. Theoretical considerations show that in some configurations, the tidal energy dissipation can be boosted by non-linear effects in dynamical tides, which are wave-like responses to tidal forcing. To probe the regime of these dynamical tides in GK stars, we searched for orbital period shortening for 6 selected hot Jupiters in systems with 0.8-1 Sun mass host stars: HATS-18, HIP 65A, TrES-3, WASP-19, WASP-43, and WASP-173A. For the hot Jupiters of our sample, we analysed transit timing data sets based on mid-transit points homogeneously determined from observations performed with the Transiting Exoplanet Survey Satellite and high-quality data available in the literature. For the TrES-3 system, we also used new transit light curves we acquired with ground-based telescopes. The mid-transit times were searched for shortening of orbital periods through statistical testing of quadratic transit ephemerides. Theoretical predictions on the dissipation rate for dynamical tides were calculated under the regimes of internal gravity waves (IGWs) undergoing wave breaking (WB) in stellar centres and weak non-linear (WNL) wave-wave interactions in radiative layers. Stellar parameters of the host stars, such as mass and age, which were used in those computations, were homogeneously redetermined using evolutionary models with the Bayesian inference. We found that transit times follow the refined linear ephemerides for all ultra-hot Jupiters of our sample. Non-detection of orbital decay allowed us to place lower constraints on the tidal dissipation rates in those planet-star systems. In three systems, HATS-18, WASP-19, and WASP-43, we reject a scenario with total dissipation of IGWs. We conclude that...

Radio pulsars allow the study of the ionised interstellar medium and its dispersive effects, a major noise source in gravitational wave searches using pulsars. In this paper, we compare the functionality and reliability of three commonly used schemes to measure temporal variations in interstellar propagation effects in pulsar-timing data. We carry out extensive simulations at low observing frequencies (100-200 MHz) by injecting long-term correlated noise processes with power-law spectra and white noise, to evaluate the robustness, accuracy and precision of the following three mitigation methods: epoch-wise (EW) measurements of interstellar dispersion; the DMX method of simultaneous, piece-wise fits to interstellar dispersion; and DMGP, which models dispersion variations through Gaussian processes using a Bayesian analysis method. We then evaluate how reliably the input signals are reconstructed and how the various methods react to the presence of achromatic long-period noise. All the methods perform well, provided the achromatic long-period noise is modeled for DMX and DMGP. The most precise method is DMGP, followed by DMX and EW, while the most accurate is EW, followed by DMX and DMGP. We also test different scenarios including simulations of L-band ToAs and realistic DM injection, with no significant variation in the obtained results. Given the nature of our simulations and our scope, we deem that EW is the most reliable method to study the Galactic ionized media. Future works should be conducted to confirm this result via more realistic simulations. We note that DM GP and DMX seem to be the most performing techniques in removing long-term correlated noise, and hence for gravitational wave studies. However, full simulations of pulsar timing array experiments are needed to support this interpretation.

We present the findings of a comprehensive and detailed analysis of merger tree data from ultra-high-resolution cosmological $N$-body simulations. The analysis, conducted with a particle mass resolution of $5 \times 10^3 h^{-1} M_{\odot}$ and a halo mass resolution of $10^7 h^{-1} M_{\odot}$, provides sufficient accuracy to suppress numerical artefacts. This study elucidates the dynamical evolution of subhaloes associated with the Milky Way-like host haloes. Unlike more massive dark matter haloes, which have been extensively studied, these subhaloes follow a distinct mass evolution pattern: an initial accretion phase, followed by a tidal stripping phase where mass is lost due to the tidal forces of the host halo. The transition from accretion to stripping, where subhaloes reach their maximum mass, occurs around a redshift of $z\simeq1$. Smaller subhaloes reach this point earlier, while larger ones do so later. Our analysis reveals that over 80 per cent of subhaloes have experienced mass loss, underscoring the universality of tidal stripping in subhalo evolution. Additionally, we derived the eccentricities and pericentre distances of subhalo orbits from the simulations and compare them with those of nearby satellite galaxies observed by the Gaia satellite. The results demonstrate a significant alignment between the orbital elements predicted by the cold dark matter model and the observed data, providing robust support for the model as a credible candidate for dark matter.

We present a novel implementation of a soft sphere, discrete elements code to simulate the dynamics of self-gravitating granular materials. The code is used to study the outcome of sub-sonic collisions between self-gravitating rubble piles with masses ranging from $\sim6\times10^{21}$ to $\sim6\times10^{22}$ g. These masses are representative of asteroids and planetesimals in the $\sim100$~km range. We simulate rubble piles composed of a range of particle sizes and analyze the collisions outcome focusing on the properties of the largest surviving fragment. We successfully test and validate the code against previous results. The results of our study show that rubble piles formed by collision of two parent rubble piles do not maintain the same particle size distribution as their parents. Rubble piles formed in low velocity collisions are characterized by a larger fraction of large particles, while the largest fragments of high-velocity collisions show a significant decrease in their mean particle size. We ascribe this effect to the fact that large particles transmit most of the forces during the collisions. In addition, we find that the mass of the largest post-collision fragment does depend on the rotation of the colliding rubble piles. This effect is especially noticeable when the pre-collision spin axes are parallel with each other and perpendicular to the relative velocity. This finding can be particularly relevant for meter to kilometer sized bodies embedded in protostellar accretion disks, where viscous stresses can efficiently align the target and projectile spin axes.

Recent cosmic ray (CR) measurements have revealed unexpected anomalies in secondary CRs, namely deviations from the predictions of the so-called standard Galactic CR paradigm regarding the composition and energy spectra of the products of interactions of primary (accelerated) CRs with interstellar gas: (i) antiparticles (positrons and antiprotons), (ii) light elements of the (Li, Be, B) group, and (iii) diffuse gamma rays. We argue that the new measurements can still be explained within the standard CR paradigm but with an additional assumption that CRs spend a significant part of their lifetime near their formation sites. The latter can be realized if CRs propagate more slowly in these localized regions than in the interstellar medium (ISM). Postulating that CRs accumulate on average energy-independent "grammage" of $0.7 \ \rm g/cm^2$ near the major contributors to galactic CRs, one can explain self-consistently the new measurements of the B/C ratio by DAMPE and the diffuse ultra-high-energy gamma-rays by LHAASO, involving a minimal number of model parameters: the energy-dependent "grammage" in the interstellar medium $\rm \lambda \approx 8 (E/10 \ GeV)^{-0.55}~\rm g/cm^{2}$ and the average CR acceleration (sourcee) spectrum, $\rm Q(E) \propto E^{-2.3}$.

We report a time-resolved analysis of the accreting X-ray pulsar Cen X-3 using observations carried out by NuSTAR, which cover approximately two binary orbits in different intensity states. The pulse profile is relatively stable over the orbital phase and shows energy dependence. It has an obvious double-peaked shape in the energy band below 15 keV -- with the second pulse peak decreasing as energy increases -- and is gradually dominated by a single peak in higher energy bands. We find that the pulse profile in the energy band of 3-5 keV at high-intensity states shows a subtle triple-peaked shape, with the main peak divided into two subpeaks. We also find a positive correlation between the pulse fraction and both energy and flux. Our spectral analysis reveals that the spectra can be well described by the continuum of Fermi-Dirac cutoff and NPEX models, and the cyclotron line is detected with the centroid energies varying from 26 keV to 29 keV, along with the iron emission line around 6.4 keV. We investigated the dependence between the cyclotron resonant scattering feature (CRSF) centroid energy and luminosity and discuss the theoretical critical luminosity. Although the variation of $E_{\rm cyc}- L_X$ is not distinct, there is a possibility that the critical luminosity lies within the range of $\sim (0.5-4)\times 10^{37}$ erg s$^{-1}$ in the band of $4-78$ keV. The photon index shows a strong positive correlation with luminosity. Our orbital-phase analysis reveals that the spectral parameters show orbital variability, and the highly variable photoelectric absorption may indicate the existence of clumpy stellar winds.

As a new generation of large-sky spectroscopic surveys comes online, the enormous data volume poses unprecedented challenges in classifying spectra. Modern unsupervised techniques have the power to group spectra based on their dominant features, circumventing the complete reliance on training data suffered by supervised methods. We outline the use of dimensionality reduction to generate a 2D map of the structure of an intermediate-resolution spectroscopic dataset. This technique efficiently separates white dwarfs of different spectral classes in the Dark Energy Spectroscopic Instrument's Early Data Release (DESI EDR), identifying spectral features that had been missed even by visual classification. By focusing the method on particular spectral regions, we identify white dwarfs with helium features at 90 per cent recall, and cataclysmic variables at 100 per cent recall, illustrating rapid selection of low-contamination samples from spectroscopic surveys. We also demonstrate the use of dimensionality reduction in a supervised manner, outlining a procedure to classify any white dwarf spectrum in comparison with those in the DESI EDR. With upcoming surveys promising tens of millions of spectra, our work highlights the potential for semi-supervised techniques as an efficient means of classification and dataset visualisation.

After over three decades of unsuccessful attempts, we report the first detection of molecular gas emission in Malin 1, the largest spiral galaxy observed to date, and one of the most iconic giant low surface brightness galaxies. Using ALMA, we detect significant $^{12}$CO(J=1-0) emission in the galaxy's central region and tentatively identify CO emission across three regions on the disc. These observations allow for a better estimate of the H$_2$ mass and molecular gas mass surface density, both of which are remarkably low given the galaxy's scale. By integrating data on its HI mass, we derive a very low molecular-to-atomic gas mass ratio. Overall, our results highlight the minimal presence of molecular gas in Malin 1, contrasting sharply with its extensive, homogeneous atomic gas reservoir. For the first time, we position Malin 1 on the Kennicutt-Schmidt (K-S) diagram, where it falls below the main sequence for normal spirals, consistent with previous upper limits but now with more accurate figures. These findings are crucial for constraining our understanding of star formation processes in environments characterized by extremely low molecular gas densities and for refining models of galaxy formation, thereby improving predictions concerning the formation, evolution, and distribution of these giant, elusive galaxies.

The innermost portions of the Milky Way's stellar halo have avoided scrutiny until recently. The lack of wide-area survey data, made it difficult to reconstruct an uninterrupted view of the density distribution of the metal-poor stars inside the Solar radius. In this study, we utilize red giant branch (RGB) stars from Gaia, with metallicities estimated using spectro-photometry from Gaia Data Release 3. Accounting for Gaia's selection function, we examine the spatial distribution of metal-poor ([M/H]<-1.3) RGB stars, from the Galactic centre (r~1 kpc) out to beyond the Solar radius (r~18 kpc). Our best-fitting single-component cored power-law model shows a vertical flattening of ~0.5 and a slope -3.4, consistent with previous studies. Motivated by the mounting evidence for two distinct stellar populations in the inner halo, we additionally test a range of two-component models. One of the components models the tidal debris from the Gaia Sausage/Enceladus merger, while the other captures the Aurora population -- stars that predate the Galactic disk formation. Our best-fit two-component model suggests that both populations contribute equally around the Solar radius, but Aurora dominates the inner halo with a steeper power-law index of -4.5, in agreement with the nitrogen-rich star distribution measured by Horta et al. (2021).

The accretion luminosity of an FU Ori disk is a fundamental system parameter, but a challenging one to estimate for all but the most well-studied systems. FU Ori objects are dynamically evolving accretion disks, especially close in time to the outburst epoch. They have a complex multi-temperature disk structure that results in distinctly shaped, broad SEDs. Detailed spectroscopic analysis is required for simultaneous constraint on relevant physical parameters such as the central stellar mass, inner disk radius, disk inclination, and disk accretion rate. However, outbursting systems that are deeply embedded and/or distant may be limited to only photometric measurement, and over only a narrow range of wavelengths. The bolometric corrections necessary to estimate accretion luminosities are not straightforward, and in particular can not be adopted from existing literature on isotropically radiating stars. We present bolometric corrections specific to astrophysical accretion disks for a variety of filters in ongoing and upcoming all-sky surveys.

Many promising explosion models for the elusive origin of Type Ia supernovae (SNe Ia) ultimately fail to completely reproduce a number of observed properties of these events. One limiting factor for many of these models is the use of the local thermodynamic equilibrium (LTE) assumption in the calculation of their synthetic observables, which has been shown to prevent the accurate prediction of a number of fundamental features of SNe Ia. The inclusion of high-accuracy non-LTE physics, however, increases computational cost and complexity such that multidimensional non-LTE calculations are often unfeasible, which can be problematic for models that are inherently multidimensional. In this work, we conduct radiative transfer calculations using 1D profiles that each correspond with a line of sight from an asymmetric, 2D SN Ia model. We find, in LTE, that the synthetic observables from these calculations efficiently reproduce those from the 2D calculation when an equivalence of bolometric luminosities between the 1D and 2D treatments is enforced. This allows for the accurate calculation of synthetic observables in 1D while still preserving multidimensional effects associated with the model. We leverage this to produce high accuracy observables from 1D non-LTE calculations, showing significantly improved agreement with observation, including a roughly 50% reduction of $B$-band decline rate into congruence with the observed Phillips relation. Additionally, our non-LTE observables show Si II $\lambda$5972 pEWs that are much more similar to observation, while spanning multiple Branch classes, suggesting that some spectral classifications of SNe Ia arise from line of sight effects.

We report high-resolution spectroscopic monitoring and long-baseline interferometric observations with the PTI of the 215-day binary system HD 174881 (K1 II-III), composed of two giant stars. The system is spatially resolved with the PTI, as well as in archival measurements with the CHARA Array. Our analysis of these observations, along with an analysis of the spectral energy distribution, have allowed us to infer accurate values for the absolute masses ($3.367^{+0.045}_{-0.041} M_{\odot}$ and $3.476^{+0.043}_{-0.043} M_{\odot}$), radii ($34.0 \pm 1.3 R_{\odot}$ and $22.7 \pm 1.8 R_{\odot}$), effective temperatures ($4620 \pm 100$ K and $4880 \pm 150$ K), and bolometric luminosities of both components, as well as other properties including the orbital parallax (distance). These provide valuable tests of stellar evolution models for evolved stars, which are still relatively uncommon compared to the situation for main-sequence stars. We find generally good agreement of all of these properties of HD 174881 with two sets of recent models (MIST, and PARSEC) at compositions near solar, for ages of 255-273 Myr. We also find evidence of an infrared excess, based largely on the flux measurements from IRAS at 60 and 100 microns.