Papers for Monday, Oct 30 2023

Tuning into the bass notes of the large-scale structure requires careful attention to geometrical effects arising from wide angles. The spherical Fourier-Bessel (SFB) basis provides a harmonic-space coordinate system that fully accounts for all wide-angle effects. To demonstrate the feasibility of the SFB power spectrum, in this paper we validate our SFB pipeline by applying it to lognormal, and both complete and realistic EZmock simulations that were generated for eBOSS DR16 LRG sample. We include redshift space distortions and the local average effect (aka integral constraint). The covariance matrix is obtained from 1000 EZmock simulations, and inverted using eigenvalue decomposition.

In this paper we present two examples of recent investigations that we have undertaken, applying Machine Learning (ML) neural networks (NN) to image datasets from outer planet missions to achieve feature recognition. Our first investigation was to recognize ice blocks (also known as rafts, plates, polygons) in the chaos regions of fractured ice on Europa. We used a transfer learning approach, adding and training new layers to an industry-standard Mask R-CNN (Region-based Convolutional Neural Network) to recognize labeled blocks in a training dataset. Subsequently, the updated model was tested against a new dataset, achieving 68% precision. In a different application, we applied the Mask R-CNN to recognize clouds on Titan, again through updated training followed by testing against new data, with a precision of 95% over 369 images. We evaluate the relative successes of our techniques and suggest how training and recognition could be further improved. The new approaches we have used for planetary datasets can further be applied to similar recognition tasks on other planets, including Earth. For imagery of outer planets in particular, the technique holds the possibility of greatly reducing the volume of returned data, via onboard identification of the most interesting image subsets, or by returning only differential data (images where changes have occurred) greatly enhancing the information content of the final data stream.

Among the sample of Galactic transient X-ray binaries (SXTs) discovered to date, ~70 have been proposed as candidates to host a black hole. Yet, only 19 have been dynamically confirmed. Such a reliable confirmation requires phase-resolved spectroscopy of their companion stars, generally feasible when the system is in a quiescent state. However, since most of the SXT population lies in the Galactic Plane, strongly affected by interstellar extinction, their optical brightness during quiescence usually falls beyond the capabilities of the current instrumentation ($R\gtrsim22$). To overcome these limitations and, thus, increase the number of confirmed Galactic black holes, a correlation between the full-width at half maximum (FWHM) of the H$\alpha$ line and the semi-amplitude of the donor's radial velocity curve ($K_2$) was presented in the past. Here, we extend the FWHM-$K_2$ correlation to the near-infrared (NIR), exploiting disc lines such as He I $\lambda$10830, Pa$\gamma$ and Br$\gamma$, in a sample of dynamically confirmed black-hole SXTs. We obtain $K_2 = 0.22(3) ~\textrm{FWHM}$, in good agreement with the optical correlation derived using H$\alpha$. The similarity of the two correlations seems to imply that the widths of H$\alpha$ and the NIR lines are consistent in quiescence. When combined with information on orbital periods, the NIR correlation will allow us to constrain the mass of the compact object of systems in quiescence by using single-epoch spectroscopy. We anticipate that this new correlation will give access to highly-reddened black-hole SXTs, which cannot be otherwise studied at optical wavelengths.

The chemical composition of the highest-energy cosmic rays, namely the atomic number $Z$ of rays with energies $E\gg1~\text{EeV}$, remains to date largely unknown. Some information on the composition can be inferred from the deflections that charged ultra-high-energy cosmic rays experience while they traverse intervening magnetic fields. Indeed, such deflections distort and suppress the original anisotropy in the cosmic rays arrival directions; thus, a measure of the anisotropy is also a measurement of the deflections, which in turn informs us on the chemical composition. In this work, we show that, by quantifying ultra-high-energy cosmic ray anisotropies through the angular, harmonic cross-correlation between cosmic rays and galaxies, we are able to exclude iron fractions $f_{\rm Fe}\leq{\cal O}(10\%)$ on a fiducial hydrogen map at $2\sigma$ level, and even smaller fractions in the reverse case of hydrogen on an iron map, going below $f_{\rm H}\lesssim10\%$ when we mask the Galactic Centre up to latitudes of $40\,\text{deg}$. This is an improvement of a factor of a few compared to our previous method, and is mostly ascribable to a new test statistics which is sensitive to each harmonic multipole individually. Our method can be applied to real data as an independent test of the recent claim that current cosmic-ray data can not be reproduced by any existing model of the Galactic magnetic field, as well as an additional handle to compare any realistic, competing, data-driven composition models.

We present the initial sample of redshifts for 3,839 galaxies in the MeerKAT DEEP2 field -- the deepest $\sim$1.4\,GHz radio field yet observed. Using a spectrophotometric technique combining coarse optical spectra with broadband photometry, we obtain redshifts with $\sigma_z \leq 0.01(1+z)$. The resulting radio luminosity functions between $0.2<z<1.3$ from our sample of 3,839 individual galaxies are in remarkable agreement with those inferred from modeling radio source counts, confirming an excess in radio-based SFRD$(z$) measurements at late times compared to those from the UV--IR. Several sources of systematic error are discussed -- with most having the potential of exacerbating the discrepancy -- with the conclusion that significant work remains to have confidence in a full accounting of the star formation budget of the universe.

Most galaxies, including the Milky Way, harbor a central supermassive black hole (SMBH) weighing millions to billions of solar masses. Surrounding these SMBHs are dense regions of stars and stellar remnants, such as neutron stars and black holes. Neutron stars and possibly black holes receive large natal kicks at birth on the order of hundreds of km s$^{-1}$. The natal kicks that occur in the vicinity of an SMBH may redistribute the orbital configuration of the compact objects and alter their underlying density distribution. We model the effects of natal kicks on a Galactic Center (GC) population of massive stars and stellar binaries with different initial density distributions. Using observational constraints from stellar orbits near the GC, we place an upper limit on the steepness of the initial stellar profile and find it to be core-like. In addition, we predict that $30-70 \%$ of compact objects become unbound from the SMBH due to their kicks and will migrate throughout the galaxy. Different black hole kick prescriptions lead to distinct spatial and kinematic distributions. We suggest that the Roman Space Telescope may be able to distinguish between these distributions and thus be able to differentiate natal kick mechanisms.

We present an archival Atacama Large Millimeter/submillimeter Array (ALMA) study of the CN N = 1 - 0 / CO J = 1 - 0 intensity ratio in nearby (z < 0.05) Ultra Luminous and Luminous Infrared Galaxies (U/LIRGs). We identify sixteen U/LIRGs that have been observed in both CN and CO lines at $\sim$ 500 pc resolution based on sixteen different ALMA projects. We measure the (CN bright)/CO and (CN bright)/(CN faint) intensity ratios at an ensemble of molecular clouds scales (CN bright = CN N = 1 - 0, J = 3/2 - 1/2; CN faint = CN N = 1 - 0, J = 1/2 - 1/2 hyperfine groupings). Our global measured (CN bright)/CO ratios range from 0.02-0.15 in LIRGs and 0.08-0.17 in ULIRGs. We attribute the larger spread in LIRGs to the variety of galaxy environments included in our sample. Overall, we find that the (CN bright)/CO ratio is higher in nuclear regions, where the physical and excitation conditions favour increased CN emission relative to the disk regions. 10 out of 11 galaxies which contain well-documented active galactic nuclei show higher ratios in the nucleus compared to the disk. Finally, we measure the median resolved (CN bright)/(CN faint) ratio and use it to estimate the total integrated CN line optical depth in ULIRGs ($\tau \sim$ 0.96) and LIRGs ($\tau \sim$ 0.23). The optical depth difference is likely due to the higher molecular gas surface densities found in the more compact ULIRG systems.

The presence of gaps or regions of small numbers of stars in the main sequence of the Hertzsprung Russell Diagram (HRD) of star clusters has been reported in literature. This is interesting and significant as it could be related to star formation and/or rapid evolution or instabilities. In this paper, using Gaia DR3 photometry and confirmed membership data, we explore the HRD of nine open clusters with reported gaps, identify them and assess their importance and spectral types.

