Papers for Wednesday, Oct 16 2024

Deep learning-based methods have been widely researched in the areas of language and vision, demonstrating their capacity to understand long sequences of data and their usefulness in numerous helio-physics applications. Foundation models (FMs), which are pre-trained on a large-scale datasets, form the basis for a variety of downstream tasks. These models, especially those based on transformers in vision and language, show exceptional potential for adapting to a wide range of downstream applications. In this paper, we provide our perspective on the criteria for designing an FM for heliophysics and associated challenges and applications using the Solar Dynamics Observatory (SDO) dataset. We believe that this is the first study to design an FM in the domain of heliophysics.

With a temperature of $\sim 285$ K WISE0855 is the coldest brown dwarf observed so far. Using the James Webb Space Telescope (JWST) we obtained observations that allow us to characterize WISE0855s atmosphere focusing on vertical variation in the water steam abundance, measuring trace gas abundances and receiving bulk parameters for this cold object. We observed the ultra cool dwarf WISE0855 using the Mid-Infrared Instrument Medium Resolution Spectrometer (MIRI/MRS) onboard JWST at a spectral resolution of up to 3750. We combined the observation with published data from the Near Infrared Spectrograph (NIRSpec) G395M and PRISM modes yielding a spectrum ranging from 0.8 to 22 um. We apply atmospheric retrievals using petitRADTRANS to measure atmospheric abundances, the pressure-temperature structure, radius and gravity of the brown dwarf. We also employ publicly available clear and cloudy self-consistent grid models to estimate bulk properties of the atmosphere such as the effective temperature, radius, gravity and metallicity. Atmospheric retrievals constrain a variable water abundance profile in the atmosphere, as predicted by equilibrium chemistry. We detect the 15NH3 isotopologue and infer a ratio of mass fraction of 14NH3/15NH3 = 332+63-43 for the clear retrieval. We measure the bolometric luminosity by integrating the presented spectrum and obtain a value of log(L/L$_{\odot}$) = -7.291+/-0.008. The detected water depletion indicates that water condenses out in the upper atmosphere due to the very low effective temperature of WISE0855. The height in the atmosphere where this occurs is covered by the MIRI/MRS data, and thus demonstrates the potential of MIRI to characterize cold gas giants atmospheres. Comparing the data to retrievals and self-consistent grid models, we do not detect signs for water ice clouds, although their spectral features have been predicted in previous studies.

We present the first seconds-long multidimensional general relativistic neutrino magnetohydrodynamic simulations of accretion-induced collapse (AIC) in rapidly rotating, magnetized white dwarfs (WDs). Our findings show that, with sufficiently strong initial magnetic fields and rotation, these systems can generate relativistic jets and neutron-rich outflows with properties consistent with long gamma-ray bursts (LGRBs) accompanied by kilonovae, such as GRB 211211A and GRB 230703A. Notably, our models reproduce the power and duration of these LGRBs and the associated kilonovae without the need for tuning parameters. Additionally, our results suggest that the AIC of WDs could be an important site of the production for heavy $r$-process elements.

Located 10.888 pc from Earth, COCONUTS-2b is a planetary-mass companion to a young (150--800 Myr) M3 star, with a wide orbital separation (6471 au) and a low companion-to-host mass ratio ($0.021\pm0.005$). We have studied the atmospheric properties of COCONUTS-2b using newly acquired 1.0--2.5~$\mu$m spectroscopy from Gemini/Flamingos-2. The spectral type of COCONUTS-2b is refined to T$9.5 \pm 0.5$ based on comparisons with T/Y dwarf spectral templates. We have conducted an extensive forward-modeling analysis, comparing the near-infrared spectrum and mid-infrared broadband photometry with sixteen state-of-the-art atmospheric model grids developed for brown dwarfs and self-luminous exoplanets near the T/Y transition. The \texttt{PH$_{3}$-free ATMO2020++}, \texttt{ATMO2020++}, and \texttt{Exo-REM} models best match the specific observations of COCONUTS-2b, regardless of variations in the input spectrophotometry. This analysis suggests the presence of disequilibrium chemistry, along with a diabatic thermal structure and/or clouds, in the atmosphere of COCONUTS-2b. All models predict fainter $Y$-band fluxes than observed, highlighting uncertainties in the alkali chemistry models and opacities. We determine a bolometric luminosity of $\log{(L_{\rm bol}/L_{\odot})}=-6.18$ dex, with a 0.5 dex-wide range of $[-6.43,-5.93]$ dex that accounts for various assumptions of models. Using thermal evolution models, we derive an effective temperature of $T_{\rm eff}=483^{+44}_{-53}$ K, a surface gravity of $\log{(g)}=4.19^{+0.18}_{-0.13}$ dex, a radius of $R=1.11^{+0.03}_{-0.04}$ R$_{\rm Jup}$, and a mass of $M=8 \pm 2$ M$_{\rm Jup}$. Various atmospheric model grids consistently indicate that COCONUTS-2b's atmosphere has sub- or near-solar metallicity and C/O. These findings provide valuable insights into COCONUTS-2b's formation history and the potential outward migration to its current wide orbit.

A fundamental ingredient in the unified model of active galactic nuclei (AGN) is the obscuring torus, whose innermost, hottest region dominates the near infrared (NIR) emission. Characterising the change in the torus properties and its interplay with the main AGN emission is key for our understanding of AGN physics, evolution and classification. Its covering factor ($CF$) is largely responsible for the classification of AGN on the basis of the detection of broad emission lines. It is still not clear whether the torus properties evolve over time and how they relate with the accretion parameters of the nucleus. In this work, we aim at investigating the evolution of the NIR properties with the redshift ($z$) and the bolometric luminosity ($L_{\rm bol}$) of the AGN. To this end, we assembled a large dataset of $\sim$36,000 Type 1 AGN between $0.5<z<2.9$ and $45.0<\log(L_{\rm bol} / (\rm erg / s))<48.0$ with UV, optical and near-infrared photometry. We produced average spectral energy distributions (SED) in different bins of the $z-L_{\rm bol}$ parameter space to estimate how the NIR SED evolves according to these parameters. We find that the NIR luminosity decreases for increasing $L_{\rm bol}$ at any redshift. At the same, time the shape of the NIR SED in our sample is consistent with a non-evolution with $z$. As a consequence, all the explored proxies for the $CF$ exhibit significant anti-correlations with $L_{\rm bol}$, but not with $z$. Additionally, the $CF$ also shows a shallower anti-correlation with the Eddington ratio ($\lambda_{\rm Edd}$), yet current systematic uncertainties, as well as the limited dynamical range, do not allow us to precisely constrain the role of the Eddington ratio. Lastly, we derived the covering factor from the ratio between the NIR and optical luminosity and we employed it to set a lower limit for the X-ray obscuration at different redshifts.

When a hot Jupiter orbits a star whose effective temperature exceeds $\sim$6100 K, its orbit normal tends to be misaligned with the stellar spin axis. Cooler stars have smaller obliquities. The latter may have been damped by hot Jupiters in resonance lock with axisymmetric stellar gravity modes (azimuthal number $m=0$), as has been recently recognized. Here we allow for resonance locks with non-axisymmetric modes, which affect both stellar obliquity and spin frequency. Obliquities damp for all modes ($-2 \leq m \leq 2$). Stars spin up for $m > 0$, and spin down for $m < 0$. We carry out a population synthesis that assumes hot Jupiters form misaligned around both cool and hot stars, and subsequently lock onto modes whose $m$-values yield the highest mode energies for given starting obliquities. Core hydrogen burning enables hot Jupiters to torque low-mass stars, but not high-mass stars, into spin-orbit alignment. Resonance locking plus stellar spin-down from magnetic braking largely reproduces observed obliquities and stellar rotation rates and how they trend with stellar effective temperature and orbital separation.

Supermassive binary black holes in galactic centers are potential multimessenger sources in gravitational waves and electromagnetic radiation. To find such objects, isolating unique electromagnetic signatures of their accretion flow is key. With the aid of three-dimensional general-relativistic magnetohydrodynamic (GRMHD) simulations that utilize an approximate, semi-analytic, super-imposed spacetime metric, we identify two such signatures for merging binaries. Both involve magnetic reconnection and are analogous to plasma processes observed in the solar corona. The first, like colliding flux tubes that can cause solar flares, involves colliding jets that form an extended reconnection layer, dissipating magnetic energy and causing the two jets to merge. The second, akin to coronal mass ejection events, involves the accretion of magnetic field lines onto both black holes; these magnetic fields then twist, inflate, and form a trailing current sheet, ultimately reconnecting and driving a hot outflow. We provide estimates for the associated electromagnetic emission for both processes, showing that they likely accelerate electrons to high energies and are promising candidates for continuous, stochastic, and/or quasi-periodic higher energy electromagnetic emission. We also show that the accretion flows around each black hole can display features associated with the magnetically arrested state. However, simulations with black hole spins misaligned with the orbital plane and simulations with larger Bondi radii saturate at lower values of horizon-penetrating magnetic flux than standard magnetically arrested disks, leading to weaker, intermittent jets due to feedback from the weak jets or equatorial flux tubes ejected by reconnecting field lines near the horizon.

We performed a chemical analysis of the asteroid-bearing white dwarf WD 1145+017 using optical and ultraviolet spectroscopic data from 25 epochs between 2015 and 2023. We present an updated gas disk model with improved opacity calculations and temperature profiles to properly account for all circumstellar absorption features. Incorporating these changes into our models, we identified at least 10 elements in the disk, including a detection of circumstellar Na. We detected 16 elements in the photosphere, including new detections of P, Co and Cu. At 16 elements, WD 1145+017 ties GD 362 as one of the most polluted white dwarfs in terms of the number of elements detected. We find that both the disk and photosphere compositions align, to first order, with CI Chondrite. Our study underscores the importance of accounting for circumstellar absorption, as neglecting them leads to significant abundance errors. Additionally, the analysis of the disk's opacity highlighted a ultraviolet flux reduction due to a pseudo-continuum due to a optically thick component. This result may affect previous analyses of other polluted white dwarfs, suggesting a need for revisiting some studies.

We present the first spectroscopic confirmation of a dust-obscured dwarf galaxy, CEERS-14821. The analysis is performed combining JWST NIRCam broad-band photometry and NIRSpec/PRISM spectroscopic data. From the detection of multiple rest-frame optical lines, we derive that CEERS-14821 is located at $z=4.883\pm0.003$. Moreover, from a secure detection of the $H_{\alpha}$ and $H_{\beta}$ we derived that the galaxy has a dust extinction ranging from Av=2.2 to Av=3.3, depending on the assumed reddening law. This value is extremely large given that we estimated a low stellar mass around log(M/Mo)=8.0-8.2. Moreover, using different metallicity tracers, we verify that the galaxy is also metal-rich, with 12+log(O/H)>8.3. This is well above the expectation from both the mass-metallicity relation and the fundamental mass-metalliticy relation. CEERS-14821 is going through a burst of star formation, there are no indications of a strong contribution from an active galactic nuclei (f(AGN)<0.5 with respect to the total dust luminosity). Based on the rest-frame optical images, this source has a size compatible with galaxies of similar stellar mass and redshift. Finally, with the current data, it seems that there are galaxies closely interacting with CEERS-14821.

