Papers for Monday, Jan 29 2024

The nearby long-duration GRB 191019A recently detected by {\em Swift} lacks an associated supernova and belongs to a host galaxy with little star formation activity, suggesting that the origin of this burst is the result of the merger of two compact objects with dynamical interactions in a high-density medium of an active galactic nucleus (AGN). By motivating of this event, it occurs in such high-density environment, the ejecta- circumstellar medium (CSM) interaction cannot be ignored to contribute a possible kilonova emission. Here, we theoretically calculate the kilonova emission by considering the contribution of the ejecta-CSM interaction in the high-density environment. We find that the contribution to the kilonova emission from the ejecta-CSM interaction will dominate at a later time, and a smaller ejecta mass will have a stronger kilonova emission from the ejecta-CSM interaction. Moreover, we try to apply it to GRB 191019A, but we find that it is difficult to identify the possible kilonova emission from the observations due to the contribution of the bright host galaxy. On the other hand, it requires a less injected mass (less than $M_{\rm ej}=2\times10^{-5}M_{\odot}$) if one can detect the kilonova emission associated with a GRB 191019A-like event in the future. The r-process-powered and spin energy contributions from the magnetar are also discussed.

Considering ten well-known relativistic mean field models, we invoke feeble interaction between hadronic matter and fermionic dark matter (DM) $\chi$ via new physics scalar ($\phi$) and vector ($\xi$) mediators in neutron star core, thereby forming DM admixed neutron stars (DMANSs). The chosen masses of the DM fermion ($m_{\chi}$) and the mediators ($m_{\phi}$ and $m_{\xi}$) are consistent with the self-interaction constraint from Bullet cluster while their respective couplings ($y_{\phi}$ and $y_{\xi}$) are also constrained by the present day relic abundance. Assuming that both $\phi$ and $\xi$ contribute equally to the relic abundance, we compute the equation of state of the DMANSs and consequently their structural properties. We found that for a particular (constant) DM density, the presence of lighter DM results in more massive DMANSs with larger radius. In the light of the various recent constraints like those from the massive pulsar PSR J0740+6620, the gravitational wave (GW170817) data and the results of NICER experiments for PSR J0030+0451 and PSR J0740+6620, we provide a bound on $m_{\chi}$ within the framework of the present work as $m_{\chi}\approx$ (0.1 $-$ 30) GeV for a wide range of fixed DM Fermi momenta $k_F^{\chi}$=(0.01 $-$ 0.06) GeV. In the case of the hadronic models that yield larger radii corresponding to the low mass neutron stars in the no-DM scenario, interaction with comparatively heavier DM fermion is necessary in order to ensure that the DMANSs obtained with such models satisfy the radius constraints from both GW170817 and NICER data for PSR J0030+0451.

Evidence for a stochastic gravitational wave (GW) background, plausibly originating from the merger of supermassive black holes (SMBHs), is accumulating with observations from pulsar timing arrays. An outstanding question is how inspiraling SMBHs get past the "final parsec" of separation, where they have a tendency to stall before GW emission alone can make the binary coalesce. We argue that dynamical friction from the dark matter (DM) spike surrounding the black holes is sufficient to resolve this puzzle, if the DM has a self-interaction cross section of order $0.2-5\,$cm$^2$/g. The same effect leads to a softening of the GW spectrum at low frequencies as suggested by the current data. For collisionless cold DM, the friction deposits so much energy that the spike is disrupted and cannot bridge the final parsec, while for self-interacting DM, the isothermal core of the halo can act as a reservoir for the energy liberated from the SMBH orbits. A realistic velocity dependence, such as generated by the exchange of a massive mediator like a dark photon, is favored to give a good fit to the GW spectrum while providing a large enough core. A similar velocity dependence has been advocated for solving the small-scale structure problems of cold DM.

We have conducted a systematic search around the Milky Way (MW) analog NGC 253 (D=3.5 Mpc), as a part of the Panoramic Imaging Survey of Centaurus and Sculptor (PISCeS) - a Magellan+Megacam survey to identify dwarfs and other substructures in resolved stellar light around MW-mass galaxies outside of the Local Group. In total, NGC 253 has five satellites identified by PISCeS within 100 kpc with an absolute V-band magnitude $M_V<-7$. We have additionally obtained deep Hubble Space Telescope imaging of four reported candidates beyond the survey footprint: Do III, Do IV, and dw0036m2828 are confirmed to be satellites of NGC 253, while SculptorSR is found to be a background galaxy. We find no convincing evidence for the presence of a plane of satellites surrounding NGC 253. We construct its satellite luminosity function, which is complete down to $M_V$$\lesssim$$-8$ out to 100 kpc and $M_V$$\lesssim$$-9$ out to 300 kpc, and compare it to those calculated for other Local Volume galaxies. Exploring trends in satellite counts and star-forming fractions among satellite systems, we find relationships with host stellar mass, environment, and morphology, pointing to a complex picture of satellite formation, and a successful model has to reproduce all of these trends.

We estimate the mass of the Large Magellanic Cloud (LMC) using the kinematics of 30 LMC globular clusters (GCs). We combine proper motions (PMs) measured with HST, Gaia, or a combination of the two, from a recent study by Bennet et al. (2022) with literature line-of-sight velocities (LOSVs) to give 3 components of motion. With these, we derive a 3D velocity dispersion anisotropy $\beta = -0.72 ^{+0.62} _{-1.07}$, consistent with the GCs forming a flattened system with significant azimuthal motion. We then apply a tracer mass estimator and measure an enclosed mass $M (<13.2 \mathrm{kpc})= 2.66^{+0.42} _{-0.36} \times 10^{10} \mathrm{M}_\odot$. This is broadly consistent with results from previous studies of the LOSVs of GCs and other luminous tracers. Assuming a cosmologically-constrained NFW distribution for the dark matter, this implies a virial mass $M_\mathrm{virial} = 1.80^{+1.05} _{-0.54} \times 10^{11} \mathrm{M}_\odot$. Despite being an extrapolation by almost an order of magnitude in radius, this result is consistent with published estimates from other methods that are directly sensitive to the LMC's total mass. Our results support the conclusion that the LMC is approximately 17$^{+10}_{-6}$% of the Milky Way's mass, making it a significant contributor to the Local Group (LG) potential.

Exoplanets with short orbit period reside very close to their host stars. They transition very rapidly between different sectors of the circumstellar space environment along their orbit, leading to large variations of the magnetic field in the vicinity of the planet on short timescales. This rapid change of the magnetic flux through the conducting and resistive layer of the planetary upper atmosphere may drive currents that dissipate in the form of Joule Heating. Here, we estimate the amount of Joule Heating dissipation in the upper atmosphere of Trappist-1e, and two hypothetical planets orbiting the Sun in close-in orbits. We find that the rapid orbital motion could drive a significant amount of atmospheric heating and could significantly affect the planetary atmosphere escape rate. Thus, the process should be accounted for when studying the long-term evolution of exoplanetary atmospheres.

Aims. In this study, we analyse magnetic field properties of a set of 15 global magnetohydrodynamics simulations of solar-type star dynamos conducted using the ASH code. Our objective is to enhance our understanding of these properties by comparing theoretical results to current observations, aiming finally to provide fresh insights into the field. Methods. We analyse rotational and magnetic properties as a function of various stellar parameters (mass, age, rotation rate) in a Sun in time approach in our extended set of 3D MHD simulations. To facilitate direct comparisons with stellar magnetism observations using various Zeeman-effect techniques, we decompose numerical data into vectorial spherical harmonics. Results. The comparison of the trends we found in our simulations set reveals a promising overall agreement with the observational context of stellar magnetism, enabling us to suggest a plausible scenario for the magneto-rotational evolution of solar-type stars. In particular, we find that the magnetic field may reach a minimum in amplitude at a transition value in Rossby number near unity. This may have important consequences on the long term evolution of solar-type stars, by impacting the relation between stellar age, rotation and magnetism. This supports the need for future observational campaigns, especially for stars in the high Rossby number regime.

