Papers for Tuesday, Apr 22 2025

Extensions to the $\Lambda\textrm{CDM}$ model prior to recombination can modify the growth of perturbations around radiation-matter equality, leaving a distinct signature in the matter power spectrum. Upcoming large-scale structure surveys will be sensitive to these features, allowing tests of early physics that are complementary to the CMB observations. In this paper, we forecast how well the combination of galaxy clustering, weak lensing and CMB lensing two point statistics, also known as $6\times2$pt analysis, will tighten constraints on extensions to the $\Lambda\textrm{CDM}$ model in the early Universe. We find significant improvements, in particular in the case of early dark energy, where the uncertainty on its density parameter could be divided by a factor of $3$ to $4$ when combining Euclid observables with Simons Observatory or CMB-S4, compared to using CMB observations alone. Testing for different scale cuts, we find that much of the constraining power comes from the largest scales which are less prone to systematic uncertainties. We take into account the most significant terms in the cross-covariance between large-scale structure tracers and CMB power spectra, which arises from gravitational lensing. Assessing the impact of this additional cross-covariance on the constraints, we find small corrections for most parameters, except for $A_s$ and $\tau$ where the lensing induced covariance leads to a more significant degradation of constraints. This forecast analysis highlights the potential of combining CMB and galaxy survey data to test the cosmological model. In particular, early Universe physics, relevant before recombination, stands out as a promising area that benefits substantially from this approach.

We present an efficient numerical algorithm for evolving self-gravitating systems of dark-matter particles that leverages the assumption of spherical symmetry to reduce the nominally six-dimensional phase space to three dimensions. It can be used to quickly determine numerically the evolution of an initially static stable self-consistent self-gravitating system if there is some additional or new physics. We illustrate here with four examples: (1) the effects of the growth of a supermassive black hole at the center; (2) the effects of stripping of the outer layers of the halo (a toy model for the effects of tidal stripping of galaxies); (3) the response of a self-gravitating system to dark matter that decays to a slightly less massive state; and (4) the effects of a slow change to Newton's constant. The approach can be extended to study dark matter with elastic and inelastic self-interactions and to study the process of virialization in spherical collapse. We describe some aspects of a code NSphere that implements this approach which we are making available.

Low-density ($\rho < 0.1 \rm{~g~cm^{-3}}$) hot Saturns are expected to quickly ($<100$ Myr) lose their atmospheres due to stellar irradiation, explaining their rarity. HAT-P-67 b seems to be an exception, with $\rho < 0.09 \rm{~g~cm^{-3}}$ and maintaining its atmosphere to well after 1 Gyr. We present a photometric and spectroscopic follow-up of HAT-P-67 b to determine how it avoided mass loss. HAT-P-67 b orbits a $V=10.1$ evolved F-type star in a 4.81 day orbit. We present new radial velocity observations of the system from the NEID spectrograph on the WIYN 3.5m Telescope from a follow-up campaign robust to stellar activity. We characterize the activity using photometry and activity indicators, revealing a stellar rotation period ($5.40\pm0.09$ d) near HAT-P-67 b's orbital period. We mitigate the stellar activity using a constrained quasi-periodic Gaussian process through a joint fit of archival ground-based photometry, TESS photometry, and our NEID observations, obtaining a planetary mass of $M_p = 0.45 \pm 0.15~M_{\rm J}$. Combined with a radius measurement of $R_p=2.140 \pm 0.025~R_{\rm J}$, this yields a density of $\rho_p = 0.061^{+0.020}_{-0.021} \rm{~g~cm^{-3}}$, making HAT-P-67 b the second lowest-density hot giant known to date. We find the recent evolution of the host star caused mass loss for HAT-P-67 b to only recently occur. The planet will be tidally disrupted in $\sim 100-250$ Myr, shortly after losing its atmosphere. With rapid atmospheric mass loss, a large, helium leading tail, and upcoming observations with the Hubble Space Telescope, HAT-P-67 b is an exceptional target for future studies, for which an updated mass measurement provides important context.

We analyze the dense gas kinematics in two Class 0/I protostellar cores, Per 30 and NGC 1333 IRAS 7, in the Perseus molecular cloud to determine whether their velocity structures are indicative of rotation. We examine the hyperfine structure of the N2H+ J=1-0 transition by combining 3" (900 AU) Atacama Large Millimeter/Submillimeter Array (ALMA) measurements with 9" (2700 AU) measurements from the Green Bank Telescope (GBT). We use the CASA Feather method to combine these data in order to maximize our sensitivity across spatial scales. We fit the N2H+ spectra to constrain the centroid velocity of the gas at each pixel and use these values to calculate the linear velocity gradient and specific angular momentum within apertures centered on each protostar with radii ranging from 5-60". Our results indicate that the velocity structure probed by the N2H+ emission is likely not a result of core rotation. These findings are consistent with other studies in the literature which indicate rotation is often not evident on scales less than 1000 AU. We instead suggest that the velocity structure we see is a result of torques caused by irregular density distributions in these protostellar systems.

Common radiative transfer methods, such as flux-limited diffusion (FLD) and the M1 closure, suffer from artificial interactions between crossing beams. In protoplanetary disks, this leads to an overestimation of the midplane temperature due to the merging of vertical inward and outward fluxes. Methods that avoid these artifacts typically require angular discretization, which can be computationally expensive. In the spirit of the two-stream approximation, we aim to remove the interaction between beams in a fixed spatial direction by introducing a half-moment (HM) closure, which integrates the radiative intensity over hemispheres. We derive a multidimensional HM closure via entropy maximization and replace it with an approximate expression that closely matches it, coinciding in the diffusion and free-streaming regimes while remaining expressible through simple operations. We implement the HM and M1 closures via implicit-explicit (IMEX) schemes, including multiple frequency groups. We test these methods in numerical benchmarks, including computing the temperature in an irradiated disk around a T Tauri star, comparing the results with Monte Carlo (MC) radiative transfer simulations. The resulting HM closure tends to the correct limit in the diffusion regime and prevents interactions between crossing fluxes in a chosen spatial direction. In disk simulations with 22 frequency groups, the M1 closure disagrees with the MC midplane temperature by up to 21%, while HM reduces this discrepancy to 6%. Even with just 3 frequency groups, HM significantly outperforms M1, with maximum departures of 8% compared to M1's 23%.

Determining how efficiently gas collapses into stars at high-redshift is key to understanding galaxy evolution in the Epoch of Reionization (EoR). Globally, this process is quantified by the gas depletion time ($t_{dep}$); on resolved scales, by the slope and normalization of the Kennicutt-Schmidt (KS) relation. This work explores the global ($\alpha_{[CII]}$) and spatially resolved ($W_{[CII]}$) [CII]-to-gas conversion factors at high-$z$ and their role in inferring reliable gas masses, surface densities, and $t_{dep}$ in the EoR. We select galaxies at 4<z<9 from the SERRA cosmological zoom-in simulation, that features on-the-fly radiative transfer and resolves interstellar medium properties down to $\approx$30 pc. The [CII] emission modelling from photodissociation regions allow us to derive global $\alpha_{ [CII]}$, and maps of $W_{[CII]}$. We study their dependence on gas metallicity (Z), density (n), Mach number (M), and burstiness parameter ($k_s$), and provide best fit relations. The $\alpha_{[CII]}$ decreases with increasing $Z$ and galaxy compactness, while the resolved $W_{[CII]}$ shows two regimes: at $Z< 0.2 Z_\odot$, it anticorrelates with n and Z, but not with $k_s$; above this threshold, it also depends on $k_s$, with more bursty regions showing lower conversion factors. This implies $W_{[CII]}\propto \Sigma_{[CII]}^{-0.5}$, as dense, metal-rich, and bursty regions exhibit higher [CII] surface brightness. Applying a constant $\alpha_{[CII]}$ overestimates $\Sigma_{gas}$ in bright $\Sigma_{[CII]}$ patches, thus flattening the KS slope and overestimating $t_{dep}$ by a factor of $\approx$4.

The Zwicky Transient Facility (ZTF), with its extensive optical monitoring capabilities, has provided an unprecedented opportunity to study the long-term variability of active galactic nuclei (AGNs). In this work, we present a comparative analysis of optical colour and brightness variability for $\mathrm{log(L_{Bol})}-z$ matched samples, consisting of 2095 Narrow-line Seyfert 1 galaxies similar to Broad-line Seyfert 1s (NLSy1_A) and 538 NLSy1_B galaxies matched with quasars (QSOs). The corresponding respective control samples of consist of 2380 Broad-line Seyfert 1 (BLSy1) galaxies and 741 QSOs. Using over six years of r-band and g-band light curves from the ZTF Data Release 22 (DR22), we characterize flux variability, fractional flux variability, and amplitude of temporal variability for each source in the samples. Our results indicate that BLSy1 galaxies exhibit significantly stronger variability compared to NLSy1_As, and similarly, QSOs show more variability than do NLSy1_Bs. To probe colour variability, we utilize quasi-simultaneous light curves, with half-hour epoch differences between $g$- and $r$-band measurements where colour index was evaluated using linear regression in magnitude-magnitude space. We find that large majorities of these sources -- 96% of NLSy1_A, 95% of BLSy1, 94% of NLSy1_B, and 91% of QSOs -- exhibit a clear "bluer-when-brighter" (BWB) trend. Furthermore, rest-frame structure function analysis reveals that BLSy1 galaxies are $1.42\pm0.06$ times more variable than NLSy1_A, while QSOs are $1.41\pm0.01$ times more variable than NLSy1_B. These results can provide valuable insights into the variability properties of AGN subclasses and their underlying physical drivers.

Hydrogen sulfide (H2S) is thought to be an important sulfur reservoir in interstellar ices. It serves as a key precursor to complex sulfur-bearing organics, and has been proposed to play a significant role in the origin of life. Although models and observations both suggest H2S to be present in ices in non-negligible amounts, its sublimation dynamics remain poorly constrained. In this work, we present a comprehensive experimental characterization of the sublimation behavior of H2S ice under astrophysically-relevant conditions. The sublimation behavior of H2S was monitored with a quadrupole mass spectrometer (QMS) during temperature-programmed desorption (TPD) experiments. These experiments are used to determine binding energies and entrapment efficiencies of H2S, which are then employed to estimate its snowline positions in a protoplanetary disk midplane. We derive mean binding energies of 3159\pm46 K for pure H2S ice and 3392\pm56 K for submonolayer H2S desorbing from a compact amorphous solid water (cASW) surface. These values correspond to sublimation temperatures of around 64 K and 69 K in the disk midplane, placing its sublimation fronts at radii just interior to the CO2 snowline. We also investigate the entrapment of H2S in water ice and find it to be highly efficient, with ~75-85% of H2S remaining trapped past its sublimation temperature for H2O:H2S mixing ratios of ~5-17:1. We discuss potential mechanisms behind this efficient entrapment. Our findings imply that, in protoplanetary disks, H2S will mostly be retained in the ice phase until water crystallizes, at radii near the water snowline, if it forms mixed into water ice. This has significant implications for the possibility of H2S being incorporated into icy planetesimals and its potential delivery to terrestrial planets, which we discuss in detail.

We investigate the presence of supermassive black hole (SMBH) binary signatures and the feasibility of identifying them through X-ray reflection spectra. The X-ray emitting region is modeled as a set of two mini-disks bound to the individual SMBHs separated by 100 $GM/c^2$ and the spectra calculated as a function of the mass, mass ratio, and total accretion rate of the binary. The X-ray reflection features are strongly influenced by the accretion-inversion phenomenon expected in SMBH binaries, which results in a wide range of ionization conditions in the two mini-disks. These are imprinted in the resulting composite spectra and the double-peaked and time-variable relativistic Fe K$\alpha$ line profiles. To test whether these features can be used as evidence for the presence of an SMBH binary, we fit mock 100\,ks observations with a single AGN model. For a $10^9$ $M_\odot$ binary targeted by Pulsar Timing Arrays (PTAs), at $z=0.1$ the single AGN model clearly fails to fit the data, while at $z=1$ the fit is acceptable but unable to converge on the SMBH spin. For a $10^6$ $M_\odot$ binary, a progenitor of a \textit{Laser Interferometer Space Antenna} (\textit{LISA}) source, spectral fitting is only possible at $z=0.1$, with the outcomes similar to the PTA binary at $z=1$. We also find that PTA binaries can be expected to show a distinct X-ray spectral variability in multi-epoch observations, whereas for \textit{LISA} precursors, orbital averaging results in the loss of spectral variability signatures.

