Papers for Tuesday, Oct 29 2024

The wolensing Python package offers a solution for gravitational wave lensing computations within the full wave-optics regime. This tool is primarily designed to calculate the gravitational lensing amplification factor including diffractive effects, an essential component for generating accurate lensed gravitational wave waveforms. These waveforms are integral to astrophysical and cosmological studies related to gravitational-wave lensing. Integrating with lensingGW (Pagano, Hannuksela, and Li 2020), wolensing provides solutions for image positions in the high-frequency regime where wave and geometrical optics converge. This functionality allows the amplification factor to be applicable across a wider frequency range. Another key feature of wolensing is its ability to plot time delay contours on the lens plane, offering researchers a visual tool to better understand the relationship between the lens system and the amplification factor. wolensing is compatible with various lens models in lenstronomy (Birrer et al. 2021). There are also built-in lens models including point mass, singular isothermal sphere (SIS),and nonsingular isothermal ellipsoid (NIE) with jax (Bradbury et al. 2018) supporting GPU computation. Users can accommodate different lens models in the code with jax. wolensing is available as an open-source package on PyPI and can be installed via pip.

The Galactic Center hosts a rotating disk of young stars between 0.05 and 0.5 pc of Sgr A*. The ``S-stars'' at a distance $<0.04$ pc, however, are on eccentric orbits with nearly isotropically distributed inclinations. The dynamical origin of the S-star cluster has remained a theoretical challenge. Using a series of $N$-body simulations, we show that a recent massive black hole merger with Sgr A* can self-consistently produce many of the orbital properties of the Galactic nuclear star cluster within 0.5 pc. A black hole merger results in a gravitational wave recoil kick, which causes the surrounding cluster to form an apse-aligned, eccentric disk. We show that stars near the inner edge of an eccentric disk migrate inward and are driven to high eccentricities and inclinations due to secular torques similar to the eccentric Kozai-Lidov mechanism. In our fiducial model, starting with a thin eccentric disk with $e = 0.3$, the initially unoccupied region within $0.04$ pc is populated with high eccentricity, high inclination S-stars within a few Myr. This dynamical channel would suggest that a black hole of mass $2^{+3}_{-1.2} \times 10^5 \ M_{\odot}$ merged with Sgr A* within the last 10 Myr.

We apply the stellar locus method to synthetic $(BP-RP)_{XPSP}$ and $(BP-G)_{XPSP}$ colors derived from corrected Gaia BP/RP (XP) spectra to obtain accurate and precise estimates of metallicity for about 100 million stars in the Milky Way (34 million giants in the color range $0.6 < (BP-RP)_0 < 1.75$ and 65 million dwarfs in the color range $0.2 < (BP-RP)_0 < 1.5$). The sub milli-magnitude precision of the derived synthetic stellar colors enables estimates of metallicity for stars as low as [Fe/H] $\sim -4$. Multiple validation tests indicate that the typical metallicity precision is between 0.05 -- 0.1 dex for both dwarfs and giants at [Fe/H] = 0 as faint as G $\sim$ 17, and decreases to 0.15 -- 0.25 dex at [Fe/H] = $-$2.0. For $-4.0 <$ [Fe/H] $ < -3.0$, the typical metallicity precision decreases to on the order of 0.4 -- 0.5 dex, based on the results from the training sample. Our achieved precision is comparable to or better than previous efforts using the entire XP spectra, and about three times better than our previous work based on Gaia EDR3 colors. This opens up new opportunities for investigations of stellar populations, the formation and chemical evolution of the Milky Way, the chemistry of stars and star clusters, and the identification of candidate stars for subsequent high-resolution spectroscopic follow-up.

The intracluster medium (ICM) today is comprised largely of hot gas with clouds of cooler gas of unknown origin and lifespan. We analyze the evolution of cool gas (temperatures $\lesssim10^{4.5}$ K) in the ICM of 352 galaxy clusters from the TNG-Cluster simulations, with present-day mass $\sim10^{14.3-15.4}\,{\rm M_\odot}$. We follow the main progenitors of these clusters over the past $\sim13$ billion years (since $z\lesssim7$) and find that, according to TNG-Cluster, the cool ICM mass increases with redshift at fixed cluster mass, implying that this cooler past of the ICM is due to more than just halo growth. The cool cluster gas at $z\lesssim0.5$ is mostly located in and around satellite galaxies, while at $z\gtrsim2$ cool gas can also accrete via filaments from the intergalactic medium. Lower-mass and higher-redshift clusters are more susceptible to cooling. The cool ICM mass correlates with the number of gaseous satellites and inversely with the central supermassive black hole (SMBH) mass. The average number of gaseous satellites decreases since $z=2$, correlating with the decline in the cool ICM mass over cosmic time, suggesting a link between the two. Concurrently, kinetic SMBH feedback shifts the ICM temperature distribution, decreasing the cool ICM mass inside-out. At $z\approx0.5$, the predicted MgII column densities are in the ballpark of recent observations, where satellites and other halos contribute significantly to the total MgII column density. Suggestively, a non-negligible amount of the ICM cool gas forms stars in-situ at early times, reaching $\sim10^{2}\,{\rm M_\odot}\,{\rm yr^{-1}}$ and an H$\alpha$ surface brightness of $\sim10^{-17}\,{\rm~erg\,s^{-1}\,cm^{-2}\,arcsec^{-2}}$ at $z\approx2$, detectable with Euclid and JWST.

Stellar-mass black holes (BHs) embedded in active galactic nuclei (AGN) may be major sources of astrophysical gravitational waves (GWs), contributing both to the observed LIGO-Virgo-KAGRA population of binary BH mergers and to future populations of LISA-band extreme mass ratio inspirals (EMRIs). The ability of these BHs to pair up into binaries, inspiral, and produce GWs will be shaped by the existence of migration traps, regions in the AGN where hydrodynamic torques vanish. Previous works have studied the existence and location of migration traps in AGN disks. Here, we investigate how individual BHs may escape such traps as an outcome of mergers, potentially suppressing hierarchical growth. We find that while GW recoil kicks are strong enough to kick merged BHs onto inclined orbits, gas drag quickly realigns them into the AGN disk. A more robust escape mechanism is gap opening: once a BH grows above a critical mass, its gravity disturbs the AGN gas sufficiently to eliminate the trap. In low-mass AGN relevant for LISA, gaps open easily and the resulting ``wet EMRI'' masses are unlikely to reflect protracted hierarchical mergers. In combination with our previous work, we find that migration traps only exist in a relatively narrow range of AGN luminosities between [10^{43.5},10^{45.5}] erg/s. We identify an even narrower AGN luminosity range for which stellar mass BHs can grow into the pair instability mass gap and beyond. This characteristic luminosity scale may assist in indirect tests of the ``AGN channel'' for binary BH mergers.

Outflows are critical components of many astrophysical systems, including accreting compact binaries and active galactic nuclei (AGN). These outflows can significantly affect a system's evolution and alter its observational appearance by reprocessing the radiation produced by the central engine. Sirocco (Simulating Ionization and Radiation in Outflows Created by Compact Objects - or "the code formerly known as Python") is a Sobolev-based Monte Carlo ionization and radiative transfer code. It is designed to simulate the spectra produced by any system with an azimuthally-symmetric outflow, from spherical stellar winds to rotating, biconical accretion disc winds. Wind models can either be parametrized or imported, e.g. from hydrodynamical simulations. The radiation sources include an optically thick accretion disc and various central sources with flexible spectra and geometries. The code tracks the "photon packets" produced by the sources in any given simulation as they traverse and interact with the wind. The code assumes radiative near-equilibrium, so the thermal and ionization state can be determined iteratively from these interactions. Once the physical properties in the wind have converged, Sirocco can be used to generate synthetic spectra at a series of observer sightlines. Here, we describe the physical assumptions, operation, performance and limitations of the code. We validate it against tardis, cmfgen and cloudy, finding good agreement, and present illustrative synthetic spectra from disc winds in cataclysmic variables, tidal disruption events, AGN and X-ray binaries. Sirocco is publicly available on GitHub, alongside its associated data, documentation and sample input files covering a wide range of astrophysical applications.

More than a decade ago, the Large Area Telescope aboard the Fermi Gamma-ray Space Telescope unveiled the existence of two gigantic gamma-ray lobes known as the Fermi bubbles. While their origin is still unknown, various studies identified intricate spectral and morphological structures within the bubbles. One peculiar region, the cocoon, has recently been associated with gamma-ray emissions from the Sagittarius dwarf spheroidal (Sgr) galaxy. We assess the validity of this claim through adaptive-template fitting and pixel-count statistical methods. Our approach introduces a substantial advancement in data interpretation by enabling a data-driven optimisation of astrophysical background models, thereby reducing the impact of background mis-modelling. We do not find evidence for gamma-ray emission from the Sgr region at the level obtained in previous work. We demonstrate that there is no pronounced difference between the source population located within the cocoon region and a reference region at similar latitudes. We examine the hypothesis that a millisecond pulsar population in Sgr causes the putative signal via dedicated simulations finding that it is unlikely that we may have detected their presence in Sgr.

By analyzing the spectral energy distributions (SEDs) of resolved stars in nearby galaxies, we can constrain their stellar properties and line-of-sight dust extinction. From the Scylla survey, we obtain ultraviolet to near-infrared photometry from Wide Field Camera 3 onboard the {\it Hubble Space Telescope} for more than 1.5 million stars in the SMC and LMC. We use the Bayesian Extinction and Stellar Tool (BEAST) to analyze the multi-band SEDs of these sources and characterize their initial masses, ages, metallicities, distances, and line-of-sight extinction properties (e.g.~$A_V$, $R_V$). We apply quality cuts and perform validation simulations to construct a catalog of over 550,000 stars with high-reliability SED fits, which we use to analyze the stellar content and extinction properties of the SMC and LMC. We detect stars with masses as low as 0.6 $M_{\odot}$. Observed stellar age distributions show a jump in stars around 6 Gyrs ago, which is in agreement with other star-formation histories. Extinctions ($A_V$) in both galaxies follow a log-normal distribution. We compare $A_V$ with ancillary gas and dust tracers like $HI$, $H_\alpha$, and far infrared (FIR) dust emission and find positive correlations on a field-by-field basis. We convert observed $A_V$ to predicted dust surface densities using the Draine et. al. (2014) model and find $A_V$-based dust surface densities are a factor of $\sim$2.5 lower than observed FIR-based dust surface densities, a correction factor similar to other studies.

Abridged: The fortunate proximity of the SN2023ixf allowed astronomers to follow its evolution from almost the moment of the collapse of the progenitor's core. SN2023ixf can be explained as an explosion of a massive star with an energy of 0.7e51 erg, however with a greatly reduced envelope mass, probably because of binary interaction. In our radiative-transfer simulations, the SN ejecta of 6 Msun interact with circumstellar material (CSM) of ~0.6 Msun extending to 1.e15 cm, which results in a light curve (LC) peak matching that of SN2023ixf. The origin of this required CSM might be gravity waves originating from convective shell burning, which could enhance wind-like mass-loss during the late stages of stellar evolution. The steeply rising, low-luminosity flux during the first hours after observationally confirmed non-detection, however, cannot be explained by the collision of the energetic SN shock with the CSM. Instead, we considered it as a precursor that we could fit by the emission from ~0.5 Msun of matter that was ejected with an energy of 1.e49 erg a fraction of a day before the main shock of the SN explosion reached the surface of the progenitor. The source of this energy injection into the outermost shell of the stellar envelope could also be dynamical processes related to the convective activity in the progenitor's interior or envelope. Alternatively, the early rise of the LC could point to the initial breakout of a highly non-spherical SN shock or of fast-moving, asymmetrically ejected matter that was swept out well ahead of the SN shock, potentially in a low-energy, nearly relativistic jet. We also discuss that pre-SN outbursts and LC precursors can be used to study or to constrain energy deposition in the outermost stellar layers by the decay of exotic particles, such as axions, which could be produced simultaneously with neutrinos in the newly formed, hot neutron star.

