Papers for Wednesday, Nov 20 2024

In the monograph, the authors systematize and generalize the results of studying solar-type activity that is characteristic of a significant part of mid- and low-mass stars of the Galaxy, outline the characteristics of such stars in the quiescent state, during the sporadic flares and variations of magnetic activity over the course of stellar evolution. The observational data taken over the whole range of the electromagnetic spectrum from decametric radio waves to X-rays are described in detail. The current models of stellar flares and stellar dynamo models are considered in two theoretical chapters. The last part of the book describes the impact of stellar activity of mid- and low-mass stars on exoplanetary environments. The Catalog of Stars with Solar-Type Activity is provided in the supplementary section. The book is intended for researchers involved in studying the physics of stars and the Sun, graduate students, and students specialized in heliophysics, astrophysics, and in the field of space physics. The cover shows the distribution of 314618 stars with solar-type activity over the celestial sphere.

Multiband observations of compact object sources offer a unique opportunity to explore their progenitors and enhance early multi-messenger alert. Recent analyses have indicated that metallicity significantly impacts the evolution of progenitors and the resulting compact objects. Using binary population synthesis, we investigate the formation of eccentric, inspiralling black hole binaries and black hole-neutron star binaries through the isolated binary evolution channel. We introduced a fiducial mass and metallicity relation for each ZAMS star. We model the stellar cluster of ZAMS stars by extending COSMIC's publicly available code. Our BPS code effectively accounts for the metallicity of each stellar object in the stellar cluster. In our analysis, we observed a significant increase in the number of inspiral binaries remaining in the stellar cluster. Instead of assuming a uniform metallicity for a stellar cluster, ZAMS stars within the cluster, characterized by diverse metallicity, evolve into more massive compact objects. The total mass of a single binary black hole inspiral varies from $\sim 9-86$ M$_\odot$; whereas for a black hole-neutron star system, this range becomes $\sim 6-32$ M$_\odot$. We compare the detectability of the characteristic strain against sub-Hz gravitational wave detectors.

Recent observations from the Juno spacecraft during its transit over flux tubes of the Galilean moons have identified sharp enhancements of particle fluxes at discrete energies. These banded structures have been suspected to originate from a bounce resonance between particles and standing Alfven waves generated by the moon-magnetospheric interaction. Here, we show that predictions from the above hypothesis are inconsistent with the observations, and propose an alternative interpretation that the banded structures are remote signals of particle absorption at the moons. In this scenario, whether a particle would encounter the moon before reaching Juno depends on the number of bounce cycles it experiences within a fixed section of drift motion determined by moon-spacecraft longitudinal separation. Therefore, the absorption bands are expected to appear at discrete, equally-spaced velocities consistent with the observations. This finding improves our understanding of moon-plasma interactions and provides a potential way to evaluate the Jovian magnetospheric models.

We present a novel density profile for halos in self-interacting dark matter (SIDM) models, which accurately captures the flat- and isothermal-core configurations. We show analytically how our density profile satisfies these conditions, with comparisons to other contemporary functional choices. We demonstrate the versatility of our profile by putting it into the context of idealized N-body simulations and show that it provides excellent representations for both density and velocity dispersion structures of the simulation data. When an estimated fitting criterion is used to approximate the general cases, such as in cosmological simulations, the resulting regressions maintain their goodness of fit in both extremes, in the initial thermalization phase and the late core-collapse regime. Our density profile provides a framework for more detailed analyses of halos in different SIDM models while serving as the basis for reducing simulation needs and constructing initial conditions for deep core-collapse simulations.

For the past decade, the high-contrast observation of disks and gas giant planets around nearby stars has been made possible on ground-based instruments using extreme adaptive optics (XAO). While these facilities produce images with a Strehl ratio larger than 90% in H-band in median observing conditions and high-flux regime, the correction leaves AO residuals which prevent the study of fainter or less massive exoplanets. Cascade AO systems with a fast second stage based on a Pyramid wavefront sensor have recently emerged as an appealing solution to reduce the atmospheric wavefront errors. Since these aberrations are expected to be small, they can also be accurately measured by a Zernike wavefront sensor (ZWFS), a well-known concept for its high sensitivity and moderate linear capture range. We propose an alternative second stage relying on the ZWFS to correct for the AO residuals. The cascade AO with a ZWFS-based control loop is implemented on the ESO's GHOST testbed to validate the scheme in monochromatic light. In median wind speed and seeing, our second-stage AO with a ZWFS and a basic integrator reduce the atmospheric residuals by a factor of 6 and increase the wavefront error stability with a gain of 2 from open to closed loop. In the presence of non common path aberrations, we also reach a contrast gain by a factor of 2 in the images with a Lyot coronagraph at short separations from the source, proving the ability of our scheme to work in cascade with an XAO loop and help imaging fainter or lighter close-in companions. In more challenging conditions, contrast improvements are also achieved by adjusting the control loop features. Our study validates the ZWFS-based second-stage AO loop as an effective solution to address small residuals left from a single-stage XAO system for the coronagraphic observations of circumstellar environments.

It is generally accepted that radio relics are the result of synchrotron emission from shock-accelerated electrons. Current models, however, are still unable to explain several aspects of their formation. In this paper, we focus on three outstanding problems: i) Mach number estimates derived from radio data do not agree with those derived from X-ray data, ii) cooling length arguments imply a magnetic field that is at least an order of magnitude larger than the surrounding intracluster medium (ICM), and iii) spectral index variations do not agree with standard cooling models. To solve these problems, we first identify typical shock conditions in cosmological simulations, using the results to inform significantly higher resolution shock-tube simulations. We apply the cosmic ray electron spectra code CREST and the emission code CRAYON+ to these, thereby generating mock observables ab-initio. We identify that upon running into an accretion shock, merger shocks generate a shock-compressed sheet, which, in turn, runs into upstream density fluctuations in pressure equilibrium. This mechanism directly gives rise to solutions to the three problems: it creates a distribution of Mach numbers at the shock-front, which flattens cosmic ray electron spectra, thereby biasing radio-derived Mach number estimates to higher values. We show that this effect is particularly strong in weaker shocks. Secondly, the density sheet becomes Rayleigh-Taylor unstable at the contact discontinuity, causing turbulence and additional compression downstream. This amplifies the magnetic field from ICM-like conditions up to microgauss levels. We show that synchrotron-based measurements are strongly biased by the tail of the distribution here too. Finally, the same instability also breaks the common assumption that matter is advected at the post-shock velocity downstream, thus invalidating laminar-flow based cooling models.

We present observations of the gamma-ray binary PSR B1259-63/LS 2883 with the Atacama Large Millimeter/submillimeter Array (ALMA) at Bands 3 (97 GHz), 6 (233 GHz), and 7 (343 GHz). PSR B1259-63/LS 2883 consists of a pulsar in a highly eccentric orbit around a massive companion star, with the pulsar passing through the circumstellar disk near periastron. Our new data were obtained over several epochs, ranging from -61 to +29 days from the periastron passage in 2024. We report an increase in flux in all bands near the periastron. The significant change in Band 3 flux suggests synchrotron emission from the interaction between the pulsar wind and the stellar wind or disk. The Band 6 flux shows an increase around periastron and a transition from thermal emission from the circumstellar disk to synchrotron emission. The Band 7 observation +24 days after periastron shows a brightening, suggesting that the pulsar's passage through the disk does not result in its immediate destruction. We discuss the implications of these results for the interaction between the pulsar wind and the circumstellar disk, such as the possible disk expansion after periastron.

The interaction between radio-jets and quasar host galaxies plays a paramount role in quasar/galaxy co-evolution. However, very little has been known so far about this interaction at very high-z. Here, we present new Atacama Large Millimeter/submillimeter Array (ALMA) observations in Band 7 and Band 3 of six radio-loud quasars' host galaxies at $z > 5$. We recover [CII] 158 $\mu$m line and underlying dust continuum emission at $>2\sigma$ for five sources, while we obtain upper limits for the CO(6-5) emission line and continuum for the remaining source. At the spatial resolution of our observations ($\sim$1.0"-1.4"), we do not recover perturbed/extended morphologies or kinematics, signatures of potential mergers. These galaxies already host large quantities of gas, with [CII]-based star formation rates of $30-400 M_{\odot} $yr$^{-1}$. Building their radio/sub-mm spectral energy distributions (SEDs), we find that in at least four cases the 1mm continuum intensity arises from a combination of synchrotron and dust emission, with an initial estimation of synchrotron contribution at 300 GHz of $\gtrsim$10%. We compare the properties of the sources inspected here with a large collection of radio-quiet sources from the literature, as well as a sample of radio-loud quasars from previous studies, at comparable redshift. We recover a potential mild decrease in $L_{\rm [CII]}$ for the radio-loud sources, which might be due to a suppression of the cool gas emission due to the radio-jets. We do not find any [CII]-emitting companion galaxy candidate around the five radio-loud quasars observed in Band 7: given the depth of our dataset, this result is still consistent with that observed around radio-quiet quasars. Further higher-spatial resolution observations, over a larger frequency range, of high-z radio-loud quasars hosts will allow for a better understanding of the physics of such sources.

The James Webb Space Telescope (JWST) is opening new frontiers of transient discovery and follow-up at high-redshift. Here we present the discovery of a spectroscopically confirmed Type Ia supernova (SN Ia; SN 2023aeax) at z=2.15 with JWST, with cadenced NIRCam observations that enable multi-band light curve fitting. SN 2023aeax lands at the edge of traditional low-z cosmology color cuts because of its blue color (peak rest-frame B-V~-0.3), but we still apply a fiducial standardization approach with the BayeSN model and find that the SN 2023aeax luminosity distance measurement is in agreement (~0.1sigma) with LambdaCDM. SN 2023aeax is only the second spectroscopically confirmed SN Ia in the dark matter-dominated Universe at z>2 (the other is SN 2023adsy), giving it rare leverage to constrain any potential evolution in SN Ia standardized luminosities. Similar to SN 2023adsy (B-V~0.8), SN 2023aeax has a fairly extreme (but opposite) color, which may be due to the small sample size or a secondary factor, such as host galaxy properties. Nevertheless, the SN 2023aeax spectrum is well-represented by normal low-z SN Ia spectra and we find no definitive evolution in SN Ia standardization with redshift. Still, the first two spectroscopically confirmed z>2 SNe Ia have peculiar colors and combine for a ~1sigma distance slope relative to LambdaCDM, the implications of which require a larger sample and dedicated host galaxy observations to investigate.

Supermassive binary black holes are a key target for the future Laser Interferometer Space Antenna, and excellent multi-messenger sources with gravitational waves. However, unique features of their electromagnetic emission that are needed to distinguish them from single supermassive black holes are still being established. Here, we conduct the first magnetohydrodynamic simulation of accretion onto eccentric binary black holes in full general relativity incorporating synchrotron radiation transport through their dual-jet. We show that the total accretion rate, jet Poynting luminosity, and the optically thin synchrotron emission exhibit periodicity on the binary orbital period, demonstrating explicitly, for the first time, that the binary accretion rate periodicity can be reflected in its electromagnetic signatures. Additionally, we demonstrate that during each periodic cycle eccentric binaries spend more time in a low emission state than in a high state. Furthermore, we find that the gravitational wave bursts from eccentric binaries are coincident with the bursts in their jet luminosity and synchrotron emission. We discuss how multimessenger observations of these systems can probe plasma physics in their jet.