We report the first X-ray and radio polarimetric results of the neutron star (NS) low-mass X-ray binary (LMXB) atoll-source 4U 1728-34 using the Imaging X-ray Polarimetry Explorer (IXPE) and Australia Telescope Compact Array (ATCA). We discovered that the X-ray source was polarized at PD = 1.9 +/- 1.0% (3-sigma errors) with a polarization angle of PA = -41 +/- 16 degree (3-sigma errors). Simultaneous Neutron Star Interior Composition Explorer (NICER) observations show that the source was in a relatively hard state, marking it as the first IXPE observation of an NS atoll source in the hard state. We do not detect any significant linear polarization (LP) in the radio band, with a 3-sigma upper limit of 2% at 5.5 GHz and 1.8% at 9 GHz. Combining the radio datasets provides the deepest upper limits on the radio polarization at < 1.5% on the linear and circular polarization (measured at 7.25 GHz). The X-ray polarimetric results suggest a source geometry with a Comptonization component possibly attributed to a boundary layer (BL) emission reflected off the disk, consistent with the other NS atoll sources.

Protoplanetary disks around luminous young A-type stars are prime observational laboratories to determine the abundances of complex organic molecules (COMs) present during planet formation. In contrast to their lower stellar mass counterparts, these warmer disks contain the sublimation fronts of complex molecules such as CH3OH on spatial scales accessible with the Atacama Large Millimeter/submillimeter Array (ALMA). We present ALMA observations of the Herbig Ae disk HD 100453 that uncover a rich reservoir of COMs sublimating from the dust cavity edge. In addition to CH3OH, we detect 13CH3OH for the first time in a Class II disk, revealing a factor of three enhancement of 13C in the disk large organics. A tentative detection of CH2DOH is also reported, resulting in a D/H of 1-2%, which is consistent with the expected deuterium enhancement from the low temperature CH3OH formation in molecular clouds and with the deuteration of CH3OH measured in comets. The detection of methyl-formate (CH3OCHO), at only a few percent level of CH3OH is an order of magnitude lower compared to claims towards other organic-rich Herbig Ae disks but is more in line with organic abundance patterns towards the earlier stages of star formation. Together these data provide multiple lines of evidence that disks, and therefore the planet and comet-forming materials, contain inherited interstellar ices and perhaps the strongest evidence to date that much of the interstellar organic ice composition survives the early stages of planet formation.

This paper presents the results of a detailed timing and spectral analysis of the TeV-detected blazar 1ES 1218+304, focused on the observations performed with the different instruments onboard the Neil Gehrels Swift Observatory in the period 2005-2024. The source showed various strengths of X-ray flaring activity and 0.3-10\,keV states differing by a factor up to 20 in brightness, exceeding a level of 2.7$\times$10$^{-10}$erg cm$^{-2}$s$^{-1}$ and representing the 3rd brightest blazar during the strongest flare. We detected tens of intraday variability instances, the majority of which occurred on sub-hour timescales and were consistent with the shock-in-jet scenario. The spectral properties were strongly and fastly variable, characterized by a frequent occurrence of very hard photon indices in the 0.3-10 keV and Fermi 0.3-300 GeV bands. The source exhibited very fast transitions of logparabolic-to-powerlaw spectra or conversely, possibly caused by changes of magnetic field properties over small spatial scales or by turbulence-driven relativistic magnetic reconnection. We detected various spectral features, which demonstrate the importance of the first-order Fermi mechanism operating by the magnetic field of changing confinement efficiencies and by the electron populations with different initial energy distributions, stochastic acceleration and cooling processes. In some periods, the source showed a softening at higher GeV-band energies, possibly due to the inverse-Compton upscatter of X-ray photons in the Klein-Nishina regime reflected in the positive correlation between X-ray and high-energy emissions.

We investigate the relationship between solar coronal holes and open-field regions using three-dimensional radiative magnetohydrodynamic (MHD) simulations combined with remote-sensing observations from the Solar Dynamics Observatory (SDO). Our numerical simulations reveal that magnetically open regions in the corona can exhibit brightness comparable to quiet regions, challenging the conventional view that open-field regions are inherently dark coronal holes. We find that the coronal brightness is primarily determined by the total energy input from photospheric magnetic activities, such as the small-scale dynamo, rather than differences in dissipative processes within the corona. Using synthesized EUV intensity maps, we show that brightness thresholds commonly used to identify coronal holes may overlook open-field regions, especially at lower spatial resolutions. Observational analysis utilizing SDO/HMI and AIA synoptic maps supports our simulation results, demonstrating that magnetic field extrapolation techniques, such as the Potential Field Source Surface (PFSS) model, are sensitive to the chosen parameters, including the source surface height. We suggest that discrepancies in estimates of open magnetic flux (the ``open flux problem'') arise both from the modeling assumptions in coronal magnetic field extrapolation and systematic biases in solar surface magnetic field observations. Our findings indicate the need for reconsidering criteria used to identify coronal holes as indicators of open-field regions to better characterize the solar open magnetic flux.

In the era of JWST, observations of hot Jupiter atmospheres are becoming increasingly precise. As a result, the signature of limb asymmetries due to temperature or abundance differences and the presence of aerosols can now be directly measured using transmission spectroscopy. Using a grid of general circulation models (GCMs) with varying irradiation temperature (1500 K - 4000 K) and prescriptions of cloud formation, we simulate 3D ingress/egress and morning/evening-limb transmission spectra. We aim to assess the impact that clouds, 3D temperature structure, and non-uniform distribution of gases have on the observed spectra, and how these inhomogeneities can be identified. A second goal is to assess the relative merits of two separate methods (ingress/egress v.s. morning/evening-limb spectroscopy) for isolating atmospheric asymmetries. From our models, it is evident that an east-west temperature difference is the leading order effect for producing ingress/egress or morning/evening-limb spectral differences. We additionally find that clouds contribute strongly to the observed limb asymmetry at moderate irradiation temperatures in our grid ($\sim 2000 \mathrm{K} < T_{\mathrm{irr}} < 3500 \mathrm{K}). At lower temperatures clouds equally dominate the optical depth on both limbs, while at higher temperatures the entire terminator region remains cloud-free. We develop limb asymmetry metrics that can be used to assess the degree of east-west asymmetry for a given planet and predict trends in these metrics with respect to irradiation temperature that are indicative of various physical processes. Our results are useful for predicting and diagnosing the signatures of limb asymmetries in JWST spectra.

We study the reheating phase following inflation in the context of single-field models, focusing on the perturbative decay of the inflaton into lighter particles. A general analytical framework is presented to compute the reheating temperature $T_{re}$ and related quantities by combining cosmological observations with model-dependent parameters. We derive expressions for $T_{re}$ for three types of interactions: gravitational, scalar, and Yukawa-type fermionic couplings, and apply these results to the class of $\alpha$-attractor inflationary models, which exhibit attractor behavior in the $(n_s, r)$ plane. The main goal of this work is to investigate how key cosmological quantities such as $T_{re}$, $N_{re}$, and $m_\phi$ among others, evolve with the scalar spectral index $n_s$ and the Yukawa coupling constant $y$, within a consistent analytical framework. Although the formulas used are approximate, they are sufficient to capture the qualitative behavior of the relevant quantities across a wide range of parameter values. Here, we are not interested in precise numerical approximations or data analysis, but rather in understanding the general trends and dependence of cosmological quantities of interest. In particular, tendencies observed in the figures, such as the sensitivity of $T_{re}$ to the coupling strength and the equation-of-state parameter $\omega_{re}$, reflect physical features that are not strongly affected by the approximations involved.

We present a catalog of Local Universe Near-Infrared Seyfert (LUNIS) \textit{K-}band integral field unit (IFU) data of 88 nearby Active Galactic Nuclei (AGN), curated from SINFONI/VLT and OSIRIS/Keck archival datasets. This catalog includes both type 1 and 2 Seyfert AGN probed down to scales of tens of parsecs with z $<$ 0.02 and spanning over five orders of magnitude in L$_{14-195keV}$ AGN luminosity. As part of this catalog we make publicly available for all galaxies the processed datacubes, a central 200 pc integrated spectrum, and two-dimensional maps of flux, velocity, and velocity dispersion for H$_2$ 1-0 S(1) 2.1218 $\mu$m, [Si VI] 1.9641 $\mu$m, and Br-$\gamma$ 2.1655 $\mu$m. The morphology and geometry of [Si VI], a tracer of AGN outflows, are reported for the 66$\%$ of galaxies with extended emission. We utilize this large sample to probe the behavior of molecular and ionized gas, identifying trends in the properties of the circumnuclear gas (surface brightness and velocity dispersion) with fundamental AGN properties (obscuration and X-ray luminosity). While there is significant variation in circumnuclear gas characteristics across the sample, we find molecular hydrogen to be less centrally concentrated and exhibit lower velocity dispersion relative to ionized gas. In addition, we find elevated molecular hydrogen surface brightness and decreased [Si VI] velocity dispersion in obscured relative to unobscured AGN. The [Si VI] and Br-$\gamma$ emission scale with L$_{14-195keV}$ X-ray luminosity, which, along with the elevated velocity dispersion compared to the molecular gas, verifies an association with AGN outflow processes.

When binaries are injected into low-angular-momentum orbits around a central supermassive black hole (SMBH), various outcomes can occur, including binary tidal breakup, double stellar disruptions and stellar collision. We use hydrodynamical simulations to study stellar collisions triggered by binary-SMBH encounters, examining both head-on and grazing collisions in deep ($\beta_b=5$) and gentle ($\beta_b=0.6$) encounters, where $\beta_b$ is the ratio of the binary tidal disruption radius to the binary pericenter distance to the SMBH. Head-on collisions consistently result in appreciable mass loss ($\sim 5\%$) and a single merger remnant. Grazing collisions in deep encounters typically leave two strongly disturbed stars with comparable mass loss, while in gentle encounters, multiple collisions eventually produce a single remnant with minimal mass loss ($\lesssim 1\%$). All merger remnants feature extended envelopes, making them susceptible to partial tidal disruptions when they return to the SMBH. The morphology and orbital energy distribution of collision-induced debris differ significantly from those of tidal disruption event (TDE) debris of single stars. Approximately half of the collision-generated debris falls back onto the SMBH, exhibiting a distinct time evolution of the fallback rate. We suggest that such mass loss and fallback can generate electromagnetic flares that mimic weak TDEs.

In this work, we present a study on the long time-scale period variations of four single-mode high-amplitude delta Scuti stars (HADS) via the classical $O-C$ analysis. The target HADS are (i) XX Cygni, (ii) YZ Bootis, (iii) GP Andromedae, and (iv) ZZ Microscopii. The newly determined times of maximum light came from the Transiting Exoplanet Survey Satellite (TESS), American Association of Variable Star Observers (AAVSO), and Bundesdeutsche Arbeitsgemeinschaft f$\ddot{\rm u}$r Ver$\ddot{\rm a}$nderliche Sterne (BAV) projects. Together with the times of maximum light obtained in the historical literature, the $O-C$ analysis was performed on these HADS, in which we obtained the linear period variation rates $\dot{P}/P$ as $(9.2 \pm 0.2) \times 10^{-9} \ \mathrm{yr^{-1}}$, $(3.2\pm 0.2)\times 10^{-9} \ \mathrm{yr^{-1}}$, $(4.22\pm 0.03) \times 10^{-8} \ \mathrm{yr^{-1}}$, and $(-2.06 \pm 0.02) \times 10^{-8} \ \mathrm{yr^{-1}}$, respectively. Based on these results and some earlier research, we also discuss the evolutionary stages and the mechanisms of the period variation of these four HADS.

In this study, we conducted a comprehensive analysis of the SX Phoenicis (SX Phe) type star CY Aquarii (CY Aqr). Our investigation included a detailed $O-C$ analysis based on a 90-year observational dataset, augmented by 1,367 newly determined times of maximum light. The $O-C$ diagram reveals that (i) the primary star of CY Aqr exhibits a linear period variation rate of $(1/P_0)(\mathrm{d} P/ \mathrm{d} t) = (2.132 \pm 0.002) \times 10^{-8} \mathrm{yr}^{-1}$ for its dominant pulsation mode; (ii) the primary star is disturbed by two companions and part of a triple system; (iii) Companion A has an orbital period of approximately 60.2 years and Companion B has an orbital period of approximately 50.8 years. It is highly probable that both Companion A and B are white dwarfs, with Companion A's elliptical orbit displaying an eccentricity of $e = 0.139 \pm 0.002$, which is the lowest confirmed value in similar binary and triple systems to date. Most notably, Companion A and B have masses that are identical within the uncertainties, with a mass ratio exceeding 0.99. Whether this is considered a coincidental event or the result of an underlying mechanism, CY Aqr is an exceptionally rare case that broadens our understanding of multiple star systems and offers a unique opportunity to delve into the enigmatic evolutionary histories of such configurations. Further intriguing characteristics of this system warrant investigation in future studies, based on additional observational data.

