Papers for Thursday, Jan 30 2025

This paper presents a design for a drop tower test to evaluate a numerical model for a structurally reconfigurable spacecraft with actuatable joints, referred to as a transformable spacecraft. A mock-up robot for a 3U-sized transformable spacecraft is designed to fit in a limited time and space of the microgravity environment available in the drop tower. The robot performs agile reorientation, referred to as nonholonomic attitude control, by actuating joints in a particular manner. To adapt to the very short duration of microgravity in the drop tower test, a successive joint actuation maneuver is optimized to maximize the amount of attitude reorientation within the time constraint. The robot records the angular velocity history of all four bodies, and the data is analyzed to evaluate the accuracy of the numerical model. We confirm that the constructed numerical model sufficiently replicates the robot's motion and show that the post-experiment model corrections further improve the accuracy of the numerical simulations. Finally, the difference between this drop tower test and the actual orbit demonstration is discussed to show the prospect.

In this report we summarize the activities of the International Astronomical Consortium for High Energy Calibration (IACHEC) from the 16th IACHEC Workshop at Parador de La Granja, Spain. Sixty-one scientists directly involved in the calibration of operational and future high-energy missions gathered during 3.5 days to discuss the status of the cross-calibration between the current international complement of X-ray observatories, and the possibilities to improve it. This summary consists of reports from the Working Groups with topics ranging across: the identification and characterization of standard calibration sources, multi-observatory cross-calibration campaigns, appropriate and new statistical techniques, calibration of instruments and characterization of background, preservation of knowledge, and results for the benefit of the astronomical community.

We examine the applicability of the initial-to-final mass relation (IFMR) for white dwarfs (WDs) in intermediate-separation binary systems (approximately 1 AU), using astrometric binaries identified in open clusters from Gaia DR3. A careful analysis of the astrometric orbits and spectral energy distributions isolates 33 main-sequence (MS) stars with highly likely WD companions. By combining cluster age estimates, dynamically measured WD masses, and, where available, WD cooling temperatures, we derive progenitor masses for 26 WD candidates. Our analysis suggests the presence of two distinct WD populations: (i) low-mass WDs, likely shaped by binary interactions during the progenitor's red-giant phase; and (ii) "spender" WDs, which experienced higher-than-expected mass loss and have progenitor masses above the IFMR predictions. The rest of the candidates, referred to as the "others", represent systems with inconclusive formation mechanisms. We suggest that at least some of these systems might be hierarchical triples, where the companion to the MS is a double WD or a double-WD merger product. However, follow-up studies are required to determine the nature of each case. These results highlight significant deviations from the IFMR derived for isolated WDs, emphasizing the role of binary evolution. Follow-up observations, particularly in the far-ultraviolet, are crucial for refining these findings and advancing our understanding of mass transfer processes and binary evolution pathways.

We report the results of a rapid follow-up campaign on the Type IIb Supernova (SN) 2022hnt. We present a daily, multi-band, photometric follow-up using the Las Cumbres Observatory, the Zwicky Transient Facility, the orbiting \textit{Swift} observatory, and the Asteroid Terrestrial-impact Last Alert System (ATLAS). A distinctive feature in the light curve of SN 2022hnt and other IIb SNe is an early narrow peak prior to the ${}^{56}$Ni peak caused by rapid shock cooling of the hydrogen envelope, which can serve as an important probe of the properties of the massive progenitor star in the moments before explosion. Using SN 2022hnt as a case study, we demonstrate a framework of considerations for the application of shock cooling models to type IIb SNe, outlining a consistent procedure for future surveys of Type IIb SNe progenitor and explosion properties. \hll{We fit several recent models of shock-cooling emission and obtain progenitor radii between $\sim50$ and $\sim100$ $R_\odot$, as well as hydrogen-enriched envelope masses between $\sim0.01$ and $\sim0.1$ $M_\odot$, both consistent with values for other IIb SNe. One of these models is the model of \cite{Morag2023}, marking the first time this model has been applied to a Type IIb SN.} We evaluate contrasting predictions between shock-cooling models to construct a fiducial parameter set which can be used for comparison to other SNe. Finally, we investigate the possibility of extended wind breakout or precursor emission captured in the earliest detections.

We present the third phase of the largest high-frequency, high-resolution imaging survey of 231 nearby, hard X-ray selected AGN, with a very high $98 \pm 1\%$ detection fraction. This survey presents VLA 22 GHz radio observations with 1" spatial resolution covering over $6$ orders of magnitude in radio luminosity in nearby AGN that span $\sim4$ orders of magnitude in black hole mass and X-ray luminosity. We identify three different radio morphologies: $44 \pm 3\%$ (102/231) are compact or unresolved, $46 \pm 3\%$ (106/231) show an extended structure (star formation, possible one-sided jets, etc.), and $8 \pm 2\%$ (19/231) have a biconical or two-sided jet-like morphology. The remaining $2 \pm 1\%$ (4/231) sources are non-detections. The radio-to-X-ray luminosity ratios of the Swift-BAT AGN ($\text{L}_R/\text{L}_{14-195 \text{keV}} \sim 10^{-5.5}$ and $\text{L}_R/\text{L}_{2-10 \text{keV}} \sim 10^{-5}$) with a scatter of $\sim0.5$ dex are similar to that of coronally active stars ($\text{L}_R/\text{L}_X \sim 10^{-5}$). For most targets, extended emission in radio-quiet objects is broadly consistent with the expectation for star formation from previous FIR observations, once the contribution from the radio core has been subtracted. Our sample represents nearby analogs of distant AGN at the peak of black hole growth, and thus the high detection fraction in our work has important implications for future high frequency AGN radio surveys with the next generation VLA (ngVLA) or Square Kilometre Array (SKA), both of which should detect large fractions of more distant AGN.

Context: Tidal tails of open clusters are the result of stellar evaporation from the cluster through the Galactic potential and internal dynamics. With the recent availability of high-precision data, tidal tails are being detected for most of the nearby open clusters. Aims: We identify the tidal tail members for all open clusters within a distance of 400 pc that are older than 100 Myr and have $>$100 members. To do this, we use model-independent methods. Methods: We used the convergent-point (CP) method to identify the co-moving stars near the open clusters using Gaia DR3 data. A new method called the self-compact convergent-point (SCCP) method was proposed and applied to some of the clusters. It performed better overall in tracing the tails. We also analysed the colour-magnitude diagrams and orbital energy to diagnose possible contamination. Results: Nineteen out of 21 clusters have tidal tails. Five of them were discovered for the first time through this work. The typical span of the tidal tails is 20--200 pc, and 30--700 member stars lie in the region inside the tidal radius and the tidal tails. Four out of 19 tidal tails are tilted away from direction of the Galactic centre. This contradicts the known theory of the tidal-tail formation. The luminosity functions of the tails and clusters are consistent with each other and with the canonical stellar interstellar mass function, but systematically higher radial velocities for the trailing tail than for the leading tail were observed for the first time. Conclusions: The CP method is useful for detecting tidal tails on a scale of approx. 100 pc for clusters closer than 400 pc. A further analysis of theoretical N-body models is required to understand the incompleteness and biases in the current sample of tidal tails.

