Papers for Wednesday, Nov 27 2024

The Euclid Wide Survey (EWS) is predicted to find approximately 170 000 galaxy-galaxy strong lenses from its lifetime observation of 14 000 deg^2 of the sky. Detecting this many lenses by visual inspection with professional astronomers and citizen scientists alone is infeasible. Machine learning algorithms, particularly convolutional neural networks (CNNs), have been used as an automated method of detecting strong lenses, and have proven fruitful in finding galaxy-galaxy strong lens candidates. We identify the major challenge to be the automatic detection of galaxy-galaxy strong lenses while simultaneously maintaining a low false positive rate. One aim of this research is to have a quantified starting point on the achieved purity and completeness with our current version of CNN-based detection pipelines for the VIS images of EWS. We select all sources with VIS IE < 23 mag from the Euclid Early Release Observation imaging of the Perseus field. We apply a range of CNN architectures to detect strong lenses in these cutouts. All our networks perform extremely well on simulated data sets and their respective validation sets. However, when applied to real Euclid imaging, the highest lens purity is just 11%. Among all our networks, the false positives are typically identifiable by human volunteers as, for example, spiral galaxies, multiple sources, and artefacts, implying that improvements are still possible, perhaps via a second, more interpretable lens selection filtering stage. There is currently no alternative to human classification of CNN-selected lens candidates. Given the expected 10^5 lensing systems in Euclid, this implies 10^6 objects for human classification, which while very large is not in principle intractable and not without precedent.

Sub-Neptunes - volatile-rich exoplanets smaller than Neptune - are intrinsically the most common type of planet known. However, the formation and nature of these objects, as well as the distinctions between sub-classes (if any), remain unclear. Two powerful tools to tease out the secrets of these worlds are measurements of (i) atmospheric composition and structure revealed by transit and/or eclipse spectroscopy, and (ii) mass, radius, and density revealed by transit photometry and Doppler spectroscopy. Here we present OrCAS, a survey to better elucidate the origins, compositions, and atmospheres of sub-Neptunes. This radial velocity survey uses a repeatable, quantifiable metric to select targets suitable for subsequent transmission spectroscopy and address key science themes about the atmospheric & internal compositions and architectures of these systems. Our survey targets 26 systems with transiting sub-Neptune planet candidates, with the overarching goal of increasing the sample of such planets suitable for subsequent atmospheric characterization. This paper lays out our survey's science goals, defines our target prioritization metric, and performs light-curve fits and statistical validation using existing TESS photometry and ground-based follow-up observations. Our survey serves to continue expanding the sample of small exoplanets with well-measured properties orbiting nearby bright stars, ensuring fruitful studies of these systems for many years to come.

We present the discovery of a large extended radio jet associated with the extremely radio-loud quasar J1601+3102 at $z\sim5$ from sub-arcsecond resolution imaging at 144 MHz with the LOFAR International Telescope. These large radio lobes have been argued to remain elusive at $z>4$ due to energy losses in the synchrotron emitting plasma as a result of scattering of the strong CMB at these high redshifts. Nonetheless, the 0.3" resolution radio image of J1601+3102 reveals a Northern and Southern radio lobe located at 9 and 57 kpc from the optical quasar, respectively. The measured jet size of 66 kpc makes J1601+3102 the largest extended radio jet at $z>4$ to date. However, it is expected to have an even larger physical size in reality due to projection effects brought about by the viewing angle. Furthermore, we observe the rest-frame UV spectrum of J1601+3102 with Gemini/GNIRS to examine its black hole properties, which results in a mass of 4.5$\times$10$^{8}$ M$_{\odot}$ with an Eddington luminosity ratio of 0.45. The BH mass is relatively low compared to the known high-$z$ quasar population, which suggests that a high BH mass is not strictly necessary to generate a powerful jet. This discovery of the first $\sim100$ kpc radio jet at $z>4$ shows that these objects exist despite energy losses from Inverse Compton scattering and can put invaluable constraints on the formation of the first radio-loud sources in the early Universe.

For aspherical interstellar dust grains aligned with their short axes preferentially parallel to the local magnetic field, the amount of extinction per grain is larger when the magnetic field is along the line of sight and smaller when in the plane of the sky. To the extent that optical extinction arises from both aligned and unaligned grain populations with different extinction properties, changes in the magnetic field orientation induces changes in its wavelength dependence, parameterized by $R_V \equiv A_V/E(B-V)$. We demonstrate that the measured total and polarized extinction curves of the diffuse Galactic interstellar medium imply $R_V$ varies from 3.21 when the magnetic field is along the line of sight ($\psi = 0$) to $R_V = 3.05$ when in the plane of the sky ($\psi = 90^\circ$). This effect could therefore account for much of the large-scale $R_V$ variation observed across the sky ($\sigma(R_V) \simeq 0.2$), particularly at high Galactic latitudes.

Magnetar Giant Flares (MGFs) are the most energetic non-catastrophic transients known to originate from stellar objects. The first discovered events were nearby. In recent years, several extragalactic events have been identified, implying an extremely high volumetric rate. We show that future instruments with a sensitivity $\lesssim 5\times 10^{-9}$ erg cm$^{-2}$ at $\sim 1$ MeV will be dominated by extragalactic MGFs over short gamma-ray bursts (sGRBs). Clear discrimination of MGFs requires intrinsic GRB localization capability to identify host galaxies. As MGFs involve a release of a sizable fraction of the neutron star's magnetic free energy reservoir in a single event, they provide us with invaluable tools for better understanding magnetar birth properties and the evolution of their magnetic fields. A major obstacle is to identify a (currently) small sub-population of MGFs in a larger sample of more energetic and distant sGRBs. We develop the tools to analyze the properties of detected events and their occurrence rate relative to sGRBs. Even with the current (limited) number of events, we can constrain the initial internal magnetic field of a typical magnetar at formation to be $B_0\approx 3\times 10^{14}-2\times 10^{15}$\,G. Larger samples will constrain the distribution of birth fields. We also estimate the contribution of MGFs to the gravitational wave (GW) stochastic background. Depending on the acceleration time of baryon-loaded ejecta involved in MGFs, their GW emission may reach beyond 10~kHz and, if so, will likely dominate over other conventional astrophysical sources in that frequency range.

Using stellar kinematic data from Gaia DR3, we revisit constraints on black hole (BH) natal kicks from observed accreting and detached BH binaries. We compare the space velocities and Galactic orbits of a sample of 12 BHs in the Galactic disk with well-constrained distances to their local stellar populations, for which we obtain proper motions and radial velocities from Gaia DR3. Compared to most previous studies, we infer lower minimum kick velocities, because our modeling accounts for the fact that most BH binaries are old and have likely been kinematically heated by processes other than kicks. Nevertheless, we find that half of the BHs have at least weak evidence for a kick, being kinematically hotter than at least 68% of their local stellar populations. At least 4 BHs are kinematically hotter than 90% of their local stellar populations, suggesting they were born with kicks of $\gtrsim 100$ km s$^{-1}$. On the other hand, 6 BHs have kinematics typical of their local populations, disfavoring kicks of $\gtrsim 50$ km s$^{-1}$. For two BHs, V404 Cyg and VFTS 243, there is strong independent evidence for a very weak kick $\lesssim 10$ km s$^{-1}$. Our analysis implies that while some BHs must form with very weak kicks, it would be wrong to conclude that most BHs do, particularly given that selection biases favor weak kicks. Although the uncertainties on most individual BHs' kicks are still too large to assess whether the kick distribution is bimodal, the data are consistent with a scenario where some BHs form by direct collapse and receive weak kicks, and others form in supernovae and receive strong kicks.

We study the influence of the ambient large-scale cold-gas vorticity on the specific star formation rate (sSFR) of all central galaxies with stellar masses of $10.0<\log\,M_{\ast}/\mathrm{M_{\odot}}<11.5$, using the IllustrisTNG-100 simulation. The cold-gas vorticity defined and calculated for gas with $T_{\rm gas} < 2\times 10^4 \mathrm{K}$ and on scales of $\sim$ 1 Mpc can well describe the angular motion of the ambient cold gas. We find crucial evidence for a clear connection between the cold gas spin/vorticity and star formation activeness, in that at any given halo mass, galaxies that live in a higher cold-gas vorticity environment are generally less actively star forming, regardless of the large-scale environment type (filament or knot) the galaxy lives in, or it being star-forming or quenched. In particular, at any fixed halo mass scale, the environmental cold-gas vorticities of galaxies in filaments are generally higher than those of galaxies in knots, naturally explaining lower the sSFRs of filament galaxies than of knot galaxies. This large-scale cold-gas vorticity is also highly connected to the orbital angular momentum of environmental galaxies up to a distance of $\sim$ 500 kpcs, indicating their common origin and/or possible angular momentum inheritance/modulation from the latter to the former. The negative modulation by the environmental vorticity to galaxy star formation is only significantly observed for the cold gas, indicating the unique role of cold-gas angular momentum.

Galaxy clusters and groups are the last link in the chain of hierarchical structure formation. Their environments can be significantly affected by outbursts from AGN, especially in groups where the medium density is lower and the gravitational potential shallower. The interaction between AGN and group weather can therefore greatly impact their evolution. We investigate the non-thermal radio emission in Abell 1213, a galaxy group which is part of a larger sample of ~50 systems (X-GAP) recently granted XMM-Newton observations. We exploit proprietary LOFAR 54 MHz and uGMRT 380 MHz observations, complementing them with 144 MHz LOFAR survey and XMM-Newton archival data. A1213 hosts a bright AGN associated with one of the central members, 4C 29.41, which was previously optically identified as a dumb-bell galaxy. Observations at 144 MHz at a resolution of 0.3'' allow us to resolve the central radio galaxy. From this source, a ~500 kpc-long tail extends North-East. Our analysis suggests that the tail likely originated from a past outburst of 4C 29.41, and its current state might be the result of the interaction with the surrounding environment. The plateau of the spectral index distribution in the Easternmost part of the tail suggests mild particle re-acceleration, that could have re-energised seed electrons from the past activity of the AGN. While we observe a spatial and physical correlation of the extended, central emission with the thermal plasma, which might hint at a mini-halo, current evidence cannot conclusively prove this. A1213 is only the first group, among the X-GAP sample, that we are able to investigate through low-frequency radio observations. Its complex environment once again demonstrates the significant impact that the interplay between thermal and non-thermal processes can have on galaxy groups.

When the core of a massive star collapses, neutrino heating can energize the stalled accretion shock, leading to a successful supernova. The critical condition that characterizes the transition from accretion to explosion is a central topic of study and is often characterized by a critical proto-neutron star (PNS) neutrino luminosity $L_\nu^{\rm crit}$, which depends on the post-collapse mass accretion rate $\dot{M}$ from the progenitor. We examine the critical condition by solving the spherically symmetric time-dependent Euler equations with a general equation of state and realistic microphysics for a range of $\dot{M}$, average neutrino energy $\left< \epsilon_{\nu}\right>$, luminosity $L_{\nu}$, PNS radius $R_{\star}$, mass $M_{\star}$, and pre-shock Mach number $\mathcal{M}$ for a fixed neutrino optical depth from the PNS surface of $2/3$. We derive $L_{\nu}^{\mathrm{crit}}$ as a function of the input parameters. We show that pressurized pre-shock flow, as parameterized by low $\mathcal{M}$, changes the normalization of the critical condition because accretion of higher entropy shells at later times after collapse leads to lower $L_{\nu}^{\mathrm{crit}}$. We connect this finding to the onset of explosion due to compositional interface accretion. Across our parameter space, we test critical conditions that have been proposed in the literature, including the ``antesonic" condition, the ``force explosion condition," and the heuristic heating-advection timescale condition. We discuss how shock oscillations impact these critical conditions. Compared to other explosion conditions, we find that the antesonic ratio shows the least variation across the model space we explore. This work is preparatory for similar experiments in 2D axisymmetry and 3D.

We present the results of applying anomaly detection algorithms to a quasar spectroscopic sub-sample from the SDSS DR16 Quasar Catalog, covering the redshift range 1.88 < z < 2.47. Principal Component Analysis (PCA) was employed for dimensionality reduction of the quasar spectra, followed by hierarchical K-Means clustering in a 20-dimensional PCA eigenvector hyperspace. To prevent broad absorption line (BAL) quasars from being identified as the primary anomaly group, we conducted the analysis with and without them, comparing both datasets for a clearer identification of other anomalous quasar types. We identified 1,888 anomalous quasars, categorized into 10 broad groups. The anomalous groups include C IV Peakers-quasars with extremely strong and narrow C IV emission lines; Excess Si IV emitters-quasars where the Si IV line is as strong as the C IV line; and Si IV Deficient anomalies, which exhibit significantly weaker Si IV emission compared to typical quasars. The anomalous nature of these quasars is attributed to lower Eddington ratios for C IV Peakers, super-solar metallicity for Excess Si IV emitters, and sub-solar metallicity for Si IV Deficient anomalies. Additionally, we identified four groups of BAL anomalies: Blue BALs, Flat BALs, Reddened BALs, and FeLoBALs, distinguished primarily by the strength of reddening in these sources. Further, among the non-BAL quasars, we identified three types of reddened anomaly groups classified as heavily reddened, moderately reddened, and plateau-shaped spectrum quasars, each exhibiting varying degrees of reddening. The detected anomalies are presented as a value-added catalog.