Gravitational-wave (GW) parameter estimation typically assumes that instrumental noise is Gaussian and stationary. Obvious departures from this idealization are typically handled on a case-by-case basis, e.g., through bespoke procedures to ``clean'' non-Gaussian noise transients (glitches), as was famously the case for the GW170817 neutron-star binary. Although effective, manipulating the data in this way can introduce biases in the inference of key astrophysical properties, like binary precession, and compound in unpredictable ways when combining multiple observations; alternative procedures free of the same biases, like joint inference of noise and signal properties, have so far proved too computationally expensive to execute at scale. Here we take a different approach: rather than explicitly modeling individual non-Gaussianities to then apply the traditional GW likelihood, we seek to learn the true distribution of instrumental noise without presuming Gaussianity and stationarity in the first place. Assuming only noise additivity, we employ score-based diffusion models to learn an empirical noise distribution directly from detector data and then combine it with a deterministic waveform model to provide an unbiased estimate of the likelihood function. We validate the method by performing inference on a subset of GW parameters from 400 mock observations, containing real LIGO noise from either the Livingston or Hanford detectors. We show that the proposed method can recover the true parameters even in the presence of loud glitches, and that the inference is unbiased over a population of signals without applying any cleaning to the data. This work provides a promising avenue for extracting unbiased source properties in future GW observations over the coming decade.

Chiral effective field theory ($\chi$EFT) has proved to be a powerful microscopic framework for predicting the properties of neutron-rich nuclear matter with quantified theoretical uncertainties up to about twice the nuclear saturation density. Tests of $\chi$EFT predictions are typically performed at low densities using nuclear experiments, with neutron star (NS) constraints only being considered at high densities. In this work, we discuss how asteroseismic quasi-normal modes within NSs could be used to constrain specific matter properties at particular densities, not just the integrated quantities to which bulk NS observables are sensitive. We focus on the crust-core interface mode, showing that measuring this mode's frequency would provide a meaningful test of $\chi$EFT at densities around half the saturation density. Conversely, we use nuclear matter properties predicted by $\chi$EFT to estimate that this mode's frequency is around 185 $\pm$ 50 Hz. Asteroseismic observables such as resonant phase shifts in gravitational-wave signals and multimessenger resonant shattering flare timings, therefore, have the potential to provide useful tests of $\chi$EFT.

X-ray polarization is a unique new probe of the particle acceleration in astrophysical jets made possible through the Imaging X-ray Polarimetry Explorer. Here we report on the first dense X-ray polarization monitoring campaign on the blazar Mrk 421. Our observations were accompanied by an even denser radio and optical polarization campaign. We find significant short-timescale variability in both X-ray polarization degree and angle, including a $\sim90^\circ$ angle rotation about the jet axis. We attribute this to random variations of the magnetic field, consistent with the presence of turbulence but also unlikely to be explained by turbulence alone. At the same time, the degree of lower-energy polarization is significantly lower and shows no more than mild variability. Our campaign provides further evidence for a scenario in which energy-stratified shock-acceleration of relativistic electrons, combined with a turbulent magnetic field, is responsible for optical to X-ray synchrotron emission in blazar jets.

We present a data reduction pipeline (DRP) for Keck/NIRSPEC built as an addition to the PypeIt Python package. The DRP is capable of reducing multi-order echelle data taken both before and after the detector upgrade in 2018. As part of developing the pipeline, we implemented major improvements to the capabilities of the PypeIt package, including manual wavelength calibration for multi-order data and new output product that returns a coadded spectrum order-by-order. We also provide a procedure for correcting telluric absorption in NIRSPEC data by using the spectra of telluric standard stars taken near the time of the science spectra. At high resolutions, this is often more accurate than modeling-based approaches.

The relationship between magnetic activity and Rossby number is one way through which stellar dynamos can be understood. Using measured rotation rates and X-ray to bolometric luminosity ratios of an ensemble of stars, we derive empirical convective turnover times based on recent observations and re-evaluate the X-ray activity-Rossby number relationship. In doing so, we find a sharp rise in the convective turnover time for stars in the mass range of $0.35-0.4\ \rm M_{\odot}$, associated with the onset of a fully convective internal stellar structure. Using $\texttt{MESA}$ stellar evolution models, we infer the location of dynamo action implied by the empirical convective turnover time. The empirical convective turnover time is found to be indicative of dynamo action deep within the convective envelope in stars with masses $0.1-1.2\ \rm M_{\odot}$, crossing the fully convective boundary. Our results corroborate past works suggesting that partially and fully convective stars follow the same activity-Rossby relation, possibly owing to similar dynamo mechanisms. Our stellar models also give insight into the dynamo mechanism. We find that empirically determined convective turnover times correlate with properties of the deep stellar interior. These findings are in agreement with global dynamo models that see a reservoir of magnetic flux accumulate deep in the convection zone before buoyantly rising to the surface.

The apparent angular positions of quasars are deflected on the sky by the gravitational field sourced by foreground matter. This weak lensing effect is measurable through the distortions it introduces in the lensed quasar spectra. Discrepancies in the statistics of the Lyman-$\alpha$ forest spectral absorption features can be used to reconstruct the foreground lensing potential. We extend the study of this method of Lyman-$\alpha$ forest weak gravitational lensing to lower angular forest spectrum source densities than previous work. We evaluate the performance of the Lyman-$\alpha$ lensing estimator of Metcalf et al. (2020) on mock data based on the angular forest source density ($50$ per square degree) and volume ($\sim$700,000 spectra total) of the DESI survey. We simulate the foreground galaxy distribution and lensing potentials with redshift evolution approximated by N-body simulation and simulate Gaussian-random Lyman-$\alpha$ forests to produce mock data for the entire DESI footprint. By correlating the foreground galaxy distribution with the potential reconstructed by the estimator, we find that a weak lensing detection with signal to noise of $\sim4$ will be possible with the full DESI data. We show that spectral surveys with low density and high volume are promising candidates for forest weak lensing in addition to the high resolution data that have been considered in previous work. We present forecasts for future spectral surveys and show that with larger datasets a detection with signal to noise $>10$ will be possible.

The main objective of this paper is to fully study 1:1 mean-motion resonance in the Solar System. We calculated stability points applying a resonant semi-analytic theory valid for any value of eccentricity or inclination. The location of each equilibrium point changes as the orbital elements of an object change, which led us to map the location of them. For the case of low inclination and low eccentricity we recovered the known L4 and L5 points. The three global types of orbits for this resonance; Tadpoles, Quasi-satellites and Horseshoes, vary as a function of the orbital elements, even disappearing for some cases. In order to build a catalog of real co-orbital objects, we filtered the NASA Horizons asteroids catalog inside the maximum resonant width and analyzed which objects are indeed in resonance. In total, we found 167 objects to be in co-orbital resonant motion with Solar System planets excluding Jupiter. We were able to recover all the already known objects and to confirm the resonant state of some new ones. Mercury remains to be the only planet with zero known co-orbitals and no L5 Earth Trojan has been discovered so far. Among the interesting identified orbits we highlight the one of 2021 FV1, a new Mars Quasi-satellite. Despite having circulating critical angles, we found some objects to be dynamically driven by the resonance such as 2012 QR50.

This study focuses on transforming galaxy images between astronomical surveys, specifically enhancing images from the Sloan Digital Sky Survey (SDSS) and the Dark Energy Camera Legacy Survey (DECaLS) to achieve quality comparable to the Hyper Suprime-Cam survey (HSC). We proposed a hybrid model called Pix2WGAN, which integrates the pix2pix framework with the Wasserstein Generative Adversarial Network with Gradient Penalty (WGAN-GP) to convert low-quality observational images into high-quality counterparts. Our model successfully transformed DECaLS images into pseudo-HSC images, yielding impressive results and significantly enhancing the identification of complex structures, such as galaxy spiral arms and tidal tails, which may have been overlooked in the original DECaLS images. Moreover, Pix2WGAN effectively addresses issues like artifacts, noise, and blurriness in both source and target images. In addition to the basic Pix2WGAN model, we further developed an advanced architecture called Cascaded Pix2WGAN, which incorporates a multi-stage training mechanism designed to bridge the quality gap between SDSS and HSC images, demonstrating similarly promising outcomes. We systematically assessed the similarity between the model-generated pseudo-HSC images and actual HSC images using various metrics, including Mean Squared Error (MSE), Peak Signal-to-Noise Ratio (PSNR), and Structural Similarity Index (SSIM), along with perceptual metrics such as Learned Perceptual Image Patch Similarity (LPIPS) and Fréchet Inception Distance (FID). The results indicate that images transformed by our model outperform both the original SDSS and DECaLS images across nearly all evaluation metrics. Our research is expected to provide significant technical support for astronomical data analysis, cross-survey image integration, and high-precision astrometry.

The magnetic braking (MB) plays an important role in driving the evolution of low-mass X-ray binaries (LMXBs). The modified MB prescription, convection and rotation boosted (CARB) model, is very successful in reproducing the detected mass-transfer rates of persistent neutron star (NS) LMXBs. In this work, we investigate whether the CARB MB prescription could account for the formation and evolution of some NS and black hole (BH) LMXBs with an observed orbital period derivative. Using the MESA code, we perform a detailed binary evolution model for six NS and three BH LMXBs. Our simulations find that the CARB MB prescription can successfully reproduce the observed donor-star masses, orbital periods, and period derivatives of four NS LMXBs and one BH LMXB. Our calculated effective temperatures are in good agreement with the detected spectral types of two NS LMXBs and one BH LMXB. However, the standard MB model is difficult to produce the observed period derivatives of those LMXBs experiencing a rapid orbital shrinkage or expansion.

Investigating line-locked phenomena within quasars is crucial for understanding the dynamics of quasar outflows, the role of radiation pressure in astrophysical flows, and the star formation history and metallicity of the early universe. We have initiated the Tracking Outflow by Line-Locking (TOLL) project to study quasar outflow by studying line-locking signatures using high-resolution high signal-to-noise ratio quasar spectra. In this paper, we present a case study of the line-locking signatures from QSO J221531-174408. The spectrum was obtained using the Very Large Telescope-UV Visual Echelle Spectrograph. We first identify associated absorbers in the spectrum using CIV, NV, and Si IV doublets and measure their velocity shifts, covering fractions, and column densities through line profile fitting technique. Then we compare the velocity separations between different absorbers, and detect nine pairs of line-locked C IV doublets, three pairs of line-locked N V doublets, and one pair of line-locked SiIV doublets. This is one of the four quasars known to possess line-locked signatures in C IV, Si IV, and N V at the same time. We also find three complex line-locked systems, where three to five absorbers are locked together through multi-ion doublets. Our study suggests that line-locking is a common phenomenon in the quasar outflows, and theoretical models involving more than two clouds and one ionic doublet are needed in the future to explain the formation of these complex line-locking signatures.