The assembly of supermassive black hole (SMBH) mass ($M_{\bullet}$) and stellar mass ($M_{*}$) in galaxies can be studied via the redshift evolution of the $M_{\bullet}-M_{*}$ relation, but the ways in which selection bias and physical assembly channels affect this evolution are uncertain. To address this, we compare the $M_{\bullet}-M_{*}$ relation for local massive ($M_{*}>10^{10.5}$M$_{\odot}$) quiescent early-type galaxies (ETGs) to that for massive ETGs hosting active galactic nuclei (AGN) at $z\sim0.8$. The restrictions on stellar mass and galaxy type limit the assembly channels that may connect the two relations. For the local sample we find $\log(M_{\bullet}) = 8.80 + 1.10(\log{M_{*}-11})$, in line with prior work. For the $z\sim0.8$ sample we find a bias-corrected relation: $\log(M_{\bullet}) = 7.80 + 1.25(\log{M_{*}-11})$. We show, however, that this relation depends on the stellar and SMBH mass functions used to compute the selection bias, the virial relation, the virial factor, and the active fraction, which together introduce uncertainty of up to $\sim0.6$\,dex in the $z\sim0.8$ relation. Adopting reasonable choices of these parameters then our $z\sim0.8$ relation lies above that for $z\sim0$ AGN by $\sim0.5$\,dex, but below our $z\sim0$ ETG relation by $0.4-1$\,dex in SMBH mass. We discuss possible sources of this offset, including further bias corrections, `downsizing" in SMBH mass assembly, and preferential SMBH growth. Our results highlight the need to reduce uncertainties from selection and measurement bias in SMBH and stellar masses at all redshifts.

Hydrodynamical simulations of protoplanetary disk dynamics are useful tools for understanding the formation of planetary systems, including our own. Approximations are necessary to make these simulations computationally tractable. A common assumption when simulating dust fluids is that of a constant Stokes number, a dimensionless number that characterizes the interaction between a particle and the surrounding gas. Constant Stokes number is not a good approximation in regions of the disk where the gas density changes significantly, such as near a planet-induced gap. In this paper, we relax the assumption of constant Stokes number in the popular FARGO3D code using semi-analytic equations for the drag force on dust particles, which enables an assumption of constant particle size instead. We explore the effect this change has on disk morphology and particle fluxes across the gap for both outward- and inward-drifting particles. The assumption of constant particle size, rather than constant Stokes number, is shown to make a significant difference in some cases, emphasizing the importance of the more accurate treatment.

Combining multiple events into population analyses is a cornerstone of gravitational-wave astronomy. A critical component of such studies is the assumed population model, which can range from astrophysically motivated functional forms to non-parametric treatments that are flexible but difficult to interpret. In practice, the current approach is to fit the data multiple times with different population models to identify robust features. We propose an alternative strategy: assuming the data have already been fit with a flexible model, we present a practical recipe to reconstruct the population distribution of a different model. As our procedure postprocesses existing results, it avoids the need to access the underlying gravitational-wave data again and handle selection effects. Additionally, our reconstruction metric provides a goodness-of-fit measure to compare multiple models. We apply this method to the mass distribution of black-hole binaries detected by LIGO/Virgo/KAGRA. Our work paves the way for streamlined gravitational-wave population analyses by fitting the data once and for all with advanced non-parametric methods and careful handling of selection effects, while the astrophysical interpretation is then made accessible using our reconstruction procedure on targeted models. The key principle is that of conceptually separating data description from data interpretation.

In this study, we analyzed the close binary open clusters CWNU 2666 and HSC 224, which are in close spatial proximity, using photometric and astrometric data from the {\it Gaia} DR3 catalog. Likely member stars were identified based on a membership probability threshold ($P \geq 0.5$), resulting in 106 and 146 members for CWNU 2666 and HSC 224, respectively. The mean proper motion components ($\mu_{\alpha}\cos\delta$, $\mu_{\delta}$) were determined to be (0.646$\pm$0.155, -0.769$\pm$0.124) mas yr$^{-1}$ for CWNU 2666, and (0.665$\pm$0.131, -0.728$\pm$0.107) mas yr$^{-1}$ for HSC 224. The isochrone distances ($d_{\rm iso}$) were estimated as 1885$\pm$44 pc for CWNU 2666 and 1866$\pm$29 pc for HSC 224. The corresponding cluster ages ($t$) were derived as 160$\pm$15 Myr and 140$\pm$15 Myr, respectively. The astrometric and fundamental astrophysical parameters derived in this study demonstrate that the two open clusters are a close pair of open clusters.

Despite the fundamental importance of the star formation history (SFH) and the initial mass function (IMF) in the description of the Milky Way, their consistent and robust derivation is still elusive. Recent and accurate astrometry and photometry collected by the Gaia satellite provide the natural framework to consolidate these ingredients in our local Galactic environment. We aim to simultaneously infer the IMF and the SFH of the Galactic disc comparing Gaia data with the mock catalog resulting from the Besançon population synthesis model (BGM). Our goal is also to estimate the impact of the systematics present in current stellar evolutionary models on this inference. We use a new implementation of the BGM Fast Approximate Simulations (BGM FASt) framework to fit the seven million star Gaia DR3 all-sky $G<13$ color-magnitude diagram (CMD) to the most updated dynamically self-consistent BGM. Our derived SFH supports an abrupt decrease of the star formation approximately 1-1.5 Gyr ago followed by a significant enhancement with a wide plateau in the range 2-6 Gyr ago. A remarkable hiatus appears around 5-7 Gyr ago with a $\sim$1 Gyr shift depending on the set of stellar models. A complex and discrepant evolution at ages older than 8 Gyr deserves further investigation. Precise but discrepant values are found for the power-law indices of the IMF. For the range 0.5-1.53 $M_\odot$ the slope takes a value of $\alpha_2 = 1.45^{+0.19}_{-0.12}$, while for masses larger than 1.53 $M_\odot$ we obtain $\alpha_3 = 1.98^{+0.13}_{-0.05}$. The current implementation of the BGM FASt framework is ready to address executions fitting all-sky Gaia data up to 14-17 apparent limiting magnitude. This will naturally allow us to derive both a reliable SFH for the early epochs of the Galactic disc evolution and a precise slope for the IMF at low masses.

The detection of the extended X-ray-emission of the intracluster medium by the first SRG/eROSITA All-Sky Survey (eRASS1), combined with optical and near-infrared follow-up, resulted in the identification of more than 12000 galaxy clusters, yielding precise constraints on cosmological parameters. However, some clusters of galaxies can be misclassified as point sources by eROSITA's source detection algorithm due to the interplay between the point-spread function, the shallow depth of the survey, compact (cool core) X-ray emission, and bright active galactic nuclei hosted in their centers or their vicinity. To identify such misclassified galaxy clusters and groups, we apply optical follow-up to the eRASS1 X-ray point sources analogously to the treatment of the extent-selected catalog. After rigorous filtering to ensure purity, we find a total of 8347 clusters of galaxies, of which 5819 are novel detections, in a redshift range $0.05 < z \lesssim 1.1$. This corresponds to a 70 % discovery rate, a fraction similar to that of the extent-selected sample. To facilitate finding new exceptional clusters such as the Phoenix cluster (which is recovered in our sample), we divide the clusters into five classes based on the optical properties of likely single-source counterparts to the X-ray emission. We further investigate potential biases in our selection process by analyzing the optical and X-ray data. With this work, we provide a catalog of galaxy clusters and groups in the eRASS1 point source catalog, including their optical and X-ray properties along with a meaningful classification.

Short-period planets provide ideal laboratories for testing star-planet interaction. Planets that are smaller than $\sim$2$R_\oplus$ are considered to be largely rocky either having been stripped of or never having acquired the gaseous envelope. Zooming in on these short-period rocky planet population, clear edges appear in the mass-period and radius-period space. Over $\sim$0.2--20 days and 0.09--1.42$M_\odot$, the maximum mass of the rocky planets stay below $\sim$10$M_\oplus$ with a hint of decrease towards $\lesssim$1 day, $\gtrsim$4 day, and $\lesssim 0.45 M_\odot$. In radius-period space, there is a relative deficit of $\lesssim$2$R_\oplus$ planets inside $\sim$1 day. We demonstrate how the edges in the mass-period space can be explained by a combination of tidal decay and photoevaporation whereas the rocky planet desert in the radius-period space is a signature of magnetic drag on the planet as it orbits within the stellar magnetic field. Currently observed catastrophically evaporating planets may have started their death spiral from $\sim$1 day with planets of mass up to $\sim$0.3$M_\oplus$ under the magnetic drag. More discoveries and characterization of small planets around mid-late M and A stars would be welcome to better constrain the stellar parameters critical in shaping the edges of rocky planet population including their UV radiation history, tidal and magnetic properties.

We present nucleosynthesis and light-curve predictions for a new site of the rapid neutron capture process ($r$-process) from magnetar giant flares (GFs). Motivated by observations indicating baryon ejecta from GFs, Cehula et al. (2024) proposed mass ejection occurs after a shock is driven into the magnetar crust during the GF. We confirm using nuclear reaction network calculations that these ejecta synthesize moderate yields of third-peak $r$-process nuclei and more substantial yields of lighter $r$-nuclei, while leaving a sizable abundance of free neutrons in the outermost fastest expanding ejecta layers. The final $r$-process mass fraction and distribution are sensitive to the relative efficiencies of $\alpha$-capture and $n$-capture freeze-outs. We use our nucleosynthesis output in a semi-analytic model to predict the light curves of novae breves, the transients following GFs powered by radioactive decay. For a baryonic ejecta mass similar to that inferred of the 2004 Galactic GF from SGR 1806-20, we predict a peak UV/optical luminosity of $\sim 10^{39}$-$10^{40}\,\rm erg\,s^{-1}$ at $\sim 10$-$15$ minutes, rendering such events potentially detectable following a gamma-ray trigger by wide-field transient monitors such as ULTRASAT/UVEX to several Mpc. The peak luminosity and timescale of the transient increase with the GF strength due to the larger ejecta mass. Although GFs likely contribute 1-10% of the total Galactic $r$-process budget, their short delay-times relative to star-formation make them an attractive source to enrich the earliest generations of stars.

We examine the origin of molecular gas heating in a sample of 42 infrared-luminous galaxies at $z<0.3$ by combining two sets of archival data. First, integrated CO line luminosities in the 1-0 and 5-4 through 13-12 transitions. Second, results from radiative transfer modelling that decompose their bolometric emission into starburst, AGN, and host galaxy components. We find that the CO 1-0 and 5-4 through 9-8 lines primarily arise via radiative heating in the starburst and the host galaxy. In contrast, the CO 10-9 through 13-12 lines may arise primarily in the starburst and AGN, with an increasing contribution from mechanical heating and shocks. For the sample as a whole, we find no evidence that AGN luminosity affects the heating of molecular gas by star formation. However, for starbursts with low initial optical depths, a more luminous AGN may reduce the efficiency of starburst heating of the CO 5-4 and above lines, consistent with negative AGN feedback.

The proliferation of exoplanet discoveries in exotic environments like the Neptune desert challenges our understanding of planetary atmospheres under intense irradiation. The unexpected discovery of LTT9779 b, an ultra-hot Neptune within this desert, offers a prime opportunity for atmospheric studies. We build on prior observations of LTT9779 b from TESS, Spitzer, and CHEOPS, incorporating new VLT/ESPRESSO data to probe its atmospheric dynamics. Preliminary analyses suggest a metal-rich atmosphere and a high day-side geometric albedo, possibly indicating silicate clouds. Minimal atmospheric escape is observed, contrasting existing models of planetary evolution under extreme irradiation. We obtained the transmission spectrum of LTT9779 b between 0.4 and 0.78 microns with ESPRESSO, addressing systematics across three transits. Our analysis focused on the sodium doublet and H-alpha, using cross-correlation with models containing Na, K, FeH, TiO, and VO. No significant atmospheric signal was detected, with metallicity limits set at [Fe/H] >= 2.25 (>= 180 times solar). The non-detection aligns with a high-metallicity, cloud-free model, implying a high mean molecular weight and reduced atmospheric scale height. We interpret this as evidence for a metal-rich atmosphere with suppressed spectral features, possibly due to high-altitude clouds or hazes. These findings are consistent with JWST observations, supporting the hypothesis of metal-rich atmospheres obscured by aerosols in extreme environments.

Plasma velocity distribution functions (VDFs) constitute a fundamental observation of numerous operational and future missions. An efficient parameterization of VDFs is crucial for (1) preserving enough information to investigate macroscopic moments along with kinetic effects, (2) producing smooth distributions whereby it is possible to perform derivatives in phase space to support numerical solvers, and (3) economic data management and its storage. Previous studies have used spherical harmonics as an efficient basis for representing electron VDFs. In this paper, we present a novel algorithm targeted towards decomposing ion VDFs measured by electrostatic analyzers onboard Magnetospheric Multiscale Mission (MMS) and Solar Orbiter (SolO) spacecrafts. We use Slepian functions, custom-designed bases providing compact support in phase space, initially developed in information theory and later used for terrestrial and planetary applications. In this paper, we choose well-studied, well-measured, and complex intervals from MMS and SolO containing a range of simpler gyrotropic and agyrotropic distributions to benchmark the robustness of our reconstruction method. We demonstrate the advantages of using Slepian functions over spherical harmonics for solar wind plasma distributions. We also demonstrate that our choice of basis representation efficiently preserves phase space complexities of a 3D agyrotropic distribution function. This algorithm shown in this study will be extended to Parker Solar Probe and future missions such as Helioswarm.

Coronal pseudostreamer flux systems have a specific magnetic configuration that influences the morphology and evolution of coronal mass ejections (CMEs) from these regions. Here we continue the analysis of the Wyper et al. (2024, ApJ 975, 168) magnetohydrodynamic simulation of a CME eruption from an idealized pseudostreamer configuration through the construction of synthetic remote-sensing and in-situ observational signatures. We examine the pre-eruption and eruption signatures in extreme ultraviolet and white-light from the low corona through the extended solar atmosphere. We calculate synthetic observations corresponding to several Parker Solar Probe-like trajectories at $\sim$10$R_\odot$ to highlight the fine-scale structure of the CME eruption in synthetic WISPR imagery and the differences between the in-situ plasma and field signatures of flank and central CME-encounter trajectories. Finally, we conclude with a discussion of several aspects of our simulation results in the context of interpretation and analysis of current and future Parker Solar Probe data.

This paper is concerned with the spread in apparent magnitudes (or absolute magnitudes) of main-sequence stars in a star cluster. I specifically consider the effect of binary stars in broadening the main sequence. I present analytic and semi-analytic expressions for the probability density function (pdf) of the magnitude at a given photometric color, including the effects of binarity. The expressions obtained employ plausible models for the pdfs of the magnitudes of the primary and secondary stars, as well as the distribution of secondary-to-primary mass ratio. A crucial parameter is the fraction of stars in the star cluster that are binaries. The resultant formulas can be used to determine the fraction of binaries in a sample of stars taken from a star cluster, or to limit long term variations in the luminosity of solar type stars.

Planetary nebulae (PNe) shown to be members of star clusters provide information on their properties and evolutionary histories that cannot be determined for PNe in the field, in particular the initial masses of their progenitor stars. Here we investigate the bipolar PN PHR J1315-6555 (hereafter PHR J1315), which lies near the open cluster AL 1 (ESO 96-SC04) on the sky. Previous work has established that the PN and cluster have similar radial velocities and amounts of interstellar reddening, and similar distances estimated using independent methods. We have obtained new images of the PN and cluster using the Hubble Space Telescope (HST). Combined with archival HST frames taken 12 years earlier, they provide high-precision proper motions (PMs) for two candidate central stars of PHR J1315. We find that the PMs of both candidates are consistent with those of cluster members, strongly confirming the PN's membership in AL 1. The candidate lying closer to the center of PHR J1315 has the color and luminosity of an early F-type dwarf, suggesting that it may be the optical primary in a close post-common-envelope binary. We used the HST data to construct a color-magnitude diagram for AL 1, which we corrected for significant foreground differential reddening. Isochrone fitting reveals that the cluster lies at a remarkably large distance of about 13 kpc, and has an age of about 1.0 Gyr. The initial mass of the progenitor of PHR J1315 was about 2.1 Msun. We suggest followup investigations that would provide tighter constraints on the object's evolution.

Collisionless shocks, essential for astrophysics, perhaps do not exist as statistically stationary solutions. If so, any quantitative statement about a collisionless shock should be qualified by the age of the shock. A theoretical description of the upstream of the 1+1 dimensional electrostatic collisionless shock is developed -- collisionless hydrodynamics. Peculiarities of collisionless hydrodynamics prevent a shock formation when a piston is driven into cold plasma. An exact self-similar solution is found instead; the spatial extent of the solution grows linearly in time. Direct numerical simulations of plasma kinetics in 1+1 dimensions confirm the hydrodynamic result -- a statistically steady collisionless shock doesn't exist. Instead, at each fixed time, there is a continuous succession in space of marginally stable velocity distribution functions. The spatial support of this continuous succession grows linearly in time.

The distribution of orbital period ratios between adjacent observed exoplanets is approximately uniform, but exhibits a strong falloff toward close orbital separations. We show that this falloff can be explained through past dynamical instabilities carving out the period ratio distribution. Our suite of numerical experiments would have required $\sim 3$ million CPU-hours through direct N-body integrations, but was achieved with only $\approx 50$ CPU-hours by removing unstable configurations using the Stability of Planetary Orbital Configurations Klassifier (SPOCK) machine learning model. This highlights the role of dynamical instabilities in shaping the observed exoplanet population, and shows that the inner part of the period ratio distribution provides a valuable observational anchor on the giant impact phase of planet formation.

Global, three-dimensional, magnetohydrodynamic simulations of Rayleigh-Taylor instabilities at the disk-magnetosphere boundary of rotating, magnetized, compact stellar objects reveal that accretion occurs in three regimes: the stable regime, the chaotic unstable regime, and the ordered unstable regime. Here we track stochastic fluctuations in the pulse period $P(t)$ and aperiodic X-ray luminosity $L(t)$ time series of 24 accretion-powered pulsars in the Small Magellanic Cloud using an unscented Kalman filter to analyze Rossi X-ray Timing Explorer data. We measure time-resolved histories of the magnetocentrifugal fastness parameter $\omega(t)$ and we connect $\omega(t)$ with the three Rayleigh-Taylor accretion regimes. The 24 objects separate into two distinct groups, with 10 accreting in the stable regime, and 14 accreting in the ordered unstable regime. None of the 24 objects except SXP 293 visit the chaotic unstable regime for sustained intervals, although several objects visit it sporadically. The Kalman filter output also reveals a positive temporal cross-correlation between $\omega(t)$ and the independently measured pulse amplitude $A(t)$, which agrees with simulation predictions regarding the pulse-forming behavior of magnetospheric funnel flows in the three accretion regimes.

The study of dark matter substructure through strong gravitational lensing has shown enormous promise in probing the properties of dark matter on sub-galactic scales. This approach has already been used to place strong constraints on a wide range of dark matter models including self-interacting dark matter, fuzzy dark matter and warm dark matter. A major source of degeneracy exists between suppression of low mass halos due to novel dark matter physics and the strength of tidal stripping experienced by subhalos. We study theoretical predictions for the statistical properties of subhalos in strong gravitational lenses using the semi-analytic galaxy formation toolkit: \galns. We present a large suite of dark matter only \gal models, spanning nearly two orders of magnitude in host halo mass (from Milky Way to group mass halos between redshifts from $0.2$ to $0.8$). Additionally, we include a smaller set of \gal runs with the potential of a central massive elliptical to complement our dark matter only suite of models. We place particular focus on quantities relevant to strong gravitational lensing; namely the projected number density of substructure near the Einstein radius as function of host stellar mass and redshift. In the innermost region in projection, we find that our \gal models agrees with N-body simulations within a factor of $\sim 2$ within the Einstein radius. We find that the addition of a central galaxy suppresses the projected number density of subhalos within in the Einstein radius by around $15\%$ relative to dark matter only simulations.

It has been suggested that star-forming galaxies may host a substantial, dark reservoir of gas in the form of planetary-mass molecular clouds that are so cold that $\text{H}_{2}$ can condense. Here we investigate the process of tidal disruption of such "snow clouds" by close passage of field stars. We construct a suite of simulations using the hydrodynamic formalism introduced by Carter and Luminet, and use it to explore the properties of the resulting tidal debris. The debris streams are tiny structures that are highly over-pressured relative to the ambient ISM. They are also unusual in their composition -- initially consisting of cold, gaseous He together with $\text{H}_{2}$ "snowballs" that may be as much as a metre in size. Each stream expands and cools and is subsequently shocked as it ploughs through the ISM; the snowballs are gradually eroded by the shocked gas. Snowballs streaming through the shocked ISM create microstructured plasma that is somewhat reminiscent of the "scattering screens" revealed by radio-wave scintillation studies. However, the tidal disruption rate is too low to account for the observed number of scattering screens if, as we assume here, the stars and clouds have no prior physical association so that disruptions occur as a result of chance encounters between stars and clouds.

Black hole X-ray binaries (BHXBs) are observed in various wavelengths from radio to GeV gamma-ray. Several BHXBs, including MAXI J1820+070 and Cygnus X-1, are also found to emit ultrahigh-energy (UHE; photon energy $>$100 TeV) gamma rays. The origin and production mechanism of the multi-wavelength emission of BHXBs are under debate. We propose a scenario where relativistic particles from magnetically arrested disks (MADs), which could form when BHXBs are in quiescent or hard states, produce UHE gamma rays, while electrons in the jets produce GeV gamma-ray emission. Specifically, magnetic turbulence in MADs heats up and accelerates electrons and protons, while magnetic reconnection in jets accelerates electrons. Sub-PeV gamma rays and neutrinos are produced when relativistic protons interact with the thermal protons and the radiation by thermal electrons in the disk. We discuss the perspectives of observing sub-PeV multi-messenger signals from individual BHXBs. Finally, we evaluate the integrated fluxes of the quiescent and hard-state BHXB population and find that BHXBs may contribute to the Galactic diffuse emission above $\sim 100$ TeV.

Wide-field surveys have markedly enhanced the discovery and study of solar system objects (SSOs). The 2.5-meter Wide Field Survey Telescope (WFST) represents the foremost facility dedicated to optical time-domain surveys in the northern hemisphere. To fully exploit WFST's capabilities for SSO detection, we have developed a heliocentric-orbiting objects processing system (HOPS) tailored for identifying these objects. This system integrates HelioLinC3D, an algorithm well suited for the WFST survey cadence, characterized by revisiting the same sky field twice on the majority of nights. In this paper, we outline the architecture and processing flow of our SSO processing system. The application of the system to the WFST pilot survey data collected between March and May 2024 demonstrates exceptional performance in terms of both temporal efficiency and completeness. A total of 658,489 observations encompassing 38,520 known asteroids have been documented, and 241 newly discovered asteroids have been assigned provisional designations. In particular, 27% of these new discoveries were achieved using merely two observations per night on three nights. The preliminary results not only illuminate the effectiveness of integrating HelioLinC3D within the SSO processing system, but also emphasize the considerable potential contributions of WFST to the field of solar system science.

Magnetohydrodynamic disk-winds play a key role in the formation of massive stars by providing the fine-tuning between accretion and ejection, where excess angular momentum is redirected away from the disk, allowing further mass growth. However, only a limited number of disk-wind sources have been detected to date. To better constrain the exact mechanism of this phenomenon, expanding the sample is critical. We performed an analysis of the disk-wind candidate G11.92-0.61 MM1 by estimating the physical parameters of the system and constraining the wind-launching mechanism. Atacama Large Millimeter/submillimeter Array (ALMA) Band 6 observations were conducted in September 2021 with ALMA's longest baselines, which provided a synthesised beam of $\sim$30 mas. We obtained high-resolution images of the CH$_3$CN ($v_8=1$ and $v=0$), CH$_3$OH, SO$_2$, and SO molecular lines, and the 1.3 mm continuum. Our molecular data allowed us to refine the parameters of the disk-outflow system in MM1. The rotating disk is resolved into two regions with distinct kinematics: the inner region ($<$300 au) is traced by high-velocity emission of high-excitation CH$_3$CN lines and shows a Keplerian rotation; the outer region ($>$300 au), traced by mid-velocity CH$_3$CN emission, rotates in a sub-Keplerian regime. The central source is estimated to be $\sim$20 $M_{\odot}$, which is about half the mass estimated in previous lower-resolution studies. A strong collimated outflow is traced by SO and SO$_2$ emission up to $\sim$3400 au around MM1a. The SO and SO$_2$ emissions show a rotation-dominated velocity pattern, a constant specific angular momentum, and a Keplerian profile that suggests a magneto-centrifugal disk-wind origin with launching radii of $\sim$50-100 au. G11.92-0.61 MM1 appears to be one of the clearest cases of molecular line-traced disk-winds detected around massive protostars.

The formation and mass distribution of the first stars depend on various environmental factors in the early universe. We compare 120 cosmological hydrodynamical simulations to explore how the baryonic streaming velocity (SV) relative to dark matter affects the formation of the first stars. We vary SV from zero to three times its cosmic root-mean-square value, $v_{\rm SV}/\sigma_{\rm SV}=0-3$, and identify 20 representative halos from cosmological simulations. For each model, we follow the evolution of a primordial star-forming cloud from the first appearance of a dense core (with gas density > $10^{6}\,{\rm cm^{-3}}$) until 2 Myr later. In each model, higher SV systematically delays the formation of primordial clouds, formed inside more massive halos ($10^{5}-10^{7}\,{\rm M}_\odot$), and promotes cloud-scale fragmentation and multiple-core formation. The number and total mass of dense cores increase with increasing SV. More than half of models with $v_{\rm SV}/\sigma_{\rm SV} \ge 1.5$ form three or more dense cores in a single halo. In extreme cases, up to 25 cores form at once, which leaves a massive first star cluster. On the other hand, models with $v_{\rm SV}/\sigma_{\rm SV} \leq 1$ form only one or two cores in a halo. In addition, HD-cooling is often enabled in models with low SV, especially in low-z, where HD-cooling is enabled in more than 50% of models. This leads to the formation of the low-mass first star. SV shapes the resulting initial mass function of the first stars and plays a critical role in setting the star-forming environment of the first galaxies.

Recent high-resolution observations have enabled detailed investigations of the circumstellar environments around Class 0/I protostars. Several studies have reported that the infall velocity of the envelope is a few times smaller than the free-fall velocity inferred from protostellar masses estimated via the observed rotational velocity of their Keplerian disks. To explore the physical origins of the slow infall, we perform a set of three-dimensional non-ideal magnetohydrodynamic simulations of the star formation process, extending to $10^5$ yr after protostar formation. Our simulations show that the infall velocity decreases markedly at the outer edge of the pseudo-disk (at radii of $\sim \! 100-1000$ au) and is much slower than the expected free-fall velocity. The degree of this reduction depends on (1) the initial magnetic field strength, (2) the alignment between the initial field and the rotation axis, and (3) the evolutionary stage of the system. Across our parameter space, the ratio of the infall velocity to the free-fall velocity is as small as $0.2-0.5$, which is consistent with the observations. We further examine the reliability of protostellar mass estimates derived from infall and rotational velocities. While the mass derived from disk rotation closely matches the true value, deviation by a factor of $0.3-2$ is found for the estimates using the infall velocity; it is underestimated due to slow infall, but could also be overestimated due to the contribution of disk mass. These findings underscore the critical role of magnetic fields in shaping star formation dynamics and highlight the uncertainties associated with protostellar mass estimates.

We propose that the core mass function (CMF) can be driven by filament fragmentation. To model a star-forming system of filaments and fibers, we develop a fractal and turbulent tree with a fractal dimension of 2 and a Larson's law exponent ($\beta$) of 0.5. The fragmentation driven by convergent flows along the splines of the fractal tree yields a Kroupa-IMF-like CMF that can be divided into three power-law segments with exponents $\alpha$ = $-0.5$, $-1.5$, and $-2$, respectively. The turnover masses of the derived CMF are approximately four times those of the Kroupa IMF, corresponding to a star formation efficiency of 0.25. Adopting $\beta=1/3$, which leads to fractional Brownian motion along the filament, may explain a steeper CMF at the high-mass end, with $\alpha=-3.33$ close to that of the Salpeter IMF. We suggest that the fibers of the tree are basic building blocks of star formation, with similar properties across different clouds, establishing a common density threshold for star formation and leading to a universal CMF.

The gravitational wave (GW) event S230518h is a potential binary neutron star-black hole merger (NSBH) event that was detected during engineering run 15 (ER15), which served as the commissioning period before the LIGO-Virgo-KAGRA (LVK) O4a observing run. Despite its low probability of producing detectable electromagnetic emissions, we performed extensive follow-up observations of this event using the GECKO telescopes in the southern hemisphere. Our observation covered 61.7\% of the 90\% credible region, a $\rm 284\:deg^2$ area accessible from the southern hemisphere, reaching a median limiting magnitude of $R=21.6$ mag. In these images, we conducted a systematic search for an optical counterpart of this event by combining a CNN-based classifier and human verification. We identified 128 transient candidates, but no significant optical counterpart was found that could have caused the GW signal. Furthermore, we provide feasible KN properties that are consistent with the upper limits of observation. Although no optical counterpart was found, our result demonstrates both GECKO's efficient wide-field follow-up capabilities and usefulness for constraining properties of kilonovae from NSBH mergers at distances of $\sim 200$ Mpc.

Early dust evolution in protoplanetary disks is dominated by sticking collisions. However, this initial phase of particle growth faces constraints - notably from destructive encounters. To find the maximum particle size achievable, we studied collisional processes during a prolonged microgravity experiment aboard a suborbital flight. Here, we specifically report an impact erosion limit. We observed individual basalt beads, each measuring 0.5 mm in diameter, colliding with and either eroding or adhering to a cluster several centimeters in size. This cluster, formed from tribocharged particles, simulates an electrostatic growth phase that surpasses the classical bouncing barrier. We find a threshold velocity of about 0.5 m/s, distinguishing between additive and erosive impacts of individual beads. Numerical simulations of grain impacts into clusters, testing both low and high charge constituents corroborate the experimental findings of surface erosion within the observed velocity range. This specific velocity threshold suggests the potential formation of pebbles several centimeters in size within protoplanetary disks. Such dimensions place these pebbles well into a regime where hydrodynamic interaction might facilitate the formation of planetesimals.

We aim to compute the impact rates for objects with a diameter of 1 km onto the regular satellites of Jupiter, Saturn and Uranus using our latest dynamical simulations of the evolution of outer solar system coupled with the best estimates of the current population of objects beyond Neptune and their size-frequency distribution. We use the outcome of the last 3.5~Gyr of evolution of the outer solar system from our database of simulations and combine this with observational constraints of the population beyond Neptune to compute the flux of objects entering the Centaur region, with uncertainties. The initial conditions resemble the current population rather than a near-circular, near-planar disc usually assumed just before the onset of giant planet migration. We obtain a better estimate of the impact probability of a Centaur with the satellites from enacting simulations of planetesimals flying past the satellites on hyperbolic orbits, which agree with literature precedents. We find that our impact rate of objects greater than 1 km in diameter with Jupiter is 0.0012/yr, which is a factor of 3--6 lower than previous estimates of 0.0044/yr from Nesvorny et al. (2023) and 0.0075/yr from Zahnle et al. (2003). On the other hand our impact probabilities with the satellites scaled to the giant planets are consistent with these earlier literature estimates, as is the leakage rate of objects from beyond Neptune into the Centaur region. However, our absolute impact probabilities with the giant planets are lower. We attribute this to our choice of initial conditions.

Cosmic-ray acceleration processes in astrophysical plasmas are often investigated with fully-kinetic or hybrid kinetic numerical simulations, which enable us to describe a detailed microphysics of particle energization mechanisms. Tracing of individual particles in such simulations is especially useful in this regard. However, visually inspecting particle trajectories introduces a significant amount of bias and uncertainty, making it challenging to pinpoint specific acceleration mechanisms. Here, we present a novel approach utilising neural networks to assist in the analysis of individual particle data. We demonstrate the effectiveness of this approach using the dataset from our recent particle-in-cell (PIC) simulations of non-relativistic perpendicular shocks that consists of 252,000 electrons, each characterised by their position, momentum and electromagnetic field at particle's position, recorded in a time series of 1200 time steps. These electrons cross a region affected by the electrostatic Buneman instability, and a small percentage of them attain high energies. We perform classification, regression, and anomaly detection algorithms on the dataset by using a convolutional neural network, a multi-layer perceptron, and an autoencoder. Despite the noisy and imbalanced dataset, all methods demonstrate the capability to differentiate between thermal and accelerated electrons with remarkable accuracy. The proposed methodology may considerably simplify particle classification in large-scale PIC and hybrid simulations.

The primary aim of this research is to evaluate several convolutional neural network-based object detection algorithms for identifying oscillation-like patterns in light curves of eclipsing binaries. This involves creating a robust detection framework that can effectively process both synthetic light curves and real observational data. The study employs several state-of-the-art object detection algorithms, including Single Shot MultiBox Detector, Faster Region-based Convolutional Neural Network, You Only Look Once, and EfficientDet besides a custom non-pretrained model implemented from scratch. Synthetic light curve images and images derived from observational TESS light curves of known eclipsing binaries with a pulsating component were constructed with corresponding annotation files using custom scripts. The models were trained and validated on established datasets, followed by testing on unseen {\it{Kepler}} data to assess their generalization performance. The statistical metrics are also calculated to review the quality of each model. The results indicate that the pre-trained models exhibit high accuracy and reliability in detecting the targeted patterns. Faster R-CNN and You Only Look Once, in particular, showed superior performance in terms of object detection evaluation metrics on the validation dataset such as mAP value exceeding 99\%. Single Shot MultiBox Detector, on the other hand, is the fastest although it shows slightly lower performance with a mAP of 97\%. These findings highlight the potential of these models to contribute significantly to the automated determination of pulsating components in eclipsing binary systems, facilitating more efficient and comprehensive astrophysical investigations.

We develop a new semi-analytic framework of Population (Pop) III and subsequent galaxy formation designed to run on dark matter halo merger trees. In our framework, we consider the effect of the Lyman-Werner flux from Pop III and II stars and the dark matter baryon streaming velocity on the critical minihalo mass for the Pop III formation. Our model incorporates the Lyman-Werner feedback in a self-consistent way, therefore, the spatial variation of Lyman-Werner feedback naturally emerges. The Pop III mass depends on the properties of a minihalo as reproducing radiative hydrodynamical simulation results. We perform statistical studies of Pop III stars by applying this framework to high-resolution cosmological N-body simulations with a maximum box size of 16 Mpc/h and enough mass resolution to resolve Pop III-forming minihalos. A top-heavy initial mass function emerges and two peaks corresponding to the H$_2$ ($20 \lesssim z \lesssim 25$) and atomic cooling halos ($z \lesssim 15$) exist in the distribution. Supermassive stars can be formed in the atomic cooling halos, and the fractions of such supermassive stars increase with the value of streaming velocity. At least an 8 Mpc/h simulation box and the self-consistent model for the Lyman-Werner feedback are necessary to correctly model the Pop III formation in the atomic cooling halos. Our model predicts one supermassive star per halo with several $10^9$ Msun at z=7.5, which is enough to reproduce a high redshift quasar.

Context. The evolution and the surrounding of stripped-envelope supernova progenitors are still under debate: some studies suggest single-star, while others prefer massive binary progenitors. Moreover, the basic physical properties of the exploding star and its interaction with circumstellar matter could significantly modify the overall light curve features of these objects. To better understand the effect of stellar evolution and circumstellar interaction, systematic hydrodynamic calculations are needed. Aims. Here, we test the hypothesis that circumstellar matter generated by an extreme episodic $\eta$ Carinae-like eruption that occurs days or weeks before the supernova explosion may explain the controversies related to the general light curve features of stripped-envelope supernovae. Methods. In this work, we present our bolometric light curve calculations of both single- and binary progenitors generated by hydrodynamic simulations via MESA and SNEC. We also studied the effect of an interaction with a close, low-mass circumstellar matter assumed to be created just a few days or weeks before the explosion. Beyond generating a model light curve grid, we compared our results with some observational data. Results. We found that merely the shape of the supernova light curve could indicate that the cataclysmic death of the massive star happened in a binary system or was related to the explosion of a single star. Moreover, our study also shows that a confined dense circumstellar matter may be responsible for the strange light curve features (bumps, re-brightening, or steeper tail) of some Type Ib/c supernovae.

We present the first results of the HI intensity mapping power spectrum analysis with the MeerKAT International GigaHertz Tiered Extragalactic Exploration (MIGHTEE) survey. We use data covering $\sim$ 4 square degrees in the COSMOS field using a frequency range 962.5 MHz to 1008.42 MHz, equivalent to HI emission in $0.4<z<0.48$. The data consists of 15 pointings with a total of 94.2 hours on-source. We verify the suitability of the MIGHTEE data for HI intensity mapping by testing for residual systematics across frequency, baselines and pointings. We also vary the window used for HI signal measurements and find no significant improvement using stringent Fourier mode cuts. Averaging in the power spectrum domain, i.e. using incoherent averaging, we calculate the first upper limits from MIGHTEE on the HI power spectrum at scales $0.5 Mpc^{-1} \lesssim k \lesssim 10 Mpc^{-1}$. We obtain the best 1$\sigma$ upper limit of 28.6 mK$^{2}$Mpc${^3}$ on $k\sim$2 Mpc$^{-1}$. Our results are consistent with the power spectrum detected with observations in the DEEP2 field with MeerKAT. The data we use here constitutes a small fraction of the MIGHTEE survey and demonstrates that combined analysis of the full MIGHTEE survey can potentially detect the HI power spectrum at $z\lesssim0.5$ in the range $0.1 Mpc^{-1} \lesssim k \lesssim 10 Mpc^{-1}$ or quasi-linear scales.

Gas-phase hydroxylamine (NH$_2$OH) has recently been detected within dense clouds in the interstellar medium. However, it is also likely present within interstellar ices, as well as on the icy surfaces of outer Solar System bodies, where it may react to form more complex prebiotic molecules such as amino acids. In this work, we aimed to provide IR spectra of NH$_2$OH in astrophysical ice analogues that will help in the search for this molecule in various astrophysical environments. Furthermore, we aimed to provide quantitative information on the stability of NH$_2$OH upon exposure to ionizing radiation analogous to cosmic rays, as well as on the ensuing chemistry and potential formation of complex prebiotic molecules. Ices composed of NH$_2$OH, H$_2$O, and CO were prepared by vapor deposition and IR spectra were acquired between $4000-500$~cm$^{-1}$ ($2.5-20 \mu$m) prior and during irradiation using 15 keV protons. In interstellar ice analogues, the most prominent IR absorption band of NH$_2$OH is that at about 1188~cm$^{-1}$, which may be a good candidate to use in searches of this species in icy space environments. Calculated effective destruction cross-sections and G-values for the NH$_2$OH-rich ices studied show that NH$_2$OH is rapidly destroyed upon exposure to ionizing radiation (more rapidly than a number of previously studied organic molecules), and that this destruction is slightly enhanced when it is mixed with other icy species. The irradiation of a NH$_2$OH:H$_2$O:CO ternary ice mixture leads to a rich chemistry that includes the formation of simple inorganic molecules such as NH$_3$, CO$_2$, OCN$^-$, and H$_2$O$_2$, as well as ammonium salts and, possibly, complex organic molecules relevant to life such as formamide, formic acid, urea, and glycine.

I argue that the mixing-heating mechanism of the intracluster medium (ICM) is compatible with the new observations by the X-ray telescope XRISM that show the dispersion velocity in the cooling flow cluster of galaxies A2029 to be 169 km/s. Past jets from the central supermassive black hole induced turbulence; the velocity dispersion value indicates that the jets were powerful, as expected in the mixing-heating mechanism. Although the kinetic energy of the ICM turbulence that XRISM finds is short of heating the ICM and counter radiative cooling, the turbulence is fast enough to mix the hot shocked jets' material with the ICM on time scales shorter than the radiative cooling time. The support of the mixing-heating mechanism from determining the turbulent velocity I claim for A2029 is similar to the conclusion from the X-ray observations by the X-ray telescope Hitomi of the Perseus cluster.

We use a sample of Type 1 active galactic nuclei (AGNs) spectra in order to investigate which atomic processes are responsible for some observed properties of the FeII emission lines and how they are connected with macroscopic physical characteristics of AGN emission regions. We especially focus on the violated relative intensities between different optical FeII lines, whose relative strengths do not follow the expected values according to atomic parameters. We investigated the connection between this effect and the ratio of optical to UV FeII lines (FeII$_{opt}$/FeII$_{UV}$). We divided the optical FeII lines into two large line groups: consistent (FeII$_{cons}$), whose relative intensities are in accordance with their atomic properties, and inconsistent (FeII$_{incons}$), whose relative intensities are significantly stronger than theoretically expected. We fitted the spectra with a flexible optical FeII model, where Fe II lines were divided into several line groups and fitted independently. We focused particularly on understanding the processes that produce strong inconsistent Fe II lines, and therefore, we investigated their correlations with FeII$_{cons}$ and with UV FeII lines. The ratios of FeII$_{incons}$/FeII$_{cons}$ and FeII$_{opt}$/FeII$_{UV}$ increase as the Eddington ratio increases and as the line widths decrease. It is possible that both ratios are affected by the process of self-absorption of stronger lines, which is responsible for the transmission of energy from the UV to the optical FeII emission lines and, analogously, from the FeII$_{cons}$ to the FeII$_{incons}$ lines. In this scenario, the high Eddington ratio causes an increase in the optical depth in FeII lines, which results in the triggering of the process of self-absorption. The observed FeII spectrum is probably a complex mixture of radiation from emission regions with different physical conditions.

Understanding the origin and evolution of supermassive black holes (SMBH) stands as one of the most important challenges in astrophysics and cosmology, with little current theoretical consensus. Improved observational constraints on the cosmological evolution of SMBH demographics are needed. Here we report results of a search via photometric variability for SMBHs appearing as active galactic nuclei (AGN) in the cosmological volume defined by the Hubble Ultra Deep Field (HUDF). This work includes particular focus on a new observation carried out in 2023 with the \textit{Hubble Space Telescope (HST)} using the WFC3/IR/F140W, which is compared directly to equivalent data taken 11 years earlier in 2012. Two earlier pairs of observations from 2009 to 2012 with WFC3/IR/F105W and WFC3/IR/F160W are also analysed. We identify 443, 149, and 78 AGN candidates as nuclear sources that exhibit photometric variability at a level of 2, 2.5 and 3~$\sigma$ in at least one filter. This sample includes 29, 14, and 9 AGN at redshifts $z>6$, when the Universe was $\lesssim900$~Myr old. After variability and luminosity function (down to $M_{\rm UV}=-17\:$mag) completeness corrections, we estimate the co-moving number density of SMBHs, $n_{\rm SMBH}(z)$. At $z = 6 - 9$, $n_{\rm SMBH}\gtrsim 10^{-2}\:{\rm cMpc^{-3}}$. At low-$z$ our observations are sensitive to AGN fainter than $M_{\rm UV}=-17 \:$mag, and we estimate $n_{\rm SMBH}\gtrsim 6\times 10^{-2}\:{\rm cMpc^{-3}}$. We discuss how these results place strong constraints on a variety of SMBH seeding theories.

Among the numerous molecular systems found in the interstellar medium (ISM), vinyl cyanide is the first identified olephinic nitrile. While it has been observed in various sources, its detection in Sgr B2 is notable as the 2$_{11}$-2$_{12}$ rotational transition exhibits maser features. This indicates that local thermodynamic equilibrium conditions are not fulfilled, and an accurate estimation of the molecular abundance in such conditions involves solving the statistical equilibrium equations taking into account the competition between the radiative and collisional processes. This in turn requires the knowledge of rotational excitation data for collisions with the most abundant species - He or H$_2$. In this paper the first three-dimensional CH$_2$CHCN - He potential energy surface is computed using explicitly correlated coupled-cluster theory [(CCSD(T)-F12] with a combination of two basis sets. Scattering calculations of the rotational (de-)excitation of CH$_2$CHCN induced by He atoms are performed with the quantum-mechanical close-coupling method in the low-energy regime. Rotational state-to-state cross sections derived from these calculations are used to compute the corresponding rate coefficients. The interaction potential exhibits a high anisotropy, with a global minimum of $-53.5$ cm$^{-1}$ and multiple local minima. Collisional cross sections are calculated for total energies up to 100 cm$^{-1}$. By thermally averaging the cross-sections, collisional rate coefficients are determined for temperatures up to 20 K. A propensity favouring the transitions with $\Delta k_a=0$ is observed.

Small-scale magnetic flux emergence in the quiet Sun is crucial for maintaining solar magnetic activity. On the smallest scales studied so far, namely within individual granules, two mechanisms have been identified: emergence in tiny magnetic loops and emergence in the form of magnetic flux sheets covering the granule. While there are abundant observations of tiny magnetic loops within granules, the evidence for the emergence of granule-covering magnetic sheets is much more limited. This work aims to statistically analyze magnetic flux sheets, quantify their frequency on the solar surface and their potential contribution to the solar magnetic budget in the photosphere, and investigate the plasma dynamics and granular-scale phenomena associated with their emergence. Using spectro-polarimetric datasets from the solar optical telescope aboard the Hinode satellite and the ground-based Swedish Solar telescope, we developed a two-step method to identify magnetic flux sheet emergence events, detecting magnetic flux patches based on the calculation of the transverse and longitudinal magnetic flux density and associating them with their host granules based on velocity field analysis. We identified 42 events of magnetic flux sheet emergence and characterized their magnetic properties and the associated plasma dynamics of their host granules. Our results align with numerical simulations, indicating a similar occurrence rate. We investigated the relationship between magnetic flux emergence and granular phenomena, finding that flux sheets often emerge in association with standard nascent granules, as well as exploding granules, or granules with granular lanes. In particular, we highlight the potential role of recycled magnetic flux from downflow regions in facilitating flux sheet emergence.

Stellar spectral classification according to the Morgan-Keenan (MK) system remains fundamental to astrophysical studies, yet modern surveys require automated, scalable tools. We present NutMaat, a Python-based package inspired by MKCLASS, designed to automate MK classification while addressing scalability and usability limitations. It employs modern computational tools for batch processing and offers a modular architecture that enables efficient, platform-independent analysis of large spectral datasets. Tested on the CFLIB and MILES libraries, NutMaat achieved spectral and luminosity classification accuracies comparable to MKCLASS, with minimal systematic offsets and a robust performance down to S/N $\approx$ 5. NutMaat successfully identified chemically peculiar stars, tested on LAMOST DR7 ACV variables, and processed the SDSS-IV MaStar library to produce a stellar catalog, demonstrating survey readiness. Future development of NutMaat will focus on extending wavelength coverage, computational acceleration via Cython, and refining peculiarity classification. As an open-source tool, NutMaat bridges traditional MK methods with modern data workflows, offering a scalable solution for current and future spectroscopic surveys.

Peaked-spectrum (PS) sources, known for their distinct peaked radio spectra, represent a type of radio-loud active galactic nuclei (AGN). Among these, megahertz-peaked spectrum (MPS) sources, which exhibit a spectral peak at a frequency of a hundred megahertz, have emerged as a potential tool for identifying high-redshift candidates. However, the potential evolutionary link between the fraction of these sources and redshift remains unclear and requires further investigation. The recent, high sensitivity Low Frequency Array (LOFAR) surveys enable statistical studies of these objects to ultra-low frequencies (< 150 MHz). In this study, we first use the multiradio data to investigate the evolution of spectral index with redshift for 1,187 quasars from the SDSS 16th quasar catalog. For each quasar, we analyze available data from the LOFAR Low Band Antenna (LBA) at 54 MHz, High Band Antenna (HBA) at 144 MHz, and the Very Large Array (VLA) the Faint Images of the Radio Sky at Twenty cm (FIRST) at 1.4 GHz. We measure the spectral index ($\alpha^{144}_{54}$ and $\alpha^{1400}_{144}$) and find no significant change in their median values with the redshift. Extended sources have steeper spectral indices than compact sources, which is consistent with previous findings. Based on the spectral indices information, we identify MPS sources using these criteria: $\rm \alpha^{144}_{54} >= 0.1$ and $\rm \alpha^{1400}_{144} < 0$, and analyze their properties. We find that the fraction of MPS sources is constant with the redshift ($0.1-4.8$), bolometric luminosity ($\rm 10^{44}-10^{48} erg/s$), and supermassive black hole mass ($\rm 10^{7}-10^{10.5} M_{\odot}$), which suggests that MPS sources have relatively stable physical conditions or formation mechanisms across various evolutionary stages and environments.

Aerosols are capable of having a huge influence on reflected, emitted and transmitted planetary spectra. The Near InfraRed Spectrograph (NIRSpec) using the PRISM mode on board of the James Webb Space Telescope (JWST) is providing valuable data of transit spectra over a wide spectral range that is able to cover the whole contribution of aerosols, potentially disentangling them from other constituents and thus allowing to constrain their properties. We aim to investigate whether NIRSpec/PRISM JWST transmission spectroscopy observations, in addition to being useful to detect and determine the abundance of gases more accurately than any previous instruments, are also capable of studying the physical properties of the aerosols in exoplanetary atmospheres. We perform nested sampling Bayesian retrievals with MultiNest library. We use the Planetary Spectrum Generator (PSG) and the Modelled Optical Properties of enSeMbles of Aerosol Particles (MOPSMAP) database as tools for the forward simulations and NIRSpec/PRISM JWST previously published observations of WASP-39b as input data. Retrievals indicate that models including an aerosol extinction weakly increasing or sharply decreasing with wavelength are decisively better than those with a flat transmission and that this increased degree of complexity is supported by the kind of data that JWST/NIRSpec can provide. Given other physical constraints from previous works, the scenario of weakly increasing particle extinction is favoured. We also find that this also has an effect on the retrieved gas abundances. JWST observations have the potential to study some physical characteristics of exoplanetary clouds, in particular their overall dependence of transmissivity with wavelength. In fact, it is important to implement more detailed aerosol models as their extinction may affect significantly retrieved molecular abundances.

We use Spitzer/IRAC deep exposure data covering two significantly larger than before sky areas to construct maps suitable for evaluating source-subtracted fluctuations in the cosmic infrared background (CIB). The maps are constructed using the self-calibration methodology eliminating artifacts to sufficient accuracy and subset maps are selected in each area containing approximately uniform exposures. These maps are clipped and removed of known sources and then Fourier transformed to probe the CIB anisotropies to new larger scales. The power spectrum of the resultant CIB anisotropies is measured from the data to >1 degree revealing the component well above that from remaining known galaxies on scales >1 arcmin. The fluctuations are demonstrated to be free of Galactic and Solar System foreground contributions out to the largest scales measured. We discuss the proposed theories for the origin of the excess CIB anisotropies in light of the new data. Out of these, the model where the CIB fluctuation excess originates from the granulation power due to LIGO-observed primordial black holes as dark matter appears most successful in accounting for all observations related to the measured CIB power amplitude and spatial structure, including the measured coherence between the CIB and unresolved cosmic X-ray background (CXB). Finally we point out the use of the data to probe the CIB-CXB cross-power to new scales and higher accuracy. We also discuss the synergy of these data with future CIB programs at shorter near-IR wavelengths with deep wide surveys and sub-arcsecond angular resolution as provided by Euclid and Roman space missions.

One of the primary objectives of the SPHEREx mission is to understand the origin of molecules such as H2O, CO2, and other volatile compounds at the early stages of planetary system formation. Because the vast majority of these compounds -- typically exceeding 95% -- exist in the solid phase rather than the gaseous phase in the systems of concern here, the observing strategy planned to characterize them is slightly unusual. Specifically, SPHEREx will target highly obscured sources throughout the Milky Way, and observe the species of concern in absorption against background illumination. SPHEREx spectrophotometry will yield ice column density measurements for millions of obscured Milky Way sources of all ages and types. By correlating those column densities with source ages, the SPHEREx mission will shed light on whether those molecules were formed in situ along with their nascent stellar systems, or whether instead they formed elsewhereand were introduced into those systems after their formation. To that end, this work describes version 7.1 of the SPHEREx Target List of Ice Sources (SPLICES) for the community. It contains about 8.6 million objects brighter than W2~12 Vega mag over much of the sky, principally within a broad strip running the length of the Milky Way midplane, but also within high-latitude molecular clouds and even the Magellanic Clouds.

We present the multi-wavelength and environmental properties of 37 variability-selected active galactic nuclei (AGNs), including 30 low luminosity AGNs (LLAGNs), using a high cadence time-domain survey (ASAS-SN) from a spectroscopic sample of 1218 nearby bright galaxies. We find that high-cadence time-domain surveys uniquely select LLAGNs that do not necessarily satisfy other AGN selection methods, such as X-ray, mid-IR, or BPT methods. In our sample, 3% of them pass the mid-infrared color based AGN selection, 18% pass the X-ray luminosity based AGN selection, and 60% pass the BPT selection. This result is supported by two other LLAGN samples from high-cadence time-domain surveys of TESS and PTF, suggesting that the variability selection method from well-sampled light curves can find AGNs that may not be discovered otherwise. These AGNs can have moderate to small amplitudes of variability from the accretion disk, but, of many of them, with no strong corona, emission lines from the central engine, or accretion power to dominate the mid-IR emission. The X-ray spectra of a sub-sample of bright sources are consistent with a power law model. Upon inspecting the environments of our sample, we find that LLAGNs are more common in denser environments of galaxy clusters in contrast with the trend established in the literature for luminous AGNs at low redshifts, which is broadly consistent with our analysis result for luminous AGNs limited by a smaller sample size. This contrast in environmental properties between LLAGN and luminous AGNs suggests that LLAGNs may have different trigger mechanisms.

Infall of interstellar material is a potential non-planetary origin of pressure bumps in protoplanetary disks. While pressure bumps arising from other mechanisms have been numerically demonstrated to promote planet formation, the impact of infall-induced pressure bumps remains unexplored. We aim to investigate the potential for planetesimal formation in an infall-induced pressure bump, starting with sub-micrometer-sized dust grains, and to identify the conditions most conducive to triggering this process. We developed a numerical model that integrates axisymmetric infall, dust drift, and dust coagulation, along with planetesimal formation via streaming instability. Our parameter space includes gas viscosity, dust fragmentation velocity, initial disk mass, characteristic disk radius, infall rate and duration, as well as the location and width of the infall region. An infall-induced pressure bump can trap dust from both the infalling material and the outer disk, promoting dust growth. The locally enhanced dust-to-gas ratio triggers streaming instability, forming a planetesimal belt inside the central infall location until the pressure bump is smoothed out by viscous gas diffusion. Planetesimal formation is favored by a massive, narrow streamer infalling onto a low-viscosity, low-mass, and spatially extended disk containing dust with a high fragmentation velocity. This configuration enhances the outward drift speed of dust on the inner side of the pressure bump, while also ensuring the prolonged persistence of the pressure bump. Planetesimal formation can occur even if the infalling material consists solely of gas. A pressure bump induced by infall is a favorable site for dust growth and planetesimal formation, and this mechanism does not require a preexisting massive planet to create the bump.

Comets, asteroids, and other small bodies are thought to be remnants of the original planetesimal population of the Solar System. As such, their physical, chemical, and isotopic properties hold crucial details on how and where they formed and how they evolved. Yet, placing precise constraints on the formation region of these bodies has been challenging. Data from spacecraft missions have a particularly high potential of addressing the question of the origin of the visited bodies. ESA's Rosetta mission to comet 67P/Churyumov-Gerasimenko returned data from the comet for two years on its journey around the Sun. This extensive data set has revolutionised our view on comets and still holds unsolved problems. We use the Rosetta/ROSINA (Rosetta Orbiter Spectrometer for Ion and Neutral Analysis) measurement of the volatile/ice composition and the Rosetta/COSIMA (COmetary Secondary Ion Mass Analyzer) measurements of the refractory composition of comet 67P. These measurements are combined using a Monte Carlo method. The refractory-to-ice ratio is a free parameter that is constrained a posteriori. Using only the composition, we constrain the refractory-to-ice ratio to $0.5<\chi<1.7$, and derive the bulk elemental abundances for 67P of H, C, N, O, Na, Mg, Al, S, K, Ar, Ca, Cr, Mn, Fe, Kr, and Xe. We find the noble gas xenon in near solar elemental abundance in comet 67P. Krypton is slightly depleted, while argon is heavily depleted. Comet 67P is enriched in all three noble gases by up to 2.5 orders of magnitude compared to CI chondrites. We show this is consistent with a formation region between 25 and 35 au in a protoplanetary disk region with temperatures between 30 and 40 K and with the trapping of dust for a long time in rings of the protoplanetary disk.

We show that primordial adiabatic curvature fluctuations generate an instability of the scalar field sourcing a kination era. We demonstrate that the generated higher Fourier modes constitute a radiation-like component dominating over the kination background after about $11$ e-folds of cosmic expansion. Current constraints on the extra number of neutrino flavors $\Delta N_{\rm eff}$ thus imply the observational bound of approximately 10 e-folds, representing the most stringent bound to date on the stiffness of the equation of state of the pre-Big-Bang-Nucleosynthesis universe.

Cosmological moduli generically come to dominate the energy density of the early universe, and thereby trigger an early matter dominated era. Such non-standard cosmological histories are expected to have profound effects on the evolution and production of axion cold dark matter and dark radiation, as well as their prospects for detection. We consider moduli-axion couplings and investigate the early history of the coupled system, considering closely the evolution of the homogeneous modulus field, the back-reaction from the axion, and the energy densities of the two fields. A particular point of interest is the enhancement of axion production from modulus decay, due to tachyonic and parametric resonant instabilities, and the implications of such production on the cosmological moduli problem, axion dark radiation, and the available parameter space for axion dark matter. Using an effective field theory approach, WKB-based semi-analytical analysis, and detailed numerical estimates of the co-evolution of the system, we evaluate the expected decay efficiency of the modulus to axions. The effects of higher-order operators are studied and implications for UV-complete frameworks such as the Large Volume Scenarios in Type IIB string theory are considered in detail.

This paper presents results for the innermost stable circular orbit in a Kerr-like spacetime. The metric employed is an approximation that combines the Kerr metric with the Erez-Rosen metric, expanded in a Taylor series. Consequently, this spacetime incorporates three relativistic multipole moments: mass, spin, and quadrupole moment. Our derivation builds upon the analysis conducted by Chandrasekhar for the Kerr metric. Utilizing the Euler-Lagrange method and Hamiltonian dynamics, we define an effective potential for the radial coordinate. This equation can be used to measure the mass quadrupole through observational methods, as it yields a quadratic polynomial for the quadrupole moment. As anticipated, the limiting cases of this equation correspond to the established cases of Kerr and Schwarzschild spacetimes.

It is commonly accepted that high energy cosmic rays up to $10^{19}$ eV can be produced in catastrophic astrophysical processes. However the source of a few observed events with higher energies remains mysterious. We propose that they may originate from decay or annihilation of ultra heavy particles of dark matter. Such particles naturally appear in some models of modified gravity related to Starobinsky inflation.

This article briefly discusses the philosophical and technical aspects of AI. It focuses on two concepts of understanding: intuition and causality, and highlights three AI technologies: Transformers, chain-of-thought reasoning, and multimodal processing. We anticipate that in principle AI could form understanding, with these technologies representing promising advancements.

We present an experimental method for the characterization of the kinetic energies of ions confined in a 22-pole radio frequency trap by inducing a small potential barrier using the surrounding ring electrodes, allowing the selective extraction of ions. Energy sampling experiments have been performed on buffer gas thermalized He$^+$ ions at trap temperatures between 10-180 K, resulting in distinct extraction curves as a function of the potential barrier, and a differentiated behavior depending on the escape time from the trap. The experiments are complemented by Monte Carlo simulations of the ion trajectories inside the calculated trap potential and allow us to investigate the properties of the sampling method, the role of ion motion coupling, and the impact of residual buffer gas collisions on the observed results. The technique has also been successfully applied to identify energetic H$_3^+$ ions produced in an exothermic reaction inside the trap. Upon calibration, this method can provide relative kinetic energy distributions or be used to filter the maximum desired kinetic energy of the ions inside the trap.

Cosmological gravitational-wave backgrounds are an exciting science target for next-generation ground-based detectors, as they encode invaluable information about the primordial Universe. However, any such background is expected to be obscured by the astrophysical foreground from compact-binary coalescences. We propose a novel framework to detect a cosmological gravitational-wave background in the presence of binary black holes and binary neutron star signals with next-generation ground-based detectors, including Cosmic Explorer and the Einstein Telescope. Our procedure involves first removing all the individually resolved binary black hole signals by notching them out in the time-frequency domain. Then, we perform joint Bayesian inference on the individually resolved binary neutron star signals, the unresolved binary neutron star foreground, and the cosmological background. For a flat cosmological background, we find that we can claim detection at $5\,\sigma$ level when $\Omega_\mathrm{ref}\geqslant 2.7\times 10^{-12}/\sqrt{T_\mathrm{obs}/\mathrm{yr}}$, where $T_\mathrm{obs}$ is the observation time (in years), which is within a factor of $\lesssim2$ from the sensitivity reached in absence of these astrophysical foregrounds.

In this work, we investigate the properties of string effective theories with scalar field(s) and a scalar potential. We first claim that in most examples known, such theories are multifield, with at least 2 non-compact field directions; the few counter-examples appear to be very specific and isolated. Such a systematic multifield situation has important implications for cosmology. Characterising properties of the scalar potential $V$ is also more delicate in a multifield setting. We provide several examples of string effective theories with $V>0$, where the latter admits an asymptotically flat direction along an off-shell field trajectory: in other words, there exists a limit $\varphi \rightarrow \infty$ for which $\frac{|\partial_{\varphi} V|}{V} \rightarrow 0$. It is thus meaningless to look for a lower bound to this single field quantity in a multifield setting; the complete gradient $\nabla V$ is then better suited. Restricting to on-shell trajectories, this question remains open, especially when following the steepest descent or more generally a gradient flow evolution. Interestingly, single field statements in multifield theories seem less problematic for $V<0$.