NASA's \textit{Kepler} primary mission observed about 116 $deg^2$ in the sky for 3.5 consecutive years to discover Earth-like exoplanets. This mission recorded pixel cutouts, known as Target Pixel Files (TPFs), of over $200,000$ targets selected to maximize the scientific yield. The Kepler pipeline performed aperture photometry for these primary targets to create light curves. However, hundreds of thousands of background sources were recorded in the TPFs and have never been systematically analyzed. This work uses the Linearized Field Deblending (LFD) method, a Point Spread Function (PSF) photometry algorithm, to extract light curves. We use Gaia DR3 as input catalog to extract $606,900$ light curves from long-cadence TPFs. $406,548$ are new light curves of background sources, while the rest are Kepler's targets. These light curves have comparable quality as those computed by the Kepler pipeline, with CDPP values $<100$ ppm for sources $G<16$. The light curve files are available as high-level science products at MAST. Files include PSF and aperture photometry, and extraction metrics. Additionally, we improve the background and PSF modeling in the LFD method. The LFD method is implemented in the \texttt{Python} library \texttt{psfmachine}. We demonstrate the advantages of this new dataset with two examples; deblending of contaminated false positive Kepler Object of Interest identifying the origin of the transit signal; and the changes in estimated transit depth of planets using PSF photometry which improves dilution when compared to aperture photometry. This new nearly unbiased catalog enables further studies in planet search, occurrence rates, and other time-domain studies.

The atmospheric dynamics of tidally-locked hot Jupiters is dominated by the equatorial winds. Understanding the interaction between global circulation and chemistry is crucial in atmospheric studies and interpreting observations. Two-dimensional (2D) photochemical transport models shed light on how the atmospheric composition depends on circulation. In this paper, we introduce the 2D photochemical transport model, VULCAN 2D, which improves on the pseudo-2D approaches by allowing for non-uniform zonal winds. We extensively validate our VULCAN 2D with analytical solutions and benchmark comparisons. Applications to HD 189733 b and HD 209458 b reveal distinct characteristics in horizontal transport-dominated and vertical mixing-dominated regimes. Motivated by the inferred carbon-rich atmosphere by Giacobbe et al. (2021), we find that HD 209458 b with super-solar carbon-to-oxygen ratio (C/O) exhibits pronounced C2H4 absorption on the morning limb but not on the evening limb, owing to horizontal transport from the nightside. We discuss when a pseudo-2D approach is a valid assumption and its inherent limitations. Finally, we demonstrate the effect of horizontal transport in transmission observations and its impact on the morning-evening limb asymmetry with synthetic spectra, highlighting the need to consider global transport when interpreting exoplanet atmospheres.

To understand giant planet formation, we need to focus on host stars close to $1.7\ \rm M_{\odot}$, where the occurrence rate of these planets is the highest. In this initial study, we carry out pebble-driven core accretion planet formation modelling to investigate the trends and optimal conditions for the formation of giant planets around host stars in the range of $1{-}2.4\ \rm M_{\odot}$. We find that giant planets are more likely to form in systems with a larger initial disk radius; higher disk gas accretion rate; pebbles of $\sim$ millimeter in size; and birth location of the embryo at a moderate radial distance of $\sim 10$ AU. We also conduct a population synthesis study of our model and find that the frequency of giant planets and super-Earths decreases with increasing stellar mass. This contrasts the observational peak at $1.7\ \rm M_{\odot}$, stressing the need for strong assumptions on stellar mass dependencies in this range. Investigating the combined effect of stellar mass dependent disk masses, sizes, and lifetimes in the context of planet population synthesis studies is a promising avenue to alleviate this discrepancy. The hot-Jupiter occurrence rate in our models is $\sim 0.7{-}0.8\%$ around $1\ \rm M_{\odot}$ - similar to RV observations around Sun-like stars, but drastically decreases for higher mass stars.

All kinds of simulations of the intergalactic medium, such as hydrodynamic simulation, N-body simulation, numerical and semi-numerical simulation, etc., have been used to realize the history of this medium. One of these simulations is 21SSD, which is specifically focused on the epoch of reionization. This simulation deepens our understanding of the physics behind the intergalactic medium by considering the free parameters related to the Wouthuysen-Field coupling fluctuations and X-ray and Lyman line transfers in the intergalactic medium, and by presenting the plots of the power spectrum, brightness temperature, etc. in different redshifts. However, due to many physical phenomena that play significant roles in this epoch, simulations of the intergalactic medium are usually extremely complex, time-consuming, and require very powerful hardware. In this work, by using the Support Vector Regression algorithm and based on the 21SSD simulation datasets, we have tried to make the machine fully understand the brightness temperature changes in terms of redshift for different astrophysical free parameters values. At first, we trained the machine with the results of the 21SSD simulation. Then, the machine was able to predict the brightness temperature in terms of redshift with very high accuracy for other interval coefficients. Although we have used this algorithm to estimate the brightness temperature, it seems that this algorithm can be easily used for other parts of cosmology and astrophysics. With its help, it is possible to save time and obtain results with extraordinary accuracy similar to complex simulations, even with normal hardware.

The primary goal of this project was to use the given data of four emission lines from two ultra-luminous infrared galaxies to calculate the molecular gas mass and dynamical mass of each galaxy. These quantities can provide valuable information about a galaxy's age, star formation properties, and molecular make-up. Ultra-luminous infrared galaxies are formed from the merger or interaction of gas-rich galaxies, and are classified by having a luminosity of greater than 1012 L?. The presence of molecular gas clouds and dust encourages formation of young, bright stars in these galaxies, but also makes it difficult to measure the total gravitational and molecular masses, since molecular hydrogen is virtually undetectable within gas clouds. Our data instead observed emissions of carbon monoxide (CO) from the two ULIRGs at transitions of J= 1 -> 0 orJ= 2 -> 1. The software Difmap was used to map the interferometer data obtained from the Plateau de Bure Interferometer (PdBI). Our data was cleaned, filtered, and subsequently used to generate UV-plots of the emissions, along with other values needed for calculations. Some of the key results include that iras 15250 is much younger as it has not used up a majority of its molecular gas while iras 17208 is older because it has used it up, the gas mass fraction can be used to estimate the amount of dark matter present in the galaxy, and that the gas content and the central surface brightness of the disk are directly correlated.

The Golden (British Columbia, Canada) meteorite fall occurred on Oct 4, 2021 at 0534 UT with the first recovered fragment (1.3 kg) landing on an occupied bed. The meteorite is an unbrecciated, low-shock (S2) ordinary chondrite of intermediate composition, typed as an L/LL5. From noble gas measurements the cosmic ray exposure age is 25 Ma while gas retention ages are all >2 Ga. Short-lived radionuclides and noble gas measurements of the pre-atmospheric size overlap with estimates from infrasound and lightcurve modelling producing a preferred pre-atmospheric mass of 70-200 kg. The orbit of Golden has a high inclination (23.5 degs) and is consistent with delivery from the inner main belt. The highest probability (60%) of an origin is from the Hungaria group. We propose that Golden may originate among the background S-type asteroids found interspersed in the Hungaria region. The current collection of 18 L and LL chondrite orbits shows a strong preference for origins in the inner main belt, suggesting multiple parent bodies may be required to explain the diversity in CRE ages and shock states.

Recently, a new wave of full modeling analyses have emerged within the Large-Scale Structure community, leading mostly to tighter constraints on the estimation of cosmological parameters, when compared with standard approaches used over the last decade by collaboration analyses of stage III experiments. However, the majority of these full-shape analyses have primarily been conducted in Fourier space, with limited emphasis on exploring the configuration space. Investigating n-point correlations in configuration space demands a higher computational cost compared to Fourier space because it typically requires an additional integration step. This can pose a limitation when using these approaches, especially when considering higher-order statistics. One avenue to mitigate the high computation time is to take advantage of neural network acceleration techniques. In this work, we present a full shape analysis of Sloan Digital Sky Survey III/BOSS in configuration space using a neural network accelerator. We show that the efficacy of the pipeline is enhanced by a time factor $10^{3}$ without sacrificing precision, making it possible to reduce the error associated with the surrogate modeling to below $10^{-2}$ percent which is compatible with the precision required for current stage IV experiments such as DESI. We find $\Omega_m=0.286\pm 0.009$, $H_0=68.8\pm 1.2$ $\mathrm{km} \mathrm{s^{-1}}\mathrm{Mpc^{-1}}$ and $A_s \times 10^9 =2.09 ^{+0.25}_{-0.29}$. Our results on public BOSS data are in good agreement with BOSS official results and compatible with other independent full modeling analyses. We explore relaxing the prior on $\omega_b$ and varying $n_s$, without significant changes in the mean values of the cosmological parameters posterior distributions, but enlarging their widths. Finally, we explore the information content of the multipoles when constraining cosmological parameters.