Recent high-cadence transient surveys have discovered rapid transients whose light curve timescales are shorter than those of typical supernovae. In this paper, we present a systematic search for rapid transients at medium-high redshifts among 3381 supernova candidates obtained from the Subaru HSC-SSP transient survey. We developed a machine learning classifier to classify the supernova candidates into four types (Type Ia, Ibc, II supernovae, and rapid transients) based on the features derived from the light curves. By applying this classifier to the 3381 supernova candidates and by further applying the quality cut, we selected 14 rapid transient samples. They are located at a wide range of redshifts ($0.34 \leq z \leq 1.85$) and show a wide range of the peak absolute magnitude ($-17 \geq M \geq -22$). The event rate of the rapid transients is estimated to be $\sim 6\times10^3~\rm{events~yr^{-1}~Gpc^{-3}}$ at $z \sim 0.74$, which corresponds to about $2$ $\%$ of the event rate of normal core-collapse supernovae at the similar redshift. Based on the luminosity and color evolution, we selected two candidates of Type Ibn supernovae at $z\sim0.75$. The event rate of Type Ibn SN candidates is more than 1 $\%$ of Type Ib SN rate at the same redshift, suggesting that this fraction of massive stars at this redshift range eruptively ejects their He-rich envelope just before the explosions. Also, two objects at $z=1.37$ and 1.85 show high luminosities comparable to superluminous supernovae. Their event rate is about 10-25 $\%$ of superluminous supernovae at $z\sim 2$.

We performed the longest numerical-relativity neutrino-radiation magnetohydrodynamics simulation for a binary neutron star merger that extends to $\approx1.5\mathrm{\,s}$ after the merger. We consider the binary model that undergoes the prompt collapse to a black hole after the merger with asymmetric mass 1.25$\,M_{\odot}$ and 1.65$\,M_{\odot}$ and SFHo equation of state. We find the Poynting flux-driven collimated outflow as well as the gravitational wave emission, neutrino emission, dynamical mass ejection, and post-merger mass ejection facilitated by magnetorotational instability-driven turbulent viscosity in a single self-consistent binary neutron star merger simulation. A magnetosphere dominated by the aligned global magnetic field penetrating the black hole develops along the black-hole spin axis after the turbulence in the remnant disk is enhanced. A jet with the Poynting flux with isotropic-equivalent luminosity of $\sim10^{49}\mathrm{\,erg/s}$ is launched, and the duration of the high luminosity is expected to be $O(1)\mathrm{\,s}$.

Accurate classification of transients obtained from spectroscopic data are important to understand their nature and discover new classes of astronomical objects. For supernovae (SNe), SNID, NGSF (a Python version of SuperFit), and DASH are widely used in the community. Each tool provides its own metric to help determine classification, such as rlap of SNID, chi2/dof of NGSF, and Probability of DASH. However, we do not know how accurate these tools are, and they have not been tested with a large homogeneous dataset. Thus, in this work, we study the accuracy of these spectral classification tools using 4,646 SEDMachine spectra, which have accurate classifications obtained from the Zwicky Transient Facility Bright Transient Survey (BTS). Comparing our classifications with those from BTS, we have tested the classification accuracy in various ways. We find that NGSF has the best performance (overall Accuracy 87.6% when samples are split into SNe Ia and Non-Ia types), while SNID and DASH have similar performance with overall Accuracy of 79.3% and 76.2%, respectively. Specifically for SNe Ia, SNID can accurately classify them when rlap > 15 without contamination from other types, such as Ibc, II, SLSN, and other objects that are not SNe (Purity > 98%). For other types, determining their classification is often uncertain. We conclude that it is difficult to obtain an accurate classification from these tools alone. This results in additional human visual inspection effort being required in order to confirm the classification. To reduce this human visual inspection and to support the classification process for future large-scale surveys, this work provides supporting information, such as the accuracy of each tool as a function of its metric.

Observational constraints on the configuration of the black hole (BH)-accretion disk-jet system are crucial for addressing key questions in black hole growth, accretion disk physics, and jet formation. The recently reported jet precession in M87 provides a novel avenue to explore these long-standing issues. This precession, attributed to the accretion disk's response to the frame-dragging effect of a spinning supermassive black hole (SMBH), indicates a non-zero spin. The relatively short precession period ($\sim$11 years) implies a compact accretion disk. In contrast to the traditional view of a strictly collimated shape, the M87 jet is inferred to curve at the innermost regions connecting to the spinning BH, which explains the unexpectedly wide innermost projected jet width.

We present JWST/NIRSpec PRISM follow-up of candidate galaxies at z=9-11 selected from deep JWST/NIRCam photometry in GLASS-JWST Early Release Science data. We spectroscopically confirm six sources with secure redshifts at z = 9.52-10.43, each showing multiple emission lines. An additional object is likely at z = 10.66, based on its Lya-break and a single emission feature, while one source is a lower redshift interloper. The sample includes the first JWST-detected candidate at z=10, GHZ1/GLASS-z10, which we confirm at z = 9.875, and the X-ray detected AGN GHZ9 confirmed at z = 10.145. Three objects in our sample, including GHZ9, have EW(CIII])>20A and occupy a region compatible with AGN emission in the EW(CIII]) vs CIV/CIII] diagram. The spectroscopic sample confirms a high abundance of galaxies at z > 9. We measure a number density of z=10 galaxies in the GLASS-JWST ERS field that is a factor of >3 higher than other JWST-based estimates at demagnified rest-frame magnitudes of -21 < Muv < -19. We find that the positions of these galaxies in redshift and angular space are not consistent with all of them being part of a unique progenitor of present-day galaxy clusters. The high density of objects in the GLASS region can be explained either by clustering on large scales or by a superposition of different forming structures of which we observe only the brightest members. By considering all the spectroscopic z=10 sources in the Abell-2744 field, we identify two potential galaxy proto-clusters centered around GHZ9 and JD1, with relative separations between their members of 1-2 pMpc. The potential AGN nature of three of the sources in our sample lends support to a scenario in which the high abundance of bright sources determined by JWST surveys at cosmic dawn may be affected by AGN contribution to their UV luminosity.

The streaming instability (SI) is currently the leading candidate for triggering planetesimal formation in protoplanetary disks. Recently, a novel variation, the `azimuthal-drift' streaming instability (AdSI), was discovered in disks exhibiting laminar gas accretion. Unlike the classical SI, the AdSI does not require pressure gradients and can concentrate dust even at low abundances. We extend previous simulations of the AdSI to explore the impact of dust abundance, accretion flow strength, pressure gradients, and grain size. For a dimensionless accretion flow strength $\alpha_{\mathrm{M}}=0.1$ and particle Stokes number $\operatorname{St}=0.1$, we find the AdSI produces dust filaments for initial dust-to-gas ratios as low as $\epsilon=0.01$. For $\epsilon\gtrsim 1$, maximum dust-to-gas ratios of order 100 are attained, which can be expected to undergo gravitational collapse. Furthermore, even in systems dominated by the classical SI, an accretion flow drives filament formation, without which the disk remains in a state of small-scale turbulence. Our results suggest that an underlying accretion flow facilitates dust concentration and may thus promote planetesimal formation.

The continuing exploration of neighboring planetary systems is providing deeper insights into the relative prevalence of various system architectures, particularly with respect to the solar system. However, a full assessment of the dynamical feasibility of possible terrestrial planets within the Habitable Zones (HZ) of nearby stars requires detailed knowledge of the masses and orbital solutions of any known planets within these systems. Moreover, the presence of as yet undetected planets in or near the HZ will be crucial for providing a robust target list for future direct imaging surveys. In this work, we quantify the distribution of uncertainties on planetary masses and semi-major axes for 1062 confirmed planets, finding median uncertainties of 11.1% and 2.2%, respectively. We show the dependence of these uncertainties on stellar mass and orbital period, and discuss the effects of these uncertainties on dynamical analyses and the locations of mean motion resonance. We also calculate the expected radial velocity (RV) semi-amplitude for a Neptune-mass planet in the middle of the HZ for each of the proposed Habitable Worlds Observatory target stars. We find that for more than half of these stars, the RV semi-amplitude is less than 1.5 m/s, rendering them unlikely to be detected in archival RV data sets and highlighting the need for further observations to understand the dynamical viability of the HZ for these systems. We provide specific recommendations regarding stellar characterization and RV survey strategies that work toward the detection of presently unseen perturbers within the HZ.

We present a sample of 1956 individual stellar clumps at redshift 0.7<z<10, detected with JWST/NIRCam in 476 galaxies lensed by the galaxy cluster Abell2744. The lensed clumps present magnifications ranging between $\mu$=1.8 and $\mu$=300. We perform simultaneous size-photometry estimates in 20 JWST/NIRCam median and broad-band filters from 0.7 to 5 $\mu$m. Spectral energy distribution (SED) fitting analyses enable us to recover the physical properties of the clumps. The majority of the clumps are spatially resolved and have effective radii in the range R_eff = 10 - 700 pc. We restrict this first study to the 1751 post-reionisation era clumps with redshift <5.5. We find a significant evolution of the average clump ages, SFR, SFR surface densities and metallicity with increasing redshift, while median stellar mass and stellar mass surface densities are similar in the probed redshift range. We observe a strong correlation between the clump properties and the properties of their host galaxies, with more massive galaxies hosting more massive and older clumps. We find that clumps closer to their host galactic center are on average more massive, while their ages do not show clear sign of migration. We find that clumps at cosmic noon sample the upper-mass end of the mass function to higher masses than at z>3, reflecting the rapid increase towards the peak of the cosmic star formation history. We conclude that the results achieved over the studied redshift range are in agreement with expectation of in-situ clump formation scenario from large-scale disk fragmentation.

Dark energy away from a cosmological constant $\Lambda$ -- like early universe inflation that ends -- can be understood in terms of well defined physical behaviors. These guide dark energy into thawing or freezing classes, with $w_0$--$w_a$ arising as a physical calibration of the phase space. Other regions of phase space -- zones of avoidance -- require violation of some basic principle. We explore these cases, drawing a direct analogy with how nonGaussianity in inflation can add physics beyond standard dynamics. We examine the physics implications if the best fit of current data is taken to be truth, outlining four properties, and investigate the reality of phantom crossing $w=-1$, finding it significantly favored.

In the Milky Way and other main-sequence galaxies, stars form exclusively in molecular gas, which is traced by CO emission. However, low metallicity dwarf galaxies are often `CO-dark' in the sense that CO emission is not observable even at the high resolution and sensitivities of modern observing facilities. In this work we use ultra high-resolution simulations of four low-metalicity dwarf galaxies (which resolve star formation down to the scale of star-forming cores, 0.01 pc) combined with a time-dependent treatment of the chemistry of the interstellar medium, to investigate the star formation environment in this previously hidden regime. By generating synthetic observations of our models we show that the galaxies have high to extremely high dark gas fractions (0.13 to 1.00 dependent on beam size and conditions), yet despite this form stars. However, when examined on smaller scales, we find that the stars still form in regions dominated by molecular gas, it is simply that these are far smaller than the scale of the beam (1.5"). Thus, while stars in CO-dark dwarf galaxies form in small molecular cores like larger galaxies, their cloud-scale environment is very different.