The stellar disc is the dominant luminous component of the Milky Way (MW). Although our understanding of its structure is rapidly expanding due to advances in large-scale stellar surveys, our picture of the MW disc remains substantially obscured by selection functions and incomplete spatial coverage of observational data. In this work, we present the comprehensive chrono-chemo-kinematic structure of the MW disc, recovered using a novel orbit superposition approach combined with data from APOGEE DR 17. We detect periodic azimuthal metallicity variations within 6-8 kpc with an amplitude of 0.05-0.1 dex peaking along the bar major axis. The radial metallicity profile of the MW also varies with azimuth, displaying a pattern typical among other disc galaxies: a decline outside the solar radius and an almost flat profile in the inner region, attributed to the presence of old, metal-poor high-{\alpha} populations, which comprise about 40% of the total stellar mass. The geometrically defined thick disc and the high-{\alpha} populations have comparable masses, with differences in their stellar population content, which we quantify using the reconstructed 3D MW structure. The well-known [{\alpha}/Fe]-bimodality in the MW disc, once weighted by stellar mass, is less pronounced at a given metallicity for the whole galaxy but distinctly visible in a narrow range of galactic radii (5-9 kpc), explaining its relative lack of prominence in external galaxies and galaxy formation simulations. Analysing a more evident double age-abundance sequence, we construct a scenario for the MW disc formation, advocating for an inner/outer disc dichotomy genetically linked to the MW's evolutionary stages. In this picture, the extended solar vicinity is a transition zone that shares chemical properties of both the inner (old age-metallicity sequence) and outer discs (young age-metallicity sequence).

The X-ray polarization observations made possible with the Imaging X-ray Polarimetry Explorer (IXPE) offer new ways of probing high-energy emission processes in astrophysical jets from blazars. Here we report on the first X-ray polarization observation of the blazar S4 0954+65 in a high optical and X-ray state. During our multi-wavelength campaign on the source, we detected an optical flare whose peak coincided with the peak of an X-ray flare. This optical-X-ray flare most likely took place in a feature moving along the parsec-scale jet, imaged at 43 GHz by the Very Long Baseline Array. The 43 GHz polarization angle of the moving component underwent a rotation near the time of the flare. In the optical band, prior to the IXPE observation, we measured the polarization angle to be aligned with the jet axis. In contrast, during the optical flare the optical polarization angle was perpendicular to the jet axis; after the flare, it reverted to being parallel to the jet axis. Due to the smooth behavior of the optical polarization angle during the flare, we favor shocks as the main acceleration mechanism. We also infer that the ambient magnetic field lines in the jet were parallel to the jet position angle. The average degree of optical polarization during the IXPE observation was (14.3$\pm$4.1)%. Despite the flare, we only detected an upper limit of 14% (at 3$\sigma$ level) on the X-ray polarization degree; although a reasonable assumption on the X-ray polarization angle results in an upper limit of 8.8% ($3\sigma$). We model the spectral energy distribution (SED) and spectral polarization distribution (SPD) of S4 0954+65 with leptonic (synchrotron self-Compton) and hadronic (proton and pair synchrotron) models. The constraints we obtain with our combined multi-wavelength polarization observations and SED modeling tentatively disfavor hadronic models for the X-ray emission in S4 0954+65.

We present a suite of zoom-in cosmological simulations of Milky Way-like galaxies with a prominent disc component and a strong bar in their centre, based on a subsample of barred galaxies from the TNG50 magneto-hydrodynamic simulation. We modify the physical models that regulate star formation, namely, supernova feedback and black hole quasar feedback, to examine how they affect the disc and bar formation. We find that the morphology remains consistent in all the galaxies that develop a massive stellar disc ($>10^{10}\Msun$), which is dominant in comparison with the bulge mass. The black hole quasar feedback models used in this work do not affect bar formation, although they can affect the properties of the bar. The energy released by the supernovae causes a delay in the time of bar formation and, in models with strong feedback, stops bar formation. This could be understood since supernova feedback affects the assembly of the disc and bulge and their structural properties, such as the mass content, size, and radial velocity dispersion. We compared our predictions to three bar instability criteria proposed in the literature. We find that galaxies with varying supernovae and black hole quasar feedback satisfy these criteria at the moment of bar formation, except in extreme cases where the galaxy does not have or has weak supernova feedback, where for some galaxies, they fail to predict the existence (or no existence) of bars. Our findings provide insights into the physical process behind bar formation while putting constraints on the analytic prescriptions that predict bar formation.

The nature of sub-Neptunes remains unknown due to degeneracies in interior structure solutions. However, a statistical set of small planets with measured masses and radii has been compiled. It can be used to test the prediction of large water reservoirs on sub-Neptunes by planet formation theory. We want to find out whether this water reservoir is included in photoevaporative winds and how much of it can partition into the rocky and metallic interior. We couple the result of a planetary formation model to evolution models which assume perfect mixing of water with H/He in the envelope or complete segregation. For the mixed envelopes, we also include fractionation during photoevaporative mass-loss. Further, the effect of equilibrium dissolution of water into an assumed magma ocean and into the metallic core is studied for the first time in coupled formation-evolution models. Out of the tested scenarios, the mass-radius relation of exoplanets is best matched under the mixed assumption without water sequestration to the interior. We quantify the radius valley location and scaling with stellar mass. Fractionation is not found to significantly alter the composition of the planets for our initial conditions due to initially massive envelopes on all planets. In contrast, water sequestration has a profound effect on the radius evolution and compositional budget of the planets. The model predicts the preservation of large quantities of water even if the gaseous envelope is lost. Planets with corresponding bulk densities are not observed in comparably large numbers. By combining formation and evolution model, we probe a parameter space favored by core accretion theory. We conclude that the dissolution of different volatiles into the planetary interior and solidification of the magma ocean are natural next steps for comprehensive treatment of atmosphere-interior interaction. (abridged)

The Generalized Needlet Internal Linear Combination (GNILC) method is a non-parametric component separation algorithm to remove the foreground contamination of the 21-cm intensity mapping data. In this work, we perform the Discrete Cosine Transform (DCT) along the frequency axis in the expanded GNILC framework (denoted eGNILC) which helps reduce the power loss in low multipoles, and further demonstrate its performance. We also calculate the eGNILC bias to modify the criterion for determining the degrees of freedom of the foreground (dof), and embed the Robust Principal Component Analysis (RPCA) in mixing matrix computation to obtain a blind component separation method. We find that the eGNILC bias is related to the averaged domain size and the dof of the foreground but not the underlying 21-cm signal. In case of no beam effect, the eGNILC bias is negligible for simple power law foregrounds outside the Galactic plane. We also examine the eGNILC performance in the SKA-MID (SKA Phase-I in mid-frequency) and BINGO (Baryon Acoustic Oscillations from Integrated Neutral Gas Observations) simulations. We show that if the adjacent frequency channels are not highly correlated, eGNILC can recover the underlying 21-cm signal with good accuracy. With the varying Airy disk beam applied to both SKA-MID and BINGO, the power spectra of 21-cm can be effectively recovered at the multipoles $\ell \in [20, 250]$ and $[20, 300]$ respectively. With no instrumental noise, the SKA-MID exhibits $\lesssim 20\%$ power loss and BINGO exhibits $\sim 10\%$ power loss. The varying Airy-disk beam only causes significant errors at large multipoles.

The next generation of cosmological spectroscopic sky surveys will probe the distribution of matter across several Gigaparsecs (Gpc) or many billion light-years. In order to leverage the rich data in these new maps to gain a better understanding of the physics that shapes the large-scale structure of the cosmos, observed matter distributions must be compared to simulated mock skies. Small mock skies can be produced using precise, physics-driven hydrodynamical simulations. However, the need to capture small, kpc-scale density fluctuations in the intergalactic medium (IGM) places tight restrictions on the necessary minimum resolution of these simulations. Even on the most powerful supercomputers, it is impossible to run simulations of such high resolution in volumes comparable to what will be probed by future surveys, due to the vast quantity of data needed to store such a simulation in computer memory. However, it is possible to represent the essential features of these high-resolution simulations using orders of magnitude less memory. We present a hybrid approach that employs a physics-driven hydrodynamical simulation at a much lower-than-necessary resolution, followed by a data-driven, deep-learning Enhancement. This hybrid approach allows us to produce hydrodynamic mock skies that accurately capture small, kpc-scale features in the IGM but which span hundreds of Megaparsecs. We have produced such a volume which is roughly one Gigaparsec in diameter and examine its relevant large-scale statistical features, emphasizing certain properties that could not be captured by previous smaller simulations. We present this hydrodynamic volume as well as a companion n-body dark matter simulation and halo catalog which we are making publically available to the community for use in calibrating data pipelines for upcoming survey analyses.