We theoretically analyze the dipole anisotropy observed in the quasar distribution from the CatWISE2020 catalog. The catalog data shows a peak around $z\approx 1$, suggesting the presence of a large-scale dipole component. We explore the possibility that this dipole could be driven by primordial density fluctuations from modes that were superhorizon at the time of CMB decoupling but have since entered the horizon and become subhorizon. In particular, we consider the impact of adiabatic modes with wavenumbers $k$ in the range $(10^{-4} - 4 \times 10^{-3})~\mathrm{Mpc}^{-1} $, corresponding to wavelength scales of several Gpc. Such modes can create large-scale density variations, likely causing anisotropies in the distribution of matter and, as a result, affecting the number density of observed quasars. We also demonstrate that a superhorizon curvature perturbations mode, with comoving wavenumber $k\lesssim0.3H_0$ can lead to a significant enhancement in the locally inferred Hubble constant. This effect offers a viable explanation for the observed discrepancy between local and CMB inferred measurements of $H_0$.

We present a comprehensive analysis of the 'canonical' soft state ($\gamma$, $\delta$, and $\phi$ spectral variability classes) of the black hole binary GRS 1915+105, using RXTE, AstroSat, and NuSTAR data from 1996 to 2017 to investigate the origin of High Frequency Quasi-periodic Oscillations (HFQPOs). Our findings reveal that HFQPOs occur only in the $\gamma$ and $\delta$ classes, with frequencies of $65.07-71.38$ Hz and are absent in the $\phi$ class. We observe an evolution of time-lag from hard-lag (1.59$-$7.55 ms) in RXTE to a soft-lag (0.49$-$1.68 ms) in AstroSat observations. Wide-band (0.7$-$50 keV) spectral modelling suggests that HFQPOs are likely observed with a higher covering fraction ($f_{cov} \gtrsim 0.5$), i.e., the fraction of seed photons being Comptonized in the corona, enhanced Comptonized flux ($\sim$ 38%), and lower optical depth ($\tau \lesssim 8.5$ ) in contrast to observations where HFQPOs are absent. We observed similar constraints for observing HFQPOs during an inter-class ($\phi \rightarrow \delta$) transition as well as in a few intra-class ($\delta \rightarrow \delta$) variations. We also find that the time-lag decreases as $\tau$ increases, indicating that a higher $\tau$ reduces Compton up-scattering, thereby decreasing the hard-lag. Interestingly, in RXTE observations, the hard-lag ($\sim$ 7 ms) gradually decreases as optical depth and Comptonization ratio increases, eventually becoming a soft-lag ($\sim$1 ms) in AstroSat observations. These constraints on spectro-temporal parameters for the likelihood of observing HFQPOs support a 'compact' coronal oscillation mechanism for generating HFQPOs, which we attempt to explain within the framework of a possible accretion scenario.

The phosphorus budget of planets is intertwined with their formation history and is thought to influence their habitability. The chemical reservoirs and volatile \emph{vs} refractory budget of phosphorus in planet-forming environments have so far eluded empirical characterisation. We employ high-resolution spectra from HST/STIS in the ultraviolet and APEX in the sub-mm to constrain the phosphorus budget in the well-characterized HD\,100546 star and protoplanetary disk system. We measure $\log{(P/H)_{\star}}=-7.50^{+0.23}_{-0.28}$ on the stellar surface, which traces the total inventory of P in accreting gas \emph{and }dust from the inner disk. The inner disk gas, inside of the main dust trap, has $\log{(P/H)_{\rm in}}\lesssim-8.70$, and the outer disk gas $\log{(P/H)_{\rm out}}\lesssim-9.30$. Phosphorus in the disk is carried by a relatively refractory reservoir, consistent with minerals such as apatite or schreibersite, or with ammonium phosphate salts, in terms of sublimation temperature. We discuss the impact this might have on the two protoplanets around HD\,100546. Our results contribute to our understanding of the chemical habitability of planetary systems and lay a foundation for future explorations, especially in the context of JWST and \emph{Ariel} which can study phosphorus in exoplanet atmospheres.

The Magellanic Clouds (MCs) are excellent locations to study stellar dust emission and its contribution to galaxy evolution. Through spectral and photometric classification, MCs can serve as a unique environment for studying stellar evolution and galaxies enriched by dusty stellar point sources. We applied machine learning classifiers to spectroscopically labeled data from the Surveying the Agents of Galaxy Evolution (SAGE) project, which involved 12 multiwavelength filters and 618 stellar objects at the MCs. We classified stars into five categories: young stellar objects (YSOs), carbon-rich asymptotic giant branch (CAGB) stars, oxygen-rich AGB (OAGB) stars, red supergiants (RSG), and post-AGB (PAGB) stars. Following this, we augmented the distribution of imbalanced classes using the Synthetic Minority Oversampling Technique (SMOTE). Therefore, the Probabilistic Random Forest (PRF) classifier achieved the highest overall accuracy, reaching ${89\%}$ based on the recall metric, in categorizing dusty stellar sources before and after data augmentation. In this study, SMOTE did not impact the classification accuracy for the CAGB, PAGB, and RSG categories but led to changes in the performance of the OAGB and YSO classes.

Circumstellar interaction of supernova (SN) ejecta is an essential process in its evolution and observations of SNe have found the signature of circumstellar interaction both in the early and late evolutionary phase of SNe. In this Letter, we show that if the SN forward shock plunges into tenuous stellar wind from dense circumstellar medium (CSM) in the vicinity of the progenitor (i.e., confined CSM), the subsequent time evolutions of the SN-CSM interaction system deviates from the prediction of self-similar solution. In this case, after all of the confined CSM is swept up by the SN forward shock (roughly $10$ days after the explosion), the propagation of the shocked shell will be driven by the freely expanding ram pressure of the confined CSM component, instead of the SN ejecta. Meanwhile, the forward shock decelerates faster than the prediction of thin-shell approximation once the confined CSM component reaches homologous expansion. This lasts until the reverse shock in the confined CSM component reaches the head of the SN ejecta, leading to the restoration of the system into the evolutionary model without confined CSM, where the SN ejecta drives the expansion of the system. We also show that this peculiar evolution will be reflected in observational signatures originating from SN-CSM interaction, taking rapid decline and rebrightening of radio emission as examples. Our results shed light on the importance of taking into account the effect of initial SN-CSM interaction even when we focus on observational properties of SNe a few years after the explosion.

Lorentz invariance violation in photons can be quantified by measuring the difference in arrival times between high- and low-energy photons originating from gamma-ray bursts (GRBs). When analyzing data, it is crucial to consider the inherent time delay in the emission of these photons at the source of the GRB. In a recent study, three distinct models were evaluated to explain the intrinsic emission times of high-energy photons by analyzing 14 multi-GeV photon events detected from 8 GRBs using the Fermi Gamma-ray Space Telescope (FGST). In this study, we examine three remarkable GRB photons recorded by different observatories: the 99.3~GeV photon from GRB 221009A observed by FGST, the 1.07~TeV photon from GRB 190114C detected by the Major Atmospheric Gamma Imaging Cherenkov (MAGIC) telescope, and the 12.2~TeV photon from GRB 221009A observed by the Large High Altitude Air-shower Observatory (LHAASO). Our analysis indicates that the newly proposed model with a linear relationship between photon energy and intrinsic emission time can offer a consistent framework to explain the behavior of all three exceptional photons with a Lorentz violation scale $E_{\rm LV}\sim 3\times 10^{17}$~GeV.

The young (23 Myr) nearby (19.4 pc) star $\beta$ Pictoris hosts an edge-on debris disk with two gas giant exoplanets in orbit around it. Many transient absorption features have been detected in the rotationally broadened stellar lines, which are thought to be the coma of infalling exocomets crossing the line of sight towards Earth. In the Solar System, the molecule cynaogen (CN) and its associated ionic species are one of the most detectable molecules in the coma and tails of comets. We perform a search for cyanogen in the spectra of $\beta$ Pictoris to detect or put an upper limit on this molecule's presence in a young, highly active planetary system. We divide twenty year's worth of HARPS spectra into those with strong exocomet absorption features, and those with only stellar lines. The high signal-to-noise stellar spectrum normalises out the stellar lines in the exocomet spectra, which are then shifted and stacked on the deepest exocomet absorption features to produce a high signal-to-noise exocomet spectrum, and search for the CN band head using a model temperature dependent cross-correlation template. We do not detect CN in our data, and place a temperature and broadening dependent 5$\sigma$ upper limit between 10$^{12}$ cm$^{-2}$ and 10$^{13}$ cm$^{-2}$, to be compared to the typical 10$^9$ - 10$^{10}$ cm$^{-2}$ expected from scaling of the values in the Solar System comets.

Understanding the origin of substructures in protoplanetary disks and their connection to planet formation is currently one of the main challenges in astrophysics. While some disks appear smooth, most exhibit diverse substructures such as gaps, rings, or inner cavities, with varying brightness and depth. As part of the Ophiuchus Disk Survey Employing ALMA (ODISEA), we previously proposed an evolutionary sequence to unify this diversity, driven by the formation of giant planets through core accretion and subsequent planet-disk interactions. By combining the disk evolution and planet formation code PLANETALP with the radiative transfer code RADMC-3D, we have now reproduced the key aspects of the proposed evolutionary sequence. Starting with a smooth disk (like e.g., WLY 2-63), we modeled the evolution of a fiducial disk with a 1 Jupiter-mass planet at 57 au. Within a few hundreds of orbits, a narrow gap forms, resembling ISO-Oph 17. By $\sim$0.1 Myr, the gap widens, and dust accumulates at the cavity edge, producing a structure similar to Elias 2-24. At $\sim$0.4 Myr, the disk evolves further into a morphology akin to DoAr 44, characterized by a smaller inner disk and a brighter inner rim. By $\sim$1 Myr, the system transitions to a single narrow ring, resembling RXJ1633.9-2442. This line of work strongly supports the planetary origin of substructures and enables the possibility of identifying a population of planets that is currently beyond the reach of more direct detection techniques.

Dusty stellar point sources are a significant stage in stellar evolution and contribute to the metal enrichment of galaxies. These objects can be classified using photometric and spectroscopic observations with color-magnitude diagrams (CMD) and infrared excesses in spectral energy distributions (SED). We employed supervised machine learning spectral classification to categorize dusty stellar sources, including young stellar objects (YSOs) and evolved stars (oxygen- and carbon-rich asymptotic giant branch stars, AGBs), red supergiants (RSGs), and post-AGB (PAGB) stars in the Large and Small Magellanic Clouds, based on spectroscopic labeled data from the Surveying the Agents of Galaxy Evolution (SAGE) project, which used 12 multiwavelength filters and 618 stellar objects. Despite missing values and uncertainties in the SAGE spectral datasets, we achieved accurate classifications. To address small and imbalanced spectral catalogs, we used the Synthetic Minority Oversampling Technique (SMOTE) to generate synthetic data points. Among models applied before and after data augmentation, the Probabilistic Random Forest (PRF), a tuned Random Forest (RF), achieved the highest total accuracy, reaching $\mathbf{89\%}$ based on recall in categorizing dusty stellar sources. Using SMOTE does not improve the best model's accuracy for the CAGB, PAGB, and RSG classes; it remains $\mathbf{100\%}$, $\mathbf{100\%}$, and $\mathbf{88\%}$, respectively, but shows variations for OAGB and YSO classes. We also collected photometric labeled data similar to the training dataset, classifying them using the top four PRF models with over $\mathbf{87\%}$ accuracy. Multiwavelength data from several studies were classified using a consensus model integrating four top models to present common labels as final predictions.

Multi-line spectropolarimetric observations allow for the simultaneous inference of the magnetic field at different layers of the solar atmosphere and provide insight into how these layers are magnetically coupled. The new upgrade of the Gregor Infrared Spectrograph (GRIS) instrument offers such a possibility, allowing for the simultaneous observation of the Ca II line at 8542 A, the Si I line at 10827 A, and the He I triplet at 10830 A in addition to some additional weaker spectral lines that can probe deeper in the photosphere. Because these spectral lines are sensitive to the plasma properties at different regions of the solar atmosphere, their combined analysis can help understand the stratification of its thermal and magnetic properties from the photosphere to the chromosphere. This work showcases recent observations of the upgraded GRIS at the active region AR13724, which shows the instrument's potential for unravelling the most minute details of solar phenomena. In particular, we analyse the spatial distribution of the polarisation signals as well as the distribution of Stokes profiles for different regimes of the magnetic field strength. We also conduct a preliminary data analysis using relatively simple and approximate methods.

