Papers for Tuesday, Dec 17 2024

Separating a stochastic gravitational wave background (SGWB) from noise is a challenging statistical task. One approach to establishing a detection criterion for the SGWB is using Bayesian evidence. If the evidence ratio (Bayes factor) between models with and without the signal exceeds a certain threshold, the signal is considered detected. We present a formalism to compute the averaged Bayes factor, incorporating instrumental-noise and SGWB uncertainties. As an example, we consider the case of power-law-shaped SGWB in LISA and generate the corresponding \textit{bayesian sensitivity curve}. Unlike existing methods in the literature, which typically neglect uncertainties in both the signal and noise, our approach provides a reliable and realistic alternative. This flexible framework opens avenues for more robust stochastic gravitational wave background detection across gravitational-wave experiments.

With ESA's PLATO space mission set for launch in December 2026, a new photometric legacy and a future of new scientific discoveries await. In this work we investigate PLATO's potential for observing pulsating stars across the Hertzsprung-Russell diagram as part of the PLATO Complimentary Science program (PLATO-CS). Specifically, a PLATO mock asteroseismic catalogue (MOCKA) of intermediate to massive stars is presented as a benchmark to highlight the asteroseismic yield of PLATO-CS in a quantitative way. MOCKA includes simulations of $\beta$~Cephei, slowly pulsating B (SPB), $\delta$~Scuti, $\gamma$~Doradus, RR Lyrae, Cepheid, hot subdwarf, and white dwarf stars. In particular, main-sequence gravity (g) mode pulsators are of interest as some of these stars form an important foundation for the scientific calibration of PLATO. MOCKA is based on a magnitude limited ($G\lesssim17$) \textit{Gaia} catalogue and is a product of realistic end-to-end \texttt{PlatoSim} simulations of stars for the first PLATO pointing field in the Southern hemisphere, which will be observed for a minimally 2-yr duration. We show that an abundant spectrum of frequencies is achievable across a wide range of magnitudes and co-pointing PLATO cameras. Within the magnitude limited regimes simulated ($G \lesssim 14$ for $\gamma$~Doradus stars and $G \lesssim 16$ for SPB stars) the dominant g-mode frequency is recovered in more than $95\%$ of the cases. MOCKA help us to understand the limits of the PLATO mission as well as highlight the opportunities to push astrophysics beyond current stellar models. All data products of this paper are made available to the community for further exploration. The key data products of MOCKA are the magnitude limited \textit{Gaia} catalogue of the first PLATO pointing field, together with fully reduced light curves from multi-camera observations for each pulsation class.

In this paper, we investigate the freezing quintessence scenario in late-time cosmic expansion using a non-linear $f(R, L_m)$ gravity model, $f(R,L_m)=\frac{R}{2}+L_m^\alpha$, where $\alpha$ is a free parameter. We consider a solution for this model using an appropriate parametrization of the scale factor, and then the model is constrained by observational datasets, including CC, Pantheon+ (SN), and CC+SN+BAO. Our analysis yields results aligning closely with observational data. The Hubble parameter, deceleration parameter, matter-energy density, and EoS parameter of our model exhibit expected trends over cosmic time, supporting its physical validity. Furthermore, the model demonstrates consistency with the $\Lambda$CDM model in late times, displaying freezing behavior in the $\omega - \omega'$ plane and stability against density perturbations. Our findings suggest that the modified $f(R, L_m)$ gravity model is a credible approach to describing the universe's accelerating phase.

Classifying variable stars is key for understanding stellar evolution and galactic dynamics. With the demands of large astronomical surveys, machine learning models, especially attention-based neural networks, have become the state-of-the-art. While achieving high accuracy is crucial, enhancing model interpretability and uncertainty estimation is equally important to ensure that insights are both reliable and comprehensible. We aim to enhance transformer-based models for classifying astronomical light curves by incorporating uncertainty estimation techniques to detect misclassified instances. We tested our methods on labeled datasets from MACHO, OGLE-III, and ATLAS, introducing a framework that significantly improves the reliability of automated classification for the next-generation surveys. We used Astromer, a transformer-based encoder designed for capturing representations of single-band light curves. We enhanced its capabilities by applying three methods for quantifying uncertainty: Monte Carlo Dropout (MC Dropout), Hierarchical Stochastic Attention (HSA), and a novel hybrid method combining both approaches, which we have named Hierarchical Attention with Monte Carlo Dropout (HA-MC Dropout). We compared these methods against a baseline of deep ensembles (DEs). To estimate uncertainty estimation scores for the misclassification task, we selected Sampled Maximum Probability (SMP), Probability Variance (PV), and Bayesian Active Learning by Disagreement (BALD) as uncertainty estimates. In predictive performance tests, HA-MC Dropout outperforms the baseline, achieving macro F1-scores of 79.8+-0.5 on OGLE, 84+-1.3 on ATLAS, and 76.6+-1.8 on MACHO. When comparing the PV score values, the quality of uncertainty estimation by HA-MC Dropout surpasses that of all other methods, with improvements of 2.5+-2.3 for MACHO, 3.3+-2.1 for ATLAS and 8.5+-1.6 for OGLE-III.

Astronomers often deal with data where the covariates and the dependent variable are measured with heteroscedastic non-Gaussian error. For instance, while TESS and Kepler datasets provide a wealth of information, addressing the challenges of measurement errors and systematic biases is critical for extracting reliable scientific insights and improving machine learning models' performance. Although techniques have been developed for estimating regression parameters for these data, few techniques exist to construct prediction intervals with finite sample coverage guarantees. To address this issue, we tailor the conformal prediction approach to our application. We empirically demonstrate that this method gives finite sample control over Type I error probabilities under a variety of assumptions on the measurement errors in the observed data. Further, we demonstrate how the conformal prediction method could be used for constructing prediction intervals for unobserved exoplanet masses using established broken power-law relationships between masses and radii found in the literature.

In our previous study, we identified a shift in the synchrotron peak frequency of the blazar B2 1308$+$326 from 10$^{12.9}$ Hz to 10$^{14.8}$ Hz during a flare, suggesting it could be a changing-look blazar (CLB). In this work, we investigate the CL behaviour of B2 1308+326 by analysing a newly acquired optical spectrum and comparing it with an archival spectrum. We find that between the two epochs, the continuum flux increased by a factor of $\sim$4.4, while the Mg II emission line flux decreased by a factor of 1.4$\pm$0.2. Additionally, the equivalent width of the Mg II line reduced from $\sim 20$ Å\ to $\sim 3$ Å, indicating an apparent shift from a flat-spectrum radio quasar (FSRQ) class to a BL Lacertae (BL Lac) class. Despite this apparent change, the ratio of accretion disk luminosity to Eddington luminosity remains $>$ 10$^{-2}$ during both epochs, indicating efficient accretion persists in B2 1308$+$326. The measured black hole mass remains consistent with an average $\log M_{\rm BH} = 8.44$ M$_{\odot}$. Our findings suggest that B2 1308$+$326 is not a genuine CLB, but rather an intrinsic FSRQ that emerges as a BL Lac during high-flux states due to enhanced non-thermal emission.

We present a spectroscopic analysis of two star-forming galaxies at z~5 observed with JWST/NIRSpec as part of the Early eXtragalactic Continuum and Emission Line Science (EXCELS) survey. The detection of the C III]$\lambda\lambda$1906,09, [O II]$\lambda\lambda$3726,29, [O III]$\lambda\lambda$4363,5007, and [N II]$\lambda$6584 nebular emission lines enables investigation of the C/O, N/O, and C/N abundance ratios using the temperature-sensitive method. The two galaxies have stellar masses of log($M_{\star}$/M$_{\odot}$ ) = 8.13$\pm$0.09 and log($M_{\star}$/M$_{\odot}$ )=8.52$\pm$0.13 and corresponding metallicities of Z~0.2Z$_{\odot}$ and Z~0.3Z$_{\odot}$. These metallicities are somewhat higher than is typical for other z>5 galaxies with similar stellar mass and are in fact comparable to high-redshift analogue galaxies at z~0. Both galaxies display evidence for N/O enhancement with respect to the z~0 sample, with log(N/O)=-1.07$\pm$0.17 and log(N/O)=-0.86$\pm$0.15 respectively. In contrast, we find low C abundances, with log(C/O)=-0.82$\pm$0.22 and log(C/O)=-1.02$\pm$0.22, consistent with the predicted yields of core-collapse supernovae. Following the trend observed in other high-redshift sources, we find that the C/N ratios are lower at fixed O/H compared to the majority of local galaxies. In contrast to the top-heavy IMF invoked in some studies to explain low C/N ratios in metal-poor galaxies, we find, via comparison to chemical evolution models, that a standard or bottom-heavy IMF better explains the observed abundance ratios in more enriched systems due to an increase in N-enrichment from intermediate mass (4-7M$_{\odot}$) stars. Our results demonstrate that robust measurements of CNO abundances with JWST can reveal unique enrichment pathways in galaxies as a function of both metallicity and redshift.

Galaxies are part of a vast cosmic ecosystem, embedded in an extensive gaseous reservoir that regulates their growth by providing the necessary fuel for star formation while preserving a fossil record of past interactions, outflows, and feedback processes. The circumgalactic medium (CGM) contains multiphase gas spanning a broad dynamic range in spatial scale, density, and temperature, with its thermodynamic and chemical properties deeply linked to the star formation histories of galaxies. As a rich laboratory for studying gas physics, the CGM offers unique insights into the processes governing gas cooling, heating, and material transfer between galaxies and their surroundings. Chemical tagging, based on the relative abundances of multiple elements, serves as a powerful timing tool to trace the origin of the gas and connect the stars in the interstellar medium (ISM) to the diffuse CGM. Developing a complete understanding of the CGM and its cosmic evolution requires multi-wavelength observational tools, ranging from X-ray, UV/optical, and sub-mm to radio, to probe this diffuse gas in both absorption and emission.

We present the polarization spectra of the nucleus of 3C 270. We confirm that the polarization angle of both the continuum and the emission lines are close to perpendicular to the jet direction after careful correction of interstellar polarization, which indicates polar scattering. The Stokes flux spectrum resembles the total flux spectrum, with no need for a broad component from the Broad Line Region. Over 94% of the Seyfert I type broad line profile would be significantly detected if present in our polarized flux spectrum. We favor that we are mostly likely observing beamed synchrotron instead of Big Blue Bump, and the innermost Narrow Line Region through reflection. This makes 3C 270 the third known case, after NGC 4258 and Centaurus A, where only narrow lines (and beamed continuum, if present) are scattered, with no evidence of an underlying Big Blue Bump.

We have performed 3D MHD simulations to study the effect of partial ionization in the process of magnetic flux emergence in the Sun. In fact, we continue previous work and we now focus: 1) on the emergence of the magnetic fields above the solar photosphere and 2) on the eruptive activity, which follows the emergence into the corona. We find that in the simulations with partial ionization (PI), the structure of the emerging field consists of arch-like fieldlines with very little twist since the axis of the initial rising field remains below the photosphere. The plasma inside the emerging volume is less dense and it is moving faster compared to the fully ionized (FI) simulation. In both cases, new flux ropes (FR) are formed due to reconnection between emerging fieldlines, and they eventually erupt in an ejective manner towards the outer solar atmosphere. We are witnessing three major eruptions in both simulations. At least for the first eruption, the formation of the eruptive FR occurs in the low atmosphere in the FI case and at coronal heights in the PI case. Also, in the first PI eruption, part of the eruptive FR carries neutrals in the high atmosphere, for a short period of time. Overall, the eruptions are relatively faster in the PI case, while a considerable amount of axial flux is found above the photosphere during the eruptions in both simulations.

Binary star systems can accrete material originating from a circumbinary disc. Since it is common for the circumbinary disc to be tilted with respect to the binary orbital plane, we test whether the accretion dynamics can be a diagnostic for binary-disc misalignment. We present hydrodynamical simulations to model the accretion flow from a circumbinary disc around an eccentric binary with initial tilts ranging from $0^\circ$ to $180^\circ$ in increments of $15^\circ$. Based on the initial tilt, the circumbinary disc will align towards three different configurations: prograde coplanar, polar, or retrograde coplanar. For discs with initial tilts evolving towards prograde coplanar alignment, the accretion rates onto the primary and secondary stars exhibit alternating preferential accretion. Circumbinary discs evolving towards polar alignment exhibit no alternating preferential accretion onto the binary unless the initial tilt is close to the critical tilt that sets the boundary between coplanar or polar alignment. Such cases cause strong disc warping, leading to disc breaking. The inner disc becomes eccentric, leading to alternating preferential accretion onto the binary. As the break propagates outward, the disc tilt damps towards a polar state and the disc eccentricity decreases. As the disc re-circularizes, the accretion rate transitions back from alternating preferential accretion to non-alternating accretion. Lastly, no alternating preferential accretion exists for discs undergoing retrograde coplanar alignment. From the summary of the accretion rates from our suite of SPH simulations, it is evident that the accretion rate evolution can be affected by the initial tilt and subsequent evolution of the circumbinary disc.

