Papers for Thursday, Jan 11 2024

Complex organic molecules (COMs) have been detected ubiquitously in protostellar systems. However, at shorter wavelengths (~0.8mm) it is more difficult to detect larger molecules than at longer wavelengths (~3mm) because of the increase of millimeter dust opacity, line confusion, and unfavorable partition function. We aim to search for large molecules (>8 atoms) in the ALMA Band 3 spectrum of IRAS 16293-2422 B. We search for more than 70 molecules and identify as many lines as possible in the spectrum. The spectral settings were set to specifically target three-carbon species such as propanol and glycerol. We identify lines of 31 molecules including many oxygen-bearing COMs such as CH3OH and c-C2H4O and a few nitrogen- and sulfur-bearing ones such as HOCH2CN and CH3SH. The largest detected molecules are gGg-(CH2OH)2 and CH3COCH3. We do not detect glycerol or propanol but provide upper limits for them which are in line with previous laboratory and observational studies. The line density in Band 3 is only ~2.5 times lower in frequency space than in Band 7. From the detected lines in Band 3 at a $\gtrsim 6\sigma$ level, ~25-30% of them could not be identified indicating the need for more laboratory data of rotational spectra. We find similar column densities and column density ratios of COMs (within a factor ~2) between Band 3 and Band 7. The effect of dust optical depth for IRAS 16293-2422 B at an off-source location on column densities and column density ratios is minimal. Moreover, for warm protostars, long wavelength spectra are not only crowded, but also take longer integration times to reach the same sensitivity limit. The 3mm search has not yet resulted in detection of larger and more complex molecules in warm sources. A full deep ALMA Band 2-3 (i.e., 3-4 mm) survey is needed to assess whether low frequency data have the potential to reveal more complex molecules in warm sources.

We investigate Lyman-alpha (Ly$\alpha$) transmission spikes at $5.2 < z < 6.8$ using synthetic quasar spectra from the ``Cosmic Reionization On Computers" simulations. We focus on understanding the relationship between these spikes and the properties of the intergalactic medium (IGM). Disentangling the complex interplay between IGM physics and the influence of galaxies on the generation of these spikes presents a significant challenge. To address this, we employ Explainable Boosting Machines, an interpretable machine learning algorithm, to quantify the relative impact of various IGM properties on the Ly$\alpha$ flux. Our findings reveal that gas density is the primary factor influencing absorption strength, followed by the intensity of background radiation and the temperature of the IGM. Ionizing radiation from local sources (i.e. galaxies) appears to have a minimal effect on Ly$\alpha$ flux. The simulations show that transmission spikes predominantly occur in regions of low gas density. Our results challenge recent observational studies suggesting the origin of these spikes in regions with enhanced radiation. We demonstrate that Ly$\alpha$ transmission spikes are largely a product of the large-scale structure, of which galaxies are biased tracers.

TWA 27B (2M1207b) is the first directly-imaged planetary-mass (MP ~ 5 MJ) companion (Chauvin et al. 2004) and was observed at 0.9--5.3 micron with JWST/NIRSpec (Luhman et al. 2023). To understand the accretion properties of TWA 27B, we search for continuum-subtracted near-infrared helium and hydrogen emission lines and measure their widths and luminosities. We detect the He I triplet at 4.3 sigma and all Paschen-series lines covered by NIRSpec (Pa alpha, Pa beta, Pa gamma, Pa delta) at 4--5 sigma. The three brightest Brackett-series lines (Br alpha, Br beta, Br gamma) as well as Pf gamma and Pf delta are tentative detections at 2--3 sigma. We provide upper limits on the other hydrogen lines, including on H alpha through Hubble Space Telescope archival data. Three lines can be reliably deconvolved to reveal an intrinsic width Delta v = 67+-9 km/s, which is 60% of the surface freefall velocity. The line luminosities seem significantly too high to be due to chromospheric activity. Converting line luminosities to an accretion rate yields Mdot ~ 5e-9 MJ/yr when using scalings relationships for planetary masses, and Mdot ~ 0.1e-9 MJ/yr with extrapolated stellar scalings. Several of these lines represent first detections at an accretor of such low mass. The weak accretion rate implies that formation is likely over. This analysis shows that JWST can be used to measure low line-emitting mass accretion rates onto planetary-mass objects, motivates deeper searches for the mass reservoir feeding TWA 27B, and hints that other young directly-imaged objects might -- hitherto unbeknownst -- also be accreting.

Sub-mm surveys toward overdense regions in the early Universe are essential to uncover the obscured star-formation and the cold gas content of assembling galaxies within massive dark matter halos. In this work, we present deep ALMA mosaic observations covering an area of $\sim 2'\times2'$ around MQN01 (MUSE Quasar Nebula 01), one of the largest and brightest Ly-$\alpha$ emitting nebulae discovered thus far surrounding a radio-quiet quasar at $z\simeq3.25$. Our observations target the 1.2- and the 3-mm dust continuum, as well as the carbon monoxide CO(4-3) transition in galaxies in the vicinity of the quasar. We identify a robust sample of eleven CO line-emitting galaxies (including a closely-separated quasar companion) which lie within $\pm 4000\,{\rm km\,s^{-1}}$ relatively to the quasar systemic redshift. A fraction of these objects are missed in previous deep rest-frame optical/UV surveys thus highlighting the critical role of (sub-)mm imaging. We also detect a total of eleven sources revealed in their 1.2-mm dust continuum with six of them having either high-fidelity spectroscopic redshift information from rest-frame UV metal absorptions, or CO line which place them in the same narrow redshift range. A comparison of the CO luminosity function (LF) and 1.2-mm number count density with that of the general fields points to a galaxy overdensity of $\delta > 10$. We find evidence of a systematic flattening at the bright-end of the CO LF with respect to the trend measured in blank fields. Our findings reveal that galaxies in dense regions at $z\sim3$ are more massive and significantly richer in molecular gas than galaxies in fields, hence enabling a faster and accelerated assembly. This is the first of a series of studies to characterize one of the densest regions of the Universe found so far at $z > 3$.