We investigate potential deviations from cold dark matter (CDM) using the latest Baryon Acoustic Oscillations (BAO) measurements from the Dark Energy Spectroscopic Instrument (DESI). Analyzing DESI data alone constrains the dark matter equation of state parameter $w_{\mathrm{dm}} = -0.042^{+0.047}_{-0.024}$, revealing a mild preference for non-cold dark matter. This preference strengthens significantly in combined analyses, but reveals a striking tension in the inferred $w_{\mathrm{dm}}$ values from observations of different nature. The DESI+DESY5 combination yields $w_{\mathrm{dm}} = -0.084 \pm 0.035$, excluding CDM ($w_{\mathrm{dm}}=0$) at 2.4$\sigma$ significance. In contrast, Planck+DESI gives $w_{\mathrm{dm}} = 0.00077\pm0.00038$, differing from concordance model at 2$\sigma$ significance. The non-vanishing $w_{\mathrm{dm}}$ preference is particularly driven by low-redshift BAO measurements ($z<1.1$), while higher redshift data remain consistent with $\Lambda$CDM. The evidence for non-cold dark matter is more pronounced in DESI compared to the previous BAO surveys. All dataset combinations show significant improvement over the $\Lambda$CDM paradigm, providing compelling evidence for non-cold dark matter scenario.

Context:Halo formation time, which quantifies the mass assembly history of dark-matter halos, directly impacts galaxy properties and evolution. Although not directly observable, it can be inferred through proxies like star formation history or galaxy spatial distributions. Recent advances in machine learning enable more accurate predictions of halo formation time using galaxy and halo properties. Aims:This study aims to investigate a machine learning-based approach to predict halo formation time-defined as the epoch when a halo accretes half of its current mass-using both halo and baryonic properties derived from cosmological simulations. By incorporating properties associated with the brightest cluster galaxy located at the cluster center, its associated intracluster light component and satellite galaxies, we aim to surpass these analytical predictions, improve prediction accuracy and identify key properties that can provide the best proxy for the halo assembly history. Methods:Using The Three Hundred cosmological simulations, we train Random Forest (RF) and Convolutional Neural Network (CNN) models on halo and baryonic properties, such as mass, concentration, stellar and gas masses, and features of the brightest cluster galaxy and intracluster light. CNN models are trained on two-dimensional radial property maps. We also construct simple linear models using only observationally accessible features. Results:RF models show median biases of 4%-9% with standard deviations of 20%. CNN models reduce median bias to <4%, although they have higher scatter. Simple linear models using a limited number of observables achieve prediction accuracy comparable to RF models. Traditional relations between halo formation time and mass/concentration are preserved.

Upcoming wide-field time-domain surveys, such as the Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST) are expected to discover up to two orders of magnitude more strongly lensed supernovae per year than have so far been observed. Of these, Type IIn supernovae have been predicted to be detected more frequently than any other supernova type, despite their small relative detection fraction amongst non-lensed supernovae. However, previous studies that predict a large population of lensed Type IIn supernova detections model their time evolving spectrum as a pure blackbody. In reality, there is a deficit in the UV flux of supernovae relative to the blackbody continuum due to line-blanketing from iron-group elements in the ejecta and scattering effects. In this work we quantify the effect of this UV suppression on the detection rates by LSST of a simulated population of strongly lensed Type IIn supernovae, relative to a pure blackbody model, using a mock LSST observing run. With a blackbody model, we predict to detect $\sim$70 lensed Type IIn supernova per year with LSST. By modelling a similar UV deficit to that seen in superluminous supernovae, we recover 60 - 80% of the detections obtained using a pure blackbody model, of which $\sim$10 detections per year are sufficiently bright ($m_\textrm{i} < 22.5$ mag) and detected early enough (> 5 observations before lightcurve peak) to enable high-cadence spectroscopic follow up.

The ngVLA is a new interferometric radio astronomy facility with transformative capabilities, being developed by the National Radio Astronomy Observatory. It combines two orders of magnitude in frequency coverage, over 1.2 - 116 GHz, with unprecedented sensitivity, spatial resolution and spatial frequency coverage, opening up new discovery space, impacting nearly every area of astrophysics. The high sensitivity that enables the path breaking science goals, which in turn lead to stringent instrument requirements, also open up new approaches to meeting them, previously only possible in limited contexts. Chief among the requirements are the image dynamic range specifications of 45 dB and 35 dB at 8 GHz and 27 GHz in single pointing and mosaiced observations. As the baseline calibration strategy to meet these requirements, we leverage the high ngVLA sensitivity through routine use of self-calibration on short time scales to counter atmospheric delay fluctuations. We recognize the broader nature of the problem - requiring a certain dynamic range, DR, at a targeted science noise level {\sigma}_science presupposes the presence of bright emission in the field at a corresponding level of detected interferometric flux of \sim DR \times {\sigma}_science, by definition. With the problem posed broadly in this way, one can derive a general limit without recourse to background source counts and the attendant Poisson fluctuations of their occurrence, and which depends only on the number of antennas in the array and the science noise level of the observation, independent of the observing band, primary beam size and antenna SEFDs. This memo formally explores these ideas as the main strategy to achieve the ngVLA image dynamic range requirements.

The ngVLA is a new interferometric radio astronomy facility with transformative capabilities, being developed by the National Radio Astronomy Observatory. It combines two orders of magnitude in frequency coverage, over 1.2 - 116 GHz, with unprecedented sensitivity, spatial resolution and spatial frequency coverage, opening up new discovery space, impacting nearly every area of astrophysics. The high sensitivity that enables the path breaking science goals, which in turn lead to stringent instrument requirements, also open up new approaches to meeting them, previously only possible in limited contexts. Chief among the requirements are the image dynamic range specifications of 45 dB and 35 dB at 8 GHz and 27 GHz in single pointing and mosaiced observations. As the baseline calibration strategy to meet these requirements, we leverage the high ngVLA sensitivity through routine use of self-calibration on short time scales to counter atmospheric delay fluctuations and pointing self-calibration to correct for pointing errors. A key benefit of leveraging self-calibration techniques, where possible, is the a reduction in system complexity of a range of subsystems, which in turn improves system reliability. Self-calibration also promises the possibility of attaining thermal noise limited dynamic range performance in some cases. This presentation provides the bases for these approaches, illustrating them to make the case for their application to the ngVLA in parallel.

Measurement of environmental parameters is one of the basic requirements for the proper operation of a telescope. This memo is intended to provide guidance for the measurement accuracy requirements in the context of the ngVLA. It relies on previous work for ALMA (Mangum, 2001) and EVLA (Butler \& Perley, 2008) and a review of the subject by Mangum \& Wallace (2015). The local operational environment can be broadly divided into two categories: electromagnetic and physical. Meteorological parameters (weather) primarily constitute the physical environmental component and radio frequency interference (RFI) is the essential element of the electromagnetic environment. This memo focuses on the weather component and does not address the RFI, safety and physical infrastructure components. Under weather, the relevant topics are (1) the correction to pointing arising from refraction in the atmosphere (2) the different delays in the arrival times of signals at different antennas due to propagation in the atmosphere (3) monitoring weather parameters to provide operations support, e.g. in determining prevalence of precision or normal conditions, dynamic scheduling and the choice of antennas to constitute a sub-array with a given set of characteristics, among others, and (4) archival. Here we restrict ourselves to the first two topics which impact the data obtained and its calibration.

In this study, we further developed and investigated the dual parameter phenomenological dark energy model (H^2 + H^-2 dark energy model) derived from Kaniadakis holographic dark energy. On the theoretical basis of the original H^2 + H^-2 dark energy model (HHDE), four types of viscosities and seven types of interactions were introduced. These were combined pairwise, and a dynamical analysis was conducted on a total of 35 Modified H^2 + H^-2 Viscous Interacting Dark Energy (MHH-VIDE) models. The advantage of the HHDE model and MHH-VIDE models is that these models can greatly relieve the Hubble tension and cicumventing the potential issue of a 'big rip', and the dark energy is Quintom-like. In this article, we performed a three-dimensional dynamical analysis of the aforementioned models with interactions and viscosity, testing their viability. The results suggest that the nature of this dark energy is closer to a property of spacetime than a cosmological component. The phase diagram analysis reveals a modified radiation-dominated epoch, a transitional matter-dominated phase, and a late-time attractor corresponding to the dark-energy-driven acceleration phase.

The red supergiant (RSG) problem refers to the observed dearth of luminous RSGs identified as progenitors of Type II supernovae (SNe II) in pre-SN imaging. Understanding this phenomenon is essential for studying pre-SN mass loss and the explodability of core-collapse SNe. In this work, we re-assess the RSG problem using late-phase spectroscopy of a sample of 50 SNe II. The [O I] $\lambda\lambda$6300,6363 emission in the spectra is employed to infer the zero-age main sequence (ZAMS) mass distribution of the progenitors, which is then transformed into a luminosity distribution via an observation-calibrated mass-luminosity relation. The resulting luminosity distribution reveals an upper cutoff at log $L/L_{\odot} = 5.21^{+0.09}_{-0.07}$ dex, and the RSG problem is statistically significant at the 2$\sigma$ to 3$\sigma$ level. Assuming single RSG progenitors that follow the mass-luminosity relation of KEPLER models, this luminosity cutoff corresponds to an upper ZAMS mass limit of $20.63^{+2.42}_{-1.64}$ $M_{\odot}$. Comparisons with independent measurements, including pre-SN imaging and plateau-phase light curve modeling, consistently yield an upper ZAMS mass limit below about 25 $M_{\odot}$, with a significance level of 1-3$\sigma$. While each individual method provides only marginal significance, the consistency across multiple methodologies suggests that the lack of luminous RSG progenitors may reflect a genuine physical problem. Finally, we discuss several scenarios to account for this issue should it be confirmed as a true manifestation of stellar physics.

Searching for primordial gravitational wave in cosmic microwave background (CMB) polarization signal is one of the key topics in modern cosmology. Cutting-edge CMB telescopes requires thousands of pixels to maximize mapping speed. Using modular design, the telescope focal plane is simplified as several detector modules. Each module has hundreds of pixels including antenna arrays, detector arrays, and readout arrays. The antenna arrays, as the beam defining component, determine the overall optical response of the detector module. In this article, we present the developments of 6-inch broadband antenna arrays from 80GHz to 170GHz for the future IHEP focal plane module. The arrays are fabricated from 42 6-inch silicon wafers including 456 antennas, 7% more pixels than usual design. The overall in-band cross polarization is smaller than -20 dB and the in-band beam asymmetry is smaller than 10%, fulfilling the requirements for primordial gravitational wave search.

We present late-time optical and infrared (IR) observations of a sample of nine extragalactic Luminous Red Novae (LRNe) discovered in the last three decades. In all these cases, the LRN survivors fade below the pre-outburst luminosity of the progenitors in the optical region. Instead, they remain visible in the near-IR (NIR), and bright in the mid-IR (MIR) domains for years. We recover AT 1997bs in $Spitzer$ images from 2004, and a residual source is visible in HST and $JWST$ NIR images 27 years after the outburst. Its spectral energy distribution (SED) is consistent with that of a red supergiant star with a photospheric temperature of 3200 K and a radius of 220 $R_{\odot}$, without a significant circumstellar dust attenuation, as the source is not detected at 4.5 $\mu$m. Another LRN, AT 2011kp, is detected by $JWST$ 12.5 years after the outburst. From its SED, we find evidence for three black-body components: a warm/hot stellar component, a colder component from cool dust at 425 K, emitting between 3 and 6 $\mu$m, and a third even colder component responsible for the excess at 8 $\mu$m, possibly due to an IR echo. We construct the $[3.6]-[4.5]$ colour curves extending up to +7 years for six LRNe, which show a similar evolution: the MIR colour is $\sim0$ before the optical maximum light, it becomes bluer ($\sim-0.5$ mag) at around +1 year, then it gradually turns to redder colours ($\sim+0.5$ mag) in the following years, before reaching $[3.6]-[4.5]\sim+1.5$ mag 7 years after the outburst. Finally, we also estimate the masses and the temperatures of newly-formed dust years after the LRN onset for different chemical compositions (silicates or graphite) and grain sizes (0.1 or 1.0 $\mu$m). We find that LRNe produce dust masses of the order of (0.1-30)$\times10^{-6}$ $M_{\odot}$ between 3.6 and 13 years after the outbursts, with temperatures in the range from 350 to 770 K.