The nearby transiting rocky exoplanet LTT 1445A b presents an ideal target for studying atmospheric retention in terrestrial planets orbiting M dwarfs. It is cooler than many rocky exoplanets yet tested for atmospheres, receiving a bolometric instellation similar to Mercury's. Previous transmission spectroscopy ruled out a light H/He-dominated atmosphere but could not distinguish between a bare-rock, a high-MMW, or a cloudy atmosphere. We present new secondary eclipse observations using JWST's MIRI/LRS, covering the 5-12 $\mu$m range. From these observations, we detect a broadband secondary eclipse depth of 41 $\pm$ 9 ppm and measure a mid-eclipse timing consistent with a circular orbit (at 1.7$\sigma$). From its emission spectrum, the planet's dayside brightness temperature is constrained to 525 $\pm$ 15 K, yielding a temperature ratio relative to the maximum average dayside temperature from instant thermal reradiation by a rocky surface $R$ = $T_{\rm day,obs}/T_{\rm max}$ = 0.952 $\pm$ 0.057, consistent with emission from a dark rocky surface. From an energy balance perspective, such a warm dayside temperature disfavors thick atmospheres, excluding $\sim$100 bar atmospheres with Bond albedo $>$ 0.08 at the 3$\sigma$ level. Furthermore, forward modeling of atmospheric emission spectra disfavor simple 100\% CO$_2$ atmospheres with surface pressures of 1, 10, and 100 bar at 4.2$\sigma$, 6.6$\sigma$, and 6.8$\sigma$ confidence, respectively. These results suggest that LTT 1445A b lacks a very thick CO$_2$ atmosphere, possibly due to atmospheric erosion driven by stellar activity. However, the presence of a moderately thin atmosphere (similar to those on Mars, Titan, or Earth) remains uncertain.

Accurately quantifying the ionizing photon production efficiency ($\xi_\rm{ion}$) of $z>6$ star-forming galaxies (SFGs) is necessary to understand their contribution to reionization. We investigate the ionizing properties of N=289 SFGs selected at $z=6.9-7.6$ from the JWST Cycle-1 imaging programmes; PRIMER and JADES. We use BAGPIPES to consistently infer the equivalent widths of their [OIII]+$\rm{H\beta}$ emission lines ($W_\lambda$) and their physical properties. To supplement our high-redshift galaxies, we measure $W_\lambda$ photometrically for a sample of N=253 $z=3.2-3.6$ SFGs selected from the VANDELS spectroscopic survey. Comparing these samples, we find a strong apparent redshift evolution in their median $W_\lambda$, increasing from $W_\lambda =310\pm25\,Å$ in VANDELS to $W_\lambda =630\pm30\,Å$ in our JWST-based sample. In the JWST sample at $z>7$, we find that $W_\lambda$ correlates with both stellar mass and UV luminosity, with high-mass, $M_{ UV}-$faint galaxies producing systematically weaker emission lines. Moreover, we discover a departure from the standard log-normal shape of the $W_\lambda$ distribution, characterised by a more pronounced tail at lower $W_\lambda$, consistent with increasingly bursty star formation. Using $W_\lambda$ as a proxy for $\xi_\rm{ion}$, and UV spectral slope as a proxy for LyC escape fraction ($f_\rm{esc}$), we find a minority of galaxies with high $\xi_\rm{ion}$ and $f_\rm{esc}$ (e.g., $\rm{log(\xi_{ion}/erg^{-1}Hz})\sim25.6$ and $f_\rm{esc}\sim0.15$). However, we find that the LyC photon budget at $z>7$ is dominated by galaxies with more moderate output, close to the sample median of $\rm{log(\xi_{ion}/erg^{-1}Hz})\sim25.3$ and $f_\rm{esc}\sim0.05$. This is consistent with estimates for the number of LyC photons required to power reionization at $z>7$, with no evidence for over or under-production.

Although the circumgalactic medium (CGM) plays a vital role in galaxy evolution, studying the emission from CGM is challenging due to its low surface brightness and the complexities involved in interpreting resonant lines like Ly$\alpha$. The near-infrared coverage, unprecedented sensitivity, and high spatial resolution of the JWST enable us to study the optical strong lines associated with the extended Ly$\alpha$ "nebulae" at redshift 2-3. These lines serve as diagnostic tools to infer the physical conditions in the massive CGM gas reservoir of these systems. In deep medium-band images taken by the JWST, we serendipitously discover the [O III] emission from the CGM surrounding a massive interacting galaxy system at redshift z ~ 2.8, known to be embedded in a bright extended (100 kpc) Ly$\alpha$ "nebula". This is the first time that the [O III] lines are detected from a Ly$\alpha$ "nebula", and the JWST images reveal that the CGM gas actually resides in narrow (~ 2.5 kpc) filamentary structures with strong [O III] emission, tracing the same extent as the Ly$\alpha$ emission. Analysis of the [O III] suggests that the emitting CGM is fully ionized and is energetically dominated by mechanical heating. We also find that the inferred density and pressure are higher than those commonly predicted by simulations of the CGM. We conclude that the observed CGM emission originates from the gas expelled by the episodic feedback processes, cooling down and enriching the CGM, while travelling to a distance of at least 60 kpc. These observations demonstrate how fierce feedback processes shape gas distribution and properties in the CGM around massive halos. While the deep high-resolution imaging opens up a new discovery space for investigating the CGM, it also challenges numerical simulations to explain and reproduce the exquisitely complex structures revealed by the observations.

The James Webb Space Telescope is revolutionising our ability to understand the host galaxies and local environments of high-z quasars. Here we obtain a comprehensive understanding of the host galaxy of the z=7.08 quasar J1120+0641 by combining NIRSpec integral field spectroscopy with NIRCam photometry of the host continuum emission. Our emission line maps reveal that this quasar host is undergoing a merger with a bright companion galaxy. The quasar host and the companion have similar dynamical masses of $\sim10^{10}M_\odot$, suggesting that this is a major galaxy interaction. Through detailed quasar subtraction and SED fitting using the NIRCam data, we obtain an estimate of the host stellar mass of $M_{\ast}\simeq2.6\times10^9M_\odot$, with $M_{*}\simeq5.0\times10^9M_\odot$ for the companion galaxy. Using the H$\beta$ Balmer line we estimate a virial black hole mass of $M_{\rm{BH}}\simeq1.4\times10^9 M_\odot$. Thus, J1120+0641 has an extreme black hole - stellar mass ratio of $M_{\rm{BH}}/M_\ast\simeq0.54$, which is ~3 dex larger than expected by the local scaling relations between black hole and stellar mass. J1120+0641 is powered by an overmassive black hole with the highest reported black hole-stellar mass ratio, in a quasar host that is currently undergoing a major merger -- these new insights highlight the power of JWST for measuring and understanding these extreme first quasars.

M dwarf systems offer a unique opportunity to study terrestrial exoplanetary atmospheres due to their smaller size and cooler temperatures. However, due to the extreme conditions these host stars impose, it is unclear whether their small, close-in rocky planets are able to retain any atmosphere at all. The Hot Rocks Survey aims to answer this question by targeting nine different M dwarf rocky planets spanning a range of planetary and stellar properties. LHS 1478 b orbits an M3-type star, has an equilibrium temperature of Teq = 585 K and experiences an instellation 21 times greater than that of Earth. We observe two secondary eclipses using photometric imaging at 15 um using the Mid-Infrared Instrument on the James Webb Space Telescope (JWST MIRI) to measure thermal emission from the dayside of the planet. We then compare these values to different atmospheric scenarios to evaluate potential heat transport and CO2 absorption signatures. We find a secondary eclipse depth of 146 +/- 56 ppm based on the first observation, while the second observation results in a non-detection due to significantly larger unexplained systematics. Based on the first observation alone, we can reject the null hypothesis of the dark (zero Bond albedo) no atmosphere bare rock model with a confidence level of 3.4 sigma. For an airless body with a Bond albedo of A=0.2, the significance decreases to 2.9 sigma. The secondary eclipse depth is consistent with the majority of atmospheric scenarios we considered, which all involve atmospheres which include different concentrations of CO2, and surface pressures from 0.1 to 10 bar. However, we stress that the two observations from our program do not yield consistent results, and more observations are needed to verify our findings. The Hot Rocks Survey serves as a relevant primer for future endeavors such as the Director's Discretionary Time (DDT) Rocky Worlds program.

Which rocky exoplanets have atmospheres, and which do not? This presumably simply question is the first that must be answered in the pathway toward understanding the prevalence of nearby habitable planets. A mere 6.9 pc from Earth, the LTT 1445A system is the closest transiting M-dwarf system, and its largest known planet, at $\rm 1.31\; R_{\oplus}$ and 424 K, is one of the best rocky targets in which to search for an atmosphere. We use HST/WFC3 transmission spectroscopy with the G280 and G141 grisms to study the spectrum of LTT 1445Ab between $\rm 0.2-1.65\;\mu m$. In doing so, we uncover a flare on the neighboring star LTT 1445C, and report one of the first simultaneous near-UV/optical spectra of an M dwarf flare. The planet spectrum is consistent with a flat line, though the infrared portion displays potential features that could be explained by known opacity sources such as HCN. Some atmospheric retrievals weakly favor ($\sim2\sigma$) an atmosphere over a flat line, but it remains challenging to discern between stellar contamination, an atmosphere, and a featureless spectrum at this time. We do, however, confidently rule out $\leq100\times$ solar metallicity atmospheres. Although stellar contamination retrievals cannot fit the infrared features well, the overall spectrum is consistent with stellar contamination from hot spots, cold spots, or both. Based on the UV/optical data, we place limits on the extent of stellar variability expected in the near-infrared ($30-40$ ppm), which will be critical for upcoming JWST observations of this highly optimal target.