Coronal loops generally trace magnetic lines of force in the upper solar atmosphere. Understanding the loop morphology and its temporal evolution has implications for coronal heating models that rely on plasma heating due to reconnection at current sheets. Simultaneous observations of coronal loops from multiple vantage points are best suited for this purpose. Here, we report a stereoscopic analysis of coronal loops in an active region based on observations from the Extreme Ultraviolet Imager on board the Solar Orbiter Spacecraft and the Atmospheric Imaging Assembly on board the Solar Dynamics Observatory. Our stereoscopic analysis reveals that coronal loops have nearly circular cross-sectional widths and that they exhibit temporally coherent intensity variations along their lengths on timescales of around 30\,minutes. Results suggest that coronal loops can be best represented as three-dimensional monolithic or coherent plasma bundles that outline magnetic field lines. Therefore, at least on the scales resolved by Solar Orbiter, it is unlikely that coronal loops are manifestations of emission from the randomly aligned wrinkles in two-dimensional plasma sheets along the line-of-sight as proposed in the `coronal veil' hypothesis.

Super sample covariance (SSC) is important when estimating covariance matrices using a set of mock catalogues for galaxy surveys. If the underlying cosmological simulations do not include the variation in background parameters appropriate for the simulation sizes, then the scatter between mocks will be missing the SSC component. The coupling between large and small modes due to non-linear structure growth makes this pernicious on small scales. We compare different methods for generating ensembles of mocks with SSC built in to the covariance, and contrast against methods where the SSC component is computed and added to the covariance separately. We find that several perturbative expansions, developed to derive background fluctuations, give similar results. We then consider scaling covariance matrices calculated for simulations of different volumes to improve the accuracy of covariance matrix estimation for a given amount of computational time. On large scales, we find that the primary limitation is from the discrete number of modes contributing to the measured power spectrum, and we propose a new method for correcting this effect. Correct implementation of SSC and the effect of discrete mode numbers allows covariance matrices created from mocks to be scaled between volumes, potentially leading to a significant saving on computational resources when producing covariance matrices. We argue that a sub-percent match is difficult to achieve because of the effects of modes on scales between the box sizes, which cannot be easily included. Even so, a 3% match is achievable on scales of interest for current surveys scaling the simulation volume by 512x, costing a small fraction of the computational time of running full-sized simulations. This is comparable to the agreement between analytic and mock-based covariance estimates to be used with DESI Y1 results.

Published analyses of very long baseline interferometry (VLBI) data for the sources included in the third International Celestial Reference Frame (ICRF3) catalog have revealed object-specific, excess astrometric variability and quasi-coherent trajectories as functions of time. A fraction of these sources show markedly elongated distributions of positions on the sky measured with diurnal observations. Here we apply a novel statistical and data-processing method to the diurnal position measurements stretching over 40 years to quantify the degree of elongation and its position angle, for each source with more than 200 data points. We find that 49\% of the examined sources have distribution elongations in excess of 1.3. Robust uncertainties of the directions of maximal astrometric dispersion are computed by the bootstrapping method, and the results are compared with a larger catalog of radio jet directions by Plavin et al. 2022. Nearly one-half of the sources with smaller position angle uncertainties are found to have astrometric position excursions from their mean positions aligned with the radio jet structures within $\pm 30\degr$.

The AU Microscopii planetary system is only 24 Myr old, and its geometry may provide clues about the early dynamical history of planetary systems. Here, we present the first measurement of the Rossiter-McLaughlin effect for the warm sub-Neptune AU\,Mic\,c, using two transits observed simultaneously with VLT/ESPRESSO, CHEOPS, and NGTS. After correcting for flares and for the magnetic activity of the host star, and accounting for transit-timing variations, we find the sky-projected spin-orbit angle of planet c to be in the range $\lambda_c=67.8_{-49.0}^{+31.7}$ degrees (1-$\sigma$). We examine the possibility that planet c is misaligned with respect to the orbit of the inner planet b ($\lambda_b=-2.96_{-10.30}^{+10.44}$ degrees), and the equatorial plane of the host star, and discuss scenarios that could explain both this and the planet's high density, including secular interactions with other bodies in the system or a giant impact. We note that a significantly misaligned orbit for planet c is in some degree of tension with the dynamical stability of the system, and with the fact that we see both planets in transit, though these arguments alone do not preclude such an orbit. Further observations would be highly desirable to constrain the spin-orbit angle of planet c more precisely.

The search for exoplanets has led to the identification of intriguing patterns in their distributions, one of which is the so-called sub-Jovian and Neptune desert. The occurrence rate of Neptunian exoplanets with an orbital period $P\lesssim 4$ days sharply decreases in this region in period-radius and period-mass space. We present a novel approach to delineating the sub-Jovian and Neptune desert by considering the incident stellar flux $F$ on the planetary surface as a key parameter instead of the traditional orbital period of the planets. Through this change of perspective, we demonstrate that the incident flux still exhibits a paucity of highly irradiated Neptunes, but also captures the proximity to the host star and the intensity of stellar radiation. Leveraging a dataset of confirmed exoplanets, we performed a systematic analysis to map the boundaries of the sub-Jovian and Neptune desert in the $(F,R_p)$ and $(F,M_p)$ diagrams, with $R_p$ and $M_p$ corresponding to the planetary radius and mass, respectively. By using statistical techniques and fitting procedures, we derived analytical expressions for these boundaries that offer valuable insights into the underlying physical mechanisms governing the dearth of Neptunian planets in close proximity to their host stars. We find that the upper and lower bounds of the desert are well described by a power-law model in the $(F,R_p)$ and $(F,M_p)$ planes. We also obtain the planetary mass-radius relations for each boundary by combining the retrieved analytic expressions in the two planes. This work contributes to advancing our knowledge of exoplanet demographics and to refining theoretical models of planetary formation and evolution within the context of the sub-Jovian and Neptune desert.

Understanding coronal structure and dynamics can be facilitated by analyzing green-line emission, which enables the investigation of diverse coronal structures such as coronal loops, streamers, coronal holes, and various eruptions in the solar atmosphere. In this study, we investigated the spatiotemporal behaviors of green-line emissions in both low and high latitudes across nine solar cycles, ranging from cycle 17 to the current cycle 25, using the Modified Homogeneous Data Set (MHDS). We employed methodologies such as cross-correlation, power spectral density (PSD), and wavelet transform techniques for this analysis. We found distinct behaviors in green line energy across various latitudinal distributions in the solar atmosphere. The trends observed at higher latitudes differ from those at lower latitudes. The emission behaviors show a close association with other solar phenomena like solar flares, sunspots, and coronal mass ejections (CMEs) throughout the solar cycles. The observed variations exhibit harmonic periods. The emission activity is significantly higher in the low latitudes, accounting for over 70 percent of the emissions, while the higher latitudes contribute less than 30 percent. The emissions exhibit asymmetric behavior between the northern and southern hemispheres, leading to a 44-year cycle of solar hemispheric dominance shifts. Various factors, such as Alfvén waves, solar magnetic fields, sunspots, differential rotation, and reconnection events, influence the observed differences in behavior between lower and higher latitudes, suggesting the existence of potential underlying phenomena contributing to deviations in properties, intensity, temporal dynamics, and spatiotemporal lifetime.

Aims. To evaluate the maximum attainable energies of electrons accelerated by means of the magneto-centrifugal mechanism. We examine how the range of maximum possible energies, as well as the primary limiting factors, vary with black hole mass. Additionally, we analyze the dependence of the maximum relativistic factor on an initial distance from the black hole. Methods. To model the acceleration of electrons on rotating magnetic field lines we apply several constraining mechanisms: the inverse Compton scattering, curvature radiation and the breakdown of the bead-on-the-wire approximation. Results. The maximal Lorentz factors for electron acceleration vary with the type of a black hole. For stellar-mass black holes, electrons can be accelerated up to the Lorentz factors 2 * 10^(6) - 2 * 10^(8) with only co-rotation constrain affecting the maximal relativistic factor; In intermediate-mass black holes, the Lorentz factors are in the interval 2 * 10^(8) - 2 * 10^(11); For the supermassive black holes the Lorentz factors range from 2.5 * 10^(10) to 2 * 10^(15); while the ultra-massive black hole located at the center of Abell 1201 can accelerate electrons up to 1.1 * 10^(13) - 6.6 * 10^(16). with both the co-rotation and curvature radiation determining the final Lorentz factor for the last three categories

Kilonova spectra provide us with the direct information of r-process nucleosynthesis in neutron star mergers. In this paper, we study the signatures of elements beyond the third r-process peak expected to be produced in neutron-rich ejecta in the photospheric spectra of kilonova. Ra II, Ac III, and Th III are our candidates because they have a small number of valence electrons and low-lying energy levels, which tend to result in strong absorption features. We systematically calculate the strength of bound-bound transitions of these candidates by constructing the line list based on the available atomic database. We find that Th III is the most promising species showing strong transitions at the NIR wavelengths. By performing radiative transfer simulations, we find that Th III produces broad absorption features at ~18000 A in the spectra when the mass ratio of actinides to lanthanides is larger than the solar r-process ratio and the mass fraction of lanthanides is $\lesssim 6\times10^{-4}$. Our models demonstrate that the Th feature may be detectable if the bulk of the ejecta in the line-forming region is dominated by relatively light r-process elements with the mixture of a small fraction of very neutron-rich material. Such conditions may be realized in the mergers of unequal-mass neutron stars or black hole-neutron star binaries. To detect the Th absorption features, the observations from the space (such as JWST) or high-altitude sites are important as the wavelength region of the Th features is overlapped with that affected by the strong telluric absorption.

Recent findings from photometric and spectroscopic JWST surveys have identified examples of high-redshift galaxies at $z \gtrsim 10$. These high-$z$ galaxies appear to form much earlier and exhibit greater UV luminosity than predicted by theoretical work. In this study, our goal is to reproduce the brightness of these sources by simulating high-redshift galaxies with virial masses $M_{\rm vir} = 10^{9} - 10^{10} M_{\odot}$ at $z > 10$. To achieve this, we conduct cosmological hydrodynamic zoom-in simulations, modifying baryonic sub-grid physics, and post-process our simulation results to confirm the observability of our simulated galaxies. Specifically, we enhanced star formation activity in high-redshift galaxies by either increasing the star formation efficiency up to 100\% or adopting a top-heavy initial mass function (IMF). Our simulation results indicate that both increasing star formation efficiency and adopting a top-heavy IMF play crucial roles in boosting the UV luminosity of high-redshift galaxies, potentially exceeding the limiting magnitude of JWST surveys in earlier epochs. Especially, the episodic starburst resulting from enhanced star formation efficiency may explain the high-redshift galaxies observed by JWST, as it evacuates dust from star-forming regions, making the galaxies more observable. We demonstrate this correlation between star formation activity and dust mass evolution within the simulated galaxies. Also, adopting a top-heavy IMF could enhance observability due to an overabundance of massive stars, although it may also facilitate rapid metal enrichment. Using our simulation results, we derive several observables such as effective radius, UV slope, and emission line rates, which could serve as valuable theoretical estimates for comparison with existing spectroscopic results and forthcoming data from the JWST NIRSpec and MIRI instruments.

We present photometric and spectroscopic observations of supernova (SN) 2014C, primarily emphasizing the initial month after the explosion at approximately daily intervals. During this time, it was classified as a Type Ib SN exhibiting a notably higher peak luminosity ($L_{\rm peak}\approx4.3\times10^{42}\rm erg\,s^{-1}$), a faster rise to brightness ($t_{\rm rise} \approx 11.6$ d), and a more gradual dimming ($\Delta m_{15}^{V} \approx 0.48$ mag) compared to typical SNe Ib. Analysis of the velocity evolution over the first $\sim$ 20 days after the explosion supports the view that the absorption near 6200Å\,is due to high-velocity H$\alpha$ in the outer layers of the ejecta, indicating the presence of a small amount of hydrogen in the envelope of progenitor before the explosion. Assuming the peak luminosity is entirely attributed to radioactive decay, we estimate that 0.14 \Msun\,of \Nifs\,was synthesized in the explosion. However, this amount of nickel could no longer maintain observed brightness approximately ten days after peak luminosity, suggesting additional energy sources beyond radioactive decay. This supplementary energy likely originates from interaction with the circumstellar medium (CSM). Consequently, the timing of the SN-CSM interaction in SN 2014C may occur much earlier than the emergence of IIn-like features during the nebular phase.

Magnetic hot stars can emit both coherent and incoherent non-thermal radio emission. Understanding the nature of these emissions and their connection to stellar rotation and magnetic field characteristics remains incomplete. The RAdio Magnetospheres of B and O stars (RAMBO) project aims to address this gap by systematically detecting and characterizing gyrosynchrotron and cyclotron maser radio emission in rapidly rotating magnetic hot stars. Using the upgraded Giant Metrewave Radio Telescope, we present the first detection of radio emission from HD55522 at 650 MHz, confirming it as a new radio-bright magnetic hot star. This supports the predictions of the Centrifugal Breakout model, furthering its application in understanding particle acceleration mechanisms in centrifugal magnetospheres of hot stars. Additionally, we report non-detections for four other targets, improving sensitivity limits by a factor of a few compared to previous observations. These findings demonstrate the potential of RAMBO to uncover the complexities of radio emission in massive stars and highlight the need for broader, multi-wavelength observations to probe magnetospheric physics comprehensively. The sensitivity of the Square Kilometre Array will enable significant advancements.

Morbidelli, Kleine & Nimmo (2024) (MKN) recently published a critical analysis on whether the terrestrial planets in the Solar System formed by rapid pebble accretion or by the classical route of multiple giant impacts between planetary embryos after the dissipation of the protoplanetary disc. They arrive at the conclusion that the terrestrial planets did not form by pebble accretion. Although we welcome debate on this topic, we want to emphasize here several points that we disagree on. We will not address in detail every claim made in MKN, but rather stick to four main points. Our conclusion is that pebble accretion remains a viable mechanism to drive significant growth of protoplanets in the protoplanetary disc, with as much as 70% of Earth formed by pebble accretion. This rapid growth phase must nevertheless have been followed by an extended period of collisional growth after the end of the protoplanetary disc phase, likely culminating with the moon-forming giant impact. We emphasize here an important recent result from Sharp & Olson (2023), namely that significant growth by pebble accretion can be reconciled with the Hf-W decay system even for a canonical moon-forming giant impact with a Mars-mass protoplanet - a near-equal mass impactor, as proposed in Johansen et al. (2023), is not necessary. Given that terrestrial planet formation naturally involves both pebble accretion and a combination of small and large impactors, this challenges the very notion of making an either/or distinction between the classical collision model and the pebble accretion model.

Extragalactic Black Hole X-ray Binaries (BH-XRBs) are the most intriguing X-ray sources as some of them are `home' to the most massive stellar-mass BHs ever found. In this work, we conduct a comprehensive study of three massive, eclipsing extragalactic BH-XRBs i.e., M33X-7, NGC300X-1, and IC10X-1 and using entire X-ray observations available from \textit{XMM-Newton} and \textit{NuSTAR} till date. Preliminary analysis using \textit{diskbb} and \textit{powerlaw} models shows that the sources have steep spectra and sub-Eddington luminosities (L<0.69 L$_{Edd}$), with major flux contribution from non-thermal component, resembling the relatively uncharted Steep Powerlaw State (SPL). To understand the accretion disc properties in this state, we explore alternate modelling scenario that reveals the presence of a `hot' ($kT_{in}=1-2$ keV) slim-disc (\textit{diskpbb}) with radial temperature profile $T(r)\propto r^{-p}$ ($p=0.5-0.66$), along with a cooler ($kT_{in}=0.1-0.2$ keV) standard thermal disc (\textit{diskbb}). We carry out the continuum-fitting method using relativistic slim-disc model (\textit{slimbh}) and estimate the mass range of M33 X-7, NGC300X-1 and IC10X-1 is to be 9$-$15 M$_{\odot}$, 9$-$28 M$_{\odot}$ and 10$-$30 M$_{\odot}$, respectively. Further, eclipse periods are determined by modelling the lightcurve, using which we estimate the size of the eclipsing bodies. Modelling of the eclipse spectra revealed the complete obscuration of soft spectral component during eclipse, implying the emission of hard component from an extended accretion region. Based on our findings, we provide an inference on geometry of accretion disc in these wind-fed systems and compare their properties with the other two extragalactic BH-XRBs.

We present an analysis of the spatial distribution of globular cluster (GC) systems of 118 nearby early-type galaxies in the Next Generation Virgo Cluster Survey (NGVS) and Mass Assembly of early-Type GaLAxies with their fine Structures (MATLAS) survey programs, which both used MegaCam on the Canada-France-Hawaii Telescope. We describe the procedure used to select GC candidates and fit the spatial distributions of GCs to a two-dimensional Sérsic function, which provides effective radii (half number radii) and Sérsic indices, and estimate background contamination by adding a constant term to the S'ersic function. In cases where a neighboring galaxy affects the estimation of the GC spatial distribution in the target galaxy, we fit two 2D Sérsic functions, simultaneously. We also investigate the color distributions of GCs in our sample by using Gaussian Mixture Modeling. For GC systems with bimodal color distributions, we divide the GCs into blue and red subgroups and fit their respective spatial distributions with Sérsic functions. Finally, we measure the total number of GCs based on our fitted Sérsic function, and calculate the GC specific frequency.