The $S_8$ parameter, which quantifies the amplitude of matter fluctuations on scales of $8h^{-1}$ Mpc, has been a source of tension between weak lensing surveys (e.g. KiDS, DES, HSC) and the Planck Cosmic Microwave Background (CMB) measurements. This discrepancy challenges the standard ${\Lambda}$CDM model and has become one of the most significant tensions in modern cosmology. The ${\Lambda'}$ model offers a potential resolution by introducing modifications to the cosmic growth history through alterations to the gravitational sector. The alterations involve including a Ricci soliton into Einstein's field equations which introduce a time dependent factor yielding a time varying cosmological constant ${\Lambda'}=(1-{\alpha(t)^2}(t)){\Lambda_{DE}}\frac{{\rho}_g}{\rho_{DE}}$ and subsequently the evolution of the cosmos. The Ricci soliton is sourced from gravitational energy density. In this study we analyze results from six surveys and compare the results for $w_a$ and $w_0$ with the ${\Lambda'}$ model. We also find ${\sigma}_8={0.750}_{-0.020}^{+0.020}$ , $S_8={0.788}$. These values are closer to some low $S_8$ measurements from weak lensing surveys (e.g DES, KiDS), which report $S_8 \approx 0.76-0.78$, suggesting that the model may alleviate the $S_8$ tension. High values of ${\alpha(t)}$ in the late universe are the cause of suppressed structure formation and low values of ${\Lambda'}$. The late universe in the ${\Lambda'}$ model is effectively or apparently 5-10% younger than in ${\Lambda}$CDM which translates to $H_0={72.734}_{-1.687}^{+1.687}$ km/s/Mpc, which is in agreement with late universe probes. ${\Lambda'}$ is classified under the dynamical dark energy models, however unlike alternatives, it does not invoke exotic particles nor phantom energy.

Semidetached binaries, distinguished by their mass transfer phase, play a crucial role in elucidating the physics of mass transfer within interacting binary systems. To identify these systems in eclipsing binary light curves provided by large-scale time-domain surveys, we have developed a methodology by training two distinct models that establish a mapping relationship between the parameters (orbital parameters and physical parameters) of semidetached binaries and their corresponding light curves. The first model corresponds to scenarios where the more massive star fills its Roche lobe, while the second model addresses situations where the less massive star does so. In consideration of the O'Connell effect observed in the light curves, we integrated a cool spot parameter into our models, thereby enhancing their applicability to fit light curves that exhibit this phenomenon. Our two-model framework was then harmonized with the Markov Chain Monte Carlo algorithm, enabling precise and efficient light-curve fitting and parameter estimation. Leveraging 2 minute cadence data from the initial 67 sectors of the Transiting Exoplanet Survey Satellite, we successfully identified 327 systems where the less massive component fills its Roche lobe, alongside three systems where the more massive component fills its Roche lobe. Additionally, we offer comprehensive fundamental parameters for these binary systems, including orbital inclination, relative radius, mass ratio, and effective temperature.

Cosmic voids are underdense regions within the large-scale structure of the Universe, spanning a wide range of physical scales - from a few megaparsecs (Mpc) to the largest observable structures. Their distinctive properties make them valuable cosmological probes and unique laboratories for galaxy formation studies. A key aspect to investigate in this context is the galaxy bias, $b$, within voids - that is, how galaxies in these underdense regions trace the underlying dark-matter density field. We want to measure the dependence of the large-scale galaxy bias on the distance to the void center, and to evaluate whether this bias profile varies with the void properties and identification procedure. We apply a void identification scheme based on spherical overdensities to galaxy data from the IllustrisTNG magnetohydrodynamical simulation. For the clustering measurement, we use an object-by-object estimate of large-scale galaxy bias, which offers significant advantages over the standard method based on ratios of correlation functions or power spectra. We find that the average large-scale bias of galaxies inside voids tends to increase with void-centric distance when normalized by the void radius. For the entire galaxy population within voids, the average bias rises with the density of the surrounding environment and, consequently, decreases with increasing void size. Due to this environmental dependence, the average galaxy bias inside S-type voids - embedded in large-scale overdense regions - is significantly higher ($\langle b\rangle_{\rm in} > 0$) at all distances compared to R-type voids, which are surrounded by underdense regions ($\langle b\rangle_{\rm in} < 0$). The bias profile for S-type voids is also slightly steeper. Since both types of voids host halo populations of similar mass, the measured difference in bias can be interpreted as a secondary bias effect.

The 21cm forest, narrow absorption features in the spectra of high redshift radio sources caused by intervening neutral hydrogen, offers a unique probe of the intergalactic medium and small-scale structures during reionization. While traditional power spectrum methods have been widely used for analyzing the 21cm forest, these techniques are limited in capturing the non-Gaussian nature of the signal. In this work, we introduce the Wavelet Scattering Transform (WST) as a novel diagnostic tool for the 21cm forest, which allows for the extraction of higher-order statistical features that power spectrum methods cannot easily capture. By decomposing simulated brightness temperature spectra into a hierarchy of scattering coefficients, the WST isolates both local intensity fluctuations (first-order coefficients) and scale-scale correlations (second-order coefficients), revealing the complex, multi-scale non-Gaussian interactions inherent in the 21cm forest. This approach enhances the power of 21cm forest in distinguishing between different cosmological models, such as Cold Dark Matter (CDM) and Warm Dark Matter (WDM), as well as scenarios with enhanced X-ray heating. Unlike traditional methods, which focus primarily on Gaussian statistics, the WST captures richer astrophysical and cosmological information. Our analysis shows that WST can significantly improve constraints on key parameters, such as the X-ray heating efficiency and the WDM particle mass, providing deeper insights into the early stages of cosmic structure formation.

We performed star counts in the region of the open cluster NGC 3532. The ranges of trigonometric parallaxes and proper motions containing all the stars of the cluster were determined using the stars with 5- and 6-parameter Gaia DR3 solutions. The estimated radius of the cluster was R_c=178+-3 arcminutes and the number of cluster stars was N_c=2200+-40. We estimate the number of stars with poor astrometric solutions that may be members of the cluster. For this purpose, we analyze the surface density distribution of stars with two-parameter Gaia DR3 solutions, stars with the parameter RUWE>1.4, and stars with large relative errors of trigonometric parallaxes in the vicinity of the cluster. We are looking for stars that fall within the area of the color-magnitude diagram occupied by probable members of the NGC 3532 cluster from the Hant&Reffert sample. The radial surface density profile plotted with such stars shows the concentration of stars toward the cluster center. An analysis of the profile yields an estimate of 2150+-230 stars that may be cluster members. Thus, nearly one half of cluster members can be lost when the probable members are selected only by exact astrometric data of Gaia DR3. Among these lost stars, there may be a significant number of unresolved binary and multiple systems.

We have discovered a substantial sodium doublet (Na D $\lambda\lambda$5890, 5896Å)-traced neutral outflow in a quenching galaxy JADES-GS-206183 at $z=1.317$ in GOODS-S field. Its JWST NIRSpec/MSA spectrum shows a significantly blueshifted and deep Na D absorption, revealing a neutral outflow with a velocity of $v_{\rm out}=828^{+79}_{-49}\,\mathrm{km\,s^{-1}}$ and a mass outflow rate of $\log(\dot{M}_{\rm out}/\mathrm{M_{\odot}\,yr^{-1}})=2.40^{+0.11}_{-0.16}$. The mass outflow rate of this outflow is higher than any of the neutral outflows identified previously beyond $z\sim1$ by the same line diagnostic and is comparable with those in local galaxies with extremely strong star formation activities or luminous AGN. Nonetheless, the best-fit SED modeling of JADES-GS-206183, based on its multi-band photometry from HST/ACS to JWST/NIRCam, suggests that the host galaxy now is quenched, and the Paschen $\alpha$ (Pa$\alpha$) emission in the FRESCO NIRCam grism spectrum confirms its current low star formation rate ($10.78\pm 0.55\,\mathrm{M_{\odot}\,yr^{-1}}$). More surprisingly, optical line ratio diagnostics indicate that the current AGN activity of JADES-GS-206183, if present, is weak. Even though we tentatively detect a broad component of the H$\alpha$ line, it is more likely tracing the ionized outflow than an AGN. The results demonstrate that the Na D outflow in JADES-GS-206183 is highly unlikely to be driven by current star formation or nuclear activity. Instead, we propose that the outflow that we are witnessing in JADES-GS-206183 may be a long-lasting fossil outflow, powered by previous AGN activity that has recently shut down.

Once per 10,000-100,000 years, an unlucky star may experience a close encounter with a supermassive black hole (SMBH), partially or fully tearing apart the star in an exceedingly brief, bright interaction called a tidal disruption event (TDE). Remnants of partial TDEs are expected to be plentiful in our Galactic Center, where at least six unexplained, diffuse, star-like "G objects" have already been detected which may have formed via interactions between stars and the SMBH. Using numerical simulations, this work aims to identify the characteristics of TDE remnants. We take 3D hydrodynamic FLASH models of partially disrupted stars and map them into the 1D stellar evolution code MESA to examine the properties of these remnants from tens to billions of years after the TDE. The remnants initially exhibit a brief, highly luminous phase, followed by an extended cooling period as they return to stable hydrogen burning. During the initial stage (< 10,000 yr) their luminosities increase by orders of magnitude, making them intriguing candidates to explain a fraction of the mysterious G objects. Notably, mild TDEs are the most common and result in the brightest remnants during this initial phase. However, most remnants exist in a long-lived stage where they are only modestly offset in temperature and luminosity compared to main-sequence stars of equivalent mass. Nonetheless, our results indicate remnants will sustain abnormal, metal-enriched envelopes that may be discernible through spectroscopic analysis. Identifying TDE survivors within the Milky Way could further illuminate some of the most gravitationally intense encounters in the Universe.

Lyman-$\alpha$ damping wings towards quasars provide a unique probe of reionization because their strength correlates strongly with the global volume-averaged neutral hydrogen (HI) fraction of the intergalactic medium (IGM). Cosmic variance in the IGM, however, is a major source of stochasticity since the local neutral environment around a quasar varies significantly even at fixed global neutral fraction. We show that the IGM damping wing carries additional information about this local ionization topology, unexploited by current analysis frameworks. We introduce a set of two new physically motivated summary statistics encoding the local information about the HI distribution in the IGM before it is altered by ionization radiation from the quasar, encompassing 1) the HI column density, weighted by a Lorentzian profile mimicking the frequency dependence of the Lyman-$\alpha$ cross section, and 2) the distance from the quasar to the first neutral patch. This description, when combined with the quasar's lifetime as a third parameter, reduces the IGM transmission scatter in the damping wing region of the spectrum to $\lesssim 1\,\%$ across the full range of physical parameter space. We introduce a simple procedure for generating synthetic HI sightlines around quasars and demonstrate that the resulting damping wing profiles are statistically indistinguishable from a realistic reionization topology. This opens the door for optimally extracting the salient local information encoded in the imprint in a model-independent fashion. In the context of a specific reionization model, measurements of these local parameters can be translated into constraints on the global timing of reionization, but in addition, they provide information about the reionization topology, hitherto unused. A marginally modified version of our framework can also be employed in the context of damping wings towards galaxies.