We observed the ultra hot Jupiter WASP-33b with the Spectro-Polarimètre Infra-Rouge on the Canada Fance Hawaii Telescope. Previous observations of the dayside of WASP-33b show evidence of CO and Fe emission indicative of a thermal inversion. We observed its nightside over five Earth-nights to search for spectral signatures of CO in the planet's thermal emission. Our three pre-transit observations and two post-transit observations are sensitive to regions near the morning or evening terminators, respectively. From spectral retrievals, we detect CO molecular absorption in the planet's emission spectrum after transit at $\sim$6.6$\sigma$. This is the strongest ground-based detection of nightside thermal emission from an exoplanet, and only the third ever. CO appearing in absorption suggests that the nightside near the evening terminator does not have a temperature inversion; this makes sense if the dayside inversion is driven by absorption of stellar radiation. On the contrary, we do not detect CO from the morning terminator. This may be consistent with heat advection by an eastward jet. Phase-resolved high-resolution spectroscopy offers an economical alternative to space-based full-orbit spectroscopic phase curves for studying the vertical and horizontal atmospheric temperature profiles of short-period exoplanets.

Carbon-rich (C-rich) stars can be found at all metallicities and evolutionary stages. They are often the result of mass-transfer from a companion, but some of the most metal-poor C-rich objects are likely carrying the imprint of the metal-free First Stars from birth. In this work, we employ a neural network to predict metallicities and carbon abundances for over 10 million stars with Gaia low-resolution XP spectra, down to [Fe/H] = -3.0 and up to [C/Fe] $\approx$ +2. We identify ~2000 high-confidence bright (G<16) carbon-enhanced metal-poor (CEMP) stars with [Fe/H] < -2.0 and [C/Fe] > +0.7. The majority of our C-rich candidates have [Fe/H] > -2.0 and are expected to be binary mass-transfer products, supported by high barium abundances in GALAH and/or their Gaia RUWE and radial velocity variations. We confirm previous findings of an increase in C-rich stars with decreasing metallicity, adopting a definition of $3\sigma$ outliers from the [C/Fe] distribution, although our frequency appears to flatten for -3.0 < [Fe/H] < -2.0 at a level of 6-7%. We also find that the fraction of C-rich stars is low among globular cluster stars (connected to their lower binary fraction), and that it decreases for field stars more tightly bound to the Milky Way. We interpret these last results as evidence that disrupted globular clusters contribute more in the inner Galaxy, supporting previous work. Homogeneous samples like these are key to understanding the full population properties of C-rich stars, and this is just the beginning.

The near-infrared (NIR) emission of the youngest protostars still needs to be characterized to better understand the evolution of their accretion and ejection activity. We analyze James Webb Space Telescope NIRSpec 1.7 -- 5.3 $\mu$m observations of two deeply embedded sources in the S68N protostellar core in Serpens. The North Central (NC) source exhibits a highly obscured spectrum (A_K ~ 4.8 mag) that is modeled with a pre-main-sequence photosphere and a hot disk component. The photospheric parameters are consistent with a young, low-mass photosphere, as suggested by the low surface gravity, log g of 1.94 $\pm$ 0.15 cm s$^{-2}$. The hot disk suggests that accretion onto the central protostellar embryo is ongoing, although prototypical accretion-tracing emission lines HI are not detected. The South Central (SC) source, which is even more embedded (A_K ~ 8 mag; no continuum is detected shortward of 3.6 $\mu$m) appears to be driving the large-scale S68N protostellar outflow, and launches a collimated hot molecular jet detected in \Ht and CO ro-vibrational lines. Shock modeling of the \Ht (ro)vibrational lines establishes that fast $C$-type shocks ($\geq$ 30 km s$^{-1}$), with high pre-shock density ($\geq$ $10^7$ cm$^{-3}$), and strong magnetic field (b ~ 3--10, where $B = b\,\times\,\sqrt{\textrm{n}_{\textrm{H}} (\textrm{cm}^{-3})}\,\mu\textrm{G}$) best match the data. The bright CO fundamental line forest suggests energetic excitation, with the contribution of non-LTE effects, ie irradiation pumping. Detected OH and CH$^{+}$ ro-vibrational lines support this hypothesis. These two Class 0 protostars seem to be in very young evolutionary stages and still have to acquire the bulk of their final stellar masses. These results demonstrate that JWST enables unprecedented diagnostics of these first stages of the protostellar evolutionary phase.

Atomic diffusion with radiative levitation is a major transport process to consider to explain abundance anomalies in Am stars. Radiative accelerations vary from one species to another, yielding different abundance anomalies at the stellar surface. Radiative accelerations can be computed using different methods: some evolution codes use an analytical approximation, while others calculate them from monochromatic opacities. We compared the abundance evolutions predicted using these various methods. Our models were computed with the Toulouse-Geneva evolution code, in which both an analytical approximation (the single-valued parameter method) and detailed calculations from Opacity Project (OP) atomic data are implemented for the calculation of radiative accelerations. The time evolutions of the surface abundances were computed using macroscopic transport processes that are able to reproduce observed Am star surface abundances in presence of atomic diffusion, namely an ad hoc turbulent model or a global mass loss. The radiative accelerations obtained with the various methods are globally in agreement for all the models below the helium convective zone, but can be much greater between the bottom of the hydrogen convective zone and that of helium. The time evolutions of the surface abundances mostly agree within the observational error, but the abundance of some elements can exceed this error for the least massive mass-loss model. The gain in computing time from using analytical approximations is significant compared to sequential calculations from monochromatic opacities for the turbulence models and for the least massive wind model; the gain is small otherwise. Test calculations of turbulence models with the tabulated OPAL opacities yield quite similar abundances as OP for most elements but in a much shorter time, so that determining Am star parameters can be done using a two-step method.

Knowledge about the internal physical structure of stars is crucial to understanding their evolution. The novel binary population synthesis code POSYDON includes a module for interpolating the stellar and binary properties of any system at the end of binary MESA evolution based on a pre-computed set of models. In this work, we present a new emulation method for predicting stellar profiles, i.e., the internal stellar structure along the radial axis, using machine learning techniques. We use principal component analysis for dimensionality reduction and fully-connected feed-forward neural networks for making predictions. We find accuracy to be comparable to that of nearest neighbor approximation, with a strong advantage in terms of memory and storage efficiency. By delivering more information about the evolution of stellar internal structure, these emulators will enable faster simulations of higher physical fidelity with large-scale simulations of binary star population synthesis possible with POSYDON and other population synthesis codes.

We analysed ALMA FIR (1.3 mm) dust continuum and CO emission of 12 starburst galaxies at $z\sim 2.1-3.6$, selected for their extreme brightness in the rest-UV with $M_{\rm UV} = -23.4$ to $-24.7$. We also analysed VLT HAWK-I $H$- and $K_{\rm s}$-band images. The galaxies are characterised by negligible dust attenuations with blue UV spectral slopes ($-2.62$ to $-1.84$), very young stellar populations of $\sim 10$ Myr, and powerful starbursts with a high mean specific star formation rate of $\rm 112~Gyr^{-1}$, placing them $\sim 1.5$~dex above the main sequence at similar redshifts and stellar masses ($M_{\rm stars} \sim (1.5-4.6)\times 10^9~M_{\odot}$). The FIR dust continuum emission revealed in 9 galaxies yields IR luminosities of $(5.9-28.3)\times 10^{11}~L_{\odot}$ and large dust masses barely produced by SNe within the 10~Myr timescale. The CO emission detected in 8 galaxies evidence large molecular gas masses with a mean molecular gas fraction of 82%. The corresponding star formation efficiencies reach $\gtrsim 40$\%, with amazingly short molecular gas depletion timescales between <13 Myr and 71 Myr. These unique properties, never reported in previously studied galaxies, highlight that these galaxies are likely caught at the very beginning of their stellar mass build-up and undergo a very efficient and fast conversion of gas into stars that can only result from the gas collapse within very short free-fall times. We find that the feedback-free starburst model seems to be able to explain the formation of these galaxies. To reconcile the co-spatial FIR dust emission with the UV-bright unattenuated emission, we speculate about radiation-driven outflows which can temporarily remove dust at the location of the starburst and expel dust at large distances in line with the measured large FIR effective radii ($\rm 1.7~kpc - 5~kpc$) in comparison to very compact stellar radii.

We review sulfur chemistry of the gas giant planets and their moons where sulfur compounds are observed. The major S-bearing gas in the upper atmospheres of the giant planets is H2S and is removed from their observable atmospheres by condensation into cloud layers (NH4SH on all four planets and additionally H2S ice on Uranus and Neptune). Any remaining H2S at higher altitudes is destroyed photochemically. Among the moons Io is the world dominated by sulfur. We summarize the sulfur cycle on Io and how pyrovolcanism is spreading sulfur across the Jovian system. Implantation of sulfur into icy surfaces of the other Galilean moons via magnetospheric transfer and radiolysis are major processes affecting the sulfur chemistry on their icy surfaces. On the icy worlds, we are literally looking at the top of the icebergs. Subsurface liquid salty bodies reveal themselves through cryovolcanism on Europa, Ganymede, and Enceladus, where salt deposits are indicated. Subsurface oceans are suspected on several other moons. We summarize the sulfur cycle for the icy Galilean moons. The occurrence of sulfates can be explained by salt exchange reactions of radiolytically produced H2SO4 with brine salts (carbonates and halides), or from a subsurface ocean that has become acidified by uptake of H2SO4 leaked from ice. In the primordial oceans of the moons that accreted with high ice rock ratios, sulfur is expected as sulfide and bisulfide anions and H2S in aqueous solution. Cosmochemical constraints suggest that pyrrhotite, tochilinite and green rusts could be important sulfide bearing compounds found with hydrous silicates such as serpentine, and magnetite on the sea floors. In N-C-rich worlds such as Titan, sulfides such as NH4SH and possibly thiazyl compounds could be important, and sulfates are unstable. Nothing is known about the sulfur chemistry on the Uranian and Neptunian moons.