During the pre-main-sequence (pre-MS) evolution stage of a star, significant amounts of stellar mass are accreted during episodic accretion events, such as multi-decade FUor-type outbursts. Here, we present a near-infrared spectroscopic follow-up study of 33 high-amplitude (most with $\Delta K_s$ > 4 mag) variable sources discovered by the Vista Variables in the Via Lactea (VVV) survey. Based on the spectral features, 25 sources are classified as eruptive young stellar objects (YSOs), including 15 newly identified FUors, six with long-lasting but EXor-like bursts of magnetospheric accretion and four displaying outflow-dominated spectra. By examining the photometric behaviours of eruptive YSOs, we found most FUor-type outbursts have higher amplitudes ($\Delta K_s$ and $\Delta W2$), faster eruptive timescales and bluer infrared colours than the other outburst types. In addition, we identified seven post-main sequence variables apparently associated with deep dipping events and an eruptive star with deep AlO absorption bands resembling those seen in the V838 Mon stellar merger.

We present an intensive multiwavelength monitoring campaign of the quasar PG 1302$-$102 with Swift and the Las Cumbres Observatory network telescopes. At $z\sim0.3$, it tests the limits of the reverberation mapping (RM) technique in probing the accretion disk around a supermassive black hole (SMBH) and extends the parameter space to high masses and high accretion rates. This is also the first time the RM technique has been applied to test disk structures predicted in the SMBH binary model that has been suggested for this source. PG 1302$-$102 was observed at a $\sim$daily cadence for $\sim 9$ months in 14 bands spanning from X-ray to UV and optical wavelengths, and it shows moderate to significant levels of variability correlated between wavelengths. We measure the inter-band time lags which are consistent with a $\tau \propto \lambda^{4/3}$ relation as expected from standard disk reprocessing, albeit with large errors. The disk size implied by the lag spectrum is consistent with the expected disk size for its black hole mass within uncertainties. While the source resembles other reverberation-mapped AGN in many respects, and we do not find evidence supporting the prevalent hypothesis that it hosts an SMBH binary, we demonstrate the feasibility of studying SMBH binaries from this novel angle and suggest possibilities for the LSST Deep Drilling Fields.

The red hypergiant VY CMa is famous for its very visible record of high mass loss events. Recent CO observations with ALMA revealed three previously unknown large scale outflows (Paper I). In this paper we use the CO maps to investigate the motions of a cluster of four clumps close to the star, not visible in the optical or infrared images. We present their proper motions measured from two epochs of ALMA images and determine the line of sight velocities of the gas in emission at the clumps. We estimate their masses and ages, or time since ejection, and conclude that all four were ejected during VY CMa's active period in the early 20th century. Together with two additional knots observed with HST, VY CMa experienced at least six massive outflows during a 30 year period with a total mass lost greater than 0.07 Msun. The position-velocity map of the $^{12}$CO emission reveals previously unnoticed attributes of the older outer ejecta. In a very narrow range of Doppler velocities, $^{12}$CO absorption and emission causes some of this outer material to be quite opaque. At those frequencies the inner structure is hidden and we see only emission from an extended outer region. This fact produces a conspicuous but illusory dark spot if one attempts to subtract the continuum in a normal way.

During the pre-main-sequence evolution, Young Stellar Objects (YSOs) assemble most of their mass during the episodic accretion process. The rarely seen FUOr-type events (FUOrs) are valuable laboratories to investigate the outbursting nature of YSOs. Here, we present multi-wavelength detection of a high-amplitude eruptive source in the young open cluster VdBH 221 with an ongoing outburst, including optical to mid-infrared time series and near-infrared spectra. The initial outburst has an exceptional amplitude of $>$6.3 mag in Gaia and 4.6 mag in $K_s$, with a peak luminosity up to 16 $L_{\odot}$ and a peak mass accretion rate of 1.4 $\times$ 10$^{-5}$ $M_\odot$ yr$^{-1}$. The optical to infrared spectral energy distribution (SED) of this object is consistent with a low-mass star (0.2$M_\odot$) with a modest extinction ($A_V < 2$ mag). A 100-d delay between optical and infrared rising stages is detected, suggesting an outside-in origin of the instability. The spectroscopic features of this object reveal a self-luminous accretion disc, very similar to FU Orionis, with a low line-of-sight extinction. Most recently, there has been a gradual increase in brightness throughout the wavelength range, possibly suggesting an enhancement of the mass accretion rate.

We have performed a comprehensive search of a VISTA Variables in the Via Lactea (VVV) database of 9.5 yr light curves for variable sources with $\Delta K_s \ge 4$ mag, aiming to provide a large sample of high amplitude eruptive young stellar objects (YSOs) and detect unusual or new types of infrared variable source. We find 222 variable or transient sources in the Galactic bulge and disc, most of which are new discoveries. The sample mainly comprises novae, YSOs, microlensing events, Long Period Variable stars (LPVs) and a few rare or unclassified sources. Additionally, we report the discovery of a significant population of aperiodic late-type giant stars suffering deep extinction events, strongly clustered in the Nuclear Disc of the Milky Way. We suggest that these are metal-rich stars in which radiatively driven mass loss has been enhanced by super-solar metallicity. Among the YSOs, 32/40 appear to be undergoing episodic accretion. Long-lasting YSO eruptions have a typical rise time of $\sim$2 yr, somewhat slower than the 6-12 month timescale seen in the few historical events observed on the rise. The outburst durations are usually at least 5 yr, somewhat longer than many lower amplitude VVV events detected previously. The light curves are diverse in nature, suggesting that multiple types of disc instability may occur. Eight long-duration extinction events are seen wherein the YSO dims for a year or more, attributable to inner disc structure. One binary YSO in NGC 6530 displays periodic extinction events (P=59 days) similar to KH 15D.

Episodic accretion is one of the competing models to explain the observed luminosity spread in young stellar clusters. These short-lived high accretion events could also have a strong impact on planet formation. Observations of high-amplitude variability in young stellar objects (YSOs) due to large changes in the accretion rate provide direct observational evidence for episodic accretion. However, there are still uncertainties in the frequency of these events and if episodic accretion is universal among YSOs. To determine the frequency of outbursts in Class I YSOs, we built a large and robust sample of objects at this evolutionary stage, and searched for high-amplitude near-infrared ($\Delta K_{\rm S}>2$~mag) variability in the VIRAC2 database of the Vista Variables in the Via Lactea (VVV) survey. By complementing with near-IR (2MASS and DENIS) and mid-IR (WISE/Neo-WISE) data, we find that from $\sim$ 7000 Class I YSOs, 97 objects can be classified as eruptive variable YSOs. The duration of the outbursts vary from a few months to longer than 9 years, and cover a similar range of amplitudes. Values of $\Delta K_{\rm S}>5$~mag, however, are only observed in outbursts with duration longer than 9 years. When considering different effects of completeness and contamination we estimate that the incidence of episodic accretion in Class I YSOs is between 2\% and 3\%. Finally, we determine a recurrence timescale of long-term outbursts (a.k.a FUors) of $\tau=1.75^{+1.12}_{-0.87}$~kyr. The latter value agrees with previous estimates and is in line with the expectations of higher frequency of FUor outbursts during younger stages of evolution.