We obtain the median, arithmetic mean, and the weighted mean-based central estimates for the distance to M87 using all the measurements collated in De Grijs et al (2020). We then reconstruct the error distribution for the residuals of the combined measurements and also splitting them based on the tracers used. We then checked for consistency with a Gaussian distribution and other symmetric distributions such as Cauchy, Laplacian, and Students-$t$ distributions. We find that when we analyze the combined data, the weighted mean-based estimates show a poor agreement with the Gaussian distribution, indicating that there are unaccounted systematic errors in some of the measurements. Therefore, the median-based estimate for the distance to M87 would be the most robust. This median-based distance modulus to M87 is given by $31.08 \pm 0.09$ mag and $31.07 \pm 0.09$ mag, with and without considering measurements categorized as "averages", respectively. This estimate agrees with the corresponding value obtained in DeGrijs et al (2020) to within $1\sigma$.

The inter-band correlations between optical/UV and X-ray luminosities of active galactic nuclei (AGN) are important for understanding the disc-coronal connection, as well as using AGN as standard candles for cosmology. It is conventional to measure the X-ray luminosity at rest frame 2 keV and compare to the UV luminosity at the rest-frame 2500 \AA, but the wavelength-dependence was never well explored. In this work, we adopt a well-defined sample of 1169 unobscured quasars in the redshift range 0.13 - 4.51, and apply the direct-correlation method to explore how the correlation with the 2 keV luminosity changes at different optical/UV wavelengths, from 1280 - 5550 \AA\ where the spectral quality is high. We find that the luminosity at all UV continuum wavelengths correlates with the X-ray luminosity similarly to that at 2500 \AA, and that these correlations are better than at the optical wavelengths. Strong self-correlation is also found in the broadband optical/UV continuum, supporting the scenario that it is dominated by the disc emission. Correlations of various emission lines are also investigated (e.g. C IV, C III], Mg II, H$\beta$, [O III]$\lambda\lambda 4959/5007$), including the Baldwin effect and correlations involving line-widths. We find the forms of these line correlations are different, and they are also different from their underlying continua, suggesting various complexities in the line-generation process. We discuss these results in the disc-wind scenario. Our study confirms that the rest-frame 2500 \AA\ is a good wavelength to represent the optical/UV continual properties of quasars, and shows the advantages of the direct-correlation method.

This study delves into the application of Large Language Models (LLMs) for Named Entity Recognition (NER) tasks in the field of astronomy literature. To enhance the zero-shot recognition capabilities of LLMs for astronomical named entities, we propose a strategy called Prompt-NER. Prompt-NER includes five prompt elements: Task Descriptions, Entity Definitions, Task Emphasis, Task Examples, and Second Conversation. To assess the effectiveness of the Prompt-NER strategy, we utilize three representative LLMs (Claude-2, GPT-3.5, and LLaMA-2-70b) to identify telescope and celestial object named entities in astronomical literature. Our experiments are conducted based on two distinct datasets. The first dataset comprises 30 original PDF documents, which we split into paragraphs in sequential order, resulting in a second dataset consisting of 30 paragraph collections. Additionally, we incorporate 30 astronomical telegrams to diversify our experiments and assess the performance of LLMs based on Prompt-NER on concise, complete texts. Our experimental results indicate that the Prompt-NER strategy enables LLMs to effectively accomplish NER tasks in the field of astronomy, even without prior astronomical knowledge during training. We carefully analyze the experimental results, including the mechanism of different prompt elements and the influence of different features of long and short texts on their respective experimental results. This research provides experience for zero-shot NER tasks in astronomical literature and suggests future work in this area.

Observations of metallic doublet emission lines, particularly Mg II 2796, 2803, provide crucial information for understanding galaxies and their circumgalactic medium. This study explores the effects of resonant scattering on the Mg II doublet lines and the stellar continuum in spherical and cylindrical geometries. Our findings show that under certain circumstances, resonance scattering can cause an increase in the doublet flux ratio and the escaping flux of the lines beyond what are expected in optically thin spherical media. As expected, the doublet ratio is consistently lower than the intrinsic ratio when the scattering medium is spherically symmetric and dusty. However, if the scattering medium has a disk shape, such as face-on disk galaxies, and is viewed face-on, the doublet ratio is predicted to be higher than two. These results may provide a valuable insight regarding the complexity of the shape and orientation of distant, spatially-unresolved galaxies. The importance of the continuum-pumped emission lines and expanding media is discussed to understand various observational aspects, including doublet flux ratios, which can be lower than 1.5 or higher than two, as well as symmetric or asymmetric line profiles. It is also discussed that the diffuse warm neutral medium would be an essential source of Mg II emission lines.

Measurements of the dependence of stellar specific angular momentum ($j_*$) on stellar mass ($M_*$) are presented for large samples of ALFALFA galaxies spanning the stellar mass range $\sim 10^8$-$10^{11}$ M$_{\odot}$. Accurate estimates of $j_*$ are generated using measurements of $I$-band effective radius and velocity width of the HI line profile. While the full sample ($N=3~607)$ of galaxies yields a $j_*$-$M_*$ relation with power-law index $\alpha = 0.404\pm 0.03$, it is shown that various sub-samples have indices that are very similar to the best literature results, yet with comparatively lower intrinsic scatters. A galaxy's mean $I$-band surface brightness within its effective radius ($<\mu_\mathrm{eff}>$) is shown to significantly correlate with $j_*$-$M_*$ scatter. A 3D plane fit to all $N=3~607$ galaxies in $\log_{10}j_*$-$\log_{10}M_*$-$<\mu_\mathrm{eff}>$ space yields $j_*\propto M_*^{0.589\pm 0.002}<\mu_\mathrm{eff}>^{0.193\pm 0.002}$ with scatter $\sigma=0.089$ dex. $<\mu_\mathrm{eff}>$-selected sub-samples of size up to $N=1~450$ yield power-law $j_*$-$M_*$ relations mostly consistent with $\alpha=0.55\pm 0.02$ from the literature and with intrinsic scatter ranging from 0.083 to 0.129 dex. Thus, this paper presents new, highly accurate measurements of the $j_*$-$M_*$ relation that can be used to better understand the important roles played by angular momentum in the formation and evolution of galaxies.

We present a novel numerical approach aiming at computing equilibria and dynamics structures of magnetized plasmas in coronal environments. A technique based on the use of neural networks that integrates the partial differential equations of the model, and called Physics-Informed Neural Networks (PINNs), is introduced. The functionality of PINNs is explored via calculation of different magnetohydrodynamic (MHD) equilibrium configurations, and also obtention of exact two-dimensional steady-state magnetic reconnection solutions (Craig & Henton 1995). Advantages and drawbacks of PINNs compared to traditional numerical codes are discussed in order to propose future improvements. Interestingly, PINNs is a meshfree method in which the obtained solution and associated different order derivatives are quasi-instantaneously generated at any point of the spatial domain. We believe that our results can help to pave the way for future developments of time dependent MHD codes based on PINNs

We use deep JWST/NIRSpec R~1000 slit spectra of 113 galaxies at 1.7 < z < 3.5, selected from the mass-complete Blue Jay survey, to investigate the prevalence and typical properties of neutral gas outflows at cosmic noon. We detect excess Na I D absorption (beyond the stellar contribution) in 46% of massive galaxies ($\log$ M$_*$/M$_\odot >$ 10), with similar incidence rates in star-forming and quenching systems. Half of the absorption profiles are blueshifted by at least 100 km/s, providing unambiguous evidence for neutral gas outflows. Galaxies with strong Na I D absorption are distinguished by enhanced emission line ratios consistent with AGN ionization. We conservatively measure mass outflow rates of 3 - 100 $M_\odot$ yr$^{-1}$; comparable to or exceeding ionized gas outflow rates measured for galaxies at similar stellar mass and redshift. The outflows from the quenching systems (log(sSFR)[yr$^{-1}$] $\lesssim$ -10) have mass loading factors of 4 - 360, and the energy and momentum outflow rates exceed the expected injection rates from supernova explosions, suggesting that these galaxies could possibly be caught in a rapid blowout phase powered by the AGN. Our findings suggest that AGN-driven ejection of cold gas may be a dominant mechanism for fast quenching of star formation at z~2.