Exploring the formation and evolution of second-generation circumbinary discs around evolved binary stars, such as post-Asymptotic Giant Branch (post-AGB) and post-Red Giant Branch (post-RGB) binaries, provides valuable insights into the complex binary interaction process that concludes the red-giant phase of evolution in these systems. Additionally, it offers a novel opportunity to investigate the formation of second-generation planets within dusty discs surrounding evolved stars. We present a pilot multi-wavelength polarimetric imaging study of the post-AGB binary system IRAS 08544-4431 using the European Southern Observatory-Very Large Telescope/SPHERE instrument. This study is focused on optical V- and I'-band ZIMPOL data to complement near-infrared H-band IRDIS data presented previously. The study aims to investigate the dust scattering properties and surface morphology of the post-AGB circumbinary disc as a function of wavelength. We successfully resolved the extended disc structure of IRAS\,08544-4431, revealing a complex disc morphology, high polarimetric disc brightness (up to ~1.5%), and significant forward scattering at optical wavelengths. Additionally, we found that the disc shows a grey polarimetric colour in both optical and near-infrared. The findings highlight similarities between post-AGB circumbinary discs and protoplanetary discs, suggesting submicron-size porous aggregates as the dominant surface dust composition, and indicating potential warping within the disc. However, further expansion of the multi-wavelength analysis to a larger sample of post-AGB binary systems, as well as high-resolution observations of dust continuum and gas emission, is necessary to fully explore the underlying structure of post-AGB circumbinary discs and associated physical mechanisms.