We conducted systematic radiative transfer (RT) modeling of the Mg II doublet line profiles for 33 low-redshift Lyman continuum (LyC) leakers, and Ly$\alpha$ modeling for a subset of six objects, using a multiphase, clumpy circumgalactic medium (CGM) model. Our RT models successfully reproduced the Mg II line profiles for all 33 galaxies, revealing a necessary condition for strong LyC leakage: high maximum clump outflow velocity ($v_{\rm MgII,\,max} \gtrsim 390\,\rm km\,s^{-1}$) and low total Mg II column density ($N_{\rm MgII,\,tot} \lesssim 10^{14.3}\,\rm cm^{-2}$). We found that the clump outflow velocity and total Mg II column density have the most significant impact on Mg II spectra and emphasized the need for full RT modeling to accurately extract the CGM gas properties. In addition, using archival HST COS/G160M data, we modeled Ly$\alpha$ profiles for six objects and found that their spectral properties do not fully align with the conventional LyC leakage criteria, yet no clear correlation was identified between the modeled parameters and observed LyC escape fractions. We inferred LyC escape fractions based on HI properties from Ly$\alpha$ RT modeling and found that LyC leakage is primarily governed by the number of optically thick HI clumps per sightline ($f_{\rm cl}$). Intriguingly, two galaxies with relatively low observed LyC leakage exhibited the highest RT-inferred LyC escape fractions due to their lowest $f_{\rm cl}$ values, driven by the strong blue peaks of their Ly$\alpha$ emission. Future high-resolution, spatially resolved observations are crucial for resolving this puzzle. Overall, our results support a "picket fence" geometry over a "density-bounded" scenario for the CGM, where a combination of high Mg II outflow velocities and low Mg II column densities may be correlated with the presence of more low-density HI channels that facilitate LyC escape.

Near-infrared (NIR) emission is less affected by dust than UV and optical emission and is therefore useful for studying the properties of dust-obscured galaxies. Although rest-frame NIR observations of high-redshift galaxies have long been made using space telescopes, their structures were unresolved due to the lack of angular resolution. This letter reports the early results from the analysis of high-resolution Pa$\beta$ imaging of the Spiderweb protocluster at $z=2.16$ with the JWST Near Infrared Camera. We investigate radial profiles of Pa$\beta$ lines and rest-frame NIR continua from luminous H$\alpha$-emitting galaxies (HAEs) in the protocluster. Particularly, we compare those of 11 HAEs (N-HAEs) on the star-forming main sequence with those of 8 HAEs (X-HAEs) with X-ray AGNs. Resultant composite Pa$\beta$ line images of N-HAEs indicate significant star formation in galactic disks. In contrast, X-HAEs are dominated by point source components rather than outer star formation, as inferred from our earlier work based on multi-wavelength SED fitting. Given their higher stellar potentials suggested from their rest-frame NIR images, the different characteristics may be driven by the impact of AGN feedback.

We study the generation of high-frequency gravitational waves (GWs) through graviton bremsstrahlung during the decay of inflaton in the post-inflationary universe, focusing on scenarios with a polynomial inflaton potential. Two main reheating channels are considered: decays into bosons (spin 0) and fermions (spin $\frac{1}{2}$). We compute the resulting GW spectra from three-body decays, where the inflaton decays into a pair of daughter particles and a graviton. We numerically compute the GW spectra for various polynomial exponents by employing two distinct approaches: one treating the inflaton as a collection of rest particles and the other treating it as a coherently oscillating classical field. In the former approach, only gravitons with energies below half the inflaton mass are produced, while the latter allows for the production of gravitons with arbitrarily high energies when the potential exponent is 4 or greater. This difference arises because the inflaton's oscillations are no longer described by a single harmonic mode but instead consist of infinitely many harmonic modes with different frequencies. As a result, the GW spectrum exhibits multiple peaks, with these peaks being less pronounced for higher powers of the potential. We also examine the dependence of the GW spectrum on the coupling constant between the inflaton and daughter particles. Our findings suggest that future GW detectors targeting GWs in the GHz band, such as resonant cavities, may have the capability to detect these signals, offering potential insights into the reheating phase of the early universe.

The fraction of dwarf galaxies hosting central, intermediate-mass black holes (IMBHs) at low redshifts is an important observational probe of black hole seeding at high redshift. Detections of nuclear accretion signatures in dwarf galaxies provides strong evidence for the presence of these IMBHs. We develop a Bayesian model to infer the black hole occupation fraction assuming a broken power law Eddington ratio distribution function. Our approach accounts for non-detections, incompleteness, and contamination from star-forming-related emission. We apply this model to galaxies with X-ray data from the Chandra Source Catalog at distances $<50$ Mpc, radio data from the VLA Sky Survey at $< 50$ Mpc, and optical variability data from the Palomar Transient Factory at $z<0.055$. We find a black hole occupation fraction of at least $90$ percent at stellar masses of $M_{\star}=10^8~M_{\odot}$ and at least $39$ percent at $M_{\star} = 10^7~M_{\odot}$ (95\% confidence intervals). We show the resulting black hole mass function. These constraints on the IMBH population have implications for the Laser Interferometer Space Antenna (LISA) mission and for cosmological models of black hole seeding and growth. We also constrain the extremely low luminosity end ($L_{\rm{bol}}\lesssim10^{40}$ erg s$^{-1}$) of the AGN luminosity functions at $z=0$. Our AGN luminosity functions are broadly consistent with an extrapolation of the shallow slope of the AGN luminosity functions from previous work.

Context. The Gaia mission has opened up a new era for the precise astrometry of stars, thus revolutionizing our understanding of the Milky Way. However, beyond a few kiloparseconds from the Sun, parallax measurements become less reliable, and even within 2 kpc, there still exist stars with large uncertainties. Aims. Our aim was to determine the distance and stellar parameters of 521,424 solar-like stars from LAMOST DR9; these stars lacked precise distance measurements (uncertainties higher than 20\% or even without any distance estimations) when checked with Gaia. Methods. We proposed a convolutional neural network (CNN) model to predict the absolute magnitudes, colors, and stellar parameters (T_eff, logg, and [FeH]) directly from low-resolution spectra. For spectra with signal-to-noise ratios at g band (S/N_g) greater than 10, the model achieves a precision of 85 K for T_eff, 0.07 dex for logg, 0.06 dex for [Fe/H], 0.25 mag for M_g, and 0.03 mag for bp-rp. The estimated distances have a median fractional error of 4% with a standard deviation of 8%. Results: We applied the trained CNN model to 521,424 solar-like stars to derive the distance and stellar parameters. Compared with other distance estimation studies and spectroscopic surveys, the results show good consistency. Additionally, we investigated the metallicity gradients of the Milky Way from a subsample, and find a radial gradient ranging from -0.05 < Delta{[Fe/H]}/Delta{R} < 0.0 dex/kpc and a vertical gradient ranging from -0.26 < Delta{[Fe/H]}/Delta{Z} < -0.07dex/kpc Conclusions. We conclude that our method is effective in estimating distances and stellar parameters for solar-like stars with limited astrometric data. Our measurements are reliable for Galactic structure studies and hopefully will be useful for exoplanet researches.

Magnetic reconnection before flux cancellation in the solar photosphere when two opposite-polarity photospheric magnetic fragments are approaching one another is usually modelled by assuming that a small so-called "floating current sheet" forms about a null point or separator that is situated in the overlying atmosphere. Here instead we consider the reconnection that is initiated as soon as the fragments become close enough that their magnetic fields interact. The resulting current sheet, which we term a "fin sheet" extends up from the null point or separator that is initially located in the solar surface. We develop here nonlinear analyses for finite-length models of both fin and floating current sheets that extend the previous models that were limited to short floating current sheets. These enable the length of the current sheet to be calculated in both cases as functions of the separation distance of the sources and the reconnection rate, as well as the rate of heating. Usually, the fin current sheet liberates more energy than a floating current sheet.

Bipolar or multipolar lobes in pre-planetary nebulae (pPNe) often exhibit intertwined outer whorled patterns, resulting from stellar wind matter accumulation during the asymptotic giant branch (AGB) phase. These structures are likely triggered by stellar or substellar companions. We regard that CW Leonis currently stands at a critical transition moment, providing a vivid illustration of the progression from an AGB star in a binary system to a pPN. We have found that CW Leonis has shown significant enhancements in its optical and near-infrared light curves over the past two decades, with the recent Hubble Space Telescope image finally revealing the long-awaited central star. Utilizing an eccentric-orbit binary model, we can reproduce the position-angle dependence of the expansion velocity in the whorled pattern around CW Leonis, suggesting a nearly face-on orbital inclination. Its contradiction to the features in the innermost circumstellar envelope, corresponding to a nearly edge-on inclination, may imply the presence of an additional companion. Our updated theoretical framework explores the complexity of the whorled pattern. Further identifying and monitoring phase-transition candidates at the tip of the AGB will provide valuable insights into the AGB-pPN transition and the role of companions in shaping the morphological evolution of these stellar objects.

Chromospheric plages are distributed between mid-latitude and the Equator and never close to the Poles. Therefore, we suspect that the inclination angle of the solar rotation axis has an impact on the observable chromospheric emission. We reproduce the solar images from any inclination in order to study the effect of the inclination axis on the solar variability by using direct observations of the Sun in the Ca II K line. More than 2700 days of observations since the beginning of the Ca II K observations with USET, in July 2012, were used in our analysis. For each observation day, we produce synoptic maps to map the entire solar surface during a full solar rotation. Then by choosing a given inclination, we generate solar-disk views, representing the segmented brightest structures of the chromosphere (plages and enhanced network), as seen under this inclination. The area fraction are extracted from the masks for each inclination and we compare the evolution of those time series to quantify the impact of the inclination angle. We find a variation of the area fraction between an Equator-on view and a Pole-on view. Our results show an important impact of the viewing angle on the detection of modulation due to the solar rotation. With the dense temporal sampling of USET data, the solar rotation is detectable up to an inclination of about $|i| = 70^{\circ}$ and the solar-cycle modulation is clearly detected for all inclinations, though with a reduced amplitude in polar views. When applying a sparse temporal sampling typical for time series of solar-like stars, the rotational modulation is no longer detected, whatever the inclination. On the other hand, we find that the activity-cycle modulation remains detectable as long as the sampling contains at least 20 observations per year and the cycle amplitude reaches at least 30\% of the solar-cycle amplitude.

This report summarizes what has happened in the mini-workshops entitled ''My Earth in 180 minutes'' organized during the lunch break at the SF2A 2024 conference in Marseille. The project showcased an innovative serious game designed to raise awareness of greenhouse gas (GHG) emissions in astronomical research laboratories. Participants, organized into teams, simulate strategies to reduce their carbon footprints by 50\%, focusing on key astronomical activities such as space instrumentation, data analysis, and laboratory work. The sessions highlight the challenges of achieving significant emissions reductions without disrupting core research activities, such as telescope observations. While the serious game facilitates important discussions on sustainable practices, the results point to the need for broader engagement, adaptation to different cultural contexts, and institutional support. The project highlights the importance of integrating climate action into the academic environment and suggests potential future directions for expanding its impact.