Protoclusters are high-$z$ overdense regions that will evolve into clusters of galaxies by $z=0$, making them ideal for studying galaxy evolution expected to be accelerated by environmental effects. However, it has been challenging to identify protoclusters beyond $z=3$ only by photometry due to large redshift uncertainties, hindering statistical study. To tackle the issue, we develop a new deep-learning-based protocluster detection model, PCFNet, which considers a protocluster as a point cloud. To detect protoclusters at $z\sim4$ using only optical broad-band photometry, we train and evaluate PCFNet with mock $g$-dropout galaxies based on the N-body simulation with the semi-analytic model. We use the sky distribution, $i$-band magnitude, $(g-i)$ color, and the redshift probability density function surrounding a target galaxy on the sky. PCFNet achieves to detect five times more protocluster member candidates while maintaining high purity (recall $=7.5\pm0.2$%, precision $=44\pm1$%) than conventional methods. Moreover, PCFNet is able to detect more progenitors ($M_\mathrm{halo}^{z=0}=10^{14-14.5}\,M_\odot$) that are less massive than supermassive clusters like the Coma cluster. We apply PCFNet to the observational photometric dataset of the HSC-SSP Deep/UltraDeep layer ($\sim17\mathrm{\,deg^2}$) and detect $121$ protocluster candidates at $z\sim4$. We find the rest-UV luminosities of our protocluster member candidates are brighter than those of field galaxies, which is consistent with previous studies. Additionally, the quenching of satellite galaxies depends on both the core galaxy's halo mass at $z\sim4$ and accumulated mass until $z=0$ in the simulation. PCFNet is very flexible and can find protoclusters at other redshifts or in future extensive surveys by Euclid, LSST, and Roman.

We present a comprehensive analysis of HI absorption around 96 lya emitters (LAEs) at $z\approx3.3$ (median lya luminosity $\approx10^{42}$ erg.s$^{-1}$). These LAEs were identified within 8 MUSE fields, each $1'\times1'$ on the sky and centered on a bright background quasar, as part of the MUSEQuBES survey. Using Voigt profile fitting for all HI absorbers detected within $\pm500$ km.$s^{-1}$ of these LAEs, we compiled a catalog of 800 HI absorption components. Our analysis shows that HI absorption is enhanced near the LAEs compared to the IGM. However, no trend is found between the column densities of HI absorbers and their impact parameters from the LAEs (spanning $\approx54$ to 260 pkpc). Additionally, all galaxies associated with Lyman-limit systems have impact parameters $>50$ pkpc from the quasar sightlines, suggesting that true absorber-hosts may be too faint to detect. The LAEs show an overall HI covering fraction (fc(HI)) of $\approx88\%$ for a threshold logN(HI)$=15$. Notably, at the same threshold, the pairs/group LAEs exhibit a $100\%$ HI covering fraction out to $\approx 250$ pkpc. In contrast, isolated LAEs consistently show a lower fc(HI) of $\approx80\%$. This environmental influence on fc(HI) is also evident up to $\approx 300$ km.$s^{-1}$ in differential bins of line-of-sight velocity. We find an anti-correlation between fc(HI) and the rest-frame lya-emission equivalent width (ew). Based on the lya-shell model, this could imply that gas-rich galaxies tend to reside in gas-rich environments or that the higher EW LAEs are more efficient at ionizing their surrounding medium.

This study presents the modeling of the gravitational wave (GW) bias parameter by bridging a connection between simulated GW sources and galaxies in low redshift galaxy surveys 2MPZ and WISExSCOS (WISC). We study this connection by creating a mock GW catalog, populating galaxy surveys with binary black holes (BBHs) for different scenarios of the GW occupation fraction (or selection function) as a function of the galaxy stellar mass. We probe the observable consequences of this connection by exploring the spatial clustering of the GW sources in terms of the GW bias parameter. We consider a phenomenological broken power law model for the selection function, with a potential turnover $M_{K}$ at high stellar mass ($10^{11}$ $M_{\odot}$ in the fiducial model) where the star formation efficiency begins to drop. We vary the parameters of the selection function and find that generically the GW bias increases as $M_{K}$ increases (and gets suppressed as $M_{K}$ decreases). The change in the GW bias parameter shows a maximum change of about $30\%$ for different scenarios explored in this work in comparison to the galaxy bias. Future measurements of the GW bias can help constrain $M_{K}$ and the slopes of the selection function and thus offer insights into the underlying astrophysical processes.

In this work, we measure the spectral index of Sagittarius A* (Sgr A*) between the $H$ (1.6 $\mu$m) and $K^\prime$ (2.2 $\mu$m) broadband filters in the near-infrared (NIR), sampling over a factor $\sim 40$ in brightness, the largest range probed to date by a factor $\sim 3$. Sgr A*-NIR is highly variable, and studying the spectral index $\alpha$ (with $F_\nu \propto \nu^{\alpha}$) is essential to determine the underlying emission mechanism. For example, variations in $\alpha$ with flux may arise from shifts in the synchrotron cutoff frequency, changes in the distribution of electrons, or multiple concurrent emission mechanisms. We investigate potential variations of $\alpha_{H-K^\prime}$ with flux by analyzing 7 epochs (2005 to 2022) of Keck Observatory imaging observations from the Galactic Center Orbits Initiative (GCOI). We remove the flux contribution of known sources confused with Sgr A*-NIR, which can significantly impact color at faint flux levels. We interpolate between the interleaved $H$ and $K^\prime$ observations using Multi-Output Gaussian Processes. We introduce a flexible empirical model to quantify $\alpha$ variations and probe different scenarios. The observations are best fit by an $\alpha_{H-K^\prime} = - 0.50 \pm 0.08 _{\rm stat} \pm 0.17_{\rm sys}$ that is constant from $\sim 1$ mJy to $\sim 40$ mJy (dereddened 2 $\mu$m flux). We find no evidence for a flux-dependence of Sgr A*'s intrinsic spectral index. In particular, we rule out a model explaining NIR variability purely by shifts in the synchrotron cutoff frequency. We also constrain the presence of redder, quiescent emission from the black hole, concluding that the dereddened 2 $\mu$m flux contribution must be $\leq 0.3$ mJy at 95% confidence level.

Recent literature on core-collapse supernovae suggests that a black hole (BH) can form within $\sim 1$ s of shock revival, while still culminating in a successful supernova. We refer to these as black hole supernovae, as they are distinct from other BH formation channels in both timescale and impact on the explosion. We simulate these events self-consistently from core-collapse until $20\text{-}50$ days after collapse using three axisymmetric models of a $60$ M$_\odot$ zero-age main sequence progenitor star and investigate how the composition of the ejecta is impacted by the BH formation. We employ Skyrme-type equations of state (EOSs) and vary the uncertain nucleonic effective mass, which affects the pressure inside the proto-neutron star through the thermal part of the EOS. This results in different BH formation times and explosion energies at BH formation, yielding final explosion energies between $0.06\text{-}0.72\times 10^{51}$ erg with $21.8\text{-}23.3$ M$_\odot$ of ejecta, of which $0\text{-}0.018$ M$_\odot$ is $^{56}$Ni. Compared to expectations from 1D simulations, we find a more nuanced EOS dependence of the explosion dynamics, the mass of the BH remnant, and the elemental composition of the ejecta. We investigate why the explosions survive despite the massive overburden and link the shape of the diagnostic energy curve and character of the ejecta evolution to the progenitor structure.

We present a novel technique using Fourier series and Laguerre polynomials to represent morphological features of disc galaxies. To demonstrate the utility of this technique, we study the evolution of asymmetry in a sample of disc galaxies drawn from the Extended Groth Strip and imaged by the JWST Cosmic Evolution Early Release Science Survey as well as archival HST observations. We measure disc asymmetry as the amplitude of the of the m = 1 Fourier harmonic for galaxies within redshift ranges of 1 < z < 4. We show that when viewed in shorter rest frame wavelengths, disc galaxies have a higher asymmetry as the flux is dominated by star forming regions. We find generally low asymmetry at rest frame infrared wavelengths, where our metric tracks asymmetry in morphological features such as bars and spiral arms. We show that higher mass galaxies have lower asymmetry and vice versa. Higher asymmetry in lower mass galaxies comes from lower mass galaxies (typically) having higher star formation rates. We measure the relation between disc galaxy asymmetry and redshift and find no conclusive relationship between them. We demonstrate the utility of the Fourier-Laguerre technique for recovering physically informative asymmetry measurements as compared to rotational asymmetry measurements. We also release the software pipeline and quantitative analysis for each galaxy.

With the advent of current and future high-resolution CMB experiments, the kinematic Sunyaev-Zel'dovich (kSZ) effect has become a unique observational probe of the distribution of baryons and velocities in the Universe. In this work, we propose a novel binned bispectrum of the form temperature-temperature-density to extract the late-time kSZ effect from cleaned CMB maps. Unlike 'kSZ tomography' methods, this estimator can use any tracer of the large-scale structure density field projected along the line-of-sight and does not require individual redshifts. With our method, we forecast signal-to-noise ratios (SNR) of $\sim$100-200 for the upcoming Simons Observatory (SO) and CMB-S4 correlated with a galaxy sample from WISE that is restricted to the linear regime. We also extend galaxy modes into the non-linear regime and explore this harmonic space to show that the SNR peaks for squeezed triangles that have a short (linear) density mode and long temperature modes in harmonic space. The existing kSZ$^{2}$-density projected-fields estimator compresses the rich information contained in this bispectrum across various scales and triangle shapes. Moreover, we find that the lensing correction to our estimator's signal is relatively small. We study the dependence of this kSZ signal on $\Lambda$CDM parameters for SO and CMB-S4 and forecast initial constraints on the sum of neutrino masses while restricting to the linear galaxy bias regime. Our work illustrates the potential of the projected-fields kSZ bispectrum as a novel probe of baryonic abundance and beyond-$\Lambda$CDM cosmology with upcoming precision measurements.

The second data release of Type Ia supernovae (SNe Ia) observed by the Zwicky Transient Facility has provided a homogeneous sample of 3628 SNe Ia with photometric and spectral information. This unprecedented sample size enables us to better explore our currently tentative understanding of the dependence of host environment on SN Ia properties. In this paper, we make use of two-dimensional image decomposition to model the host galaxies of SNe Ia. We model elliptical galaxies as well as disk/spiral galaxies with or without central bulges and bars. This allows for the categorisation of SN Ia based on their morphological host environment, as well as the extraction of intrinsic galaxy properties corrected for both cosmological and atmospheric effects. We find that although this image decomposition technique leads to a significant bias towards elliptical galaxies in our final sample of galaxies, the overall results are robust. By successfully modelling 728 host galaxies, we find that the photometric properties of SNe Ia found in disks and in elliptical galaxies, correlate fundamentally differently with their host environment. We identified strong linear relations between light-curve stretch and our model-derived galaxy colour for both the elliptical (16.8$\sigma$) and disk (5.1$\sigma$) subpopulations of SNe Ia. Lower stretch SNe Ia are found in redder environments, which we identify as an age/metallicity effect. Within the subpopulation of SNe Ia found in disk containing galaxies, we find a significant linear trend (6.1$\sigma$) between light-curve stretch and model-derived local $r$-band surface brightness, which we link to the age/metallicity gradients found in disk galaxies. SN Ia colour shows little correlation with host environment as seen in the literature. We identify a possible dust effect in our model-derived surface brightness (3.3$\sigma$), for SNe Ia in disk galaxies.

Space-based astronomy of hard X-rays and gamma rays covers more than seven orders of magnitude in photon energy, from 10 keV to several hundred GeV. Detecting cosmic photons in this energy range is a challenge, due to the relatively low probability of interaction of high-energy photons with matter and the high background noise generated in space detectors by environmental charged particles and radiation. However, the development of new detection technologies is constantly improving the performance of space-based X- and gamma-ray telescopes. This chapter presents the different detectors used in this field of astronomy, their configuration within space telescopes and some proposals for new instruments.