Fermi balls are non-topological solitons that can naturally form in an early universe containing a dark sector with heavy fermions and an attractive interaction mediated by a light scalar field. We compute the Fermi ball mass and radius scaling relations when the potential of the scalar field $\varphi$ has a non-negligible quartic coupling $\lambda\varphi^4$. The resulting Fermi balls reach `saturation' very rapidly, even when their radius is much smaller than the effective Yukawa force range. These objects can therefore grow by mergers or by accretion of ambient dark fermions, until they become so dense that they fall within their Schwarzschild radius and collapse to black holes. This setup, therefore, provides an example of a rather natural and economical dark sector scenario for the formation of primordial black holes.

Due to the lack of experimental data on extremely neutron-rich nuclei, theoretical values derived from nuclear physics models are essential for the rapid neutron capture process ($r$-process). Metal-poor stars enriched by the $r$-process offer valuable cases for studying the impact of nuclear physics models on $r$-process nucleosynthesis. This study analyzes four widely used nuclear physics models in detail: Finite-Range Droplet Model, Hartree-Fock-Bogoliubov, Duflo-Zuker, and Weizs$\ddot{\rm a}$cker-Skyrme (WS4). Theoretical values predicted by the WS4 model are found to be in good agreement with experimental data, with deviations significantly smaller than those predicted by other models. The heavy element abundances observed in $r$-process enhanced metal-poor stars can be accurately reproduced by $r$-process nucleosynthesis simulations using the WS4 model, particularly for the rare earth elements. This suggests that nuclear data provided by nuclear physics model like WS4 are both essential and crucial for $r$-process nucleosynthesis studies.

We present a study on the evolution of transitional dwarf galaxies, specifically dwarf lenticulars (dS0s) and early-type dwarfs with blue cores (ETdG(bc)s), driven by environmental processes in the Virgo cluster utilizing the Extended Virgo Cluster Catalog. We investigated the morphological fraction and stellar mass of transitional dwarf galaxies in relation to the clustercentric distance, compared to dwarf elliptical galaxies (dEs) and dwarf irregular galaxies (dIrrs). We found that dS0s beyond 0.7R_vir exhibit a similar trend in the morphology-clustercentric distance relation to dEs, demonstrating a decreasing fraction with clustercentric distance, whereas ETdG(bc)s display an opposite trend to dS0s. The spatial distributions of transitional dwarf galaxies and dEs correlate with the mass, in which fractions of bright, massive galaxies increase towards the central region of the Virgo cluster. In the mass-clustercentric distance plane, dS0s exhibit a skewed distribution that favors more massive galaxies than dEs at a given clustercentric distance. In the projected phase-space diagram, dS0s are scarce in the stripped region, whereas ETdG(bc)s are absent in both the stripped and virialized regions. In addition, the dS0s in the virialized region are predominantly brighter and more massive than the dEs, indicating that the transformation of dS0s into dEs depends on the stellar mass. We propose that the majority of observed dS0s constitute a population that has settled into the Virgo cluster, whereas ETdG(bc)s represent a recently accreted population. We discuss the impact of ram pressure stripping effects on mass-dependent morphological evolution.

The 21-cm forest, composed of spectral absorption features from high-redshift background radio sources, provides a unique probe for studying small-scale structures during the epoch of reionization. It is particularly sensitive to detecting small-scale structures and early heating processes. Despite the rich information contained in the 21-cm forest signal, the complexity of directly modeling the signal has led to a lack of effective analytical models. However, the one-dimensional (1D) power spectrum of the 21-cm forest contains valuable information about the matter power spectrum, making analytical modeling feasible. This work employs an analytical modeling approach based on the halo model, which links the distribution of matter to dark matter halos, allowing for effective predictions of cosmic structure formation and its impact on the 21-cm signal. By considering various parameter scenarios within the halo model framework, particularly different dark matter particle masses and varying levels of cosmic heating, we can capture the complexities of small-scale structures and make the 1D power spectrum modeling applicable across a wide range of parameters. This method not only enhances our understanding of the 21-cm forest signal but also provides theoretical support for future observational data. Observing the 21-cm forest with large radio telescopes, such as the Square Kilometre Array, is anticipated to enable simultaneous exploration of dark matter properties and the heating history of the early universe.

The existence of the ultracompact object XTE J1814-338, with an inferred mass and radius of $M$ = 1.21 $\pm$ 0.05 $M_\odot$ and R = 7.0 $\pm$ 0.4 km, presents a great challenge for the theory of neutron stars. Within this context, we revisit the theory of dark-matter-admixed neutron stars and infer the physical properties of this compact object, such as the Fermi momentum of dark matter necessary to compress the star at such low radius, its equation of state, speed of sound, and some macroscopic properties, such as the inertia moment and the dimensionless tidal parameter. We also compare the physical properties of the XTE J1814-338 with other pulsars, such as the canonical 1.4 M$_\odot$, the PSR J0740 + 6620, and the HESS J1731-347.

European Space Agency's Rosetta mission is the only space mission that performed long-term monitoring of comet at close distances. Its over two years' rendezvous with comet 67P/Churyumov-Gerasimenko revealed diverse evolutionary processes of the cometary nucleus. One of the most striking events is the migration of a 30-m boulder in the southern hemisphere region of Khonsu. Previous works found the boulder's 140-m displacement occurred during the three months from August to October 2015, and several triggering mechanisms were proposed, including outburst at the boulder site, seismic vibrations from nearby activities, or surface erosion of the slope beneath the boulder. In this work, we further analyze this impressive event by analysing imaging data from Rosetta's OSIRIS camera. We constrained the boulder's migration time to within 14 hours and derived a detailed timeline of the boulder migration event and local dust activities. High-resolution thermophysical modelling shows significant dichotomy in the thermal history of the boulder's southern and northern sides, which could have triggered or facilitated its migration via its own volatile activity.

While protoplanetary disks (PPDs) are generally thought to dissipate within several Myr, recent observations have revealed gas in debris disks. The origin of this gas remains uncertain, with one possibility being the unexpectedly long survival of PPDs (the primordial-origin scenario). To explore the plausibility of this scenario, we conduct 1D disk evolution simulations, varying parameters like stellar mass, disk mass, turbulent stress, and magnetohydrodynamic winds, while incorporating stellar evolution to account for time-varying photoevaporation rates. Our focus is on disks where small grains are depleted, as these are potentially long-lived due to reduced far-ultraviolet photoevaporation. Our results show that gas in these disks can survive beyond 10 Myr regardless of the stellar mass, provided they are initially massive ($M_{\mathrm{disk}}\approx 0.1M_*$) with relatively weak turbulent stress ($\alpha \ll 10^{-2}$). The longest lifetimes are consistently found for $M_* = 2 M_{\odot}$ across a wide parameter space, with gas typically persisting at $\sim 10$--$10^3 \mathrm{au}$. Roughly estimated CO masses for these disks fall within the observed range for the most massive gas-rich debris disks around early A stars. These alignments support the plausibility of the primordial-origin scenario. Additionally, our model predicts that accretion persists for as long as the disk survives, which could explain the accretion signatures detected in old disks hosted by low-mass stars, including Peter Pan disks. Our finding also suggests that ongoing accretion may exist in gas-rich debris disks. Thus, searching for accretion signatures could be a key factor to identifying the origin of gas in debris disks.

Modern radio telescopes generate large amounts of data, with the next generation Very Large Array (ngVLA) and the Square Kilometre Array (SKA) expected to feed up to 292 GB of visibilities per second to the science data processor (SDP). However, the continued exponential growth in the power of the world's largest supercomputers suggests that for the foreseeable future there will be sufficient capacity available to provide for astronomers' needs in processing 'science ready' products from the new generation of telescopes, with commercial platforms becoming an option for overflow capacity. The purpose of the current work is to trial the use of commercial high performance computing (HPC) for a large scale processing task in astronomy, in this case processing data from the GASKAP-HI pilot surveys. We delineate a four-step process which can be followed by other researchers wishing to port an existing workflow from a public facility to a commercial provider. We used the process to provide reference images for an ongoing upgrade to ASKAPSoft (the ASKAP SDP software), and to provide science images for the GASKAP collaboration, using the joint deconvolution capability of WSClean. We document the approach to optimising the pipeline to minimise cost and elapsed time at the commercial provider, and give a resource estimate for processing future full survey data. Finally we document advantages, disadvantages, and lessons learned from the project, which will aid other researchers aiming to use commercial supercomputing for radio astronomy imaging. We found the key advantage to be immediate access and high availability, and the main disadvantage to be the need for improved HPC knowledge to take best advantage of the facility.

This study examines the average X-ray properties of massive halos at z< 0.2, covering the largest halo mass range to date, from Milky Way-like halos to massive clusters. The analysis is based on stacking in the eFEDS area of the GAMA galaxy group sample, validated with synthetic data that mimic observed eROSITA X-ray and GAMA optical data using Magneticum lightcones. Stacking was conducted in halo mass bins and tested for AGN and X-ray binary contamination, systematics in the halo mass proxy, and uncertainties in optical group centers. The study provides average X-ray surface brightness profiles in six mass bins, spanning Milky Way-like systems to poor clusters. The scatter in the X-ray luminosity-mass (LX-M) relation is attributed to gas concentration: low X-ray luminosity systems at fixed halo mass exhibit lower central gas concentrations than high-luminosity systems, consistent with Magneticum predictions. However, discrepancies in dark matter concentration arise, with Magneticum predicting undetected groups as older and more relaxed, while observations suggest the opposite. New LX-M relations are presented covering three decades of halo mass. These relations fit a single power law, aligning with previous studies. Magneticum matches observed gas distributions across all masses, whereas IllustrisTNG, EAGLE, Simba, and FLAMINGO exhibit significant discrepancies at various mass scales. Simulations calibrated on local galaxy properties accurately reproduce central galaxies but fail to capture gas properties. Conversely, simulations like Magneticum excel in gas predictions but produce overly massive central galaxies. Further exploration of gas and dark matter distributions and their effects on galaxy properties is critical to comprehending the role of gravitational forces and feedback in shaping large-scale structure and galaxy evolution.

In the detection of gravitational waves in space, during the science phase of the mission, the point ahead angle mechanism (PAAM) serves to steer a laser beam to compensate for the angle generated by the relative motion of the two spacecrafts (SCs) during the approximately 10 seconds of flight time a laser beam will take from one SC to reach a distant SC of three million kilometers away. The common practice for pointing stability control of a laser beam is to first do a coarse tracking by the PAAM to steer a laser beam to compensate for the relative motion between two SCs, to be followed by a fine pointing stability control. In the present work, by exploiting the near-circular orbit structure of individual SC in the triangular constellation, the feasibility of inserting an adaptive Kalman filter (AEKF) into the PAAM control loop is investigated. By adopting a colored measurement noise model that closely resembles the prospective on orbit situation, numerical simulation suggests that the dynamic range of the PAAM may be reduced to the level of nano-radians using the prediction of the pointing head angle (PAA) by the AEKF. This will cut down on the TTL coupling noise and the position noise budget allocated to the PAAM. This in turn reduces the dynamic range of the fine pointing control and leaves room to improve its accuracy, thereby offers the prospect of reduction of the position noise budget allocated to the laser pointing instability as a whole.

We continue to study the structure and kinematics of HH flows. Herbig-Haro (HH) flows exhibit large variety of morphological and kinematical structures. Both proper motion (PM) and radial velocity investigations are essential to understand the physical nature of such structures. We investigate the kinematics and PM of spectrally separated structures in the PV Cep HH flow HH 215. We present the observational results obtained with a 6 m telescope (Russia) using the SCORPIO multi-mode focal reducer with scanning Fabry-Perot interferometer. Two epochs of the observations of the PV Cep region in H$\alpha$ and [SII] emission (2003 and 2020-2021) allowed us to study the morphology of HH 215 jet in detail and to measure the PM and the radial velocities for its inner structures. Already known emission knots in the HH 215 flow and new features were studied. Moreover, a newly-formed HH knot was revealed, presumably formed during the large maximum of PV Cep star in 1976-1977. We found the high-velocity inner channel in the HH 215 ionized outflow, oriented accordingly to the mean direction of the whole HH outflow and the axis of the symmetry of the reflection nebula. The HH-knots located along the axis of the high-velocity channel have a position angle coinciding with its axis (abut 325$^{\circ}$), however other ones have completely different value (about 25$^{\circ}$), which supports the idea that those knots are formed by oblique shocks. We derived the value of i $\approx$ 30$^{\circ}$$\pm$ 5$^{\circ}$ for the inclination angle between the flow axis and the line of sight. The total length of HH 215 outflow should be about 0.2 pc, and the full length of the bipolar outflow from PV Cep (HH 315 + HH 215) can be estimated as 3.6 pc, assuming that it more or less keeps the same inclination angle.

The circumgalactic medium (CGM) is the largest baryon reservoir around galaxies, but its extent, mass, and temperature distribution remain uncertain. We propose that cold gas in the CGM resides primarily in $\sim 100 \hbox{--} 10^4$ cloud complexes (CCs), each containing a mist of tiny cold cloudlets dispersed in a warm/hot medium ($\sim 10^5 \hbox{--} 10^6$~K). Modeling CCs as uniform and misty simplifies the calculation of observables like ion absorption columns compared to resolving tiny individual cloudlets. Using Monte Carlo simulations, we explore how CC properties affect the observed spread in column densities. A power-law distribution of CCs ($dN_{\rm CC}/dR \propto R^{-1.2}$) reproduces MgII equivalent width (EW) and column density distribution trends with impact parameter ($R_\perp$). We show that the area-averaged MgII column density, combined with the covering fraction of CCs, provides a robust proxy for estimating the cold CGM mass, independent of other model parameters. Modeling individual CCs demonstrates that turbulent broadening blends cloudlet absorption lines, allowing CCs to approximate the observational effects of their constituent cloudlets analytically. Direct simulations of cloudlets within multiple CCs confirm the computational challenges of fully resolving the mist-like structure. Comparing modeled MgII absorption with observations, we estimate the cold CGM mass of Milky Way-like galaxies to be $\sim 10^{10} \, M_\odot$, about $10\%$ of the total CGM mass. This work provides a practical framework for connecting CGM models with observations, shedding light on the cold gas distribution in galaxy halos and its role in the galactic baryon cycle.

Black hole (BH) mergers are natural sources of gravitational waves (GWs) and are possibly associated with electromagnetic events. Such events from a charged rotating BH with an accretion on to it could be more energetic and ultra-short-lived if the magnetic force dominates the accretion process because the attraction of ionized fluid with a strong magnetic field around the rotating BH further amplifies the acceleration of the charged particle via a gyromagnetic effect. Thus a stronger magnetic field and gravitational pull will provide an inward force to any fluid displaced in the radial direction and move it toward the axis of rotation with an increasing velocity. After many twists during rotation and the existence of restoring agents, Such events could produce a narrow intense jet starts in the form of Poynting flux along the axis of rotation resembling the Blandford-Znajek (BZ) mechanism. We investigated a charged rotating BH and obtained characteristic results (e.g., the remnant mass, magnetic field strength, luminosity, opening angle, viewing angle, and variation of viewing angle on the SGRB luminosity detection) that have a nice coincidence with rare events having GW associated with EM counterparts. This study gives a new insight into events with a strongly magnetized disk dominating the accretion process of energy extraction.

Understanding the processes within compact stars hinges on astrophysical observations. A recent study reported on the central object in the HESS J1731-347 supernova remnant (SNR), estimating a mass of $M=0.77_{-0.17}^{+0.20} \ M_{\odot}$ and a radius of $R=10.40_{-0.78}^{+0.86} \ \rm{km}$, making it the lightest neutron star ever observed. Conventional models suggest that neutron stars form with a minimum gravitational mass of about $1.17M_{\odot}$, raising the question: is this object a typical neutron star, or could it be our first encounter with an "exotic" star? To explore this, we employ the Color-Flavored Locked (CFL) equation of state (EoS), aiming to constrain it by integrating data from the HESS J1731-347 event with pulsar observations and gravitational wave detections. Additionally, we model hybrid EoS by combining the MDI-APR1 (hadronic) and CFL (quark) EoS, incorporating phase transitions via Maxwell construction. Our analysis indicates that CFL quark matter adequately explains all measurements, including the central compact object of HESS J1731-347. In contrast, hybrid models featuring CFL quark phases fail to account for the masses of the most massive observed pulsars.

The short gamma-ray burst (sGRB), GRB~170817A, is often considered a rare event. However, its inferred event rate, $\mathcal{O}(100s)\ \text{Gpc}^{-3}\ \text{yr}^{-1}$, exceeds cosmic sGRB rate estimates from high-redshift samples by an order of magnitude. This discrepancy can be explained by geometric effects related to the structure of the relativistic jet. We first illustrate how adopting a detector flux threshold point estimate rather than an efficiency function, can lead to a large variation in rate estimates. Simulating the Fermi-GBM sGRB detection efficiency, we then show that for a given a universal structured jet profile, one can model a geometric bias with redshift. Assuming different jet profiles, we show a geometrically scaled rate of GRB~170817A is consistent with the cosmic beaming uncorrected rate estimates of short $\gamma$-ray bursts (sGRBs) and that geometry can boost observational rates within $\mathcal{O}(100s)$\,Mpc. We find an apparent GRB~170817A rate of $303_{-300}^{+1580}$ $\mathrm{Gpc}^{-3}\, \mathrm{yr}^{-1} $ which when corrected for geometry yields $6.15_{-6.06}^{+31.2}$ $\mathrm{Gpc}^{-3}\, \mathrm{yr}^{-1} $ and $3.34_{-3.29}^{+16.7}$ $\mathrm{Gpc}^{-3}\, \mathrm{yr}^{-1} $ for two different jet profiles, consistent with pre-2017 estimates of the isotropic sGRB rate. Our study shows how jet structure can impact rate estimations and could allow one to test structured jet profiles. We finally show that modelling the maximum structured jet viewing angle with redshift can transform a cosmic beaming uncorrected rate to a representative estimate of the binary neutron star merger rate. We suggest this framework can be used to demonstrate parity with merger rates or to yield estimates of the successful jet fraction of sGRBs.

Galaxies evolve hierarchically through merging with lower-mass systems and the remnants of destroyed galaxies are a key indicator of the past assembly history of our Galaxy. However, accurately measuring the properties of the accreted galaxies and hence unraveling the Milky Way's (MW) formation history is a challenging task. Here we introduce CASBI (Chemical Abundance Simulation Based Inference), a novel inference pipeline for Galactic Archeology based on Simulation-based Inference methods. CASBI leverages on the fact that there is a well defined mass-metallicity relation for galaxies and performs inference of key galaxy properties based on multi-dimensional chemical abundances of stars in the stellar halo. Hence, we recast the problem of unraveling the merger history of the MW into a SBI problem to recover the properties of the building blocks (e.g. total stellar mass and infall time) using the multi-dimensional chemical abundances of stars in the stellar halo as observable. With CASBI we are able to recover the full posterior probability of properties of building blocks of Milky Way like galaxies. We highlight CASBI's potential by inferring posteriors for the stellar masses of completely phase mixed dwarf galaxies solely from the 2d-distributions of stellar abundance in the iron vs. oxygen plane and find accurate and precise inference results.

This paper introduces a hybrid numerical scheme for the fuzzy dark matter model: It combines a wave-based approach to solve the Schrödinger equation using Fourier continuations with Gram polynomials and a fluid-based approach to solve the Hamilton-Jacobi-Madelung equations. This hybrid scheme facilitates zoom-in simulations for cosmological volumes beyond the capabilities of wave-based solvers alone and accurately simulates the full nonlinear dynamics of fuzzy dark matter. We detail the implementation of a Hamilton-Jacobi-Madelung solver, the methodology for phase matching at fluid-wave boundaries, the development of a local pseudospectral wave solver based on Fourier continuations, new grid refinement criteria for both fluid and wave solvers, an interpolation algorithm based on Fourier continuations, and the integration of these building blocks into the adaptive mesh refinement code GAMER. The superiority of the scheme is demonstrated through various performance and accuracy tests, tracking the linear power spectrum evolution in a 10 Mpc/h box, and a hybrid cosmological simulation in a 5.6 Mpc/h box. The corresponding code is published as part of the GAMER project on this https URL.

We present the discovery of a new Odd Radio Circle (ORC J0219--0505) in 1.2~GHz radio continuum data from the MIGHTEE survey taken with the MeerKAT telescope. The radio-bright host is a massive elliptical galaxy, which shows extended stellar structure, possibly tidal tails or shells, suggesting recent interactions or mergers. The radio ring has a diameter of 35", corresponding to 114~kpc at the host galaxy redshift of $z_{\rm spec} = 0.196$. This MIGHTEE ORC is a factor 3--5 smaller than previous ORCs with central elliptical galaxies. The discovery of this MIGHTEE ORC in a deep but relatively small-area radio survey implies that more ORCs will be found in deeper surveys. While the small numbers currently available are insufficient to estimate the flux density distribution, this is consistent with the simplest hypothesis that ORCs have a flux density distribution similar to that of the general population of extragalactic radio sources.

We present measurements of the specific angular momentum $j_\star$ of 41 star-forming galaxies at $1.5<z<2.5$. These measurements are based on radial profiles inferred from near-IR \textit{HST} photometry, along with multi-resolution emission-line kinematic modelling using integral field spectroscopy (IFS) data from KMOS, SINFONI, and OSIRIS. We identified 24 disks (disk fraction of $58.6\pm 7.7\%$) and used them to parametrize the $j_\star$ \textit{vs} stellar mass $M_\star$ relation (Fall relation) as $j_\star\propto M_\star^{\beta}$. We measure a power-law slope $\beta=0.25\pm0.15$, which deviates by approximately $3\sigma$ from the commonly adopted local value $\beta = 0.67$, indicating a statistically significant difference. We find that two key systematic effects could drive the steep slopes in previous high-redshift studies: first, including irregular (non-disk) systems due to limitations in spatial resolution and second, using the commonly used approximation $\tilde{j}_\star\approx k_n v_s r_\mathrm{eff}$, which depends on global unresolved quantities. In our sample, both effects lead to steeper slopes of $\beta=0.48\pm0.21$ and $\beta=0.61\pm0.21$, respectively. To understand the shallow slope, we discuss observational effects and systematic uncertainties and analyze the retention of $j_\star$ relative to the angular momentum of the halo $j_h$ (angular momentum retention factor $f_j =j_\star/j_h$). For the $M_\star$ range covered by the sample $9.5 <\log_{10} (M_\star/M_\odot) < 11.5$ (halo mass $11.5 < \log_{10} (M_h/M_\odot) < 14$), we find large $f_j$ values ($>1$ in some cases) in low-mass haloes that decrease with increasing mass, suggesting a significant role of efficient angular momentum transport in these gas-rich systems, aided by the removal of low-$j_\star$ gas via feedback-driven outflows in low-mass galaxies.

We report on the results of X-ray observations (XMM-Newton, INTEGRAL and Swift) of two hard X-ray sources, IGR J17503-2636 and IGR J17507-2647, whose nature is not fully elucidated in the literature. Three XMM-Newton observations covered the field of IGR J17503-2636, in 2020 and twice in 2023. The analysis of the two XMM-Newton observations performed in September 2023, six days apart, did not detect IGR J17503-2636, allowing us to pose the most stringent 3sigma upper limit on the source flux to date (~9.5x10^-14 erg/cm2/s, 2-10 keV, flux corrected for absorption). This value implies that the amplitude of the X-ray flux variability exceeds a factor of ~2100, compared with the discovery outburst in 2018. A candidate X-ray periodicity at 0.335397(3) seconds has been barely detected (significance of about 3.8sigma) from IGR J17503-2636 with XMM-Newton (pulsed fraction of (10+/-1) per cent). The new data, put into the context of previous literature, allow us to propose a new classification of IGR J17503-2636 as a symbiotic X-ray binary, rather than a candidate supergiant fast X-ray transient. IGR J17507-2647 was formerly reported below 10 keV only during Chandra observations performed in 2009. We report here on two XMM-Newton observations that serendipitously covered the source field in 2020 and in 2023, finding a stable X-ray emission, both in X-ray flux and spectral shape. The long-term, persistent X-ray emission has also been probed by several Swift/XRT short observations and by INTEGRAL data spanning several years. We have detected an iron line in the emission (with centroid energy in the range of 6.3-6.6 keV), never reported before in the IGR J17507-2647 spectrum. The source properties favor the identification with a cataclysmic variable.

The stability criteria of rapid mass transfer and common-envelope evolution are fundamental in binary star evolution. They determine the mass, mass ratio, and orbital distribution of many important systems, such as X-ray binaries, type Ia supernovae, and merging gravitational-wave sources. In the limit of extremely rapid mass transfer, the response of a donor star in an interacting binary becomes asymptotically one of adiabatic expansion. We built the adiabatic mass-loss model and systematically surveyed the thresholds for dynamical timescale mass transfer over the entire span of possible donor star evolutionary states. Many studies indicate that new mass transfer stability thresholds play an essential role in the formation and properties of double compact object populations and the progenitors of SNe Ia and detectable GW sources. For example, our studies show that the mass transfer in the red giant and the asymptotic giant branch stars and the massive stars can be more stable than previously believed. Consequently, detailed binary population synthesis studies, using updated unstable mass transfer criteria, predicate the non-conservative stable mass transfer may dominate the formation channel of double stellar-mass black holes and can explain the population of the large mass ratio double stellar-mass black holes. Using our updated mass transfer thresholds, binary population thesis studies by Li et al. show that Ge et al.'s results support the observational double white dwarfs merger rate distribution per Galaxy and the space density of double white dwarfs in the Galaxy.

The Mid-infrared ELT Imager and Spectrograph (METIS) is a first-generation instrument for the Extremely Large Telescope (ELT), Europe's next-generation 39 m ground-based telescope for optical and infrared wavelengths. METIS will offer diffraction-limited imaging, low- and medium-resolution slit spectroscopy, and coronagraphy for high-contrast imaging between 3 and 13 microns, as well as high-resolution integral field spectroscopy between 3 and 5 microns. The main METIS science goals are the detection and characterisation of exoplanets, the investigation of proto-planetary disks, and the formation of planets. The Single-Conjugate Adaptive Optics (SCAO) system corrects atmospheric distortions and is thus essential for diffraction-limited observations with METIS. Numerous challenging aspects of an ELT Adaptive Optics (AO) system are addressed in the mature designs for the SCAO control system and the SCAO hardware module: the complex interaction with the telescope entities that participate in the AO control, wavefront reconstruction with a fragmented and moving pupil, secondary control tasks to deal with differential image motion, non-common path aberrations and mis-registration. A K-band pyramid wavefront sensor and a GPU-based Real-Time Computer (RTC), tailored to the needs of METIS at the ELT, are core components. This current paper serves as a natural sequel to our previous work presented in Hippler et al. (2018). It includes updated performance estimations in terms of several key performance indicators, including achieved contrast curves. We outline all important design decisions that were taken, and present the major challenges we faced and the main analyses carried out to arrive at these decisions and eventually the final design. We also elaborate on our testing and verification strategy, and, last not least, comprehensively present the full design, hardware and software.

Interior compositions are key for our understanding of Earth-like exoplanets. The composition of the core can influence the presence of a magnetic dynamo and the strength of gravity on the planetary surface, both of which heavily impact thermal and possible biological processes and thus the habitability for life and its evolution on the planet. However, detailed measurements of the planetary interiors are extremely challenging for small exoplanets, and existing data suggest a wide diversity in planet compositions. Hitherto, only certain photospheric chemical abundances of the host stars have been considered as tracers to explain the diversity of exoplanet compositions. Here we present a homogeneous analysis of stars hosting rocky exoplanets, with ages between 2 and 14 Gyr, revealing a correlation between rocky exoplanet compositions and the ages of the planetary systems. Denser rocky planets are found around younger stars. This suggests that the compositional diversity of rocky exoplanets can be linked to the ages of their host stars. We interpret this to be a result of chemical evolution of stars in the Milky Way, which modifies the material out of which stars and planets form. The results imply that rocky planets which form today, at similar galactocentric radii, may have different formation conditions, and thus different properties than planets which formed several billion years ago, such as the Earth.

Context. Relativistic jets from Active Galactic Nuclei are observed to be collimated on the parsec scale. When the pressure between the jet and the ambient medium is mismatched, recollimation shocks and rarefaction shocks are formed. Previous numerical simulations have shown that instabilities can destroy the recollimation structure of jets. Aims. In this study, we aim to study the instabilities of non-equilibrium over-pressured relativistic jets with helical magnetic fields. Especially, we investigate how the magnetic pitch affects the development of instabilities. Methods. We perform three-dimensional relativistic magnetohydrodynamic simulations for different magnetic pitches, as well as a two-dimension simulation and a relativistic hydrodynamic simulation served as comparison groups Results. In our simulations, Rayleigh-Taylor Instability (RTI) is triggered at the interface between the jet and ambient medium in the recollimation structure of the jet. We found that when the magnetic pitch decreases the growth of RTI becomes weak but interestingly, another instability, the CD kink instability is excited. The excitement of CD kink instability after passing the recollimation shocks can match the explanation of the quasi-periodic oscillations observed in BL Lac qualitatively.

In an ongoing effort to understand planet formation the link between the chemistry of the protoplanetary disk and the properties of resulting planets have long been a subject of interest. These connections have generally been made between mature planets and young protoplanetary disks through the carbon-to-oxygen (C/O) ratio. In a rare number of systems, young protoplanets have been found within their natal protoplanetary disks. These systems offer a unique opportunity to directly study the delivery of gas from the protoplanetary disk to the planet. In this work we post-process 3D numerical simulations of an embedded Jupiter-massed planet in its protoplanetary disk to explore the chemical evolution of gas as it flows from the disk to the planet. The relevant dust to this chemical evolution is assumed to be small, co-moving grains with a reduced dust-to-gas ratio indicative of the upper atmosphere of a protoplanetary disk. We find that as the gas enters deep into the planet's gravitational well, it warms significantly (up to $\sim 800$ K), releasing all of the volatile content from the ice phase. This change in phase can influence our understanding of the delivery of volatile species to the atmospheres of giant planets. The primary carbon, oxygen, and sulfur carrying ices: CO$_2$, H$_2$O, and H$_2$S are released into the gas phase and along with the warm gas temperatures near the embedded planets lead to the production of unique species like CS, SO, and SO$_2$ compared to the protoplanetary disk. We compute the column densities of SO, SO$_2$, CS, and H$_2$CS in our model and find that their values are consistent with previous observational studies.

In kinetic theory, the classic $n \Sigma v$ approach calculates the rate of particle interactions from local quantities: the number density of particles $n$, the cross-section $\Sigma$, and the average relative speed $v$. In stellar dynamics, this formula is often applied to problems in collisional (i.e. dense) environments such as globular and nuclear star clusters, where blue stragglers, tidal capture binaries, binary ionizations, and micro-tidal disruptions arise from rare close encounters. The local $n \Sigma v$ approach implicitly assumes the ergodic hypothesis, which is not well motivated for the densest star systems in the Universe. In the centers of globular and nuclear star clusters, orbits close into 1D ellipses because of the degeneracy of the potential (either Keplerian or harmonic). We find that the interaction rate in perfectly Keplerian or harmonic potentials is determined by a global quantity -- the number of orbital intersections -- and that this rate can be far lower or higher than the ergodic $n \Sigma v$ estimate. However, we find that in most astrophysical systems, deviations from a perfectly Keplerian or harmonic potential (due to e.g. granularity or extended mass) trigger sufficient orbital precession to recover the $n \Sigma v$ interaction rate. Astrophysically relevant failures of the $n \Sigma v$ approach only seem to occur for tightly bound stars orbiting intermediate-mass black holes, or for the high-mass end of collisional cascades in certain debris disks.

To test the scenario that outflows accelerated by active galactic nuclei (AGN) have a major impact on galaxy-wide scales, we have analysed deep VLT/MUSE data for the type-2 quasar/ultraluminous infrared galaxy F13451+1232 - an object that represents the major mergers considered in models of galaxy evolution. After carefully accounting for the effects of atmospheric seeing that had smeared the emission from known compact nuclear outflows across the MUSE field of view, we find that the large-scale kinematics in F13451+1232 are consistent with gravitational motions that are expected in a galaxy merger. Therefore, the fast ($\mathrm{W_{80}}>500$ km s$^{-1}$) warm-ionised AGN-driven outflows in this object are limited to the central $\sim$100 pc of the galaxy, although we cannot rule out larger-scale, lower-velocity outflows. Moreover, we directly demonstrate that failure to account for the beam-smearing effects of atmospheric seeing would have led to the mass outflow rates and kinetic powers of spatially-extended emission being overestimated by orders of magnitude. We also show that beam-smeared compact-outflow emission can be significant beyond radial distances of 3.5 arcseconds (more than eight times the radius of the seeing disk), and support the argument that some previous claims of large-scale outflows in active galaxies were likely the result of this effect rather than genuine galaxy-wide ($r>5$ kpc) outflows. Our study therefore provides further evidence that warm-ionised AGN-driven outflows are limited to the central kiloparsecs of galaxies and highlights the critical importance of accounting for atmospheric seeing in ground-based observational studies of active galaxies.

Planet formation in the solar system was started when the first planetesimals were formed from the gravitational collapse of pebble clouds. Numerical simulations of this process, especially in the framework of streaming instability, produce various power laws for the initial mass function for planetesimals. While recent advances have shed light on turbulence and its role in particle clustering, a comprehensive theoretical framework linking turbulence characteristics to particle cluster properties and planetesimal mass function remains incomplete. Recently, a kinetic field theory for ensembles of point-like classical particles in or out of equilibrium has been applied to cosmic structure formation. This theory encodes the dynamics of a classical particle ensemble by a generating functional specified by the initial probability distribution of particles in phase space and their equations of motion. Here, we apply kinetic field theory to planetesimal formation. A model for the initial probability distribution of dust particles in phase space is obtained from a quasi-initial state for a three-dimensional streaming-instability simulation that is a particle distribution with velocities for gas and particles from the Nakagawa relations. The equations of motion are chosen for the simplest case of freely streaming particles. We calculate the non-linearly evolved density power spectrum of dust particles and find that it develops a universal $k^{-3}$ tail at small scales, suggesting scale-invariant structure formation below a characteristic and time-dependent length scale. Thus, the KFT analysis indicates that the initial state for streaming instability simulations does not impose a constraint on structure evolution during planetesimal formation.

There are hosts of surveys that will provide excellent data to search for and locate Fast Radio Bursts (FRBs) at cosmological distances. The BINGO project is one such surveys, and this collaboration has already estimated a FRB detection rate that the project will yield with the main telescope helped by a set of outrigger stations. This paper aims to simulate and estimate the potential of FRBs in constraining our current cosmological model. We present a forecast of the future constraints that the BINGO FRB detections and localizations will have when added to other current cosmological datasets. We quantify the dispersion measure (DM) as a function of redshift ($z$) for the BINGO FRB mock sample. Furthermore, we use current datasets (Supernovae, Baryonic Acoustic Oscillations, and Cosmic Microwave Background data) prior to assessing the efficacy of constraining dark energy models using Monte Carlo methods. Our results show that spatially located BINGO FRB dataset will provide promising constraints on the population of host galaxies intrinsic DM and be able to measure the nuisance parameters present within a FRB cosmological analysis. Still, they will also provide alternative estimates on other parameters such as the Hubble constant and the dark energy equation of state. In particular, we should see that BINGO FRB data can put constraints on the degenerate $w-H_0$ plane, which the CMB is incapable of measuring, allowing FRBs to be a viable alternative to BAO to constrain the dark energy equation of state. We conclude that FRBs remain a promising future probe for cosmology and that the FRBs detected by the BINGO project will contribute significantly to our knowledge of the current cosmological model.

The results of the analysis of 205 brightest sources ( $S>15$ mJy), which were found in the sky survey at the declination of the pulsar in the Crab Nebula, are presented. The survey was conducted at a frequency of 4.7~GHz using a three-beam radiometer complex installed in the focus of the Western Sector of the RATAN-600 radio telescope in 2018-2019. Based on the measurements and data collected in the database of astrophysical catalogs CATS built radio spectra of objects. For a quarter of all detected sources, data at a frequency higher than 4~GHz were obtained for the first time, and for the rest, they were appended. The variability of radiation sources on the scales of the year, from days to months, was studied. The greatest change in the radio flux was found in the blazar B2~1324+22. The search for daily variability was carried out for 26 the brightest sources with an average value of $S_{4.7} \sim 250$ mJy. All sources are identified with objects from optical and infrared catalogs. Radio luminosity is calculated for 112 objects with a known redshift.

The dramatic dimming episode of the red supergiant Betelgeuse in 2019/2020, caused by a partial darkening of the stellar disk, has highlighted gaps in the understanding of the evolution of massive stars. We analyzed numerical models to investigate the processes behind the formation of dark surface patches and the associated reduction in the disk-integrated stellar light. With the CO5BOLD code, we performed global 3D radiation-hydrodynamical simulations of evolved stars, including convection in the stellar interior, self-excited pulsations, and the resulting atmospheric dynamics with strong radiative shocks. We attribute dimming phenomena to obscuring clouds of cool gas in the lower atmosphere, forming according to three different scenarios. One process transports material outward in a strong shock, similar to what occurs in 1D simulations of radially pulsating AGB stars. Another mechanism is triggered by a large convective upflow structure, in combination with exceptionally strong radial pulsations. This induces Rayleigh-Taylor instabilities, causing plumes of material to be sent outward into the atmosphere. The third and rarest scenario involves large-amplitude convective fluctuations, leading to enhanced flows in deep downdrafts, which rebound and send material outward. In all cases, the dense gas above the stellar surface cools and darkens rapidly in visible light. AGB stars show localized dark patches regularly during intermediate phases of their large-amplitude pulsations, while more massive stars will only intermittently form such patches during luminosity minima. The episodic levitation of dense gas clumps above the stellar surface, followed by the formation of complex molecules in the cooling gas and possibly dust grains at a later stage, can account for the dark patches and strong dimming events of supergiant stars such as Betelgeuse.

The helium abundances in the multiple populations that are now known to comprise all closely studied Milky Way globular clusters are often inferred by fitting isochrones generated from stellar evolutionary models to globular cluster photometry. It is therefore important to build stellar models that are chemically self-consistent in terms of their structure, atmosphere, and opacity. In this work we present the first chemically self-consistent stellar models of the Milky Way globular cluster NGC 2808 using MARCS model atmospheres, OPLIB high-temperature radiative opacities, and AESOPUS low-temperature radiative opacities. These stellar models were fit to the NGC 2808 photometry using Fidanka , a new software tool that was developed to optimally fit cluster photometry to isochrones and for population synthesis. Fidanka can determine, in a relatively unbiased way, the ideal number of distinct populations that exist within a dataset and then fit isochrones to each population. We achieve this outcome through a combination of Bayesian Gaussian Mixture Modeling and a novel number density estimation algorithm. Using Fidanka and F275W-F814W photometry from the Hubble UV Globular Cluster Survey we find that the helium abundance of the second generation of stars in NGC 2808 is higher than the first generation by $15\pm3\%$. This is in agreement with previous studies of NGC 2808. This work, along with previous work by Dotter et al. (2015) focused on NGC 6752, demonstrates that chemically self-consistent models of globular clusters do not significantly alter inferred helium abundances, and are therefore unlikely to be worth the significant additional time investment.

Estimating the true background in an astronomical image is fundamental to detecting faint sources. In a typical low-photon count astronomical image, such as in the far and near-ultraviolet wavelength range, conventional methods relying on the 3-sigma clipping and median or mode estimation often fail to capture the true background level accurately. As a consequence, differentiating true sources from noise peaks remains a challenging task. Additionally, in such images, effectively identifying and excluding faint sources during the background estimation process remains crucial, as undetected faint sources could contaminate the background. This results in overestimating the true background and obscuring the detection of very faint sources. To tackle this problem, we introduce a geometric approach based on the method of steepest descent to identify local minima in an astronomical image. The proposed algorithm based on the minima statistics effectively reduces the confusion between sources and background in the image; thereby ensuring a better background estimation and enhancing the reliability of faint source detection. Our algorithm performs well compared to conventional methods in estimating the background even in crowded field images. In low-photon count, less crowded images, our algorithm recovers the background within 10\%, while traditional methods drastically underestimate it by a few orders of magnitude. In crowded fields, the conventional methods overestimates the background by $\sim 200\%$ whereas our algorithm recovers the true background within $\sim 14\%$. We provide a simple prescription to create a background map using our algorithm and discuss its application in large astronomical surveys.

The axion were proposed as a result to a solution to the Strong CP Problem in quantum chromodynamics (QCD) and is now considered a leading candidate for dark matter. Direct axion dark matter detection experiments are challenging due to the axion's weak interaction with electromagnetism. Recent work has suggested the possibility of an enhancement of astrophysical axion-to-photon decay through parametric resonance. We explore here the feasibility of using parametric resonance to enhance the signal in direct axion-like particle dark matter detectors.

We analyse the large-scale clustering of the Luminous Red Galaxy (LRG) and Quasar (QSO) sample from the first data release (DR1) of the Dark Energy Spectroscopic Instrument (DESI). In particular, we constrain the primordial non-Gaussianity (PNG) parameter $f_{\rm NL}^{\rm loc}$ via the large-scale scale-dependent bias in the power spectrum using $1,631,716$ LRGs ($0.6 < z < 1.1$) and $1,189,129$ QSOs ($0.8 < z < 3.1$). This new measurement takes advantage of the enormous statistical power at large scales of DESI DR1 data, surpassing the latest data release (DR16) of the extended Baryon Oscillation Spectroscopic Survey (eBOSS). For the first time in this kind of analysis, we use a blinding procedure to mitigate the risk of confirmation bias in our results. We improve the model of the radial integral constraint proposing an innovative correction of the window function. We also carefully test the mitigation of the dependence of the target selection on the photometry qualities by incorporating an angular integral constraint contribution to the window function, and validate our methodology with the blinded data. Finally, combining the two samples, we measure $f_{\rm NL}^{\rm loc} = {-3.6}_{-9.1}^{+9.0}$ at $68\%$ confidence, where we assume the universality relation for the LRG sample and a recent merger model for the QSO sample about the response of bias to primordial non-Gaussianity. Adopting the universality relation for the PNG bias in the QSO analysis leads to $f_{\rm NL}^{\rm loc} = 3.5_{-7.4}^{+10.7}$ at $68\%$ confidence. This measurement is the most precise determination of primordial non-Gaussianity using large-scale structure to date, surpassing the latest result from eBOSS by a factor of $2.3$.

High-energy neutrinos from the blazar TXS~0506+056 are usually thought to arise from the relativistic jet pointing to us. However, the composition of jets of active galactic nuclei (AGNs), whether they are baryon dominated or Poynting flux dominated, is largely unknown. In the latter case, no comic rays and neutrinos are expected from the AGN jets. In this work, we study whether the neutrino emission from TXS~0506+056 could be powered by the accretion flow of the supermassive black hole. Protons could be accelerated by magnetic reconnection or turbulence in the inner accretion flow. To explain the neutrino flare of TXS~0506+056 in the year of 2014-2015, a super-Eddington accretion is needed. During the steady state, a sub-Eddington accretion flow could power a steady neutrino emission that may explain the long-term steady neutrino flux from TXS 0506+056. We consider the neutrino production in both magnetically arrested accretion (MAD) flow and the standard and normal evolution (SANE) regime of accretion. In the MAD scenario, due to a high magnetic field, a large dissipation radius is required to avoid the cooling of protons due to the synchrotron emission.

When an exoplanet passes in front of its host star, the resulting eclipse causes an observable decrease in stellar flux, and when multiple such transits are detected, the orbital period of the exoplanet can be determined. Over the past six years, NASA's Transiting Exoplanet Survey Satellite (TESS) has discovered thousands of potential planets by this method, mostly with short orbital periods, although some have longer reported values over one hundred days. These long orbital periods, however, are note easy to confirm due to frequent lengthy data gaps. Here we show that while the majority of these long period candidates likely have periods much shorter than reported, there are a sizable number of TESS candidates with true long periods. These candidates generally only have two reported transits, but the periods of duo-transits like this, and even candidates with three or more transits, can be confirmed if the data rules out all possible shorter period aliases. Using TESS data, we confirm long orbital periods for nine candidate planets, and present five others that are likely long period. Due to their long periods, these planets will have relatively cool equilibrium temperatures, and may be more likely to host exomoons or rings. We present these TOIs, along with a variety of small corrections to other TESS orbital periods and three planet candidates with possible transit timing variations, with the goal of refining the TESS data set and enabling future research with respect to cool transiting planets.

Full-disk measurements of the solar magnetic field by the Helioseismic and Magnetic Imager (HMI) are often used for magnetic field extrapolations, but its limited spatial and spectral resolution can lead to significant errors. We compare HMI data with observations of NOAA 12104 by the Hinode Spectropolarimeter (SP) to derive a scaling curve for the magnetic field strength, $B$. The SP data in the \ion{Fe}{i} lines at 630\,nm were inverted with the SIR code. We find that the Milne-Eddington inversion of HMI underestimates $B$ and the line-of-sight flux, $\Phi$, in all granulation surroundings by an average factor of 4.5 in plage and 9.2 in the quiet Sun in comparison to the SP. The deviation is inversely proportional to the magnetic fill factor, $f$, in the SP results. We derived a correction curve to match the HMI $B$ with the effective flux $B\,f$ in the SP data that scaled HMI $B$ up by 1.3 on average. A comparison of non-force-free field extrapolations over a larger field of view without and with the correction revealed minor changes in connectivity and a proportional scaling of electric currents and Lorentz force ($\propto B \sim 1.3$) and free energy ($\propto B^2\sim 2$). Magnetic field extrapolations of HMI vector data with large areas of plage and quiet Sun will underestimate the photospheric magnetic field strength by a factor of 5--10 and the coronal magnetic flux by at least 2. An HMI inversion including a fill factor would mitigate the problem.

Nova super-remnants (NSRs) are substantially extended structures (>100 parsecs) encompassing recurrent novae. NSRs grow as a result of frequent nova eruptions transporting vast quantities of the locally surrounding interstellar medium away from the binary system over many millennia into a thin high-density shell, as the central white dwarf grows towards the Chandrasekhar limit. The prototypical NSR, first identified as such in 2014, is situated in the Andromeda Galaxy and belongs to the annually erupting nova, M31N 2008-12a (or '12a'). In this short review, modelling of evolving NSRs (including the 12a NSR) will be outlined as motivation towards searching for more of these phenomena in the Galaxy and beyond. The latest developments in this upcoming subfield of nova research will then be presented including the discovery of two new Galactic nova super-remnants (and their consequent modelling) and the first survey undertaken with the sole purpose of finding NSRs in the Andromeda Galaxy and the Large Magellanic Cloud.

\texttt{PSpectCosmo} is a high-performance \texttt{C++} program developed to investigate early-universe cosmological dynamics, with a specific emphasis on the inflationary epoch. Utilizing a Fourier-space pseudo-spectral method, \texttt{PSpectCosmo} enables the precise evolution of interacting scalar fields and gravitational waves, ensuring accurate representation of the power spectrum during inflation. This approach overcomes key limitations of finite difference methods, particularly in maintaining consistency between effective and lattice wave vectors. The code employs the adaptive step size velocity-Verlet algorithm for time integration, offering a balance of numerical stability and high precision. Additionally, \texttt{PSpectCosmo} incorporates a robust mechanism to compute convergent energy density, effectively resolving the issue of divergent energy density at the onset of inflation. These capabilities establish \texttt{PSpectCosmo} as a reliable and versatile tool for probing non-linear cosmological phenomena and inflationary dynamics with exceptional accuracy. The code is publicly available at \href{this https URL}{this https URL}.

The relationship between B-field orientation and density structure in molecular clouds is often assessed using the Histogram of Relative Orientations (HRO). We perform a plane-of-the-sky geometrical analysis of projected B-fields, by interpreting HROs in dense, spheroidal, prestellar and protostellar cores. We use James Clerk Maxwell Telescope (JCMT) POL-2 850 $\mu$m polarisation maps and Herschel column density maps to study dense cores in the Ophiuchus molecular cloud complex. We construct two-dimensional core models, assuming Plummer column density profiles and modelling both linear and hourglass B-fields. We find high-aspect-ratio ellipsoidal cores produce strong HRO signals, as measured using the shape parameter $\xi$. Cores with linear fields oriented $< 45^{\circ}$ from their minor axis produce constant HROs with $-1 < \xi < 0$, indicating fields are preferentially parallel to column density gradients. Fields parallel to the core minor axis produce the most negative value of $\xi$. For low-aspect-ratio cores, $\xi \approx 0$ for linear fields. Hourglass fields produce a minimum in $\xi$ at intermediate densities in all cases, converging to the minor-axis-parallel linear field value at high and low column densities. We create HROs for six dense cores in Ophiuchus. $\rho$ Oph A and IRAS 16293 have high aspect ratios and preferentially negative HROs, consistent with moderately strong-field behaviour. $\rho$ Oph C, L1689A and L1689B have low aspect ratios, and $\xi \approx 0$. $\rho$ Oph B is too complex to be modelled using a simple spheroidal field geometry. We see no signature of hourglass fields, agreeing with previous findings that dense cores generally exhibit linear fields on these size scales.

There is an ever-growing need in the gravitational wave community for fast and reliable inference methods, accompanied by an informative error bar. Nested sampling satisfies the last two requirements, but its computational cost can become prohibitive when using the most accurate waveform models. In this paper, we demonstrate the acceleration of nested sampling using a technique called posterior repartitioning. This method leverages nested sampling's unique ability to separate prior and likelihood contributions at the algorithmic level. Specifically, we define a `repartitioned prior' informed by the posterior from a low-resolution run. To construct this repartitioned prior, we use a $\beta$-flow, a novel type of conditional normalizing flow designed to better learn deep tail probabilities. $\beta$-flows are trained on the entire nested sampling run and conditioned on an inverse temperature $\beta$. Applying our methods to simulated and real binary black hole mergers, we demonstrate how they can reduce the number of likelihood evaluations required for convergence by up to an order of magnitude, enabling faster model comparison and parameter estimation. Furthermore, we highlight the robustness of using $\beta$-flows over standard normalizing flows to accelerate nested sampling. Notably, $\beta$-flows successfully recover the same posteriors and evidences as traditional nested sampling, even in cases where standard normalizing flows fail.