Missing frequency channels pose a problem for estimating $P(k_\perp,k_\parallel)$ the redshifted 21-cm power spectrum (PS) from radio-interferometric visibility data. This is particularly severe for the Murchison Widefield Array (MWA), which has a periodic pattern of missing channels that introduce spikes along $k_\parallel$. The Tracking Tapered Gridded Estimator (TTGE) overcomes this by first correlating the visibilities in the frequency domain to estimate the multi-frequency angular power spectrum (MAPS) $C_\ell(\Delta\nu)$ that has no missing frequency separation $\Delta\nu$. We perform a Fourier transform along $\Delta\nu$ to estimate $P(k_\perp,k_\parallel)$. Considering our earlier work, simulations demonstrate that the TTGE can estimate $P(k_\perp,k_\parallel)$ without any artifacts due to the missing channels. However, the spikes were still found to persist for the actual data, which is foreground-dominated. The current work presents a detailed investigation considering both simulations and actual data. We find that the spikes arise due to a combination of the missing channels and the strong spectral dependence of the foregrounds. Based on this, we propose and demonstrate a technique to mitigate the spikes. Applying this, we find the values of $P(k_\perp,k_\parallel)$ in the region $0.004 \leq k_\perp \leq 0.048\,{\rm Mpc^{-1}}$ and $k_\parallel > 0.35 \,{\rm Mpc^{-1}}$ to be consistent with zero within the expected statistical fluctuations. We obtain the $2\sigma$ upper limit of $\Delta_{\rm UL}^2(k)=(918.17)^2\,{\rm mK^2}$ at $k=0.404\,{\rm Mpc^{-1}}$ for the mean squared brightness temperature fluctuations of the $z=8.2$ epoch of reionization (EoR) 21-cm signal. This upper limit is from just $\sim 17$ minutes of observation for a single pointing direction. We expect tighter constraints when we combine all $162$ different pointing directions of the drift scan observation.

We use PyDynamicaLC, a model using the least number of, and the least correlated, degrees of freedom needed to derive a photodynamical model, to describe some of the smallest -- and lowest TTV (transit timing variations) amplitude -- of the Kepler planets. We successfully analyze 64 systems containing 218 planets, for 88 of which we were able to determine significant masses (to better than $3\sigma$). We demonstrate consistency with literature results over two orders of magnitude in mass, and for the planets that already had literature mass estimations, we were able to reduce the relative mass error by $\sim22\%$ (median value). Of the planets with determined masses 23 are new mass determinations with no previous significant literature value, including a planet smaller and lighter than Earth (KOI-1977.02 / Kepler-345 b). These results demonstrate the power of photodynamical modeling with the appropriately chosen degrees of freedom. This will become increasingly more important as smaller planets are detected, especially as the TESS mission gathers ever longer-baseline light curves and for the analysis of the future PLATO mission data

Electrons accelerated by solar flares and observed as type III solar radio bursts are not only a crucial diagnostic tool for understanding electron transport in the inner heliosphere but also a possible early indication of potentially hazardous space weather events. The electron beams travelling in the solar corona and heliosphere along magnetic field lines generate Langmuir waves and quasilinearly relax towards a plateau in velocity space. The relaxation of the electron beam over the short distance in contrast to large beam-travel distances observed is often referred to as Sturrok's dilemma. Here, we develop a new electron transport model with quasilinear distance/time self-consistently changing in space and time. The model results in a nonlinear advection-diffusion equation for the electron beam density with nonlinear diffusion term that inversely proportional to the beam density. The solution predicts slow super-diffusive (ballistic) spatial expansion of a fast propagating electron beam. The model also provides the evolution of the spectral energy density of Langmuir waves, which determines brightness temperature of plasma radiation in solar bursts. The model solution is consistent with the results of numerical simulation using kinetic equations and can explain some characteristics of type III solar radio bursts.

The precise Doppler method for measuring stellar radial velocities (RV) is a fundamental technique in modern astronomy. This method records a star's spectrum and detects periodic Doppler shifts in its spectral features, which indicate the gravitational influences induced by orbiting companions. The Doppler technique has yielded remarkable successes in exoplanet detection, uncovering a diverse array of planetary systems ranging from hot Jupiters to Neptune-mass planets and super-Earths. Having led to the discovery of over 1100 exoplanets, the RV method is the most effective approach for measuring orbital geometries and minimum masses, which are of fundamental importance when accessing planet formation and evolution scenarios. However, exoplanet detection via precise RV measurements poses significant challenges, including managing various sources of interference, such as instrumental errors, and mitigating spurious Doppler shifts induced by phenomena like stellar activity. Key to this technique's advancement is instrumental calibration methods, notably precise calibration methods and ultra-stable spectrographs. This technique holds promise in systematically exploring the domain of Jovian analogs, rocky and icy planets within the habitable zones of their parent stars, and providing crucial follow-up observations for transiting candidates detected by space missions. The synergy between transit and Doppler measurements of exoplanets, when feasible, has provided a comprehensive set of orbital and physical parameters for exoplanets, such as the dynamical mass and mean planet density, which is enhancing our understanding of their internal composition. Additionally, ongoing efforts aim to improve the RV technique further by developing more stable calibration techniques aimed at detecting Earth-like analogs around Solar-type stars that require cm/s RV precision.

There are about 6000 stars, that can be seen with the naked eye and have been observed for centuries for various purposes. More modern investigations using advanced telescopes show that our Milky Way, a quite common galaxy, consists of about 100 -- 400 billion stars. And, it is estimated that there are between 200 billion to 2 trillion galaxies in the observable universe -- all of them consist mostly of stars, and sending observable signals which also represents nothing more than a superposition of the light of individual stars. So we can conclude that the most common observable objects in the Universe are $\textit{stars}$. In this chapter, we focus on the long history of the observation of stars (compared to studies in other fields of science) to find out more about the nature of these objects.

Strong gravity in the immediate vicinity of compact objects (e.g., black holes, neutron stars) necessitates inclusion of general relativistic effects. Traditionally, pseudo-Newtonian potential representation of gravity were favored to simulate the fluid motion in this region since that reduced the calculation complexity. However, with the advent of easily implementable, reliable numerical algorithms and computer hardware, more and more research groups are shifting towards the numerical solutions of general relativistic fluid dynamics equations. In this work, we report our progress on the development of such simulation tool and present results of sub-Keplerian accretion flow onto black holes.

Brown dwarfs provide a unique opportunity to study atmospheres and their physical and chemical processes with high precision, especially in temperature ranges relevant to exoplanets. In this study, we performed high-resolution ($R \sim 70,000$) spectroscopy using Subaru/IRD of Gl 229 B, the first discovered T-type (T7.0p) brown dwarf, which orbits an M1V host star at a separation of 33 au. We conducted atmospheric retrieval on the reduced $H$-band spectrum using the high-resolution spectrum model compatible with automatic differentiation and GPU, ExoJAX. In contrast to previous retrieval studies on medium-resolution spectra, we obtained a C/O ratio consistent with that of the host star, aligning with the expected formation process for such a massive brown dwarf. Additionally, based on the strong constraint on temperature from the high-resolution spectrum and previously measured photometric magnitude, our analysis indicates that Gl 229 B is a binary, which was also proposed by Brandt et al. (2021). Finally, we validated current molecular line lists by leveraging the obtained high-precision, high-resolution spectrum of this warm ($\sim 900$ K) atmosphere. This study highlights the importance of observing companion brown dwarfs as benchmark objects for establishing characterization techniques for low-mass objects and enhancing our understanding of their atmospheres, given the wealth of available information and the relative ease of observation.

We present the unplanned detection of a white dwarf companion to the star HD 218261 in mid-infrared (10-16 $\mu$m) observations with JWST/MIRI. This star was observed as a calibrator for coronagraphic observations of the exoplanet host HR 8799. HD 218261 B has only previously been detected by Gaia, and only in visible light. We confidently detect the companion in the mid-infrared, where it is less luminous than the primary by a factor of ~10$^4$. The visible and mid-infrared photometry are consistent with a white dwarf of $T_\text{eff}\approx10000$ K, $M\approx0.8 M_\odot$, though observation of its optical spectrum is required to precisely constrain its physical parameters. These results demonstrate that precise mid-infrared photometry of white dwarf companions to bright stars can be obtained with MIRI, opening up new possibilities for studying white dwarfs in close binaries.

Stellar ultraviolet (UV) emission serves as a crucial indicator for estimating magnetic activity and evaluating the habitability of exoplanets orbiting stars. In this paper, we present a straightforward method to derive stellar photospheric UV emission for F to M main-sequence stars. By using PARSEC models, we establish relations between near-UV (NUV) and far-UV (FUV) magnitudes from the Galaxy Evolution Explorer (GALEX), NUV magnitudes from the China Space Station Telescope, and stellar effective temperatures and Gaia BP$-$RP color for different metallicities. Together with the observed sample, we find that for NUV emission, the photospheric contribution to the observed flux is less than 20% for M stars, around 10% to 70% for G stars, and ranges from 30% to 85% for G and F stars. For FUV emission, the photospheric contribution is less than $10^{-6}$ for M stars, below $10^{-4}$ for K stars, around $10^{-4}$ to 10% for G stars, and between 6% and 50% for F stars. Our work enables the simple and effective determination of stellar excess UV emission and the exploration of magnetic activity.

A detection of the 21-cm signal power spectrum from the Epoch of Reionization is imminent, thanks to consistent advancements from telescopes such as LOFAR, MWA, and HERA, along with the development of SKA. In light of this progress, it is crucial to expand the parameter space of simulations used to infer astrophysical properties from this signal. In this work, we explore the role of cosmological parameters such as the Hubble constant $H_0$ and the matter clustering amplitude $\sigma_8$, whose values as provided by measurements at different redshifts are in tension. We run $N$-body simulations using GADGET-4, and post-process them with the reionization simulation code POLAR, that uses L-GALAXIES to include galaxy formation and evolution properties and GRIZZLY to execute 1-D radiative transfer of ionizing photons in the intergalactic medium (IGM). We compare our results with the latest JWST observations and explore which astrophysical properties for different cosmologies are necessary to match the observed UV luminosity functions at redshifts $z = 10$ and 9. Additionally, we explore the impact of these parameters on the observed 21-cm signal power spectrum, focusing on the redshifts within the range of LOFAR 21-cm signal observations ($z \approx 8.5-10$). Despite differences in cosmological and astrophysical parameters, the 21-cm power spectrum at these redshifts agrees with presently observed upper limits. This suggests the need for broader physical parameter spaces for inference modeling to account for all models that agree with observations. However, we also propose stronger constraining power by using a combination of galactic and IGM observables.