Most research on astrophysical lensing has been conducted using the geometric optics framework, where there exists a clear concept of lensing images. However, wave optics effects can be important for coherent sources, e.g. pulsars, fast raio bursts, and gravitational waves observed at long wavelengths. There, the concept of lensing images needs an extension. We introduce the concept of the `lensing point-spread function' (LPSF), the smoothed flux density distribution of a coherent point source after being lensed, as a generalization of the lensing image concept at finite frequencies. The frequency-dependent LPSF captures the gradual change of the flux density distribution of the source from discrete geometric images at high frequencies to a smooth distribution at low frequencies. It complements other generalizations of lensing images, notably the imaginary images and the Lefschetz thimbles. Being a footprint of a lensing system, the LPSF is useful for theoretical studies of lensing. Using the LPSF, we identify a frequency range with non-trivial wave effects, where both geometric optics and perturbative wave optics fail, and determine this range to be $|\kappa|^{-1} \lesssim \nu \lesssim 10$, with $\kappa$ and $\nu$ being the dimensionless lens amplitude and the reduced observing frequency, respectively. Observation of LPSFs with non-trivial wave effects requires either very close-by lenses or very large observing wavelengths. The potential possibilities are the lensing of gravitational waves, the plasma lensing of Milky Way pulsars, and lensing by the solar gravitational lens.

Context. The PLAnetary Transits and Oscillations of stars (PLATO) mission will observe the same area of the sky continuously for at least two years in an effort to detect transit signals of an Earth-like planet orbiting a solar-like star. Aims. We aim to study how short-term solar-like variability caused by oscillations and granulation would affect PLATO's ability to detect and size Earth if PLATO were to observe the Solar System itself. Methods. We injected Earth-like transit signals onto real solar data taken by the Helioseismic and Magnetic Imager (HMI) instrument. We isolated short-term stellar variability by removing any variability with characteristic timescales longer than five hours. We then added a noise model for a variety of different stellar magnitudes computed by PlatoSim assuming an observation by all 24 normal cameras. We first compared four different commonly used treatments of correlated noise in the time domain. We then tried to recover pairs of transit signals. Finally, we performed transit fits using realistic priors on planetary and stellar parameters. Results. We find that short-term solar-like variability affects the correct retrieval of Earth-like transit signals in PLATO data. Variability models accounting for variations with typical timescales at the order of one hour are sufficient to mitigate these effects. For bright targets (8.5 - 10.5 mag), the transit signal of an Earth analogue can reliably be detected in PLATO data. For faint targets the results of transit search algorithms have to be verified by transit-fitting algorithms to avoid false positive detections being flagged. For bright targets (V-mag $\leq$ 9.5), the radius of an Earth-like planet orbiting a solar-like star can be correctly determined at a precision of 3% or less, assuming that at least two transit events are observed and the characteristics of the host star are well understood.

Building a comprehensive catalog of galaxy clusters is a fundamental task for the studies on the structure formation and galaxy evolution. In this paper, we present COSMIC (Cluster Optical Search using Machine Intelligence in Catalogs), an algorithm utilizing machine learning techniques to efficiently detect galaxy clusters. COSMIC involves two steps, including the identification of the brightest cluster galaxies and the estimation of the cluster richness. We train our models on the galaxy data from the Sloan Digital Sky Survey and WHL galaxy cluster catalog. Validated to a test data in the region of northern Galactic cap, COSMIC algorithm demonstrates a high completeness when cross-matching with previous cluster catalogs. Richness comparison with previous optical and X-ray measurements also demonstrated a tight correlation. Our methodology showcases robust performance in galaxy cluster detection and holds promising prospects for applications in upcoming large-scale surveys.

The neutron capture rates and temperature dependent stellar beta decay rates of Mo isotopes are investigated within the framework of the statistical code (Talys v1.96) and proton neutron quasi particle random phase approximation (pnQRPA) models. The Maxwellian average cross section (MACS) and neutron capture rates for 95 98Mo(n,{\gamma}) 96 99Mo radiative capture processes are analyzed within the framework of statistical code Talys v1.96 based on the phenomenological nuclear level density (NLD) parameters and gamma strength functions (GSFs). The current model based computed data for MACS provide a good comparison with the existing measured data. The temperature dependent stellar \b{eta} decay rates of 95 98Mo are investigated at different densities within the framework of the pn-QRPA model. For the considered temperature range, the neutron capture rates were found to be higher, by orders of magnitude, than the stellar \b{eta} rates.

The geometrically thick dusty torus structure is believed to exist in the nuclear region of galaxies (especially in active galactic nuclei, AGNs). The debris stream from a tidal disruption event (TDE) will possibly collide with the dusty torus and produce a transient flare. We perform three-dimensional hydrodynamic simulations to model the dynamical evolution of the interaction between unbound debris and dusty torus. During the continuous interaction, the shocked material will be spilled out from the interaction region and form an outflow. We calculate the temporal evolution of synchrotron emission by assuming that the shock accelerates a fraction of electrons in the outflow into a non-thermal distribution. We find that radio emission from the debris-torus collision generates a steep-rise and slow-decline radio light curve due to the sharp edge and dense gas of dusty torus, where the radio outburst delays the main optical/X-ray outburst by several years or even several tens of years. We apply our model to a TDE that happened in a narrow-line Seyfert I (PS16dtm), where both the radio spectrum and the light curve can be roughly reproduced. Future high-sensitivity, wide-field-of-view radio surveys have the opportunity to detect more such radio flares.

We aim to gain new insights into the controversial dynamical status of MACS J0329-0211 (MACS0329), a massive cluster at z=0.4503, with a new analysis using a large sample of member galaxies as kinematic tracers. Our analysis is based on extensive spectroscopic data for more than 1700 galaxies obtained with the VIMOS and MUSE spectrographs at the Very Large Telescope (VLT), in combination with B and Rc Suprime-Cam photometry from the Subaru archive. According to our member selection procedure, we define a sample of 430 MACS0329 galaxies within 6 Mpc, corresponding to about 3 times the virial radius. We estimate the global velocity dispersion, sigmaV=841 km/s, and present the velocity dispersion profile. We estimate a mass M200=9.2E14 in units of solar masses, using 227 galaxies within R200=1.71 Mpc, for which sigmaV200=1018 km/s. The spatial distribution of the red galaxies traces a SE-NW elongated structure, without signs of a velocity gradient. This structure likely originates from the main phase of cluster assembly. The distribution of the blue galaxies is less concentrated, more rounded and shows signs of substructure, all characteristics indicating a recent infall of groups from the field. We detect two loose clumps of blue galaxies in the south and southwest at a distance about R200 from the cluster center. The strong spatial segregation among galaxy populations is not accompanied by kinematical difference. Thanks to our extensive catalog of spectroscopic redshift, we are able to study galaxy systems that are intervening along the line of sight. We can identify two foreground galaxy systems (GrG1 at z=0.31 and GrG2 at z=0.38) and one background system (GrG3 at z=0.47). We point out that the second brightest galaxy projected onto the MACS0329 core is in fact the dominant galaxy of the foreground group GrG2.

Identifying the dominant ionizing sources in galaxies is essential for understanding their formation and evolution. Traditionally, spectra are classified based on their dominant ionizing source using strong emission lines and Baldwin, Phillips, \& Terlevich (BPT) diagrams. The ionizing source is traditionally determined by the emission line ratios using the BPT diagrams. Low-Ionization Nuclear Emission-line Regions (LINERs) are a class of ionizing mechanisms that is observationally identified but with a poorly understood origin, unlike the case of star forming regions and active galactic nuclei (AGN). LINERs, typically found in early-type galaxies, are often associated with low-luminosity AGN activity but may also be powered by aging stellar populations, particularly post-Asymptotic Giant Branch (p-AGB) stars. In this study, we employ a machine-learning-based encoder, Spender, to analyze the full MaNGA IFU spectra and identify key spectral features of LINERs. By examining the continuum and line emission of these spaxels, our approach aims to uncover hidden patterns and better understand the dominant ionizing sources. We show in this work that the neural network-based encoder was able identify LINER sources from the stellar continuum alone. The characteristics of the stellar population underlying LINER regions are consistent with evolved low mass stars implying that the source driving LINER emission is probably p-AGB stars rather than AGN activity.

With the arrival of the next generation of ultra-deep optical imaging surveys reaching $\mu_V$$\sim$30 mag/arcsec$^2$ (3$\sigma$; 10"$\times$10"), the removal of scattered light due to the point spread function (PSF) effect remains a critical step for the scientific exploitation of the low surface brightness information contained in these data. Because virtually all pixels in the ground-based images are affected by an unwanted screen of light with a brightness greater than $\mu_V$$\sim$29 mag/arcsec$^2$, the characterization of the extended PSF (R$>$5 arcmin) is mandatory. We describe the procedure used to construct the extended PSFs of the LIGHTS survey in the g- and r-band images taken with the Large Binocular Cameras (LBCs) of the Large Binocular Telescope (LBT). We produce PSFs with a radial extension of 7.5 arcmins. These are later extended to 30 arcmins following an empirically motivated power-law extrapolation of their behaviour in their outermost regions. As an example of the application of our methodology, we subtract the scattered light around the galaxy NGC3198. The result of this subtraction clearly shows the outermost parts of the galaxy's disc, which have been obscured by the influence of nearby bright stars. Following a reproducible science philosophy, we make all the PSF models and the scripts used to do the analysis of this paper publicly available.

We explore red stellar populations toward the W3 giant molecular cloud through the use of optical-to-infrared (IR) photometry and Gaia DR 3 data, simultaneously characterizing stellar content and properties of dust in the molecular medium. We use a Rayleigh-Jeans Color Excess (RJCE) method modified to de-redden stellar observations of both red giants (RGs) and OB stars, and construct an IR Hertzsprung-Russell diagram validated against the Besanccon Galactic model. Taking advantage of the near-universal IR interstellar extinction law and precise Gaia measurements, we develop a method for obtaining the spectral classification, foreground extinction, and distance moduli of stars, validated by spectroscopically-confirmed OB stars. We constrain the observed parallax and proper motion of OB stars in W3, demonstrating the importance of considering systematic effects in the parallax bias, and assign parallax- and proper motion-based cloud membership to our stellar samples. While it has been assumed that all spectroscopic OB stars are inside the W3 cloud, we find evidence of seven background B stars and three potential runaway OB stars. The methods developed here based on known stellar populations enable us to identify 82 new OB candidates that are confidently within the cloud. We quantify several dust-to-dust empirical correlations, in particular the IR color excess $E(H-[4.5])$ and the optical depth $\tau_1$ of submillimeter dust emission at 1 THz using RGs behind W3, measuring a best fit of $E(H-[4.5]) = (1.07 \pm 0.04) \times 10^3 \, \tau_{1,\,\mathrm{HOTT}} + (0.00 \pm 0.02)$ mags.

The number density and redshift evolution of optically selected galaxy clusters offer an independent measurement of the amplitude of matter fluctuations, $S_8$. However, recent results have shown that clusters chosen by the redMaPPer algorithm show richness-dependent biases that affect the weak lensing signals and number densities of clusters, increasing uncertainty in the cluster mass calibration and reducing their constraining power. In this work, we evaluate an alternative cluster proxy, outskirt stellar mass, $M_{\textrm{out}}$, defined as the total stellar mass within a $[50,100]$ kpc envelope centered on a massive galaxy. This proxy exhibits scatter comparable to redMaPPer richness, $\lambda$, but is less likely to be subject to projection effects. We compare the Dark Energy Survey Year 3 redMaPPer cluster catalog with a $M_{\textrm{out}}$ selected cluster sample from the Hyper-Suprime Camera survey. We use weak lensing measurements to quantify and compare the scatter of $M_{\textrm{out}}$ and $\lambda$ with halo mass. Our results show $M_{\textrm{out}}$ has a scatter consistent with $\lambda$, with a similar halo mass dependence, and that both proxies contain unique information about the underlying halo mass. We find $\lambda$-selected samples introduce features into the measured $\Delta \Sigma$ signal that are not well fit by a log-normal scatter only model, absent in $M_{\textrm{out}}$ selected samples. Our findings suggest that $M_{\textrm{out}}$ offers an alternative for cluster selection with more easily calibrated selection biases, at least at the generally lower richnesses probed here. Combining both proxies may yield a mass proxy with a lower scatter and more tractable selection biases, enabling the use of lower mass clusters in cosmology. Finally, we find the scatter and slope in the $\lambda-M_{\textrm{out}}$ scaling relation to be $0.49 \pm 0.02$ and $0.38 \pm 0.09$.

We present an analysis of 126 new radial velocity measurements from the MAROON-X spectrograph to investigate the TOI-1266 system, which hosts two transiting sub-Neptunes at 10.8 and 18.8 days. We measure masses of $M_{b}=4.01~\pm~0.55~M_{\oplus}$ for TOI-1266 b and $M_{c}=2.00~\pm~0.72~M_{\oplus}$ for TOI-1266 c. Our mass measurements agree with existing HARPS-N observations which we combined using a weighted average yielding masses for TOI-1266 b, and c of $M_{b}=4.10~\pm~0.43~M_{\oplus}$, $M_{c}=2.4~\pm~0.54~M_{\oplus}$ respectively. The combined dataset enabled a $\approx30\%$ improvement in mass precision. With bulk densities of $\rho_{b}$ = 1.25 $\pm$ 0.36 g cm$^{-3}$ and $\rho_{c}$ = 1.36 $\pm$ 0.31 g cm$^{-3}$, the planets are among the lowest density sub-Neptunes orbiting an M dwarf. They are both consistent with rocky cores surrounded by hydrogen helium envelopes. TOI-1266 c may also be consistent with a water-rich composition, but we disfavor that interpretation from an Occam's razor perspective.

The recent mass ($0.77 \pm ^{0.20}_{0.17}M_{\odot}$) and radius ($10.4\pm^{0.86}_{0.78} \text{km}$) measurement of HESS J1731-347 made it one of the most fascinating object if it is indeed a neutron star. In this work, we examine the current status of the dense matter equation of states in the context of this compact object being a neutron star. We use three sets of equation of states corresponding to the three classes - neutron stars, strange stars, and hybrid stars and perform Bayesian model selection on them. Our results show that for hadronic models, the EoS is preferred to be stiff at the intermediate densities. This makes the Brueckner-Hartree-Fock approximation and models based on Effective-interactions deviate from current astrophysical observations on the inclusion of HESS J1731-347. Furthermore, for the strange star family, the equation of states composed of three flavor quarks prefers relatively smaller bag parameters. Analyzing the hybrid family of equation of states consisting of a first-order phase transition revealed preferences for early first-order phase transition. Comparing all the preferred equations of state among each family, it was found that the current astrophysical constraints most prefer the hybrid equation of states.

The possibility that the vacuum energy density (VED) could be time dependent in the expanding Universe is intuitively more reasonable than just a rigid cosmological constant for the entire cosmic history. The framework of the running vacuum model (RVM) is a salient example derived from QFT in curved spacetime, wherein the VED appears as a power series of the Hubble rate, $H(t)$, and its derivatives. The RVM contributes to alleviate the cosmological tensions and at a more fundamental level it also helps to smooth out certain hardcore aspects of the cosmological constant problem. Composite dark energy (DE) extensions of the RVM are possible, in which the DE is a mixed fluid made out of running vacuum and an entity X called ``phantom matter'' which, I should stress, is radically different from phantom DE, since the former produces positive pressure like ordinary matter (therein its name). The prototype is the old $\Lambda$XCDM model from 2006. Recently a simplified version of the latter, called the $w$XCDM, has proven capable to yield a truly outstanding fit to the cosmological data with dramatic implications for mitigating the cosmological tensions. The analysis is sensitive to the kind of BAO data used, transversal (2D) or anisotropic (3D). For both BAO types the overall fit quality is substantially better than for the standard $\Lambda$CDM model, and the growth tension becomes alleviated. Nevertheless, only for 2D BAO the $H_0$-tension can be completely cut down. Most notably, this scenario favours (at $3.3\sigma$ c.l.) quintessence-like behavior around our time, as recently reported by the DESI Collaboration. In what follows I will summarize in very simple phenomenological terms the essentials of the RVM idea and its composite extensions, as well as the remarkable practical implications ensuing from such scenarios.

Observations of the Jovian upper atmosphere at high latitudes in the UV, IR and mm/sub-mm all indicate that the chemical distributions and thermal structure are broadly influenced by auroral particle precipitations. Mid-IR and UV observations have shown that several light hydrocarbons (up to 6 carbon atoms) have altered abundances near Jupiter's main auroral ovals. Ion-neutral reactions influence the hydrocarbon chemistry, with light hydrocarbons produced in the upper stratosphere, and heavier hydrocarbons as well as aerosols produced in the lower stratosphere. One consequence of the magnetosphere-ionosphere coupling is the existence of ionospheric jets that propagate into the neutral middle stratosphere, likely acting as a dynamical barrier to the aurora-produced species. As the ionospheric jets and the background atmosphere do not co-rotate at the same rate, this creates a complex system where chemistry and dynamics are intertwined. The ion-neutral reactions produce species with a spatial distribution following the SIII longitude system in the upper stratosphere. As these species sediment down to the lower stratosphere, and because of the progressive dynamical decoupling between the ionospheric flows and the background atmosphere, the spatial distribution of the auroral-related species progressively follows a zonal distribution with increasing pressures that ultimately produces a system of polar and subpolar hazes that extends down to the bottom of the stratosphere. This paper reviews the most recent work addressing different aspects of this environment.

The existence of a mass gap of $3-5\,M_{\odot}$ between the heaviest neutron stars (NSs) and the lightest black holes (BHs), inferred from the BH mass distribution in low mass X-ray binaries (LMXBs), has been suggested for decades. The recently reported gravitational-wave source GW230529 has been confidently identified as a NSBH merger, with the BH mass falling within this lower mass gap. This detection provides strong evidence against the existence of the latter and introduces new implications for the coalescing NSBH population, including a revised BH mass distribution and an updated local merger rate. In this study, we employ POSYDON, a binary population synthesis code that integrates detailed single- and binary-star models, to investigate coalescing NSBH binaries formed through isolated binary evolution. In particular, we focus on the BH mass distribution of the intrinsic NSBH merger population. We find that, with a high common-envelope efficiency of $\alpha_{\rm{CE}} =2 $, the BH masses in NSBH mergers concentrate in the lower mass gap, aligning more closely with observations. However, after accounting for the constraints of the selection bias against mass-gap BHs in LMXBs, which suggests that the maximum NS birth mass is below $\simeq 2\,M_{\odot}$, we find that introducing a high $\alpha_{\rm{CE}}$ is not required to match observations. Additionally, we explore the impact of core-collapse supernova kicks. Finally, we present the property distributions of observable NSBH mergers from our simulation and find that they match well with the observations. We find that the fraction of electromagnetic counterparts in observable populations is $\approx 4-30\%$, depending on different NS equations of state. Future detections of coalescing NSBH binaries would provide invaluable insights into SN mechanisms, common envelope evolution, and NS physics.

There are some traces of the existence of internal ocean in some icy moons, such as the vapor plumes of Europa and Enceladus. This implies a region of liquid water beneath the surface ice shell. Since liquid water would be essential for the origin of life, it is important to understand the development of these internal oceans, particularly their temperature distribution and evolution. The balance between tidal heating and radiative cooling is believed to sustain liquid water beneath an icy moon's surface. We aim to simulate the tidal heating of an internal ocean in an icy moon using 3-dimensional numerical fluid calculations with the Smoothed Particle Hydrodynamics (SPH) method. We incorporated viscosity and thermal conduction terms into the governing equations of SPH. However, we encountered two issues while calculating rigid body rotation using SPH with a viscous term: (1) conventional viscosity formulations generated unphysical forces that hindered rotation, and (2) there was artificial internal energy partitioning within the layered structure, which was due to the standard SPH formulations. To address the first issue, we modified the viscosity this http URL the second, we adopted Density Independent SPH (DISPH) developed in previous studies to improve behavior at discontinuous surfaces. Additionally, we implemented radiative cooling using an algorithm to define fluid surfaces via the particle method. We also introduced an equation of state accounting for phase transitions. With these modifications, we have refined the SPH method to encompass all necessary physical processes for simulating the evolution of icy moons with internal oceans.

Unimolecular gas phase chemical reactions could be activated by both infrared (IR) radiation and inter-molecular collision in the interstellar environment. Understanding the interplay and competition between the radiation and collision activation mechanisms is crucial for assessing accurate reaction rate constants with an appropriate model. In this work, guided by an extended version of Lindemann theory that considers the contribution of both the radiation-activation and collision-activation to the rate constant of unimolecular reactions, we show that the relative importance of the two mechanisms can be measured by a dimensionless number $PR$ that is the ratio of the collision frequency to the radiation absorption rate of the molecule. The reaction kinetic is dominated by collision-activation or by radiation activation depending on whether $PR$ is larger or smaller than a reference value $PR^*$, which is determined to be $PR^* \approx 10$ based on magnitudes of molecular properties and is verified by detailed calculations of a number of typical interstellar unimolecular reactions. This method of evaluating the relative importance of the two mechanisms is checked against master equation calculations of two interstellar reactions: the dissociation reaction of silicilic acid around the asymptotic giant branch (AGB) star and the methyl association in Titan's atmosphere, and the validity is verified. The method can be used in the future to help determine the appropriate and effective modeling approach for chemical reactions in astrophysical environments.

The Local Ultraviolet to Infrared Treasury (LUVIT) is a Hubble Space Telescope program that combines newly acquired data in the near ultraviolet (NUV), optical, and near infrared (NIR) with archival optical and NIR imaging to produce multiband panchromatic resolved stellar catalogs for 23 pointings in 22 low-mass, star-forming galaxies ranging in distance from the outskirts of the Local Group to ~3.8 Mpc. We describe the survey design, detail the LUVIT broadband filter observations and the archival datasets included in the LUVIT reductions, and summarize the simultaneous multiband data reduction steps. The spatial distributions and color-magnitude diagrams (CMDs) from the resulting stellar catalogs are presented for each target, from the NUV to the NIR. We demonstrate in which regions of the CMDs stars with NUV and optical, optical and NIR, and NUV through NIR detections reside. For each target, we use the results from artificial star tests to measure representative completeness, bias, and total photometric uncertainty as a function of magnitude in each broadband filter. We also assess which LUVIT targets have significant spatial variation in the fraction of stars recovered at a given magnitude. The panchromatic LUVIT stellar catalogs will provide a rich legacy dataset for a host of resolved stellar population studies.

Dust-obscured galaxies (DOGs) containing central supermassive black holes (SMBHs) that are rapidly accreting (i.e., having high Eddington ratios, $\lambda_\mathrm{Edd}$) may represent a key phase closest to the peak of both the black-hole and galaxy growth in the coevolution framework for SMBHs and galaxies. In this work, we present a 68 ks XMM-Newton observation of the high-$\lambda_\mathrm{Edd}$ DOG J1324+4501 at \mbox{$z \sim 0.8$}, which was initially observed by Chandra. We analyze the XMM-Newton spectra jointly with archival Chandra spectra. In performing a detailed \mbox{X-ray} spectral analysis, we find that the source is intrinsically \mbox{X-ray} luminous with $\log (L_\mathrm{X}$/erg s$^{-1}) = 44.71^{+0.08}_{-0.12}$ and heavily obscured with $\log (N_\mathrm{H}/\mathrm{cm}^{-2}) = 23.43^{+0.09}_{-0.13}$. We further utilize UV-to-IR archival photometry to measure and fit the source's spectral energy distribution (SED) to estimate its host-galaxy properties. We present a supplementary comparison sample of 21 \mbox{X-ray} luminous DOGs from the XMM-SERVS survey with sufficient ($> 200$) $0.5-10$ keV counts to perform a similarly detailed X-ray spectral analysis. Of the X-ray luminous DOGs in our sample, we find that J1324+4501 is the most remarkable, possessing one of the highest \mbox{X-ray} luminosities, column densities, and star-formation rates. We demonstrate that J1324+4501 is in an extreme evolutionary stage where SMBH accretion and galaxy growth are at their peaks.

High-cadence high-resolution spectroscopic observations of infant Type II supernovae (SNe) represent an exquisite probe of the atmospheres and winds of exploding red-supergiant (RSG) stars. Using radiation hydrodynamics and radiative transfer, we study the gas and radiation properties during and after the phase of shock breakout, considering RSG progenitors enshrouded within a circumstellar material (CSM) that varies in extent, density, and velocity profile. In all cases, the original, unadulterated CSM structure is probed only at the onset of shock breakout, visible in high-resolution spectra as narrow, often blueshifted emission, possibly with an absorption trough. As the SN luminosity rises during breakout, radiative acceleration of the unshocked CSM starts, leading to a broadening of the ``narrow'' lines by ~100 and up to ~1000km/s, depending on CSM properties. This acceleration is maximum close to the shock, where the radiative flux is greater, and thus typically masked by optical-depth effects. Generally, narrow-line broadening is greater for more compact, tenuous CSM because of the proximity to the shock where the flux is born, and smaller in denser and more extended CSM. Narrow-line emission should show a broadening that slowly increases first (the line forms further out in the original wind), then sharply rises (the line forms in a region that is radiatively accelerated), before decreasing until late times (the line forms further away in regions more weakly accelerated). Radiative acceleration should inhibit X-ray emission during the early, IIn phase. Although high spectral resolution is critical at the earliest times to probe the original slow wind, radiative acceleration and associated line broadening may be captured with medium resolution allowing a simultaneous view of narrow, Doppler-broadened as well as extended, electron-scattering broadened line emission.

We present PySCo, a fast and user-friendly Python library designed to run cosmological $N$-body simulations across various cosmological models, such as $\Lambda$CDM and $w_0w_a$CDM, and alternative theories of gravity, including $f(R)$, MOND and time-dependent gravitational constant parameterisations. PySCo employs Particle-Mesh solvers, using multigrid or Fast Fourier Transform (FFT) methods in their different variations. Additionally, PySCo can be easily integrated as an external library, providing utilities for particle and mesh computations. The library offers key features, including an initial condition generator based on up to third-order Lagrangian Perturbation Theory (LPT), power spectrum estimation, and computes the background and growth of density perturbations. In this paper, we detail PySCo's architecture and algorithms and conduct extensive comparisons with other codes and numerical methods. Our analysis shows that, with sufficient small-scale resolution, the power spectrum at redshift $z = 0$ remains independent of the initial redshift at the 0.1\% level for $z_{\rm ini} \geq$ 125, 30, and 10 when using first, second, and third-order LPT, respectively. Although the seven-point Laplacian method used in multigrid also leads to power suppression on small scales, this effect can largely be mitigated when computing ratios. In terms of performance, PySCo only requires approximately one CPU hour to complete a Newtonian simulation with $512^3$ particles (and an equal number of cells) on a laptop. Due to its speed and ease of use, PySCo is ideal for rapidly generating vast ensemble of simulations and exploring parameter spaces, allowing variations in gravity theories, dark energy models, and numerical approaches. This versatility makes PySCo a valuable tool for producing emulators, covariance matrices, or training datasets for machine learning.

We investigate the cyclotron resonant scattering features (CRSFs) of the accreting X-ray pulsar Cen X-3 and significantly detect the 29 keV cyclotron line features in the hard X-ray averaged spectroscopy studies based on the recent Insight-HXMT observations in 2022, when Cen X-3 has X-ray luminosity $L_{\rm X} > \sim 5 \times 10^{37}$ erg\ s$^{-1}$ in the bands of 2 -- 60 keV. We do not find a harmonic line in the average spectra based on different continuum models. We showed that the CRSF energies have no correlation with time or luminosity in the average spectra. In addition, by performing a pulse phase-dependent spectral analysis, we revealed the fundamental line with the centroid energy ranging from 25 to 29 keV with a high significance over the spin phases. The evolution of the cyclotron line centroid energies over pulse phase is similar to the shape of pulse profiles, illustrating a positive correlation between the energy of CRSFs and the pulse phase flux.

An accretion disk can be formed around a secondary star in a binary system when the primary companion leaves the Main sequence and starts to lose mass at an enhanced rate. We study the accretion disk evolution and planetary migration in wide binaries. We use a numerical model of a non-stationary alpha-disk with a variable mass inflow. We take into account that the low-mass disk has an extended region that is optically thin along the rotation axis. We consider irradiation by both the host star and the donor. The migration path of a planet in such a disk is determined by the migration rate varying during the disk evolution. Giants may open/close the gap several times over the disk lifetime. We identify the new type of migration specific to parts of the growing disk with a considerable radial gradient of an aspect ratio. Its rate is enclosed between the type 2 and the fast type 1 migration rates being determined by the ratio of time and radial derivatives of the disk aspect ratio. Rapid growth of the wind rate just before the envelope loss by the donor leads to the formation of a zone of decretion, which may lead to substantial outward migration. In binaries with an initial separation $a\lesssim 100$ AU migration becomes most efficient for planets with 60-80 Earth masses. This results in approaching a short distance from the host star where tidal forces become non-negligible. Less massive Neptune-like planets at the initial orbits $r_\mathrm{p} \lesssim 2$ AU can reach these internal parts in binaries with $a \lesssim 30$ AU. In binaries, mass loss by the primary component at late evolutionary stages can significantly modify the structure of a planetary system around the secondary component resulting in mergers of relatively massive planets with a host star.

The origin of magnetic fields in white dwarfs (WDs) remains mysterious. Magnetic WDs are traditionally associated with field strengths $\gtrsim1\,\mathrm{MG}$, set by the sensitivity of typical spectroscopic magnetic field measurements. Informed by recent developments in red giant magnetoasteroseismology, we revisit the use of WD pulsations as a seismic magnetometer. WD pulsations primarily probe near-surface magnetic fields, whose effect on oscillation mode frequencies is to asymmetrize rotational multiplets and, if strong enough, suppress gravity-mode propagation altogether. The sensitivity of seismology to magnetic fields increases strongly with mode period and decreases quickly with the depth of the partial ionization-driven surface convective zone. We place upper limits for magnetic fields in $24$ pulsating WDs: $20$ hydrogen-atmosphere (DAV) and three helium-atmosphere (DBV) carbon-oxygen WDs, and one extremely low-mass (helium-core) pulsator. These bounds are typically $\sim1$-$10\,\mathrm{kG}$, although they can reach down to $\sim10$-$100\,\mathrm{G}$ for DAVs and helium-core WDs in which lower-frequency modes are excited. Seismic magnetometry may enable new insights into the formation and evolution of WD magnetism.

Synchrotron X-ray emission has been detected from nearly a dozen young supernova remnants (SNRs). X-rays of synchrotron origin exhibit linear polarization in a regular, non-randomly oriented magnetic field. The significant polarized X-ray emission from four such SNRs has already been reported on the basis of observations with the Imaging X-ray Polarimetry Explorer (IXPE). The magnetic-field structure as derived from IXPE observations is radial for Cassiopeia A, Tycho's SNR, and SN 1006, and tangential for RX J1713.7-3946. The latter together with the recent detection of a tangential magnetic field in SNR 1E 0102.2-7219 by the Australia Telescope Compact Array in the radio band shows that tangential magnetic fields can also be present in young SNRs. Thus, the dichotomy in polarization between young and middle-aged SNRs (radial magnetic fields in young SNRs, but tangential magnetic fields in middle-aged SNRs), previously noticed in the radio band, deserves additional attention. The present analysis of IXPE observations determines, for the first time, a magnetic-field structure in the northwestern rim of Vela Jr, also known as RX J0852.0-4622, and provides a new example of a young SNR with a tangential magnetic field.

The new generation of upcoming deep photometric and spectroscopic surveys will allow us to measure the astrophysical properties of faint galaxies in massive clusters. This would demand to produce simulations of galaxy clusters with better mass resolution than the ones available today if we want to make comparisons between the upcoming observations and predictions of cosmological models. But producing full-physics hydrodynamical simulations of the most massive clusters is not an easy task. This would involve billions of computational elements to reliably resolve low mass galaxies similar to those measured in observations. On the other hand, dark matter only simulations of cluster size halos can be done with much larger mass resolution but at the cost of having to apply a model that populate galaxies within each of the subhalos in these simulations. In this paper we present the results of a new set of dark matter only simulations with different mass resolutions within the THE THREE HUNDRED project. We have generated catalogs of galaxies with stellar and luminosity properties by applying the SAGE Semi-Analytical Model of galaxy formation. To obtain the catalogs consistent with the results from hydrodynamical simulations, the internal physical parameters of SAGE were calibrated with the Particle Swarm Optimization method using a subset of full-physics runs with the same mass resolution than the dark matter only ones.

We report the first observation of the vertical and temporal structure of the H3+ emission at the auroral footprint of Io, as observed by Juno/JIRAM. The brightness vertical profile shows a maximum at 600 km above 1 bar, with no apparent difference between the Main Alfvén Wing spot emission and the tail of the footprint. This observation is more compatible with a broadband energy distribution of the precipitating electrons, than a monoenergetic one. The temporal profile of H3+ column density has been observed after the passage of the MAW and shows a hyperbolic decrease. A model of H3+ decay is proposed, which takes into account the second-order kinetic of dissociative recombination of H3+ ions with electrons. The model is found to be in very good agreement with Juno observation. The conversion factor from radiance to column density has been derived, as well as the half-life for H3+, which is not constant but inversely proportional to the H3+ column density. This explains the wide range of H3+ lifetimes proposed before.

This study examines the relationship between early solar coronal mass ejection (CME) propagation, the associated filament eruption, and coronal dimming in the rare event observed on March 28, 2022, which featured a three-part CME in the low corona of active region AR 12975, including a bright core/filament, dark cavity, and bright front edge. We employ 3D filament and CME shock reconstructions using data from SolO, STEREO-A, and SDO to track the filament's path, height, and kinematics. Our analysis across three viewpoints shows the outer front in SolO/EUI 304 Å aligns with shock structures in STEREO-A/EUVI 195 Å, showing a full 3D EUV wave dome, later matching the outer CME front in STEREO-A COR2. We introduce the method ATLAS-3D (Advanced Technique for single Line-of-sight Acquisition of Structures in 3D) and validate it against traditional approaches to reconstruct CME shock using SOLO data exclusively. Additionally, we estimate early CME propagation characteristics based on coronal dimming evolution with the DIRECD method. Results indicate that the filament height increased from 28 to 616 Mm (0.04 to 0.89 Rs) within 30 minutes (11:05 to 11:35 UT), reaching peak velocity of around 648 km/s and acceleration of around 1624 m/s$^2$. At 11:45 UT, the filament deflected by 12° to a height of 841 Mm (1.21 Rs), while the CME shock expanded from 383 to 837 Mm (0.55 to 1.2 Rs) over 10 minutes. Key parameters include a CME direction inclined by 6°, a 21° half-width, and a 1.12 Rs cone height at the dimming's impulsive phase end. This event demonstrates that expanding dimming correlates with early CME development, with the DIRECD method linking 2D dimming to 3D CME evolution. These insights underscore the value of multi-viewpoint observations and advanced reconstructions for improving space weather forecasting.

Planet engulfment has been identified as one of the mechanisms for enhancing lithium abundance in stars. However, comprehensive investigations into lithium signatures following such events remain limited. Stars born together, sharing a common origin and stellar characteristics, provide a unique opportunity to study these signatures and compare lithium abundances. We demonstrate that the distinctive signature of planet engulfment in lithium abundance is only discernible among highly similar stellar twins. We present lithium abundance measurements for 125 co-moving pairs of stars, representing the largest sample to date with a single, homogeneous assessment of high-precision lithium abundance. While lithium abundance enhancements in pairs showing planet engulfment signatures are within 0.35 dex, we find that even at fixed stellar parameters (temperature and age), the intrinsic scatter in lithium abundance is typically 0.35 dex for G/F dwarfs and can be as large as 0.6 dex for older and cooler stars due to internal stellar evolution processes. Since the planet engulfment signature from lithium can be masked by stellar intrinsic scatter, our findings raise questions about relying solely on lithium as an indicator for planet engulfment or attributing lithium-richness in stars primarily to planet engulfment events.

The James Webb Space Telescope (JWST) has uncovered an abundant population of compact, extremely red, and X-ray weak objects at $z\gtrsim4$, knows as ``Little Red Dots" (LRDs). These objects exhibit spectral energy distributions that resemble both active galactic nuclei (AGN) and stellar population templates. However, whether dominated by AGN activity or compact star formation, the high redshifts and masses/luminosities of LRDs, coupled with their significant abundance, present potential challenges to the standard $\Lambda$CDM model. In this work, we proposes a novel cosmic interpretation of this anomaly, suggesting that these LRDs are likely massive galaxies seeded by supermassive primordial black holes (SMPBHs) came into being in the very early universe. We analyze 434 known LRDs from the 0.54 ${\rm deg}^2$ COSMOS-Web survey and test the hypothesis that they originated from SMPBHs assuming sub-Eddington accretion. According to our result, SMPBHs actually could lead to the existence of more LRDs, even at higher redshifts ($z>8$).

We investigate the influence of inelastic neutrino microphysics in general-relativistic magnetohydrodynamics simulations of a hypermassive neutron star. In particular, we include species/energy groups coupled neutrino-matter interactions, such as inelastic neutrino-electron scattering and electron-positron annihilation kernels, into simulations up to 50 ms. Neutrino-electron inelastic scattering is known to have effective neutrino-matter energy exchange. We show that, with neutrino-electron inelastic scattering, simulations predict 75% higher disc mass with slightly different mass-averaged compositions, and 18% more ejected mass with similar distributions. The enhancement of the mass of the disc and the ejecta results in stronger baryon pollution, leading to less favourable jet launching environments. Furthermore, neutrino luminosities are about 50, 40, and 30% higher for electron neutrino, electron anti-neutrino, and heavy-lepton neutrinos. In contrast, we do not see any significant impacts due to electron-positron annihilation.

In this paper we present a statistical study of the [C II] 158 $\mu$m line and the CO(1-0) emission for a sample of $\sim$200 local and high-$z$ (32 sources with $z>1$) galaxies with much different physical conditions. We explore the correlation between the luminosities of [C II] and CO(1-0) lines, and obtain a strong linear relationship, confirming that [C II] is able to trace total molecular gas mass, with a small difference between (U)LIRGs and less-luminous galaxies. The tight and linear relation between [C II] and CO(1-0) is likely determined by the average value of the observed visual extinction $A_V$ and the range of $G_0/n$ in galaxies. Further investigations into the dependence of $L_{\mathrm{[CII]}}/L_{\mathrm{CO(1-0)}}$ on different physical properties show that $L_{\mathrm{[CII]}}/L_{\mathrm{CO(1-0)}}$ (1) anti-correlates with $\Sigma_{\mathrm{IR}}$, and the correlation becomes steeper when $\Sigma_{\mathrm{IR}} \gtrsim 10^{11}$ $L_\odot\,\mathrm{kpc}^{-2}$; (2) correlates positively with the distance from the main sequence $\Delta(\mathrm{MS})$ when $\Delta(\mathrm{MS})\lesssim 0$; and (3) tends to show a systematically smaller value in systems where the [C II] emission is dominated by ionized gas. Our results imply that caution needs to be taken when applying a constant [[C II]-to-$M_{\mathrm{H_2}}$ conversion factor to estimate the molecular gas content in extreme cases, such as galaxies having low-level star formation activity or high SFR surface density.

Here we report our study on a low-mass ratio contact binary system NSVS 3198272. We subjected our multi-filter ground-based photometry to the light curve numerical modeling using a modified Wilson-Devinney code. We present three scenarios fitting to the data best: a simple, non-spotted model, a model with a circumpolar spot and a model with a marginal third light. The first two models return similar physical properties which are comparable to the previously reported results on this object. In discussion section we are arguing for advocating the model with a circumpolar spot. At the and we are considering the impact of specific scenarios on the calculated physical parameters of contact binaries. In addition, we report a discovery of a new bright variable star, found by chance as a field star.

We present a joint analysis of high-resolution CO(2-1) and Paschen-$\alpha$ emission lines to trace gas dynamics and spatially resolved star formation in ASPECS-LP.3mm.06, a $z=1.1$ main sequence galaxy. Utilizing data from the ALMA and JWST NIRCam Wide Field Slitless Spectroscopy (WFSS), we explore both ionized gas and molecular gas within this galaxy. With a substantial molecular gas fraction (f$_\mathrm{mol}$ = 0.44 $\pm$ 0.02), ASPECS-LP.3mm.06 remains on the star-forming main sequence and adheres to the Kennicutt-Schmidt (KS) relation, indicating typical gas-to-star conversion efficiency. Our analysis reveals extended structures across multiple wavelengths, suggesting regulated star formation within a stable disk. The spatially resolved star formation efficiency (SFE) and kinematic analysis indicate that ASPECS-LP.3mm.06 features a smooth mass assembly process across bulge and disk. Additionally, the galaxy exhibits modest dust extinction (A$_\mathrm{V}$ = 0.8), potentially linked to self-regulation during bulge formation. These findings position ASPECS-LP.3mm.06 as a prototypical galaxy, offering valuable insights into the mechanisms governing normal disk galaxy growth at z$\sim$1.

Gravitational waves (GWs) from stellar binary black hole (sBBH) mergers can be strongly gravitational lensed by intervening galaxies/galaxy clusters. Only a few works investigated the cluster-lensed sBBH mergers by adopting oversimplified models, while galaxy-lensed ones were intensively studied. In this paper, we estimate the detection rate of cluuster-lensed sBBH mergers with the third-generation GW detectors and its dependence on the lens models. We adopt detailed modeling of galaxy cluster lenses by using the mock clusters in the Synthetic Sky Catalog for Dark Energy Science with LSST (CosmoDC2) and/or approximations of the pseudo-Jaffe profile or an eccentric Navarro-Frenk-White dark matter halo plus a bright central galaxy with singular isothermal sphere profile. Considering the formation of sBBH mergers dominates by the channel of evolution of massive binary stars (EMBS), we find that the detection rate of cluster-lensed sBBHs is $\sim5-84$ yr$^{-1}$, depending on the adopted lens model and uncertainty in the merger rate density, and it is about $\sim{13_{-2.0}^{+28}}$yr$^{-1}$ if adopting relatively more realistic galaxy clusters with central main and member galaxies in the CosmoDC2 catalog, close to the estimated detection rate of sBBH mergers lensed by galaxies. In addition, we also consider the case that the production of sBBH mergers dominated by the dynamical interactions in dense stellar systems. We find that the detection rate of cluster-lensed sBBHs if from the dynamical channel is about $1.5$ times larger than that from the EMBS channel and the redshift distribution of former peaking at a higher redshift ($\sim3$) compared with that from latter ($\sim2$).

Novalike variables are a subgroup of cataclysmic variables (CVs) that -- unlike dwarf novae -- do not exhibit strong brightenings in their long-term light curves. Variations over time scales of weeks, months or years are mostly restricted to irregular low-amplitude modulations. However, some of them occasionally suffer from so-called stunted outbursts, that is, small-scale brightenings of less than a magnitude lasting for a couple of days to weeks. There is no consensus about the physical mechanisms behind these outbursts. Here I discuss the common properties of a group of novalike variables (which I call AH~Pictoris stars after its most prominent member) that exhibit a continuous train of successive stunted outbursts over their entire observational history, or at least for several years. The outburst amplitudes are stable in a given system, always ranging between 0.5 and 1~mag in the visual band. The outburst intervals, at an overall range between 12 and 30 days, and the outburst profiles can gradually evolve, but no sudden changes are observed. On shorter time scales the orbital waveforms are not only surprisingly similar, but also evolve in the same way over the outburst cycle. All AH~Pic stars have absolute visual magnitudes in the overlap region between absolute magnitudes of all novalike variables and of quiescent dwarf novae above the CV period gap. So far, I identified seven novalike variables with the consistent photometric behavior which may be termed the AH~Pic syndrome. Several more systems may be related objects. The relationship to the anomalous Z~Cam stars is discussed.

The magnetic orientation of coronal mass ejections (CMEs) is of great importance to understand their space weather effects. Although many evidences suggest that CMEs can undergo significant rotation during the early phases of evolution in the solar corona, there are few reports that CMEs rotate in the interplanetary space. In this work, we use multi-spacecraft observations and a numerical simulation starting from the lower corona close to the solar surface to understand the CME event on 2021 December 4, with an emphatic investigation of its rotation. This event is observed as a partial halo CME from the back side of the Sun by coronagraphs, and reaches the BepiColombo spacecraft and the MAVEN/Tianwen-1 as a magnetic flux rope-like structure. The simulation discloses that in the solar corona the CME is approximately a translational motion, while the interplanetary propagation process evidences a gradual change of axis orientation of the CME's flux rope-like structure. It is also found that the downside and the right flank of the CME moves with the fast solar wind, and the upside does in the slow-speed stream. The different parts of the CME with different speeds generate the nonidentical displacements of its magnetic structure, resulting in the rotation of the CME in the interplanetary space. Furthermore, at the right flank of the CME exists a corotating interaction region (CIR), which makes the orientation of the CME alter, and also deviates from its route due to the CME. These results provide new insight on interpreting CMEs' dynamics and structures during their travelling through the heliosphere.

Understanding the complex ionization structure and chemical composition of \hii\ regions poses a significant challenge in astrophysics. The abundance discrepancy problem, characterized by inconsistencies between abundances derived from recombination lines (RLs) and collisionally excited lines (CELs), has long been a puzzle in the field. In this theoretical study, we present novel photoionization models that incorporate temperature, density, and chemical inhomogeneities within a single cloud to comprehensively address this discrepancy. By accounting for the intricate interplay between ionization, excitation, and chemistry, our models successfully reproduce both observed RLs and CELs with with an average difference between our models and the observations of 25% -- within uncertainties inherent in Galactic archival long-slit and new SDSS-V Local Volume Mapper observations. Through comparisons between generic inhomogeneous model predictions and observations, demonstrating the ability of our theoretical framework to analyze the abundance discrepancy problem within \hii\ regions. Our results highlight the importance of incorporating spatially resolved temperature, density, and chemical structures when interpreting the physical processes governing emission line spectra in these astrophysical environments.

We present synthetic light curves and spectra from three-dimensional (3D) Monte Carlo radiative transfer simulations based on a 3D core-collapse supernova explosion model of an ultra-stripped $3.5\,\mathrm{M}_{\odot}$ progenitor. Our calculations predict a fast and faint transient with $\Delta m_{15} \sim 1\texttt{-} 2\,\mathrm{mag}$ and peak bolometric luminosity between $-15.3\,\mathrm{mag}$ and $-16.4\,\mathrm{mag}$. Due to a large-scale unipolar asymmetry in the distribution of $^{56}\mathrm{Ni}$, there is a pronounced viewing-angle dependence with about $1\,\mathrm{mag}$ difference between the directions of highest and lowest luminosity. The predicted spectra for this rare class of explosions do not yet match any observed counterpart. They are dominated by prominent Mg~II lines, but features from O, C, Si, and Ca are also found. In particular, the O~I line at \wl{7}{774} appears as a blended feature together with Mg~II emission. Our model is not only faster and fainter than the observed Ib/c supernova population, but also shows a correlation between higher peak luminosity and larger $\Delta m_{15}$ that is not present in observational samples. A possible explanation is that the unusually small ejecta mass of our model accentuates the viewing-angle dependence of the photometry. We suggest that the viewing-angle dependence of the photometry may be used to constrain asymmetries in explosion models of more typical stripped-envelope supernova progenitors in future.

We investigated the pulsating behavior of KIC 10855535 using Kepler 4-year long cadence data. Two independent frequencies were detected: a pulsation frequency F0 = 17.733260(5)d-1 and a low frequency f8=0.412643(8)d-1 We identify F0 as the fundamental frequency, at which a equidistant quintuplet is centered, suggesting that the star orbits in a binary system. The fitted orbital parameters align well with those reported in previous literature. Long-term phase modulation caused by binarity has been confirmed by considering TESS light curve. Through adjusting light time via removing the light time effect, we measured a linear change in period of order $\dot{P}/P \simeq 1.44\times 10^{-7}yr^{-1}$, a value that could be indicative of stellar evolution. The star also exhibits a gradual and stable amplitude growth, thereby raising the possibility of structural changes during its evolution. We attributed f8 and its two harmonics to rotation and surface spots, with further analysis suggesting evolving characteristics over time. Based on the hypothesis, KIC 10855535 may rotate slowly for its type, with a speed of 37(2)km/s. Overall, KIC 10855535 presents an exceptionally clean spectrum and a relatively slow rotation as a {\delta} Sct pulsator, exhibiting a single pulsation mode that undergoes both amplitude and phase modulation.

In 2023, the 12th edition of Global Trajectory Competition was organised around the problem referred to as "Sustainable Asteroid Mining". This paper reports the developments that led to the solution proposed by ESA's Advanced Concepts Team. Beyond the fact that the proposed approach failed to rank higher than fourth in the final competition leader-board, several innovative fundamental methodologies were developed which have a broader application. In particular, new methods based on machine learning as well as on manipulating the fundamental laws of astrodynamics were developed and able to fill with remarkable accuracy the gap between full low-thrust trajectories and their representation as impulsive Lambert transfers. A novel technique was devised to formulate the challenge of optimal subset selection from a repository of pre-existing optimal mining trajectories as an integer linear programming problem. Finally, the fundamental problem of searching for single optimal mining trajectories (mining and collecting all resources), albeit ignoring the possibility of having intra-ship collaboration and thus sub-optimal in the case of the GTOC12 problem, was efficiently solved by means of a novel search based on a look-ahead score and thus making sure to select asteroids that had chances to be re-visited later on.

Powerful, galactic outflows driven by Active Galactic Nuclei (AGNs) are commonly considered as a main mechanism to regulate star formation in massive galaxies. Ultra- and hyper-luminous IR galaxies (U/HyLIRGs) are thought to represent a transition phase of galaxies from a rapidly growing period to a quiescent status as gas swept out by outflows, providing a laboratory to investigate outflows and their feedback effects on the hosts. In this paper we report recent Gemini and ALMA observations of a HyLIRG, J1126 at $z=0.46842$, which has been identified with a puzzling co-existence of a fast ionized outflow ($>2000$ km s$^{-1}$) and an intense starburst (star formation rate of 800 $M_{\odot}$ yr$^{-1}$). The Gemini observation shows the fast ionized outflow is extended to several kpc with a mass-loss rate of 180 $M_{\odot}$ yr$^{-1}$. A massive molecular outflow with a high mass-loss rate (2500 $M_{\odot}$ yr$^{-1}$) is revealed by ALMA. The multi-phase outflows show large factors of momentum boost and loading of kinetic power, indicating a driving by thermal pressure of a nuclear hot wind and/or radiation pressure of a highly obscured AGN. In addition to ejection of kinetic energy, it is also found that the powerful outflow can induce an ionizing shock in the galaxy disk and enhance the excitation and dissociation of molecular gas. The powerful outflow probably results in an instantaneous negative feedback and shows potential to regulate the host growth in a long term.

The solar magnetized corona is responsible for various manifestations with a space weather impact, such as flares, coronal mass ejections (CMEs) and, naturally, the solar wind. Modeling the corona's dynamics and evolution is therefore critical for improving our ability to predict space weather In this work, we demonstrate that generative deep learning methods, such as Denoising Diffusion Probabilistic Models (DDPM), can be successfully applied to simulate future evolutions of the corona as observed in Extreme Ultraviolet (EUV) wavelengths. Our model takes a 12-hour video of an Active Region (AR) as input and simulate the potential evolution of the AR over the subsequent 12 hours, with a time-resolution of two hours. We propose a light UNet backbone architecture adapted to our problem by adding 1D temporal convolutions after each classical 2D spatial ones, and spatio-temporal attention in the bottleneck part. The model not only produce visually realistic outputs but also captures the inherent stochasticity of the system's evolution. Notably, the simulations enable the generation of reliable confidence intervals for key predictive metrics such as the EUV peak flux and fluence of the ARs, paving the way for probabilistic and interpretable space weather forecasting. Future studies will focus on shorter forecasting horizons with increased spatial and temporal resolution, aiming at reducing the uncertainty of the simulations and providing practical applications for space weather forecasting. The code used for this study is available at the following link: this https URL

Millimetron is a space observatory for millimeter and sub-millimeter observations planned for launch around 2030. The 10-meter diameter space unfolded telescope will be cooled down to 10~K and operated in the vicinity of Lagrange point L2. Mission lifetime is 10 years and it includes astronomical observations in two modes: as a space-ground interferometer and as a single-dish telescope. This paper presents the results of the Millimetron space observatory orbit design that takes into account technical and scientific requirements and constraints. It covers the computation of suitable operational orbits, the selection of an orbit, and the transfer from Earth. Furthermore, scientific objectives are demonstrated through VLBI visibility simulation and image reconstruction. Unlike previous works that used analytical methods, this work employs numerical integration for orbital design. Based on the orbital design methods developed in this work, we calculate the exact dates for departure, halo formation, and trajectory correction. Additionally, we investigate the existence of the short baseline projection and its specific dates for VLBI mode and show the feasibility of scientific objectives through VLBI visibility simulation and image reconstruction.

The redshift and size distributions of galaxy scale strong lenses depend on the evolution of early-type galaxies (ETGs) in the redshift range 0.2<z<1. We use this dependence to constrain the velocity dispersion function (VDF) evolution from the Strong Lensing Legacy Survey (SL2S) sample of lenses. Our modeling of the lens population includes lens identifiability given survey parameters, and constrains the evolution of the VDF based on the redshift distributions of sources and lenses as well as the distribution of Einstein radii. We consider five different assumptions for the reference VDF at redshift zero and two sets of scaling relations for the VDF. We find that in all cases the observed lens sample favors a slow evolution of both the VDF normalization factor and the VDF characteristic velocity with redshift which is consistent with a VDF that is constant in redshift for z<1.

Relativistic magnetic reconnection is one of the most fundamental mechanisms considered responsible for the acceleration of relativistic particles in astrophysical jets and magnetospheres of compact objects. Understanding the properties of the dissipation of magnetic fields and the formation of non-ideal electric fields is of paramount importance to quantify the efficiency of reconnection at energizing charged particles. Recent results from particle-in-cell (PIC) simulations suggest that the fundamental properties of how magnetic fields dissipate in a current sheet might be captured by an ``effective resistivity'' formulation, which would locally enhance the amount of magnetic energy dissipated and favor the onset of fast reconnection. Our goal is to assess this ansatz quantitatively by comparing fluid models of magnetic reconnection with a non-constant magnetic diffusivity and fully-kinetic models. We perform 2D resistive relativistic magnetohydrodynamic (ResRMHD) simulations of magnetic reconnection combined to PIC simulations using the same initial conditions (namely a Harris current sheet). We explore the impact of crucial parameters such as the plasma magnetization, its mass density, the grid resolution, and the characteristic plasma skin depth. Our ResRMHD models with effective resistivity can quantitatively reproduce the dynamics of fully-kinetic models of relativistic magnetic reconnection. In particular, they lead to reconnection rates consistent with PIC simulations, while for constant-resistivity fluid models the reconnection dynamics is generally 10 times slower. Even at modest resolutions the adoption of an effective resistivity can qualitatively capture the properties of kinetic reconnection models and produce reconnection rates compatible with collisionless models, i.e. of the order of $\sim10^{-1}$.

The increasing congestion in the near-Earth space environment has amplified the need for robust and efficient conjunction analysis techniques including the computation of the minimum distance between orbital paths in the presence of perturbations. After showing that classical Minimum Orbit Intersection Distance (MOID) computation schemes are unsuitable to treat Earth orbiting objects, the article presents an analytical approach to provide a more accurate estimate of the true distance between perturbed trajectories by incorporating the effect of zonal harmonics of arbitrary order. Cook's linear secular theory for the motion of the eccentricity vector is extended to include higher order eccentricity effects and applied to the computation of the minimum and maximum radii attained by two orbits at their mutual nodes, which can be employed to estimate the true distance between the two orbital paths and to establish an efficient algorithm for determining or excluding potential conjunctions. Extensive testing and validation are conducted using a high-fidelity propagator and a comprehensive dataset of resident space objects. The results demonstrate an accuracy below the km level for the orbit distance computation in 99\% of cases, which enables high-efficiency conjunction filtering.

Tidal forces in close binary systems have diverse impacts on magnetic activity. The synchronicity characteristic of close systems counteracts magnetic braking, thereby sustaining rapid rotation-a key factor in increased levels of magnetic activity. Tidal effects can also work against the slowing down of rotation during stellar evolution, when the star inflates into a red giant. A notable manifestation of the effect of binarity on activity in such systems is the appearance of active longitudes, which are thought to arise from the excitation of non-axisymmetric dynamo modes. Through some recent examples, the dynamo operation in RS CVn and BY Dra type systems is briefly reviewed in terms of spot cycles, active longitudes, flare activity, and differential rotation.

This study continues our investigation of early-type binaries using high-precision broad-band polarimetry, focusing on HD 165052, a massive O+O-type binary in the young cluster NGC 6530. Our aim was to monitor linear polarization variations and independently determine the orbital period through polarization data analysis. By examining the phase-locked Stokes parameters, we estimated the orbital inclination at about 55+5/-55 deg, and the orientation at 148+20/-22 deg. The binary's rotation direction is clockwise. Using the Dipol-2 polarimeter with the 60 cm KVA and Tohoku T60 telescopes, we achieved polarization measurement accuracy around 0.01 percent in the B, V, and R passbands. Period analysis identified a strong periodic signal with a cycle of about 1.48 days, corresponding to half of the orbital period of about 2.96 days. The second Fourier harmonics suggest a symmetric scattering geometry, likely due to electron scattering in interacting stellar winds. Our estimates indicate a mass-loss rate of approximately 4.0 x 10^(-7) solar masses per year. Observations of neighboring stars in NGC 6530 revealed complex interstellar polarization behaviors within the cluster, enhancing our understanding of its polarimetric environment.

The evolution of protoplanetary disks in regions with massive OB stars is influenced by externally driven winds that deplete the outer parts of disks. These winds have previously been studied via forbidden oxygen emission lines, which also arise in isolated disks in low-mass star forming-regions (SFRs) with weak external UV fields in photoevaporative or magnetic (internal) disk winds. It is crucial to determine how to disentangle external winds from internal ones. Here, we report a proxy for unambiguously identifying externally driven winds with a forbidden line of neutral atomic carbon, [C i] 8727 A. We compare for the first time the spatial location of the emission in the [O i] 5577 A, [O i] 6300 A, and [C i] 8727 A lines traced by VLT/MUSE-NFM, with the ALMA Band 7 continuum disk emission in a sample of 12 proplyds in the Orion Nebula Cluster (ONC). We confirm that the [O i] 5577 A emission is co-spatial with the disk emission, whereas the [O i] 6300 A is emitted both on the disk surface and on the ionization front of the proplyds. We show for the first time that the [C i] 8727 A line is also co-spatial with the disk surface in proplyds, as seen in the MUSE and ALMA data comparison. To verify whether the [C i] 8727 A line is detected in regions where external photoevaporation is not expected, we examine VLT/X-Shooter spectra for young stars in low-mass SFRs. Although the [O i] lines are well detected in all these targets, there is <<10% detection rate in the case of the [C i] 8727 A line. This number increases substantially to a ~40% detection rate in sigma-Orionis, a region with intermediate UV radiation. The spatial location of the [C i] 8727 A line emission and the lack of its detection in isolated disks in low-mass SFRs strongly suggest that this line is a tell-trace tracer of externally driven photoevaporative winds, which agrees with recent excitation models.

Synchrotron maser emission is a leading candidate to explain the coherent emission from Fast Radio Bursts (FRBs). This mechanism requires a population inversion in order to operate. We show that nonresonant interactions between Alfvén waves and a relativistic plasma result in the formation of population inversions across a wide range of magnetizations, $\sigma\gtrsim10^{-4}$, and temperatures, $10^{-2} \leq k_bT/mc^2 \leq 3$, spanning the parameters expected in FRB environments. We calculate the fraction of energy contained in the inversion across the whole of this parameter space for the first time and we show that energy fractions of $f_{inv}\sim 0.3$ are achieved for high magnetizations $\sigma >1$. The population inversion forms on time-scales compatible with the typical dynamical time-scales of magnetars for all magnetizations. Furthermore, we provide physical explanations for the behaviour of the interaction in different magnetization regimes, and identify the important characteristic values at which this behaviour changes. We also show that the mechanism is capable of producing an FRB signal at GHz frequencies in a relativistic magnetar wind close to the light cylinder and that this signal can escape the magnetar environment without significant damping.

We assess the impact of QSOs on the high redshift (z > 4) Intergalactic Medium using Monte Carlo realisations of QSO populations and the HeIII regions they generate, applied to the Sherwood-Relics simulations, allowing for uncertainties in the QSO luminosity function, its evolution, and QSO spectra and ages. While QSO luminosity functions based on optical-infrared selection are unable to reproduce the broadening HI Lyman-alpha optical depth distributions at z > 5, much broader distributions are found for the higher numbers of QSOs based on x-ray selection, suggesting a large QSO contribution to the ultra-violet background at z > 5 may offer an alternative to late reionization models to account for the broad HI Lyman-alpha optical depth distributions. Realisations using QSOs based on the higher QSO counts also much better recover the measured pixel flux auto-correlation function at z > 5. The HeIII regions from QSO sources according to both types of luminosity function suppress the pixel flux power spectrum on small scales, k > 0.02 s/km, while enhancing it on larger, both by amounts of up to tens of percent at z > 4, with the magnitude increasing with the intergalactic HeIII filling factor and the boost in temperature within the HeIII regions.

Human exploration of the moon is expected to resume in the next decade, following the last such activities in the Apollo programme time. One of the major objectives of returning to the Moon is to continue retrieving geological samples, with a focus on collecting high-quality specimens to maximize scientific return. Tools that assist astronauts in making informed decisions about sample collection activities can maximize the scientific value of future lunar missions. A lunar rock classifier is a tool that can potentially provide the necessary information for astronauts to analyze lunar rock samples, allowing them to augment in-situ value identification of samples. Towards demonstrating the value of such a tool, in this paper, we introduce a framework for classifying rock types in thin sections of lunar rocks. We leverage the vast collection of petrographic thin-section images from the Apollo missions, captured under plane-polarized light (PPL), cross-polarised light (XPL), and reflected light at varying magnifications. Advanced machine learning methods, including contrastive learning, are applied to analyze these images and extract meaningful features. The contrastive learning approach fine-tunes a pre-trained Inception-Resnet-v2 network with the SimCLR loss function. The fine-tuned Inception-Resnet-v2 network can then extract essential features effectively from the thin-section images of Apollo rocks. A simple binary classifier is trained using transfer learning from the fine-tuned Inception-ResNet-v2 to 98.44\% ($\pm$1.47) accuracy in separating breccias from basalts.

Context. Pulsar wind nebulae (PWNe) are a source of very high energy (VHE) radiation that can reach up to TeV and even PeV energies. Our work uses the pulsar tree, a graph theory tool recently presented to analyze the pulsar population to select the candidates of interest. Aims. We aim to discover possible detectable PWNe. We also aim to test to what extent the pulsar tree groups detectable PWNe despite it considering only pulsars' intrinsic properties. Methods. We select four pulsars as candidates for TeV PWNe based on their positions in the pulsar tree. Using observed and assumed ranges of values for relevant parameters, we anticipate possible spectral energy distributions (SEDs) of the PWNe of four pulsars (PSR J1208-6238, J1341-6220, J1838-0537, and J1844-0346) via a detailed time-dependent, leptonic model that was already found appropriate to describe almost all other detected nebulae. Results. We estimate the likelihood of detection for the four candidates studied by comparing the TeV fluxes predicted by the possible models with the sensitivities of different observatories. In doing so, we provide context for analyzing the advantages and caveats of the pulsar tree position as a marker for properties that go beyond the intrinsic features of pulsars that are considered in producing it.

We present a novel deep learning framework using Long Short-Term Memory (LSTM) networks to predict the spectra of galactic cosmic rays by leveraging historical solar activity data, addressing the limitations of traditional transport models. By incorporating multiple solar parameters, such as the heliospheric magnetic field, solar wind speed, and sunspot numbers, our model achieves accurate short-term and long-term predictions of cosmic-ray flux. The inclusion of historical cosmic-ray data enhances prediction accuracy, making the model highly effective for space weather forecasting. Moreover, it provides reliable one-day-ahead predictions of full spectra of cosmic rays for different species. Our approach surpasses traditional physics-based methods, providing a scalable, data-driven solution for reliable daily and long-term forecasts. This work paves the way for advanced models that can integrate broader observational data, with significant implications for space weather monitoring and mission planning.

We present an accelerated pipeline, based on high-performance computing techniques and normalizing flows, for joint Bayesian parameter estimation and model selection and demonstrate its efficiency in gravitational wave astrophysics. We integrate the Jim inference toolkit, a normalizing flow-enhanced Markov chain Monte Carlo (MCMC) sampler, with the learned harmonic mean estimator. Our Bayesian evidence estimates run on $1$ GPU are consistent with traditional nested sampling techniques run on $16$ CPU cores, while reducing the computation time by factors of $5\times$ and $15\times$ for $4$-dimensional and $11$-dimensional gravitational wave inference problems, respectively. Our code is available in well-tested and thoroughly documented open-source packages, ensuring accessibility and reproducibility for the wider research community.

The advent of large astronomical surveys has made available large and complex data sets. However, the process of discovery and interpretation of each potentially new astronomical source is, many times, still handcrafted. In this context, machine learning algorithms have emerged as a powerful tool to mine large data sets and lower the burden on the domain expert. Active learning strategies are specially good in this task. In this report, we used the PineForest algorithm to search for superluminous supernova (SLSN) candidates in the Zwicky Transient Facility. We showcase how the use of previously confirmed sources can provide important information to boost the convergence of the active learning algorithm. Starting from a data set of $\sim$14 million objects, and using 8 previously confirmed SLSN light curves as priors, we scrutinized 120 candidates and found 8 SLSN candidates, 2 of which have not been reported before (AT 2018moa and AT 2018mob). These results demonstrate how existing spectroscopic samples can be used to improve the efficiency of active learning strategies in searching for rare astronomical sources.

Demixing properties of planetary major constituents influence the interior structure and evolution of planets. Comparing experimental and computational data on the miscibility of hydrogen and water to adiabatic profiles suggests phase separation between these components occurs in the ice giants Uranus and Neptune. We aim to predict the atmospheric water abundance and transition pressure between the water-poor outer envelope and the water-rich deep interior in Uranus and Neptune. We construct seven H2-H2O phase diagrams from the available experimental and computational data. We compute interior adiabatic structure models and compare these to the phase diagrams to infer whether demixing is occurring. We obtain a strong water depletion in the top layer due to rain-out of water and find upper limits on the atmospheric water mass fraction Z_atm of 0.21 for Uranus and 0.16 for Neptune. The transition from the water-poor to the water-rich layer is sharp and occurs at pressures P_Z between 4 and 11 GPa. Using these constraints on Z_atm and P_Z, we find that the observed gravitational harmonics J2 and J4 can be reproduced if P_Z > 10 GPa in Uranus and > 5 GPa in Neptune, and if the deep interior has a high primordial water mass fraction of 0.8, unless rocks are also present. The agreement with J4 is improved if rocks are confined deeper than P_Z, for instance below a rock cloud level at 2000 K (20-30 GPa). These findings confirm classical few-layer models and suggest that a layered structure may result from a combination of primordial mass accretion and subsequent phase separation. Reduced observational uncertainty in J4 and its dynamic contribution, atmospheric water abundance measurements from an Orbiter with a Probe mission to Uranus (UOP) or Neptune, and better understanding of the mixing behaviour of constituents are needed to constrain the interiors of ice giants.

Non-circular (NC) motions represent the imprints of non-axisymmetric structures in galaxies, providing opportunities to study the physical properties of gas departing from circular rotation. In this work, we have conducted a systematic study of the non-circular motions in a sample of 1624 gas-rich disk galaxies from the MaNGA MPL-11. By using the H$\alpha$ velocity as a tracer of the disk rotation, we find indications that the amplitude of the non-circular motions is related to the stellar mass, with the low mass and late-type galaxies the most affected. In our sample, we find ratios of circular to non-circular rotation ranging from 5% to 20%. By implementing harmonic models to include NC motions associated with spiral arms and stellar bars, we find that the rotational curves traced with H$\alpha$ are barely affected by the NC induced by these structures. Consequently, in our sample, we do not find evidence that NC motions contribute to the scatter of the stellar Tully-Fisher relation. Our results suggest that non-circular motions might have a more localized effect in galaxies rather than a global one.

The kinematic properties of the Sco-Cen association have been studied using the spatial velocities of young stars. New kinematic age estimates for the three components of the association with the age of UCL and LCC being $17.7\pm2.4$ Myrs and the age of US being $6.4\pm1.7$ Myrs have been obtained. The parameters of the residual velocities US, UCL, and LCC ellipsoid have been estimated.

The origin of planetary mass objects (PMOs) wandering in young star clusters remains enigmatic, especially when they come in pairs. They could represent the lowest-mass object formed via molecular cloud collapse or high-mass planets ejected from their host stars. However, neither theory fully accounts for their abundance and multiplicity. Here, we show via hydrodynamic simulations that free-floating PMOs have a unique formation channel via the fragmentation of tidal bridge between encountering circumstellar disks. This process can be highly productive in density clusters like Trapezium forming metal-poor PMOs with disks. Free-floating multiple PMOs also naturally emerge when neighboring PMOs are caught by mutual gravity. PMOs may thus form a distinct population different from stars and planets.

We present the Variable Star PhasE Curve (VSPEC) Collection, a set of Python packages for simulating combined-light spectroscopic observations of 3-dimensional exoplanet atmospheres in the presence of stellar variability and inhomogeneity. VSPEC uses the Planetary Spectrum Generator's Global Emission Spectra (PSG/GlobES) application along with a custom-built multi-component time-variable stellar model based on a user-defined grid of stellar photosphere models to produce spectroscopic light curves of the planet-host system. VSPEC can be a useful tool for modeling observations of exoplanets in transiting geometries (primary transit, secondary eclipse) as well as orbital phase curve measurements, and is built in a modular and flexible configuration for easy adaptability to new stellar and planetary model inputs. We additionally present a set of codes developed alongside the core VSPEC modules, including the stellar surface model generator vspec-vsm, the stellar spectral grid interpolation code GridPolator, and a Python interface for PSG, libpypsg.

The origin of Galactic Cosmic Rays (GCRs) and the potential role of Supernova Remnants (SNRs) as cosmic-ray (CR) accelerators remain subjects of ongoing debate. To shed more light on this topic, we have studied the spectral shape of two SNRs, RX J1713.7-3946 and HAWC J2227+610, performing simulations for the Cherenkov Telescope Array Observatory (CTAO). The previous multi-wavelength (MWL) analysis on these two sources showed an important hadronic contribution at high energies. The interaction of the GCRs accelerated by the SNRs with the medium around the accelerator leads to a process of pion decay (PD) that produces gamma-rays ($\gamma$-rays). These emissions, detectable by CTAO, offer an indirect means of pinpointing the CR source. Two scenarios have been considered: the contribution of heavy CRs and different cut-off sharpnesses ($\beta$) of the particle spectra. The simulations were performed by using different CR composition distributions (protons, CNO, Fe) and different sharpness values ranging from $\beta$=0.5 to $\beta$=1.5. The results show that, in the cases studied here, CTAO will increase the sensitivity to the spectral shape of $\gamma$-rays. This allows us to distinguish protons from heavy CRs and obtain information on $\beta$ values and therefore on different acceleration scenarios.

Optical observations of Blue Compact Dwarf galaxies (BCDs) show they typically have high specific star formation rates and low metallicites. A subset of these galaxies (those with the lowest gas phase metallicities) display cometary optical morphologies similar to those found at high redshift. Whether this combination of properties predominantly arises from interactions with neighbours or accretion from the cosmic web or something else remains unclear. We used high resolution HI mapping to gain insights into the processes driving the observed properties of a sample of extremely metal poor (XMP) BCDs. We present Very Large Array B- and C-configuration HI mapping of four BCDs. For three of the targeted BCDs we also detect and map the HI in their nearby companions. In these three cases there is HI evidence for a recent flyby interaction between the BCD and a nearby companion. The HI evidence for recent interactions for these three BCDs is corroborated by our analysis of the tidal forces exerted on the BCDs by companions with available spectroscopic redshifts. For J0204-1009 we had sufficient spatial resolution to determine that it is dark matter dominated and estimate its DM halo mass to be in the range 1.2 x 10^11 to 5.2 x 10^11 solar masses. It is the most isolated BCD in our small sample, J0301-0052, which shows one of the most asymmetric HI morphologies. J0301-0052 has a similar cometary HI morphology to the BCD's optical morphology, although the HI column density maximum is projected at the end of the of the optical tail. Our HI observations suggest J0301-0052 may be undergoing a merger, while the other BCDs show evidence of a recent tidal interaction with a near neighbour. While our selection criteria favoured BCDs with companions our results are consistent with the earlier finding by other authors that most BCDs are associated with either mild tidal interactions or mergers.

We use a suite of the most recent cosmological observations to test models of dynamical dark energy motivated by quantum gravity. Specifically, we focus on hilltop quintessence scenarios, able to satisfy theoretical constraints from quantum gravity. We discuss their realisation based on axions, their supersymmetric partners, and Higgs-like string constructions. We also examine a specific parameterisation for dynamical dark energy suitable for hilltop quintessence. We then perform an analysis based on Markov Chain Monte-Carlo to assess their predictions against CMB, galaxy surveys, and supernova data. We show to what extent current data can distinguish amongst different hilltop set-ups, providing model parameter constraints that are complementary to and synergetic with theoretical bounds from quantum gravity conjectures, as well as model comparisons across the main dark energy candidates in the literature. However, all these constraints are sensitive to priors based on theoretical assumptions about viable regions of parameter space. Consequently, we discuss theoretical challenges in refining these priors, with the aim of maximizing the informative power of current and forthcoming cosmological datasets for testing dark energy scenarios in quantum gravity.

The thermal Sunyaev-Zel'dovich effect (tSZ) is a sensitive probe of cosmology, as it traces the abundance of galaxy clusters and groups in the late-time Universe. Upcoming cosmic microwave background experiments such as the Simons Observatory (SO) and CMB-S4 will provide low-noise and high-resolution component-separated tSZ maps covering a large sky fraction. The tSZ signal is highly non-Gaussian; therefore, higher-order statistics are needed to optimally extract information from these maps. In this work, we study the cosmological constraining power of several tSZ statistics -- Minkowski functionals (MFs), peaks, minima, and moments -- that have yielded promising results in capturing non-Gaussian information from other cosmological data. Using a large suite of halo-model-based tSZ simulations with varying $\Omega_{c}$ and $\sigma_{8}$ (154 cosmologies and over $800, 000$ maps, each $10.5\times10.5$ deg$^{2}$), we show that by combining these observables, we can achieve $\approx 29\times$ tighter constraints compared to using the tSZ power spectrum alone in an idealized noiseless case, with the MFs dominating the constraints. We show that much of the MF constraining power arises from halos below the detection threshold of cluster surveys, suggesting promising synergies with cluster-count analyses. Finally, we demonstrate that these statistics have the potential to deliver tight constraints even in the presence of noise. For example, using post-component-separation tSZ noise expected for SO, we obtain $\approx1.6\times$ and $\approx1.8\times$ tighter constraints than the power spectrum with MFs and all statistics combined, respectively. We show that the constraints from MFs approach the noiseless case for white-noise levels $\lesssim 1 \,\, \mu$K-arcmin.

Galactic-scale simulations rely on sub-grid models to provide prescriptions for the coupling between supernova (SN) feedback and the interstellar medium (ISM). Many of these models are computed in 1-D to allow for an efficient way to account for the variability of properties of their local environment. However, small-scale simulations revealed that the release of energy from SNe within molecular clouds can be highly asymmetrical. This is largely due to the presence of pre-SN feedback, such as ionizing radiation, that are able to carve cavities and channels around the progenitors prior to their detonation. Being partially confined, the SN energy escapes into the outer ISM preferentially through these channels, departing from the spherically symmetric 1-D descriptions. To understand by how much the feedback output could differ, we present a theoretical model for a semi-confined SN. The problem concerns a SN expanding into an evolved HII region, bounded by a molecular cloud with pre-existing vents. With the aid of simple 3-D hydrodynamical simulations, we show that this mode of energy release increases the dynamical impact of the outflows, and extends the timescales over which the SN is energetically coupled to the surrounding matter. We also show that the amount of small-scale solenoidal turbulence driven by semi-confined SNe may be amplified.