Observations indicate dust populations vary between galaxies and within them, suggesting a complex life cycle and evolutionary history. Here we investigate the evolution of galactic dust populations across cosmic time using a suite of cosmological zoom-in simulations from the Feedback in Realistic Environments (FIRE) project, spanning $M_{\rm vir}=10^{9-12}M_{\odot};\,M_{*}=10^{6-11}\,M_{\odot}$. Our simulations incorporate a dust evolution model that accounts for the dominant sources of dust production, growth, and destruction and follows the evolution of specific dust species with set chemical compositions. All galactic dust populations in our suite exhibit similar evolutionary histories, with gas-dust accretion being the dominant producer of dust mass for all but the most metal-poor galaxies. The onset of efficient gas-dust accretion occurs above a `critical' metallicity threshold ($Z_{\rm crit}$). Due to this threshold, our simulations reproduce observed trends between galactic D/Z and metallicity and element depletion trends in the ISM. Furthermore, $Z_{\rm crit}$ varies between dust species due to differences in key element abundances, dust physical properties, and life cycle processes resulting in $Z_{\rm crit} \sim 0.05 Z_{\odot},\,0.2 Z_{\odot},\,0.5 Z_{\odot}$ for metallic iron, silicates, and carbonaceous dust, respectively. These variations could explain the lack of small carbonaceous grains observed in the Magellanic Clouds. We also find a delay between the onset of gas-dust accretion and when a dust population reaches equilibrium, which we call the equilibrium timescale ($\tau_{\rm eq}$). The relation between $\tau_{\rm eq}$ and the metal enrichment timescale of a galaxy, determined by its recent evolutionary history, can contribute to the scatter in the observed relation between galactic D/Z and metallicity.

We present a novel approach to measuring the expansion rate and the geometry of the Universe, which combine time-delay cosmography in lens galaxy clusters with pure samples of 'cosmic chronometers' (CCs) by probing the member galaxies. The former makes use of the measured time delays between the multiple images of time-varying sources strongly lensed by galaxy clusters, while the latter exploits the most massive and passive cluster member galaxies to measure the differential time evolution of the Universe. We applied two different statistical techniques, adopting realistic errors on the measured quantities, to assess the accuracy and the gain in precision on the values of the cosmological parameters. We demonstrate that the proposed combined method allows for a robust and accurate measurement of the value of the Hubble constant. In addition, this provides valuable information on the other cosmological parameters thanks to the complementarity between the two different probes in breaking parameter degeneracies. Finally, we showcase the immediate observational feasibility of the proposed joint method by taking advantage of the existing high-quality spectro-photometric data for several lens galaxy clusters.

We present high-definition observations with the James Webb Space Telescope of >1000 Cepheids in a geometric anchor of the distance ladder, NGC4258, and in 5 hosts of 8 SNe~Ia, a far greater sample than previous studies with JWST. These galaxies individually contain the largest samples of Cepheids, an average of >150 each, producing the strongest statistical comparison to those previously measured with the Hubble Space Telescope in the NIR. They also span the distance range of those used to determine the Hubble constant with HST, allowing us to search for a distance-dependent bias in HST measurements. The superior resolution of JWST negates crowding noise, the largest source of variance in the NIR Cepheid Period-Luminosity relations (Leavitt laws) measured with HST. Together with the use of two-epochs to constrain Cepheid phases and three filters to remove reddening, we reduce the dispersion in the Cepheid PL relations by a factor of 2.5. We find no significant difference in the mean distance measurements determined from HST and JWST, with a formal difference of -0.01+/-0.03 mag. This result is independent of zeropoints and analysis variants including metallicity dependence, local crowding, choice of filters, and relation slope. We can reject the hypothesis of unrecognized crowding of Cepheid photometry from HST that grows with distance as the cause of the ``Hubble Tension'' at 8.2 sigma, i.e., greater confidence than that of the Hubble Tension itself. We conclude that errors in photometric measurements of Cepheids across the distance ladder do not significantly contribute to the Tension.

Understanding how galaxies quench their star formation is crucial for studies of galaxy evolution. Quenching is related to the cold gas decrease. In the first paper we showed that the dust removal timescale in early-type galaxies (ETGs) is about 2.5 Gyr. Here we present carbon monoxide (CO) and 21 cm hydrogen (H I) line observations of these galaxies and measure the timescale of removal of the cold interstellar medium (ISM). We find that all the cold ISM components (dust, molecular and atomic gas) decline at similar rates. This allows us to rule out a wide range of potential ISM removal mechanisms (including starburst-driven outflows, astration, a decline in the number of asymptotic giant branch stars), and artificial effects like stellar mass-age correlation, environmental influence, mergers, and selection bias, leaving ionization by evolved low-mass stars and ionization/outflows by supernovae Type Ia or active galactic nuclei as viable mechanisms. We also provide evidence for an internal origin of the detected ISM. Moreover, we find that the quenching of star formation in these galaxies cannot be explained by a reduction in gas amount alone, because the star formation rates (SFRs) decrease faster (on a timescale of about 1.8 Gyr) than the amount of cold gas. Furthermore, the star formation efficiency of the ETGs (SFE = SFR/MH2) is lower than that of star-forming galaxies, whereas their gas mass fractions (fH2 = MH2/M*) are normal. This may be explained by the stabilization of gas against fragmentation, for example due to morphological quenching, turbulence, or magnetic fields.

The tip of the red giant branch (TRGB) allows for the measurement of precise and accurate distances to nearby galaxies, based on the brightest ascent of low-mass red giant branch stars before they undergo the helium flash. With the advent of JWST, there is great promise to utilize the technique to measure galaxy distances out to at least 50 Mpc, significantly further than HST's reach of 20 Mpc. However, with any standard candle, it is first necessary to provide an absolute reference. Here we use Cycle 1 data to provide an absolute calibration in the F090W filter. F090W is most similar to the F814W filter commonly used for TRGB measurements with HST, which had been adopted by the community due to minimal dependence from the underlying metallicities and ages of stars. The imaging we use was taken in the outskirts of NGC 4258, which has a direct geometrical distance measurement from the Keplerian motion of its water megamaser. Utilizing several measurement techniques, we find $M_{TRGB}^{F090W}$ = -4.362 $\pm$ 0.033 (stat) $\pm$ 0.045 (sys) mag (Vega) for the metal-poor TRGB. We also perform measurements of the TRGB in two Type Ia supernova hosts, NGC 1559, and NGC 5584. We find good agreement between our TRGB distances and previous distance determinations to these galaxies from Cepheids ($\Delta$ = 0.01 $\pm$ 0.06 mag), with these differences being too small to explain the Hubble tension ($\sim$0.17 mag). As a final bonus, we showcase the serendipitous discovery of a faint dwarf galaxy near NGC 5584.

We study stars in the J-regions of the asymptotic giant branch (JAGB) of near-infrared color magnitude diagrams in the maser host NGC 4258 and 4 hosts of 6 Type Ia supernovae (SN~Ia): NGC 1448, NGC 1559, NGC 5584, and NGC 5643. These clumps of stars are readily apparent near $1.0< F150W-F277W <1.5$ and $m_{F150W}$=22-25 mag with \textit{James Webb Space Telescope} NIRCam photometry. Various methods have been proposed to assign an apparent reference magnitude for this recently proposed standard candle, including the mode, median, sigma-clipped mean or a modeled parameter of its luminosity function. We test the consistency of these by measuring intra-host variations, finding differences of up to $\sim$0.2 mag that significantly exceed statistical uncertainties. Brightness differences appear intrinsic, and are further amplified by the non-uniform shape of the JAGB LF, an issue already apparent in comparing the LMC and SMC. We follow a ``many methods'' approach to consistently measure JAGB magnitudes and distances to the SN Ia host sample calibrated by NGC 4258. We find broad agreement with the distances measured from Cepheids, tip of the red giant branch (TRGB) and Miras. However, the SN host mean distance estimated via the JAGB method and necessary to estimate $H_0$, differs by $\sim$ 0.18 mag amongst the above definitions, a result of the different levels of asymmetry of the JAGB luminosity functions. The methods yield a full range of $72-78$ km s$^{-1}$ Mpc$^{-1}$, i.e., a fiducial result of $H_0=74.7 \pm 2.2$ (stat) $\pm$ 2.3 (sys) ($\pm$ 3.2 if combined in quadrature) km s$^{-1}$ Mpc$^{-1}$, independent of Cepheids, TRGB, or Miras, with systematic errors limited by the differences in methods. Future work may seek to further standardize and refine this promising tool, making it more competitive with established distance indicators.

Modified Newtonian Dynamics (MOND), postulating a breakdown of Newtonian mechanics at low accelerations, has considerable success at explaining galaxy kinematics. However, the quadrupole of the gravitational field of the Solar System (SS) provides a strong constraint on the way in which Newtonian gravity can be modified. In this paper we assess the extent to which modified gravity formulations of MOND are capable of accounting simultaneously for the Radial Acceleration Relation (RAR) -- encapsulating late-type galaxy dynamics -- the Cassini measurement of the SS quadrupole and the kinematics of wide binaries in the Solar neighbourhood. We achieve this by extending the method of Desmond (2023) to infer the location and sharpness of the MOND transition from the SPARC RAR under broad assumptions for the behaviour of the interpolating function and external field effect. We constrain the same quantities from the SS quadrupole, finding that it requires a significantly sharper transition between the deep-MOND and Newtonian regimes than is allowed by the RAR (an 8.7$\sigma$ tension under fiducial model assumptions). This may be relieved by allowing additional freedom in galaxies' mass-to-light ratios -- which also provides a better RAR fit -- and more significantly by removing galaxies with bulges. We show that the SS quadrupole constraint implies, to high precision, no deviation from Newtonian gravity in wide binaries in the Solar neighbourhood, and speculate on possible resolutions of this incompatibility between SS and galaxy data within the MOND paradigm.

As the closest major galaxy merger and home to thousands of super star clusters (SSCs), the Antennae Galaxies (NGC 4038 and NGC 4039) are an important location to study the molecular clouds at sites of vigorous star formation. We cataloged giant molecular clouds (GMCs) in the region where the two galaxies overlap using high-resolution (~0.1" ~10 pc) Atacama Large Millimeter/submillimeter Array (ALMA) observations of the 12CO(2-1) and 13CO(2-1)} emission lines. Of the 72 individual GMCs identified in the overlap region, 17 are within uncertainties of having the necessary mass, pressure, and size needed to form super star clusters (SSCs). Of those 17 GMCs, only one has significant ionizing radiation, indicating that the birth environments are likely still intact in the 16 other GMCs. We compared the physical properties calculated from 12CO(2-1) GMC data with observations of 10 other galaxies obtained using the same emission line and similar resolution. Compared to other sources in this sample, the GMCs from the Antennae, as well as in other starbursts and in the centers of galaxies, have the highest luminosities, surface densities, and turbulent pressures. The GMCs in starbursts and at the centers of galaxies also have large linewidths, although the linewidths in the Antennae are among the widest. These comparative results also indicate that the Antennae GMCs have the highest virial parameters despite their high densities.

We present a data reduction pipeline written in Python for data obtained with the near-infrared cross-dispersed echelle spectrograph, WINERED, which yields a 0.91$-$1.35 $\mu$m spectrum with the resolving power of $R_{\text{max}} \equiv \lambda / \Delta \lambda = 28,000$ or 70,000 depending on the observing mode. The pipeline was developed to efficiently extract the spectrum from the raw data with high quality. It comprises two modes: the calibration and the science mode. The calibration mode automatically produces the flat-fielding image, bad pixel map, echellogram distortion map and the dispersion solution from the set of the calibration data. Using calibration images and parameters, the science data of astronomical objects can be reduced automatically using the science mode. The science mode is also used for the real-time quick look at the data during observations. An example of the spectra reduced with WARP is presented. The effect of the highly inclined slit image on the spectral resolution is discussed.

The age and mass of red giants are essential for understanding the structure and evolution of the Milky Way. Traditional isochrone methods for these estimations are inherently limited due to overlapping isochrones in the Hertzsprung-Russell diagram, while asteroseismology, though more precise, requires high-precision, long-term observations. In response to these challenges, we developed a novel framework, Spectral Transformer (SPT), to predict the age and mass of red giants aligned with asteroseismology from their spectra. A key component of SPT, the Multi-head Hadamard Self-Attention mechanism, designed specifically for spectra, can capture complex relationships across different wavelength. Further, we introduced a Mahalanobis distance-based loss function to address scale imbalance and interaction mode loss, and incorporated Monte Carlo dropout for quantitative analysis of prediction uncertainty.Trained and tested on 3,880 red giant spectra from LAMOST, the SPT achieved remarkable age and mass estimations with average percentage errors of 17.64% and 6.61%, respectively, and provided uncertainties for each corresponding prediction. The results significantly outperform those of traditional machine learning algorithms and demonstrate a high level of consistency with asteroseismology methods and isochrone fitting techniques. In the future, our work will leverage datasets from the Chinese Space Station Telescope and the Large Synoptic Survey Telescope to enhance the precision of the model and broaden its applicability in the field of astronomy and astrophysics.

Using James Webb Space Telescope near-infrared data of the inner Orion Nebula, \citet{Pearson_McCaughrean_2023} detected 40 Jupiter-Mass Binary Objects (JuMBOS). These systems are not associated with stars and their components have masses of giant Jupiter-like planets and separations in the plane of the sky of order $\sim$100 au. The existence of these wide free-floating planetary mass binaries was unexpected in our current theories of star and planet formation. Here we report the radio continuum (6.1 and 10.0 GHz) Karl G.\ Jansky Very Large Array detection of a counterpart to JuMBO\,24. The radio emission appears to be steady at a level of $\sim$50 $\mu$Jy over timescales of days and years. We set an upper limit of $\simeq15$~km~s$^{-1}$ to the velocity of the radio source in the plane of the sky. As in the near-infrared, the radio emission seems to be coming from both components of the binary.

Massive star-forming galaxies in the high-redshift universe host large reservoirs of cold gas in their circumgalactic medium (CGM). Traditionally, these reservoirs have been linked to diffuse H I Lyman-$\alpha$ (Ly$\alpha)$ emission extending beyond $\approx 10$ kpc scales. In recent years, millimeter/submillimeter observations are starting to identify even colder gas in the CGM through molecular and/or atomic tracers such as the [C II] $158\,\mu$m transition. In this context, we study the well-known J1000+0234 system at $z=4.54$ that hosts a massive dusty star-forming galaxy (DSFG), a UV-bright companion, and a Ly$\alpha$ blob. We combine new ALMA [C II] line observations taken by the CRISTAL survey with data from previous programs targeting the J1000+0234 system, and achieve a deep view into a DSFG and its rich environment at a 0.2" resolution. We identify an elongated [C II]-emitting structure with a projected size of 15 kpc stemming from the bright DSFG at the center of the field, with no clear counterpart at any other wavelength. The plume is oriented $\approx 40^{\circ}$ away from the minor axis of the DSFG, and shows significant spatial variation of its spectral parameters. In particular, the [C II] emission shifts from 180 km/s to 400 km/s between the bottom and top of the plume, relative to the DSFG's systemic velocity. At the same time, the line width starts at 400-600 km/s but narrows down to 190 km/s at top end of the plume. We discuss four possible scenarios to interpret the [C II] plume: a conical outflow, a cold accretion stream, ram pressure stripping, and gravitational interactions. While we cannot strongly rule out any of these with the available data, we disfavor the ram pressure stripping scenario due to the requirement of special hydrodynamic conditions.

Recent deep and wide radio surveys extend the studies for radio-excess active galactic nuclei (radio-AGNs) to lower luminosities and higher redshifts, providing new insights into the abundance and physical origin of radio-AGNs. Here we focus on the cosmic evolution, physical properties and AGN-host galaxy connections of radio-AGNs selected from a sample of ~ 500,000 galaxies at 0 < z < 4 in GOODS-N, GOODS-S, and COSMOS fields. Combining deep radio data with multi-band, de-blended far-infrared (FIR) and sub-millimeter data, we identify 1162 radio-AGNs through radio excess relative to the FIR-radio relation. We study the cosmic evolution of 1.4 GHz radio luminosity functions (RLFs) for star-forming galaxies (SFGs) and radio-AGNs, which are well described by a pure luminosity evolution of $L_*\propto (1+z)^{-0.31z+3.41}$ and a pure density evolution of $\Phi_*\propto (1+z)^{-0.80z+2.88}$, respectively. We derive the turnover luminosity above which the number density of radio-AGNs surpasses that of SFGs. This crossover luminosity increases as increasing redshift, from $10^{22.9}$ W Hz$^{-1}$ at z ~ 0 to $10^{25.2}$ W Hz$^{-1}$ at z ~ 4. At full redshift range (0 < z < 4), we further derive the probability ($p_{radio}$) of SFGs and quiescent galaxies (QGs) hosting a radio-AGN as a function of stellar mass ($M_*$), radio luminosity ($L_R$), and redshift (z), which yields $p_{radio}\propto (1+z)^{3.54}M_*^{1.02}L_R^{-0.90}$ for SFGs, and $p_{radio}\propto (1+z)^{2.38}M_*^{1.39}L_R^{-0.60}$ for QGs, respectively. It indicates that radio-AGNs in QGs prefer to reside in more massive galaxies with larger $L_R$ than those in SFGs, and radio-AGN fraction increases towards higher redshift in both SFGs and QGs with a more rapid increase in SFGs. Further, we find that the radio-AGN fraction depends on accretion states of BHs and redshift in SFGs, while in QGs it also depends on BH (or galaxy) mass.

Quasi-periodic oscillations (QPOs) have been widely observed in black hole X-ray binaries (BHBs), which often exhibit significant X-ray variations. Extensive research has explored the long-term evolution of the properties of QPOs in BHBs. In contrast, such evolution in active galactic nuclei (AGNs) has remained largely unexplored due to limited observational data. By using the 10 new XMM-Newton observations for the narrow-line Seyfert 1 galaxy RE J1034+396 from publicly available data, we analyze the characteristics of its X-ray QPOs and examine their long-term evolution. The hard-band (1--4 keV) QPOs are found in all 10 observations and the frequency of these QPOs evolves ranging at $(2.47\text{--}2.83)\times10^{-4}\rm\ Hz$. Furthermore, QPO signals in the soft (0.3--1 keV) and hard bands exhibit strong coherence, although, at times, the variations in the soft band lead those in the hard band (the hard-lag mode), while at other times, it is the reverse (the soft-lag mode). The observations presented here serendipitously captured two ongoing lag reversals within about two weeks, which are first seen in RE J1034+396 and also among all AGNs. A transition in QPO frequency also takes place within a two-week timeframe, two weeks prior to its corresponding lag reversal, indicating a possible coherence between the transitions of QPO frequency and lag mode with delay. The diagram of time lag versus QPO frequency clearly evidences this interconnected evolution with hysteresis, which is, for the first time, observed among AGNs.

The JWST North Ecliptic Pole (NEP) Time Domain Field (TDF) is a $>$14 arcmin diameter field optimized for multi-wavelength time-domain science with JWST. Its location within JWST's northern continuous viewing zone offers optimal conditions for time-domain studies with JWST, ensuring year-round observability, absence of bright stars, and minimal Zodiacal foreground interference. It has been observed across the electromagnetic spectrum both from the ground and from space, including with the Hubble Space Telescope (HST). As part of HST observations over 3 cycles (the "TREASUREHUNT" program), deep images were obtained with ACS/WFC in F435W and F606W that cover almost the entire JWST NEP TDF. Many of the individual pointings of these programs partially overlap, allowing an initial assessment of the potential of this field for time-domain science with HST and JWST. The cumulative area of overlapping pointings is ~88 arcmin$^2$, with time intervals between individual epochs that range between 1 day and 4$+$ years. We present the discovery of 12 transients and ~100 variable sources to a depth of $m_{\rm AB} \simeq$ 29.5 mag (F606W). The majority of the transients will be supernovae, although at least two are likely quasars. Most variable sources are AGN, where we find ~0.35% (0.04%) of the general $z < 6$ field galaxy population to vary at the 3$\sigma$ (5$\sigma$) level. Based on a 5-year timeframe, this translates into a random supernova areal density of up to ~0.07 transients per arcmin$^2$ (~245 deg$^{-2}$) per epoch, and a variable AGN areal density of ~1.4 variables per arcmin$^2$ (~5184 deg$^{-2}$) to these depths.

We present results from the analysis of 88 carbon stars selected from Hamburg/ESO (HES) survey using low-resolution spectra (R$\sim$1330 \& 2190). The spectra were obtained with the Himalayan Faint Object Spectrograph Camera (HFOSC) attached to the 2-m Himalayan Chandra Telescope (HCT). Using a well-defined spectral criteria based on the strength of carbon molecular bands, the stars are classified into different groups. In our sample, we have identified 53 CH stars, four C-R stars, and two C-N type stars. Twenty-nine stars could not be classified due to the absence of prominent C$_{2}$ molecular bands in their spectra. We could derive the atmospheric parameters for 36 stars. The surface temperature is determined using photometric calibrations and synthesis of the H-alpha line profile. The surface gravity log g estimates are obtained using parallax estimates from the Gaia DR3 database whenever possible. Microturbulent velocity ($\zeta$) is derived using calibration equation of log g \& ${\zeta}$. We could determine metallicity for 48 objects from near-infrared Ca II triplet features using calibration equations. The derived metallicity ranges from $-$0.43$\leq$[Fe/H]$\leq$$-$3.49. Nineteen objects are found to be metal-poor ([Fe/H] $\leq$$-$1), 14 very metal-poor ([Fe/H]$\leq$$-$2), and five extremely metal-poor ([Fe/H]$\leq$$-$3.0) stars. Eleven objects are found to have a metallicity in the range $-$0.43 $\leq$[Fe/H]$\leq$$-$0.97. We could derive the carbon abundance for 25 objects using the spectrum synthesis calculation of the C$_{2}$ band around 5165\AA. The most metal-poor objects found will make important targets for follow-up detailed chemical composition studies based on high-resolution spectroscopy, that are likely to provide insight into the Galactic chemical evolution.

The satellite galaxy Crater II of the Milky Way is extremely cold and exceptionally diffuse. These unusual properties are challenging to understand in the standard model of cold dark matter. We investigate the formation of Crater II in self-interacting dark matter (SIDM), where dark matter particles can scatter and thermalize. We conduct a series of controlled N-body simulations to model the tidal evolution of Crater II, varying the self-interacting cross section, orbit parameters, and initial stellar distribution. Dark matter self-interactions lead to halo core formation and the distribution of stars expands accordingly. A cored SIDM halo also boosts tidal mass loss, allowing for a high orbit. Our simulations show that SIDM halos with a $1~{\rm kpc}$ core can simultaneously explain the low stellar velocity dispersion and the large half-light radius of Crater II, remaining robust to the initial distribution of stars. For the orbit motivated by the measurements from Gaia Early Data Release 3, the favored self-interacting cross section is approximately $60~{\rm cm^2/g}$ on the mass scale of Crater II.

We present an analysis of the outer Galaxy giant molecular filament (GMF) G214.5-1.8 (G214.5) using IRAM 30m data of $^{12}$CO, $^{13}$CO and C$^{18}$O. We find that the $^{12}$CO (1-0) and (2-1) derived excitation temperatures are near identical and are very low, with a median of 8.2 K, showing that the gas is extremely cold across the whole cloud. Investigating the abundance of $^{13}$CO across G214.5, we find that there is a significantly lower abundance along the entire 13 pc spine of the filament, extending out to a radius of $\sim 0.8$ pc, corresponding to $A_v \gtrsim 2$ mag and $T_{dust} \lesssim 13.5$ K. Due to this, we attribute the decrease in abundance to CO freeze-out, making G214.5 the largest scale example of freeze-out yet. We construct an axisymmetric model of G214.5's $^{13}$CO volume density considering freeze-out and find that to reproduce the observed profile significant depletion is required beginning at low volume densities, $n\gtrsim2000$ cm$^{-3}$. Freeze-out at this low number density is possible only if the cosmic-ray ionisation rate is $\sim 1.9 \times 10^{-18}$ s$^{-1}$, an order of magnitude below the typical value. Using timescale arguments, we posit that such a low ionisation rate may lead to ambipolar diffusion being an important physical process along G214.5's entire spine. We suggest that if low cosmic-ray ionisation rates are more common in the outer Galaxy, and other quiescent regions, cloud-scale CO freeze-out occurring at low column and number densities may also be more prevalent, having consequences for CO observations and their interpretation.

The prompt emission, X-ray plateau, and X-ray flares of GRB are thought to be from internal dissipation, and the magnetar as the central engine with propeller-fallback-accretion is proposed to interpret the observed phenomena of GRB. In this paper, by systematically searching for X-ray emission observed by Swift/XRT, we find that seven robust GRBs include both X-ray flares and plateau emissions with measured redshift. More interestingly, the X-ray flare/bump for those seven GRBs are simultaneously observed in $\gamma-$ray band. By adopting the propeller-fallback-accretion model to fit the observed data, it is found that the free parameters of two GRBs (140512A and 180329B) can be constrained very well, while the other five cases, more or less, are not all sufficiently constrained. On the other hand, it requires that the conversion efficiency of propeller is two or three times higher than that of spin-down dipole radiation of magnetar. If this is the case, it is contradictory to the expectation from the propeller model, namely, a dirtier ejecta should be less efficient in producing $\gamma-$ray emissions. Our results hint that at least the magnetar central engine with propeller-fallback-accretion model cannot interpret very well both the GRB X-ray flares simultaneously observed in $\gamma-$ray band and the X-ray flares of GRB with a high Lorentz factor.

Supernova remnants interacting with molecular/atomic clouds are interesting X-ray sources to study broadband nonthermal emission. X-ray line emission in these systems can be produced by different processes, e.g. low energy cosmic rays interacting with the cloud and fast ejecta fragments moving in the cloud. The paper aims at studying the origin of the non-thermal X-ray emission of the southwestern limb of SN 1006 beyond the main shock, in order to distinguish if the emission is due to low energy cosmic rays diffusing in the cloud or to ejecta knots moving into the cloud. We analyzed the X-ray emission of the southwestern limb of SN 1006, where the remnant interacts with an atomic cloud, with three different X-ray telescopes ({NuSTAR, Chandra and XMM-Newton) and performed a combined spectro-imaging analysis of this region. The analysis of the non thermal X-ray emission of the southwestern limb of SN 1006, interacting with an atomic cloud, has shown the detection of an extended X-ray source in the atomic cloud, approximately $2$ pc upstream of the shock front. The source is characterized by a hard continuum (described by a power law with photon index $\Gamma\sim1.4$) and by Ne, Si and Fe emission lines. The observed flux suggests that the origin of the X-ray emission is not associated with low energy cosmic rays interacting with the cloud. On the other hand, the spectral properties of the source, together with the detection of an IR counterpart visible with \textit{Spitzer}-MIPS at 24 $\mu$m are in good agreement with expectations for a fast ejecta fragment moving within the atomic cloud. We detected X-ray and IR emission from a possible ejecta fragment, with radius approximately 1$\times10^{17}$ cm, and mass approximately $10^{-3}M_\odot$ at about 2 pc out of the shell of SN 1006, in the interaction region between the southwestern limb of the remnant and the atomic cloud.

The jet composition of GRB plays an important role in understanding the energy dissipation and radiation mechanisms in GRB physics, but it is poorly constrained from the observational data. Recently, an interesting long-duration GRB 230307A with redshift $z=$0.065 has attracted great attention. The lack of detected thermal emission and mini-structure of prompt emission lightcurve of this burst suggest that the outflow is Poynting-flux-dominated and point towards the ICMART model. In this paper, we invoke two independent methods to investigate the jet composition of GRB 230307A. The high magnetization parameter ($\sigma>7$ or ever large) for$R_0=10^{10}$ cm that is used to suppress thermal component, strongly suggests that a significant fraction of the outflow energy is likely in a Poynting flux entrained with the baryonic matter. Moreover, it is found that the radiation efficiency of this burst for typical values $\epsilon_e=0.1$ and $\epsilon_B=0.01$ can reach as high as $~50\%$ which disfavors the internal shock model, but is consistent with ICMART model. Finally, a possible unified picture to produce GRB 230307A originated from a compact star merger is also discussed.

Estimating the magnetic field strength in the solar corona is crucial for understanding different physical processes happening over diverse spatio-temporal scales. However, the high temperatures and low density of the solar corona make this task challenging. The coronal magnetic field is too weak to produce a measurable splitting of the spectral lines using the Zeeman effect, and high temperature causes spectral lines to become weak and broad, making it difficult to detect the small Zeeman splitting. Coronal magneto-seismology, which combines the theoretical and observed properties of magnetohydrodynamic (MHD) waves, can be used to infer the magnetic field strength of oscillating structures in the solar corona, which are otherwise difficult to estimate. In this work, we use the Doppler velocity and density data obtained from the Coronal Multichannel Polarimeter (CoMP) on 2016 October 14 to obtain the global map of the coronal magnetic field using Bayesian inference. Two priors are used for plasma density, viz Gaussian and uniform distributions. Bayesian inference provides us with the probability distribution for the magnetic field strength at each location from 1.05 to 1.35 $R_\odot$. A comparison between the magnetic field obtained using simple inversion and Bayesian inference is also drawn. We find that the values obtained using simple inversion do not always match the maximum posterior estimates obtained using Bayesian inference. We find that the inferred values follow a power-law function for the radial variation of the coronal magnetic field, with the power-law indices for simple and Bayesian inversion being similar.

The spatial distribution of Galactic stars with different chemical abundances encodes information on the processes that drove the formation and evolution of the Milky Way. Survey selection functions are indispensable for analysing astronomical catalogues produced by large-scale surveys. The use of these selection functions in data modelling is more complex when data from different surveys are to be modelled simultaneously. We introduce a procedure for constructing the selection function of a sample of red clump stars that have parallaxes and elemental abundances from the Gaia mission. We separately constructed the selection function of the APOGEE DR17 red clump stars, which depends on very different observables and has a very different spatial coverage. We combined the two surveys and accounted for their joint selection function to provide strong constraints on the radial and vertical density distribution of mono-abundance populations, with Gaia offering a dense coverage of the solar neighbourhood, while APOGEE reaches larger distances near the Galactic plane. We confirm that the radial density profile steepens with increasing metallicity. The combined sample also indicates a metallicity-dependent flaring of the alpha-poor disc. We provide the code for constructing the Gaia selection function we used in this study through the GaiaUnlimited Python package.

Prevailing $N$-body planet formation models typically start with lunar-mass embryos and show a general trend of rapid migration of massive planetary cores to the inner Solar System in the absence of a migration trap. This setup cannot capture the evolution from a planetesimal to embryo, which is crucial to the final architecture of the system. We aim to model planet formation with planet migration starting with planetesimals of $\sim10^{-6}$ -- $10^{-4}M_\oplus$ and reproduce the giant planets of the Solar System. We simulated a population of 1,000 -- 5,000 planetesimals in a smooth protoplanetary disc, which was evolved under the effects of their mutual gravity, pebble accretion, gas accretion, and planet migration, employing the parallelized $N$-body code SyMBAp. We find that the dynamical interactions among growing planetesimals are vigorous and can halt pebble accretion for excited bodies. While a set of results without planet migration produces one to two gas giants and one to two ice giants beyond 6 au, massive planetary cores readily move to the inner Solar System once planet migration is in effect. Dynamical heating is important in a planetesimal disc and the reduced pebble encounter time should be considered in similar models. Planet migration remains a challenge to form cold giant planets in a smooth protoplanetary disc, which suggests an alternative mechanism is required to stop them at wide orbits.

We observe the magnetic field morphology towards a nearby star-forming filamentary cloud, G202.3+2.5, by the JCMT/POL-2 850 {\mu}m thermal dust polarization observation with an angular resolution of 14.4" (~0.053 pc). The average magnetic field orientation is found to be perpendicular to the filaments while showing different behaviors in the four subregions, suggesting various effects from filaments' collision in these subregions. With the kinematics obtained by N2H+ observation by IRAM, we estimate the plane-of-sky (POS) magnetic field strength by two methods, the classical Davis-Chandrasekhar-Fermi (DCF) method and the angular dispersion function (ADF) method, B_{pos,dcf} and B_{pos,adf} are ~90 {\mu}G and ~53 {\mu}G. We study the relative importance between the gravity (G), magnetic field (B) and turbulence (T) in the four subregions, find G > T > B, G >= T > B, G ~ T > B and T > G > B in the north tail, west trunk, south root and east wing, respectively. In addition, we investigate the projection effect on the DCF and ADF methods based on a similar simulation case and find the 3D magnetic field strength may be underestimated by a factor of ~3 if applying the widely-used statistical B_{pos}-to-B_{3D} factor when using DCF or ADF method, which may further underestimate/overestimate related parameters.

Planetesimals are believed to form by the gravitational collapse of aerodynamically concentrated clumps of pebbles. Many properties of the objects in the cold classical Kuiper belt -- such as binarity, rotation, and size distribution -- are in agreement with this gravitational collapse model. Further support comes from the pebble-pile structure inferred for comet nuclei. For this study, we simulated the final assembly of a planetesimal from the gravitational collapse of a rotating clump of pebbles. We implemented a numerical method from granular dynamics to follow the collapse that includes the transition from a pebble swarm to solid cells at a high density. We compared the shapes of the simulated planetesimals with the shapes of the lobes of contact binaries and bilobed Solar System objects. We find that the gravitational collapse of slowly rotating pebble clouds naturally explains the formation of flattened ellipsoidal bodies. This result agrees well with the flattened structure of the bilobed planetesimal Arrokoth and the shapes of the components of bilobed comets.

Electron-capture supernovae (EC-SNe) provide an alternative channel for producing neutron stars (NSs). They play an important role in the formation of double NS (DNS) systems and the chemical evolution of galaxies, and contribute to the NS mass distribution in observations. It is generally believed that EC-SNe originate from $e$-captures on $\rm^{24}Mg$ and $\rm^{20}Ne$ in the massive degenerate oxygen-neon (ONe) cores with masses close to the Chandrasekhar limit ($M_{\rm Ch}$). However, the origin of EC-SNe is still uncertain. In this paper, we systematically studied the EC-SNe in NS+He star systems by considering the explosive oxygen burning that may occur in the near-$M_{\rm Ch}$ ONe core. We provided the initial parameter spaces for producing EC-SNe in the initial orbital period $-$ initial He star mass (log$P_{\rm orb}^{\rm i}-M_{\rm He}^{\rm i}$) diagram, and found that both $M_{\rm He}^{\rm i}$ and minimum $P_{\rm orb}^{\rm i}$ for EC-SNe increase with metallicity. Then, by considering NS kicks added to the newborn NS, we investigated the properties of the formed DNS systems after the He star companions collapse into NSs, such as the orbital periods, eccentricities and spin periods of recycle pulsars ($P_{\rm spin}$), etc. The results show that most of the observed DNS systems can be produced by NS kicks of $\lesssim50\rm\,km\,s^{-1}$. In addition, we found that NSs could accrete more material if the residual H envelope on the He star companions is considered, which can form the mildly recycled pulsars ($P_{\rm spin}\sim20\,$ms) in DNS systems.

A recent study by Hon et al. reported that a close-in planet around the red clump star, 8 UMi, should have been engulfed during the expansion phase of the stellar evolution. They explained the survival of this exoplanet through a binary merger channel for 8 UMi. The key to testing this formation scenario is to derive the true age of this star: an old "imposter" resulting from a binary merger or a genuinely young red clump giant? To accomplish this, we derive kinematic and chemical properties for 8 UMi using astrometric data from $Gaia$ DR3 and the element-abundance pattern measured from a high-resolution ($R \sim 75,000$) spectrum taken by SOPHIE. Our analysis shows that 8 UMi is a normal thin-disk star with orbital rotation speed of $\it{V}_\mathrm{\phi}=\mathrm{244.96 km s^{-1}}$, and possesses a Solar metallicity ([Fe/H] $= -0.01 \pm 0.07$) and $\alpha$-element abundance ratio ([$\alpha$/Fe] $= 0.00 \pm 0.03$). By adopting well-established relationships between age and space velocities/elemental abundances, we estimate a kinematic age of $3.50^{+3.00}_{-2.00}$ Gyr, and a chemical age of $2.75^{+1.50}_{-1.00}$ Gyr from [C/N] and $3.47 \pm 1.96$ Gyr from [Y/Mg] for 8 UMi, respectively. These estimates are consistent with the isochrone-fitting age ($1.85^{+0.55}_{-0.30}$ Gyr) of 8 UMi, but are all much younger than 8.6 Gyr required in a binary merger scenario. This result challenges the binary merger model; the existence of such a closely orbiting exoplanet around a giant star remains a mystery yet to be resolved.

Individual subpulses of pulsars are regarded as the basic emission components, providing invaluable information to understand the radio emission process in the pulsar magnetosphere. Nevertheless, subpulses are overlapped with each other along the rotation phase for most pulsars, making it difficult to study the statistical properties of subpulses. Among the pulsars observed by the Five-hundred-meter Aperture Spherical radio Telescope, PSR B1916+14 has a large number of isolated well-resolved subpulses in the high time resolution observations, having a typical width of 0.15 ms and a high linear polarization. We find that the number distribution of subpulses contributes dominantly to the mean profile. According to the emission geometry, these emission units come from a region roughly 155 km above the polar cap in the pulsar magnetosphere, and the length scale of basic emission units is approximately 120 m. The deviations of polarization position angles for these single subpulses from the standard S-shaped curve are closely related to their fractional linear and circular polarization, and the large deviations tend to come from drifting subpulses.

Similarity solutions to the nonlinear non-equilibrium Marshak wave problem with a time dependent radiation driving source are presented. The radiation transfer model used is the gray, non-equilibrium diffusion approximation in the supersonic regime. These solutions constitute an extension of existing non-equilibrium supersonic Marshak wave solutions which are linear, to the nonlinear regime, which prevails in realistic high energy density systems. The generalized solutions assume a material model with power law temperature dependent opacities and a material energy density which is proportional to the radiation energy density, as well as a surface radiation temperature drive which obeys a temporal power-law. The solutions are analyzed in detail and it is shown that they take various qualitatively different forms according to the values of the opacity exponents. The solutions are used to construct a set of standardized benchmarks for supersonic non-equilibrium radiative heat transfer, which are nontrivial but straightforward to implement. These solutions are compared in detail to implicit Monte-Carlo and discrete-ordinate transport simulations as well gray diffusion simulations, showing a good agreement, which demonstrates the usefulness of these solutions as a code verification test problem.

A search for pulsed radiation at a frequency of 111 MHz in the direction of 116 RRAT candidates was carried out. For the search, archival data obtained on a meridian 128-beam radio telescope, a Large Phased Array (LPA), was used. For each candidate, about six days of observations were accumulated over an interval of eight years. Eleven new RRATs have been discovered. It was possible to estimate periods for six of them, and to construct average profiles for four of them. Some of the candidates turned out to be known pulsars observed in the side lobes of the radio telescope and interference. For the part of the candidates could not find pulses with a signal-to-noise ratio of more than seven, and their nature remains unknown.

Massive stars condition the evolution of the interstellar medium by the amount of energy released during their lives and especially by their deaths as supernova explosions. The vast amounts of spectroscopic data for massive stars provided by previous and existing instruments on ground-based and space-based telescopes have already saturated our capability to process them by the use of human routines. As a consequence, there is a pressing need for machine-assisted tools to help handle incoming data. To this end, we present the development of a massive star spectroscopic multiwavelength interactive database designed for scientific research and a fully automatic stellar parameter determination tool. Here we show the preliminary results of the application of these tools to optical spectra of O-type stars.

Methods are studied to compute the boosting effects produced by the observer motion that modifies and transfers to higher l the isotropic monopole frequency spectrum of the cosmic background. Explicit analytical solutions for spherical harmonic coefficients are presented and applied to various background spectra, alleviating computational effort. High l frequency spectra are led by higher order derivatives of the spectrum. Tabulated frequency spectra are computed with a relatively poor frequency resolution in comparison with the Doppler shift, calling for interpolation. They are affected by uncertainties due to intrinsic inaccuracies in modelling, observational data or limited computation accuracy, propagate and increase with the derivative order, possibly preventing a trustworthy computation to higher l and of the observed monopole. We filter the original function and the multipole spectra to derive reliable predictions of the harmonic coefficients. For spectra expressed in Taylor series, we derive explicit solutions for the harmonic coefficients up to l=6 in terms of spectra derivatives. We consider filters and study the quality of these methods on suitable analytical approximations, polluted with simulated noise. We consider the extragalactic sources microwave background from radio loud AGN and the 21cm line superimposed to the CMB. Gaussian pre-filtering coupled to a real space filtering of derivatives allows accurate predictions up to l=6, while log-log polynomial representation gives accurate solutions at any l. Describing the 21 cm model variety is difficult, so it is relevant to relax assumptions. Pre-filtering gives accurate predictions up to l=3-4, while further filtering or boosting amplification/deamplification method improves the results allowing reasonable estimations. The methods can extend the range of realistic background models manageable with a fast computation.

In this work, we revisit the relationship between [O III] line width $w_{\rm 90}$ (as the indicator of AGN outflow velocity) and the radio emission in RQQs by employing a large sample of Type I quasars ($\sim 37,000$) selected from the Sloan Digital Sky Survey (SDSS) Data Release Sixteen. By median stacking the radio images (to include the dominant fraction of individually radio non-detected RQQs) of Karl G. Jansky Very Large Array (VLA) Sky Survey (VLASS) for subsamples of RQQs with different $w_{\rm 90}$, our study demonstrates that, the correlation between $w_{\rm 90}$ and radio emission in our SDSS RQQs is significant, and remains solid after controlling the effects of black hole mass, quasar luminosity, Eddington ratio and redshift. This intrinsic link supports that the [O III] outflows in quasars, most likely resulted from wide-angled sub-relativistic quasar winds launched from the accretion disc, could make a dominant contribution to radio emission in the general RQQs. Alternatively, the correlation may be attributed to low-power jets in RQQs if they are ubiquitous and could efficiently enhance the [O III] width through interacting with the ISM. Meanwhile, the star-formation rates traced by the flux ratio of [Ne V]/[O II] emission lines display no dependence on $w_{\rm 90}$ after controlling the effects of black hole mass, quasar luminosity, Eddington ratio and redshift. This suggests that the stronger radio emission in RQQs with larger $w_{\rm 90}$ could not be attributed to outflow enhanced (positive feedback) star formation in the hosts. However, this also indicates the outflows, though exhibiting robust correlation with radio power, produce neither positive nor negative feedback to the star formation in their hosts.

So far no diffuse emissions in dwarf spheroidal satellites of the Milky Way have ever been observed. Given that dwarf galaxies are predominantly composed of Dark Matter, the discovery of these signals could offer valuable insights into understanding the nature of Dark Matter. We present "diffSph", a Python tool which in its present version provides fast predictions of such diffuse signals in radio frequencies. It also features a very comprehensive module for the computation of "J" and "D" factors that are relevant for indirect Dark Matter detection using gamma rays. Routines are coupled to parton-shower algorithms and Dark Matter halo mass functions from state-of-the-art kinematic fits. This code is also useful for testing generic hypotheses (not necessarily associated with any Dark Matter candidate) about the cosmic-ray electron/positron sources in the dwarf galaxies. The diffSph tool has already been employed in searches for diffuse signals from dwarf spheroidal galaxies using the LOw Frequency ARray (LOFAR).

We present the first JWST MIRI and NIRCam observations of the prominent debris disk around Beta Pictoris. Coronagraphic observations in 8 filters spanning from 1.8 to 23~$\mu$m provide an unprecedentedly clear view of the disk at these wavelengths. The objectives of the observing program were to investigate the dust composition and distribution, and to investigate the presence of planets in the system. In this paper, we focus on the disk components, providing surface brightness measurements for all images and a detailed investigation of the asymmetries observed. A companion paper by Kammerer et al. will focus on the planets in this system using the same data. We report for the first time the presence of an extended secondary disk in thermal emission, with a curved extension bent away from the plane of the disk. This feature, which we refer to as the ``cat's tail", seems to be connected with the previously reported CO clump, mid-infrared asymmetry detected in the southwest side, and the warp observed in scattered light. We present a model of this secondary disk sporadically producing dust that broadly reproduces the morphology, flux, and color of the cat's tail, as well as other features observed in the disk, and suggests the secondary disk is composed largely of porous, organic refractory dust grains.

Radial profiles play a crucial role in the analysis and interpretation of astronomical data, facilitating the extraction of spatial information. However, highly customizable (for different scenarios) measurements over each elliptical annulus can be challenging. In response, we present 'astscript-radial-profile', which is part of Gnuastro from version 0.15 and has an extensive documentation. A convenient feature of this program is its capability to make the measurements with different operators (mean, median, sigma-clipping, and many more) over ellipses, very quickly and directly on the command-line with minimal dependencies. This research note is reproducible with Maneage, on the Git commit 104aad5.

We present a detailed view of cluster formation (CF) to trace the evolution and interaction history of the Magellanic Clouds (MCs) in the last 3.5 Gyr. Using the \textit{Gaia} DR3 data, we parameterized 1710 and 280 star clusters in the Large Magellanic Cloud (LMC) and the Small Magellanic Cloud (SMC), where 847 and 113 clusters are newly characterized in the outer LMC and SMC, respectively. We estimated the age-extinction-metallicity-distance parameters using an automated fitting of the color-magnitude diagram (CMD) after field star removal, followed by an MCMC technique. We report a first-time detection of two synchronized CF peaks in the MCs at 1.5$\pm$0.12 Gyr and 800$\pm$60 Myr. We recommend that the choice of the metallicity ($Z$) values of isochrones for clusters with age $\le$ 1 - 2 Gyr are Z$_{\text{LMC}}$ = 0.004 - 0.008 and Z$_{\text{SMC}}$ = 0.0016 - 0.004 for the LMC and SMC, respectively. We found evidence for spiral arms in the LMC, as traced by the cluster count profiles over the last 3.5 Gyr. The density maps provide evidence of ram-pressure stripping in the North-East of the LMC, a severe truncation of CF in the South of the LMC, and a radial shrinkage of CF in the SMC in the last 450 Myr. The last SMC-LMC interaction ($\sim$ 150 Myr) resulted in a substantial CF in the North and Eastern SMC, with a marginal impact on the LMC. This study provides important insights into the CF episodes in the MCs and their connection to the LMC-SMC-MW interactions.

Multi-spacecraft observations of solar energetic particle (SEP) events not only enable a deeper understanding and development of particle acceleration and transport theories, but also provide important constraints for model validation efforts. However, because of computational limitations, a given physics-based SEP model is usually best-suited to capture a particular phase of an SEP event, rather than its whole development from onset through decay. For example, magnetohydrodynamic (MHD) models of the heliosphere often incorporate solar transients only at the outer boundary of their so-called coronal domain -- usually set at a heliocentric distance of 20-30 $R_{\odot}$. This means that particle acceleration at CME-driven shocks is also computed from this boundary onwards, leading to simulated SEP event onsets that can be many hours later than observed, since shock waves can form much lower in the solar corona. In this work, we aim to improve the modelled onset of SEP events by inserting a "fixed source" of particle injection at the outer boundary of the coronal domain of the coupled WSA-Enlil 3D MHD model of the heliosphere. The SEP model that we employ for this effort is SEPMOD, a physics-based test-particle code based on a field line tracer and adiabatic invariant conservation. We apply our initial tests and results of SEPMOD's fixed-source option to the 2021 October 9 SEP event, which was detected at five well-separated locations in the inner heliosphere -- Parker Solar Probe, STEREO-A, Solar Orbiter, BepiColombo, and near-Earth spacecraft.

The relative importance of magnetic fields, turbulence, and gravity in the early phases of star formation is still not well understood. We report the first high-resolution dust polarization observations at 850 $\mu$m around the most massive clump, located at the hub of the Giant Molecular Cloud G148.24+00.41, using SCUBA-2/POL-2 at the James Clerk Maxwell Telescope. We find that the degree of polarization decreases steadily towards the denser portion of the cloud. Comparing the intensity gradients and local gravity with the magnetic field orientations, we find that local gravity plays a dominant role in driving the gas collapse as the magnetic field orientations and gravity vectors seem to point towards the dense clumps. We also find evidence of U-shaped magnetic field morphology towards a small-scale elongated structure associated with the central clump, hinting at converging accretion flows towards the clump. Our observation has resolved the massive clump into multiple substructures. We study the magnetic field properties of two regions, central clump (CC) and northeastern elongated structure (NES). Using the modified Davis-Chandrasekhar Fermi method, we determine that the magnetic field strengths of CC and NES are $\sim$24.0 $\pm$ 6.0 $\mu$G and 20.0 $\pm$ 5.0 $\mu$G, respectively. The mass-to-flux ratios are found to be magnetically transcritical/supercritical, while the Alfv$\acute{\text{e}}$n Mach number indicates a trans-Alfv$\acute{\text{e}}$nic state in both regions. These results, along with Virial analysis, suggest that at the hub of G148.24+00.41, gravitational energy has an edge over magnetic and kinetic energies.

In this work, we compare star formation histories of massive (10.5 $< \log(\mathrm{M_*/M_{\odot}}) <$ 12) galaxies in the UniverseMachine model to those measured from the Large Early Galaxy Astrophysics Census (LEGA-C) at $0.6<z<1$. Following the LEGA-C study, we investigate how 50% ($t_{50}$) and 90% ($t_{90}$) formation timescales depend on total stellar mass. We find good agreement between the observed and model timescales for the star-forming population $\Delta\,t_{SF}\lesssim1\,\mathrm{Gyr}$ across the full mass range. In contrast, the observed age-mass correlation is weaker for the quiescent population compared to UniverseMachine models ($\Delta t_{Q}\lesssim2\,\mathrm{Gyr}$), especially at the high-mass end. This indicates continued star formation or additional processes in the most massive quiescent galaxies, a behavior not accounted for in the UniverseMachine model.

We present the second data release (DR2) of the Far-Infrared Polarimetric Large-Area CMZ Exploration (FIREPLACE) survey. This survey utilized the Stratospheric Observatory for Infrared Astronomy (SOFIA) High-resolution Airborne Wideband Camera plus (HAWC+) instrument at 214 $\mu$m (E-band) to observe dust polarization throughout the Central Molecular Zone (CMZ) of the Milky Way. DR2 consists of observations that were obtained in 2022 covering the region of the CMZ extending roughly from the Brick to the Sgr C molecular clouds (corresponding to a roughly 1$^{\circ}$ $\times$ 0.75$^{\circ}$ region of the sky). We combine DR2 with the first FIREPLACE data release covering the Sgr B2 region to obtain full coverage of the CMZ (a 1.5$^{\circ}$ $\times$0.75$^{\circ}$ region of the sky). After applying total and polarized intensity significance cuts on the full FIREPLACE data set we obtain $\rm\sim$65,000 Nyquist-sampled polarization pseudovectors. The distribution of polarization pseudovectors confirms a bimodal distribution in the CMZ magnetic field orientations, recovering field components that are oriented predominantly parallel or perpendicular to the Galactic plane. These magnetic field orientations indicate possible connections between the previously observed parallel and perpendicular distributions. We also inspect the magnetic fields toward a set of prominent CMZ molecular clouds (the Brick, Three Little Pigs, 50 km s$\rm^{-1}$, Circum-nuclear Disk, CO 0.02-0.02, 20 km s$\rm^{-1}$, and Sgr C), revealing spatially varying magnetic fields that generally trace the morphologies of the clouds. We find evidence that compression from stellar winds and shear from tidal forces are prominent mechanisms influencing the structure of the magnetic fields observed within the clouds.

The construction of high-resolution shock-capturing schemes is vital in producing highly accurate gravitational waveforms from neutron star binaries. The entropy based flux limiting (EFL) scheme is able to perform fast converging binary neutron star merger simulations reaching up to fourth-order convergence in the gravitational waveform phase. Here, we extend the applicability of the EFL method beyond special/general relativistic hydrodynamics to scalar conservation laws and show how to treat systems without a thermodynamic entropy. This is an indication that the method has universal applicability to any system of partial differential equations that can be written in conservation form. We also present some further very challenging special/general relativistic hydrodynamics applications of the EFL method and use it to construct eccentricity reduced initial data for a specific neutron star binary and show up to optimal fifth-order convergence in the gravitational waveform phase for this simulation.

In our galaxy, the white dwarfs (WDs) will inevitably capture the dark matter (DM) particles streaming through them, if there exist interactions between DM particles and nucleons/electrons. At the same time, these DM particles can also be evaporated by the nucleons/electrons in a WD if they have proper mass. The evaporation of DM particles will lead to a faster cooling evolution than that predicted by stellar evolution theory. In this letter, we ascribe the faster cooling evolution of 3 observed WDs to the capture and evaporation of DM particles, and get possible DM signals as follows: for $F(q) = 1$, $40\ \mathrm{MeV}/c^{2} \lesssim m_{\chi} \lesssim 70\ \mathrm{MeV}/c^{2}$ and $10^{-57} \mathrm{cm}^{2} \lesssim \sigma_{\chi,e} \lesssim 10^{-55} \mathrm{cm}^{2}$; for $F(q) = (\alpha m_{e})^{2}/q^{2}$, $30\ \mathrm{MeV}/c^{2} \lesssim m_{\chi} \lesssim 60\ \mathrm{MeV}/c^{2}$ and $10^{-53} \mathrm{cm}^{2} \lesssim \sigma_{\chi,e} \lesssim 10^{-51} \mathrm{cm}^{2}$. These results should be cross checked by more novel scenarios in the future.

The radiation damage to a silicon photomultiplier (SiPM) set on a satellite orbit increases energy threshold for scintillator detectors. We confirmed that 1 krad of radiation increases the energy threshold by approximately a factor of 10, which is worst for our system. Using one or two SiPMs damaged by proton irradiation and a plastic scintillator, we performed the following three experiments in our attempt to lower the energy threshold of radiation-damaged SiPMs to the greatest extent: (1) measurements using a current waveform amplifier rather than a charge-sensitive amplifier, (2) coincidence measurements with two radiation-damaged SiPMs attached to one scintillator and summing up their signals, and (3) measurements at a low temperature. Our findings confirmed that the use of a current waveform amplifier, as opposed to a charge-sensitive amplifier and a shaping amplifier, could lower the energy threshold to approximately 65% (from 198 keV to 128 keV). Furthermore, if we set the coincidence width appropriately and sum up the signals of the two SiPMs in the coincidence measurement, the energy threshold could be lowered to approximately 70% (from 132 keV to 93 keV) with little loss of the acquired signal, compared to that of use of only one scintillator. Finally, if we perform our measurements at a temperature of -20 {\deg}C, we could lower the energy threshold to approximately 34% (from 128 keV to 43 keV) compared to that of at 20 {\deg}C. Accordingly, we conclude that the energy threshold can be lowered to approximately 15% by using a combination of these three methods.

In this work, we analyze in detail the problem of piston driven shock waves in planar media. Similarity solutions to the compressible hydrodynamics equations are developed, for a strong shock wave, generated by a time dependent pressure piston, propagating in a non-homogeneous planar medium consisting of an ideal gas. Power law temporal and spatial dependency is assumed for the piston pressure and initial medium density, respectively. The similarity solutions are written in both Lagrangian and Eulerian coordinates. It is shown that the solutions take various qualitatively different forms according to the value of the pressure and density exponents. We show that there exist different families of solutions for which the shock propagates at constant speed, accelerates or slows down. Similarly, we show that there exist different types of solutions for which the density near the piston is either finite, vanishes or diverges. Finally, we perform a comprehensive comparison between the planar shock solutions and Lagrangian hydrodynamic simulations, by setting proper initial and boundary conditions. A very good agreement is reached, which demonstrates the usefulness of the analytic solutions as a code verification test problem.

Flavor-dependent neutrino transport is described by a well-known kinetic equation for occupation-number matrices in flavor space. However, as an overlooked theoretical problem, we show that in the inhomogeneous case, neutrino-neutrino refractive energy is not conserved. We derive the missing gradient terms in the fast flavor limit (vanishing neutrino masses), and prove that the missing refractive energy is traded with the huge reservoir of neutrino kinetic energy through gradients of the weak interaction potential. Even small changes of the kinetic energy accommodate the refractive energy gained or lost. Flavor evolution alone is negligibly affected by the new terms.

We obtain bounds on dark matter annihilation using 14 years of publicly available Fermi-LAT data from a set of 54 dwarf spheroidal galaxies, using spectral information from 16 energy bins. We perform this analysis using our updated and publicly available code MADHATv2, which can be used to test a variety of models for dark matter particle physics and astrophysics in an accessible manner. In particular, we note that including Carina III in the analysis strengthens constraints on $s$-wave annihilation into two-body Standard Model final states by a factor of $\sim 3$ but broadens the error on the constraint due to the large uncertainty of its $J$-factor. Our findings illustrate the importance of verifying if Carina III is in fact a dwarf spheroidal galaxy and measuring more precisely its $J$-factor. More generally, they highlight the significance of forthcoming discoveries of nearby ultra-faint dwarfs for dark matter indirect detection.