The Hf-W isotopic system is the reference chronometer for determining the chronology of Earth's accretion and differentiation. However, its results depend strongly on uncertain parameters, including the extent of metal-silicate equilibration and the siderophility of tungsten. Here we show that a multistage core-formation model based on N-body accretion simulations, element mass balance and metal-silicate partitioning, largely eliminates these uncertainties. We modified the original model of Rubie et al. (2015) by including (1) smoothed particle hydrodynamics estimates of the depth of melting caused by giant impacts and (2) the isotopic evolution of 182W. We applied two metal-silicate fractionation mechanisms: one when the metal delivered by the cores of large impactors equilibrates with only a small fraction of the impact-induced magma pond and the other when metal delivered by small impactors emulsifies in global magma oceans before undergoing progressive segregation. The latter is crucial for fitting the W abundance and 182W anomaly of Earth's mantle. In addition, we show, for the first time, that the duration of magma ocean solidification has a major effect on Earth's tungsten isotope anomaly. We re-evaluate the six Grand Tack N-body simulations of Rubie et al. (2015). Only one reproduces epsilon182W=1.9+/-0.1 of Earth's mantle, otherwise accretion is either too fast or too slow. Depending on the characteristics of the giant impacts, results predict that the Moon formed either 143-183 Myr or 53-62 Myr after the start of the solar system. Thus, independent evaluations of the Moon's age provide an additional constraint on the validity of accretion simulations.