The Far- and Lyman-Ultraviolet Imaging Demonstrator (FLUID) is a rocket-borne arcsecond-level ultraviolet (UV) imaging instrument covering four bands between 92 -- 193 nm. FLUID will observe nearby galaxies to find and characterize the most massive stars that are the the primary drivers of the chemical and dynamical evolution of galaxies, and the co-evolution of the surrounding galactic environment. The FLUID short wave channel is designed to suppress efficiency at Lyman-$\alpha$ (121.6 nm), while enhancing the reflectivity of shorter wavelengths. Utilizing this technology, FLUID will take the first ever images of local galaxies isolated in the Lyman ultraviolet (90 -- 120 nm). As a pathfinder instrument, FLUID will employ and increase the TRL of band-selecting UV coatings, and solar-blind UV detector technologies including microchannel plate and solid state detectors; technologies prioritized in the 2022 NASA Astrophysical Biennial Technology Report. These technologies enable high throughput and high sensitivity observations in the four co-aligned UV imaging bands that make up the FLUID instrument. We present the design of FLUID, status on the technology development, and results from initial assembly and calibration of the FLUID instrument.

This paper explores the effects of numerical algorithms on global magnetohydrodynamics (MHD) simulations of solar wind (SW) in the inner heliosphere. To do so, we use sunRunner3D, a 3-D MHD model that employs the boundary conditions generated by CORHEL and the PLUTO code to compute the plasma properties of the SW with the ideal MHD approximation up to 1.1 AU in the inner heliosphere. Mainly, we define three different combinations of numerical algorithms based on their diffusion level. This diffusion level is related to the way of solving the MHD equations using the finite volume formulation, and, therefore, we set in terms of the divergence-free condition methods, Riemann solvers, variable reconstruction schemes, limiters, and time-steeping algorithms. According to the simulation results, we demonstrate that sunRunner3D reproduces global features of Corotating Interaction Regions (CIRs) observed by Earth-based spacecraft (OMNI) for a set of Carrington rotations that cover a period that lays in the late declining phase of solar cycle 24, independently of the numerical algorithms. Moreover, statistical analyses between models and in-situ measurements show reasonable agreement with the observations, and remarkably, the high diffusive method matches better with in-situ data than low diffusive methods.

Self-organised critical avalanche models are a class of cellular automata that, despite their simplicity, can be applied to the modeling of solar (and stellar) flares and generate robust power-law distributions in event size measures. However, bridging the conceptual gap to both magnetohydrodynamics and real flare observations continues to prove challenging. In this paper, we focus on a specific, key aspect of this endeavour, namely the definition of magnetic energy and its consequences for the model's internal dynamics and energy release statistics. We show that the dual requirement of releasing energy and restoring local stability demands that the instability criterion and boundary conditions be set in a manner internally consistent with a given energy definition, otherwise unphysical behavior ensues, e.g., negative energy release. Working with three energy definitions previously used in the literature, we construct such internally consistent avalanche models and compare/contrast their energy release statistics. Using the same set of models, we also explore a recent proposal by Farhang et al. (2018, 2019), namely that avalanches/flares should maximize the amount of energy released by the lattice when instabilities are triggered. This tends to produce avalanches of shorter duration but higher peak energy release, but adding up to similar total energy release. For the three energy definition we tested, such avalanche models exhibit almost identical distributions of event size measures. Our results indicate that the key to reproduce solar-like power-law slopes in these size measures is lattice configurations in which most nodes remain relatively far from the instability threshold.

A sample of SiO-masing late-type stars located in the inner Galaxy is analyzed with the goal of better constraining their obscuration. This reference sample allows us to define mathematical relations between their dereddened infrared colours and the observed colours (e.g. \Ks-[8], \Ks-[24]). The derived equations define a property (a locus) of these late-type stars. Therefore, they enable us to derive the interstellar extinction. With estimated spectral types, it is possible to decompose the total extinction in the two components (interstellar and envelope extinction). These relations are very useful to classify extremely obscured late-type stars located in the inner Galaxy. Estimating the two extinction components is performable on an individual late-type star, independently of its surrounding, and also when a few mid-infrared measurements are available.

Tight relationships exist in the local universe between the central stellar properties of galaxies and the mass of their supermassive black hole. These suggest galaxies and black holes co-evolve, with the main regulation mechanism being energetic feedback from accretion onto the black hole during its quasar phase. A crucial question is how the relationship between black holes and galaxies evolves with time; a key epoch to probe this relationship is at the peaks of star formation and black hole growth 8-12 billion years ago (redshifts 1-3). Here we report a dynamical measurement of the mass of the black hole in a luminous quasar at a redshift of 2, with a look back time of 11 billion years, by spatially resolving the broad line region. We detect a 40 micro-arcsecond (0.31 pc) spatial offset between the red and blue photocenters of the H$\alpha$ line that traces the velocity gradient of a rotating broad line region. The flux and differential phase spectra are well reproduced by a thick, moderately inclined disk of gas clouds within the sphere of influence of a central black hole with a mass of 3.2x10$^{8}$ solar masses. Molecular gas data reveal a dynamical mass for the host galaxy of 6x10$^{11}$ solar masses, which indicates an under-massive black hole accreting at a super-Eddington rate. This suggests a host galaxy that grew faster than the supermassive black hole, indicating a delay between galaxy and black hole formation for some systems.

We present the luminosity functions and stellar mass functions of supernova (SN) host galaxies and test if they differ from the functions of normal field galaxies. We utilize homogeneous samples consisting of 273 SNe Ia ($z\leq0.3$) and 44 core-collapse (CC) SNe ($z \leq 0.1$) from the Sloan Digital Sky Survey (SDSS) II Supernova Survey and the high-signal-to-noise-ratio photometry of galaxies from the Hyper Suprime-Cam Subaru Strategic Program (HSC SSP). SN hosts are classified into star-forming and passive galaxy groups based on the spectral energy distribution (SED) fitting. We find that the SN host luminosity functions and stellar mass functions deviate from those of normal field galaxies. Star-forming galaxies dominate the low-mass end of the SN Ia host mass function, while passive galaxies dominate the high-mass end. CC SNe are predominantly hosted by star-forming galaxies. In addition, intermediate-mass hosts produce CC SNe with the highest efficiency, while the efficiency of producing SNe Ia monotonically increases as the hosts become more massive. Furthermore, We derive the pseudo mass normalized SN rates (pSNuM) based on the mass functions. We find that the star-forming component of pSNuM$_{Ia}$ is less sensitive to the changes in stellar mass, in comparison with the total rate. The behavior of pSNuM$_{CC}$ suggests that the CC rate is proportional to the star-forming rate.

We investigate the diffuse light (DL) content of dark matter haloes in the mass range $11.5\leq \log M_{halo}\leq13$, a range that includes also the dark matter halo of the Milky-Way, taking advantage of a state-of-the-art semi-analytic model run on the merger trees extracted from a set of high-resolution cosmological simulations. The fraction of DL in such relatively small haloes is found to progressively decrease from the high to the low mass end, in good agreement with analytic (\citealt{purcell2007}) and numerical results from simulations (\citealt{proctor2023,ahvazi2023}), in good agreement also with the fraction of the DL observed in the Milky-Way (\citealt{deason2019}) and M31 (\citealt{harmsen2017}). Haloes with different masses have a different efficiency in producing DL: $\log M_{halo} \simeq 13$ is found to be the characteristic halo mass where the production of DL is the most efficient, while the overall efficiency decreases at both larger (\citealt{contini2024}) and smaller scales (this work). The DL content in this range of halo mass is the result of stellar stripping due to tidal interaction between satellites and its host (95\%) and mergers between satellites and the central galaxy (5\%), with pre-processed material, sub-channel of mergers and stripping and so already included in the 100\%, that contributes no more than 8\% on average. The halo concentration is the main driver of the DL formation: more concentrated haloes have higher DL fractions that come from stripping of more massive satellites in the high halo mass end, while dwarfs contribute mostly in the low halo mass end.