[O III]$\lambda\lambda$4960,5008 doublet are often the strongest narrow emission lines in starburst galaxies and quasi-stellar objects (QSOs), and thus are a promising probe to possible variation of the fine-structure constant $\alpha$ over cosmic time. Previous such studies using QSOs optical spectra were limited to $z<1$. In this work, we constructed a sample of 40 spectra of Ly$\alpha$ emitting galaxies (LAEs) and a sample of 46 spectra of QSOs at $1.09<z<3.73$ using the VLT/X-Shooter near-infrared spectra publicly available. We measured the wavelength ratios of the two components of the spin-orbit doublet and accordingly calculated $\alpha(z)$ using two methods. Analysis on all of the 86 spectra yielded $\Delta\alpha/\alpha=(-3\pm6)\times10^{-5}$ with respect to the laboratory $\alpha$ measurements, consistent with no variation over the explored time interval. If assuming a uniform variation rate, we obtained $\alpha^{-1}{\rm d}\alpha/{\rm d}t = (-3\pm6)\times10^{-15}$ yr$^{-1}$ within the last 12 Gyrs. Extensive tests indicate that $\alpha$ variation could be better constrained using starburst galaxies' spectra than using QSO spectra in future studies.

In the pursuit of understanding the large-scale structure of the Universe, the synergy between complementary cosmological surveys has proven to be a powerful tool. Using multiple tracers of the large-scale structure can significantly improve the constraints on cosmological parameters. We explore the potential of combining the Square Kilometre Array Observatory (SKAO) and the Dark Energy Spectroscopic Instrument (DESI) spectroscopic surveys to enhance precision on the growth rate of cosmic structures. We employ a multitracer Fisher analysis to estimate precision on the growth rate when using pairs of mock surveys that are based on SKAO and DESI specifications. The pairs are at both low and high redshifts. For SKA-MID, we use the HI galaxy and the HI intensity mapping samples. In order to avoid the complexities and uncertainties at small scales, we confine the analysis to scales where linear perturbations are reliable. The consequent loss of signal in each individual survey is mitigated by the gains from the multi-tracer. After marginalising over cosmological and nuisance parameters, we find a significant improvement in the precision on the growth rate.

Context. Traditionally, supersonic turbulence is considered to be one of the most likely mechanisms to slow down the gravitational collapse in dense clumps, thereby enabling the formation of massive stars. However, several recent studies have raised differing points of view based on observations carried out with sufficiently high spatial and spectral resolution. These studies call for a re-evaluation of the role turbulence plays in massive star-forming regions. Aims. Our aim is to study the gas properties, especially the turbulence, in a sample of massive star-forming regions with sufficient spatial and spectral resolution, which can both resolve the core fragmentation and the thermal line width. Methods. We observed NH3 metastable lines with the Very Large Array (VLA) to assess the intrinsic turbulence. Results. Analysis of the turbulence distribution histogram for 32 identified NH3 cores reveals the presence of three distinct components. Furthermore, our results suggest that (1) sub- and transonic turbulence is a prevalent (21 of 32) feature of massive star-forming regions and those cold regions are at early evolutionary stage. This investigation indicates that turbulence alone is insufficient to provide the necessary internal pressure required for massive star formation, necessitating further exploration of alternative candidates; and (2) studies of seven multi-core systems indicate that the cores within each system mainly share similar gas properties and masses. However, two of the systems are characterized by the presence of exceptionally cold and dense cores that are situated at the spatial center of each system. Our findings support the hub-filament model as an explanation for this observed distribution

Radio observations of the neutral hydrogen signal from the Cosmic Dawn and Epoch of Reionisation have helped to provide constraints on the properties of the first stars and galaxies. Since this global 21-cm cosmological signal from the Cosmic Dawn is effectively constant on observing timescales and since effects resulting from systematics will vary with time, the effects of these systematics can be mitigated without the need for a model of the systematic. We present a method to account for unmodelled time-varying systematics in 21-cm radio cosmology experiments using a squared-exponential Gaussian process kernel to account for correlations between time bins in a fully Bayesian way. We find by varying the model parameters of a simulated systematic that the Gaussian process method improves our ability to recover the signal parameters by widening the posterior in the presence of a systematic and reducing the bias in the mean fit parameters. When varying the amplitude of a model sinusoidal systematic between 0.25 and 2.00 times the 21-cm signal amplitude and the period between 0.5 and 4.0 times the signal width, we find on average a 5% improvement in the root mean squared error of the fitted signal. We can use the fitted Gaussian process hyperparameters to identify the presence of a systematic in the data, demonstrating the method's utility as a diagnostic tool. Furthermore, we can use Gaussian process regression to calculate a mean fit to the residuals over time, providing a basis for producing a model of the time-varying systematic.

We use observations from the Catalina Sky Survey (CSS) to determine the bias-corrected population of small members in four very young families down to sizes equivalent to several hundred meters. Using the most recent catalog of known asteroids, we identified members from four young families for which the population has grown appreciably over recent times. A large fraction of these bodies have also been detected by CSS. We used synthetic populations of asteroids, with their magnitude distribution controlled by a small number of parameters, as a template for the bias-corrected model of these families. Applying the known detection probability of the CSS observations, we could adjust these model parameters to match the observed (biased) populations in the young families. In the case of three families, Datura, Adelaide, and Rampo, we find evidence that the magnitude distribution transitions from steep to shallow slopes near $300$ to $400$ meters. Conversely, the Hobson family population may be represented by a single power-law model. The Lucascavin family has a limited population; no new members have been discovered over the past two decades. We consider a model of parent body rotational fission with the escaping secondary tidally split into two components (thereby providing three members within this family). In support of this idea, we find that no other asteroid with absolute magnitude $H\leq 18.3$ accompanies the known three members in the Lucascavin family. A similar result is found for the archetypal asteroid pair Rheinland--Kurpfalz.

We test whether we can predict optical spectra from deep-field photometry of distant galaxies. Our goal is to perform a comparison in data space, highlighting the differences between predicted and observed spectra. The Large Early Galaxy Astrophysics Census (LEGA-C) provides high-quality optical spectra of thousands of galaxies at redshift $0.6<z<1$. Broad-band photometry of the same galaxies, drawn from the recent COSMOS2020 catalog, is used to predict the optical spectra with the spectral energy distribution (SED) fitting code Prospector and the MILES stellar library. The observed and predicted spectra are compared in terms of two age and metallicity-sensitive absorption features (H$\delta_\mathrm{A}$ and Fe4383). The global bimodality of star-forming and quiescent galaxies in photometric space is recovered with the model spectra. But the presence of a systematic offset in the Fe4383 line strength and the weak correlation between the observed and modeled line strength imply that accurate age or metallicity determinations cannot be inferred from photometry alone. For now we caution that photometry-based estimates of stellar population properties are determined mostly by the modeling approach and not the physical properties of galaxies, even when using the highest-quality photometric datasets and state-of-the-art fitting techniques. When exploring a new physical parameter space (i.e. redshift or galaxy mass) high-quality spectroscopy is always needed to inform the analysis of photometry.

Observations indicate that the accretion process in star formation may occur through accretion outbursts. This phenomenon has also now been detected in a few young massive (proto)stars (>8 Msun). The recent outburst at radio wavelengths of the massive (proto)star S255 NIRS3 has been interpreted by us as expansion of a thermal jet, fed by the infalling material. To follow up on our previous study and confirm our interpretation, we monitored the source for more than 1 yr in six bands from 1.5 GHz to 45.5 GHz and, after ~1.5 yr, with the Atacama Large Millimeter/submillimeter Array at two epochs, which made it possible to detect the proper motions of the jet lobes. The prediction of our previous study is confirmed by the new results. The radio jet is found to expand, while the flux, after an initial exponential increase, appears to stabilise and eventually decline. The radio flux measured during our monitoring is attributed to a single NE lobe, However, from 2019 a second lobe has been emerging to the SW, probably powered by the same accretion outburst, although with a delay of at least a couple of years. Flux densities at >6 GHz were satisfactorily fitted with a jet model, whereas those below 6 GHz are clearly underestimated by the model. This indicates that non-thermal emission becomes dominant at long wavelengths. Our results suggest that thermal jets can be a direct consequence of accretion events, when yearly flux variations are detected. The end of the accretion outburst is mirrored in the radio jet, as ~1 yr after the onset of the radio outburst, the inner radius of the jet began to increase while the jet mass stopped growing, as expected if the powering mechanism of the jet is quenched. Our findings support a tight connection between accretion and ejection in massive stars, consistent with a formation process involving a disk-jet system similar to that of low-mass stars.