Searches for impulsive, astrophysical transients are often highly computationally demanding. A notable example is the dedispersion process required for performing blind searches for Fast Radio Bursts (FRBs) in radio telescope data. We introduce a novel approach - Efficient Summation of Arbitrary Masks (ESAM) - which efficiently computes 1-D convolution of many arbitrary 2-D masks, and can be used to carry out dedispersion over thousands of dispersion trials efficiently. Our method matches the accuracy of the traditional brute force technique in recovering the desired Signal-to-Noise ratio (S/N) while reducing computational cost by around a factor of 10. We compare its performance with existing dedispersion algorithms, such as the Fast Dispersion Measure Transform (FDMT) algorithm, and demonstrate how ESAM provides freedom to choose arbitrary masks and further optimise computational cost versus accuracy. We explore the potential applications of ESAM beyond FRB searches.

We describe some of our group's recent work on common envelopes. Our goal is to understand the onset and outcomes of dynamically unstable mass transfer, including the properties of the binaries left behind and the outflows during the common envelope stage. We have also started thinking about light curves of common envelope events. During a talk at StanFest, a meeting in honour of Stan Owocki's retirement held in Leuven in July, 2024, the first author reported on some of the results we recently obtained and briefly outlined future prospects.

Planetary bodies are formed by coagulation of solid dust grains in protoplanetary disks. Therefore, it is crucial to constrain the physical and chemical properties of the dust grains. In this study, we measure the dust albedo at mm-wavelength, which depends on dust properties at the disk midplane. Since the albedo and dust temperature are generally degenerate in observed thermal dust emission, it is challenging to determine them simultaneously. We propose to break this degeneracy by using multiple optically-thin molecular lines as a dust-albedo independent thermometer. In practice, we employ pressure-broadened CO line wings that provide an exceptionally high signal-to-noise ratio as an optically thin line. We model the CO $J=2-1$ and $3-2$ spectra observed by the Atacama Large Millimeter/sub-millimeter Array (ALMA) at the inner region ($r<6\ {\rm au}$) of the TW Hya disk and successfully derived the midplane temperature. Combining multi-band continuum observations, we constrain the albedo spectrum at $0.9-3$ mm for the first time without assuming a dust opacity model. The albedo at these wavelengths is high, $\sim0.5-0.8$, and broadly consistent with the Ricci et al. (2010), DIANA, and DSHARP dust models. Even without assuming dust composition, we estimate the maximum grain size to be $\sim 340\ \mu m$, the power law index of the grain size distribution to be $>-4.1$, and porosity to be $<0.96$. The derived dust size may suggest efficient fragmentation with the threshold velocity of $\sim 0.08\ {\rm m\ s^{-1}}$. We also note that the absolute flux uncertainty of $\sim10\%$ ($1\sigma$) is measured and used in the analysis, which is approximately twice the usually assumed value.

Using nearly simultaneous radio, near-infrared, optical, and ultraviolet data collected since 2009, we constructed 106 spectral energy distributions (SEDs) of the blazar OJ 287. These SEDs were well-fitted by a log-parabolic model. By classifying the data into `flare' and `quiescent' segments, we found that the median flux at peak frequency of the SEDs during flare segments was 0.37$\pm$0.22 dex higher compared to quiescent segments, while no significant differences were observed in the median values of the curvature parameter $b$ or the peak frequency $\log \nu_{\mathrm{p}}$. A significant bluer-when-brighter trend was confirmed through a relation between $V$ magnitude and $B-V$ color index, with this trend being stronger in the flare segments. Additionally, a significant anti-correlation was detected between $\log \nu_{\mathrm{p}}$ and $b$, with a slope of 5.79 in the relation between $1/b$ and $\log \nu_{\mathrm{p}}$, closer to the prediction from a statistical acceleration model other than a stochastic acceleration interpretation, though a notable discrepancy persists. This discrepancy indicates that additional factors, such as deviations from idealized conditions or radiative contributions-such as thermal emission from the accretion disk in the optical-UV range during quiescent states-may play a role in producing the observed steeper slope. Within the framework of statistical acceleration mechanism, lack of correlation between change in peak intensity and change in peak frequency suggests that change in electron energy distribution is unlikely to be responsible for the time-dependent SED changes. Instead, changes in Doppler boosting or magnetic fields may have a greater influence.

We propose the HI 21-cm power spectrum from the post-reionization epoch as a probe of a cosmological model with decaying dark matter particles. The unstable particles are assumed to undergo a 2-body decay into a massless and massive daughter. We assume, that a fraction $f$ of the total dark matter budget to be, unstable and quantify the decay using the life-time $\Gamma^{-1}$ and the relative mass splitting $\epsilon$ between the parent and the massive daughter. The redshift space anisotropic power spectrum of the post-reionization 21-cm signal brightness temperature, as a tracer of the dark matter clustering, imprints the decaying dark matter model through its effect on background evolution and the suppression of power on small this http URL find that with an idealized futuristic intensity mapping experiment with a SKA-I Mid like radio-array, $\epsilon$ and $\Gamma$ can be measured at $3.1\%$ and $4.64\%$ around their fiducial values of $\epsilon = 0.01 $ and $\Gamma = 0.074 {\rm Gyr}^{-1}$ respectively.

By analyzing data from DESI Legacy Imaging Survey of the dwarf galaxies in the Arecibo Legacy Fast Alfa Survey, we have identified five ultra-diffuse galaxies (UDGs) featuring central pseudo-bulges. These UDGs display blue pseudo-bulges with Sérsic indices $n<2.5$ and effective radii spanning 300-700 pc, along with bluer thin stellar disks exhibiting low surface brightness and expansive effective radii that align with the UDG definition. The rotation velocities of these UDGs, determined using HI line widths and optical inclinations, exceed those of most dwarf galaxies of similar mass, suggesting the high halo spins or substantial dark matter halos. We propose that these UDGs likely formed through mergers of dwarf galaxies lacking old stars in their progenitors, resulting in the development of central bulge-like structures during starbursts triggered by the mergers, while also enhancing their halo spin. Subsequent gas accretion facilitated the formation of extended stellar disks. It is also worth noting the possibility that these UDGs could alternatively represent ``failed $L^{\star}$ galaxies'' with massive dark matter halos but reduced star formation efficiencies. If future high-resolution HI observations confirm the presence of massive halos around these UDGs, they may have formed due to intense AGN feedback in the early universe, and may be the descendants of ``little red dots'' observed by the James Webb Space Telescope, which are characterized by heightened central black hole masses and intensified accretion and feedback processes in the early universe.

Due to its peculiar and highly variable nature, the blazar 3C 454.3 has been extensively monitored by the WEBT team. Here, we present for the first time these long-term optical flux and color variability results using data acquired in B, V, R, and I bands over a time span of $\sim$ 2 decades. We include data from WEBT collaborators and public archives such as SMARTS, Steward Observatory, and ZTF. The data are binned and segmented to study the source over this long term when more regular sampling was available. During our study, the long-term spectral variability reveals a redder when brighter (RWB) trend, which, however, stabilizes at a particular brightness cutoff $\sim$ 14.5 mag in the I-band, after which it saturates and evolves into a complex state. This trend indicates increasing jet emission dominance over accretion disk emission until jet emission completely dominates. Plots of the spectral index variation (following $F_{\nu} \propto \nu^{-\alpha}$) reveal a bimodal distribution using a one-day binning. These correlate with two extreme phases of 3C 454.3, an outburst or high flux state and quiescent or low flux state, which are respectively jet and accretion disk dominated. We have also conducted intra-day variability studies of nine light curves and found that six of them are variable. Discrete Correlation Function (DCF) analysis between different optical waveband pairs peak at zero lags, indicating co-spatial emission in different optical bands.

Low-density and unbound stellar groups, OB associations have been historically delineated through their bright and massive members. They have been analysed for decades, but the arrival of Hipparcos, and more recently of Gaia led to a change of paradigm by allowing the identification of more reliable members using parallaxes and proper motions. This renewed interest offers an opportunity to emphasize the role of OB associations across many areas of astronomy. In this review, I highlight their importance across multiple scales: how OB associations constitute suitable sites to study massive stars and stellar multiplicity, their relation with star clusters, their interactions with the interstellar medium through the feedback of their massive members, and how they shape the structure and evolution of the Milky Way and beyond.

We examine lengthy radio light curves of the flat spectrum radio galaxy 3C 454.3 for possible quasi-periodic oscillations (QPOs). The data used in this work were collected at five radio frequencies, 4.8, 8.0, 14.5, 22.0, and 37.0 GHz between 1979--2013 as observed at the University of Michigan Radio Astronomical Observatory, Crimean Astrophysical Observatory, and Aalto University Mets{ä}hovi Radio Observatory. We employ generalized Lomb-Scargle periodogram and weighted wavelet transform analyses to search for periodicities in these light curves. We confirm a QPO period of $\sim$ 2000 day to be at least 4$\sigma$ significant using both methods at all five radio frequencies between 1979 and 2007, after which a strong flare changed the character of the light curve. We also find a $\sim$~600 day period which is at least 4$\sigma$ significant, but only in the 22.0 and 37.0 GHz light curves. We briefly discuss physical mechanisms capable of producing such variations.

The reverberation mapping (RM) technique has seen wide applications in probing geometry and kinematics of broad-line regions (BLRs) and measuring masses of supermassive black holes (SMBHs) in active galactic nuclei. However, the key quantities in RM analysis like emissivity, responsivity, transfer functions, and mean and root-mean-square (RMS) spectra are fragmentally defined in the literature and largely lack a unified formulation. Here, we establish a rigorous framework for BLR RM and include a locally dependent responsivity according to photoionization calculations. The mean and RMS spectra are analytically expressed with emissivity- and responsivity-weighted transfer functions, respectively. We demonstrate that the RMS spectrum is proportional to the responsivity-weighted transfer function only when the continuum variation timescale is much longer than the typical extension in time delay of the BLR, otherwise, biases arise in the obtained RMS line widths. The long-standing phenomenon as to the different shapes between mean and RMS spectra can be explained by a radial-increasing responsivity of BLRs. The debate on the choice of emission line widths for SMBH mass measurements is explored and the virial factors are suggested to also depend on the luminosity states, in addition to the geometry and kinematics of BLRs.

We propose that the drastic photometric and spectroscopic changes affecting the symbiotic star V694 Mon since 2018, are due to its transition from the accreting-only state to steady hydrogen-burning on the surface of the white dwarf, closely mimicking the pattern followed by V4368 Sgr. The phase of peak optical brightness and weakest emission lines has probably been reached in early 2024. The high-velocity absorptions powered by jet-ejection and the wild flickering, which dominated the century-long quiescence, should not reappear as long has nuclear burning will hold (time-scale of decades). The 3500 Lsun burning luminosity suggests a mass of 0.60 Msun for the WD.

O-stars generally show clear signs of strong line-broadening (in addition to rotational broadening) in their photospheric absorption lines (typically referred to as 'macroturbulence'), believed to originate in a turbulent sub-surface zone associated with enhanced opacities due to recombination of iron-group elements (at T~ 150-200 kK). O-stars with detected global magnetic fields also display such macroturbulence; the sole exception to this is NGC 1624-2, which also has the strongest (by far) detected field of the known magnetic O-stars. It has been suggested that this lack of additional line-broadening is because NGC 1624-2's exceptionally strong magnetic field might be able able to suppress the turbulent velocity field generated in the iron opacity peak zone. For moderately strong magnetic cases (~1 kG) the simulated atmospheres are highly structured characterised by large root-mean-square velocities, and our results are qualitatively similar to those found in previous non-magnetic studies. By contrast, we find that a strong horizontal magnetic field in excess of 10 kG can indeed suppress the large velocity fluctuations and thus stabilise (and thereby also inflate) the atmosphere of a typical early O-star in the Galaxy. On the other hand, an equally strong radial field is only able to suppress horizontal motions, and as a consequence these models exhibit significant radial fluctuations. Our simulations provide an overall physical rationale as to why NGC 1624-2 with its strong ~20 kG dipolar field lacks the large macroturbulent line broadening that all other known slowly rotating early O-stars exhibit. However, they also highlight the importance of field geometry for controlling the atmospheric dynamics in massive and luminous stars that are strongly magnetic, tentatively suggesting latitudinal dependence of macroturbulence and basic photospheric parameters.

Fuzzy (wave) dark matter (FDM), the dynamical model underlying an ultralight bosonic dark matter species, produces a rich set of non-gravitational signatures that distinguishes it markedly from the phenomenologically related warm (particle) dark matter (WDM) scenario. The emergence of extended interference fringes hosted by cosmic filaments is one such phenomenon reported by cosmological simulations, and a detailed understanding of such may strengthen existing limits on the boson mass but also break the degeneracy with WDM, and provide a unique fingerprint of interference in cosmology. In this paper, we provide initial steps towards this goal. In particular, we show in a bottom-up approach, how the presence of interference in an idealised filament population can lead to a non-suppressive feature in the matter power spectrum -- an observation supported by fully-cosmological FDM simulations. To this end, we build on a theoretically motivated and numerically observed steady-state approximation for filaments and express the equilibrium dynamics of such in an expansion of FDM eigenstates. We optimise the size of the expansion by incorporating classical phase-space information. Ellipsoidal collapse considerations are used to construct a fuzzy filament mass function which, together with the reconstructed FDM wave function, allow us to efficiently compute the one-filament power spectrum. We showcase our non-perturbative interference model for a selection of boson masses and confirm our approach is able to produce the matter power boost observed in fully-cosmological FDM simulations. More precisely, we find an excess in correlation between the spatial scale associated with the FDM ground state and the quantum pressure scale. We speculate about applications of this effect in data analysis.