High-latitude ($|b|>30^{\circ}$) molecular clouds have virial parameters that exceed 1, but whether these clouds can form stars has not been studied systematically. Using JCMT SCUBA-2 archival data, we surveyed 70 fields that target high-latitude Planck galactic cold clumps (HLPCs) to find dense cores with density of $10^{5}$-$10^{6}$ cm$^{-3}$ and size of $<0.1$ pc. The sample benefits from both the representativeness of the parent sample and covering densest clumps at the high column density end ($>1\times10^{21}$ cm$^{-2}$). At an average noise rms of 15 mJy/beam, we detected Galactic dense cores in only one field, G6.04+36.77 (L183), while also identifying 12 extragalactic objects and two young stellar objects. Compared to the low-latitude clumps, dense cores are scarce in HLPCs. With synthetic observations, the densities of cores are constrained to be $n_c\lesssim10^5$ cm$^{-3}$, should they exist in HLPCs. Low-latitude clumps, Taurus clumps, and HLPCs form a sequence where a higher virial parameter corresponds to a lower dense core detection rate. If HLPCs were affected by the Local Bubble, the scarcity should favor turbulence-inhibited rather than supernova-driven star formation. Studies of the formation mechanism of the L183 molecular cloud are warranted.

The Tsinghua University-Ma Huateng Telescope for Survey (TMTS) has been constantly monitoring the northern sky since 2020 in search of rapidly variable stars. To find variable white dwarfs (WDs), the TMTS catalog is cross-matched with the WD catalog of Gaia EDR3, resulting in over 3000 light curves of WD candidates. The WD TMTS J17184064+2524314 (hereafter J1718) is the second ZZ~Ceti star discovered among these common sources. Based on the light curves from TMTS, follow-up photometric observations, and TESS, 10 periods and 3 combination periods are detected. A rotation period of $25.12\pm0.18$ hr is derived, according to the identified rotational splitting. Our spectroscopic observation indicates that this WD belongs to DA type with $T_{\rm eff}=11,670\pm604$ K, log $g=8.16\pm0.36$, $M = 0.70\pm0.23$ M$_{\odot}$, and age=$0.51\pm0.34$ Gyr. Based on core-parameterized asteroseismological model grids ($\geqslant$ 14 million), we derive a best-fit solution of $T_{\rm eff}=11,640\pm20$ K, log $g=8.267\pm0.008$, and $M = 0.750\pm0.005$ M$_{\odot}$ for J1718, consistent with the spectral fitting results. For this WD, the corresponding carbon and oxygen abundances in the core are 0.43 and 0.57, respectively. The distance derived from the intrinsic luminosity given by asteroseismology is $64\pm15$ pc, in accord with the distance of $70.1\pm0.2$ pc from Gaia DR3 within the uncertainties.

This paper presents a general analytical method to describe the center manifolds of collinear libration points in the Restricted Three-body Problem (RTBP). It is well-known that these center manifolds include Lissajous orbits, halo orbits, and quasihalo orbits. Previous studies have traditionally tackled these orbits separately by iteratively constructing high-order series solutions using the Lindstedt-Poincar\'e method. Instead of relying on resonance between their frequencies, this study identifies that halo and quasihalo orbits arise due to intricate coupling interactions between in-plane and out-of-plane motions. To characterize this coupling effect, a novel concept, coupling coefficient $\eta$, is introduced in the RTBP, incorporating the coupling term $\eta \Delta x$ in the $z$-direction dynamics equation, where $\Delta$ represents a formal power series concerning the amplitudes. Subsequently, a uniform series solution for these orbits is constructed up to a specified order using the Lindstedt-Poincar\'e method. For any given paired in-plane and out-of-plane amplitudes, the coupling coefficient $\eta$ is determined by the bifurcation equation $\Delta = 0$. When $\eta$ = 0, the proposed solution describes Lissajous orbits around libration points. As $\eta$ transitions from zero to non-zero values, the solution describes quasihalo orbits, which bifurcate from Lissajous orbits. Particularly, halo orbits bifurcate from planar Lyapunov orbits if the out-of-plane amplitude is zero. The proposed method provides a unified framework for understanding these intricate orbital behaviors in the RTBP.

M dwarfs are ubiquitous in our galaxy, and the rate at which they host stellar companions, and the properties of these companions, provides a window into the formation and evolution of the star(s), and of any planets that they may host. The Pervasive Overview of 'Kompanions' of Every M dwarf in Our Neighborhood (POKEMON) speckle survey of nearby M dwarfs is volume-limited from M0V through M9V out to 15 pc, with additional targets at larger distances. In total, 1125 stars were observed, and 455 of these are within the volume-limited, 15-pc sample of M-dwarf primaries. When we combine the speckle observations with known companions from the literature, we find that the stellar multiplicity rate of M dwarfs within 15 pc is 23.7% plus or minus 2.0%, and that the companion rate is 29.0% plus or minus 2.1%. We also find that the projected separation distribution for multiples that are known to host planets peaks at 198 au, while the distribution for multiples that are not yet known to host planets peaks at 5.57 au. This result suggests that the presence of close-in stellar companions inhibits the formation of M-dwarf planetary systems, similar to what has been found for FGK stars.

DEASA (Dayalbagh Educational Air Shower Array) consists of eight plastic scintillators each with an area of 1 square meter. The cosmic ray showers have been simulated in CORSIKA [1] for the different primary particles in the energy range of 1014- 1015 eV. The longitudinal and lateral profiles have been studied for Agra. The real-life applications of cosmic ray particles in space have been studied to protect the astronaut from the galactic cosmic rays [2]. A plastic scintillation detector is simulated in Geant4 to study applications in hadron and carbon ion therapy [3]. The proton and carbon beam are simulated through the tumour region to study the stopping power and depth dose distribution for different organs. The energy range for each study is optimized and the Bragg curve is then interpreted with Bragg peak position and range.

We conducted an investigation of the Coma cluster of galaxies by running a series of constrained hydrodynamic simulations with GIZMO-SIMBA and GADGET-3, based on initial conditions reconstructed from the SDSS survey volume in the ELUCID project. We compared simulation predictions and observations for galaxies, ICM and IGM in and around the Coma cluster to constrain galaxy formation physics. Our results demonstrate that this type of constrained investigation allows us to probe in more detail the implemented physical processes, because the comparison between simulations and observations is free of cosmic variance and hence can be conducted in a ''one-to-one'' manner. We found that an increase in the earlier star formation rate and the supernova feedback of the original GIZMO-SIMBA model is needed to match observational data on stellar, ISM and ICM metallicity. The simulations without AGN feedback can well reproduce the observational ICM electron density, temperature, and entropy profiles, ICM substructures, and the IGM temperature-density relation, while the ones with AGN feedback usually fail. However, one requires something like AGN feedback to reproduce a sufficiently large population of quiescent galaxies, particularly in low-density regions. The constrained simulations of the Coma cluster thus provide a test bed to understand processes that drive galaxy formation and evolution.