We present the R-band lightcurves of the Flora family asteroid (12499) 1998 FR47, obtained in 2022 at two different astronomical sites: Bulgarian National Astronomical Observatory Rozhen (MPC Code 071) and Astronomical Station Vidojevica (MPC Code C89). The quadramodal lightcurves reveal a rotation period of 6.172+/-0.003 h and an amplitude of about 0.44 mag. Using the lightcurve inversion method, with the combination of our dense lightcurves and the sparse data from Gaia DR3, we found the sidereal period, an indication of a retrograde rotation of (12499) and its low-convex resolution shape. Nonetheless, the unusual shape of the quadramodal lightcurve and its additional analysis reveals two possible periods, 3.0834+/-0.0085 h and 4.1245+/-0.0151 h, making the suspect that the asteroid might be a non-synchronised wide binary system. Spectral analysis of the asteroid using data from the GAIA DR3 shows that it is either an M- or an L-type object and maybe a piece of the first planetesimals to form in the Solar System protoplanetary disk. On the other hand, (12499)'s dynamical properties indicate a significantly shorter lifetime. The asteroid lies exactly on the chaotic border of the 7:4 mean motion resonance with Mars (7M:4), alternating between being in and out of it for almost 190 Myrs. During 200 Myrs of integration, (12499) visited other resonances in the Flora family, but it never became a Near Earth Object (NEO). Additional integration of fictive objects from the 7M:4 resonance showed a possibility of transportation to the NEO region already at about 20 Myrs.

Light echoes offer a means of studying protoplanetary disks, including their geometry and composition, even when they are not spatially resolved. We present a test of this approach applied specifically to optically thick, geometrically flared disks around active stars. Here we adopt stellar parameters of an active M dwarf to calculate light echoes for disks and rings with radii that would produce time delays consistent with TESS short cadence (about 2 minutes) time bins. Our results show successful fits to disk parameters, highlighting the potential effectiveness of this method in the search for protoplanetary disks.

JWST has revealed a large population of compact, red galaxies at $z>4$ known as Little Red Dots (LRDs). We analyze the spectral energy distributions (SEDs) of 95 LRDs from the JWST PRIMER survey with complete photometric coverage from $1-18\ \mu$m using NIRCam and MIRI imaging, representing the most extensive SED analysis on a large LRD sample with long-wavelength MIRI data. We examine SED models in which either galaxy or active galactic nucleus (AGN) emission dominates the rest-frame UV or optical continuum, extracting physical properties to explore each scenario's implications. In the galaxy-only model, we find massive, dusty stellar populations alongside unobscured, low-mass components, hinting at inhomogeneous obscuration. The AGN-only model indicates dusty, luminous AGNs with low hot dust fractions compared to typical quasars. A hybrid AGN and galaxy model suggests low-mass, unobscured galaxies in the UV, with stellar mass estimates spanning $\sim$2 dex across the different models, underscoring the need for caution in interpreting LRD stellar masses. With MIRI photometry, the galaxy-only model produces stellar masses within cosmological limits, but extremely high stellar mass densities are inferred. The black hole and stellar masses inferred from the hybrid model lead to highly overmassive black holes even exceeding those in recently reported high-redshift AGNs, hinting at a partial AGN contribution to the rest-optical continuum or widespread super-Eddington accretion. Our findings highlight the extreme conditions required for both AGN or galaxy dominated scenarios in LRDs, supporting a mixed contribution to the red continuum, or novel scenarios to explain the observed emission.

Dark matter is theorised to form massive haloes, which could be further condensed into so-called spikes when a black hole grows at the centre of such a halo. The existence of these spikes is instrumental for several dark matter detection schemes such as indirect detection and imprints on gravitational wave inspirals, but all previous work on their formation has been (semi-)analytical. We present fully numerically simulated cold dark matter spikes using the SWIFT code. Based on these results, we propose a simple empirical density profile - dependent on only a single mass-ratio parameter between the black hole and total mass - for dark matter spikes grown in Hernquist profiles. We find that the radius of the spike scales differently compared to theoretical predictions, and show a depletion of the outer halo that is significant for high mass-ratio systems. We critically assess approximations of the spike as used in the field, show that our profile significantly deviates, and contextualise the potential influence for future dark matter detections by simulating binary black hole inspirals embedded in our profile.

The estimation of uncertainties in cosmological parameters is an important challenge in Large-Scale-Structure (LSS) analyses. For standard analyses such as Baryon Acoustic Oscillations (BAO) and Full Shape, two approaches are usually considered. First: analytical estimates of the covariance matrix use Gaussian approximations and (nonlinear) clustering measurements to estimate the matrix, which allows a relatively fast and computationally cheap way to generate matrices that adapt to an arbitrary clustering measurement. On the other hand, sample covariances are an empirical estimate of the matrix based on en ensemble of clustering measurements from fast and approximate simulations. While more computationally expensive due to the large amount of simulations and volume required, these allow us to take into account systematics that are impossible to model analytically. In this work we compare these two approaches in order to enable DESI's key analyses. We find that the configuration space analytical estimate performs satisfactorily in BAO analyses and its flexibility in terms of input clustering makes it the fiducial choice for DESI's 2024 BAO analysis. On the contrary, the analytical computation of the covariance matrix in Fourier space does not reproduce the expected measurements in terms of Full Shape analyses, which motivates the use of a corrected mock covariance for DESI's Full Shape analysis.

Clouds and hazes are abundant in the thin and cold atmospheres of Triton and Pluto, where they are thought to be produced by interactions between atmospheric gases and ultraviolet photons from the Sun and those scattered by the local interstellar medium. These interactions lead to a rich network of chemical reactions that produces higher order hydrocarbons and nitriles that condense out to form ice clouds, and ultimately complex haze particles that rain down onto the surface that impact the atmospheric thermal structure, gas chemistry, and surface evolution. In this chapter, we will review the observational evidence for clouds and hazes in the atmospheres of Triton and Pluto and theoretical interpretations thereof, and the emerging set of experiments aiming to produce Triton and Pluto clouds and hazes in the lab to learn about them in detail.

This study investigates the structural properties of strange quark stars (SQS) using a Quantum Chromodynamics (QCD) perturbative model combined with the latest Particle Data Group dataset. Given the energy scale present in compact stars, QCD perturbation theory alone may not fully explain their structure. To account for non-perturbative contributions, we incorporate a density-dependent effective bag parameter, $B$, and derive the equation of state (EOS) for strange quark matter (SQM). We start by demonstrating the limitations of EOSs with a constant $B$ in describing massive objects with $ M_{TOV}> 2M_{\odot} $. Subsequently, we show that considering $B$ as a density-dependent function significantly changes the results. Our definition of $B$ includes two parameters determined by both theoretical and observational constraints. We demonstrate that incorporating a density-dependent $B$ into the perturbative EOS can yield SQSs with masses exceeding $2M_{\odot}$, while complying with gravitational wave constraints such as tidal deformability, and thermodynamic considerations, including stability conditions and speed of sound behavior. Specifically, we show that massive compact objects like PSR J0952-0607, PSR J2215+5135, PSR J0740+6620, and the secondary mass of GW190814 can be SQSs. Additionally, we compare our EOS with the EOS of the authors who use a generalized polytropic form with adjustable parameters and obtain an interesting result.

Gamma-ray bursts (GRBs) are the most energetic explosions in the universe, and their afterglow emission provides an opportunity to probe the physics of relativistic shock waves in an extreme environment. Several key pieces for completing the picture of the GRB afterglow physics are still missing, including jet properties, emission mechanism, and particle acceleration. Here we present a study of the afterglow emission of GRB 221009A, the most energetic GRB ever observed. Using optical, X-ray, and gamma-ray data up to approximately two days after the trigger, we trace the evolution of the multi-wavelength spectrum and the physical parameters behind the emission process. The broadband spectrum is consistent with the synchrotron emission emitted by relativistic electrons with its index of $p = 2.29\pm 0.02$. We identify a break energy at keV and an exponential cutoff at GeV in the observed multi-wavelength spectrum. The break energy increases in time from $16.0_{-4.9}^{+7.1}$ keV at 0.65 days to $46.8_{-15.5}^{+25.0}$ keV at 1.68 days, favoring a stellar wind-like profile of the circumburst medium with $k=2.4\pm0.1$ as in $\rho (r) \propto r^{-k}$. The high-energy attenuation at around 0.4 to 4 GeV is attributed to the maximum of the particle acceleration in the relativistic shock wave. This study confirms that the synchrotron process can explain the multi-wavelength afterglow emission and its evolution.

A consequence of a non-zero occupation fraction of massive black holes (MBHs) in dwarf galaxies is that these MBHs can become residents of larger galaxy halos via hierarchical merging and tidal stripping. Depending on the parameters of their orbits and original hosts, some of these MBHs will merge with the central supermassive black hole in the larger galaxy. We examine four cosmological zoom-in simulations of Milky Way-like galaxies to study the demographics of the black hole mergers which originate from dwarf galaxies. Approximately half of these mergers have mass ratios less than 0.04, which we categorize as intermediate mass ratio inspirals, or IMRIs. Inspiral durations range from 0.5 - 8 Gyr, depending on the compactness of the dwarf galaxy. Approximately half of the inspirals may become more circular with time, while the eccentricity of the remainder does not evolve. Overall, IMRIs in Milky Way-like galaxies are a significant class of black hole merger that can be detected by LISA, and must be prioritized for waveform modeling.

This paper presents the first analysis of the contact binary TYC 3801-1529-1. We observed four sets of multiple bands complete light curves and one set of radial velocity curve of the primary component. Based on a simultaneous investigation of our observed and TESS light curves and the radial velocity curve, we found that TYC 3801-1529-1 is an extremely low-mass-ratio, medium contact binary with $q=0.0356$, with the contribution of the third light at a level of about 10\%. Its mass ratio is lower than V1187 Her, making TYC 3801-1529-1 the lowest mass-ratio contact binary ever found in the universe. The light curves observed in 2022 are asymmetric, which is aptly explained by a hot spot on the primary component. A 16-year eclipse timings analysis indicates a secular increase orbital period with a rate of dp/dt$=7.96(\pm0.35)\times10^{-7}$ d yr$^{-1}$. We studied the stability of this target and identified that not only the value of $J_{spin}/J_{orb}$, but also the mass ratio surpass the unstable boundary. Hence, TYC 3801-1529-1 presents a challenge to theoretical research and ought to be considered a progenitor of a contact binary merger.

A planetary growth rate, a.k.a., the mass accretion rate, is a fundamental parameter in planet formation, as it determines a planet's final mass. Planetary mass accretion rates have been estimated using hydrogen lines, based on the models originally developed for accreting stars, known as the accretion flow model. Recently, Aoyama et al. (2018) introduced the accretion shock model as an alternative mechanism for hydrogen line emission. However, it remains unclear which model is more appropriate for accreting planets and substellar objects. To address this, we applied both models to archival data consisting of 96 data points from 76 accreting brown dwarfs and very low-mass stars, with masses ranging from approximately 0.02 to 0.1 $M_\sun$, to test which model best explains their accreting properties. The results showed that the emission mechanisms of 15 data points are best explained by the shock model, while 55 data points are best explained by the flow model. For the 15 data points explained by the planetary shock model, the shock model estimates up to several times higher mass accretion rates than the flow model. As this trend is more pronounced for planetary mass objects, it is crucial to determine which emission mechanism is dominant in individual planets. We also discuss the physical parameters that determine the emission mechanisms and the variability of line ratios.