Supersymmetry-based hybrid inflation models (referred to as `spontaneously broken supersymmetry' by the Planck collaboration) are attractive for several reasons, including the appealing feature that inflation is associated with local gauge symmetry breaking in the early universe. Following the Planck collaboration's notation, the inflationary potential with sub-Planckian inflaton field values is given by: $V = \Lambda^4 [1 + \alpha_h \log(\phi / M_{Pl})] - m_{3/2} \Lambda^2 \phi + \Lambda^4 O((\phi / M_{Pl})^4)$. Here, $\Lambda = \sqrt{\kappa} M$ denotes the energy scale of inflation, $M$ is the gauge symmetry breaking scale, $\kappa$ is a dimensionless parameter that sets the inflaton mass ($\sqrt{2} \kappa M$), and $\alpha_h$ is determined from quantum corrections in terms of $\kappa$ and the underlying gauge group. A soft supersymmetry-breaking term proportional to the gravitino mass $m_{3/2}$ (~10 TeV) and linear in the inflaton field $\phi$ is also present during inflation. The final term in $V$ represents the leading-order supergravity correction. (Note that the last two terms were not taken into account in the Planck analysis.) We provide estimates for the parameters $\kappa$ (and $\alpha_h$) that yield a scalar spectral index $n_s$ in the range 0.96 to 0.98, which is fully consistent with recent P-ACT-LB measurements presented by the Atacama Cosmology Telescope, as well as earlier measurements by Planck. We recall that in the absence of the soft SUSY-breaking term proportional to $m_{3/2}$ in $V$, the spectral index $n_s = 1 - 1/N = 0.98$, where $N = 50$ denotes the number of e-foldings. The tensor-to-scalar ratio $r$ in this minimal model is tiny, but it can reach potentially observable values ($r \lesssim 0.01$) in non-minimal models.

Hypervelocity stars (HVSs) represent a unique class of objects capable of escaping the gravitational pull of the Milky Way due to extreme acceleration events, such as close encounters with the supermassive black hole at the Galactic center (GC), supernova explosions in binary systems, or multi-body dynamical interactions. Finding and studying HVSs are crucial to exploring these ejection mechanisms, characterizing central black holes, probing the GC environment, and revealing the distribution of dark matter in our galaxy. The Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) spectroscopic surveys have so far identified four B-type unbound HVSs. To expand this sample with the second-phase LAMOST survey that started in 2018, we conducted a systematic search for early-type HVSs using the LAMOST Data Release 10. We identified 125 early-type high-velocity candidates with total velocities exceeding 300 km\,s$^{-1}$. Among them, we report ten new unbound B- and A-type hypervelocity star (HVS) candidates (designated LAMOST-HVS5 through LAMOST-HVS14), tripling the number of unbound HVSs previously identified by LAMOST. Kinematic analyses suggest that these newly discovered HVS candidates likely originated either from the Galactic Center or via dynamical interactions. Future high-resolution follow-up observations promise to refine the stellar parameters, distances, and elemental abundances of these candidates, thereby providing deeper insights into their origins and broadening their potential applications across astrophysics.

The thorough understanding on the initiation of coronal mass ejections (CMEs), which is manifested as a slow rise of pre-eruptive structures before the impulsive ejection in kinematics, is the key for forecasting the solar eruptions. In our previous work, we showed that the slow rise of a hot flux rope with coronal mass density is caused by the moderate magnetic reconnection occurring in the hyperbolic flux tube (HFT) combined with the torus instability. However, it remains unclear how the initiation process varies when a filament is present in the pre-eruptive flux rope. In this work, we reveal the complete initiation route of a CME containing filament mass with a state-of-the-art full-magnetohydrodynamics simulation. The comprehensive analyses show that the filament mass has an important impact on the CME initiation through triggering and driving the slow rise of flux rope with its drainage, besides the contributions of HFT reconnection and torus instability. Finally, in combination with our previous work, we propose that the enhanced drainage of filament mass and various features related to the HFT reconnection, such as, the split of pre-eruptive structure and the pre-flare loops and X-ray emissions, can serve as the precursors of CME initiation in observations.

Despite its scientific importance, the low-surface-brightness universe has yet to be fully explored due to various systematic uncertainties that affect the achievable surface-brightness limit. Reducing these uncertainties requires very accurate data processing. The dark-sky flat is a widely used calibration frame for accurate flat-field correction, generated by combining the sky background from science images. However, the night sky will likely contain complex local fluctuations, thus may still lead to photometric errors in data calibrated with dark-sky flats. To address this concern, we conduct mock observations with semi-realistic sky simulation data and evaluate observation strategies to mitigate the impact of the fluctuating sky background. Our experiments consider two representative sky conditions (clear and dirty) and perform intensive comparative analysis on two observation methods (offset and rolling). Our findings suggest that the rolling dithering method, which incorporates the operation of camera rotation into conventional dithering, can provide more accurate dark-sky flats. Finally, we discuss the broader implications of this method through additional experiments examining several factors that may affect the imaging quality of observational data.

Recently, Sawala et al. 2025 claimed to refute the cosmological significance of the Giant Arc based on their analysis of the FLAMINGO-10K simulation data. In our paper here, we highlight several shortcomings of the authors' analysis. We then perform an enhanced analysis on the FLAMINGO-10K simulation data with applications of: the Single-Linkage Hierarchical Clustering (SLHC), the Convex Hull of Member Spheres (CHMS), and the Minimal Spanning Tree (MST) algorithms. Using the full $2.8^3$ Gpc$^3$ FLAMINGO-10K box, with subhaloes at $z=0.7$, and $100$ random realisations (from random subset selections) we find no gigaparsec structures in FLAMINGO-10K, and only a few ultra-large large-scale structures (uLSSs, structures exceeding a maximum pairwise separation of $370$ Mpc). Somewhat surprisingly, we found that the large-scale aspects of the FLAMINGO-10K data could be adequately represented by a Poisson point distribution. The enhanced analysis presented here further supports the remarkable nature of the Giant Arc as a cosmologically-significant structure. Of course, the Giant Arc is also accompanied by a second uLSS, the Big Ring. The analysis presented here builds on the work presented by Sawala et al., but amends the application of their statistical assessments. We do not yet know why there appears to be such a large discrepancy between the FLAMINGO-10K data and the observed LSS in MgII absorbers. Perhaps the results presented here might suggest that the GA, and especially the GA + BR, presents a more direct challenge to $\Lambda$CDM. In contrast to the conclusion of Sawala et al. that `gigaparsec patterns abound in a $\Lambda$CDM universe' we find that they are nowhere to be seen.

We study the G013.313+0.193 G013.313 region, a complex environment characterized by molecular cloud interactions indicative of cloud-cloud collision (CCC). Observations of the NH3(1,1) and (2,2) inversion transitions were obtained using the Nanshan 26 m radio telescope, while HCO+ (1-0), 12CO, 13CO, and C18O(1-0) transitions from the Purple Mountain Observatory Delingha 14 m telescope. Archival data are also included. We identified key observational signatures of CCC, including complementary spatial distributions, U-shaped structures, bridge features, and V-shaped velocity distributions. The position-velocity diagrams (P-V) reveal clear indications of gas interaction between two velocity components, suggesting an ongoing collision at an estimated angle of approximately 45 degree to the line of sight. The estimated collision timescale is 0.35-1.03 Myr, aligned with the inferred ages of young stellar objects (YSOs) in the region, supporting the hypothesis of collision-induced star formation. Hub-filament system (HFS) are identified in the compressed gas region, where filaments converge toward a dense hub, suggesting the CCC as a potential driver of HFS formation and massive star formation. The high column density suggests favorable conditions for the formation of massive stars. Although alternative kinematic drivers such as longitudinal collapse and shear motion are considered, CCC remains the most plausible explanation for the observed features. Our findings contribute to our understanding of the mechanisms of cloud dynamics and massive star formation in turbulent molecular environments.

Laboratory searching for dark matter is crucial for understanding several fundamental conundrums in physics and cosmology. Most cavity-based haloscope searches focus on the frequency range below 10 GHz, while the parameter space with higher frequency remains rarely explored, due to the challenges lying in the fabrication of microwave cavities. Here we report the first Q-band haloscope searching for dark photons with a 33.141 GHz cavity. A novel coupling tuning structure separated from the cavity was designed so as not to degrade the quality factor of the cavity. We have established the most stringent constraints $\chi<2.5\times10^{-12}$ at a confidence level of 90$\%$ in the frequency range from 33.139 GHz to 33.143 GHz, corresponding to the mass of dark photons ranging from 137.05 $\mu$eV to 137.07 $\mu$eV. The results surpass the previous astronomical constraints by nearly three orders of magnitude. This work has demonstrated the feasibility of dark matter haloscopes at Q band. In the future, the constraints can be further improved by more than one order of magnitude through low-temperature experiments, and the setup can be extended to search for axions, axion-like particles, and high-frequency gravitational waves.

We use the highest-resolution EAGLE simulation, Recal-L025N0752, to study the properties and formation of ultra-diffuse galaxies (UDGs). We identify 181 UDGs and find their properties closely match observations. The total masses of EAGLE UDGs range from ${\sim}5\times 10^{8}~M_{\odot}$ to ${\sim}2\times 10^{11}~M_{\odot}$, indicating that they are dwarf galaxies rather than failed $L_\star$ galaxies. EAGLE UDGs are not a distinct population, but rather a subset of dwarf galaxies, as their properties generally form a continuous distribution with those of normal dwarf galaxies. Unlike the situations in previous studies, the extended sizes of field UDGs in EAGLE are not driven by high halos spin or by supernova-induced stellar expansion, but instead largely arise from high spins in their star-forming gas and thus the newly formed stars at large radii. This might be attributed to galactic fountains, by which star-forming gas are launched to large halo-centric distances and acquire additional angular momentum through interactions with the circumgalactic this http URL satellite UDGs, ${\sim} 60 \%$ of them were already UDGs before falling into the host galaxy, while the remaining ${\sim} 40\%$ were normal galaxies prior to infall and subsequently transformed into UDGs due to tidal effects after infall.

Interstellar hydrogen atoms (H atoms) penetrate into the heliosphere through the region of the solar wind interaction with the interstellar plasma due to their large mean free path. Resonant charge exchange of H atoms with protons has been considered as the main interaction process between the components. In the majority of models, other processes like elastic H-H and H-p collisions are not included. Moreover, it has been assumed that the velocities of the colliding particles remain unchanged during charge exchange. This corresponds to the scattering on the angle of $\pi$ in the centre mass rest frame. The goal of this paper is to explore effects of the elastic H-H and H-p collisions as well as the angular scattering during charge exchange on the distribution of the interstellar atoms in the heliosphere and at its boundary. We present results of simple (and therefore, easily repeatable) kinetic model of the interstellar atom penetration through the region of the solar and interstellar winds interaction into the heliosphere. As a result of the model we compute the distribution function of the interstellar atoms at different heliospheric distances. Further, this distribution function is used to compute its moments and potentially observable features such as absorption and backscattered spectra in the Lyman-alpha line. Results show that there are differences in the behavior of the distribution function when considering elastic collisions and the changes in the moments of the distribution achieve 10%. Therefore, in cases where precise calculation of H atom parameters is essential, such as in the modeling of backscattered Lyman-$\alpha$ emission, elastic collisions must be considered.

The eruption of solar prominences can eject substantial mass and magnetic field into interplanetary space and cause geomagnetic storms. However, various questions about prominences and their eruption mechanism remain unclear. In particular, what causes the intriguing Doppler bullseye pattern in prominences has not yet been solved, despite some preliminary studies proposing that they are probably associated with counterstreaming mass flows. Previous studies are mainly based on single-angle and short timescale observations, making it difficult to determine the physical origin of Doppler bullseye patterns in prominences. Here, taking advantage of stereoscopic observations taken by the Solar Dynamics Observatory and the Solar Terrestrial Relations Observatory and a three-dimensional numerical simulation, we investigate the origin of prominence Doppler bullseye pattern by tracing a long-lived transequatorial filament/prominence from July 23 to August 4, 2012. We find that repeated coronal jets at one end of the prominence can launch the Doppler bullseye pattern. It is evidenced in our observations and simulation that during the forward traveling of jet plasma along the helical magnetic field structure of the prominence, part of the ejecting plasma can not pass through the apex of the prominence due to the insufficient kinetic energy and therefore forms a backward-moving mass flow along the same or neighboring magnetic field lines. This process finally forms counterstreaming mass flows in on-disk filaments. When the on-disk filament rotates to the solar limb to be a prominence, the counterstreaming mass flows are naturally observed as a Doppler bullseye pattern.

The prompt emission of gamma-ray bursts (GRBs) is supposed to be released from the relativistic jet launched from the central engine. Apart from the non-thermal nature of the spectra in a majority of GRBs, there is evidence for the presence of quasi-thermal components in the prompt emission of a few GRBs according to observations by Fermi satellite. On the other hand, the GRB jet has been revealed as structured in recent research. The theoretical observed spectra of photosphere emissions by an off-axis observer and the dependence of the spectra on the viewing angle under the assumption of structured jets remain unexplored. In this paper, we numerically calculate the instantaneous photosphere spectra by different viewing angles from a structured jet, from which relevant temporal and spectral characteristics are derived. Moreover, we address the necessity of proper treatment of the outflow boundary in the photosphere emission scenario. Furthermore, our calculations suggest that the Einstein Probe and Space-based multi-band astronomical Variable Object Monitor will have the capability to detect the short GRBs similar to GRB 170817A up to a luminosity distance of 200Mpc if the off-axis viewing angle is less than 10 degrees.

We report on the spectral energy distributions (SEDs) of infrared-bright dust-obscured galaxies (DOGs) with $(i - [22])_{\rm AB} \geq 7.0$. Using photometry from the deep and wide Ultraviolet Near-Infrared Optical Northern Survey, combined with near-IR and mid-IR data from the UKIRT Infrared Deep Sky Survey and the Wide-field Infrared Survey Explorer, we successfully identified 382 DOGs in $\sim$ 170 deg$^2$. Among them, the vast majority (376 DOGs) were classified into two subclasses: bump DOGs (132/376) and power-law (PL) DOGs (244/376), which are dominated by star formation and active galactic nucleus (AGN), respectively. Through the SED analysis, we found that roughly half (120/244) of the PL DOGs show ``broken'' power-law SEDs. The significant red slope from optical to near-IR in the SEDs of these ``broken power-law DOGs'' (BPL DOGs) probably reflects their large amount of dust extinction. In other words, BPL DOGs are more heavily obscured AGNs, compared to PL DOGs with non-broken power-law SEDs.

Slow magnetoacoustic waves (SMAWs) have been considered in the past as a possible candidate for chromospheric heating. This study analyzes 20 active regions observed between 2012 and 2016 to examine the amplitude and energy flux variation of SMAWs in the umbral atmosphere. Six different wavelength channels from the Atmospheric Imaging Assembly onboard the Solar Dynamics Observatory, covering regions from the photosphere to the low corona, were utilized for this purpose. The wave amplitude estimations show a gradual increase in 3-minute oscillation amplitude, peaking between 700--900 km, followed by a steady decrease; at altitudes greater than 1800 km, the amplitude appears to increase and then decrease again. The corresponding energy flux, on the other hand, displays a steady and monotonous decrease with a significant reduction in value from approximately $3.32 \pm 0.50~\mathrm{kW,m^{-2}}$ near the photosphere to about $(6.47 \pm 3.16) \times 10^{-4}~\mathrm{W,m^{-2}}$ at an altitude of 2585 km. This decay may be attributed to radiative damping and shock dissipation in the lower altitudes, and thermal conduction and viscosity in the higher altitudes. The missing flux is a factor of 3--15 lower than that required to counterbalance the chromospheric radiative losses.

$\delta$ Scuti ($\delta$ Sct) stars are potential distance tracers for studying the Milky Way structure. We conduct a comprehensive analysis of the period-luminosity (PL) and period-luminosity-metallicity (PLZ) relation for $\delta$ Sct stars, integrating data from the Zwicky Transient Facility (ZTF), the Transiting Exoplanet Survey Satellite (TESS), Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST), Apache Point Observatory Galactic Evolution Experiment (APOGEE), and Gaia. To mitigate the impact of the Gaia parallax zero point offset, we applied a correction method, determining the optimal zero point value to be $zp_\varpi = 35 \pm 2 \, \mu\text{as}$. Using the three best bands, by varying the parallax error threshold, we found that the total error of the PLR zero point was minimized to 0.9\% at a parallax error threshold of 6\%. With this threshold, we derived the PL and PLZ relation for nine bands (from optical to mid-infrared) and five Wesenheit bands. Through our analysis, we conclude that the influence of metallicity on the PLR of $\delta$ Sct stars is not significant, and the differences across various bands are minimal.