Hot rocky super-Earths are thought to be sufficiently irradiated by their host star to melt their surface and thus allow for long-lasting magma oceans. Some processes have been proposed for such planets to have retained primordial hydrogen captured during their formation while moving inward in the planetary system. The new generation of space telescopes such as the JWST may provide observations precise enough to characterize the atmospheres and perhaps the interiors of such exoplanets. We use a vaporization model that calculates the gas-liquid equilibrium between the atmosphere (including hydrogen) and the magma ocean, to compute the elemental composition of a variety of atmospheres for different quantities of hydrogen. The elemental composition is then used in a steady-state atmospheric model to compute the atmospheric structure and generate synthetic emission spectra. With this method, we confirm previous results showing that silicate atmospheres exhibit a thermal inversion, with notably an emission peak of SiO at 9~$\mu m$. We compare our method to the literature on the inclusion of hydrogen in the atmosphere, and show hydrogen reduces the thermal inversion, because of the formation of H2O which has a strong greenhouse potential. However planets that are significantly irradiated by their host star are sufficiently hot to dissociate H2O and thus also maintain a thermal inversion. The observational implications are twofold: 1) H2O is more likely to be detected in colder atmospheres; 2) Detecting a thermal inversion in hotter atmospheres does not a priori exclude the presence of H (in its atomic form). Due to the impact of H on the overall chemistry and atmospheric structure, and therefore observations, we emphasize the importance of including volatiles in the calculation of the gas-liquid equilibrium. Finally, we provide a criterion to determine potential targets for observation.

Solar-type stars form with a wide range of rotation rates. A wide range persists until a stellar age of 0.6 Gyr, after which solar-type stars exhibit Skumanich spin-down. Rotational evolution models incorporating polytropic stellar winds struggle to simultaneously reproduce these two regimes, namely the initially wide range and the Skumanich spin-down without imposing an a-priori cap on the wind mass-loss rate. We show that a three-dimensional wind model driven by Alfv\'en waves and observational data yields wind torques that agree with the observed age distribution of rotation rates. In our models of the Sun and twenty-seven open cluster stars aged from 0.04 to 0.6 Gyr that have observationally derived surface magnetic maps and rotation rates, we find evidence of exponential spin-down in young stars that are rapid rotators and Skumanich spin-down for slow rotators. The two spin-down regimes emerge naturally from our data-driven models. Our modelling suggests that the observed age distribution of stellar rotation rates arises as a consequence of magnetic field strength saturation in rapid rotators.

Starting from the Hamiltonian representation of the dynamics in \cite{rosengren2015chaos,colombo2019long}, this work proposes an innovative procedure to design fully-analytical maneuvers for post-mission disposal of HEOs satellites, exploiting the third-body perturbations. The Hamiltonian representation has been selected to include the external perturbing effects and to obtain a phase space representation. Notably, the orbit evolution can be described through the variation of double-averaged orbital elements over the orbital periods of the spacecraft and the perturbing bodies around the central planet. this work conveys a two-dimensional Hamiltonian representation under the third-body perturbations and the central planet's oblateness. The effect of solar radiation pressure has been neglected in this analysis.

W50 is a radio nebula around hyper-accreting Galactic microquasar SS433. In this letter, we focus on one peculiar feature of W50 - a pair of so-called "extended X-ray jets" (EXJs). These "jets" have a size of $\sim20\, {\rm pc}$, a sharp inner boundary, and their spectra are well represented by a featureless X-ray continuum. We argue that EXJ could be an outcome of a powerful {\it anisotropic} wind produced by a super-critical accretion disk. The wind itself consists of two components. The first component is a nearly isotropic outflow that subtends most of the solid angle as seen from the compact source. The second component is a more collimated wind aligned with the binary system rotation axis (polar wind). The termination shock of the former component recollimates the latter, giving it an appearance of an extended X-ray structure. In this model, the EXJ continuum spectrum is due to synchrotron emission of electrons accelerated in the recollimated polar wind. At variance with many other studies, in this model, the EXJ structures are not directly related to the highly collimated $0.26\;\!c$ baryonic jets. Instead, the EXJ and the W50's ears are produced by the part of the wind with an Eddington level kinetic luminosity confined to a half opening angle of $\sim 10$ degrees, which is not necessarily a recollimated version of the jets.

There are currently two open questions in white dwarf physics: why are massive dwarfs observed less often in astronomical surveys, and why have not any super-Chandrasekhar white dwarfs been found despite the discovery of more than a dozen peculiar, overly-luminous type Ia supernovae in about a couple of decades? According to different research, magnetic fields appear to somewhat resolve these issues, but stability remains a concern. For the first time, we investigate how modified gravity affects the specific heat of electrons and ions, the crystallization process, and the cooling mechanism in white dwarfs. We demonstrate it for the Ricci-based gravity. We show that massive white dwarfs fade faster and conclude that it could be a physical reason, apart from the presence of high magnetic fields, both for finding fewer massive white dwarfs and the lack of direct detection of super-Chandrasekhar white dwarfs.

Based on the horizon-scale magnetofluid model developed in [arXiv:2309.13304], we investigate the millimeter-wave images of a geometrically thick accretion disk or a funnel wall jet around a Kerr black hole. By employing the numerical method to solve the null geodesic and radiative transfer equations, we obtain the optical appearances at various observational angles and frequencies, generated by the thermal synchrotron radiation within the magnetofluid. For the thick disk, we specifically examine the impact of emission anisotropy on images, concluding that anisotropic synchrotron radiation could play an important role in the observability of the photon ring. For the funnel wall jet, we find that both the outflow and inflow funnel walls exhibit annular structures on the imaging plane. The outflow jet yields a brighter primary image than the photon ring, whereas the inflow jet does not. Based on our investigation, the inflow funnel wall model can not be ruled out by current observations of M87*.

Hierarchical black hole (BH) mergers are one of the most straightforward mechanisms to produce BHs inside and above the pair-instability mass gap. Here, we investigate the impact of globular cluster (GC) evolution on hierarchical mergers, and we account for the uncertainties related to BH mass pairing functions on the predicted primary BH mass, mass ratio and spin distribution. We find that the evolution of the host GC quenches the hierarchical BH assembly already at the third generation, mainly due to cluster expansion powered by a central BH sub-system. Hierarchical mergers match the primary BH mass distribution from GW events for $m_1 > 50 \, \mathrm{M_{\odot}}$, regardless of the assumed BH pairing function. At lower masses, however, different pairing functions lead to dramatically different predictions on the primary BH mass merger rate density. We find that the primary BH mass distribution evolves with redshift, with a larger contribution from mergers with $m_1 \geq 30 \, \mathrm{M_{\odot}}$ for $z\geq{}1$. Finally, we calculate the mixing fraction of BBHs from GCs and isolated binary systems. Our predictions are very sensitive to the spins, which favor a large fraction ($>0.6$) of BBHs born in GCs, in order to reproduce misaligned spin observations.