The development of technologies for creating various types of solid-state detectors for optical astronomy is reviewed. The principles of designing astronomical photodetecting systems with large-format sensors based on charge-coupled device (CCD) and complementary metal oxide semiconductor (CMOS) structures are analyzed. Examples of the most advanced projects to which they have been applied are given. The history of the creation of optical detectors for telescopes operated in Russia is described, and a brief description and characteristics of the developed systems are provided. The results of testing in real research are displayed. The prospects for creating large-format systems based on CCD and CMOS detectors manufactured in Russia and abroad are discussed.

We aim to investigate profile changes at the cyclotron line energy of the accreting X-ray pulsar V 0332+53 by means of the analysis of its energy-resolved pulse profile behaviour, using the full set of available NuSTAR observations. We apply a tailored pipeline to study the energy dependence of the pulse profiles and to build the pulsed fraction spectra (PFS) for the different observations. We study the profile changes also using cross-correlation and lag spectra. We re-analyse the energy spectra to search for links between the local features observed in the PFS and spectral emission components associated with the shape of the fundamental cyclotron line. In the PFS data, with sufficiently high statistics, we observe a consistent behaviour around the cyclotron line energy. Specifically, two Gaussian-shaped features appear symmetrically on either side of the putative cyclotron line. These features exhibit minimal variation with source luminosity, and their peak positions consistently remain on the left and right of the cyclotron line energy. We interpret these features as evidence for cyclotron emission wings (also referred to as shoulders), as predicted by theoretical models of line formation for resonant cyclotron absorption and its propagation along the observer's line of sight. A phase-resolved analysis of the pulse in the energy bands surrounding these features enables us to determine both the spectral shape and the intensity of the photons responsible for these peaks in the PFS. Assuming these features correspond to a spectral component, we use their shapes as priors for the corresponding emission components finding a statistically satisfactorily description of the spectra. To explain these results, we propose that our line of sight is close to the direction of the spin axis, while the magnetic axis is likely orthogonal to it.

We analyze the radio emission from the $\beta$ Cep star V2187 Cyg using archive data from the Jansky Very Large Array. The observations were made in ten epochs at 1.39 and 4.96 GHz in the highest angular resolution A configuration. We determine a spectral index of of $\alpha = 0.6\pm0.2$ ($S_{\nu} \propto \nu^\alpha$), consistent with an ionized wind or a partially optically-thick synchrotron or gyrosynchrotron source. The emission is spatially unresolved at both frequencies. The 4.96 GHz data shows a radio pulse with a duration of about one month that can be modeled in terms of an internal shock in the stellar wind produced by a sudden increase in the mass-loss rate and the terminal velocity. The quiescent radio emission of V2187 Cyg at 4.96 GHz (with a flux density of $\simeq 150~\mu Jy$), cannot be explained in terms of an internally (by V2187 Cyg) or externally (by a nearby O star) photoionized wind. We conclude that, despite the spectral index suggestive of free-free emission from an ionized wind, the radio emission of V2187 Cyg most likely has a magnetic origin, a possibility that can be tested with a sensitive search for circular polarization in the radio, as expected from gyro-synchrotron radiation, and also by trying to measure the stellar magnetic field, that is expected to be in the range of several kGauss.