Potential vibrational modes associated with diffuse interstellar bands (DIBs) could be discerned by examining energy differences between correlated DIBs. Consequently, $\approx 10^3$ higher correlated DIB pairs ($r-\sigma_r \ge 0.8$, $\ge 12$ sightlines) were extracted from the Apache Point Observatory DIB catalog, and their energy spacings computed. In this first macro exploratory step, a histogram possibly reveals chemical bond signatures of C$\equiv$C, C$\equiv$N, S$-$H, C$-$O, C$=$O, Si$-$H, N$-$H, C$-$H (aliphatic), C$\mathbf{^{\underline{...}}}$C (in-ring), and aromatics (C$-$H stretch, C$\mathbf{^{\underline{...}}}$C in-ring, oop C$-$H bending, and overtones). Continued research is required to (in)validate the histogram approach, mitigate noise, scrutinize maxima, break degeneracies, and converge upon an optimal framework.

Filaments are special plasma phenomena embedded in the solar atmosphere, characterized by unique thermodynamic properties and magnetic structures. Magnetohydrodynamic (MHD) simulations are useful to investigate the eruption mechanisms of filaments. We conduct a data-constrained zero-$\beta$ MHD simulation in spherical coordinates to investigate a C3.5 class flare triggered by an eruptive filament on 2022 August 15 in a decaying weak active region 13079. We reconstruct the three-dimensional coronal magnetic field using vector magnetograms and synoptic maps from the Solar Dynamics Observatory/Helioseismic and Magnetic Imager (SDO/HMI). We transform vector magnetic field into Stonyhurst heliographic spherical coordinates combined with a synoptic map and constructed a potential field source surface (PFSS) model with a magnetic flux rope (MFR) embedded using the Regularized Biot--Savart Laws (RBSL). Subsequently, we conduct a spherical zero-$\beta$ MHD simulation using the Message Passing Interface Adaptive Mesh Refinement Versatile Advection Code (MPI-AMRVAC) and replicated the entire dynamic process of the filament eruption consistent with observations. With the calculation of time-distance profile, Qusai-Separatrix Layers (QSL), and synthetic radiation from simulated current density, we find a good agreement between our simulation and observations in terms of dynamics and magnetic topology. Technically, we provide a useful method of advanced data-constrained simulation of weak active regions in spherical coordinates. Scientifically, the model allows us to quantitatively describe and diagnose the entire process of filament eruption.

The 21-cm forest offers a powerful cosmological probe of the thermal history and small-scale structure of the intergalactic medium during the Epoch of Reionization (EoR). Its success, however, critically depends on the availability of high-redshift radio-loud quasars (HzRLQs) as background sources. In this work, we investigate the configuration requirements for a Moon-based low-frequency radio interferometer aimed at maximizing the detection of HzRLQs for future 21-cm forest studies. Building upon a previously developed quasar luminosity function (QLF), we forecast HzRLQ abundances under various array configurations. Assuming a total survey area of $10^4\,\mathrm{deg}^2$ and 1 year of observation, we compare continuum surveys with 10 MHz bandwidth and 21-cm forest surveys with 5 kHz resolution. Our results show that a minimum collecting area of $\sim$6 500 m$^2$ enables detection at $z \sim 6$, while SKA-like arrays ($N_{\mathrm{st}} = 512$) extend the detection limit to $z \sim 10$ for 21-cm forest survey and $z \sim 16$ for continuum survey. Larger arrays with $N_{\mathrm{st}} = 2048$ can reach $z \sim 11$ in 21-cm forest mode. We also explore configurations that maintain fixed collecting areas while increasing the number to enhance survey efficiency. This boosts source detection but significantly increases the data volume and computational demands. These results underscore the importance of optimizing array design for different survey goals and balancing sensitivity, spectral resolution, and data management. A well-designed Moon-based array could open a new observational window on reionization and early cosmic structure formation.

We present multiwavelength analysis of the pre-flare phase and onset of the powerful X4.9 near-limb eruptive solar flare on February 25, 2014, revealing the tether-cutting (TC) geometry. We aim at determining relationship between the region of pre-flare energy release with the regions where the flare started to develop, and to investigate a detailed chronology of energy release during the pre-flare time interval and the beginning of the impulse phase. Using X-ray, ultraviolet and radio microwave data we found that the pre-flare energy release site was compact and localized in the vicinity of TC interaction of magnetic structures near the polarity inversion line. The analysis indicates that a pre-flare current sheet (CS) could be in this region. Good correspondence between the location of the pre-flare and flare emission sources visible at the very beginning of the impulsive phase is shown. We found relationship between dynamics of the energy release in the pre-flare CS and formation of the future flare eruptive structure. The growth of the magnetic flux rope was associated with activation of plasma emissions, flows and an increase of UV radiation fluxes from the region where the pre-flare CS was located. The eruptive flux rope gradually grew due to feeding by magnetized plasma ejected from the reconnecting pre-flare CS. Finally, it is shown that the most probable trigger of the eruption was a local fast microflare-like magnetic reconnection in the pre-flare CS. Some local instability in the pre-flare sheet could lead to a transition from the slow to fast reconnection regime. As a result, an ejection from the sheet was initiated and the eruptive flux rope lost its stability. Then, the eruptive flux rope itself initiated formation of the main reconnecting flare CS as in the Standard Flare Model during its movement, and intense emissions associated with the impulsive phase were observed.

A search was carried out for pulsars with periods (P) from 2 to 90 s in daily observations carried out over an interval of 5 years in a area measuring 6300 this http URL. The data was obtained on a Large Phased Array (LPA) at a frequency of 111 MHz. The periodograms calculated using the Fast Folding Algorithm (FFA) were used for the search. To increase the sensitivity, the periodograms obtained in different observation sessions were added together. Of the 14 known pulsars that entered the study area, with periods of P>2 s and dispersion measures (DM) less than 200 pc/cm$^3$, 9 were detected. 2 new pulsars have been found. The mean profiles of pulsars are obtained and estimates of their flux densities are given. The open pulsar J1951+28, with a period of P = 7.3342 s and DM = 3.5 pc/cm$^3$, turned out to be one of the pulsars closest to the Sun. The absence of new pulsars with periods of tens of seconds with a search sensitivity of 1 mJy outside the Galactic plane indicates a low probability of the existence of pulsars with extremely long periods. Most likely, the recently found sources of periodic radiation with periods from a minute to tens of minutes are white dwarfs.

Fast Radio Bursts (FRBs) have emerged as a powerful tool for cosmological studies, particularly through the dispersion measure-redshift ($\mathrm{DM}-z$) relation. This work proposes a novel calibration method for FRBs using the Yang-Li-Zhang (YLZ) empirical relation, which links the rotation measure (RM) of FRBs to the luminosity of their associated persistent radio sources (PRS). We demonstrate that this approach provides independent constraints on cosmological parameters, bypassing limitations inherent to traditional $\mathrm{DM}-z$ method. Utilizing the current sample of four YLZ-calibrated FRBs, we derive a Hubble constant measurement of $H_0 = 86.18_{-14.99}^{+18.03}\ \mathrm{km\ s^{-1}\ Mpc^{-1}}$ (68\% CL). Monte Carlo simulations indicate that a future catalog of 400 FRB-PSR systems could reduce the relative uncertainty of $H_0$ to 4.5\%. Combining YLZ-calibrated FRBs with $\mathrm{DM}-z$ sample reveals critical synergies: joint analysis of equalized samples ($N=100$ for both methods) reduces the relative uncertainty of $H_0$ to 2.9\%, mainly because the incorporation of PRS observations substantially mitigates the degeneracy between the parameters such as IGM baryon mass fraction ($f_{\rm IGM}$) and other cosmological parameters inherent to the $\mathrm{DM}-z$ relation.

This work investigates a potential time dependence of the absolute magnitude of Type Ia Supernovae (SN Ia). Employing the Gaussian Process approach, we obtain the SN Ia absolute magnitude and its derivative as a function of redshift. The data set considered in the analysis comprises measurements of apparent magnitude from SN Ia, Hubble rate from cosmic chronometers, and the ratio between angular and radial distances from Large-Scale Structure data (BAO and voids). Our findings reveal good compatibility between the reconstructed SN Ia absolute magnitudes and a constant value. However, the mean value obtained from the Gaussian Process reconstruction is $M=-19.456\pm 0.059$, which is $3.2\sigma$ apart from local measurements by Pantheon+SH0ES. This incompatibility may be directly associated to the $\Lambda$CDM model and local data, as it does not appear in either model-dependent or model-independent estimates of the absolute magnitude based on early universe data. Furthermore, we assess the implications of a variable $M$ within the context of modified gravity theories. Considering the local estimate of the absolute magnitude, we find $\sim3\sigma$ tension supporting departures from General Relativity in analyzing scenarios involving modified gravity theories with variations in Planck mass through Newton's constant.

We present a Simulation-Based Inference (SBI) framework for cosmological parameter estimation via void lensing analysis. Despite the absence of an analytical model of void lensing, SBI can effectively learn posterior distributions through forward modeling of mock data. We develop a forward modeling pipeline that accounts for both cosmology and the galaxy-halo connection. By training a neural density estimator on simulated data, we infer the posteriors of two cosmological parameters, $\Omega_m$ and $S_8$. Validation tests are conducted on posteriors derived from different cosmological parameters and a fiducial sample. The results demonstrate that SBI provides unbiased estimates of mean values and accurate uncertainties. These findings highlight the potential to apply void lensing analysis to observational data even without an analytical void lensing model.

Recent data from DESI, in combination with other data, provide moderate evidence of dynamical dark energy, $w\neq-1$. In the $w_0, w_a$ parametrization of $w$, there is a preference for a phantom crossing, $w<-1$, at redshift $z\sim0.5$. In general relativity, the phantom equation of state is unphysical. Thus it is important to check whether phantom crossing is present in other physically self-consistent models of dark energy that have equivalent evidence to the $w_0, w_a$ parametrization. We find that thawing quintessence with nonzero cosmic curvature can fit the recent data slightly better than $w_0, w_a$ in a flat background. The phantom crossing may be a spurious artifact of a parametrization that is not based on a physical model.

The latest ACT data release disfavors the attractor $n_s=1-2/N$. In inflationary models with nonminimal coupling, such attractors typically arise in the strong coupling limit. To align with observational constraints, we focus on nonminimal coupling models with small coupling constants. For the model with the coupling function $\Omega(\phi) = 1 + \xi f(\phi)$ and the potential $V(\phi) = \lambda^2 f^2(\phi)$, we find that observational data constrain the parameters as $0.1 \lesssim \xi \lesssim 35$ and $0 \lesssim k \lesssim 1.5$ for $f(\phi) = \phi^k$ at the $1\sigma$ confidence level. With the help of the nonmiminal coupling $\Omega(\phi) = 1 + \xi \phi^2$, the hilltop inflation and power-law inflation models with power indices $2/3$ and $1/3$ can be consistent with observational data within the $1\sigma$ range. We also give the viable parameter regions for $\xi$ for these three models.