We performed the temperature diagnostics using inversions of data from nine spectroscopic observations obtained with the IRIS spectrograph in the MgII h & k lines. The sensitivity to the temperature of the emission peaks of these lines was exploited to determine the temperature of the coronal rain plasma using inversions of the spectroscopic profiles. Additional relationships between different spectral features of these lines, derived from the use of 3D radiative transfer line synthesis applied to simulations, were employed in order to derive the line-of-sight (LoS) velocities in different parts of the coronal rain plasma. For the first time, spectroscopic inversions of coronal rain were successfully performed. Temperatures derived from the inversions yield coronal rain clump temperatures at the formation height of the emission peaks of the MgII h & k lines in the range between 5000 and 7000 K. This narrow range of values remains consistent among all the different observations used in this work. We obtained LoS velocities of up to 40 km/s, which are consistent with the motion of the plasma being mostly constrained to the plane of the sky, as the coronal rain was mostly detected shortly after its formation and the observations took place in the disc. Furthermore, velocity diagnostics led to the detection of larger velocities at higher layers of the coronal rain plasma in some cases. This increased velocity seems to indicate that at some point (at least) during the fall of coronal rain clumps towards the chromosphere, the material in the upper part of the coronal rain plasma is falling with greater velocity than the material below it. The conditions of the temperature and density of the coronal rain plasma where the Mg II h line forms appear to be slightly different that those of the Mg II k line, with the former found at slightly colder and denser parts of the plasma.

High redshift blazars are among the most powerful non-explosive sources in the Universe and play a crucial role in understanding the evolution of relativistic jets. To understand these bright objects, we performed a detailed investigation of the multiwavelength properties of 79 $\gamma$-ray blazars with redshifts ranging from z = 2.0 to 2.5, using data from Fermi LAT, Swift XRT/UVOT, and NuSTAR observations. In the $\gamma$-ray band, the spectral analysis revealed a wide range of flux and photon indices, from $5.32 \times 10^{-10}$ to $3.40 \times 10^{-7}$ photons cm$^{-2}$ s$^{-1}$ and from 1.66 to 3.15, respectively, highlighting the diverse nature of these sources. The detailed temporal analysis showed that flaring activities were observed in 31 sources. Sources such as 4C+71.07, PKS 1329-049, and 4C+01.02, demonstrated significant increase in the $\gamma$-ray luminosity and flux variations, reaching peak luminosity exceeding $10^{50}$ erg s$^{-1}$. The temporal analysis extended to X-ray and optical/UV bands, showed clear flux changes in some sources in different observations. The time-averaged properties of high redshift blazars were derived through modeling the spectral energy distributions with a one-zone leptonic scenario, assuming the emission region is within the broad-line region (BLR) and the X-ray and $\gamma$-ray emissions are due to inverse Compton scattering of synchrotron and BLR-reflected photons. This modeling allowed us to constrain the emitting particle distribution, estimate the magnetic field inside the jet, and evaluate the jet luminosity, which is discussed in comparison with the disk luminosity derived from fitting the excess in the UV band.

The development of a corrector for low-frequency variations in the star image at the input of ESPriF, the echelle spectropolarimeter of the BTA primary focus, is reported. New technical solutions have made it possible to extend the operating frequency range to 10Hz for stars brighter than 13^m.

We present an investigation of evolved stars in the nearby star-forming galaxy WLM, using NIRCam imaging from the JWST resolved stellar populations early-release science (ERS) program. We find that various combinations of the F090W, F150W, F250M, and F430M filters can effectively isolate red supergiants (RSGs) and thermally-pulsing asymptotic giant branch (TP-AGB) stars from one another, while also providing a reasonable separation of the primary TP-AGB subtypes: carbon-rich C-type stars and oxygen-rich M-type stars. The classification scheme we present here agrees very well with the well-established Hubble Space Telescope (HST) medium-band filter technique. The ratio of C to M-type stars (C/M) is 0.8$\pm$0.1 for both the new JWST and the HST classifications, which is within one sigma of empirical predictions from optical narrow-band CN and TiO filters. The evolved star colors show good agreement with the predictions from the PARSEC$+$COLIBRI stellar evolutionary models, and the models indicate a strong metallicity dependence that makes stellar identification even more effective at higher metallicity. However, the models also indicate that evolved star identification with NIRCam may be more difficult at lower metallicies. We test every combination of NIRCam filters using the models and present additional filters that are also useful for evolved star studies. We also find that $\approx$90\% of the dusty evolved stars are carbon-rich, suggesting that carbonaceous dust dominates the present-day dust production in WLM, similar to the findings in the Magellanic Clouds. These results demonstrate the usefulness of NIRCam in identifying and classifying dust-producing stars without the need for mid-infrared data.

There is a common assumption in the particulate disc community that the pressure in particulate discs is essentially zero and that the disc streamlines follow Keplerian orbits, in the absence of self-gravity or external perturbations. It is also often assumed that the fluid description of particulate discs is not valid in the presence of crossing orbits (e.g. from nonzero free eccentricities). These stem from the misconception that fluid pressure arises due to the (typically rare) collisions between particles and that the velocity of particles in fluids are single-valued in space. In reality, pressure is a statistical property of the particle distribution function which arises precisely because there is a distribution of velocities at a given position. In this letter we demonstrate, with simple examples, that pressure in particulate discs is non-zero and is related to the inclination and free eccentricity distributions of the constituent particles in the discs. This means many common models of debris discs implicitly assume a nonzero, and potentially quite significant, dust pressure. We shall also demonstrate that the bulk motion of the dust is not the same as the particle motion and that the presence of pressure gradients can lead to strong departures from Keplerian motion.

Some post-common-envelope binaries are binary stars with short periods that exhibit significant period variations over long observational time spans. These eclipse timing variations (ETVs) are most likely to be accounted for by the presence of an unseen massive companion, potentially of planetary or substellar nature, and the light-travel time (LTT) effect. In this study, our main objective is to describe the diversity of compatible nontransit companions around PCEBs and explore the robustness of the solutions by employing tools for uncertainty estimation. We select the controversial data of the QS Vir binary star, which previous studies have suggested hosts a planet. We employ a minimizing strategy, using genetic algorithms to explore the global parameter space followed by refinement of the solution using the simplex method. We evaluate errors through the classical MCMC approach and discuss the error range for parameters. Our results highlight the strong dependence of ETV models for close binaries on the dataset used, which leads to relatively loose constraints on the parameters of the unseen companion. We find that the shape of the $O-C$ curve is influenced by the dataset employed. We propose an alternative method to evaluate errors on the orbital fits based on a grid search surrounding the best-fit values, obtaining a wider range of plausible solutions that are compatible with goodness-of-fit statistics. We also analyze how the parameter solutions are affected by the choice of the dataset, and find that this system continuously changes the compatible solutions as new data are obtained from eclipses. The best-fit parameters for QS Vir correspond to a low-mass stellar companion (57.71 $M_{jup}$ ranging from 40 to 64 $M_{jup}$) on an eccentric orbit ($e=0.91^{+0.07}_{-0.17}$) with a variety of potential periods ($P = 16.69 ^{+0.47}_{-0.42}$ yr.)

Context. The modelling of chemical transport mechanisms is crucial for accurate stellar characterizations. Atomic diffusion is one of these processes and it is commonly included in stellar models. However, it is usually neglected for F-type or more massive stars because it produces surface abundance variations that are unrealistic. Additional mechanisms to counteract atomic diffusion must therefore be considered. It has been demonstrated that turbulent mixing can prevent the surface abundance over-variations, and can also be calibrated to mimic the effects of radiative accelerations on iron. Aims. We aim to evaluate the effect of a calibrated turbulent mixing on the characterisation of a sample of F-type stars, and how the estimates compare with those obtained when the chemical transport mechanisms are neglected. Methods. We selected stars from two samples - one from the Kepler LEGACY sample and the other from a sample of Kepler planet-hosting stars. We inferred their stellar properties using two grids. The first grid considers atomic diffusion only in models that do not show chemical over-variations at the stellar surface. The second grid includes atomic diffusion in all the stellar models and the calibrated turbulent mixing to avoid unrealistic surface abundances. Results. Comparing the derived results from the two grids, we found that the results for the more massive stars in our sample will have higher dispersion in the inferred values of mass, radius and age, due to the absence of atomic diffusion in one of the grids. This can lead to relative uncertainties for individual stars of up to 5% for masses, 2% for radii and 20% for ages. Conclusions. This work shows that a proper modelling of the microscopic transport processes is key for an accurate estimation of their fundamental properties not only for G-type stars, but also for F-type stars.