We present JWST-NIRCam narrow-band, 4.05 $\mu$m Brackett-$\alpha$ images of the Sgr C HII region, located in the Central Molecular Zone (CMZ) of the Galaxy. Unlike any HII region in the Solar vicinity, the Sgr C plasma is dominated by filamentary structure in both Brackett-$\alpha$ and the radio continuum. Some bright filaments, which form a fractured arc with a radius of about 1.85 pc centered on the Sgr C star-forming molecular clump, likely trace ionization fronts. The brightest filaments form a `$\pi$-shaped' structure in the center of the HII region. Fainter filaments radiate away from the surface of the Sgr C molecular cloud. The filaments are emitting optically thin free-free emission, as revealed by spectral index measurements from 1.28 GHz (MeerKAT) to 97 GHz (ALMA). But, the negative in-band 1 to 2 GHz spectral index in the MeerKAT data alone reveals the presence of a non-thermal component across the entire Sgr C HII region. We argue that the plasma flow in Sgr C is controlled by magnetic fields, which confine the plasma to rope-like filaments or sheets. This results in the measured non-thermal component of low-frequency radio emission plasma, as well as a plasma $\beta$ (thermal pressure divided by magnetic pressure) below 1, even in the densest regions. We speculate that all mature HII regions in the CMZ, and galactic nuclei in general, evolve in a magnetically dominated, low plasma $\beta$ regime.

The gamma-ray emission in blazars can be attributed to the leptonic Synchrotron Self-Compton (SSC) model, photo-hadronic interactions, or a combination thereof. While evidence supports both models, their specific contributions remain uncertain. One supportive piece of evidence for the SSC model is the correlation between synchrotron and SSC fluxes in some blazar's Spectral Energy Distribution (SED), indicating the relative contributions of leptonic and hadronic mechanisms. Observational studies of the HBL blazar Markarian 421 over several years, spanning TeV gamma rays and X-rays, have reported a linear correlation across various timescales, which breaks at the highest gamma-ray fluxes. Extending this analysis to four High synchrotron peaked BL Lac (HBL) blazars -- Markarian 501, 1ES 1959+650, PKS 2155-304 and 1ES 2344+514.-- we utilize multiwavelength data from ground-based Imaging Atmospheric Cherenkov Telescopes (IACTs) for gamma rays and satellite observations for X-rays. Our long-term study confirms a linear correlation between fluxes across these energy bands, except for Markarian 501, which shows a correlation index of $1.45 \pm 0.01$. Notably, the exceptional flaring episode of PKS 2155-304 exhibits a correlation index of 2 with extreme values of gamma-ray fluxes. We observe outliers with high gamma-ray fluxes, suggesting the involvement of another mechanism, either of hadronic or leptonic origin. Finally, all other correlations exhibit alignment with a general correlation, suggesting a common acceleration mechanism among them with slight variations likely due to individual magnetic field strengths.

The standard model for planet formation is a bottom-up process in which the origin of rocky and gaseous planets can be traced back to the collision of micron-sized dust grains within the gas-rich environment of protoplanetary disks. Key milestones along the way include disk formation, grain growth, planetesimal formation, core growth, gas accretion, and planetary system evolution. I provide an introductory overview of planet formation, emphasizing the main ideas and reviewing current theoretical understanding. Many of the phases of planet formation have a well-developed physical understanding, though the complexity of the problem means that few can be quantitatively modeled with complete confidence. Transformative advances in disk imaging and exoplanet detection provide the first direct information on the initial conditions for planet formation, and have motivated new formation models that are faster, more efficient, and lead to a more diverse set of architectures than their Solar System inspired forebears. Much remains to be learned, and I close with a personal, incomplete list, of open problems.

IGR J17591--2342, a transient accretion-powered millisecond X-ray pulsar, was discovered during its 2018 outburst. Here, we present a timing and spectral analysis of the source using {\it AstroSat} data of the same outburst. From the timing analysis, we obtain updated values of binary orbital parameters, which reveal an average pulsar spin frequency of 527.4256984(8) Hz. The pulse profiles can be fit well with four harmonically related sinusoidal components with fractional amplitudes of fundamental and second, third, and fourth harmonics as $\sim13$\%, $\sim$6\%, $\sim$0.9\%, $\sim$0.2\%, respectively. The energy-dependent study of pulse profiles in the range of $3-20$ keV shows that the fractional amplitude of both the fundamental and first overtone is consistent with being constant across the considered energy band. Besides, a decaying trend has been observed for both the fundamental and first overtone in the phase-delay versus energy relation resulting in soft X-ray (2.8-3.3 keV) phase lags of $\sim$0.05 and $\sim$0.13 with respect to $\leq 15$ keV photons, for the fundamental and first overtone, respectively. The combined spectra from the Large Area X-ray Proportional Counters and the Soft X-ray Telescope aboard {\it AstroSat} in the $1-18$ keV range can be fit well with an absorbed model consisting of a Comptonization, a blackbody and a Gaussian emission line component yielding as best-fit parameters a blackbody seed photon temperature $kT_{\rm bb}$ $\sim 0.95 \pm 0.03$ keV, and an electron temperature $kT_{\rm e}$ $\sim 1.54 \pm0.03$ keV. The spectral aspects suggest the scattering of photons from the accretion disc or the neutron star's surface.

A core sample of 59 unobscured type 1 AGNs with simultaneous XMM-Newton X-ray and UV observations is compiled from archive to probe the nature of soft X-ray excess (SE). In the first paper of this series, our focus centers on scrutinizing the spectral profile of the soft excess. Of the sources, $\approx$ 71% (42/59) exhibit powerlaw-like (po-like) soft excess, while $\approx$ 29% (17/59) exhibit blackbody-like (bb-like) soft excess. We show a cut-off powerlaw could uniformly characterize both types of soft excesses, with median Ecut of 1.40 keV for po-like and 0.14 keV for bb-like. For the first time, we report a robust and quantitative correlation between the SE profile and SE strength (the ratio of SE luminosity to that of the primary powerlaw continuum in 0.5 - 2.0 keV), indicating that stronger soft excess is more likely to be po-like, or effectively has a higher Ecut. This correlation cannot be explained by ionized disk reflection alone, which produces mostly bb-like soft excess (Ecut $\sim$ 0.1 keV) as revealed by relxilllp simulation. Remarkably, we show with simulations that a toy hybrid scenario, where both ionized disk reflection (relxilllp, with all reflection parameters fixed at default values except for ionization of the disk) and warm corona (compTT, with temperature fixed at 1 keV) contribute to the observed soft excess, can successfully reproduce the observed correlation. This highlights the ubiquitous hybrid nature of the soft X-ray excess in AGNs, and underscores the importance of considering both components while fitting the spectra of soft excess.

We continue to present the results of a Byurakan Narrow Band Imaging Survey (BNBIS). In this work we present the results of the search and further detailed investigation of the objects, found in the course of the BNBIS survey in the southern part of the Mon R2 association. For the search of HH objects the narrow band images, obtained with the 1-m Schmidt telescope of the Byurakan Observatory, were used. Newly found objects were imaged in optical and near-IR range with the Apache Point Observatory 3.5 meter telescope, and observed spectrally with long-slit spectrograph and scanning Fabry-Perot interferometer on 6 m telescope of Special Astrophysical Observatory of the Russian Academy of Sciences using SCORPIO-2. We found three new HH groups: HH 1233, HH 1234 and HH 1235, two of them represent extended collimated flows. HH 1233 is the C-shape bipolar outflow system associated with the 2MASS 06084223$-$0657385 source star. HH 1234 is the helical chain of HH knots near the star V963 Mon. HH 1235 is a separate compact knot, connected with the visible only in mid- and far-IR source WISE J060856.57$-$070103.5. We found also several molecular hydrogen outflows, one of which coincides with HH 1233 and two other are associated with the deeply embedded IR sources in the same field. One more probable bipolar H$_2$ outflow is related to WISE J060856.57$-$070103.5. The emission spectra and spectral energy distributions of the source stars were analyzed. According to them they should be under rather early evolutional stage.

Globular clusters host large populations of millisecond pulsars (MSPs) due to their high gravitational encounter rates, producing many binary systems and thus MSPs via the recycling process. Seven pulsars with spin periods ranging from 3 ms to 134 ms have been discovered in Terzan 1, which was targeted for pulsar searches with the Green Bank Telescope after Australia Telescope Compact Array imaging revealed steep-spectrum point sources in the cluster core. We have obtained timing observations over seven years, for the first seven Green Bank Telescope (GBT) discoveries (Terzan 1 A through G), using the GBT and Murriyang, CSIRO's Parkes radio telescope. All seven pulsars are isolated, consistent with Terzan 1's classification as a core-collapsed cluster (core collapse is predicted to disrupt, or ionize, binaries). With these timing solutions, we measured the positions and observed period derivatives, dP/dt, for each pulsar. The measured dP/dt values are composed of intrinsic spin-down and accelerations experienced by the pulsars (primarily from the cluster's gravitational potential), and they can be used to infer line-of-sight accelerations. We attempted to constrain the radius and density of the cluster core using these inferred accelerations. A wide range of radii and densities are possible, pointing to the need for continued timing as well as new discoveries to better constrain these cluster properties. We additionally find that Ter 1 A may be younger than the cluster and thus may have formed via a formation channel other than a core-collapse supernova. Theoretical formation mechanisms such as electron-capture supernovae from accretion- or merger-induced collapse of white dwarfs could potentially explain these pulsars' origins. It may therefore be a member of a small but growing class of globular cluster pulsars that appear to be significantly younger than their host clusters.

We present the first $\textit{NuSTAR}$ X-ray observation of EF Eri, a well-known polar system. The $\textit{NuSTAR}$ observation was conducted in conjunction with $\textit{NICER}$, shortly after EF Eri entered a high accretion state following an unprecedented period of low activity lasting 26 years since 1997. $\textit{NuSTAR}$ detected hard X-ray emission up to 50 keV with an X-ray flux of $1.2\times10^{-10}$ ergs s$^{-1}$ cm$^{-2}$ ($3\rm{-}50 keV$). Folded X-ray lightcurves exhibit a single peak with $\sim65\%$ spin modulation throughout the $3\rm{-}50$ keV band. We found no evidence of QPO signals at $\nu = 0.1\rm{-}100$ Hz with an upper limit on the QPO amplitude below $5\%$ ($90\%$ CL) at $\nu \sim 0.5$ Hz where the optical QPO was previously detected. Our 1-D accretion column model, called ${\tt MCVSPEC}$, was fitted to the $\textit{NuSTAR}$ spectral data, yielding an accurate WD mass measurement of $M = (0.55\rm{-}0.58) M_\odot$. $\texttt{MCVSPEC}$ accounts for radiative cooling by thermal bremsstrahlung and cyclotron emission, X-ray reflection off the WD surface, and a previously constrained range of the accretion column area. The derived WD mass range is in excellent agreement with the previous measurement of $M = (0.55\rm{-}0.60) M_\odot$ in the optical band. This demonstrates a combination of broadband X-ray spectral analysis and the ${\tt MCVSPEC}$ model that can be employed in our ongoing $\textit{NuSTAR}$ observation campaign of other polars to determine their WD masses accurately.

The question of what is the total entropy of the universe, how it compares to the maximal entropy of de Sitter space, and how it is distributed across the universe's components, bears considerable importance for a number of reasons. Here, we first update the computation of the entropy associated with various sectors of the observed universe, including in the diffuse cosmic and late-time gamma-ray and neutrino backgrounds, in baryonic matter both in diffuse components, in stars and stellar remnants, and in cosmic rays; we then update, crucially, the estimate of entropy in stellar-mass and super-massive black holes, whose abundance and mass function has come into increasingly sharp definition with recent observations and with the rapidly growing statistics of black-hole-black-hole mergers observed with gravity wave detectors. We also provide a new, corrected estimate of the potential entropy associated with a stochastic gravitational wave background, with dark sector radiations, and with several dark matter models. Finally, we utilize the similarly recently updated constraints on the abundance of hypothetical primordial black holes - black holes, that is, of non-stellar origin - to assess the maximal amount entropy they could store. We find that if supermassive primordial black holes exist, they can dominate the entropy budget of the universe consistently with current constraints on their abundance and mass function, to a level potentially not distant from the posited entropy associated with the cosmic event horizon of de Sitter spacetime. The same conclusion holds for certain dark sector models featuring a large number of dark degrees of freedom.

Active galactic nuclei (AGN) are powerful sources of panchromatic radiation. All AGN emit in X-rays, contributing around $\sim 5-10\%$ of the AGN bolometric luminosity. The X-ray emitting region, popularly known as the corona, is geometrically and radiatively compact with a size typically $\lesssim 10 \, R_{\rm G}$ (gravitational radii). The rapid and extreme variability in X-rays also suggest that the corona must be a dynamic structure. Decades of X-ray studies have shed much light on the topic, but the nature and origin of AGN corona are still not clearly understood. This is mostly due to the complexities involved in several physical processes at play in the high-gravity, high-density and high-temperature region in the vicinity of the supermassive black hole (SMBH). It is still not clear how exactly the corona is energetically and physically sustained near a SMBH. Since corona is ubiquitously found in AGN, is there something fundamental about the accretion process that produces it? In this review we discuss the X-ray observational properties of the corona in radio quiet AGN.

Detecting a pulsar associated with a supernova remnant (SNR) and/or pulsar wind nebula (PWN) is crucial for unraveling its formation history and pulsar wind dynamics, yet the association with a radio pulsar is observed only in a small fraction of known SNRs and PWNe. In this paper, we report the discovery of a young pulsar J1631$-$4722, associated with the Galactic SNR G336.7$+$0.5 using Murriyang, CSIRO's Parkes radio telescope. It is also potentially associated with a PWN revealed by the Rapid ASKAP (Australian Square Kilometre Array Pathfinder) Continuum Survey (RACS). This 118 ms pulsar has a high dispersion measure of 873 $\mathrm{pc\,cm^{-3}}$ and a rotation measure of $-$1004 $\mathrm{rad\,m^{-2}}$. Because of such a high DM, at frequencies below 2 GHz, the pulse profile is significantly scattered, making it effectively undetectable in previous pulsar surveys at $\sim$1.4 GHz. Follow-up observations yield a period derivative of $\dot{P} = 3.6 \times 10^{-15}$, implying a characteristic age, $\tau_{c} = 33\,$kyr, and spin-down luminosity, $\dot{E} = 1.3\times10^{36}\,$erg$\,s^{-1}$. PSR$\,$J1631$-$4722, with its high spin-down luminosity and potential link to a PWN, stands out as a promising source of the high-energy $\gamma$-ray emission observed in the region.

There has been a rapid increase in the known fast radio burst (FRB) population, yet the progenitor(s) of these events have remained an enigma. A small number of FRBs have displayed some level of quasi-periodicity in their burst profile, which can be used to constrain their plausible progenitors. However, these studies suffer from the lack of polarisation data which can greatly assist in constraining possible FRB progenitors and environments. Here we report on the detection and characterisation of FRB 20230708A by the Australian Square Kilometre Array Pathfinder (ASKAP), a burst which displays a rich temporal and polarimetric morphology. We model the burst time series to test for the presence of periodicity, scattering and scintillation. We find a potential period of T = 7.267 ms within the burst, but with a low statistical significance of 1.77$\sigma$. Additionally, we model the burst's time- and frequency-dependent polarisation to search for the presence of (relativistic and non-relativistic) propagation effects. We find no evidence to suggest that the high circular polarisation seen in FRB 20230708A is generated by Faraday conversion. The majority of the properties of FRB 20230708A are broadly consistent with a (non-millisecond) magnetar model in which the quasi-periodic morphology results from microstructure in the beamed emission, but other explanations are not excluded.

We report a rare case where an elliptical radio-loud quasar host, 3C 59, rejuvenates star formation activity through minor mergers with its nearby satellite galaxies. The inferred star formation history of 3C 59 shows significant star formation rejuvenation within the past 500 Myr, before which remains rather quiescent for most of the cosmic time. Three nearest satellite galaxies of 3C 59 exhibit significant morphological disturbances, and two of them present strong tidal tails pointing towards 3C 59. In addition, all the satellite galaxies within a projected distance of 200 kpc show low star formation activities. They also have systematically lower effective radius ($R_{\rm e}$) than local late-type galaxies, while 3C 59 has significantly larger $R_{\rm e}$ than both early- and late-type galaxies. All these features suggest that ongoing minor mergers between 3C 59 and its nearby satellites could be causing gas to flow into 3C 59, which induces the star formation rejuvenation, and possibly also triggers the quasar activity. The enormous power from the large-scale radio jet of 3C 59 may in turn help keep the halo hot, prevent gas cooling, and further reduce star formation in its satellite galaxies. These results provide important insights into the mass and size growth of central galaxies and star formation quenching of satellite galaxies in galaxy groups.

We present measurements of the neutral atomic hydrogen (HI) mass function (HIMF) and cosmic HI density ($\Omega_{\rm HI}$) at $0 \leq z \leq 0.088$ from the Looking at the Distant Universe with MeerKAT Array (LADUMA) survey. Using LADUMA Data Release 1 (DR1), we analyze the HIMF via a new "recovery matrix" (RM) method that we benchmark against a more traditional Modified Maximum Likelihood (MML) method. Our analysis, which implements a forward modeling approach, corrects for survey incompleteness and uses extensive synthetic source injections to ensure robust estimates of the HIMF parameters and their associated uncertainties. This new method tracks the recovery of sources in mass bins different from those in which they were injected and incorporates a Poisson likelihood in the forward modeling process, allowing it to correctly handle uncertainties in bins with few or no detections. The application of our analysis to a high-purity subsample of the LADUMA DR1 spectral line catalog in turn mitigates any possible biases that could result from the inconsistent treatment of synthetic and real sources. For the surveyed redshift range, the recovered Schechter function normalization, low-mass slope, and "knee" mass are $\phi_\ast = 3.56_{-1.92}^{+0.97} \times 10^{-3}$ Mpc$^{-3}$ dex$^{-1}$, $\alpha = -1.18_{-0.19}^{+0.08}$, and $\log(M_\ast/M_\odot) = 10.01_{-0.12}^{+0.31}$, respectively, which together imply a comoving cosmic HI density of $\Omega_{\rm HI}=3.09_{-0.47}^{+0.65}\times 10^{-4}$. Our results show consistency between RM and MML methods and with previous low-redshift studies, giving confidence that the cosmic volume probed by LADUMA, even at low redshifts, is not an outlier in terms of its HI content.

Event-based sensors detect only changes in brightness across a scene, each pixel producing an asynchronous stream of spatial-temporal data, rather than recording frames of overall illumination like a traditional frame-based sensor. This is advantageous for implementing into a wavefront sensor, which benefits from high temporal resolution and high dynamic range. The determination of tip-tilt in particular is still a problem in laser guide star adaptive optics as there is no current technological capabilities to measure it. This study characterised the behaviour of an event-based sensor in the context of tip-tilt sensing,investigating if the high temporal resolution of the event streams could address these challenges. Different conditions of tip-tilt and background illumination levels are explored and found to be a strong contender for tip-tilt sensing with laser guide stars.

We studied magnetically arrested disks (MAD) around rotating black holes (BH), under the influence of radiative cooling. We introduce a critical value of the mass accretion rate $\dot M_{\rm crit}$ for which the cooling by the synchrotron process efficiently radiates the thermal energy of the disk. We find $\dot M_{\rm crit} \approx 10^{-5.5} \dot M_{\rm Edd}$, where $\dot M_{\rm Edd}$ is the Eddington mass accretion rate. The normalization constant depends on the saturated magnetic flux and on the ratio of electron to proton temperatures, but not on the BH mass. We verify our analytical estimate using a suite of general relativistic magnetohydrodynamic (GRMHD) simulations for a range of black hole spin parameters $a \in \{ -0.94, -0.5, 0, 0.5, 0.94 \}$ and mass accretion rates ranging from $10^{-7}\dot M_{\rm Edd}$ to $10^{-4}\dot M_{\rm Edd}$. We numerically observe that the MAD parameter and the jet efficiency vary by a factor of $\approx 2$ as the mass accretion rate increases above $\dot M_{\rm crit}$, which confirms our analytical result. We further detail how the forces satisfying the quasi-equilibrium of the disk change, with the magnetic contribution increasing as the thermal contribution decreases.

We report the results of quasi-simultaneous multiwavelength (near-infrared, optical, UV, and X-ray) observations of the Galactic X-ray black hole binary MAXI J1820+070 performed in 2019 May 10-13, $\sim 60$ days after the onset of the first rebrightening phase. It showed a much larger optical-to-X-ray luminosity ratio ($\sim 8$) than in the initial outburst epoch. The primary components of the spectral energy distribution (SED) can be best interpreted by radiatively inefficient accretion flow (RIAF) spectrum showing a luminosity peak in the optical band. By comparison with theoretical calculations, we estimate the mass accretion rate to be $\dot{M}/(8 L_{\rm Edd}/c^2) \sim 10^{-3}$, where $c$ is the light speed and $L_{\rm Edd}$ is the Eddington luminosity. In addition to the RIAF emission, a blue power-law component is detected in the optical-UV SED, which is most likely synchrotron radiation from the jet. The optical spectrum taken at the Seimei telescope shows a weak and narrow H$\alpha$ emission line, whose emitting region is constrained to be $\gtrsim 2 \times 10^{4}$ times the gravitational radius. We suggest that the entire disk structure cannot be described by a single RIAF solution but cooler material responsible for the H$\alpha$ emission must exist at the outermost region.

To understand the evolution of global accretion disk structure in the ``rebrightening'' phase of MAXI J1820$+$070, we perform a comprehensive analysis of its near infrared/optical/UV to X-ray spectral energy distribution (SED) utilizing data obtained by OISTER, Las Cumbres Observatory (LCO), Swift, NICER, and NuSTAR in 2019. Optical spectra observed with Seimei telescope in 2019 and 2020 are also analyzed. On the basis of the optical and X-ray light curves and their flux ratios, we divide the whole phase into 3 periods, Periods I (flux rise), II (decay), and III (dim). In the first 2 periods, the source stayed in the low/hard state (LHS), where the X-ray (0.3--30 keV) and optical/UV SED can be both fitted with power-law models. We interpret that the X-ray emission arises from hot corona via Comptonization, whereas the optical/UV flux is dominated by synchrotron radiation from the jets, with a partial contribution from the irradiated disk. The optical/UV power-law component smoothly connects to a simultaneous radio flux, supporting its jet origin. Balmer line profiles in the optical spectra indicate that the inner radius of an irradiated disk slightly decreased from $\sim 2\times 10^5 r_{\rm g}$ (Period I) to $\sim 1\times 10^5 r_{\rm g}$ (Period II), where $r_{\rm g}$ is the gravitational radius, implying a change of the hot corona geometry. In Period III, the SED can be reproduced by an advection-dominated accretion flow and jet emission. However, the double-peaked H$\alpha$ emission line indicates that a cool disk remained at large radii.

The Reuven Ramaty High Energy Solar Spectrocopy Imager (RHESSI) $\gamma$-ray observations of the extraordinary GOES X25 flare SOL2003-10-28T11:10 are revisited to investigate previously reported conclusions that flare-accelerated electrons and protons precipitate along spatially separated flare loops. In contrast to previous works which reconstructed 2.223 MeV images over extended time periods ($\sim$20 minutes), we selected shorter integration times of the order of 2 to 3 minutes. Using simulations of the 2.223 MeV profile by Murphy et al. (2003) in combination with observations of the prompt $\gamma$-ray lines from the INTEGRAL mission, we obtain two separated integration time ranges representing the peak of the flare and the start of the decay, respectively. The resulting $\gamma$-ray images are then compared to GONG white-light (WL) observations to identify where along the flaring ribbons electrons and protons precipitation occurs. We point out that previously reported results comparing RHESSI hard X-ray (HXR) and $\gamma$-ray images only hold if the relative time evolution in the two energy ranges is the same. As the decay times for the 28 October 2003 is different at the considered two energy ranges (200-300 keV and around 2.223 MeV), the previously published conclusion that electrons and protons precipitate at different locations is an overstatement. Using shorter integration times reveals that the $\gamma$-ray and HXR sources spatially coincide with the WL flare ribbons. Our key conclusion is that electron and proton precipitation sites coincide with the flare ribbons, suggesting that the electron and proton precipitation sites are the same, at least within RHESSI's imaging capabilities. This result solves the twenty-years-long mystery around the previously reported different electron and proton precipitation sites.

Observations of protoplanetary disks have revealed the presence of both crescent-shaped and ring-like structures in dust continuum emission. These crescents are thought to arise from dust-trapping vortices generated by the Rossby Wave Instability (RWI), which induces density waves akin to those caused by planets. These vortices have the potential to create gaps and rings within the disk, resulting from the dissipation of their density waves. We carry out 2D hydrodynamic simulations in the shearing box to investigate vortex-disk interaction. We find that long-lived vortices can produce dust rings and gaps in inviscid discs detectable by ALMA, and a more elongated vortex produces rings at larger separations. Vortex-induced density waves carry over two orders of magnitude higher angular momentum flux compared to planet-induced ones that shock at the same location, making the former much more effective at producing dust gaps and rings far away.

The X-ray light curves of gamma-ray bursts (GRBs) display complex features, including plateaus and flares, that challenge theoretical models. Here, we study the properties of flares that are observed in the early afterglow phase (up to a few thousand of seconds). We split the sample into two groups: bursts with and without X-ray plateau. We find that the distributions of flare properties are similar in each group. Specifically, the peak time ($t_{\rm pk}$) of the flares and the ratio of the flare width to the flare peak time ($w/t_{\rm pk}$) which is found to be $\approx 1$, regardless of the presence of a plateau. We discuss these results in view of the different theoretical models aimed at explaining the origin of the plateau. These results are difficult to explain by viewing angle effects or late-time energy injection, but do not contradict the idea that GRBs with X-ray plateau have low Lorentz factor, of the order of tens. For these GRBs, the dissipation processes that produce the flares naturally occur at smaller radii compared to GRBs with higher Lorentz factors, while the flares maintain a similar behaviour. Our results therefore provide an independent support for the idea that many GRBs have a Lorentz factor of a few tens rather than a few hundreds.

Temperature-dependent biological productivity controls silicate weathering and thereby extends the potential habitable timespan of Earth. Models and theoretical considerations indicate that the runaway greenhouse on Earth-like exoplanets is generally accompanied by a dramatic increase in atmospheric H$_2$O and CO$_2$, which might be observed with the upcoming generation of space telescopes. If an active biosphere extends the habitable timespan of exoplanets similarly to Earth, observing the atmospheric spectra of exoplanets near the inner edge of the habitable zone could then give insights into whether the planet is inhabited. Here, we explore this idea for Earth-like stagnant-lid planets. We find that while for a reduced mantle, a surface biosphere extends the habitable timespan of the planet by about 1 Gyr, for more oxidising conditions, the biologically enhanced rate of weathering becomes increasingly compensated for by an increased supply rate of CO$_2$ to the atmosphere. Observationally, the resulting difference in atmospheric CO$_2$ near the inner edge of the habitable zone is clearly distinguishable between biotic planets with active weathering and abiotic planets that have experienced a runaway greenhouse. For an efficient hydrological cycle, the increased bioproductivity also leads to a CH$_4$ biosignature observable with JWST. As the planet becomes uninhabitable, the H$_2$O infrared absorption bands dominate, but the 4.3-micron CO$_2$ band remains a clear window into the CO$_2$ abundances. In summary, while the effect of life on the carbonate-silicate cycle leaves a record in the atmospheric spectrum of Earth-like stagnant-lid planets, future work is needed especially to determine the tectonic state and composition of exoplanets and to push forward the development of the next generation of space telescopes.

Context. The Stardust flyby mission to Jupiter-family comet (JFC) 81P/Wild 2 (hereafter, 81P) captured its dense quasicircular depressions. Nevertheless, the formation mechanism remains a subject of ongoing debate. Aims. We aimed to study how cometary activity contributed to the formation and enlargement of these depressions by analyzing Stardust flyby images and ground-based observation data. Methods. We calculated the time-dependent water production rate of 81P inside the snow line (<3 au) and compared it with the observational data. In addition, we estimated the fallback debris mass using an observation-based model, where a dust ejection from 81P was considered to reproduce ground-based observations of the dust tail. We compared the total excavated volume of water and dust with the total depression volume derived using the 81P shape model. Results. We find that the total excavated volume after 81P was injected into the inner Solar System accounts for up to only 30 % of the depression volume. This insufficiency suggests that a large portion (>70 %) of the depressions had already existed before the comet was injected into the current orbit. In addition, we estimate the dust-to-ice mass ratio for 81P to be 2-14. Conclusions. We suggest that most depressions were formed in the source region of the Kuiper-belt objects.

We review the main topics in the field of blazar multiwavelength variability as a tool to understand the physics and structure of extragalactic jets and their central engine. We address issues such as the cross-correlation between flux variations at different frequencies, the mechanisms to explain the long-term and short-term variability, the size of the jet emitting regions, the polarization behaviour, and the periodicities detected in the multiwavelength light curves.

A recent investigation highlighted peculiar trends between the radii derived from surface brightness-colour relations (SBCRs) combined with Gaia DR3 parallaxes with respect to asteroseismic scaling relation radii from K2 data. [...] We investigated on the robustness of the results based on Kepler data. We cross-matched asteroseismic and astrometric data for over 12,000 red giant branch and red clump stars from the end-of-mission Kepler catalogue with the Gaia DR3 and TIC v8.2 to obtain precise parallaxes, V- and K-band magnitudes, and E(B - V) colour excesses. Two well-tested SBCRs from the literature were adopted to estimate stellar radii. The analysis confirmed that SBCR and asteroseismic radii agree very well. The overall differences are only 1-2% depending on the adopted SBCR. The dispersion of 7% was about two-thirds of what was found for K2-based data. As a difference from the K2-based investigation, the ratio of SBCRs-to-asteroseismic radii did not depend on the metallicity [Fe/H]. Moreover, the intriguing decreasing trend with [$\alpha$/Fe] of the radius ratio for massive stars that was observed in K2 data was absent in Kepler data. The SBCR radii are systematically higher than asteroseismic estimates by 5% for stars with masses below 1.0 $M_{\odot}$. The SBCRs have proven to be a highly effective tool for estimating radii with a precision comparable to that obtained from asteroseismology, but at a significantly lower observational cost. Moreover, the superior concordance of Kepler-derived radii with SBCR measurements and the absence of the discrepancies observed in the K2-derived radii suggest the existence of underlying systematic errors that impact specific mass and metallicity regimes within the K2 dataset.

The Euclid telescope, launched from Cape Canaveral on July 1st, 2023, is dedicated to studying dark matter and dark energy from its orbit at the Sun-Earth Lagrangian point L2. It is equipped with two instruments: the visual imager (VIS) and the Near-Infrared Spectrometer and Photometer (NISP). The Euclid Wide Survey (Scaramella et al. 2022) will cover approximately 14500 deg2 of the extragalactic sky, and will be complemented by the Euclid Deep Survey that will cover about 40 deg2 and achieve two magnitudes fainter depth. Several ground-test campaigns were performed to assess NISP instrument basic functionalities and performances, some highlights will be reviewed in this thesis. The analysis comprises two key aspects: the evaluation of NISP image data acquired during ground-tests and during commissioning, and the verification of instrument performance in spectroscopic mode through the use of pixel-level simulated images. The analysis includes the characterisation of NISP detectors, from instrumental background to point-spread function (PSF) evaluations, confirming that NISP performance is well within the Euclid requirement. We introduce a calibration methodology to estimate the filter wheel assembly position from PSF distortions that have been successfully applied during the commissioning and will continue to be used during the mission for calibration purpose. Additionally, we investigated noise events that occur during image acquisition, with a particular focus on ŠsnowballsŠ and cosmic ray tracks. Finally, we present a spectroscopic pixel-level simulation campaign of the NISP instrument, where spectral energy distributions were generated and processed by the Euclid spectroscopic channel simulator.

Spider pulsars constitute a distinct subset within the domain of radio millisecond pulsars, divided further into the categories of black widows and redbacks. Evident across multiple wavelengths, these pulsars manifest periodic variations and reside within binary systems. Investigating and discovering additional spider-type pulsars carries significant implications for comprehending the evolution of high-mass stars. Particularly crucial is the validation of the "Recycling" theory of millisecond pulsar genesis. In this investigation, we systematically explore spider pulsar binary systems utilizing time-domain variability data from the Zwicky Transient Facility, in conjunction with Fermi unassociated gamma-ray sources sourced from the 4FGL-DR3 catalog. We have implemented a time-domain data processing pipeline utilizing the Lomb-Scargle Periodogram algorithm, integrated with the wget data crawling technology. This approach has led to the identification of 194 ellipsoidal variables and irradiation-type binary stars. Subsequent refinement through the Gaia Hertzsprung-Russell diagram has culled a selection of 24 spider pulsar gold sample candidates. By incorporating the 4FGL 95\% confidence error ellipse, the pool was narrowed down to 19 gold sample candidates. Utilizing the Gaia color-reduced proper motion diagram further refined the selection to 9 gold sample candidates. These newly identified spider pulsar candidates will inform subsequent observational campaigns across radio, X-ray, and optical spectroscopy, thereby facilitating a deeper validation of their physical characteristics.

The first JWST/MIRI photometric observations of TRAPPIST-1 b allowed for the detection of the thermal emission of the planet at 15 $\mu m$, suggesting that the planet could be a bare rock with a zero albedo and no redistribution of heat. These observations at 15 $\mu m$ were acquired as part of GTO time that included a twin program at 12.8 $\mu m$ in order to have a measurement in and outside the CO$_2$ absorption band. Here we present five new occultations of TRAPPIST-1 b observed with MIRI in an additional photometric band at 12.8 $\mu m$. We perform a global fit of the 10 eclipses and derive a planet-to-star flux ratio and 1-$\sigma$ error of 452 $\pm$ 86 ppm and 775 $\pm$ 90 ppm at 12.8 $\mu m$ and 15 $\mu m$, respectively. We find that two main scenarios emerge. An airless planet model with an unweathered (fresh) ultramafic surface, that could be indicative of relatively recent geological processes fits well the data. Alternatively, a thick, pure-CO2 atmosphere with photochemical hazes that create a temperature inversion and result in the CO2 feature being seen in emission also works, although with some caveats. Our results highlight the challenges in accurately determining a planet's atmospheric or surface nature solely from broadband filter measurements of its emission, but also point towards two very interesting scenarios that will be further investigated with the forthcoming phase curve of TRAPPIST-1 b.

The Solar Ultraviolet Imaging Telescope (SUIT) on board the Aditya-L1 mission is designed to observe the Sun across 200-400 nm wavelength. The telescope used 16 dichroic filters tuned at specific wavelengths in various combinations to achieve its science goals. For accurate measurements and interpretation, it is important to characterize these filters for spectral variations as a function of spatial location and tilt angle. Moreover, we also measured out-of-band and in-band transmission characteristics with respect to the inband transmissions. In this paper, we present the experimental setup, test methodology, and the analyzed results. Our findings reveal that the transmission properties of all filters meet the expected performance for spatial variation of transmission and the transmission band at a specific tilt angle. The out-of-band transmission for all filters is below 1% with respect to in-band, except for filters BB01 and NB01. These results confirm the capabilities of SUIT to effectively capture critical solar features in the anticipated layer of the solar atmosphere.

The structural information of spiral galaxies such as the spiral arm number, offer valuable insights into their formation processes and physical roles in galaxy evolution. We developed classifiers based on CNNs using variants of the EfficientNet architecture with different transfer learning techniques and pre-trained weights to categorise spiral galaxies by their number of spiral arms. A dataset from GZ2, comprising 11718 images filtered based on appropriate criteria is used for training and evaluation. Both the EfficientNetV2M model fine-tuned on ImageNet and the EfficientNetB0 model with Zoobot pre-trained weights achieved high accuracy on the down-sampled dataset, with most performance metrics exceeding 0.8 across all classes, except for galaxies with 4 arms due to the limited number of samples in this category. Merging higher-arm-number classes (more than 4 arms) improved the V2M model's accuracy significantly for 4-arm galaxies, as this approach allowed the model to focus on more distinct features with a more balanced class distribution. GradCAM++ and SmoothGrad highlight the networks' effectiveness in classifying galaxies, through the distinction of the galaxy structures and the extraction of the spiral arms, with the V2M model showing better capabilities in both tasks. Lower-arm galaxies tend to be misclassified as "can't tell" when their spiral arms are not clearly visible, while higher-arm galaxies tend to be misclassified as having fewer arms when their features are only partially detected. The study also found that galaxies with 3 arms tend to have lower stellar masses, and this tendency is reduced in the model predictions. The models' mispredictions between 2-arm and 1/3-arm are likely resulting from external interference and dynamic nature of spiral arms. The V2M model prediction also shows a slight tendency towards higher stellar mass in high-arm galaxies.

Subdwarf B stars are a well-known class of hot, low-mass stars thought to be formed through interactions in stellar binary systems. While different formation channels for subdwarf B stars have been studied through a binary population synthesis approach, it has also become evident that the characteristics of the found populations depend on the initial set of assumptions that describe the sometimes poorly constrained physical processes, such as common envelope episodes or angular momentum loss during mass transfer events. In this work we present a parameter study of subdwarf B populations, including a novel analytic prescription that approximates the evolution of subdwarf B stars with hydrogen-rich outer shells, an element previously overlooked in rapid binary population synthesis. We find that all studied parameters strongly impact the properties of the population, with the possibility of igniting helium below the expected core-mass value near the tip of the red giant branch strongly affecting the total number of subdwarf B candidates. Critically, our newly proposed prescription for the evolution of subdwarf B stars with hydrogen-shells helps to reconcile theoretical predictions of surface gravity and effective temperature with observational results. Our prescription is useful in the context of rapid binary population synthesis studies and can be applied to other rapid binary population synthesis codes' output.

We study the properties of 26 PNe with PG1159-type central stars known till date and compare them with the properties of PNe having [WR], $wels$ and hydrogen-rich central stars published earlier. We use archival photometric measurements of $2MASS$ for near-IR analysis and $WISE$ and $IRAS$ data for mid- and far-IR analysis and derive the IR properties of PG1159-PNe. We analyze the IR colour-colour diagrams of PG1159-PNe and compare them with the other three groups of PNe. Similar to the [WR]-PNe, many PG1159-PNe also show large amount of near-IR emission from the hot-dust component but their AGB dust is relatively cooler. We also report here the dust colour temperatures, dust masses, dust-to-gas mass ratios, IR luminosities and IR excess of PG1159-PNe and plot them against their surface H$\beta$ brightness (age) and compare them with the distribution of other groups of PNe. The IR luminosity and dust temperature show strong correlation with surface H$\beta$ brightness, however, the dust-to-gas mass ratio and IR excess do not show any trend. While the mean dust mass has a lower value for PG1159-PNe, in compared to other groups, the average dust-to-gas mass ratio is found to be marginally larger for PG1159-PNe. An analysis of the number distribution of different groups of PNe against surface H$\beta$ brightness shows that a) younger [WR]-, $wels$- and normal-PNe have a similar distribution indicating that they all have evolved from the AGB in a similar way, b) while there is an overlap of surface H$\beta$ brightness between [WR]- and PG1159-PNe, showing an evolutionary connection between them, there exists a significant gap between the values derived for $wels$- and PG1159-PNe.

We present $\sim8-40\,\mu$m SOFIA-FORCAST images of seven regions of ``clustered" star formation as part of the SOFIA Massive (SOMA) Star Formation Survey. We identify a total of 34 protostar candidates and build their spectral energy distributions (SEDs). We fit these SEDs with a grid of radiative transfer models based on the Turbulent Core Accretion (TCA) theory to derive key protostellar properties, including initial core mass, $M_c$, clump environment mass surface density, $\Sigma_{\rm cl}$, and current protostellar mass, $m_*$. We also carry out empirical graybody (GB) estimation of $\Sigma_{\rm cl}$, which allows a case of restricted SED fitting within the TCA model grid. We also release version 2.0 of the open-source Python package \emph{sedcreator}, designed to automate the aperture photometry and SED building and fitting process for sources in clustered environments, where flux contamination from close neighbors typically complicates the process. Using these updated methods, SED fitting yields values of $M_c\sim30-200\:M_{\odot}$, $\Sigma_{\text{cl,SED}}\sim0.1-3\:{\rm{g\:cm}}^{-2}$, and $m_*\sim4-50\:M_{\odot}$. The graybody fitting yields smaller values of $\Sigma_{\text{cl,GB}}\lesssim1\:{\rm{g\:cm}}^{-2}$. From these results, we do not find evidence for a critical $\Sigma_{\rm{cl}}$ needed to form massive ($\gtrsim 8\:M_\odot$) stars. However, we do find tentative evidence for a dearth of the most massive ($m_*\gtrsim30\:M_\odot$) protostars in the clustered regions suggesting a potential impact of environment on the stellar initial mass function.

Our current understanding is that intermediate- to high-mass stars form in a way similar to low-mass stars, that is, through disk accretion. However, the physical conditions that play a role in disk formation, evolution, and the possibility of (sub)stellar companion formation, are significantly different. We search for the mm counterparts of four intermediate- to high-mass (4-10 Solar mass) young stellar objects (YSOs) in the giant Hii region M17 at a distance of 1.7 kpc. These objects expose their photospheric spectrum such that their location on the pre-main-sequence (PMS) is well established. They have a circumstellar disk that is likely remnant of the formation process. With ALMA we have detected, for the first time, these four YSOs in M17, in Band 6 and 7, as well as four other serendipitous objects. Besides the flux measurements, the source size and spectral index provide important constraints on the physical mechanism(s) producing the observed emission. We apply different models to estimate the dust and gas mass contained in the disks. All our detections are spatially unresolved, constraining the source size to <120 au, and have a spectral index in the range 0.5-2.7. The derived (upper limits on the) disk dust masses are on the order of a few Earth masses and estimations of the upper limits on the gas mass vary between $10^{-5}$ and $10^{-3}$ Solar mass. In two objects (B331 and B268) free-free emission indicates the presence of ionized material around the star. The four serendipitous detections are likely (low-mass) YSOs. We compare the derived disk masses of our M17 targets to those obtained for YSOs in low-mass star-forming regions (SFRs) and Herbig stars, as a function of stellar mass, age, luminosity, and outer disk radius. The M17 sample, though small, is both the most massive and the youngest sample, yet has the lowest mean disk mass.

In modern astronomy, the quantity of data collected has vastly exceeded the capacity for manual analysis, necessitating the use of advanced artificial intelligence (AI) techniques to assist scientists with the most labor-intensive tasks. AI can optimize simulation codes where computational bottlenecks arise from the time required to generate forward models. One such example is PHOEBE, a modeling code for eclipsing binaries (EBs), where simulating individual systems is feasible, but analyzing observables for extensive parameter combinations is highly time-consuming. To address this, we present a fully connected feedforward artificial neural network (ANN) trained on a dataset of over one million synthetic light curves generated with PHOEBE. Optimization of the ANN architecture yielded a model with six hidden layers, each containing 512 nodes, provides an optimized balance between accuracy and computational complexity. Extensive testing enabled us to establish ANN's applicability limits and to quantify the systematic and statistical errors associated with using such networks for EB analysis. Our findings demonstrate the critical role of dilution effects in parameter estimation for EBs, and we outline methods to incorporate these effects in AI-based models. This proposed ANN framework enables a speedup of over four orders of magnitude compared to traditional methods, with systematic errors not exceeding 1\%, and often as low as 0.01\%, across the entire parameter space.

Chemical clocks based on [s-process elements/alpha-elements] ratios are widely used to estimate ages of Galactic stellar populations. However, the [s/alpha] vs. age relations are not universal, varying with metallicity, location in the Galactic disc, and specific s-process elements. Current Galactic chemical evolution models struggle to reproduce the observed [s/alpha] increase at young ages. We provide chemical evolution models for the Milky Way disc to identify the conditions required to reproduce the observed [s/H], [s/Fe], and [s/alpha] vs. age relations. We adopt a multi-zone chemical evolution model including state-of-the-art nucleosynthesis prescriptions for neutron-capture elements (AGB stars, rotating massive stars, neutron star mergers, magneto-driven supernovae). We explore variations in gas infall, AGB yield dependencies on progenitor stars, and rotational velocity distributions for massive stars. Results are compared with open cluster data from the Gaia-ESO survey. A three-infall scenario for disc formation captures the rise of [s/alpha] with age in the outer regions but fails in the inner ones, especially for second s-process peak elements. Ba production in the last 3 Gyr of chemical evolution would need to increase by half to match observations. S-process contributions from low-mass AGB stars improve predictions but require increases not supported by nucleosynthesis calculations, even with potential i-process contribution. Variations in the metallicity dependence of AGB yields show inconsistent effects across elements. Distributions of massive star rotational velocities fail to improve results due to balanced effects on elements. We confirm that there is no single relationship [s/alpha] vs. age, but that it varies along the MW disc. Current prescriptions for neutron-capture element yields cannot fully capture the complexity of evolution, particularly in the inner disc.

Gradual solar energetic particle (SEP) events are generally attributed to the particle acceleration in shock waves driven by coronal mass ejections (CMEs). Space-weather effects of such events are important, so there has been continuous effort to develop models able to forecast their various characteristics. Here we present the first version of a new such model with the primary goal to address energetic storm particle (ESP) events. The model, PARASOL, is built upon the PArticle Radiation Asset Directed at Interplanetary Space Exploration (PARADISE) test-particle simulation model of SEP transport, but includes a semi-analytical description of an inner (i.e., near the shock) part of the foreshock region. The semi-analytical foreshock description is constructed using simulations with the SOLar Particle Acceleration in Coronal Shocks (SOLPACS) model, which simulates proton acceleration self-consistently coupled with Alfven wave generation upstream of the shock, and subsequent fitting of the simulation results with suitable analytical functions. PARASOL requires input of solar wind and shock magnetohydrodynamic (MHD) parameters. We evaluate the performance of PARASOL by simulating the 12 July 2012 SEP event, using the EUropean Heliospheric FORecasting Information Asset (EUHFORIA) MHD simulation of the solar wind and CME in this event. The PARASOL simulation has reproduced the observed ESP event ($E \lesssim 5$ MeV) in the close vicinity of the shock within one order of magnitude in intensity.

We present a detailed investigation into the abundance and morphology of high redshift quenched galaxies at $3 < z < 7$ using James Webb Space Telescope data in the NEP, CEERS and JADES fields. Within these fields, we identify 90 candidate passive galaxies using specific star formation rates modelled with the BAGPIPES SED fitting code, which is more effective at identifying recently quenched systems than the classical UVJ method. With this sample of galaxies, we find number densities broadly consistent with other works and a rapidly evolving passive fraction of high mass galaxies ($\log_{10}{(M_{\star}/M_{\odot})} > $ 9.5) between $3 < z < 5$. We find that the fraction of galaxies with low star formation rates and mass 9.5 $ < \log_{10}{(M_{\star}/M_{\odot})} < $ 10.5 decreases from $\sim$25% at $3 < z < 4$ to $\sim$2% at $5 < z < 7$. Our passive sample of galaxies is shown to exhibit more compact light profiles compared to star-forming counterparts and some exhibit traces of AGN activity through detections in either the X-ray or radio. At the highest redshifts ($z > 6.5$) passive selections start to include examples of 'little red dots' which complicates any conclusions until their nature is better understood.

The recent announcement of evidence for a stochastic background of gravitational waves (GWB) in pulsar timing array (PTA) data has piqued interest across the scientific community. A combined analysis of all currently available data holds the promise of confirming the announced evidence as a solid detection of a GWB. However, the complexity of individual pulsar noise models and the variety of modeling tools used for different types of pulsars present significant challenges for a truly unified analysis. In this work we propose a novel approach to the analysis of PTA data: first a posterior distribution over Fourier modes is produced for each pulsar individually. Then, in a global analysis of all pulsars these posterior distributions can be re-used for a GWB search, which retains all information regarding the signals of interest without the added complexity of the underlying noise models or implementation differences. This approach facilitates combining radio and gamma-ray pulsar data, while reducing the complexity of the model and of its implementations when carrying out a GWB search with PTA data.

This study aims to test a potential application of lognormal seminumerical simulations to recover the thermal parameters and Jeans length. This could be suitable for generating large number of synthetic spectra with various input data and parameters, and thus ideal for interpreting the high-quality data obtained from QSO absorption spectra surveys. We use a seminumerical approach to simulate absorption spectra of quasars at redshifts $ 3 \leq z \leq 5$. These synthetic spectra are compared with the 1D flux power spectra and using the Markov Chain Monte Carlo analysis method we determine the temperature at mean density, slope of the temperature-density relation and Jeans length. Our best-fit model is also compared with the evolution of the temperature of the intergalactic medium from various UVB models. We show that the lognormal simulations can effectively recover thermal parameters and Jeans length. Besides, by comparing the synthetic flux power spectra with observations from Baryon Oscillation Spectroscopy Survey we found, that such an approach can be also used for the cosmological parameter inference.

Primordial Black Holes (PBHs) have not been experimentally detected so far, but their existence would provide important insights about the early Universe and serve as one of the possible candidates of dark matter (DM). In this work, we explore the accretion of radiation and matter by PBHs, with relevance for the growth of PBH seeds to form early Supermassive Black Holes; the emission from accreting PBHs; and constraints from gravitational wave observations, among others. We study the growth of PBH masses in the early Universe due to the accretion of radiation, highlighting uncertainties which arise from estimates of the PBH formation time. For baryonic accretion, we review the traditional Bondi-Hoyle-Lyttleton (BHL) and its refined version known as the Park-Ricotti (PR) model, which also includes radiative feedback. We find that in the BHL model, PBHs heavier than $\sim 100 \,\mathrm{M_{\odot}}$ can grow in mass by several orders of magnitude by $z \lesssim 10$, though only when surrounded by DM halos and only when the accretion efficiency is large. By contrast, the inclusion of radiation feedback in the PR model can drastically suppress the baryonic accretion rate of PBHs, leading to a negligible change in PBH mass over cosmic time. Furthermore our calculations show that the accretion rate depends sensitively on the modelling of various parameters such as the speed of sound in the baryonic gas and the velocity of PBHs. These findings highlight the uncertainties associated with accretion onto PBHs, and we find that a large increase in the PBH mass due to accretion is by no means guaranteed.

Ultrahigh-energy cosmic rays (UHECRs) are the highest energy messenger from space, with energies exceeding 1 EeV. Although UHECRs were discovered over 60 years ago, their origin still remains a mystery. Pinpointing sources of UHECRs is crucial for understanding the extreme astrophysical processes that accelerate particles to such extraordinary energies. We searched for UHECR multiplets via analyzing 17 years of data with energies greater than 40 EeV from the Pierre Auger Observatory. A spatial association is found between a multiplet of $25.7^{+6.2}_{-7.0}$ cosmic rays and the Sombrero galaxy with a local (global) significance of $4.5~\sigma~(3.3~\sigma)$. The Sombrero galaxy hosts a supermassive central black hole with a mass of $\sim1\times 10^9 M_{\odot}$ and exhibits large-scale radio lobes and jets. Our finding provides critical evidence on active supermassive black holes as the source of the highest-energy cosmic rays.

Atmospheric energy transport is central to the cooling of primordial magma oceans. Theoretical studies of atmospheres on lava planets have assumed that convection is the only process involved in setting the atmospheric temperature structure. This significantly influences the ability for a magma ocean to cool. It has been suggested that convective stability in these atmospheres could preclude permanent magma oceans. We develop a new 1D radiative-convective model in order to investigate when the atmospheres overlying magma oceans are convectively stable. Using a coupled interior-atmosphere framework, we simulate the early evolution of two terrestrial-mass exoplanets: TRAPPIST-1 c and HD 63433 d. Our simulations suggest that the atmosphere of HD 63433 d exhibits deep isothermal layers which are convectively stable. However, it is able to maintain a permanent magma ocean and an atmosphere depleted in H2O. It is possible to maintain permanent magma oceans underneath atmospheres without convection. Absorption features of CO2 and SO2 within synthetic emission spectra are associated with mantle redox state, meaning that future observations of HD 63433 d may provide constraints on the geochemical properties of a magma ocean analogous with the early Earth. Simulations of TRAPPIST-1 c indicate that it is expected to have solidified within 100 Myr, outgassing a thick atmosphere in the process. Cool isothermal stratospheres generated by low molecular-weight atmospheres can mimic the emission of an atmosphere-less body. Future work should consider how atmospheric escape and chemistry modulates the lifetime of magma oceans, and the role of tidal heating in sustaining atmospheric convection

We study the effects of hyperons, delta baryons, and quark matter phase transitions on f-mode oscillations in neutron stars. Using the density-dependent relativistic mean-field model (DDME2) for the hadronic phase and the density-dependent quark mass (DDQM) model for the quark phase, we construct hadronic and hybrid equations of state (EoSs) consistent with astrophysical constraints. Including hyperons and delta baryons soften the EoS, reducing maximum masse, while phase transition to the quark matter further softens the EoS, decreasing the speed of sound and hence the maximum mass. f-mode frequencies, calculated using both the Cowling approximation and the general relativistic (GR) frameworks, reveal a significant overestimation by the Cowling method of about 10-30%, with discrepancies decreasing for more massive stars. We derive universal relations connecting the frequencies of the f-mode to the average density, compactness, and tidal deformability, finding significant deviations due to hyperons and delta baryons. Empirical relations for mass-scaled and radius-scaled frequencies are also provided, highlighting the importance of GR calculations for accurate modeling. These findings highlight the potential of gravitational wave asteroseismology to constrain neutron star EoSs and internal structure.

Isolated neutron stars are thought to receive a natal kick velocity at birth nearly aligned with their spin axis. Direct observational confirmation of this alignment has been limited to a single source in a supernova remnant (PSR J0538+2817) whose three-dimensional velocity has been well-constrained. Pulsar polarisation statistical properties indicate the presence of a spin-kick correlation, but aligned and orthogonal cases remain plausible. However, if the three-dimensional velocities of radiopulsars are indeed predominantly aligned with their spin axes, a systematic difference in the observed transverse velocities of pulsars with small and large magnetic obliquities would be expected. In particular, due to projection effects, weakly oblique rotators should show systematically smaller and less scattered transverse velocities. In contrast, transverse velocities of pulsars with large obliquities should be close to their actual three-dimensional velocities. This study analyzed samples of 13 weakly and 25 strongly oblique pulsars with known distances and proper motions. We find their peculiar velocities being distributed differently with the statistical confidence of 0.007 and 0.016 according to Anderson-Darling and Kolmogorov-Smirnov tests, respectively. We performed a detailed population synthesis of the isolated pulsars, considering the evolution of their viewing geometry in both isotropic and spin-aligned kick scenarios. The observed split in the transverse velocity distributions and its amplitude are consistent with the spin-aligned kick model but not the isotropic case. At the same time, an orthogonal kick predicts a similar effect but of the opposite sign. This provides robust support for pulsar spin-kick alignment based on their statistics and independent of their polarization properties.

We study a class of homogeneous but anisotropic cosmologies within the family of shift-symmetric Horndeski theories, where the scalar field features an inhomogeneous profile but it preserves a translational symmetry that is realised as a combination of spatial translations and internal shifts. The spatial gradient of the scalar field introduces a preferred direction, so the resulting cosmologies are of the axisymmetric Bianchi I type. The momentum density of these configurations exhibits a universal evolution and an additional component with non-vanishing momentum density is required to have non-trivial effects. We show the relation of these scenarios with cosmologies of non-comoving components and, in particular, we explain how they provide a specific realisation of moving dark energy models. Among the class of shift-symmetric Horndeski theories, we analyse in more detail the case of Kinetic Gravity Braiding with emphasis on its application to moving dark energy models and its effects on large scale dark flows as well as the CMB dipole and quadrupole.

In this paper, solar cycles 21 to 24 were compared using complex network analysis. A network was constructed for these four solar cycles to facilitate the comparison. In these networks, the nodes represent the active regions of the Sun that emit flares, and the connections correspond to the sequence of solar flares over time. This resulted in a directed network with self-connections allowed. The model proposed by Abe and Suzuki for earthquake networks was followed. The incoming degree for each node was calculated, and the degree distribution was analyzed. It was found that for each solar cycle, the degree distribution follows a power law, indicating that solar flares tend to appear in correlated active zones rather than being evenly distributed. Additionally, a variation in the characteristic exponent {\gamma} for each cycle was observed, with higher values in even cycles compared to odd cycles. A more detailed analysis was performed by constructing 11-year networks and shifting them in one-year intervals. This revealed that the characteristic exponent shows a period of approximately 22 years coincident with the Hale cycle, suggesting that the complex networks provide information about the solar magnetic activity.

New radio interferometers such as MeerKAT, SKA, ngVLA, and DSA-2000 drive advancements in software for two key reasons. First, handling the vast data from these instruments requires subdivision and multi-node processing. Second, their improved sensitivity, achieved through better engineering and larger data volumes, demands new techniques to fully exploit it. This creates a critical challenge in radio astronomy software: pipelines must be optimized to process data efficiently, but unforeseen artefacts from increased sensitivity require ongoing development of new techniques. This leads to a trade-off among (1) performance, (2) flexibility, and (3) ease-of-development. Rigid designs often miss the full scope of the problem, while temporary research code is unsuitable for production. This work introduces a framework for developing radio astronomy techniques while balancing the above trade-offs. It prioritizes flexibility and ease-of-development alongside acceptable performance by leveraging Open Source data formats and software. To manage growing data volumes, data is distributed across multiple processors and nodes for parallel processing, utilizing HPC and cloud infrastructure. We present two Python libraries, Dask-MS and Codex Africanus, which enable distributed, high-performance radio astronomy software with Dask. Dask is a lightweight parallelization and distribution framework that integrates with the PyData ecosystem, addressing the "Big Data" challenges of radio astronomy.

Several groups have recently suggested that small planets orbiting very closely around white dwarf stars could be promising locations for life to arise, even after stellar death. There are still many uncertainties, however, regarding the existence and habitability of these worlds. Here, we consider the retention of water during post-main-sequence evolution of a Sun-like star, and during the subsequent migration of planets to the white dwarf's habitable zone. This inward migration is driven by dynamical mechanisms such as planet-planet interactions in packed systems, which can excite planets to high eccentricities, setting the initial conditions for tidal migration into short-period orbits. In order for water to persist on the surfaces of planets orbiting white dwarfs, the water must first survive the AGB phase of stellar evolution, then avoid being lost due to photoevaporation due to X-ray and extreme ultraviolet (XUV) radiation from the newly-formed white dwarf, and then finally survive the tidal migration of the planet inwards to the habitable zone. We find that while this journey will likely desiccate large swaths of post-main-sequence planetary systems, planets with substantial reservoirs of water may retain some surface water, especially if their migration occurs at later white dwarf cooling ages. Therefore, although stellar evolution may pose a challenge for the retention of water on exoplanet surfaces, it is possible for planets to retain surface oceans even as their host stars die and their orbits evolve.

The neglect of modelling both stellar and nebular emission significantly affects the derived physical properties of galaxies, particularly those with high star formation rates. While this issue has been studied, it has not been established a clear threshold for a significant impact on the estimated physical properties of galaxies due to accounting for both stellar and nebular emission. We analyse galaxies from SDSS-DR7 across a wide range of star-forming activity levels, comparing the results obtained from two spectral fitting tools: FADO (which considers both stellar and nebular continuum) and STARLIGHT (only considers the stellar continuum). A strong linear correlation is found between the rest-frame H$\alpha$ and H$\beta$ equivalent widths (EWs) and the optical nebular contribution, identifying these as reliable tracers. The results show that when the nebular contribution exceeds 8% (corresponding to EW(H$\alpha$)$\simeq$500 Å and EW(H$\beta$)$\simeq$110 Å), there is a significant impact on the estimation of galaxy properties, namely stellar mass, age and metallicity. Our results highlight the importance of taking into account both the stellar and nebular continuum when analysing the optical spectra of star-forming galaxies. In particular, this is a fundamental aspect for galaxies with a rest-frame EW(H$\alpha$)$\gtrsim$500 Å (or the scaled value of 375 Å for pseudo-continuum measures). At low redshifts, this mostly impacts extreme emission line galaxies, while at higher redshifts it becomes a dominant aspect given the higher star-forming activity in the younger Universe. In light of current JWST observations and future instruments designed for high-redshift observations, such as MOONS, this reveals as a critical issue to take into consideration.

We measure the cross-correlation between cosmic shear from the third-year release of the Dark Energy Survey, thermal Sunyaev-Zel'dovich (tSZ) maps from Planck, and X-ray maps from ROSAT. We investigate the possibility of developing a physical model able to jointly describe both measurements, simultaneously constraining the spatial distribution and thermodynamic properties of hot gas. We find that a relatively simple model is able to describe both sets of measurements and to make reasonably accurate predictions for other observables (the tSZ auto-correlation, its cross-correlation with X-rays, and tomographic measurements of the bias-weighted mean gas pressure). We show, however, that contamination from X-ray AGN, as well as the impact of non-thermal pressure support, must be incorporated in order to fully resolve tensions in parameter space between different data combinations. We obtain simultaneous constraints on the mass scale at which half of the gas content has been expelled from the halo, $\mathrm{log}_{10}(M_c)=14.83^{+0.16}_{-0.23}$, on the polytropic index of the gas, $\Gamma=1.144^{+0.016}_{-0.013}$, and on the ratio of the central gas temperature to the virial temperature $\alpha_T=1.30^{+0.15}_{-0.28}$.

The DarkNESS (Dark Matter Nano-satellite Equipped with Skipper Sensors) mission aims to deploy a skipper-CCD CubeSat Observatory to search for dark matter (DM) from Low Earth Orbit. This mission will employ novel skipper-CCDs to investigate O(keV) X-rays from decaying DM, as well as electron recoils from strongly-interacting sub-GeV DM. The DarkNESS mission will be the first space deployment of skipper-CCDs, and the DarkNESS team is developing a skipper-CCD instrument that is compatible with the CubeSat platform. DarkNESS has recently progressed from laboratory validation to a Critical Design Review (CDR) phase, with a launch opportunity anticipated in late 2025. The implementation of the DarkNESS skipper-CCD payload on the CubeSat platform will pave the way for future demonstrators of space-based imagers for X-ray and single-electron counting applications.

It has been suggested that merging black hole (BH) binaries in active galactic nucleus (AGN) discs formed through two-body scatterings via the gas-capture process may explain a significant fraction of BH mergers in AGN and a non-negligible contribution to the observed rate from LIGO-VIRGO-KAGRA. We perform Monte Carlo simulations of BH and binary BH formation, evolution and mergers across the observed AGN mass function using a novel physically motivated treatment for the gas-capture process derived from hydrodynamical simulations of BH-BH encounters in AGN and varying assumptions on the AGN disc physics. The results suggest that gas-captured binaries could result in merger rates of 0.73 - 7.1Gpc$^{-3}$yr$^{-1}$. Most mergers take place near the outer boundary of the accretion disk, but this may be subject to change when migration is considered. The BH merger rate in the AGN channel in the Universe is dominated by AGN with supermassive BH masses on the order of 10$^{7} M_\odot$ , with 90% of mergers occurring in the range 10$^{6} M_\odot$ - 10$^{8} M_\odot$ . The merging mass distribution is flatter than the initial BH mass power law by a factor $\Delta \xi$ = 1.1 to 1.2, as larger BHs can align with the disc and successfully form binaries more efficiently. Similarly, the merging mass ratio distribution is flatter, therefore the AGN channel could easily explain the high mass and unequal mass ratio detections such as GW190521 and GW190814. When modelling the BH binary formation process using a simpler dynamical friction treatment, we observe very similar results, where the primary bottleneck is the alignment time with the disk. We find the most influential parameters on the rates are the anticipated number of BHs and their mass function. We conclude that AGN remain an important channel for consideration, particularly for gravitational wave detections involving one or two high mass BHs.

We show that the Milgromian acceleration of MOND and the cosmological constant can be understood and quantified as the effects of quantum fluctuations of spin connection which are described by precanonical quantum gravity put forward by one of us earlier. We also show that a MOND-like modification of Newtonian dynamics at small accelerations emerges from this picture in the non-relativistic approximation.

The pulsar timing technique, which compares the observed arrival times of electromagnetic radiation from a pulsar with the predicted arrival times derived from a theoretical model of the pulsar system, is used in pulsar astronomy to infer a multitude of physical information and to constrain possible corrections to General Relativity (GR). The propagation delay is usually computed using formulas based on a post-Newtonian approach, for both the light trajectory and the orbital motion. However, evidence has recently emerged that this approximation may no longer be sufficient when the companion object is a supermassive black hole; deviations from a full GR computation of the propagation delay can reach a few seconds. In this paper, we analyze the case of binary pulsars with a stellar or intermediate black hole companion, whose discovery and timing are key goals of SKA. With a numerical algorithm, we have found that in this case, the full GR value depends only on the semi-major axis of the relative orbit and on the mass of the black hole companion. If the mass of the latter is sufficiently large ($100 M_{\odot}$), the maximum difference between the two approaches is significant ($\sim10^{-7}$ s) even for large binaries ($\sim10^{16}$ cm), and increases up to $\sim 10^{-4}$ s when the mass is $10^5 M_{\odot}$. We also consider relativistic corrections to the orbital motion, and discover that they can strongly affect the value of the propagation delay. We conclude that in the future, post-Newtonian formulas should be replaced with a more accurate approach in these systems, especially in view of future discoveries made by new large telescopes such as SKA.

Space weather modelling has been gaining importance due to our increasing dependency on technology sensitive to space weather effects, such as satellite services, air traffic and power grids. Improving the reliability, accuracy and numerical performance of space weather modelling tools, including global coronal models, is essential to develop timely and accurate forecasts and to help partly mitigate the space weather threat. Global corona models, however, require accurate boundary conditions, for the formulations of which we have very limited observational data. Unsuitable boundary condition prescriptions may lead to inconsistent features in the solution flow field and spoil the code's accuracy and performance. In this paper, we develop an adjustment to the inner boundary condition of the COCONUT global corona model to better capture the dynamics over and around the regions of stronger magnetic fields by constraining the plasma \b{eta} and the Alfvén speed. Using data from solar observations and solar atmospheric modelling codes such as Bifrost, we find that the baseline homogeneous boundary condition formulations for pressure and density do not capture the plasma conditions physically accurately. We develop a method to adjust these prescribed pressure and density values by placing constraints on the plasma \b{eta} and the Alfvén speed that act as proxies. We demonstrate that we can remove inexplicable fast streams from the solution by constraining the maximum Alfvén speed and the minimum plasma \b{eta} on the boundary surface. We also show that the magnetic topology is not significantly affected by this treatment otherwise. The presented technique shows the potential to ease the modelling of solar maxima, especially removing inexplicable features while, at the same time, not significantly affecting the magnetic field topology around the affected regions.

The Standard Model extended by a real scalar singlet $S$ with an approximate $\mathbb{Z}_2$ symmetry offers a minimal framework for realizing electroweak baryogenesis (EWBG) during a first-order electroweak phase transition. In this work, we explore a novel mechanism where spontaneous $\mathbb{Z}_2$ breaking enables EWBG via domain walls separating two distinct phases of the $S$ field. These domain walls feature restored (or weakly broken) EW symmetry in their cores and sweep through space, generating the baryon asymmetry below the temperature of EW symmetry breaking. We identify the key conditions for the existence of EW-symmetric domain wall cores and analyze the dynamics required for wall propagation over sufficient spatial volumes. Additionally, we outline the CP-violating sources necessary for baryogenesis under different regimes of domain wall evolution. The parameter space accommodating this mechanism spans singlet masses from sub-eV to 15 GeV, accompanied by a non-vanishing mixing with the Higgs boson. Unlike the standard realization of EWBG in the minimal singlet-extended SM, which is notoriously difficult to test, our scenario can be probed by a wide range of existing and upcoming experiments, including fifth force searches, rare meson decays, and EDM measurements.

This paper reports the first search for stellar-origin binary black holes within the LISA Data Challenges (LDC). The search algorithm and the Yorsh LDC datasets, both previously described elsewhere, are only summarized briefly; the primary focus here is to present the results of applying the search to the challenge of data. The search employs a hierarchical approach, leveraging semi-coherent matching of template waveforms to the data using a variable number of segments, combined with a particle swarm algorithm for parameter space exploration. The computational pipeline is accelerated using GPU hardware. The results of two searches using different models of the LISA response are presented. The most effective search finds all five sources in the data challenge with injected signal-to-noise ratios $\gtrsim 12$. Rapid parameter estimation is performed for these sources.

We propose a novel detection method for axion dark matter using the Rabi oscillation of neutron spins in beam-based measurements. If axions couple to neutron spins, a background oscillating axion dark matter field would drive transitions between spin-up and spin-down neutron states in a magnetic field when the axion particle energy matches the energy gap between the spin states. The transition can be detected in a double-Stern-Gerlach-type apparatus, with the first splitter producing a pure spin-polarized neutron beam and the second splitter selecting spin-flipped signals. Our approach offers enhanced detection capability for axions within the $10^{-12} - 10^{-10} \,$eV mass window with the capability to surpass the sensitivity of current laboratory experiments.

Compact binary sources that emit gravitational waves (GW) are expected to be both spinning and have eccentric orbits. To date, there has been no closed-form expression for the phasing of GWs that contain information from both spin and eccentricity. The introduction of eccentricity can slow waveform generation, as obtaining closed-form expressions for the waveform phase is unattainable due to the complexity of the differential equations involved, often requiring slower numerical methods. However, closed-form expressions for the waveform phase can be obtained when eccentricity is treated as a small parameter, enabling quick waveform generation. In this paper, closed-form expressions for the GW phasing in the form of Taylor approximants up to the eighth power in initial eccentricity $(e_0)$ are obtained while also including aligned spins up to the third post-Newtonian order. The phasing is obtained in both time and frequency domains. The TaylorT2 phasing is also resummed for usage in scenarios where initial eccentricities are as high as 0.5. Finally, a waveform is constructed using the $e_{0}^2$ expanded TaylorF2 phasing for aligned-spin systems added to TaylorF2Ecc. We perform mismatch computation between this model and TaylorF2Ecc. The findings indicate that for eccentricities $\gtrsim 0.15$ (defined at 10 Hz) and small spins $(\sim 0.2 )$, the mismatches can be higher than 1%. This leads to an overall loss in signal-to-noise ratio and lower detection efficiency of GWs coming from eccentric spinning compact binary inspirals if the combined effects of eccentricity and aligned spins are neglected in the waveforms.

Symmetric instability has broad applications in geophysical fluid dynamics. It plays a crucial role in the formation of mesoscale rainbands at mid-latitudes on Earth, instability in the ocean's mixed layer, and slantwise convection on gas giants and in the oceans of icy moons. Here, we apply linear instability analysis to an arbitrary zonally symmetric Boussinesq flow on a rotating spherical planet, with applicability to planetary atmospheres and icy moon oceans. We characterize the instabilities into three types: (1) gravitational instability, occurring when stratification is unstable along angular momentum surfaces, (2) inertial instability, occurring when angular momentum shear is unstable along buoyancy surfaces, and (3) a mixed ``PV'' instability, occurring when the potential vorticity has the opposite sign as planetary rotation. We note that $N^2>0$, where $N$ is the Brunt-Väisälä frequency, is neither necessary nor sufficient for stability. Instead, $b_z \sin{\theta}>0$, where $b_z$ is the stratification along the planetary rotation axis and $\theta$ is latitude, is always necessary for stability and also sufficient in the low Rossby number limit. In the low Rossby number limit, applicable to convection in the oceans of icy moons and in the atmospheres of gas giants, the most unstable mode is slantwise convection parallel to the planetary rotation axis.

In this work, we derive a generalized modified Friedmann equation based on an entropy-area relation that incorporates established modifications, such as volumetric, linear, and logarithmic terms, in addition to novel entropic modifications that might yield to relevant cosmological implications at different stages of the evolution of the Universe. Some of these modifications are capable of mimicking the effects of dark energy and describing the current state of accelerated expansion of the Universe. We study particular cases of the generalized Friedmann equation and constrain the free parameters using observational datasets, including Hubble parameter measurements, baryon acoustic oscillations, and strong lensing systems. Our findings indicate that the proposed models align well with current observational data, particularly in low-redshift regimes; furthermore, these models are compatible with the value of $H_0$ obtained by the SH0ES program.

Parameter estimation of gravitational wave data is often computationally expensive, requiring simplifying assumptions such as circularisation of binary orbits. Although, if included, the sub-dominant effects like orbital eccentricity may provide crucial insights into the formation channels of compact binary mergers. To address these challenges, we present a pipeline strategy leveraging minimally modelled waveform reconstruction to identify the presence of eccentricity in real time. Using injected signals, we demonstrate that ignoring eccentricity ($e_{\rm 20Hz} \gtrsim 0.1$) leads to significant biases in parameter recovery, including chirp mass estimates falling outside the 90% credible interval. Waveform reconstruction shows inconsistencies increase with eccentricity, and this behaviour is consistent for different mass ratios. Our method enables low-latency inferences of binary properties supporting targeted follow-up analyses and can be applied to identify any physical effect of measurable strength.

With the increasing use of Machine Learning (ML) algorithms in scientific research comes the need for reliable uncertainty quantification. When taking a measurement it is not enough to provide the result, we also have to declare how confident we are in the measurement. This is also true when the results are obtained from a ML algorithm, and arguably more so since the internal workings of ML algorithms are often less transparent compared to traditional statistical methods. Additionally, many ML algorithms do not provide uncertainty estimates and auxiliary algorithms must be applied. Conformal Prediction (CP) is a framework to provide such uncertainty quantifications for ML point predictors. In this paper, we explore the use and properties of CP applied in the context of glitch classification in gravitational wave astronomy. Specifically, we demonstrate the application of CP to the Gravity Spy glitch classification algorithm. CP makes use of a score function, a nonconformity measure, to convert an algorithm's heuristic notion of uncertainty to a rigorous uncertainty. We use the application on Gravity Spy to explore the performance of different nonconformity measures and optimise them for our application. Our results show that the optimal nonconformity measure depends on the specific application, as well as the metric used to quantify the performance.

We study fuzzy axion dark matter in type IIB string theory, for axions descending from the Ramond-Ramond four-form in compactifications on orientifolds of Calabi-Yau hypersurfaces. Such models can be tested by cosmological measurements if a significant relic abundance of fuzzy dark matter arises, which we argue is most common in models with small numbers of axions. We construct a topologically exhaustive ensemble of more than 350,000 Calabi-Yau compactifications yielding up to seven axions, and in this setting we perform a systematic analysis of misalignment production of fuzzy dark matter. In typical regions of moduli space, the fuzzy axion, the QCD axion, and other axions have comparable decay constants of $f_a \approx 10^{16}$ GeV. We find that overproduction of heavier axions is problematic, except at special loci in moduli space where decay constant hierarchies can occur: without a contrived reheating epoch, it is necessary to fine-tune initial displacements. The resulting dark matter is typically a mix of fuzzy axions and heavier axions, including the QCD axion. Dark photons are typically present as a consequence of the orientifold projection. We examine the signatures of these models by simulating halos with multiple fuzzy axions, and by computing new cosmological constraints on ultralight axions and dark radiation. We also give evidence that cosmic birefringence is possible in this setting. Our findings determine the phenomenological correlates of fuzzy axion dark matter in a corner of the landscape.

The Hellings-Downs (HD) correlation, which characterizes the signature of a stochastic gravitational wave background measured via Pulsar Timing Arrays (PTA), is derived using a harmonic formalism. This approach closely follows the theoretical framework traditionally employed to compute correlations of temperature fluctuations in the Cosmic Microwave Background (CMB). This parallel enables a direct comparison between the correlations observed in PTA and those in CMB experiments. In particular, we clarify that what has been referred to as cosmic variance in previous literature is not an intrinsic limitation for PTA measurements. Instead, when an optimal estimator is used, this variance can be mitigated by accumulating more observation time or improving the cadence of pulsar monitoring. Therefore, unlike CMB angular correlations, where cosmic variance represents an irreducible constraint, it can be reduced in PTA measurements and will continue to diminish in future experiments. Finally, we show that if the primordial power spectrum of tensor fluctuations was very blue with $n_T>4$, the CMB angular correlation due to these tensor modes would also exhibit a HD correlation.