The dispersion measure (DM) of fast radio bursts (FRBs) is sensitive to the electron distribution in the Universe, making it a promising probe of cosmology and astrophysical processes such as baryonic feedback. However, cosmological analyses of FRBs require knowledge of the contribution to the observed DM coming from the FRB host. The size and distribution of this contribution is still uncertain, thus significantly limiting current cosmological FRB analyses. In this study, we extend the baryonification (BCM) approach to derive a physically-motivated, analytic model for predicting the host contribution to FRB DMs. By focusing on the statistical properties of FRB host DMs, we find that our simple model is able to reproduce the probability distribution function (PDF) of host halo DMs measured from the CAMELS suite of hydrodynamic simulations, as well as their mass- and redshift dependence. Furthermore, we demonstrate that our model allows for self-consistent predictions of the host DM PDF and the matter power spectrum suppression due to baryonic effects, as observed in these simulations, making it promising for modelling host-DM-related systematics in FRB analyses. In general, we find that the shape of the host DM PDF is determined by the interplay between the FRB and gas distributions in halos. Our findings indicate that more compact FRB profiles require shallower gas profiles (and vice versa) in order to match the observed DM distributions in hydrodynamic simulations. Furthermore, the analytic model presented here shows that the shape of the host DM PDF is highly sensitive to the parameters of the BCM. This suggests that this observable could be used as an interesting test bed for baryonic processes, complementing other probes due to its sensitivity to feedback on galactic scales. We further discuss the main limitations of our analysis, and point out potential avenues for future work.