Neutron stars are compact and dense celestial objects that offer the unique opportunity to explore matter and its interactions under conditions that cannot be reproduced elsewhere in the Universe. Their extreme gravitational, rotational and magnetic energy reservoirs fuel the large variety of their emission, which encompasses all available multi-messenger tracers: electromagnetic and gravitational waves, neutrinos, and cosmic rays. However, accurately measuring global neutron-star properties such as mass, radius, and moment of inertia poses significant challenges. Probing internal characteristics such as the crustal composition or superfluid physics is even more complex. This article provides a comprehensive review of the different methods employed to measure neutron-star characteristics and the level of reliance on theoretical models. Understanding these measurement techniques is crucial for advancing our knowledge of neutron-star physics. We also highlight the importance of employing independent methods and adopting a multi-messenger approach to gather complementary data from various observable phenomena as exemplified by the recent breakthroughs in gravitational-wave astronomy and the landmark detection of a binary neutron-star merger. Consolidating the current state of knowledge on neutron-star measurements will enable an accurate interpretation of the current data and errors, and better planning for future observations and experiments.

Precision measurements of astrometry and photometry require stringent control of systematics such as those arising from imperfect correction of sensor effects. In this work, we develop a parametric method to model the wavelength dependence of photo-response non-uniformity (PRNU) for a laser annealed backside-illuminated charge-coupled device. The model accurately reproduces the PRNU patterns of flat-field images taken at nine wavelengths from 290nm to 950nm, leaving the root mean square (RMS) residuals no more than 0.2% in most cases. By removing the large-scale non-uniformity in the flat fields, the RMS residuals are further reduced. This model fitting approach gives more accurate predictions of the PRNU than cubic-spline interpolation does with fewer free parameters. It can be applied to make PRNU corrections for individual objects according to their spectral energy distribution to reduce photometry errors.

The prompt emission mechanism of gamma-ray bursts (GRBs) is still unclear, and the time-resolved spectral analysis of GRBs is a powerful tool for studying their underlying physical processes. We performed a detailed time-resolved spectral analysis of 78 bright long GRB samples detected by Fermi/Gamma-ray Burst Monitor (GBM). A total of 1490 spectra were obtained and their properties were studied using a typical Band-shape model. Firstly, the parameter distribution of the time-resolved spectrum given as follows: the low-energy spectral index $\alpha \sim -0.72$, high-energy spectral index $\beta \sim -2.42$, the peak energy $E_{\rm p} \sim 221.69 \,\rm{keV}$, and the energy flux $F \sim 7.49\times 10^{-6} \rm{\, erg\,cm^{-2}\,s^{-1}}$. More than 80\% of the bursts exhibit the hardest low-energy spectral index $\alpha_{\rm max}$ exceeding the synchrotron limit (-2/3). Secondly, the evolution patterns of $\alpha$ and $E_{\rm p}$ were statistically analyzed. The results show that for multi-pulse GRBs the intensity-tracking pattern is more common than the hard-to-soft pattern in the evolution of both $E_{\rm p}$ and $\alpha$. The hard-to-soft pattern is generally shown in single-pulse GRBs or in the initial pulse of multi-pulse GRBs. Finally, we found a significant positive correlation between $F$ and $E_{\rm p}$, with half of the samples exhibiting a positive correlation between $F$ and $\alpha$. We discussed the spectral evolution of different radiation models. The diversity of spectral evolution patterns indicates that there may be more than one radiation mechanism occurring in the gamma-ray burst radiation process, including photospheric radiation and synchrotron radiation. However, it may also involve only one radiation mechanism, but more complicated physical details need to be considered.

The geometry of the neutral gas in and around galaxies is a key regulator of the escape of ionizing photons. We present the first statistical study aiming at linking the neutral and ionized gas distributions to the Lyman continuum (LyC) escape fraction (fesc(LyC)) in a sample of 22 confirmed LyC leakers and non-leakers at z~0.35 using the Keck Cosmic Web Imager (Keck/KCWI) and the Low Resolution Spectrograph 2 (HET/LRS2). Our integral field unit data enable the detection of neutral and low-ionization gas, as traced by Mg II, and ionized gas, as traced by [O II], extending beyond the stellar continuum for 7 and 10 objects, respectively. All but one object with extended Mg II emission also shows extended [O II] emission; in this case, Mg II emission is always more extended than [O II] by a factor 1.3 on average. Most of the galaxies with extended emission are non or weak LyC leakers (fesc(LyC) < 5%), but we find a large diversity of neutral gas configurations around these weakly LyC-emitting galaxies. Conversely, the strongest leakers (fesc(LyC) > 10%) appear uniformly compact in both Mg II and [O II] with exponential scale lengths <1 kpc. We also find a trend between fesc(LyC) and the spatial offsets of the nebular gas and the stellar continuum emission. Moreover, we find significant anti-correlations between the spatial extent of the neutral gas and the [O III]/[O II] ratio, and H$\beta$ equivalent width, as well as positive correlations with metallicity and UV size, suggesting that galaxies with more compact neutral gas sizes are more highly ionized. The observations suggest that strong LyC emitters do not have extended neutral gas halos and ionizing photons may be emitted in many directions. Combined with high ionization diagnostics, we propose the Mg II, and potentially [O II], spatial compactness are indirect indicators of LyC emitting galaxies at high-redshift.

In galaxy cluster collisions, the gas can be separated from dark matter halos. Abell $56$ displays signatures of a dissociative bullet-like merger with a possible high inclination angle between the plane of orbit and the sky. Our objective is to provide a comprehensive description of the features observed in the collision scenario of Abell $56$. Additionally, we aim to apply a potential weak lensing mass bias correction attributed to the merger to evaluate its impact on our findings. To investigate this, we perform tailored hydrodynamical $N$-body simulations, varying the impact parameter. We initially identified an early scenario at $0.12$ Gyr after the central passage that reproduces some observational features. However, the mean temperature of $9.8$ keV exceeded the observed value. Despite applying a mass bias correction due to the merger process, the new mean temperature of $8.4$ keV remained higher than the observed value. Our best model corresponds to the late scenario at $0.52$ Gyr after the pericenter, reproducing observed features of Abell $56$, with an inclination of $58^\circ$. These features include the offset between the main gas density peak and the south dark matter density peak of $103$ kpc, gas morphology, a line of sight relative velocity of $184$ km s$^{-1}$, and a mean temperature of $6.7$ keV. In the Abell $56$ collision scenario, the weak lensing mass bias did not significantly impact the overall dynamics of the cluster merger. The correction only resulted in a slight decrease in the final mean temperature.