The search for Galactic PeVatrons - astrophysical accelerators of cosmic rays to PeV energies - has entered a new phase in recent years with the discovery of the first Ultra-High-Energy (UHE, $E>100$ TeV) gamma-ray sources by the HAWC and LHAASO experiments. Establishing whether the emission is leptonic or hadronic in nature, however, requires multiwavelength data and modelling studies. Among the currently known UHE sources, LHAASO J2108+5157 is an enigmatic source without clear association to a plausible accelerator, yet spatially coincident with molecular clouds. We investigate the scenario of a molecular cloud illuminated by cosmic rays accelerated in a nearby supernova remnant (SNR) as an explanation for LHAASO J2108+5157. We aim to constrain the required properties of the SNR as well as which of the clouds identified in the vicinity is the most likely association. We use a model for cosmic ray acceleration in SNRs, their transport through the interstellar medium and subsequent interaction with molecular material, to predict the corresponding gamma-ray emission. The parameter space of SNR properties is explored to find the most plausible parameter combination that can account for the gamma-ray spectrum of LHAASO J2108+5157. In the case that a SNR is illuminating the cloud, we find that it must be young ($<10$ kyr) and located within $40-60$ pc of the cloud. A SN scenario with a low Sedov time is preferred, with a maximum proton energy of 3 PeV assumed. No SNRs matching these properties are currently known, although an as yet undetected SNR remains feasible. The galactic CR sea is insufficient to solely account for the observed flux, such that a PeVatron accelerator must be present in the vicinity.

Recently, brown dwarfs have emerged as a new topic for the astrophysical studies. These objects are intermediate between solar-type stars and giant gaseous planets. In this article, the analogies between brown dwarfs and the planet Jupiter are considered with a focus on the surrounding plasma. I consider the magnetohydrodynamic version of the Rayleigh-Taylor instability (or so called ``interchange instability'') as a minimal model of the expansion of the plasma disc surrounding Jupiter. By comparing the theoretical prediction for the radial expansion rate of the disc with the observations I quantitatively confirm the existing qualitative result, which predicts that the Rayleigh-Taylor instability provides too quick expansion. Therefore, in the realistic plasma disc yet another mechanism must operate which slows down the expansion. I suggest that similar mechanisms take place in the observed radiation belts of brown dwarfs.

Galaxy clusters and surrounding medium, can be studied using X-ray bremsstrahlung emission and Sunyaev Zel'dovich (SZ) effect. Both astrophysical probes, sample the same environment with different parameters dependance. The SZ effect is relatively more sensitive in low density environments and thus is useful to study the filamentary structures of the cosmic web. In addition, observations of the matter distribution require high angular resolution in order to be able to map the matter distribution within and around galaxy clusters. MISTRAL is a camera working at 90GHz which, once coupled to the Sardinia Radio Telescope, can reach $12''$ angular resolution over $4'$ field of view (f.o.v.). The forecasted sensitivity is $NEFD \simeq 10-15mJy \sqrt{s}$ and the mapping speed is $MS= 380'^{2}/mJy^{2}/h$. MISTRAL was recently installed at the focus of the SRT and soon will take its first photons.

We use synthetic model spectra to investigate the potential of near-ultraviolet (3000-4050 \r{A}) observations of massive O-type stars. We highlight the He I $\lambda$3188 and He II $\lambda$3203 pair as a potential temperature diagnostic in this range, supported by estimates of gravity using the high Balmer series lines. The near-ultraviolet also contains important metallic lines for determinations of chemical abundances (oxygen in particular) and estimates of projected rotational velocities for O-type spectra. Using the model spectra we present performance estimates for observations of extragalactic massive stars with the Cassegrain U-Band Efficient Spectrograph (CUBES) now in construction for the Very Large Telescope. The high efficiency of CUBES will open-up exciting new possibilities in the study of massive stars in external galaxies. For instance, CUBES will provide new insights into the physical properties of O-type stars, including oxygen abundances, in metal-poor irregular galaxies at ~1 Mpc from integrations of just 2-3 hrs. Moreover, CUBES will bring quantitative spectroscopy of more distant targets within reach for the first time, such as the O-type star (V~21.5 mag) in Leo P (at 1.6 Mpc) in only half a night of observations.

As demonstrated by Planck, SPT, and ACT, the abundance of Sunyaev-Zeldovich-detected galaxy clusters across mass and redshift is a powerful cosmological probe. Upcoming experiments such as the Simons Observatory (SO) will detect over an order of magnitude more objects than what previous experiments have found, thereby providing an unprecedented constraining potential. However, in order for this potential to be realised, the cluster detection and analysis pipelines will have to be built and understood to a much higher level of accuracy than has been demonstrated to date. Here we discuss ongoing efforts towards the accurate modelling of tSZ cluster counts, focusing on several improvements regarding optimisation bias, covariance estimation, and foreground deprojection, which are implemented in the publicly-available SZiFi package. Next, we briefly discuss the application of these improved cluster detection methods to Planck data. Finally, we introduce cosmocnc, a new cluster number count likelihood code that will be publicly available soon.

The Centaur (60558) Echeclus was discovered on March 03, 2000, orbiting between the orbits of Jupiter and Uranus. After exhibiting frequent outbursts, it also received a comet designation, 174P. If the ejected material can be a source of debris to form additional structures, studying the surroundings of an active body like Echeclus can provide clues about the formation scenarios of rings, jets, or dusty shells around small bodies. Stellar occultation is a handy technique for this kind of investigation, as it can, from Earth-based observations, detect small structures with low opacity around these objects. Stellar occultation by Echeclus was predicted and observed in 2019, 2020, and 2021. We obtain upper detection limits of rings with widths larger than 0.5 km and optical depth of $\tau$ = 0.02. These values are smaller than those of Chariklo's main ring; in other words, a Chariklo-like ring would have been detected. The occultation observed in 2020 provided two positive chords used to derive the triaxial dimensions of Echeclus based on a 3D model and pole orientation available in the literature. We obtained $a = 37.0\pm0.6$ km, $b = 28.4 \pm 0.5$ km, and $c= 24.9 \pm 0.4$ km, resulting in an area-equivalent radius of $30.0 \pm 0.5$ km. Using the projected limb at the occultation epoch and the available absolute magnitude ($\rm{H}_{\rm{v}} = 9.971 \pm 0.031$), we calculate an albedo of $p_{\rm{v}} = 0.050 \pm 0.003$. Constraints on the object's density and internal friction are also proposed.

The Gaia mission has provided us full astrometric solutions for over $1.5$B sources. However, only the brightest 34M of those have radial velocity measurements. As a proof of concept, this paper aims to close that gap, by obtaining radial velocity estimates from the low-resolution BP/RP spectra that Gaia now provides. These spectra are currently published for about 220M sources, with this number increasing to the full $\sim 2$B Gaia sources with Gaia Data Release 4. To obtain the radial velocity measurements, we fit Gaia BP/RP spectra with models based on a grid of synthetic spectra, with which we obtain the posterior probability on the radial velocity for each object. Our measured velocities show systematic biases that depend mainly on colours and magnitudes of stars. We correct for these effects by using external catalogues of radial velocity measurements. We present in this work a catalogue of about $6.4$M sources with our most reliable radial velocity measurements and uncertainties $<300$ km s$^{-1}$ obtained from the BP/RP spectra. About 23% of these have no previous radial velocity measurement in Gaia RVS. Furthermore, we provide an extended catalogue containing all 125M sources for which we were able to obtain radial velocity measurements. The latter catalogue, however, also contains a fraction of measurements for which the reported radial velocities and uncertainties are inaccurate. Although typical uncertainties in the catalogue are significantly higher compared to those obtained with precision spectroscopy instruments, the number of potential sources for which this method can be applied is orders of magnitude higher than any previous radial velocity catalogue. Further development of the analysis could therefore prove extremely valuable in our understanding of Galactic dynamics.

The luminous narrow line Seyfert galaxy PG1211+143 was the first non-BAL AGN to reveal a powerful ionized wind, based on early observations with ESA's XMM X-ray Observatory. Subsequent observations, mainly with XMM and the Japanese SUZAKU Observatory, found such winds to be a common feature of luminous AGN. Typical outflow velocities of v $\sim 0.1$c and flow momenta $m v \sim L_{\rm Edd} /c$ are consistent with winds being launched by continuum driving from a disc when the local mass accretion rate is super-Eddington. Here we report the launch of a new, ultra-fast outflow component in PG1211+143, near the end of a 5-week XMM observing campaign, and discuss its origin in an ultra-fast {\it inflow} of similar velocity detected some 3 weeks earlier. We find that the inflow lasted for at least 3 days and delivered some 10 Earth mass of fresh material into the innermost region of the source. While this mass by itself is insufficient to cause a complete inner disc restructuring, we show that it is sufficient to disrupt the X-ray emitting corona of the disc. We conclude that it is this coronal re-arrangement of the inner tens gravitational radii in PG1211+143 that subsequently caused the launch of a new wind.