We obtain CO(1-0) molecular gas measurements with ATCA on a sample of 43 spectroscopically confirmed H$\alpha$ emitters in the Spiderweb protocluster at $z=2.16$ and investigate the relation between their star formation and cold gas reservoirs as a function of environment. We achieve a CO(1-0) detection rate of $\sim23\pm12\%$ with 10 dual CO(1-0) and H$\alpha$ detections at $10<\log M_{*}/M_\odot<11.5$. In addition, we obtain upper limits for the remaining sources. In terms of total gas fractions ($F_{gas}$), our sample is divided into two different regimes with a steep transition at $\log M_{*}/M_\odot\approx10.5$. Galaxies below that threshold have gas fractions that in some cases are close to unity, indicating that their gas reservoir has been replenished by inflows from the cosmic web. However, objects at $\log M_{*}/M_\odot>10.5$ display significantly lower gas fractions and are dominated by AGN (12 out of 20). Stacking results yield $F_{gas}\approx0.55$ for massive emitters excluding AGN, and $F_{gas}\approx0.35$ when examining only AGN candidates. Furthermore, depletion times show that most H$\alpha$ emitters may become passive by $1<z<1.6$, concurrently with the surge and dominance of the red sequence in the most massive clusters. Our analyses suggest that galaxies in the outskirts of the protocluster have larger molecular-to-stellar mass ratios and lower star formation efficiencies than in the core. However, star formation across the protocluster remains consistent with the main sequence, indicating that evolution is primarily driven by the depletion of the gas reservoir towards the inner regions. We discuss the relative importance of in-/outflow processes in regulating star formation during the early phases of cluster assembly and conclude that a combination of feedback and overconsumption may be responsible for the rapid cold gas depletion these objects endure.

An updated analysis of pulsar glitch parameters was conducted using a sample of 215 pulsars, encompassing 677 recorded glitch events, each glitch at least once between 1968 and 2024. The glitch rates were estimated and plotted against various pulsar parameters, including rotational frequency, spin-down rate, and characteristic age. These relationships were analyzed using Pearson's correlation coefficient. The results indicate linear and weak relationships between glitch rate and the pulsar parameters. However, the relationship between glitch rate and characteristic age is inversely weak.

Observational studies are identifying stars thought to be remnants from the earliest stages of the hierarchical mass assembly of the Milky Way, referred to as the proto-Galaxy. We use red giant stars with kinematics and [$\alpha$/M] and [M/H] estimates from Gaia DR3 data to investigate the relationship between azimuthal velocity and metallicity, aiming to understand the transition from a chaotic proto-Galaxy to a well-ordered, rotating (old) disc-like population. To analyse the structure of the data in [M/H]-v$_\phi$ space for both high- and low-$\alpha$ samples with carefully defined $\alpha$-separation, we develop a model with two Gaussian components in v$_\phi$: one representing a disc-like population and the other a halo-like population. This model is designed to capture the conditional distribution P(v$_\phi$ $\mid$ [M/H]) with a 2-component Gaussian Mixture Model with fixed azimuthal velocities means and standard deviations. To quantify the spin-up of the high-$\alpha$ disc population, we extend this two-component model by allowing the mean velocity and velocity dispersion to vary between the spline knots across the metallicity range used. We also compare our findings with existing literature using traditional Gaussian Mixture Modelling in bins of [M/H] and investigate using orbital circularity instead of azimuthal velocity. Our findings show that the metal-poor high-$\alpha$ disc gradually spins up across [M/H] $\sim$ -1.7 to -1.0, while the low-$\alpha$ sample exhibits a sharp transition at [M/H] $\sim$ -1.0. This latter result is due to the accreted debris dominating the metal-poor end, underscoring the critical role of [$\alpha$/M] selection in studying the (old) disc evolution of the Milky Way. These results indicate that the proto-Galaxy underwent a slow, monotonic spin-up phase rather than a rapid, dramatic spin-up at [M/H] $\sim$ -1.0, as previously inferred.

I present CMBAnalysis, a state-of-the-art Python framework designed for high-precision analysis of Cosmic Microwave Background (CMB) radiation data. This comprehensive package implements parallel Markov Chain Monte Carlo (MCMC) techniques for robust cosmological parameter estimation, featuring adaptive integration methods and sophisticated error propagation. The framework incorporates recent advances in computational cosmology, including support for extended cosmological models, detailed systematic error analysis, and optimized numerical algorithms. I demonstrate its capabilities through analysis of Planck Legacy Archive data, achieving parameter constraints competitive with established pipelines while offering significant performance improvements through parallel processing and algorithmic optimizations. Notable features include automated convergence diagnostics, comprehensive uncertainty quantification, and publication-quality visualization tools. The framework's modular architecture facilitates extension to new cosmological models and analysis techniques, while maintaining numerical stability through carefully implemented regularization schemes. My implementation achieves excellent computational efficiency, with parallel MCMC sampling reducing analysis time by up to 75\% compared to serial implementations. The code is open-source, extensively documented, and includes a comprehensive test suite, making it valuable for both research applications and educational purposes in modern cosmology.

We estimate halo spins for HI-rich galaxies in the Arecibo Legacy Fast Alfa Survey using a semi-analytic approach, examining the relationship between halo spin and stellar surface density. Our findings reveal an inverse correlation in both low- and high-mass galaxy samples, with stellar surface density decreasing as halo spin increases. This trend highlights the pivotal role of halo spin in galaxy evolution and suggests a universal formation scenario: high-spin halos, accompanied by high-spin accreted gas, retain angular momentum, preventing gas from efficiently condensing in the galactic center and thus suppressing star formation. Consequently, weak feedback redistributes gas to the halo outskirts without significant expulsion. The shallower central gravitational potential in high-spin halos promotes outward stellar migration, leading to more extended stellar distributions and lower stellar surface densities.

Leveraging the semi-analytic method, we compute halo spins for a substantial sample of HI-bearing galaxies observed in the Arecibo Legacy Fast Alfa Survey. Our statistical analysis reveals a correlation between halo spin and environment, although the trend is subtle. On average, galaxies exhibit a decreasing halo spin tendency in denser environments. This observation contrasts with previous results from $N$-body simulations in the Lambda cold dark matter framework. The discrepancy may be attributed to environmental gas stripping, leading to an underestimation of halo spins in galaxies in denser environments, or to baryonic processes that significantly alter the original dark matter halo spins, deviating from previous $N$-body simulation findings.

Utilizing ALFALFA HI data, we investigate the relationship between specific star formation rate (sSFR) and halo spin across various star-forming galaxies. Our analysis reveals no significant correlation between sSFR and halo spin, irrespective of the galactic environment. Previous research suggests that high-spin halos tend to harbor extended, low-density stellar distributions due to suppressed gas cooling and star formation. However, unlike galaxy size and density, sSFR may primarily reflect the current star-forming state rather than long-term history, indicating potential independence from halo spin.

Interferometric closure invariants, constructed from triangular loops of mixed Fourier components, capture calibration-independent information on source morphology. While a complete set of closure invariants is directly obtainable from measured visibilities, the inverse transformation from closure invariants to the source intensity distribution is not established. In this work, we demonstrate a deep learning approach, Deep learning Image Reconstruction with Closure Terms (DIReCT), to directly reconstruct the image from closure invariants. Trained on both well-defined mathematical shapes (two-dimensional gaussians, disks, ellipses, $m$-rings) and natural images (CIFAR-10), the results from our specially designed model are insensitive to station-based corruptions and thermal noise. The median fidelity score between the reconstruction and the blurred ground truth achieved is $\gtrsim 0.9$ even for untrained morphologies, where a unit score denotes perfect reconstruction. In our validation tests, DIReCT's results are comparable to other state-of-the-art deconvolution and regularised maximum-likelihood image reconstruction algorithms, with the advantage that DIReCT does not require hand-tuned hyperparameters for each individual prediction. This independent approach shows promising results and offers a calibration-independent constraint on source morphology, ultimately complementing and improving the reliability of sparse VLBI imaging results.

In this manuscript, an oversimplified model is proposed for the first time to explain the different variability trends in the observed broad H$\alpha$ emission line luminosity $L_{H\alpha}(t)$ and in the TDE model determined bolometric luminosity $L_{bol}(t)$ in the known TDE ASASSN-14li. Considering broad emission line regions (BLRs) lying into central accretion disk related to accreted materials onto central black hole in TDEs, mass evolution $M_{BLRs}(t)$ of central BLRs can be determined by the maximum mass $M_{BLRs,0}$ of central BLRs minus the corresponding accreted mass in a TDE. Meanwhile, through the simple linear dependence of Broad Balmer emission line luminosity on mass of BLRs, the mass evolution $M_{BLRs}(t)$ of central BLRs can be applied to describe the observed $L_{H\alpha}(t)$. Although the proposed model is oversimplified with only one free model parameter $M_{BLRs,0}$, the model with $M_{BLRs,0}\sim0.02{\rm M_\odot}$ can be applied to well describe the observed $L_{H\alpha}(t)$ in the TDE ASASSN-14li. Meanwhile, the oversimplified model can also be applied to roughly describe the observed $L_{H\alpha}(t)$ in the TDE ASASSN-14ae. The reasonable descriptions to the observed $L_{H\alpha}(t)$ in ASASSN-14li and ASASSN-14ae indicate the oversimplified model with only one free model parameter is probably efficient enough to describe mass evolutions of $M_{BLRs}$ related to central accreted debris in TDEs.

We present a reconstruction of the line-of-sight motions of the local interstellar medium (ISM) based on the combination of a state-of-the-art model of the three-dimensional dust density distribution within 1.25 kpc from the Sun and the HI and CO line emission within Galactic latitudes $|b| < 5^{\circ}$. We use the histogram of oriented gradient (HOG) method to match the plane-of-the-sky 3D dust distribution across distances with the atomic and molecular line emission. We identify a significant correlation between the 3D dust model and the line emission. We employ this correlation to assign line-of-sight velocities to the dust across density channels and produce a face-on map of the local ISM radial motions with respect to the local standard of rest (LSR). We find that most of the material in the 3D dust model follows the large-scale pattern of Galactic rotation; however, we also report local departures from the rotation pattern with standard deviations of 12.1 and 6.1 km/s for the HI and CO line emission, respectively. The mean kinetic energy densities of these streaming motions are around 0.68 and 0.18 eV/cm$^{3}$ from either gas tracer. Assuming homogeneity and isotropy in the velocity field, these values are within a factor of a few of the total kinetic energy density. These kinetic energy values are roughly comparable to other energy densities, thus confirming the near-equipartition introduced by the feedback loops connecting the physical processes in the local ISM. Yet, we find energy and momentum overdensities of around a factor of ten concentrated in the Radcliffe Wave, the Split, and other local density structures. Although we do not find evidence of the Local Spiral Arm impact in these overdensities, their distribution suggests the influence of other large-scale effects that, in addition to supernova feedback, shape the energy distribution in the Solar neighborhood.