We search for radio emission from star-planet interactions in the M-dwarf system GJ~486, which hosts an Earth-like planet. We observed the GJ~486 system with the upgraded Giant Metrewave Radio Telescope (uGMRT) from 550 to 750 MHz in nine different epochs, between October 2021 and February 2022, covering almost all orbital phases of GJ~486 b from different orbital cycles. We obtained radio images and dynamic spectra of the total and circularly polarized intensity for each individual epoch We do not detect any quiescent radio emission in any epoch above 3$\sigma$. Similarly, we do not detect any bursty emission in our dynamic spectra. While we cannot completely rule out that the absence of a radio detection is due to time variability of the radio emission, or to the maximum electron-cyclotron maser emission being below our observing range, this seems unlikely. We discuss two possible scenarios: an intrinsic dim radio signal, or alternatively, that the anisotropic beamed emission pointed away from the observer. If the non-detection of radio emission from star-planet interaction in GJ~486 is due to an intrinsically dim signal, this implies that, independently of whether the planet is magnetized or not, the mass-loss rate is small (\dot{M}_\star $\lesssim$ 0.3 \dot{M}_\sun) and that, concomitantly, the efficiency of the conversion of Poynting flux into radio emission must be low ($\beta \lesssim 10^{-3}$). Free-free absorption effects are negligible, given the high value of the coronal temperature. Finally, if the anisotropic beaming pointed away from us, this would imply that GJ~486 has very low values of its magnetic obliquity and inclination.

Certain holographic dark energy (HDE) models allow for the possibility of a ``long freeze,'' in which the scale factor evolves to a constant in the long-time limit. Here we extend previous calculations by adding a nonrelativistic matter component. The addition of a matter component tends to destroy the long freeze behavior, driving the universe to recollapse. Long freeze evolution is still possible, but only for a limited set of HDE models.

Solitons in relativistic field theories are not necessarily topologically charged. In particular, non-topological solitons -- known as Q-balls -- arise naturally in nonlinear field theories endowed with attractive interactions and internal symmetries. Even without stabilizing internal symmetries, quasi-solitons known as oscillons, which are long-lived, can also exist. Both Q-balls and oscillons have significant applications in cosmology and particle physics. This review is an updated account of the intriguing properties and dynamics of these non-topological solitons and quasi-solitons, as well as their important roles in early-universe scenarios and particle physics models.

We present a Python module for simulating Silicon Photo-Multipliers, Avalanche Photo-Diodes, and Multi-Pixel Photon Counters. This module allows users to perform noise analyses: Dark Count Rate, crosstalk, and afterpulsing. Furthermore, the simulation framework novelty is the capability of simulating assemblies of SiPM arrays (MPPCa) for large area detectors like Ring Imaging Cherenkov detectors, Cherenkov Telescopes, Positron Emission Tomography, and any detector using SiPM arrays. Users can simulate ring- or shower-like-shaped signals based on the expected number of photons generated by the source. We validate the performance of the simulation module with data from four different SiPM: Broadcom AFBR-S4N66P024M, Hamamatsu S14160-636050HS, Onsemi MICROFC-60035, and FBK NUV-HD3.

Viscous fluids can dissipate and alter the propagation of gravitational waves, as well as modify the relaxation and stability properties of self-gravitating fluids. This is particularly relevant in order to understand the relaxation to equilibrium of neutron stars, and their gravitational wave emission. Here we study the linearized theory of perturbations of spherically symmetric self-gravitating fluids. Dissipative effects are included through the hydrodynamics theory of Bemfica, Disconzi, Noronha, and Kovtun (BDNK). This theory has been shown to be causal and stable, despite involving only first order gradients. We show how the problem reduces to two coupled wave equations in the axial sector, one of them associated to a novel viscous mode, and including explicitly dissipative terms. In the polar sector, we reduce the problem to five coupled wave equations and one additional constraint. We comment on their causal structure, and recover the causality bounds of the BDNK theory.

We study the dynamical response of viscous materials to gravitational waves, in the context of a fully relativistic theory of fluid dynamics. For the first time, we calculate oscillation modes and scattering properties of viscous stars. Viscous stars absorb high frequency radiation, following a dispersion relation introduced by Press. In the extremely large viscosity regime, stars would become reflectors of waves, but this regime appears to be forbidden by causality bounds. In the context of black hole mimickers, we show how maximally viscous stars on the threshold of stability mimic the absorption of a black hole with the same mass. Our results suggest that rotating viscous stars will amplify incoming radiation.

In a cosmological context, the Einstein-Gauss-Bonnet theory contains, in $d+4$ dimensions, a dynamical compactification scenario in which the additional dimensions settle down to a configuration with a constant radion/scale factor. Sadly however this work demonstrates that such a quite appealing framework is plagued by instabilities, either from the background configuration's unsteadiness or the ghostly behaviors of the tensorial perturbations. New and stable solutions are found by relaxing one of the hypotheses defining the original compactification scenario. However, such configurations do not respect the current bounds on the speed of propagation of gravitational waves, and thus have to be discarded. Those results thus advocate for a comprehensive study of compactification scenarios in the Gauss-Bonnet framework, their stability, and the effects of matter inclusion.

The final stage of black hole evaporation is a potent probe of physics beyond the Standard Model: Hawking-Bekenstein radiation may be affected by quantum gravity "memory burden effects", or by the presence of "dark", beyond-the-Standard-Model degrees of freedom in ways that are testable with high-energy gamma-ray observations. We argue that information on either scenario can best be inferred from measurements of the evaporation's lightcurve and by correlating observations at complementary energies. We offer several new analytical insights in how such observations map on the fundamental properties of the evaporating black holes and of the possible exotic particles they can evaporate into.

Ultralight bosons, proposed as candidates for dark matter, are predicted by various new physics models. In the presence of bosons with suitable masses, superradiant (SR) instability can naturally transform a spinning black hole (BH) into a gravitational atom (GA). Here we study the dynamics of intermediate mass-ratio inspirals (IMRIs) around a GA formed by ultralight vector field saturated in its SR ground state. We employ a perturbative model at the leading Newtonian order to consistently account for both the conservative effect of cloud gravity and the dissipative effect of cloud ionization. We find the cloud can make a sizable negative contribution to the secular periastron precession at binary separations comparable to the gravitational Bohr radius. Meanwhile, the backreaction of ionization could significantly accelerate the process of orbital decay and circularization. Considering reasonably small vector boson masses, we investigate the adiabatic orbital evolution and gravitational waveforms of eccentric inspirals. The results indicate that vector GAs might be detectable through observations of low-frequency IMRIs by future space-based gravitational-wave detectors, such as LISA and Taiji.

Using the Skyrme model Sk$\chi$450 constrained by the chiral effective field theory and the ground-state energies of doubly-magic nuclei, we explore the macroscopic static energy spectrum of dense matter. Structure of the matter is idealized as 1-, 2- or 3-dimensional periodic Coulomb lattices (pasta phases) at average baryon density in the range of $0.005n_0-0.5n_0$ with $n_0\sim0.16$ fm$^{-3}$, corresponding to the inner crust of neutron stars. In the earlier work, the bulk and the surface nuclear properties in this scenario were described on the basis of different nuclear interactions. As a result, predictions for the inner crust structure based on those numerical results are not necessarily self-consistent. In this work, we solve this problem by describing the nuclear properties and proton-proton Cooper pairing in a unified manner, \emph{i.e.} based on the same nuclear interaction. We calculate the surface tension to the leading order (planar surface) and to the next-to-leading order (from nonzero principal curvature) at the Extended Thomas-Fermi level. Next, we develop and solve a system of equations within the compressible liquid drop model (LDM), which provides us with all necessary information about the macroscopic static energy spectrum of the nuclear structure. We find that the curvature corrections change the ground state in a relevant way. We also find that, typically, the energy differences between the different pasta phases are less than the thermal energy for temperatures $\sim10^8-10^9$ K, which implies that the real nuclear pasta phases are likely polymorphic. Finally, for 1-dimensional pasta phase we evaluate the superconducting coherence length of protons, the London penetration depth and the superconducting energy gap. Our results offer a preliminary insight into rich magnetic properties of the pasta phases.

We show that, in the framework of quantum field theory in curved spacetime, the semiclassical energy-momentum tensor of the neutrino flavor vacuum fulfills the equation of state of dust and cold dark matter. We consider spherically symmetric spacetimes, and we demonstrate that, within the weak field approximation, the flavor vacuum contributes as a Yukawa correction to the Newtonian potential. This corrected potential may account for the flat rotation curves of spiral galaxies. In this perspective, neutrino mixing could contribute to dark matter

The gauge singlet right-handed neutrinos are one of the essential fields in neutrino mass models that explain tiny masses of active neutrinos. We consider the effective field theory of the Standard Model extended with these fields under the assumption that neutrinos are Dirac particles. In this framework, we provide a comprehensive study for the phenomenological consequences of various dimension six interactions employing various high and low energy observables such as the $pp\to\ell^\pm+\not{E}_T$, $pp\to j+\not{E}_T$, decays of proton, meson, tau and top, as well as the cosmological parameter $N_{\rm eff}$.

The proton-rich nucleus $^{22}$Si is studied using Nuclear Lattice Effective Field Theory with high-fidelity chiral forces. Our results indicate that $^{22}$Si is more tightly bound than $^{20}$Mg, thereby excluding the possibility of two-proton emission. The $Z = 14$ shell closure in $^{22}$Si is supported by the evolution of the $2^+$ state in the neighboring nuclei. We then focus on the charge radius and spatial distribution information of $^{22}$Si, considering the novel phenomena that may emerge due to the small two-proton separation energy and the shell closure. We present the distribution of the $14$ protons and $8$ neutrons obtained from our lattice simulation, revealing insights into the spatial arrangement of the nucleons. Moreover, the spatial localization of the outermost proton and neutron suggests that $^{22}$Si is a doubly magic nucleus. Furthermore, we develop the pinhole method based on the harmonic oscillator basis, which gives insight into the nuclear structure in terms of the shell model picture from lattice simulations. Our calculated occupation numbers support that $Z = 14$ and $N = 8$ are the shell closures and show that the $\pi 1s_{1/2}$ orbital component is minor in $^{22}$Si.

The interaction between a thin foil target and a circularly polarized laser light injected along an external magnetic field is investigated numerically by particle-in-cell simulations. A standing wave appears at the front surface of the target, overlapping the injected and partially reflected waves. Hot electrons are efficiently generated at the standing wave due to the relativistic two-wave resonant acceleration if the magnetic field amplitude of the standing wave is larger than the ambient field. A bifurcation occurs in the gyration motion of electrons, allowing all electrons with non-relativistic velocities to acquire relativistic energy through the cyclotron resonance. The optimal conditions for the highest energy and the most significant fraction of hot electrons are derived precisely through a simple analysis of test-particle trajectories in the standing wave. Since the number of hot electrons increases drastically by many orders of magnitude compared to the conventional unmagnetized cases, this acceleration could be a great advantage in laser-driven ion acceleration and its applications.