We study the differences between models computed with Ledoux and Schwarzschild criteria on the internal structure, evolutionary track in the Hertzsprung-Russell diagram (HRD), lifetimes, evolution of the surface abundances and velocities, and masses of the He and CO cores. We investigate the consequences on the nature of the supernova (SN) progenitors and the type of SN events, as well as on the yields of light elements. We also study the impact on the outputs of population synthesis models. Models with initial masses between 7 and 120 M$_\odot$ at solar metallicity ($Z$=0.014) and with an initial rotation equal to 0 or 0.4 times the critical velocity at the zero-age main sequence were computed with either the Schwarzschild or the Ledoux criterion until the end of the C-burning phase. Models with initial masses between 15 and 32 M$_\odot$ computed with the Schwarzschild criterion show larger intermediate convective zones attached to the H-burning shell than models computed with the Ledoux criterion. Their CO cores and outer convective zones in the red supergiant (RSG) phase are also smaller. This impacts many outputs of stars during the core He-burning phase. Schwarzschild models have smaller CO cores and outer convective zones in the RSG phase, and their blue-to-red supergiant ratio is much higher than for Ledoux models. They also produce longer crossings of the Hertzsprung gap and favour blue loops. The upper luminosity of RSGs is little affected by the change in the convective criterion. The maximum luminosity of RSG progenitors for type II-P SN events is lowered from 5.2 to 4.95 when the Ledoux criterion is used instead of the Schwarzschild criterion in non-rotating models. The Schwarzschild criterion predicts longer-lasting, less nitrogen-enriched, and faster-rotating Cepheids. Rotational mixing decreases the differences between Schwarzschild and Ledoux models.

We investigate the detectability of single-event coalescing black hole binaries with total mass of $100-600 M_{\odot}$ at cosmological distances ($5 \lesssim z \lesssim 20$) with the next generation of terrestrial gravitational wave observatories, specifically Einstein Telescope and Cosmic Explorer. Our ability to observe these binaries is limited by the low-frequency performance of the detectors. Higher-order Multipoles of the gravitational wave signal are observable in these systems, and detection of such multipoles serves to both b the mass range over which black hole binaries are observable and improve the recovery of their individual masses and redshift. For high redshift systems of $\sim 200 M_{\odot}$ we will be able to confidently infer that the redshift is at least $z=12$, and for systems of $\sim 400 M_{\odot}$ we can infer a minimum redshift of at least $z=8$. We discuss the impact that these observations will have in narrowing uncertainties on the existence of the pair-instability mass-gap, and their implications on the formation of the first stellar black holes that could be seeds for the growth of supermassive black holes powering high-$z$ quasars.

Future generation of astronomical imaging spectrometers are targeting the far infrared wavelengths to close the THz astronomy gap. Similar to lens antenna coupled Microwave Kinetic Inductance Detectors (MKIDs), lens absorber coupled MKIDs are a candidate for highly sensitive large format detector arrays. However, the latter is more robust to misalignment and assembly issues at THz frequencies due to its incoherent detection mechanism while requiring a less complex fabrication process. In this work, the performance of such detectors is investigated. The fabrication and sensitivity measurement of several lens absorber coupled MKID array prototypes operating at 6.98 and 12 THz central frequencies is on-going.

Studies of the dynamics of globular clusters assume different values of bar parameters (mass, velocity, size) and analyse the results of orbit classifications over the range of the chosen values. It is also a usual thing that a spherical bulge component is converted into the bar to obtain a non-axisymmetric potential from an axisymmetric one. The choice of bar parameters and the way the bar is converted from the bulge introduce systematics into the orbit classifications that we explore in the present study. We integrate orbits of 30 bulge globular clusters residing in the inner area of the Galaxy ($R \lesssim 5$ kpc) backwards in time for three different potentials, two of which are obtained by fitting the rotation curve, and one is taken from the surrogate $N$-body model representing our Galaxy. We analyse each orbit in terms of dominant frequencies obtained from its coordinate spectra. We find that the bar pattern speed is a key factor in orbital classification. With an increase of it, frequencies deviate more and more from the "bar" frequency ratio 2:1. The mass and size of the bar play a smaller role. We also find that, in the $N$-body potential, the fraction of orbits that follow the bar is higher than in those obtained from fitting the rotation curve.

This work explores the relationships between galaxy sizes and related observable galaxy properties in a large volume cosmological hydrodynamical simulation. The objectives of this work are to both develop a better understanding of the correlations between galaxy properties and the influence of environment on galaxy physics in order to build an improved model for the galaxy sizes, building off of the {\it fundamental plane}. With an accurate intrinsic galaxy size predictor, the residuals in the observed galaxy sizes can potentially be used for multiple cosmological applications, including making measurements of galaxy velocities in spectroscopic samples, estimating the rate of cosmic expansion, and constraining the uncertainties in the photometric redshifts of galaxies. Using projection pursuit regression, the model accurately predicts intrinsic galaxy sizes and have residuals which have limited correlation with galaxy properties. The model decreases the spatial correlation of galaxy size residuals by a factor of $\sim$ 5 at small scales compared to the baseline correlation when the mean size is used as a predictor.

Halo Occupation Distribution (HOD) models help us to connect observations and theory, by assigning galaxies to dark matter haloes. In this work we study one of the components of HOD models: the probability distribution function (PDF), which is used to assign a discrete number of galaxies to a halo, given a mean number of galaxies. For satellite galaxies, the most commonly used PDF is a Poisson Distribution. PDFs with super-Poisson variances have also been studied, allowing for continuous values of variances. This has not been the case for sub-Poisson variances, for which only the Nearest Integer distribution, with a single variance, has been used in the past. In this work we propose a distribution based on the binomial one, which provides continuous sub-Poisson variances. We have generated mock galaxy catalogues from two dark-matter only simulations, UNIT and OUTERIM, with HOD models assuming different PDFs. We show that the variance of the PDF for satellite galaxies affects the one-halo term of the projected correlation function, and the Count-In-Cells (CIC) one point statistics. We fit the clustering of eBOSS Emission Line Galaxies, finding a preference for a sub-poissonian PDF, when we only vary the parameter controlling the PDF variance and the fraction of satellites. Using a mock catalogue as a reference, we have also included both the clustering and CIC to constrain the parameters of the HOD model. CIC can provide strong constraints to the PDF variance of satellite galaxies.

We carried out a uniform and systematic analysis of a sample of 112 nearby bright Seyfert 1 type AGN, the observations of which were carried out by the {\it Nuclear Spectroscopic Telescope Array (NuSTAR)} between August 2013 and May 2022. The main goal of this analysis is to investigate the nature of the X-ray corona in Seyfert 1 galaxies. From the physical model that fits the {\it NuSTAR} spectra, we could constrain the high energy cut-off ($\rm{E_{cut}}$) for 73 sources in our sample. For those 73 sources, we fitted the Comptonization model to estimate the temperature ($\rm{kT_{e}}$) of their corona. $\rm{kT_{e}}$ could be constrained in 42 sources. We investigated for possible correlations between various properties of the corona obtained from physical model fits to the observed spectra and between various coronal parameters and physical properties of the sources such as Eddington ratio and black hole mass. We found (a) a strong correlation between $\rm{E_{cut}}$ and the photon index and (b) a significant negative correlation between $\rm{kT_{e}}$ and the optical depth.