Characterising stellar jet asymmetries is key to setting robust constraints on jet launching models and improving our understanding of the underlying mechanisms behind jet launching. We aim to characterise the asymmetric properties of the bipolar jet coming from the Classical T Tauri Star Th 28. We combined data from integral field spectroscopy with VLT/MUSE and high-resolution spectra from VLT/X-shooter to map the optical emission line ratios in both jet lobes. We carried out a diagnostic analysis of these ratios to compare the density, electron temperature, and ionisation fraction within both lobes. The mass accretion rate was derived from the emission lines at the source and compared with the mass outflow rate derived for both lobes, using the estimated densities and measured [O I]6300 and [S II ]6731 luminosities. The blue-shifted jet exhibits a significantly higher electron temperature and moderately higher ionisation fraction than the red-shifted jet. In contrast to previous studies, we also estimated higher densities in the blue-shifted jet by a factor of ~2. These asymmetries are traced to within 160 au of the source in the line ratio maps. We estimate the mass accretion rate onto the central star and compare this with estimates of the mass outflow rates through each side of the jet. The emission line maps and diagnostic results suggest that the jet asymmetries originate close to the source and are likely to be intrinsic to the jet. Furthermore, the combined dataset offers access to a broad array of accretion tracers. In turn, this enables a more accurate estimation of the mass accretion rate, revealing a value of Macc that is higher by a factor >350 than would otherwise be determined. Supplementary figures and tables are available via a public Zenodo repository (doi:https://doi.org/10.5281/zenodo.13373809).

For high-energy cosmic-ray physics, it is imperative to determine the mass and energy of the cosmic ray that initiated the air shower in the atmosphere. This information can be extracted from the longitudinal profile of the air shower. In radio-metric observations, this profile is customarily determined through an extensive fitting procedure where calculated radio intensity is fitted to data. Beamforming the measured signals offers a promising alternative to bypass the cumbersome fitting procedure and to determine the longitudinal profile directly. Finite aperture effects in beamforming hamper the resolution with which this profile can be determined. We present a comprehensive investigation of the beamforming resolution in radiometric observations of air showers. There are two, principally different, approaches possible in air-shower beamforming, one where the total beamforming intensity is determined and an alternative where the beamforming trace is cross-correlated with a known response function. The effects due to a finite aperture (size of antenna array and bandwidth) are large for both approaches. We argue that it is possible to correct for the aperture corrections using an unfolding procedure. We give an explicit expression for the folding function, the kernel. Being able to calculate the folding function allows for unfolding the finite aperture effects from the data. We show that, in a model-to-model comparison, this allows for an accurate reconstruction of the current profile as the shower develops in the atmosphere. We present also an example where we reconstruct the longitudinal current profile of a shower developing under thunderstorm conditions where the atmospheric electric fields greatly alter the orientation of the transverse current in the shower front.

Context. HD140283, or the Methuselah star, is a well-known reference object in stellar evolution. Its peculiar chemical composition, proximity and absence of reddening makes it an interesting case-study of Pop II stars. Thanks to recent observational efforts, we now have precise interferometric and spectroscopic constraints, as well as revised astrometric parallaxes from the Gaia mission. Aims. We aim at determining the age of HD140283 with these lastest constraints, as well as quantifying the impact of systematics from physical inaccuracies in the stellar evolution models. Methods. Using recent spectroscopic abundances from the literature, including 3D non-LTE values for C, O, and Fe, we compute opacity tables specific to HD140283. We then use them in grids of stellar evolution models coupled to a Markov Chain Monte Carlo tool to determine the age of HD140283. Results. With our tailored models we find an age of 12.3Gy. Using a solar-scaled mixture instead results in an age value of 14Gy, in tension with the age of the universe ($13.77\pm0.06$Gy). We also find that reducing the mixing length parameter from its solar calibrated value will lead to an even lower age, in agreement with other recent studies. However, we find no direct evidence to favour a lower mixing length parameter value from our modelling. Conclusions. Taking into account the specific elemental abundances is crucial for the modelling of HD140283, as it leads to significant differences in the inferred age. However, this effect is degenerate with a lowering of the mixing length parameter. In this respect, asteroseismic constraints might play a key role in accurately deriving the mass of HD140283, therefore strongly constraining its age.

We detail a platform for partial g environment and an experiment for simulated impacts on asteroid surfaces based on it. The partial g environment is created by a two stage approach: First, create microgravity using the ZARM drop tower. Second, convert microgravity to partial gravity by steady acceleration of experiment volume on linear drive inside microgravity environment. The experiment we conducted on this platform simulates low-velocity impacts into a simulated asteroid surface. To recreate the asteroid environment, in addition to the partial gravity, a vacuum chamber is used. We explain requirements, setup and operation of partial gravity platform and experiment and discuss its performance. Finally, we are open for requests for external experiments which might benefit from our platform with $9.3\,$s of controlled partial gravity down to the mm/s$^2$ range with low g-jitter.

We investigate chemistry in the compression layer behind the interstellar shock waves, where molecular cloud formation starts. We perform three-dimensional magnetohydrodynamics simulations of converging flows of atomic gas with shock parameters of inclination between the interstellar magnetic field and the shock wave, pre-shock density, and shock velocity. Then we derive 1D mean-flow models, along which we calculate a detailed gas-grain chemical reaction network as a post process with various chemical parameters, i.e. cosmic-ray ionization rate, abundances of PAHs, and metals in the gas phase. While carbon chains reach their peak abundances when atomic carbon is dominant in the pseudo-time-dependent models of molecular clouds, such behavior is less significant in our models since the visual extinction of the compression layer is low ($\lesssim 1$ mag) when atomic carbon is abundant. Carbon chains, CN, and HCN increase at $A_V \gtrsim 1$ mag, where the gas-phase C/O ratio increases due to water ice formation. Shock parameters affect the physical structure and the evolutional timescale of the compression layer, and thus molecular evolution. Carbon chains are more abundant in models with higher post-shock density and slower gas accumulation. We calculate molecular column densities in the compression layer and compare them with the observations of diffuse and translucent clouds, which show reasonable agreement for water ice, carbon chains, and HCO$^+$. The observed variation of their column densities could be due to the difference in shock parameters and chemical parameters. The column density of CN is overestimated, for which we discuss possible reasons.

Inferring the coupling of different atmospheric layers requires observing spectral lines sensitive to the atmospheric parameters, particularly the magnetic field vector, at various heights. The best way to tackle this goal is to perform multi-line observations simultaneously. For instance, the new version of the Gregor Infrared Spectrograph instrument offers the possibility to observe the spectral lines at 8542 and 10830 A simultaneously for the first time. The first spectral window contains the Ca II 8542 A spectral line, while the Si I 10827 A transition and He I 10830 A triplet infrared lines can be found in the second spectral window. As the sensitivity to the atmospheric parameters and the height of formation of those transitions is different, combining them can help understand the properties of the solar photosphere and chromosphere and how they are magnetically coupled. Traditionally, the analysis of the Si I 10827 A transition assumes local thermodynamic equilibrium (LTE), which is not the best approximation to model this transition. Hence, in this work, we examine the potential of performing non-LTE (NLTE) inversions of the full Stokes vector of the Si I 10827 A spectral line. The results indicate that we properly infer the atmospheric parameters through an extended range of atmospheric layers in comparison with the LTE case (only valid for the spectral line wings, i.e., the low photosphere), with no impact on the robustness of the solution and just a minor increase in computational time. Thus, the NLTE assumption will help to accurately constrain the photospheric physical parameters when performing combined inversions with, e.g., the Ca II 8542 A spectral line.

Understanding the processes that drive the morphology and kinematics of molecular gas in galaxies is crucial for comprehending star formation and, ultimately, galaxy evolution. Using data obtained with the James Webb Space Telescope (JWST) and the Atacama Large Millimeter/submillimeter Array (ALMA), we study the behavior of the warm molecular gas at temperatures of hundreds of Kelvin and the cold molecular gas at tens of Kelvin in the galaxy MCG$-$05$-$23$-$16, which hosts an active galactic nucleus (AGN). Hubble Space Telescope (HST) images of this spheroidal galaxy, classified in the optical as S0, show a dust lane resembling a nuclear spiral and a surrounding ring. These features are also detected in CO(2$-$1) and H2, and their morphologies and kinematics are consistent with rotation plus local inward gas motions along the kinematic minor axis in the presence of a nuclear bar. The H2 transitions 0-0 S(3), 0-0 S(4), and 0-0 S(5), which trace warmer and more excited gas, show more disrupted kinematics than 0-0 S(1) and 0-0 S(2), including clumps of high-velocity dispersion (of up to $\sim$ 160 km/s), in regions devoid of CO(2$-$1). The kinematics of one of these clumps, located at $\sim$ 350 pc westward from the nucleus, are consistent with outflowing gas, possibly driven by localized star formation traced by Polycyclic Aromatic Hydrocarbon (PAH) emission at 11.3 ${\mu}$m. Overall, we observe a stratification of the molecular gas, with the colder gas located in the nuclear spiral, ring, and connecting arms, while most warmer gas with higher velocity-dispersion fills the inter-arm space. The compact jet, approximately 200 pc in size, detected with Very Large Array (VLA) observations, does not appear to significantly affect the distribution and kinematics of the molecular gas, possibly due to its limited intersection with the molecular gas disc.

The ongoing physical and chemical processes in planet-forming disks set the stage for planet formation. The asymmetric disk around the young star Oph-IRS 48 has one of the most well-characterised chemical inventories, showing molecular emission from a wide variety of species at the dust trap. One of the explanations for the asymmetric structure is dust trapping by a perturbation-induced vortex. We aim to constrain the excitation properties of the molecular species SO$_2$, CH$_3$OH, and H$_2$CO. We further characterise the extent of the molecular emission, through the determination of important physical and chemical timescales at the location of the dust trap. We also investigate whether the potential vortex can influence the observable temperature structure of the gas. Through a pixel-by-pixel rotational diagram analysis, we create rotational temperature and column density maps for SO$_2$ and CH$_3$OH, while temperature maps for H$_2$CO are created using line ratios. We find temperatures of $T\sim$55 K and $T\sim$125 K for SO$_2$ and CH$_3$OH, respectively, while the line ratios point towards temperatures of T$\sim$150-300 K for H$_2$CO. The rotational diagram of CH$_3$OH is dominated by scatter and subsequent non-LTE RADEX calculations suggest that both CH$_3$OH and H$_2$CO must be sub-thermally excited. The temperatures suggest that SO$_2$ comes from a layer deep in the disk, while CH$_3$OH and H$_2$CO originate from a higher layer. While a potential radial gradient is seen in the temperature map of SO$_2$, we do not find any hints of a vortex influencing the temperature structure. The determined turbulent mixing timescale is not able to explain the emitting heights of the molecules, but the photodissociation timescales are able to explain the wider azimuthal extents of SO$_2$ and H$_2$CO compared to CH$_3$OH, where a secondary, gas-phase formation reservoir is required for H$_2$CO.

In the astrophysics community it is common practice to model collisionless dust, entrained in a gas flow, as a pressureless fluid. However a pressureless fluid is fundamentally different from a collisionless fluid - the latter of which generically possess a non-zero anisotropic pressure or stress tensor. In this paper we derive a fluid model for collisionless dust, entrained in a turbulent gas, starting from the equations describing the motion of individual dust grains. We adopt a covariant formulation of our model to allow for the geometry and coordinate systems prevalent in astrophysics, and provide a closure valid for the accretion disc context. We show that the continuum mechanics properties of a dust fluid corresponds to a higher-dimensional anisotropic Maxwell fluid, after the extra dimensions are averaged out, with a dynamically important rheological stress tensor. This higher-dimensional treatment has the advantage of keeping the dust velocity and velocity of the fluid seen, and their respective moments, on the same footing. This results in a simplification of the constitutive relation describing the evolution of the dust Rheological stress.

The Itô and Stratonovich approaches are two ways to integrate stochastic differential equations. Detailed knowledge of the origin of the stochastic noise is needed to determine which approach suits a particular problem. I discuss this topic pedagogically in stochastic inflation, where the noise arises from a changing comoving coarse-graining scale or, equivalently, from `zooming in' into inflating space. I introduce a zoom-in scheme where deterministic evolution alternates with instantaneous zoom-in steps. I show that this alternating zoom-in scheme is equivalent to the Itô approach in the Markovian limit, while the Stratonovich approach doesn't have a similar interpretation. In the full non-Markovian setup, the difference vanishes. The framework of zoom-in schemes clarifies the relationship between computations in stochastic inflation, linear perturbation theory, and the classical $\Delta N$ formalism. It informs the numerical implementation of stochastic inflation and is a building block for a first-principles derivation of the stochastic equations.

As mixed with real pulsations, the reflection of super-Nyquist frequencies (SNFs) pose a threat to asteroseismic properties. Although SNFs have been studied in several pulsating stars, a systematic survey remains scarcely explored. Here we propose a method to identify SNFs from Kepler and TESS photometry by characterizing their periodic frequency modulations using a sliding Fourier transform. After analyzing long cadence photometry in the Kepler legacy, we have identified 304 SNFs in 56 stars from 45607 frequencies in $\sim600$ $\gamma$ Doradus stars, corresponding to a fraction of approximately $0.67\%$ and $9.2\%$, respectively. Most SNFs are detected in the frequency range of pressure mode over 120 $\mu$Hz and the fraction of SNF detection increases as frequency up to $\sim7\%$. We barely found two potential SNFs mixed with gravity modes in two $\gamma$ Doradus stars. These findings indicate that SNFs have a negligible impact on the global seismic properties, such as those derived from period spacing in $\gamma$ Doradus stars. However, we stress that SNFs must be carefully and systematically examined by this method in other pulsating stars, particularly $\delta$ Scuti and hot B subdwarf stars, to establish a solid foundation for precise asteroseismolgy of various types of pulsators.

Stellar evolution is driven by the changing composition of a star from nuclear reactions. At the late stages of evolution and during explosive events, the timescale can be short and drive strong hydrodynamic flows, making simulations of astrophysical reacting flows challenging. Over the past decades, the standard approach to modeling reactions in simulation codes has been operator splitting, using implicit integrators for reactions. Here we explore some of the assumptions in this standard approach and describe some techniques for improving the efficiency and accuracy of astrophysical reacting flows.

The gamma-ray binary HESS J0632+057 consists of a Be star and an undetected compact object in a $\sim$317 day orbit. The interpretation of the emission from this system is complicated by the lack of a clear orbital solution, as two different and incompatible orbital solutions were obtained by previous radial velocity studies of this source. In order to address this, we report on 24 new observations, covering $\sim$60 per cent of the orbit which we have undertaken with the Southern African Large Telescope (SALT). We obtained new radial velocity measurements from cross-correlation of the narrower spectral features, and by fitting Voigt profiles to the wings of the Balmer emission lines. Additionally, we find an indication of orbital variability in the equivalent widths and V/R of the Balmer lines. Using the combined data from this study, as well as archival data from the earlier radial velocity studies, we have derived updated orbital solutions. Using reported H $\alpha$ emission radial velocities - previously not considered for the orbital solution - along with the new SALT data, a solution is obtained where the brighter peak in the X-ray and gamma-ray light curves is closer to periastron. However, continuing sparse coverage in the data around the expected phases of periastron indicates that the orbital solution could be improved with further observation.

Around 50 years ago, the famous bet between Stephen Hawking and Kip Thorne on whether Cyg X-1 hosts a stellar-mass black hole became a well-known story in the history of black hole science. Today, Cyg X-1 is widely recognised as hosting a stellar-mass black hole with a mass of approximately 20 solar masses. With the advancement of X-ray telescopes, Cyg X-1 has become a prime laboratory for studies in stellar evolution, accretion physics, and high-energy plasma physics. In this review, we explore the latest results from X-ray observations of Cyg X-1, focusing on its implications for black hole spin, its role in stellar evolution, the geometry of the innermost accretion regions, and the plasma physics insights derived from its X-ray emissions. This review primarily focuses on Cyg X-1; however, the underlying physics applies to other black hole X-ray binaries and, to some extent, to AGNs.

Gravitational waves from core-collapse supernovae are a promising yet challenging target for detection due to the stochastic and complex nature of these signals. Conventional detection methods for core-collapse supernovae rely on excess energy searches because matched filtering has been hindered by the lack of well-defined waveform templates. However, numerical simulations of core-collapse supernovae have improved our understanding of the gravitational wave signals they emit, which enables us, for the first time, to construct a set of templates that closely resemble predictions from numerical simulations. In this study, we investigate the possibility of detecting gravitational waves from core-collapse supernovae using a matched-filtering methods. We construct a theoretically-informed template bank and use it to recover a core-collapse supernova signal injected into real LIGO-Virgo-KAGRA detector data. We evaluate the detection efficiency of the matched-filtering approach and how well the injected signal is reconstructed. We discuss the false alarm rate of our approach and investigate the main source of false triggers. We recover 88\% of the signals injected at a distance of 1 kpc and 50% of the signals injected at 2 kpc. For more than 50% of the recovered events, the underlying signal characteristics are reconstructed within an error of 15%. We discuss the strengths and limitations of this approach and identify areas for further improvements to advance the potential of matched filtering for supernova gravitational-wave detection. We also present the open-source Python package SynthGrav used to generate the template bank.

Consistent stellar evolution and nonlinear radial stellar pulsation calculations were carried out for models of asymptotic giant branch stars with initial masses $1.5M_\odot\le M_\mathrm{ZAMS}\le 3M_\odot$ and initial metal abundance $Z=0.006$. All the models are shown to be either the fundamental mode or the first overtone pulsators. The lower limit of the first overtone period increases with increasing mass of the Mira model from $\Pi_{1,\min}\approx 80$ days for $M=1.3M_\odot$ to $\Pi_{1,\min}\approx 120$ days for $M=2.6M_\odot$. The upper limit of the first overtone period and lower limit of the fundamental mode period depend on the stellar structure during mode switching and range from $\Pi_{1,\max}=130$, $\Pi_{0,\min}=190$ days for $M=0.96M_\odot$ to $\Pi_{1,\max}=210$, $\Pi_{0,\min}=430$ days for $M=2.2M_\odot$. The slope of the theoretical period--luminosity relation of Mira variables perceptibly increases with decreasing $Z$. Fourier spectra of the kinetic energy of twelve hydrodynamic models show a split of the fundamental mode maximum into several equidistant components. Frequency intervals between split components fall within the range $0.03 \le \Delta\nu/\nu_0 \le 0.1$. The superposition of radial oscillations with the fundamental mode splitting leads to the long-term amplitude variations with the cycle length from 10 to 30 times longer than the fundamental mode period. A more thorough analysis of hydrodynamic models is required for understanding the origin of the principal pulsation mode splitting.

Radio halos of edge-on galaxies are crucial for investigating cosmic ray propagation and magnetic field structures in galactic environments. We present VLA C-configuration S-band (2--4 GHz) observations of the spiral galaxy NGC 3556, a target from the Continuum Halos in Nearby Galaxies - an EVLA Survey (CHANG-ES). We estimate the thermal contribution to the radio emission from a combination of the H$\alpha$ and mid-IR data, and employ Rotation Measure Synthesis to reveal the magnetic field structures. In our data, NGC 3556 exhibits a box-like radio halo extending nearly 7 kpc from the galactic plane. The scale height of the total S-band intensity in the halo is $1.68\pm 0.29$ kpc, while that of the non-thermal intensity is $1.93\pm 0.28$ kpc. Fitting the data to a 1-D cosmic-ray transport model, we find advection to describe the cosmic-ray propagation within the halo better than diffusion, with advection speeds of $245 \pm 15$ km s$^{-1}$ and $205 \pm 25$ km s$^{-1}$ above and below the disk, respectively. The magnetic field is detected patchily across the galaxy, displaying a toroidal configuration in the rotation measure map. The mean equipartition magnetic field strength is approximately $8.3\ \mu$G in the disk and $4.5\ \mu$G in the halo. In addition, a bubble-like structure extends nearly 3~kpc into the southern halo, aligned with the polarized intensity and H$\alpha$ image, suggestive of superwinds generated by recent star formation feedback in the nuclear region.

The velocities of Ic-BL supernovae can be determined using two techniques (spline fitting and template fitting), sometimes resulting in different velocities for the same event. This work compares and contrasts both methods, identifying sources of error which are not accounted for by most authors and quantifying their impact on the final velocity measurement. Finally, it identifies the cause of velocity discrepancies for events measured using both methods. We quantified the impact of pre-smoothing the spectra prior to use of both methods using two well-sampled cases. To identify the source of velocity discrepancies, two cases were measured and directly compared. Additional sources of error for template fitting arise due to the choice of phase of the template spectrum ($\sim$1000 km/s) and smoothing of the input spectrum ($\sim$500 km/s). The impact of phase shifts is minimised at peak time. The spline fitting method tends to underestimate uncertainties by around 1000 km/s. This method can also be impacted by fine tuning of the smoothing parameters ($\sim$500-1000 km/s). Optimum smoothing parameters for different cases are presented along with suggestions for best practice. Direct comparison of both methods showed that velocity discrepancies are not always present, debunking the claim that the template fitting method always handles blending better than spline fitting. Spline fitting seems to struggle to handle blending only in cases where the Fe II features are superimposed on a red continuum, which biases the minimum of this feature towards the bluest line of the triplet, creating an artificially higher velocity. This situation may be relatively rare among Ic-BLs, based on typical temperature evolution. Both methods can be applied under certain circumstances with similar results. The morphology of the velocity evolution of an SN appears to be the same regardless of the method used.

The gravitational potential decay rate (DR) is caused by the cosmic acceleration of the universe, providing a direct probe into the existence of dark energy (DE). We present measurements of DR and explore its implications for DE models using the Data Release 9 galaxy catalog of DESI imaging surveys and the Planck cosmic microwave background maps. Our analysis includes six redshift bins within the range of $0.2\le z<1.4$ and achieves a total significance of 3.1$\sigma$, extending the DR measurements to a much higher redshift comparing to Dong et al. (2022), which focused on $0.2\le z<0.8$. Other improvements involve addressing potential systematics in the DR-related measurements of correlation functions, including imaging systematics and magnification bias. We explore the constraining power of DR both the $w$CDM model and the $w_0w_a$CDM model. We find that, the addition of DR can significantly improves DE constraints, over Sloan Digital Sky Survey baryon acoustic oscillation (BAO) data alone or PantheonPlus supernovae (SNe) compilation alone, although it shows only a modest improvement for DESI BAO. In the $w$CDM model, all three probes-DR, DESI BAO and SNe-favor $w=-1$. For the $w_0w_a$CDM, while DESI BAO prefers $w_0>-1$ and $w_a<0$, SNe Ia and DR data constrain $w_0=-0.94^{+0.11}_{-0.13}$ and $w_a=-0.22^{+0.57}_{-0.97}$. Namely SNe Ia and DR data has no preference on dynamical dark energy over $\Lambda$.

Gamma rays measured by the Fermi-LAT satellite tell us a lot about the processes taking place in high-energetic astrophysical objects. The fluxes coming from these objects are, however, extremely variable. Hence, gamma-ray light curves optimally use adaptive bin sizes in order to retrieve most information about the source dynamics and to combine gamma-ray observations in a multi-messenger perspective. However, standard adaptive binning approaches are slow, expensive and inaccurate in highly populated regions. Here, we present a novel, powerful, deep-learning-based approach to estimate the necessary time windows for adaptive binning light curves in Fermi-LAT data using raw photon data. The approach is shown to be fast and accurate. It can also be seen as a prototype to train machine-learning models for adaptive binning light curves for other astrophysical messengers.

Observations have revealed unique temperature profiles in hot Jupiter atmospheres. We propose that the energy transport by vertical mixing could lead to such thermal features. In our new scenario, strong absorbers, TiO and VO are not necessary. Vertical mixing could be naturally excited by atmospheric circulation or internal gravity wave breaking. We perform radiative transfer calculations by taking into account the vertical mixing driven energy transport. The radiative equilibrium (RE) is replaced by radiative-mixing equilibrium (RME). We investigate how the mixing strength, $K_{\rm zz}$, affects the atmospheric temperature-pressure profile. Strong mixing can heat the lower atmosphere and cool the upper atmosphere. This effect has important effects on the atmosphere thermal features that would form without mixing. In certain circumstances, it can induce temperature inversions in scenarios where the temperature monotonically increases with increasing pressure under conditions of lower thermal band opacity. Temperature inversions show up as $K_{\rm zz}$ increases with altitude due to shear interaction with the convection layer. The atmospheric thermal structure of HD~209458b can be well fitted with $K_{\rm zz} \propto (P/1\ {\rm bar})^{-1/2}\ {\rm cm}^{2} \ {\rm s}^{-1}$. Our findings suggest vertical mixing promotes temperature inversions and lowers $K_{\rm zz}$ estimates compared to prior studies. Incorporating chemical species into vertical mixing will significantly affect the thermal profile due to their temperature sensitivity.

Misaligned circumbinary disks will produce dust traffic jams during alignment or anti-alignment to the binary orbital plane. We conduct a hydrodynamical simulation of an initially misaligned circumbinary disk undergoing polar alignment with multiple dust species. Due to differential precession between the gas and dust components, multiple dust traffic jams are produced within the disk during polar alignment. The radial locations of the dust traffic jams depend on the Stokes number of the grains, which depends on grain size. We compute the dust temperature structure using post-processing radiative transfer to produce continuum images at cm-wavelengths. Multiple distinct rings emerge in the continuum images, corresponding to the dust traffic jams. The angular resolution of upcoming observations from SKA and ngVLA will be sufficient to detect centimeter-sized grains in protoplanetary disks and resolve the widths of dust traffic jams. Therefore, dust traffic jams resulting from the differential precession of gas and dust in misaligned circumbinary disks will be a prime target for more extended wavelength observations.

Galaxies grow and evolve in dark matter halos. Because dark matter is not visible, galaxies' halo masses ($\rm{M}_{\rm{halo}}$) must be inferred indirectly. We present a graph neural network (GNN) model for predicting $\rm{M}_{\rm{halo}}$ from stellar mass ($\rm{M}_{*}$) in simulated galaxy clusters using data from the IllustrisTNG simulation suite. Unlike traditional machine learning models like random forests, our GNN captures the information-rich substructure of galaxy clusters by using spatial and kinematic relationships between galaxy neighbour. A GNN model trained on the TNG-Cluster dataset and independently tested on the TNG300 simulation achieves superior predictive performance compared to other baseline models we tested. Future work will extend this approach to different simulations and real observational datasets to further validate the GNN model's ability to generalise.

A large sample of planet-planet scattering events for three planet systems with different orbital separations and masses is analyzed with a multiple regression model. The dependence of the time for the onset of instability on the masses of the planets and on their initial orbital separations is modeled with a quadratic function. The same analysis is applied to the timespan of the chaotic evolution dominated by mutual close encounters. The configurations with the less massive planet on an outside orbit are stable over longer timescales. The same configuration leads to shorter chaotic evolution times before the ejection of one planet. In about 70\% of the cases the lighter planet is the one escaping from the system. If a different separation is assumed between the inner and outer planet pairs, then the dominant effect on the instability time is due to the pair with the smaller separation, as a first approximation.

General Relativistic Magnetohydrodynamics (GRMHD) simulations are an indispensable tool in studying accretion onto compact objects. The Event Horizon Telescope (EHT) frequently uses libraries of ideal GRMHD simulations to interpret polarimetric, event-horizon-scale observations of supermassive black holes at the centers of galaxies. In this work, we present a library of ten non-radiative, ideal GRMHD simulations that were utilized by the EHT Collaboration in their analysis of Sagittarius A*. The parameter survey explores both low (SANE) and high (MAD) magnetization states across five black hole spins $a_{*}=-15/16,-1/2,0,+1/2,+15/16$ where each simulation was run out to $30,000\hspace{0.1cm}\mathrm{GM/c}^{3}$. We find the angular momentum and energy flux in SANE simulations closely matches the thin-disk value, with minor deviations in prograde models due to fluid forces. This leads to spin equilibrium around $a_{*}\sim0.94$, consistent with previous studies. We study the flow of conserved quantities in our simulations and find mass, angular momentum, and energy transport in SANE accretion flows to be primarily inward and fluid-dominated. MAD models produce powerful jets with outflow efficiency $>1$ for $a_{*}=+0.94$, leading to black hole spin-down in prograde cases. We observe outward directed energy and angular momentum fluxes on the horizon, as expected for the Blandford-Znajek mechanism. MAD accretion flows are sub-Keplerian and exhibit greater variability than their SANE counterpart. They are also hotter than SANE disks within $r\lesssim 10\hspace{0.1cm}\mathrm{GM/c}^{2}$. This study is accompanied by a public release of simulation data at \url{this http URL}.

It is widely assumed that dust opacities in molecular clouds follow a power-law profile with an index, $\beta$. Recent studies of the Orion Molecular Cloud (OMC) 2/3 complex, however, show a flattening in the spectral energy distribution (SED) at $ \lambda > 2$ mm implying non-constant indices on scales $\gtrsim$ 0.08 pc. The origin of this flattening is not yet known but it may be due to the intrinsic properties of the dust grains or contamination from other sources of emission. We investigate the SED slopes in OMC 2/3 further using observations of six protostellar cores with NOEMA from 2.9 mm to 3.6 mm and ALMA-ACA in Band 4 (1.9 -- 2.1 mm) and Band 5 (1.6 -- 1.8 mm) on core and envelope scales of $\sim 0.02 - 0.08$ pc. We confirm flattened opacity indices between 2.9 mm and 3.6 mm for the six cores with $\beta \approx -0.16 - 1.45$, which are notably lower than the $\beta$ values of $> 1.3$ measured for these sources on $0.08$ pc scales from single-dish data. Four sources have consistent SED slopes between the ALMA data and the NOEMA data. We propose that these sources may have a significant fraction of emission coming from large dust grains in embedded disks, which biases the emission more at longer wavelengths. Two sources, however, had inconsistent slopes between the ALMA and NOEMA data, indicating different origins of emission. These results highlight how care is needed when combining multi-scale observations or extrapolating single-band observations to other wavelengths.

Thermal changes in coronal loops are well-studied, both in quiescent active regions and in flaring scenarios. However, relatively little attention has been paid to loop emission in the hours before the onset of a solar flare; here, we present the findings of a study of over 50 off-limb flares of GOES class C5.0 and above. We investigated the integrated emission variability for Solar Dynamics Observatory Atmospheric Imaging Assembly channels 131, 171, 193, and 304 Ångstroms for 6 hours before each flare, and compared these quantities to the same time range and channels above active regions without proximal flaring. We find significantly increased emission variability in the 2-3 hours before flare onset, particularly for the 131 and 304 channels. This finding suggests a potential new flare prediction methodology. The emission trends between the channels are not consistently well-correlated, suggesting a somewhat chaotic thermal environment within the coronal portion of the loops that disturbs the commonly-observed heating and cooling cycles of quiescent active region loops. We present our approach, the resulting statistics, and discuss the implications for heating sources in these pre-flaring active regions.

We present the in-lab and on-sky performance for the upgraded 90 GHz focal plane of the Cosmology Large Angular Scale Surveyor (CLASS), which had four of its seven detector wafers updated during the austral winter of 2022. The update aimed to improve the transition-edge-sensor (TES) stability and bias range and to realize the high optical efficiency of the sensor design. Modifications included revised circuit terminations, electrical contact between the TES superconductor and the normal metal providing the bulk of the bolometer's heat capacity, and additional filtering on the TES bias lines. The upgrade was successful: 94% of detectors are stable down to 15% of the normal resistance, providing a wide overlapping range of bias voltages for all TESs on a wafer. The median telescope efficiency improved from $0.42^{+0.15}_{-0.22}$ to $0.60^{+0.10}_{-0.32}$ (68% quantiles). For the four upgraded wafers alone, median telescope efficiency increased to $0.65^{+0.06}_{-0.06}$. Given our efficiency estimate for the receiver optics, this telescope efficiency implies a detector efficiency exceeding $0.90$. The overall noise-equivalent temperature of the 90 GHz focal plane improved from 19 $\mu$K$\sqrt{s}$ to 11.3 $\mu$K$\sqrt{s}$.

Understanding and predicting the structure and evolution of coronal mass ejections (CMEs) in the heliosphere remains one of the most sought-after goals in heliophysics and space weather research. A powerful tool for improving current knowledge and capabilities consists of multi-spacecraft observations of the same event, which take place when two or more spacecraft fortuitously find themselves in the path of a single CME. Multi-probe events can not only supply useful data to evaluate the large-scale of CMEs from 1D in-situ trajectories, but also provide additional constraints and validation opportunities for CME propagation models. In this work, we analyse and simulate the coronal and heliospheric evolution of a slow, streamer-blowout CME that erupted on 23 September 2021 and was encountered in situ by four spacecraft approximately equally distributed in heliocentric distance between 0.4 and 1 au. We employ the Open Solar Physics Rapid Ensemble Information (OSPREI) modelling suite in ensemble mode to predict the CME arrival and structure in a hindcast fashion and to compute the "best-fit" solutions at the different spacecraft individually and together. We find that the spread in the predicted quantities increases with heliocentric distance, suggesting that there may be a maximum (angular and radial) separation between an inner and an outer probe beyond which estimates of the in-situ magnetic field orientation (parameterised by flux rope model geometry) increasingly diverge. We discuss the importance of these exceptional observations and the results of our investigation in the context of advancing our understanding of CME structure and evolution as well as improving space weather forecasts.

The atmospheres and surfaces of planets show tremendous amount of spatial variation, which has a direct effect on the spectrum of the object, even if this may not be spatially resolved. Here, we apply hyper realistic radiative simulations of Earth as an exoplanet comprising thousands of simulations and study the unresolved spectrum. The GlobES module on the Planetary Spectrum Generator was used, and we parameterized the atmosphere as described in the modern earth retrospective analysis for research and applications, MERRA2, database. The simulations were made into high spatial resolution images and compared to space based observations from the DSCOVR EPIC, at L1, and Himawari8, geostationary, satellites, confirming spatial variations and the spectral intensities of the simulations. The DSCOVR EPIC camera only functions in narrow wavelength bands, but strong agreement is demonstrated. It is shown that aerosols and small particles play an important role in defining Earths reflectance spectra, contributing significantly to its characteristic blue color. Subsequently, a comprehensive noise model is employed to constrain the exposure time required to detect O2, O3 and H2O as a function of varying ground and cloud cover for several concept observatories, including the habitable worlds observatory. Cloud coverage enhances the detectability of planets in reflected light, with important consequences for the design of the future HWO. The HWO concept would require between 3 to 10 times longer to observe the studied features than LUVOIR A but performs better than the HabEx without a starshade. The codes, routines, and the noise models are made publicly available.

An accurate estimation of the continuum excess emission from accretion spots and inner circumstellar disk regions is crucial for a proper derivation of fundamental stellar parameters in accreting systems. However, the presence of starspots can make disentangling the complicated multi-component emission in these systems challenging. Subtraction of a single-temperature spectral template is insufficient to account for the composite stellar emission, as we demonstrated in a recent campaign involving Weak-Lined T Tauri Stars. Here, we model the moderate resolution near-infrared spectra of Classical T Tauri Stars, presenting new spectral models that incorporate spotted stars plus emission from accretion hot-spots and a warm inner disk, allowing us to simultaneously reconstruct the entire 0.8-2.4 micrometer spectrum of our sixteen targets. Using these models, we re-derive the continuum excess emission. Our results indicate that accounting for starspots resolves the need to include a previously proposed intermediate temperature component in the IYJ excess, and highlights the importance of a proper treatment of starspots in studies of accreting low-mass stars.

The aim of this study is to identify quiescent galaxies in the 2-deg$^2$ COSMOS field at $z \sim 3.1$ and analyze their environment. Using data from the ODIN survey and COSMOS2020 catalog, we identify 24 massive quiescent galaxies (MQGs) with stellar masses $\geq 10^{10.6}$ and derive their star formation histories and quenching timescales using SED fitting with BAGPIPES. Voronoi-based density maps trace local and large-scale environments using Lyman-$\alpha$ Emitters and photometric galaxies. Results indicate uniformly short quenching timescales ($<$500 Myr) independent of environmental density, suggesting rapid internal mechanisms such as AGN feedback dominate over environmental factors. MQGs do not correlate with protoclusters or filaments, although some are near gas-rich filaments but show no rejuvenation. These findings suggest quenching at high redshift is driven primarily by internal processes rather than environmental interactions.

Empirical exoplanet mass-radius relations have been used to study the demographics and compositions of small exoplanets for many years. However, the heterogeneous nature of these measurements hinders robust statistical analysis of this population, particularly with regard to the masses of planets. For this reason, we perform a homogeneous and consistent re-analysis of the radial velocity (RV) observations of 85 small exoplanets using publicly available HARPS RV data and the fitting toolkit Pyaneti. For the entire sample, we run 12 different models to investigate the impact of modelling choices, including the use of multi-dimensional Gaussian Processes (GPs) to mitigate stellar activity. We find that the way orbital eccentricity is modelled can significantly impact the RV amplitude found in some cases. We also find that the addition of a GP to mitigate stellar activity does impact the RV amplitude found - though if the GP is modelled on activity indicators as well as the RVs the results are more robust. The RV amplitude found for every planet in our sample using all the models is made available for other groups to perform demographics studies. Finally, we provide a list of recommendations for the RV community moving forward.

Recent works have challenged our canonical view of RR Lyrae (RRL) stars as tracers of exclusively old populations ($\gtrsim10$~Gyr) by proposing a fraction of these stars to be of intermediate ages ($\sim$2-5~Gyr). Since it is currently not possible to infer stellar ages directly for individual RRL stars, our goal in this work is to search for these in association to intermediate-age clusters whose reliable ages can then be safely be attributed to the RRL. We used the Gaia DR3 Specific Object Study and OGLE IV public catalogues to search for RRL stars around stellar clusters older than 1~Gyr in the Large and Small Magellanic Clouds. Modelling membership probabilities based on proper motion and photometric distance we obtained a list of 302 RRL stars associated with Magellanic clusters. Of these, 23 RRL are likely members of 10 intermediate-age clusters: 3 and 7 in the Small and Large Magellanic Clouds, respectively. By modelling the inferred expectation values of the number of RRL stars per cluster, we inferred the delay time distribution of the RRL in three age ranges. For the old population ($>8$~Gyr) we find $2.6^{+0.4}_{-0.3}$ RRL$/10^5 M_\odot$. For the young (1-2 Gyr) and intermediate age (2-8 Gyr) populations we find rates of $0.9^{+0.3}_{-0.2}$ and $0.27^{+0.1}_{-0.09}$ RRL$/10^5 M_\odot$, respectively. While radial velocities are necessary for definitively confirming cluster memberships, the high-probability list of intermediate-age RRL stars presented here offers a promising opportunity for the first direct confirmation of these enigmatic stars.

In this thesis, we present a comprehensive and pedagogical overview of dark matter (DM). Chapter 1 discusses the main evidences for its existence, its properties, and potential candidates. We then explore major detection strategies, with Chapter 2 specifically dedicated to indirect detection. In the following chapters, we study the emission of secondary photons resulting from the interaction between DM products and the Galactic environment. Chapters 3 and 4 focus on DM as sub-GeV particles, analysing how the DM-produced electrons and positrons interact with ambient photons to generate X-rays through inverse Compton scattering. Comparing the predicted spectra with data from X-ray observatories yields strong constraints on sub-GeV DM. Chapter 5 extends these techniques to the case of primordial black hole (PBH) evaporation, imposing significant limits on PBHs as potential DM candidates.

The radiation chemistry and physics of solid N2O have been increasingly studied due to its potential presence on the surfaces of cold, outer Solar System bodies. However, to date, no study has investigated systematically the influence of temperature on this chemistry and physics. In this present study, crystalline N2O ices were irradiated using 2 keV electrons at five different temperatures in the 20-60 K range and the radiolytic dissociation of the molecular solid (as well as the radiolytic formation of seven product molecules) was quantified through the G-value. Our results indicate that temperature does indeed play a role in the radiolytic destruction of crystalline N2O, with higher temperatures being associated with higher destruction G-values. The formation G-values of NO, NO2, N2O2, N2O3, N2O4, N2O5, and O3 were also noted to vary with temperature, with each product molecule exhibiting a distinct trend. The applications of our experimental results to further understanding solid-phase radiation chemistry in the outer Solar System are discussed.

If the dark sector possesses long-range self-interactions, these interactions can source dramatic collective instabilities even in astrophysical settings where the collisional mean free path is long. Here, we focus on the specific case of dark matter halos composed of a dark $U(1)$ gauge sector undergoing a dissociative cluster merger. We study this by performing the first dedicated particle-in-cell plasma simulations of interacting dark matter streams, tracking the growth, formation, and saturation of instabilities through both the linear and nonlinear regimes. We find that these instabilities give rise to local (dark) electromagnetic inhomogeneities that serve as scattering sites, inducing an effective dynamic collisional cross-section. Mapping this effective cross-section onto existing results from large-scale simulations of the Bullet Cluster, we extend the limit on the dark charge-to-mass ratio by over ten orders of magnitude. Our results serve as a simple example of the rich phenomenology that may arise in a dark sector with long-range interactions and motivate future dedicated study of such ``dark plasmas.''

The Hayward regular BH solution attempted to resolve the curvature singularity issue by entering the domain of non-singular spacetimes. Recently, Dutta Roy and Kar (Phys. Rev. D 106, 044028) expanded this solution to encompass a broader range of spacetimes. These spacetimes are constructed based on the Damour-Solodukhin prescription, which involves introducing different metric parameters in the $g_{tt}$ and $g_{rr}$ components of the original Hayward line element, and are characterized by two parameters ($\sigma, \kappa$). This generalization gives rise to both known and novel regular/singular BHs as well as various types of wormhole spacetimes. In this work, we explore the spacetimes that emerge for different values of ($\sigma, \kappa$) from the generalized Hayward metric, particularly focusing on their shadows in vacuum and when surrounded by plasma. Intriguingly, we observe the presence of both photon and anti-photon spheres for certain regular spacetimes. Our study highlights the differences in the shadows of different types of regular spacetime compared to those of the singular BH derived from the generalized Hayward metric and also sheds light on the impact of plasma on the shadow radius.

We show that nickel oxide, which is already a very promising target to look for sub-MeV dark matter scattering, can be employed to hunt axion dark matter, with masses in the meV range and couplings to electrons allowing them to potentially be QCD axions. We describe the interactions between axions and the collective excitations of nickel oxide in terms of a universal effective field theory, built solely out of symmetry arguments. The processes of conversion into one or two excitations provide, respectively, a narrowband and a broadband channel for the axion search, and the possibility of varying an external magnetic field up to a phase transition point allows to cover a large portion of a yet unexplored parameter space, reaching axion masses down to few fractions of an meV. Our results underline nickel oxide as an ideal candidate for a multi-purpose target for light dark matter searches.

Motivated by known facts about effective field theory and non-Abelian gauge theory, we argue that the post-Newtonian approximation might fail even in the limit of weak fields and small velocities under certain conditions. Namely, the post-Newtonian approximation might break down for wide extended bodies with angular momentum, where angular momentum spans significant spacetime curvature. We construct a novel dimensionless quantity that samples this breakdown, and we evaluate it by means of existing analytical solutions of rotating extended bodies and observational data. We give estimates for galaxies and binary systems, as well as our home in the Cosmos, Laniakea. We thus propose that a novel effective field theory of general relativity is needed to account for the onset of nonlocal angular momentum effetcs, with significant consequences for gravitational physics and cosmology at large.

We extend the recently developed quantum van der Waals quarkyonic matter to non-zero isospin asymmetries by utilizing the two-component van der Waals equation with a generalized excluded volume prescription. The isospin dependence of van der Waals interaction parameters is determined by constraints on the symmetry energy, slope of the symmetry energy, and nuclear ground state properties. We find that the speed of sound has a peak for all values of the asymmetry parameter, signifying a transition to quarkyonic matter. The quarkyonic matter onset density is found to have a mild dependence on isospin asymmetry, with specific details influenced by the isospin dependence of the repulsive interactions. We also incorporate leptonic degrees of freedom and explore the neutron star matter equation of state, calculating mass-radius relations and tidal properties of neutron stars. We find that quarkyonic matter supports heavy neutron stars with a maximum mass of at least 2.6 solar masses. We observe quantitatively different behavior for the excluded volume cases of isospin-blind ($b_{n}=b_{pn}$) and isospin-dependent ($b_{n} \neq b_{pn}$) repulsion, the latter being preferred by observational constraints.

The presence of a massive body between the Earth and a gravitational-wave source will produce the so-called gravitational lensing effect. In the case of strong lensing, it leads to the observation of multiple deformed copies of the initial wave. Machine-learning (ML) models have been proposed for identifying these copies much faster than optimal Bayesian methods, as will be needed with the detection rate of next-generation detector. Most of these ML models are based on a time-frequency representation of the data that discards the phase information. We introduce a neural network that directly uses the time series data to retain the phase, limit the pre-processing time and keep a one-dimensional input. We show that our model is more efficient than the base model used on time-frequency maps at any False Alarm Rate (FPR), up to $\sim 5$ times more for an FPR of $10^{-4}$. We also show that it is not significantly impacted by the choice of waveform model, by lensing-induced phase shifts and by reasonable errors on the merger time that induce a misalignment of the waves in the input.

We present a single-fluid approach for the simulation of partially-ionized plasmas (PIPs) which is designed to capture the non-ideal effects introduced by neutrals while remaining close in computational efficiency to single-fluid MHD. This is achieved using a model which treats the entire partially-ionized plasma as a single mixture, which renders internal ionization/recombination source terms unnecessary as both the charged and neutral species are part of the mixture's conservative system. Instead, the effects of ionization and the differing physics of the species are encapsulated as material properties of the mixture. Furthermore, the differing dynamics between the charged and neutral species is captured using a relative-velocity quantity, which impacts the bulk behavior of the mixture in a manner similar to the treatment of the ion-electron relative-velocity as current in MHD. Unlike fully-ionized plasmas, the species composition of a PIP changes rapidly with its thermodynamic state. This is captured through a look-up table referred to as the tabulated equation of state (TabEoS), which is constructed prior to runtime using empirical physicochemical databases and efficiently provides the ionization fraction and other material properties of the PIP specific to the thermodynamic state of each computational cell. Crucially, the use of TabEoS also allows our approach to self-consistently capture the non-linear feedback cycle between the PIP's macroscopic behavior and the microscopic physics of its internal particles, which is neglected in many fluid simulations of plasmas today.

In this study, we explore the back reaction of phase transitions in the spectator sector on the inflaton field during slow-roll inflation. Due to the significant excursion of the inflaton field, these phase transitions are likely to occur and can induce substantial non-Gaussian correlations in the curvature perturbation. Our results suggest that these correlations could be detectable by future observations of the cosmic microwave background radiation and large-scale structure surveys. Furthermore, we demonstrate that in certain parameter spaces, a scaling non-Gaussian signal can be produced, offering deeper insights into both the inflaton and spectator sectors. Additionally, phase transitions during inflation can generate gravitational wave signals with distinctive signatures, potentially explaining observations made by pulsar timing array experiments. The associated non-Gaussian correlations provide collateral evidence for these phase transitions.