A fundamental question in cosmology is whether dark energy evolves over time, a topic that has gained prominence since the discovery of cosmic acceleration. Recently, the DESI collaboration has reported increasing evidence for evolving dark energy using combinations of cosmic microwave background (CMB), type Ia supernova (SN), and their new measurements of baryon acoustic oscillations (BAO). However, our analysis reveals that these combinations are problematic due to clear tensions among the CMB, BAO and SN datasets. Consequently, DESI's claim of dynamical dark energy (DDE) is not robust. A more reliable approach involves constraining the evolution of dark energy using each dataset independently. Through a statistical comparison for each dataset, on average, we find that DDE is strongly preferred over the $\Lambda$CDM model. This suggests that DDE likely exists, although its real parameter space remains elusive due to weak constraints on the dark energy equation of state and inconsistencies among the datasets. Interestingly, when considering DDE, none of the individual datasets -- including CMB, DESI DR2, Pantheon+, Union3, and DESY5 -- can independently detect cosmic acceleration at a significant level. Our findings not only clarify the current understanding of the nature of dark energy but also challenge the established discovery of cosmic acceleration and the long-held notion that dark energy exerts negative pressure. Both individual and combined datasets suggest that the ultimate fate of the universe is likely to be dominated by matter rather than dark energy.

We characterize the warm circumgalactic medium (CGM) of a dwarf galaxy pair with properties similar to the Magellanic Clouds in the \textsc{Hestia} cosmological simulations. The system consists of a massive dwarf ($M_{\rm halo} \sim 10^{11.5} M_{\odot}$) and a lower-mass companion ($M_{\rm halo} \sim 10^{10} M_{\odot}$), dynamically evolving in isolation before infall into a Milky Way-mass halo. The massive dwarf hosts a warm coronal gas envelope with a temperature of $T \sim 3 \times 10^5$ K, consistent with expectations for virialized CGM in dwarf halos. Tidal interactions produce a neutral gas stream that extends over $\sim 150$ kpc, with an \ion{H}{1} mass of $M_{\rm HI} \sim 10^8 M_{\odot}$, similar to the Magellanic Stream. Furthermore, in the \textsc{Hestia} simulation suite, we find that coronal gas is ubiquitous in all halos with $M_{\rm halo} > 10^{11} M_{\odot}$, implying that massive dwarfs generically develop extended gaseous envelopes prior to accretion. This result has significant implications for the survival of neutral tidal structures, and suggests that current and future high-ion UV absorption-line observations are indicative of warm coronae surrounding LMC-mass dwarfs, independent of their environment.

Compact Obscured Nuclei (CONs) are heavily obscured infrared cores that have been found in local (ultra)luminous infrared galaxies (U/LIRGs). They show bright emission from vibrationally excited rotational transitions of HCN, known as HCN-vib, and are thought to harbor Compton Thick (CT, $N_{\text{H}} \geq 10^{24}$ cm$^{-2}$) active galactic nuclei (AGN) or extreme compact starbursts. We explore the potential evolutionary link between CONs and CT AGN by searching for CONs in hard X-ray-confirmed CT AGN from the Great Observatories All-sky LIRG Survey (GOALS). Here, we present new Atacama Large Millimeter/submillimeter Array Band 6 observations that targeted HCN-vib emission in four hard X-ray-confirmed CT AGN. We analyze these objects together with literature HCN-vib measurements of five additional hard X-ray-confirmed CT AGN from the GOALS sample. We do not detect any CONs in this combined sample of nine CT AGN. We then explore a proposed evolutionary sequence in which CONs evolve into X-ray-detectable CT AGN once outflows and feedback reduce the column densities of the enshrouding gas. We find, however, no evidence of well-developed dense molecular outflows in the observed CT AGN. While this could suggest that CT AGN are not universally linked to CONs, it could also be explained by a short duty cycle for molecular outflows.

Weak gravitational lensing is the slight distortion of galaxy shapes caused primarily by the gravitational effects of dark matter in the universe. In our work, we seek to invert the weak lensing signal from 2D telescope images to reconstruct a 3D map of the universe's dark matter field. While inversion typically yields a 2D projection of the dark matter field, accurate 3D maps of the dark matter distribution are essential for localizing structures of interest and testing theories of our universe. However, 3D inversion poses significant challenges. First, unlike standard 3D reconstruction that relies on multiple viewpoints, in this case, images are only observed from a single viewpoint. This challenge can be partially addressed by observing how galaxy emitters throughout the volume are lensed. However, this leads to the second challenge: the shapes and exact locations of unlensed galaxies are unknown, and can only be estimated with a very large degree of uncertainty. This introduces an overwhelming amount of noise which nearly drowns out the lensing signal completely. Previous approaches tackle this by imposing strong assumptions about the structures in the volume. We instead propose a methodology using a gravitationally-constrained neural field to flexibly model the continuous matter distribution. We take an analysis-by-synthesis approach, optimizing the weights of the neural network through a fully differentiable physical forward model to reproduce the lensing signal present in image measurements. We showcase our method on simulations, including realistic simulated measurements of dark matter distributions that mimic data from upcoming telescope surveys. Our results show that our method can not only outperform previous methods, but importantly is also able to recover potentially surprising dark matter structures.

Recent measurements of baryon acoustic oscillations (BAO) by the Dark Energy Spectroscopic Instrument (DESI) exhibit a mild-to-moderate tension with cosmic microwave background (CMB) and Type Ia supernova (SN) observations when interpreted within a flat $\Lambda$CDM framework. This discrepancy has been cited as evidence for dynamical dark energy (DDE). Given the profound implications of DDE for fundamental physics, we explore whether the tension can instead be resolved by modifying the physics of recombination. We find that a phenomenological model of modified recombination can effectively reconcile the BAO and CMB datasets and, unlike DDE, also predicts a higher Hubble constant $H_0$, thereby partially alleviating the Hubble tension. A global fit to BAO, CMB, and calibrated SN data clearly favors modified recombination over DDE, suggesting that current claims of a DDE detection may be premature.

The Roman Space Telescope Galactic Bulge Time Domain Survey (GBTDS) is expected to detect ~10^5 transiting planets. Many of these planets will have short orbital periods and are thus susceptible to tidal decay. We use a catalog of simulated transiting planet detections to predict the yield of orbital decay detections in the Roman GBTDS. Assuming a constant stellar tidal dissipation factor, Q^{'}_{*}, of 10^6, we predict ~ 5 - 10 detections. We additionally consider an empirical period-dependent parameterization of Q^{'}_{*} \propto P^{-3} and find a substantially suppressed yield. We conclude that Roman will provide constraints on the rate of planet engulfment in the Galaxy and probe the physics of tidal dissipation in stars.

Galaxy formation models, particularly semi-analytic models (SAMs), rely on differential equations with free parameters to describe the physical mechanisms governing galaxy formation and evolution. Traditionally, most SAMs calibrate these parameters manually to match observational data. However, this approach fails to fully explore the multidimensional parameter space, resulting in limited robustness and inconsistency with some observations. In contrast, the L-Galaxies SAM features a unique Markov Chain Monte Carlo (MCMC) mode, enabling robust model calibration. Using this functionality, we address a long-standing tension in galaxy formation models: simultaneously reproducing the number densities of dusty star-forming galaxies (DSFGs) and high-redshift massive quiescent galaxies (MQs). We test nine combinations of observational constraints - including stellar mass functions, quiescent fractions, neutral hydrogen mass functions, and DSFG number densities - across different redshifts. We then analyze the resulting galaxy property predictions and discuss the underlying physical mechanisms. Our results identify a model that reasonably matches the number density of DSFGs while remaining consistent with observationally-derived lower limits on the number density of high-redshift MQs. This model requires high star formation efficiencies in mergers and a null dependency of supermassive black hole (SMBH) cold gas accretion on halo mass, facilitating rapid stellar mass and SMBH growth. Additionally, our findings highlight the importance of robust calibration procedures to address the significant degeneracies inherent to multidimensional galaxy formation models.

How momentum, energy, and magnetic fields are transported in the presence of macroscopic gradients is a fundamental question in plasma physics. Answering this question is especially challenging for weakly collisional, magnetized plasmas, where macroscopic gradients influence the plasma's microphysical structure. In this paper, we introduce thermodynamic forcing, a new method for systematically modeling how macroscopic gradients in magnetized or unmagnetized plasmas shape the distribution functions of constituent particles. In this method, we propose to apply an anomalous force to those particles inducing the anisotropy that would naturally emerge due to macroscopic gradients in weakly collisional plasmas. We implement thermodynamic forcing in particle-in-cell (TF-PIC) simulations using a modified Vay particle pusher and validate it against analytic solutions of the equations of motion. We then carry out a series of simulations of electron-proton plasmas with periodic boundary conditions using TF-PIC. First, we confirm that the properties of two electron-scale kinetic instabilities -- one driven by a temperature gradient and the other by pressure anisotropy -- are consistent with previous results. Then, we demonstrate that in the presence of multiple macroscopic gradients, the saturated state can differ significantly from current expectations. This work enables, for the first time, systematic and self-consistent transport modeling in weakly collisional plasmas, with broad applications in astrophysics, laser-plasma physics, and inertial confinement fusion.

We report the bulk soil electrical conductivity and relative permittivity at a site in the Canadian High Arctic (79.37980 degrees N, 90.99885 degrees W). The soil parameters are determined using impedance measurements of a dipole antenna mounted horizontally 52 cm above the surface. The antenna is part of the Mapper of the IGM Spin Temperature (MIST) radio cosmology experiment. The measurements were conducted on July 17-28, 2022, every 111 minutes, and in the frequency range 25-125 MHz. To estimate the soil parameters, we compare the impedance measurements with models produced from numerical electromagnetic simulations of the antenna, considering single- and two-layer soil models. Our best-fit soil model corresponds to a two-layer model in which the electrical parameters are consistent with unfrozen soil at the top and frozen soil underneath. The best-fit parameters further agree with measurements done at other Arctic sites with more traditional techniques, such as capacitively-coupled resistivity, electrical resistivity tomography, and ground-penetrating radar.

Black holes are the sources of the strongest gravitational fields that can be found today in the Universe and are ideal laboratories for testing Einstein's theory of General Relativity in the strong field regime. In this letter, I show that the possibility of an interstellar mission to send small spacecrafts to the nearest black hole, although very speculative and extremely challenging, is not completely unrealistic. Certainly we do not have the necessary technology today, but it may be available in the next 20-30 years. The mission may last 80-100 years, but we would be able to obtain very valuable information about black holes and General Relativity that could be unlikely obtained in other ways.

As a high-level discipline, the development of remote sensing depends on the contribution of many other basic and applied disciplines and technologies. For example, due to the close relationship between remote sensing and photogrammetry, remote sensing would inevitably integrate disciplines such as optics and color science. Also, remote sensing integrates the knowledge of electronics in the conversion from optical signals to electrical signals via CCD (Charge-Coupled Device) or other image sensors. Moreover, when conducting object identification and classification with remote sensing data, mathematical morphology and other digital image processing technologies are used. These examples are only the tip of the iceberg of interdisciplinary integration of remote sensing. This work briefly reviews the interdisciplinary integration of remote sensing with four examples - ecology, mathematical morphology, machine learning, and electronics.

We study the possibility of probing leptogenesis via stochastic gravitational waves (GW) arising from a dark sector assisted first-order electroweak phase transition. The same dark sector, with non-trivial transformation under an unbroken $Z_2$ symmetry is also responsible for providing the only source of CP asymmetry via one-loop interference with the tree level decay of a heavy right-handed neutrino into lepton and Higgs doublets. The new Yukawa and scalar portal couplings enhance the CP asymmetry allowing TeV scale leptogenesis without any resonant enhancement. Light neutrino masses arise from a combination of type-I and one-loop contributions with vanishing lightest neutrino mass. While the new degrees of freedom in sub-TeV range keep the detection prospects at terrestrial experiments promising, the new scalars enhance the strength of the electroweak phase transition keeping the GW signals within reach of near future experiments like LISA.

Superheavy dark matter has been attractive as a candidate of particle dark matter. We propose a "natural" particle model, in which the dark matter serves as the inflaton in natural inflation, while decaying to high-energy particles at energies of $10^{9}-10^{13} \, \text{GeV}$ from the prediction of the inflation. A scalar field responsible to dilute the dark matter abundance revives the natural inflation. Since the dark matter must be a spin zero scalar, we study carefully the galactic dark matter 3-body decay into fermions and two body decays into a gluon pair, and point out relevant multi-messenger bounds that constrain these decay modes. Interestingly, the predicted energy scale may coincide with the AMATERASU event and/or the KM3NeT neutrino event, KM3-230213A. We also point out particle models with dark baryon to further alleviate $\gamma$-ray bounds. This scenario yields several testable predictions for the UHECR observations, including the highest-energy neutrons that are unaffected by magnetic fields, the tensor-to-scalar ratio, the running of spectral indices, $\alpha_s\gtrsim\mathcal{O}(0.001)$, and the existence of light new colored particles that could be accessible at future collider experiments. Further measurements of high energy cosmic rays, including their components and detailed coordinates may provide insight into not only the origin of the cosmic rays but also inflation.