Here we present the cloud population extracted from M51, following the application of our new high-resolution dust extinction technique to the galaxy (Faustino Vieira et al. 2023). With this technique, we are able to image the gas content of the entire disc of M51 down to 5 pc (0.14"), which allows us to perform a statistical characterisation of well-resolved molecular cloud properties across different large-scale dynamical environments and with galactocentric distance. We find that cloud growth is promoted in regions in the galaxy where shear is minimised; i.e. clouds can grow into higher masses (and surface densities) inside the spiral arms and molecular ring. We do not detect any enhancement of high-mass star formation towards regions favourable to cloud growth, indicating that massive and/or dense clouds are not the sole ingredient for high-mass star formation. We find that in the spiral arms there is a significant decline of cloud surface densities with increasing galactocentric radius, whilst in the inter-arm regions they remain relatively constant. We also find that the surface density distribution for spiral arm clouds has two distinct behaviours in the inner and outer galaxy, with average cloud surface densities at larger galactocentric radii becoming similar to inter-arm clouds. We propose that the tidal interaction between M51 and its companion (NGC 5195) - which heavily affects the nature of the spiral structure - might be the main factor behind this.

Coronal mass ejections (CMEs) are primary drivers of space weather and studying their evolution in the inner heliosphere is vital to prepare for a timely response. Solar wind streams, acting as background, influence their propagation in the heliosphere and associated geomagnetic storm activity. This study introduces SWASTi-CME, a newly developed MHD-based CME model integrated into the Space Weather Adaptive SimulaTion (SWASTi) framework. It incorporates a non-magnetized elliptic cone and a magnetized flux rope CME model. To validate the model's performance with in-situ observation at L1, two Carrington rotations were chosen: one during solar maxima with multiple CMEs, and one during solar minima with a single CME. The study also presents a quantitative analysis of CME-solar wind interaction using this model. To account for ambient solar wind effects, two scenarios of different complexity in solar wind conditions were established. The results indicate that ambient conditions can significantly impact some of the CME properties in the inner heliosphere. We found that the drag force on the CME front exhibits a variable nature, resulting in asymmetric deformation of the CME leading edge. Additionally, the study reveals that the impact on the distribution of CME internal pressure primarily occurs during the initial stage, while the CME density distribution is affected throughout its propagation. Moreover, regardless of the ambient conditions, it was observed that after a certain propagation time (t), the CME volume follows a non-fractal power-law expansion ($\propto t^{3.03-3.33}$) due to the attainment of a balanced state with ambient.

Through observational tests of strong lensing galaxy clusters, we can test simulation derived structure predictions that follow from $\Lambda$ Cold Dark Matter ($\Lambda$CDM) cosmology. The shape and centroid deviations between the total matter distribution, stellar matter distributions, and hot intracluster gas distribution serve as an observational test of these theoretical structure predictions. We measure the position angles, ellipticities, and locations/centroids of the brightest cluster galaxy (BCG), intracluster light (ICL), the hot intracluster medium (ICM), and the core lensing mass for a sample of strong lensing galaxy clusters from the SDSS Giant Arcs Survey (SGAS). We utilize HST WFC3/IR imaging data to measure the shapes/centroids of the ICL and BCG distributions and use Chandra ACIS-I X-ray data to measure the shapes/centroids of ICM. Additionally, we measure the concentration parameter c and asymmetry parameter A to incorporate cluster dynamical state into our analysis. Using this multicomponent approach, we attempt to constrain the astrophysics of our strong lensing cluster sample and evaluate the different components in terms of their ability to trace out the DM halo of clusters in various dynamical states.

We test the hypothesis that (3) Juno is a parent body of the H chondrites with dynamical modeling of an asteroid-family-forming impact and comparison to current observational data. Using a dynamical model that includes the Yarkovsky force on a simulated Juno family and a simplified cosmic ray exposure age model we examine the expected distribution of Juno family members in both the main belt and near-Earth orbits over 300 Myrs and the cosmic ray exposure distribution for fragments exiting the main belt via the 3:1J, 5:2J, and 8:3J mean motion resonances. We find that the smallest modeled ($D<$10 m) family members of (3) Juno cannot be directly responsible for the observed H chondrite flux and that the breakup of larger family members creates an CRE distribution that resembles the measured H chondrite CRE distribution but is still unable to adequately explain the significant number of H chondrites with CRE ages of 6-8 Myrs. A similar model was performed for the asteroid (6) Hebe, another parent body candidate, and produced a CRE age distribution that is inconsistent with the measured H chondrite CRE ages. We also find from our dynamical models that we can expect $<$7 km-scale Juno family members in near-Earth orbits in the present day, consistent with the recent discovery of the shock-darkened H chondrite-like asteroid (52768) 1998 OR$_{2}$.

We analyze the role of the general relativity (GR) on the nodal librations of test particles located at the Habitable Zone (HZ) around a solar-mass star, which evolve under the influence of an eccentric planetary-mass perturber with a semimajor axis of 0.1 au. Based on a secular Hamiltonian up to quadrupole level, we derive analytical criteria that define the nodal libration region of a HZ particle as a function of its eccentricity $e_2$ and inclination $i_2$, and the mass $m_1$ and the eccentricity $e_1$ of the perturber. We show that a HZ particle can experience nodal librations with orbital flips or purely retrograde orbits for any $m_1$ and $e_1$ by adopting a suitable combination of $e_2$ and $i_2$. For $m_1 <$ 0.84 M$_\textrm{Jup}$, the greater the $m_1$ value, the smaller the $e_2$ value above which nodal librations are possible for a given $e_1$. For $m_1 >$ 0.84 M$_{\textrm{Jup}}$, a HZ test particle can undergo nodal librations for any $e_2$ and appropriate values of $e_1$ and $i_2$. The same correlation between $m_1$ and $e_2$ is obtained for nodal librations with orbital flips, but a mass limit for $m_1$ of 1.68 M$_{\textrm{Jup}}$ is required in this case. Moreover, the more massive the inner perturber, the greater the nodal libration region associated with orbital flips in the ($e_1$, $i_2$) plane for a given value of $e_2$. Finally, we find good agreements between the analytical criteria and results from N-body simulations for values of $m_1$ ranging from Saturn-like planets to super-Jupiters.

The recent discovery of astrophysical neutrinos from the Seyfert galaxy NGC 1068 suggests the presence of non-thermal protons within a compact "coronal" region close to the central black hole. The acceleration mechanism of these non-thermal protons remains elusive. We show that a large-scale magnetic reconnection layer, of the order of a few gravitational radii, may provide such a mechanism. In such a scenario, rough energy equipartition between magnetic fields, X-ray photons, and non-thermal protons is established in the reconnection region. Motivated by recent three-dimensional particle-in-cell simulations of relativistic reconnection, we assume that the spectrum of accelerated protons is a broken power law, with the break energy being constrained by energy conservation (i.e., the energy density of accelerated protons is at most comparable to the magnetic energy density). The proton spectrum is $dn_p/dE_p\propto E_p^{-1}$ below the break, and $dn_p/dE_p\propto E_p^{-s}$ above the break, with IceCube neutrino observations suggesting $s \simeq 3$. Protons above the break lose most of their energy within the reconnection layer via photohadronic collisions with the coronal X-rays, producing a neutrino signal in good agreement with the recent observations. Gamma-rays injected in photohadronic collisions are cascaded to lower energies, sustaining the population of electron-positron pairs that makes the corona moderately Compton thick.

The Milky Way (MW) is by far the best-studied galaxy and has been regarded as an ideal laboratory for understanding galaxy evolution. However, direct comparisons of Galactic and extra-galactic observations are marred by many challenges, including selection effects and differences in observations and methodology. In this study, we present a novel code GalCraft to address these challenges by generating mock integral-field spectrograph data cubes of the MW using simple stellar population models and a mock stellar catalog of the Galaxy derived from E-Galaxia. The data products are in the same format as external galaxies, allowing for direct comparisons. We investigate the ability of pPXF to recover kinematics and stellar population properties for an edge-on mock observation of the MW. We confirm that pPXF can distinguish kinematic and stellar population differences between thin and thick disks. However, pPXF struggles to recover star formation history, where the SFR is overestimated in the ranges between 2-4 and 12-14 Gyr compared to the expected values. This is likely due to the template age spacing, pPXF regularization algorithm, and spectral similarities in old population templates. Furthermore, we find systematic offsets in the recovered kinematics, potentially due to insufficient spectral resolution and the variation of line-of-sight velocity with [M/H] and age through a line-of-sight. With future higher resolution and multi-[$\alpha$/Fe] SSP templates, GalCraft will be useful to identify key signatures such as [$\alpha$/Fe]-[M/H] distribution at different $R$ and $|z|$ and potentially measure radial migration and kinematic heating efficiency to study detailed chemodynamical evolution of MW-like galaxies.

We present an investigation into a hitherto unexplored systematic that affects the accuracy of galaxy cluster mass estimates with weak gravitational lensing. Specifically, we study the covariance between the weak lensing signal, $\Delta\Sigma$, and the "true" cluster galaxy number count, $N_{\rm gal}$, as measured within a spherical volume that is void of projection effects. By quantifying the impact of this covariance on mass calibration, this work reveals a significant source of systematic uncertainty. Using the MDPL2 simulation with galaxies traced by the SAGE semi-analytic model, we measure the intrinsic property covariance between these observables within the 3D vicinity of the cluster, spanning a range of dynamical mass and redshift values relevant for optical cluster surveys. Our results reveal a negative covariance at small radial scales ($R \lesssim R_{\rm 200c}$) and a null covariance at large scales ($R \gtrsim R_{\rm 200c}$) across most mass and redshift bins. We also find that this covariance results in a $2-3\%$ bias in the halo mass estimates in most bins. Furthermore, by modeling $N_{\rm gal}$ and $\Delta\Sigma$ as multi-(log)-linear equations of secondary halo properties, we provide a quantitative explanation for the physical origin of the negative covariance at small scales. Specifically, we demonstrate that the $N_{\rm gal}$-$\Delta\Sigma$ covariance can be explained by the secondary properties of halos that probe their formation history. We attribute the difference between our results and the positive bias seen in other works with (mock)-cluster finders to projection effects. These findings highlight the importance of accounting for the covariance between observables in cluster mass estimation, which is crucial for obtaining accurate constraints on cosmological parameters.

The gas in the circumgalactic medium (CGM) of galaxies is being probed by absorbers in the spectrum of background objects in observations, yet the physical properties of the medium remain unconstrained. We use the cosmological hydrodynamical simulation TNG50 to statistically trace the origins of HI Ly-$\alpha$ absorbers around galaxies at $z = 0.5$ with stellar masses ranging from 10$^8$ to 10$^{11}$ M$_\odot$. By including all gas with a line of sight velocity within $\pm 500$ \kms\ of the central, we mimic observational studies of the CGM to quantitatively assess the impact of other galaxy haloes and overdense gas in the IGM that intersect sightlines. We find that 75 per cent of HI absorbers with column densities log(NHI) $> 16.0$ trace the central galaxy within $\pm150$ (80) \kms\ of $M_* = 10^{10} (10^8)$ M$_\odot$ central galaxies. The impact of satellites to the total absorber fraction is most significant at impact parameters $0.5 R_{\rm vir} < b < R_{\rm vir}$ and satellites with masses below typical detection limits ($M_* < 10^8$ M$_\odot$) account for 10 (40) per cent of absorbers that intersect any satellite bound to $10^{10}$ and $10^{11}$ $(10^9)$ M$_\odot$ centrals. After confirming outflows are more dominant along the minor axis, we additionally show that at least 20 per cent of absorbers are found without strong radial motions, highlighting that absorbers can also trace quasi-static gas. The metallicity of absorbers also depends on the azimuthal angle, but this signal is largely driven by enriched inflowing and quasi-static gas. Our work shows that determining the stellar mass of galaxies at $z_{\rm abs}$ is essential to constrain the physical origin of the gas traced in absorption, which in turn is key to characterising the kinematics and distribution of gas and metals in the CGM.

Magnetic materials are particularly favorable targets for detecting axions interacting with electrons because the collective excitation of electron spins, the magnon, can be excited through the axion-magnon conversion process. It is often assumed that only the zero-momentum uniformly precessing magnetostatic (Kittel) mode of the magnon is excited. This is justified if the de Broglie wavelength of the axion is much longer than the size of the target magnetic material. However, if the de Broglie wavelength is shorter, finite-momentum magnon modes can also be excited. We systematically analyze the target material size dependence of the axion-magnon conversion rate. We discuss the importance of these effects in the detection of relativistic axions as well as in the detection of axion dark matter of relatively heavy mass with large material size.

This work investigates static and dynamical quark star properties within a $D_3-D_7$ holographic model. We solve the Tolman-Oppenheimer-Volkoff equations for the quark matter equation of state obtained from the brane configuration and determine the range of model parameters in which the quark star family mass-radius diagram are compatible with recent NICER observational data for the pulsars PSR J$0030+0451$ and PSR J$0740+6620$. We show that the model supports stable configurations with maximum masses higher than $2$ Solar masses, in line with the inferred masses of the pulsars PSR J$1614-2230$, PSR J$0348+0432$ and PSR J$0740+6620$. Furthermore, we show that for the constituent quark mass ranging from $m=330$ MeV to $m=350$ MeV, the tidal deformability parameter for each component of the binary star system is consistent with the GW170817 event detected by the LIGO-Virgo collaboration.

Terrestrial gravity perturbations caused by seismic fields produce the so-called Newtonian noise in gravitational-wave detectors, which is predicted to limit their sensitivity in the upcoming observing runs. In the past, this noise was seen as an infrastructural limitation, i.e., something that cannot be overcome without major investments to improve a detector's infrastructure. However, it is possible to have at least an indirect estimate of this noise by using the data from a large number of seismometers deployed around a detector's suspended test masses. The noise estimate can be subtracted from the gravitational-wave data; a process called Newtonian-noise cancellation (NNC). In this article, we present the design and implementation of the first NNC system at the Virgo detector as part of its AdV+ upgrade. It uses data from 110 vertical geophones deployed inside the Virgo buildings in optimized array configurations. We use a separate tiltmeter channel to test the pipeline in a proof-of-principle. The system has been running with good performance over months.

Earlier, it was a standard assumption that the entire core of neutron stars is superconducting. However, the matter contents in the inner core has been unknown even qualitatively, because the density of matter in that region is expected to be higher than the nuclear saturation density 0.16 $\mathrm{fm}^{-3}$. As a consequence, no reliable model exists that would describe the neutron star matter in the inner core of neutron stars. Thus, a possibility of presence of normal, nonsuperconducting, plasma in the inner core cannot be excluded as of today. This point is supported by the numerical calculations performed in [1]. The calculations are based on the equation of state and the proton Cooper pairing gap energy derived from the chiral effective field theory. The numerical results show that the superconducting gap goes to zero beyond the depth about 1 km below the crust-core boundary. Given that the stellar radius is of the order of 12 km, therefore the superconducting proton matter is expected to exist only in a thin layer at the tip of the outer core. Recently it has been realized that the symmetry of superconductor is anisotropic in the lasagna region of the pasta phases located at the bottom of the crust. However the question of whether this symmetry is continuous or discreet was unsolved. The numerical calculations performed in [1] have shown that the tunneling rate between the adjacent slabs in the entire range of the corresponding densities is negligibly small. Thus, a discreet model is necessary for the description of the lasagna region. Uncertainties and future directions of the research are discussed.

We use previously obtained experimental results by neutron interferometry to effectively constrain the parameter space of several prominent dark energy models. This investigation encompasses the environment-dependent dilaton field, a compelling contender for dark energy that emerges naturally within the strong coupling limit of string theory, alongside symmetron and chameleon fields. Our study presents substantial improvements over previous constraints of the dilaton and symmetron fields, improving parameter constraints by several orders of magnitude. However, the analysis does not yield any new constraints on the chameleon field. Furthermore, we establish constraints for the projected neutron split interferometer, which has recently concluded a decisive proof-of-principle demonstration. Our symmetron simulations reveal that depending on the parameter values there are multiple static solutions with increasing number of nodes and increasing energy inside a cylindrical vacuum chamber. This agrees with results obtained earlier in the literature for infinitely parallel plates. Interestingly, while these multiple solutions can correspond to domain walls forming inside the vacuum chamber, we also find solutions that do not reach their vacuum expectation value inside the vacuum chamber, but display multiple nodes nonetheless.

We consider a nearly collisionless plasma consisting of a species of `test particles' in 1D-1V, stirred by an externally imposed stochastic electric field. The mean effect on the particle distribution function is stochastic heating. Accompanying this heating is the generation of fine-scale structure in the distribution function, which we characterize with the collisionless (Casimir) invariant $C_2 \propto \iint dx dv \, \langle f^2 \rangle$. We find that $C_2$ is transferred from large scales to small scales in both position and velocity space via a phase-space cascade enabled by both particle streaming and nonlinear interactions between particles and the stochastic electric field. We compute the steady-state fluxes and spectrum of $C_2$ in Fourier space, with $k$ and $s$ denoting spatial and velocity wavenumbers, respectively. Whereas even the linear phase mixing alone would lead to a constant flux of $C_2$ to high $s$ (towards the collisional dissipation range) at every $k$, the nonlinearity accelerates this cascade by intertwining velocity and position space so that the flux of $C_2$ is to both high $k$ and high $s$ simultaneously. Integrating over velocity (spatial) wavenumbers, the $k$-space ($s$-space) flux of $C_2$ is constant down to a dissipation length (velocity) scale that tends to zero as the collision frequency does, even though the rate of collisional dissipation remains finite. The resulting spectrum in the inertial range is a self-similar function in the $(k,s)$ plane, with power-law asymptotics at large $k$ and $s$. We argue that stochastic heating is made irreversible by this entropy cascade and that, while collisional dissipation accessed via phase mixing occurs only at small spatial scales rather than at every scale as it would in a linear system, the cascade makes phase mixing even more effective overall in the nonlinear regime than in the linear one.