The extreme environments of transiting close-in exoplanets in highly-eccentric orbits are ideal for testing exo-climate physics. Spectroscopically resolved phase curves not only allow for the characterization of their thermal response to irradiation changes but also unveil phase-dependent atmospheric chemistry and dynamics. We observed a partial phase curve of the highly-eccentric close-in giant planet HAT-P-2b ($e=0.51,M=9M_{\rm{Jup}}$) with the Wide Field Camera 3 aboard the Hubble Space Telescope. Using these data, we updated the planet's orbital parameters and radius, and retrieved high-frequency pulsations consistent with the planet-induced pulsations reported in Spitzer data. We found that the peak in planetary flux occurred at $6.7\pm0.6$ hr after periastron, with a heating and cooling timescales of $9.0^{+3.5}_{-2.1}$ hr, and $3.6^{+0.7}_{-0.6}$ hr, respectively. We compare the light-curve to various 1-dimensional and 3-dimensional forward models, varying the planet's chemical composition. The strong contrast in flux increase and decrease timescales before and after periapse indicates an opacity term that emerges during the planet's heating phase, potentially due to more H$^{-}$ than expected from chemical equilibrium models. The phase-resolved spectra are largely featureless, that we interpret as indicative an inhomogeneous dayside. However, we identified an anomalously high flux in the spectroscopic bin coinciding with the hydrogen Paschen $\beta$ line and that is likely connected to the planet's orbit. We interpret this as due to shock heating of the upper atmosphere given the short timescale involved, or evidence for other star-planet interactions.

Hot subdwarf (SD) stars are the stripped cores of red giant stars in transition to the white dwarf sequence. The B-type subdwarfs (sdB) are powered by helium fusion in the core, more evolved ones (sdO) by shell burning. Low mass SDs may evolve through this stage without any support by nuclear fusion. Because the loss of the giants' envelopes is likely caused by mass transfer in binaries, hot SDs are test beds for close-binary evolution through stable and unstable Roche lobe overflow, common envelope formation and ejection as well as mergers. Many classes of hot SDs can be identified according to surface composition, binarity, magnetism, pulsation characteristics and population membership, including members of globular clusters. Observed binaries show a wide spread of orbital periods from 20 minutes to more than 1,000 days with white dwarf or main sequence companions. The closest systems qualify as type Ia supernova progenitors and LISA detectable gravitational wave sources. High-precision light curves from Kepler and TESS combined with radial velocity curves are used to derive masses, while asteroseismology adds information on the internal structure, slow rotation, and synchronization. Gaia's parallax measurements now allow us to place the stars in the Hertzsprung-Russell diagram and derive stellar parameters by combining them with multi-band photometry. The stellar radius can be determined to high precision this way. Newton's law can then be used to derive masses if accurate surface gravities are available. Large-scale spectroscopic surveys will provide atmospheric parameters for large samples of stars, allowing the mass distributions for the diverse subtypes to be established. These are crucial for testing binary synthesis models and constraining poorly known parameters such as the common envelope efficiency as well as the critical threshold mass-ratio for mass transfer stability.

Having a detailed knowledge of the intracluster medium (ICM) to infer the exact cluster physics such as the cluster dynamical state is crucial for cluster-based cosmological studies. This knowledge limits the accuracy and precision of mass estimation, a key parameter for such studies. In this paper, we conduct an in-depth analysis of cluster ACT-CL J0240.0+0116 using a multi-wavelength approach, with a primary focus on high angular resolution Sunyaev-Zeldovich (SZ) thermal component observations obtained under the NIKA2 Sunyaev-Zeldovich Large Programme (LPSZ). We create composite images using NIKA2, X-ray, and optical galaxy number density maps. The results reveal distinct signs of disturbance within the cluster with the distributions of gas and member galaxies that do not overlap. We also find suggestions of an inflow of matter onto the cluster from the southwestern direction. Ultimately, we classify the cluster as disturbed, using morphological indicators derived from its SZ, X-ray, and optical image. The cluster SZ signal is also contaminated by a strong central point source. We adopt different approaches to handling this contaminant and find the estimates of our pressure and hydrostatic mass profiles robust to the point source mitigation model. The cluster hydrostatic mass is estimated at $4.25^{+0.50}_{-0.45\, } \times 10^{14} \,\mathrm{M}_{\odot}$ for the case where the point source was masked. These values are consistent with the mass estimated using only X-ray data and with those from previous SZ studies of the Atacama cosmology telescope (ACT) survey, with improved precision on the mass estimate. Our findings strongly suggest that ACT-CL J0240.0+0116 is a disturbed cluster system, and the detailed observations and derived values serve as a compelling case study for the capabilities of the LPSZ in mapping the cluster ICM with high precision.

To study the dust dynamics in the dust trapping vortices in the protoplanetary disk around MWC~758, we analyzed the 1.3 mm continuum images of the MWC~758 disk obtained with the Atacama Large Millimeter/submillimeter Array (ALMA) in 2017 and 2021. We detect proper motions of 22 mas and 24 mas in the two dust clumps at radii of 0\farcs32 and 0\farcs54 in the disk on the plane of the sky, respectively. On the assumption that the dust clumps are located in the disk midplane, the velocities of the observed proper motions along the azimuthal direction of the inner and outer dust clumps are sub- and super-Keplerian, respectively, and both have angular velocities corresponding to the Keplerian angular velocity at a radius of $0\farcs46\pm0\farcs04$. This deviation from the Keplerian motion is not expected in the conventional theory of vortices formed by the Rossby wave instability. The observed non-Keplerian proper motions of the dust clumps are unlikely due to the disk warp and eccentricity, nor be associated with any predicted planets. The two dust clumps are likely spatially coincident with the infrared spirals. In addition, we detect the changes in the intensity profiles of the dust clumps over the four-year span. Therefore, we suggest that the observed proper motions are possibly due to changes in the density distributions in the dust clumps caused by their interaction with the spirals in the disk.

Observations of the high-$z$ Universe from JWST have revealed a new population of bright, early galaxies. A robust statistical interpretation of this data requires fast forward models that account for uncertainties in galaxy evolution and incorporate observational systematic effects. We present a probabilistic framework for population synthesis of high-$z$ galaxies and inference of their properties. Our framework is based on the semi-analytic galaxy-formation model Galacticus. To infer the astrophysical parameters governing high-$z$ galaxy evolution, we analyze JWST data from the CEERS and NGDEEP surveys and calculate the likelihood of observing individual objects in apparent magnitude--redshift space, for $z\geq8.5$ galaxy candidates. We include observational selection effects due to limited survey volume and depth, as well as photometric redshift uncertainties. We recover the posterior probability distribution for parameters describing star formation and outflow rates. We place an upper limit on the star formation timescale of $500~\mathrm{Myr}$ at a disk velocity of $50~\mathrm{km\ s}^{-1}$, and we infer a characteristic velocity at which the outflow mass-loading factor is $\sim 1$ of $150^{+280}_{-60}~\mathrm{km\ s}^{-1}$, both at $95\%$ confidence. Marginalizing over our astrophysical model, we find that galaxies in CEERS and NGDEEP data occupy halos with virial masses $10^{10\pm 0.5}~M_{\mathrm{\odot}}$ at $8.5\leq z\leq 12$, at $95\%$ confidence. The star formation timescale preferred by our fit is relatively short compared to typical values at lower redshifts, consistent with previous findings. The modeling and analysis framework presented here can enable systematic tests of high-$z$ galaxies' dust content, initial mass functions, and star formation burstiness in the future.

Recent studies showed that observations of line emission from shocks in YSO jets require a substantial amount of ionization of the pre-shock matter. Photoionization from X-ray emitted close to the central source may be responsible for the initial ionization fraction. The aim of our work is to study the effect of X-ray photoionization, coming from the vicinity of the central star, on the ionization fraction inside the jet that can be advected at large distances. For this purpose we have performed axisymmetric MHD jet launching simulations including photoionization and optically thin losses using PLUTO. For typical X-ray luminosities in classical T-Tauri stars, we see that the photoionization is responsible for ionizing to 10 % -20 % the jet close to the star.

Scylla is a deep Hubble Space Telescope survey of the stellar populations, interstellar medium and star formation in the LMC and SMC. As a pure-parallel complement to the Ultraviolet Legacy Library of Young Stars as Essential Standards (ULLYSES) survey, Scylla obtained 342 orbits of ultraviolet (UV) through near-infrared (IR) imaging of the LMC and SMC with Wide Field Camera 3. In this paper, we describe the science objectives, observing strategy, data reduction procedure, and initial results from our photometric analysis of 96 observed fields. Although our observations were constrained by ULYSSES primary exposures, we imaged all fields in at least two filters (F475W and F814W), and 64% of fields in at least three and as many as seven WFC3 filters spanning the UV to IR. Overall, we reach average 50% completeness of $m_{\rm F225W}=26.0$, $m_{\rm F275W}=26.2$, $m_{\rm F336W}=26.9$, $m_{\rm F475W}=27.8$, $m_{\rm F814W}=25.5$, $m_{\rm F110W}=24.7$, and $m_{\rm F160W}=24.0$ Vega magnitudes in our photometric catalogs, which is faintward of the ancient main sequence turnoff in all filters. The primary science goals of Scylla include characterizing the structure and properties of dust in the MCs, as well as their spatially-resolved star formation and chemical enrichment histories. Our images and photometric catalogs, which represent the widest-area coverage of MCs with HST photometry to date, are available as a high-level science product at the Barbara A. Mikulski Archive for Space Telescopes.

Ultra-hot Jupiters (UHJs) are gas giant planets orbiting close to their host star, with equilibrium temperatures exceeding 2000 K, and among the most studied planets in terms of their atmospheric composition. Thanks to a new generation of ultra-stable high-resolution spectrographs, it is possible to detect the signal from the individual lines of the species in the exoplanetary atmospheres. We employed two techniques in this study. First, we used transmission spectroscopy, which involved examining the spectra around single lines of FeII. Then we carried out a set of cross-correlation studies for two UHJs: KELT-9b and KELT-20b. Both planets orbit fast-rotating stars, which resulted in the detection of the strong Rossiter-McLaughlin (RM) effect and center-to-limb variations in the transmission spectrum. These effects had to be corrected to ensure a precise analysis. Using the transmission spectroscopy method, we detected 21 single lines of FeII in the atmosphere of KELT-9b. All of the detected lines are blue-shifted, suggesting strong day-to-night side atmospheric winds. The cross-correlation method leads to the detection of the blue-shifted signal with a signal-to-noise ratio (S/N) of 13.46. Our results are in agreement with models based on non-local thermodynamical equilibrium (NLTE) effects, with a mean micro-turbulence of nu_mic = 2.73 +/- 1.5 km/s and macro-turbulence of nu_mac = 8.22 +/- 3.85 km/s. In the atmosphere of KELT-20b, we detected 17 single lines of FeII. Considering different measurements of the systemic velocity of the system, we conclude that the existence of winds in the atmosphere of KELT-20b cannot be determined conclusively. The detected signal with the cross-correlation method presents a S/N of 11.51. The results are consistent with NLTE effects, including means of nu_mic = 3.04 +/- 0.35 km/s and nu_mac = 6.76 +/- 1.17 km/s.