We present observational evidence of late-time interaction between the ejecta of the hydrogen-poor Type Ib supernova (SN) 2019yvr and hydrogen-rich circumstellar material (CSM), similar to the Type Ib SN 2014C. A narrow H{\alpha} emission line appears simultaneously with a break in the light-curve decline rate at around 80-100 d after explosion. From the interaction delay and the ejecta velocity, under the assumption that the CSM is detached from the progenitor, we estimate the CSM inner radius to be located at ~6.5-9.1 {\times} 10^{15} cm. The H{\alpha} emission line persists throughout the nebular phase at least up to +420 d post-explosion, with a full width at half maximum of ~2000 km/s. Assuming a steady mass-loss, the estimated mass-loss rate from the luminosity of the H{\alpha} line is ~3-7 {\times} 10^{-5} M_\odot yr^{-1}. From hydrodynamical modelling and analysis of the nebular spectra, we find a progenitor He-core mass of 3-4 M{_\odot}, which would imply an initial mass of 13-15 M{_\odot}. Our result supports the case of a relatively low-mass progenitor possibly in a binary system as opposed to a higher mass single star undergoing a luminous blue variable phase.

We revisit the state of the light element abundances from big bang nucleosynthesis in early 2024 with particular focus on the derived baryon abundance. We find that the largest differences between the final baryon abundances are typically driven by the assumed Deuterium burning rates, characterized in this work by the underlying code. The rates from theoretical ab-initio calculations favor smaller baryon abundances, while experimentally-determined rates prefer higher abundances. Through robust marginalization over a wide range of nuclear rates, the recently released $\mathtt{PRyMordial}$ code allows for a conservative estimate of the baryon abundance at $\Omega_b h^2 = 0.02218 \pm 0.00055$ (using PDG-recommended light element abundances) in $\Lambda$CDM and $\Omega_b h^2 = 0.02196 \pm 0.00063$ when additional ultra-relativistic relics are considered ($\Lambda$CDM + $N_\mathrm{eff}$). These additional relics themselves are constrained to $\Delta N_\mathrm{eff} = -0.10 \pm 0.21$ by light element abundances alone.

We present the results obtained from detailed X-ray timing and spectral studies of X-ray pulsar Swift J0243.6+6124 during its giant and normal X-ray outbursts between 2017 and 2023 observed by the Neutron star Interior Composition Explorer (NICER). We focused on the timing analysis of the normal outbursts. A distinct break is found in the power density spectra of the source. The corresponding break frequency and slope of power-laws around the break vary with luminosity, indicating the change in accretion dynamics with mass accretion rate. Interestingly, we detected quasi-periodic oscillations within a specific luminosity range, providing further insights into the underlying physical processes. We also studied the neutron star spin period evolution and a luminosity variation in pulse profile during the recent 2023 outburst. The spectral analysis was conducted comprehensively for the giant and all other normal outbursts. We identified a double transition at luminosities of $\approx$7.5$\times$10$^{37}$ and 2.1$\times$10$^{38}$ erg s$^{-1}$ in the evolution of continuum parameters like photon index and cutoff energy with luminosity. This indicates three distinct accretion modes experienced by the source mainly during the giant X-ray outburst. A soft blackbody component with a temperature of 0.08-0.7 keV is also detected in spectra. The observed temperature undergoes a discontinuous transition when the pulsar evolves from a sub- to super-Eddington state. Notably, in addition to an evolving 6-7 keV iron line complex, a 1 keV emission line was observed during the super-Eddington state of the source, implying the X-ray reflection from the accretion disc or outflow material.

Primordial black hole formation has been discussed widely, when density perturbations in the early universe cause matter to collapse gravitationally, giving rise to these ultra-compact objects. We propose here that such a gravitational collapse also gives rise to primordial naked singularities, that would play an important role in the observable features of the present universe. We consider two types of collapse scenarios which would give rise to event-like and object-like visible singularities. We briefly discuss implications of primordial naked singularities, including those on dark matter, vis-a-vis primordial black holes.

Gravitational wave physics can probe theories that extend beyond General Relativity. Motivated by recent attention on the Kalb-Ramond field as a dark matter candidate, in this work, we study a parity violating dimension four operator which couples the dual Riemann curvature to the 2-form field. After mapping the equations of motion for the right- and left-handed gravitational wave amplitudes to the novel parameterization recently presented in \cite{Jenks:2023pmk}, we discuss various constraints on the model parameters in light of the coincident electromagnetic/gravitational wave signal of GW170817 and the GWTC-3 dataset.

The primary challenge in detecting ultrahigh energy (UHE) neutrinos with energies exceeding $10^{16}$ eV is to instrument a large enough volume to detect the extremely low flux, which falls as $\sim E^{-2}$. We explore in this article the feasibility of using the forest as a detector. Trees have been shown to be efficient broadband antennas, and may, without damage to the tree, be instrumented with a minimum of apparatus. A large scale array of such trees may be the key to achieving the requisite target volumes for UHE neutrino astronomy.

A set of equations are developed that extend the macroscale magnetic reconnection simulation model kglobal to include particle ions. The extension from earlier versions of kglobal, which included only particle electrons, requires the inclusion of the inertia of particle ions in the fluid momentum equation, which was not required in the electron-only model. The new equations will facilitate the exploration of the simultaneous non-thermal energization of ions and electrons during magnetic reconnection in macroscale systems. Numerical tests of the propagation of Alfv\'en waves in a plasma with anisotropic electron and ion pressure are presented to benchmark the new model.

We present a theoretical investigation of the expected experimental signals produced by freely falling atoms with time oscillating mass and transition frequency. These oscillations could be produced in a variety of models, in particular, models of scalar dark matter (DM) non universally coupled to the standard matter (SM) such as axion-like particles (ALP) and dilatons. Performing complete and rigorous calculations, we show that, on one hand, two different atomic species would accelerate at a different rate, and on the other hand, they would produce a non-zero differential phase shift in atom interferometers (AI). The former would produce observable signals in equivalence principle tests like the recent MICROSCOPE mission, and we provide a corresponding sensitivity estimate, showing that MICROSCOPE can reach beyond the best existing searches in the ALP case. We also compare the expected sensitivity of two future AI experiments, namely the AION-10 gradiometer and an isotope differential AI considered for MAGIS-100, that we will refer to as SPID. We show that the SPID setup would be more sensitive to these dark matter fields compared to the gradiometer one, assuming equivalent experimental parameters.

We propose a diagnostic tool for future analyses of stochastic gravitational wave background signals of extra-galactic origin in LISA data. Next-generation gravitational wave detectors hold the capability to track unresolved gravitational waves bundled into a stochastic background. This composite background contains cosmological and astrophysical contributions, the exploration of which offers promising avenues for groundbreaking new insights into very early universe cosmology as well as late-time structure formation. In this article, we develop a full end-to-end pipeline for the extraction of extra-galactic signals, based on kinematic anisotropies arising from the galactic motion, via full-time-domain simulations of LISA's response to the gravitational wave anisotropic sky. Employing a Markov-Chain-Monte-Carlo map-making scheme, multipoles up to $\ell=2$ are recovered for scale-free spectra that support an interpretation as signals originating from cosmic strings in the case of a high signal-to-noise ratio. We demonstrate that our analysis is consistently beating cosmic variance and is robust against statistical and systematic errors. The impact of instrumental noise on the extraction of kinematic anisotropies is investigated, and we establish a detection threshold of $\Omega_{GW}\gtrsim 5\times 10^{-8}$ in the presence of instrument-induced noise. Potential avenues for improvement in our methodology are highlighted.

We perform a numerical simulation of compression of a highly porous dust aggregate of monodisperse spheres. We find that the average interparticle normal force within the aggregate is inversely proportional to both the filling factor and the average coordination number, and we also derive this relation theoretically. Our findings would be applicable for granular matter of arbitrary structures, as long as the constituent particles are monodisperse spheres.