In this work, we investigate the properties of a population of cool clouds in simulated galaxy outflows. Using data from the CGOLS isolated galaxy simulations, we generate catalogues of $\sim 10^5$ clouds. We describe the impact of two different supernova feedback models -- a centrally concentrated starburst and disk-wide distributed star formation -- on the resulting cloud population. In both cases we find that the mass distribution function $dN/dM \propto M^{-2}$, in good agreement with model predictions of turbulent fragmentation. We explore how cloud properties change with distance from the galaxy and find no qualitative distinction between the two feedback modes, although significant quantitative differences exist in attributes such as the total number of clouds, their densities, etc. We further show that both internal cloud velocities and cloud-cloud relative velocities are described well by properties of turbulent motion, despite significant bulk radial velocities. Finally, we investigate the distribution of cloud sizes in the context of recent theoretical arguments about cloud survival in winds. We find that proposed cloud survival criteria are a good predictor of cloud survival, in both the case where clouds are primarily destroyed and the case where cloud growth occurs in the outflow.

The mechanism responsible for the prompt gamma-ray emission of a gamma-ray burst continues to remain an enigma. The detailed analysis of the spectrum of GRB 131014A observed by the $Fermi$ gamma ray burst monitor and Large Area Telescope has revealed an unconventional spectral shape that significantly deviates from the typical Band function. The spectrum exhibits three distinctive breaks and an extended power law at higher energies. Furthermore, the lower end of the spectrum aligns with power-law indices greater than -0.5, and in the brightest region of the burst, these values approach +1. The lowest spectral break is thereby found to be consistent with a blackbody. These observed spectral characteristics strongly suggest the radiation process to be inverse Compton scattering in an optically thin region. Applying the empirical fit parameters for physical modeling, we find that the kinetic energy of the GRB jet of bulk Lorentz factor, $\Gamma \sim 400$, gets dissipated just above the photosphere, approximately at a radius of $\sim 10^{14}$ cm. The electrons involved in this process are accelerated to a power-law index of $\delta = -1.5$, and the minimum electron Lorentz factor, $\gamma_{min}$, is approximately $3$. In summary, this study provides a comprehensive identification and detailed modeling of optically thin inverse Compton scattering in the prompt emission of GRB 131014A.

We present the all-sky catalogue of absolute positions and estimates of correlated flux density of 21,942 compact radio sources determined from processing interferometric visibility data of virtually all VLBI observing sessions at 2 to 23 GHz from 72 programs suitable for absolute astrometry collected for 30 years. We used a novel technique of generation of a dataset of fused observables that allowed us to incorporate all available data in our analysis. The catalogue is the most complete and most precise to date. It forms the foundation and reference for positional astronomy, space geodesy, space navigation, population analysis of active galactic nuclei (AGNs), and provides calibrators for phase-referencing for differential astrometry and VLBI astrophysical observations. Its accuracy was evaluated through a detailed accounting systematic errors, rigorous decimation tests, comparison of different datasets, and comparison with other catalogues. The catalogue samples preferentially the AGNs with a strong contemporary parsec-scale synchrotron emission. Its milliarcsecond level positional accuracy allows association of these AGNs with detections in a wide range of the electromagnetic spectrum from low-frequency radio to gamma-rays and high energy neutrinos. We describe the innovative data processing and calibration technique in full detail, report the in depth analysis of random and systematic positional errors, and provide the list of associations with large surveys at different wavelengths.

Recently, possible hints of parity violation have been observed in the connected galaxy four-point correlation function. Although the true origin of the signal from the analysis has been debated, should they have a physical origin, they might point to primordial non-Gaussianity and would be evidence of new physics. In this work, we examine the single-field slow-roll model of inflation within chiral scalar-tensor theories of modified gravity. These theories extend the Chern-Simons one by including parity-violating operators containing first and second derivatives of the non-minimally coupled scalar (inflaton) field. This model is capable of imprinting parity-violating signatures in late-time observables, such as the galaxy four-point correlation function. We perform an analysis of the graviton-mediated scalar trispectrum of the gauge-invariant curvature perturbation $\zeta(t,\mathbf{x})$. We estimate that for a set of parameters of the theory it is possible to produce a signal-to-noise ratio for the parity-violating part of the trispectrum of order one without introducing modifications to the single-field slow-roll setup. Even if the signal found in the analysis turns out to be spurious or if no parity violation is ever detected in the galaxy four-point correlation function, our analysis can be used to constrain the free parameters of these theories.

We present the first results of the JWST Emission Line Survey (JELS). Utilising the first NIRCam narrow-band imaging at 4.7$\mu$m, over 63 arcmin$^{2}$ in the PRIMER/COSMOS field, we have identified 609 emission line galaxy candidates. From these, we robustly selected 35 H$\alpha$ star-forming galaxies at $z \sim 6.1$, with H$\alpha$ star-formation rates ($\rm{SFR_{H\alpha}}$) of $\sim0.9-15\ \rm{M_{\odot} yr^{-1}}$. Combining our unique H$\alpha$ sample with the exquisite panchromatic data in the field, we explore their physical properties and star-formation histories, and compare these to a broad-band selected sample at $z\sim 6$ to offer vital new insight into the nature of high-redshift galaxies. UV-continuum slopes ($\beta$) are considerably redder for our H$\alpha$ sample ($\langle\beta\rangle\sim-1.92$) compared to the broad-band sample ($\langle\beta\rangle\sim-2.35$). This is not due to dust attenuation as our H$\alpha$ sample is relatively dust-poor (median $A_V=0.23$); instead, we argue that the reddened slopes could be due to nebular continuum. We compared $\rm{SFR_{H\alpha}}$ and the UV-continuum-derived $\rm{SFR_{UV}}$ to SED-fitted measurements averaged over canonical timescales of 10 and 100 Myr ($\rm{SFR_{10}}$ and $\rm{SFR_{100}}$). We find an increase in recent SFR for our sample of H$\alpha$ emitters, particularly at lower stellar masses ($<10^9 \rm{M_{\odot}}$). We also find that $\rm{SFR_{H\alpha}}$ strongly traces SFR averaged over 10 Myr timescales, whereas the UV-continuum over-predicts SFR on 100 Myr timescales at low stellar masses. These results point to our H$\alpha$ sample undergoing `bursty' star formation. Our F356W $z \sim 6$ sample show a larger scatter in $\rm{SFR_{10}/SFR_{100}}$ across all stellar masses, highlighting how narrow-band photometric selections of H$\alpha$ emitters are key to quantifying the burstiness of star-formation activity.

The dynamical properties of the bodies in the outer Solar System hold information regarding the planets' orbital histories. In early Solar System numerical simulations, where chaos is a primary driver, it is difficult to explore parameter space in a systematic way. In such simulations, stable configurations are hard to come by, and often require special fine-tuning. In addition, it is infeasible to run suites of well-resolved, realistic simulations with massive particles to drive planetary evolution where enough particles remain to represent the transneptunian populations to robustly statistically compare with observations. To complement state of the art full N-body simulations, we develop a method to artificially control each planet's orbital elements independently from each other, which when carefully applied, can be used to test a wider suite of models. We modify two widely used publicly available N-body integrators: (1) the C code, REBOUND and (2) the FORTRAN code, Mercury. We show how the application of specific fictitious forces within numerical integrators can be used to tightly control planetary evolution to more easily explore migration and orbital excitation and damping. This tool allows us to replicate the impact a massive planetesimal disk would have on the planets, without actually including the massive planetesimals, thus decreasing the chaos and simulation runtime. We highlight an appropriate application that shows the impact of Neptune's eccentricity damping and radial outward migration on Uranus' eccentricity

We present the first spatially-resolved infrared images of supernova remnants (SNRs) in M33 with the unprecedented sensitivity and resolution of JWST. We analyze 43 SNRs in four JWST fields: two covering central and southern M33 with separate NIRCam (F335M, F444W) and MIRI (F560W, F2100W) observations, one $\sim$5 kpc-long radial strip observed with MIRI F770W, and one covering the giant HII region NGC 604 with multiple NIRCam and MIRI broad/narrowband filters. Of the 21 SNRs in the MIRI field, we found three clear detections (i.e., identical infrared and \ha morphologies), and six partial-detections, implying a detection fraction of 43\% in these bands. In contrast, only one SNR (out of 16) is detectable in the NIRCam field. One of the SNRs, L10-080, is a potential candidate for having freshly-formed ejecta dust, based on its size and centrally-concentrated 21 \mum emission. Two SNRs near NGC 604 have strong evidence of molecular (H$_2$) emission at 4.7 \mum, making them the farthest known SNRs with visible molecular shocks. Five SNRs have F770W observations, with the smaller younger objects showing tentative signs of emission, while the older, larger ones have voids. Multi-wavelength data indicate that the clearly-detected SNRs are also among the smallest, brightest at other wavelengths (\ha, radio and X-ray), have the broadest line widths (H$\alpha$ FWHM$\sim$250-350 km/s), and the densest environments. No strong correlation with star-formation histories are seen, with the clearly-detected SNRs having both high-mass ($\sim$35 \Msun) and low-mass ($\lesssim$10 \Msun) progenitors.

Binary population synthesis (BPS) is an essential tool for extracting information about massive binary evolution from gravitational-wave (GW) detections of compact object mergers. It has been successfully used to constrain the most likely permutations of evolution models among hundreds of alternatives, while initial condition models, in contrast, have not yet received the same level of attention. Here, we introduce BOSSA, a detailed initial sampling code including a set of 192 initial condition permutations for BPS that capture both "invariant" and "varying" models, the latter accounting for a possible metallicity- and star formation rate (SFR)-dependence of the initial mass function (IMF); as well as correlations between the initial primary mass, orbital period, mass ratio and eccentricity of binaries. We include 24 metallicity-specific cosmic star formation history (cSFH) models and propose two alternate models for the mass-dependent binary fraction. We build a detailed pipeline for time-evolving BPS, such that each binary has well-defined initial conditions, and we are able to distinguish the contributions from populations of different ages. We discuss the meaning of the IMF for binaries and introduce a refined initial sampling procedure for component masses. We also discuss the treatment of higher-order multiple systems when normalizing a binary sample. In particular, we argue for how a consistent interpretation of the IMF implies that this is not the distribution from which any set of component masses should be independently drawn, and show how the individual IMF of primaries and companions is expected to deviate from the full IMF.