Papers for Monday, Jun 03 2024

The Hubble Space Telescope (HST) has been our most prolific tool to study exoplanet atmospheres. As the age of JWST begins, there is a wealth of HST archival data that is useful to strengthen our inferences from JWST. Notably, HST/STIS and its 0.3-1 $\mu$m wavelength coverage extends past JWST's 0.6 $\mu$m wavelength cutoff and holds an abundance of potential information: alkali (Na, K) and molecular (TiO, VO) species opacities, aerosol information, and the presence of stellar contamination. However, time series observations with HST suffer from significant instrumental systematics and can be highly dependent on choices made during the transit fitting process. This makes comparing transmission spectra of planets with different data reduction methodologies challenging, as it is difficult to discern if an observed trend is caused by differences in data reduction or underlying physical processes. Here, we present the Sculpting Hubble's Exoplanet Legacy (SHEL) program, which aims to build a consistent data reduction and light curve analysis methodology and associated database of transmission spectra from archival HST observations. In this paper, we present the SHEL analysis framework for HST/STIS and its low-resolution spectroscopy modes, G430L and G750L. We apply our methodology to four notable hot Jupiters: WASP-39 b, WASP-121 b, WASP-69 b, and WASP-17 b, and use these examples to discuss nuances behind analysis with HST/STIS. Our results for WASP-39 b, WASP-121 b, and WASP-17 b are consistent with past publications, but our analysis of WASP-69 b differs and shows evidence of either a strong scattering slope or stellar contamination. The data reduction pipeline and tutorials are available on Github.

Primordial black holes (PBHs) in the mass range $10^{-16}-10^{-11}~M_\odot$ may constitute all the dark matter. We show that gravitational microlensing of bright x-ray pulsars provide the most robust and immediately implementable opportunity to uncover PBH dark matter in this mass window. As proofs of concept, we show that the currently operational NICER telescope can probe this window near $10^{-14}~M_\odot$ with just two months of exposure on the x-ray pulsar SMC-X1, and that the forthcoming STROBE-X telescope can probe complementary regions in only a few weeks. These times are much shorter than the year-long exposures obtained by NICER on some individual sources. We take into account the effects of wave optics and the finite extent of the source, which become important for the mass range of our PBHs. We also provide a spectral diagnostic to distinguish microlensing from transient background events and to broadly mark the PBH mass if true microlensing events are observed. In light of the powerful science case, i.e., the imminent discovery of dark matter searchable over multiple decades of PBH masses with achievable exposures, we strongly urge the commission of a dedicated large broadband telescope for x-ray microlensing. We derive the microlensing reach of such a telescope by assuming sensitivities of detector components of proposed missions, and find that with hard x-ray pulsar sources PBH masses down to a few $10^{-17}~M_\odot$ can be probed.

In many cases accretion proceeds from disks onto planets, stars, white dwarfs, and neutron stars via a boundary layer, a region of intense shear where gas transitions from a near-Keplerian speed to that of the surface. These regions are \textit{not} susceptible to the common magnetorotational and Kelvin-Helmholtz instabilities, and instead global modes generated by supersonic shear instabilities are a leading candidate to govern transport in these regions. This work investigates the dynamics of these systems under a range of thermodynamic conditions, surveying both disk sound speeds and cooling rates. Very fast or very slow cooling has little effect on wave dynamics: in the fast-cooling limit, waves propagate in an effectively isothermal manner, and in the slow limit wave propagation is effectively adiabatic. However, when the cooling timescale is comparable to the wave period, wave damping becomes extreme. In cases with intermediate cooling rates, mass and angular momentum transport can be suppressed by orders of magnitude compared to isothermal and uncooled cases. Cooling in accretion disks leads to a preference for wavenumbers near and below the Mach number of the disk; the corresponding lower frequencies can (in non-isothermal systems) couple to gravity modes within the star, potentially causing low-frequency variability such as dwarf nova and quasi-periodic oscillations in accreting systems.

We developed a model for the star formation history (SFH) of super-early galaxies and applied it to GS-z14-0, the most distant galaxy known, located at $z=14.32$ (294 million years after the Big Bang). The SFH, starting at $z=26.7$, is complex. Initially ($z>18$), the galaxy experiences feedback-regulated phases that are bursty, relatively faint (reaching $M_{\rm UV}=-18.4$), and unattenuated. When dust shielding allows for a smooth star formation rate (SFR), the galaxy quickly becomes heavily obscured. During this obscured phase, which lasts for approximately 20% of the total star-forming time, 70% of the observed stars are formed. Super-early galaxies in this phase should be detectable by ALMA. Twenty-six million years before observation, as the galaxy becomes super-Eddington, a powerful radiation-driven outflow clears most of the dust and significantly reduces the SFR by a factor of seven, from $100 \to 15\ M_\odot \rm yr^{-1}$. The galaxy transitions into a "blue monster" dominating the bright end of the UV luminosity function. When the outflow ceases due to decreased dust opacity, the galaxy relaxes into a post-starburst phase, in which it is currently observed. Our model accurately reproduces all the observed and inferred properties of the galaxy. The analysis of this extreme system opens exciting opportunities for studying the beginnings of the luminous Universe.

Folk wisdom dictates that a lower bound on the dark matter particle mass, $m$, can be obtained by demanding that the de Broglie wavelength in a given galaxy must be smaller than the virial radius of the galaxy, leading to $m\gtrsim 10^{-22}\text{ eV}$ when applied to typical dwarf galaxies. This lower limit has never been derived precisely or rigorously. We use stellar kinematical data for the Milky Way satellite galaxy Leo II to self-consistently reconstruct a statistical ensemble of dark matter wavefunctions and corresponding density profiles. By comparison to a data-driven, model-independent reconstruction, and using a variant of the maximum mean discrepancy as a statistical measure, we determine that a self-consistent description of dark matter in the local Universe requires $m>2.2 \times 10^{-21}\,\mathrm{eV}\;\mathrm{(CL>95\%)}$. This lower limit is free of any assumptions pertaining to cosmology, microphysics (including spin), or dynamics of dark matter, and only assumes that it is predominantly composed of a single bosonic particle species.

Triple body systems are prevalent in nature, from planetary to stellar to supermassive black hole scales. In a hierarchical triple system, oscillations of the inner orbit's eccentricity and inclination can be induced on secular timescales. Over many cycles, the octupole-level terms in the secular equations of motion can drive the system to extremely high eccentricities via the Eccentric Kozai-Lidov (EKL) mechanism. The overall decrease in the inner orbit's pericenter distance has potentially dramatic effects for realistic systems, such as tidal disruption. We present an analytical approximation in the test particle limit to describe the step-wise eccentricity evolution of the inner orbit. We also integrate the equations of motion and present approximations in both the test particle limit and in the general case to describe the overall octupole-level time evolution of the eccentricity. The analytical approximations are compared to numerical simulations to show that the models accurately describe the form and timescale of the secular descent from large distances to a close-encounter distance (e.g., the Roche limit). By circumventing the need for numerical simulations to obtain the long-term behavior, these approximations can be used to readily estimate descent timescales for populations of systems. We demonstrate by calculating rates of EKL-driven migration for Hot Jupiters in stellar binaries.

We have performed a systematic search for galaxy-scale strong lenses using Hyper Suprime-Cam imaging data, focusing on lenses in overdense environments. To identify these lens candidates, we exploit our neural network from HOLISMOKES VI, which is trained on realistic gri mock-images as positive examples, and real images as negative examples. Compared to our previous work, we lower the i-Kron radius limit to >0.5". This results in an increase by around 73 million sources to more than 135 million images. During our visual multi-stage grading of the network candidates, we now inspect simultaneously larger stamps (80"x80") to identify large, extended arcs cropped in the 10"x10" cutouts, and classify additionally their overall environment. Here we also reinspect our previous lens candidates and classify their environment. Using these 546 visually identified lens candidates, we further define various criteria by exploiting extensive and complementary photometric redshift catalogs, to select the candidates in overdensities. In total, we identified 24 grade-A and 138 grade-B candidates with either spatially-resolved multiple images or extended, distorted arcs in the new sample. Furthermore, with our different techniques, we identify in total 237/546 lens candidates in a cluster-like or overdense environment, containing only 49 group- or cluster-scale re-discoveries. These results demonstrate the feasibility of downloading and applying network classifiers to hundreds of million cutouts, necessary in the upcoming era of big data from deep, wide-field imaging surveys like Euclid and the Rubin Observatory Legacy Survey of Space and Time, while leading to a sample size that can be inspected by humans. These networks, with false-positive rates of ~0.01%, are very powerful tools to identify such rare galaxy-scale strong lensing systems, while also aiding in the discovery of new strong lensing clusters.

Narrowband galaxy surveys have recently gained interest as a promising method to achieve the necessary accuracy on the photometric redshift estimate of individual galaxies for stage-IV cosmological surveys. One key advantage is the ability to provide higher spectral resolution information about galaxies that should allow a more accurate and precise estimation of galaxy stellar population properties. However, the impact of adding narrow-band photometry on the stellar population properties estimate is largely unexplored. The scope of this work is two-fold: on one side, leveraging the predictive power of broad-band and narrow-band data to infer galaxy physical properties such as stellar masses, ages, star formation rates and metallicities. On the other hand, evaluating the improvement of performance in estimating galaxy properties when we use narrow-band data instead of broad-band. In this work we measure the stellar population properties of a sample of galaxies in the COSMOS field for which both narrowband and broadband data are available. In particular, we employ narrowband data from PAUS and broad-band data from CFHTLS. We use two different spectral energy distribution fitting codes to measure galaxy properties, namely CIGALE and Prospector. We find that the increased spectral resolution of narrow-band photometry does not yield a substantial improvement on constraining galaxy properties using spectral energy distribution fitting. Still we find that we obtain a more diverse distribution of metallicities and dust optical depths with cigale when employing the narrowband data. The effect is not as prominent as expected, which we relate this to the low narrowband SNR of a majority of the galaxies, the respective drawbacks of both codes as well as the coverage only in the optical regime. The measured properties are afterwards compared to the COSMOS2020 catalogue, showing good agreement.

The Parameterised Post-Newtonian (PPN) approach is the default framework for performing precision tests of gravity in nearby astrophysical systems. In recent works we have extended this approach for cosmological applications, and in this paper we use observations of the anisotropies in the Cosmic Microwave Background to constrain the time variation of the PPN parameters $\alpha$ and $\gamma$ between last scattering and the present day. We find their time-averages over cosmological history should be within $\sim 20\%$ of their values in GR, with $\bar{\alpha}=0.89^{+0.08}_{-0.09}$ and $\bar{\gamma}=0.90^{+0.07}_{-0.08}$ at the $68\%$ confidence level. We also constrain the time derivatives of these parameters, and find that their present-day values should be within a factor of two of the best Solar System constraints. Many of these results have no counter-part from Solar System observations, and are entirely new constraints on the gravitational interaction. In all cases, we find that the data strongly prefer $\bar{\alpha}\simeq \bar{\gamma}$, meaning that observers would typically find local gravitational physics to be compatible with GR, despite considerable variation of $\alpha$ and $\gamma$ being allowed over cosmic history. This study lays the groundwork for future precision tests of gravity that combine observations made over all cosmological and astrophysical scales of length and time.

Large Language Models (LLMs) are shifting how scientific research is done. It is imperative to understand how researchers interact with these models and how scientific sub-communities like astronomy might benefit from them. However, there is currently no standard for evaluating the use of LLMs in astronomy. Therefore, we present the experimental design for an evaluation study on how astronomy researchers interact with LLMs. We deploy a Slack chatbot that can answer queries from users via Retrieval-Augmented Generation (RAG); these responses are grounded in astronomy papers from arXiv. We record and anonymize user questions and chatbot answers, user upvotes and downvotes to LLM responses, user feedback to the LLM, and retrieved documents and similarity scores with the query. Our data collection method will enable future dynamic evaluations of LLM tools for astronomy.

We use COSMIC, a galaxy population synthesis code, to investigate how metallicity affects the rate of formation of massive stars with a closely orbiting compact object companion, the suggested progenitors of radio loud long gamma-ray bursts. We present the evolution time of these systems at different metallicities, and how the formation rates of these systems are anti-correlated with metallicity. In particular, these systems occur about 10 times more frequently in at metallicities between $Z = 2\times 10^{-4}$ and $2 \times 10^{-3}$, compared to those between $Z = 2\times 10^{-3}$ and $2 \times 10^{-2}$. This work serves as a prerequisite to predicting the global rates of these systems as a function of redshift, ultimately giving crucial insight into our understanding of the progenitors of long gamma-ray bursts and their evolution over cosmic time.

The Primordial Inflation Explorer (PIXIE) is an Explorer-class mission concept to measure the energy spectrum and linear polarization of the cosmic microwave background (CMB). A single cryogenic Fourier transform spectrometer compares the sky to an external blackbody calibration target, measuring the Stokes I, Q, U parameters to levels ~200 Jy/sr in each 2.65 degree diameter beam over the full sky, in each of 300 frequency channels from 28 GHz to 6 THz. With sensitivity over 1000 times greater than COBE/FIRAS, PIXIE opens a broad discovery space for the origin, contents, and evolution of the universe. Measurements of small distortions from a CMB blackbody spectrum provide a robust determination of the mean electron pressure and temperature in the universe while constraining processes including dissipation of primordial density perturbations, black holes, and the decay or annihilation of dark matter. Full-sky maps of linear polarization measure the optical depth to reionization at nearly the cosmic variance limit and constrain models of primordial inflation. Spectra with sub-percent absolute calibration spanning microwave to far-IR wavelengths provide a legacy data set for analyses including line intensity mapping of extragalactic emission and the cosmic infrared background amplitude and anisotropy. We describe the PIXIE instrument sensitivity, foreground subtraction, and anticipated science return from both the baseline 2-year mission and a potential extended mission.

SNe Ia are used to determine the distance-redshift relation and build the Hubble diagram. Neglecting their host-galaxy peculiar velocities (PVs) may bias the measurement of cosmological parameters. The smaller the redshift, the larger the effect is. We use realistic simulations of SNe Ia observed by the Zwicky Transient Facility (ZTF) to investigate the effect of different methods to take into account PVs. We study the impact of neglecting galaxy PVs and their correlations in an analysis of the SNe Ia Hubble diagram. We find that it is necessary to use the PV full covariance matrix computed from the velocity power spectrum to take into account the sample variance. Considering the results we have obtained using simulations, we determine the PV systematic effects in the context of the ZTF DR2 SNe Ia sample. We determine the PV impact on the intercept of the Hubble diagram, $a_B$, which is directly linked to the measurement of $H_0$. We show that not taking into account PVs and their correlations results in a shift of the $H_0$ value of about $1.0$km.s$^{-1}$.Mpc$^{-1}$ and a slight underestimation of the $H_0$ error bar.

Among the nearly 30,000 known near-Earth asteroids (NEAs), only tens of them possess Earth co-orbital characteristics with semi-major axes $\sim$1 au. In particular, 469219 Kamo`oalewa (2016 HO3), upcoming target of China's Tianwen-2 asteroid sampling mission, exhibits a meta-stable 1:1 mean-motion resonance with Earth. Intriguingly, recent ground-based observations show that Kamo`oalewa has spectroscopic characteristics similar to space-weathered lunar silicates, hinting at a lunar origin instead of an asteroidal one like the vast majority of NEAs. Here we use numerical simulations to demonstrate that Kamo`oalewa's physical and orbital properties are compatible with a fragment from a crater larger than 10--20 km formed on the Moon in the last few million years. The impact could have ejected sufficiently large fragments into heliocentric orbits, some of which could be transferred to Earth 1:1 resonance and persist today. This leads us to suggest the young lunar crater Giordano Bruno (22 km diameter, 1--10 Ma age) as the most likely source, linking a specific asteroid in space to its source crater on the Moon. The hypothesis will be tested by the Tianwen-2 mission when it returns a sample of Kamo`oalewa. And the upcoming NEO Surveyor mission will possibly help us to identify such a lunar-derived NEA population.

New longwave HgCdTe detectors are critical to upcoming plans for ground-based infrared astronomy. These detectors, with fast-readouts and deep well-depths, will be key components of extremely large telescope instruments and therefore must be well understood prior to deployment. We analyze one such HgCdTe detector, a Teledyne Imaging Sensors GeoSnap, at the University of Michigan. We find that the properties of the GeoSnap are consistent with expectations from analysis of past devices. The GeoSnap has a well-depth of 2.75 million electrons per pixel, a read noise of 360 e-/pix, and a dark current of 330,000 e-/s/pix at 45 K. The device experiences 1/f noise which can be mitigated relative to half-well shot noise with modest frequency image differencing. The GeoSnap's quantum efficiency is calculated to be 79.7 +- 8.3 % at 10.6 microns. Although the GeoSnap's bad pixel fraction, on the order of 3%, is consistent with other GeoSnap devices, close to a third of the bad pixels in this detector are clustered in a series of 31 "leopard" spots spread across the detector plane. We report these properties and identify additional analyses that will be performed on future GeoSnap detectors.

Cosmic filaments are the main transport channels of matter in the Megaparsec universe, and represent the most prominent structural feature in the matter and galaxy distribution. Here we describe and define the physical and dynamical nature of cosmic filaments. It is based on the realization that the complex spatial pattern and connectivity of the cosmic web are already visible in the primordial random Gaussian density field, in the spatial pattern of the primordial tidal and deformation eigenvalue field. The filaments and other structural features in the cosmic web emerging from this are multistream features and structural singularities in phase-space. The caustic skeleton formalism allows a fully analytical classification, identification, and treatment of the nonlinear cosmic web. The caustic conditions yield the mathematical specification of weblike structures in terms of the primordial deformation tensor eigenvalue and eigenvector fields, in which filaments are identified -- in 2D -- with the so-called cusp caustics. These are centered around points that are maximally stretched as a result of the tidal force field. The resulting mathematical conditions represent a complete characterization of filaments in terms of their formation history, dynamics, and orientation. We illustrate the workings of the formalism on the basis of a set of constrained $N$-body simulations of protofilament realizations. These realizations are analyzed in terms of spatial structure, density profiles, and multistream structure and compared to simpler density or potential field saddle point specifications. The presented formalism, and its 3D generalization, will facilitate the mining of the rich cosmological information contained in the observed weblike galaxy distribution, and be of key significance for the analysis of cosmological surveys such as SDSS, DESI, and Euclid.

Hot Dust-Obscured Galaxies (Hot DOGs) are a rare population of hyper-luminous infrared galaxies discovered by the WISE mission. Despite the significant obscuration of the AGN by dust in these systems, pronounced broad and blue-shifted emission lines are often observed. Previous work has shown that 8 Hot DOGs, referred to as Blue-excess Hot DOGs (BHDs), present a blue excess consistent with type 1 quasar emission in their UV-optical SEDs, which has been shown to originate from the light of the obscured central engine scattered into the line of sight. We present an analysis of the rest-frame optical emission characteristics for 172 Hot DOGs through UV-MIR SED modeling and spectroscopic details, with a particular focus on the identification of BHDs. We find that while the optical emission observed in Hot DOGs is in most cases dominated by a young stellar population, 26% of Hot DOGs show a significant enough blue excess emission to be classified as BHDs. Based on their broad CIV and MgII lines, we find that the $M_{\rm BH}$ in BHDs range from $10^{8.7}$ to $10^{10} \ M_{\odot}$. When using the same emission lines in regular Hot DOGs, we find the $M_{\rm BH}$ estimates cover the entire range found for BHDs while also extending to somewhat lower values. This agreement may imply that the broad lines in regular Hot DOGs also originate from scattered light from the central engine, just as in BHDs, although a more detailed study would be needed to rule out an outflow-driven nature. Similar to $z\sim 6$ quasars, we find that Hot DOGs sit above the local relation between stellar and black hole mass, suggesting either that AGN feedback has not yet significantly suppressed the stellar mass growth in the host galaxies, or that they will be outliers of the relation when reaching $z$=0.

Galactic white dwarf binaries (WDBs) and black hole binaries (BHBs) will be gravitational wave (GW) sources for LISA. Their detection will provide insights into binary evolution and the evolution of our Galaxy through cosmic history. Here, we make predictions of the expected WDB and BHB population within our Galaxy. We combine predictions of the compact remnant binary populations expected by stellar evolution by using the detailed Binary Population and Spectral Synthesis code (BPASS) with a Milky Way analogue galaxy model from the Feedback In Realistic Environment (FIRE) simulations. We use \textsc{PhenomA} and \textsc{LEGWORK} to simulate LISA observations. Both packages make similar predictions that on average four Galactic BHBs and 673 Galactic WDBs above the signal-to-noise ratio (SNR) threshold of 7 after a four-year mission. We compare these predictions to earlier results using the Binary Star Evolution (BSE) code with the same FIRE model galaxy. We find that BPASS predicts a few more LISA observable Galactic BHBs and a twentieth of the Galactic WDBs. The differences are due to the different physical assumptions that have gone into the binary evolution calculations. These results indicate that the expected population of compact binaries that LISA will detect depends very sensitively on the binary population synthesis models used and thus observations of the LISA population will provide tight constraints on our modelling of binary stars. Finally, from our synthetic populations we have created mock LISA signals that can be used to test and refine data processing methods of the eventual LISA observations.

Abell 1758 (z~0.278) is a galaxy cluster composed of two structures: A1758N and A1758S, separated by ~2.2 Mpc. The northern cluster is itself a dissociative merging cluster that has already been modelled by dedicated simulations. Recent radio observations revealed the existence of a previously undetected bridge connecting A1758N and A1758S. New simulations are now needed to take into account the presence of A1758S. We wish to evaluate which orbital configuration would be compatible with a bridge between the clusters. Using N-body hydrodynamical simulations that build upon the previous model, we explore different scenarios that could have led to the current observed configuration. Five types of orbital approaches were tested: radial, tangential, vertical, post-apocentric, and outgoing. We found that the incoming simulated scenarios are generally consistent with mild enhancements of gas density between the approaching clusters. The mock X-ray images exhibit a detectable bridge in all cases. Compared to measurements of Chandra data, the amplitude of the X-ray excess is overestimated by a factor of ~2--3 in the best simulations. The scenario of tangential approach proved to be the one that best matches the properties of the profiles of X-ray surface brightness. The scenarios of radial approach of vertical approach are also marginally compatible.

High-precision isotopic measurements of meteorites revealed that they are classified into non-carbonaceous (NC) and carbonaceous (CC) meteorites. One plausible scenario for achieving this grouping is the early formation of Jupiter because massive planets can create gaps that suppress the mixing of dust across the gap in protoplanetary disks. However, the efficiency of this suppression by the gaps depends on dust size and the strength of turbulent diffusion, allowing some fraction of the dust particles to leak across the Jovian gap. In this study, we investigate how isotopic ratios of NC and CC meteorites are varied by the dust leaking across the Jovian gap in the solar nebula. To do this, we constructed a model to simulate the evolution of the dust size distribution and the $^{54}$Cr-isotopic anomaly $\varepsilon^{54}$Cr in isotopically heterogeneous disks with Jupiter. Assuming that the parent bodies of NC and CC meteorites are formed in two dust-concentrated locations inside and outside Jupiter's orbit, referred to as the NC reservoir and CC reservoir, we derive the temporal variation of $\varepsilon^{54}$Cr at the NC and CC reservoir. Our results indicate that substantial contamination of CC materials occurs at the NC reservoir in the fiducial run. Nevertheless, the values of $\varepsilon^{54}$Cr at the NC reservoir and the CC reservoir in the run are still consistent with those of NC and CC meteorites formed around 2 Myrs after the formation of calcium-aluminum-rich inclusions. Moreover, this dust leakage causes a positive correlation between the $\varepsilon^{54}$Cr value of NC meteorites and the accretion ages of their parent bodies.

Gravitational-wave observations of binary neutron-star (BNS) mergers have the potential to revolutionize our understanding of the nuclear equation of state (EOS) and the fundamental interactions that determine its properties. However, Bayesian parameter estimation frameworks do not typically sample over microscopic nuclear-physics parameters that determine the EOS. One of the major hurdles in doing so is the computational cost involved in solving the neutron-star structure equations, known as the Tolman-Oppenheimer-Volkoff (TOV) equations. In this paper, we explore approaches to emulating solutions for the TOV equations: Multilayer Perceptrons (MLP), Gaussian Processes (GP), and a data-driven variant of the reduced basis method (RBM). We implement these emulators for three different parameterizations of the nuclear EOS, each with a different degree of complexity represented by the number of model parameters. We find that our MLP-based emulators are generally more accurate than the other two algorithms whereas the RBM results in the largest speedup with respect to the full, high-fidelity TOV solver. We employ these emulators for a simple parameter inference using a potentially loud BNS observation, and show that the posteriors predicted by our emulators are in excellent agreement with those obtained from the full TOV solver.

We present confusion-limited SCUBA-2 450-$\mu$m observations in the COSMOS-CANDELS region as part of the JCMT Large Program, SCUBA-2 Ultra Deep Imaging EAO Survey (STUDIES). Our maps at 450 and 850 $\mu$m cover an area of 450 arcmin$^2$. We achieved instrumental noise levels of $\sigma_{\mathrm{450}}=$ 0.59 mJy beam$^{-1}$ and $\sigma_{\mathrm{850}}=$ 0.09 mJy beam$^{-1}$ in the deepest area of each map. The corresponding confusion noise levels are estimated to be 0.65 and 0.36 mJy beam$^{-1}$. Above the 4 (3.5) $\sigma$ threshold, we detected 360 (479) sources at 450 $\mu$m and 237 (314) sources at 850 $\mu$m. We derive the deepest blank-field number counts at 450 $\mu$m, covering the flux-density range of 2 to 43 mJy. These are in agreement with other SCUBA-2 blank-field and lensing-cluster observations, but are lower than various model counts. We compare the counts with those in other fields and find that the field-to-field variance observed at 450 $\mu$m at the $R=6^\prime$ scale is consistent with Poisson noise, so there is no evidence of strong 2-D clustering at this scale. Additionally, we derive the integrated surface brightness at 450 $\mu$m down to 2.1 mJy to be $57.3^{+1.0}_{-6.2}$~Jy deg$^{-2}$, contributing to (41$\pm$4)\% of the 450-$\mu$m extragalactic background light (EBL) measured by COBE and Planck. Our results suggest that the 450-$\mu$m EBL may be fully resolved at $0.08^{+0.09}_{-0.08}$~mJy, which extremely deep lensing-cluster observations and next-generation submillimeter instruments with large aperture sizes may be able to achieve.

We investigate the influence of star formation and instantaneous AGN feedback processes on the ionized gas velocity dispersion in a sample of 1285 emission-line galaxies with stellar masses $\log\,(M_*/M_{\odot}) \geq 9$ from the integral-field spectroscopy SAMI Galaxy Survey. We fit both narrow and broad emission line components using aperture spectra integrated within one effective radius, while ensuring the elimination of velocity differences between the spectra of individual spaxels. Our analysis reveals that 386 (30%) galaxies can be adequately described using a single emission component while 356 (28%) galaxies require two (broad and narrow) components. Galaxies characterized by high mass, elevated star formation rate surface density, or type-2 AGN-like emissions tend to feature an additional broad emission-line component, leading to their classification as double-component galaxies. We explore the correlations between $M_*$ and gas velocity dispersions, highlighting that the prominence of the broad component significantly contributes to elevating the gas velocity dispersion. Galaxies displaying AGN-like emission based on optical definitions show enhanced gas velocity dispersions. In star-forming galaxies, both stellar mass and star-formation rate surface density substantially contribute to the velocity dispersion of the narrow component. Increased star-forming activity appears to elevate the velocity dispersion of the narrow component. The broad component exhibits a weaker dependence on stellar mass and is primarily driven by galactic outflows. We suggest that strong star forming activity leads to the formation of a broad emission-line component, but the impact on inflating gas velocity dispersion is moderate. On the other hand, AGN-driven outflows appear to be a more important contributor to the elevated velocity dispersion of the ionized gas.

Radio observations from normal pulsars indicate that the coherent radio emission is excited by curvature radiation from charge bunches. In this review we provide a systematic description of the various observational constraints on the radio emission mechanism. We have discussed the presence of highly polarized time samples where the polarization position angle follow two orthogonal well defined tracks across the profile, that closely match the rotating vector model in an identical manner. The observations also show the presence of circular polarization, with both the right and left handed circular polarization seen across the profile. Other constraints on the emission mechanism is provided by the detailed measurements of the spectral index variation across the profile window, where the central part of the profile, corresponding to the core component, has a steeper spectrum than the surrounding cones. Finally, the detailed measurements of the subpulse drifting behaviour can be explained by considering the presence of non-dipolar field on the stellar surface and the formation of the Partially screened Gap (PSG) above the polar cap region. The PSG gives rise to a non-stationary plasma flow, that has a multi-component nature, consisting of highly energetic primary particles, secondary pair plasma and iron ions discharged from the surface, with large fragmentation resulting is dense plasma clouds and lower density inter-cloud regions. The physical properties of the outflowing plasma and the observational constraints lead us to consider coherent curvature radiation as the most viable explanation for the emission mechanism in normal pulsars, where propagation effects due to adiabatic walking and refraction are largely inconsequential.

Low-mass galaxies at high redshifts are the building blocks of more massive galaxies at later times and are thus key populations for understanding galaxy formation and evolution. We have made deep narrow-band observations for two protoclusters and the general field in COSMOS at $z$ $\sim$ 2. In a clumpy young protocluster, USS1558$-$003, at $z$ = 2.53, we find many star-forming galaxies well above the star-forming main sequence of field galaxies at the low-mass end ($M_{\star}/\mathrm{M_{\odot}}<10^{8.9}$). This suggests that some environmental effects may be at work in low-mass galaxies in high-density regions to enhance their star formation activities. In the core of this protocluster, we also find that enhanced star formation activity of middle-mass galaxies ($10^{8.9} < M_{\star}/\mathrm{M_{\odot}} < 10^{10.2}$) while such trends are not observed in a more mature protocluster, PKS1138$-$262 at $z$ = 2.16. We expect these activities to be mainly due to galaxy mergers/interactions and differences in the amount of cold gas accretion. As one piece of evidence, we show that the star formation activity within individual galaxies in the protoclusters is more centrally concentrated than those in the field. This is probably due to the enhanced interactions between galaxies in the protocluster, which can reduce the angular momentum of the gas, drive the gas towards the galaxy center, and lead to a central starburst.

We investigate the role of carbon monoxide ice in the chemical evolution of prestellar cores using astrochemical rate equation models. We constrain the ratios of the binding energies on CO ice and H$_{2}$O ice for a series of adsorbates deemed important in diffusive chemistry on H$_{2}$O ices. We later include these ratios in our chemical reaction network model, where the binding and diffusion energies of icy species vary as a function of the surface composition. When the surface coverage of CO increases, the model shows an enhancement of O-bearing complex organic molecules, especially those formed from the intermediate products of CO hydrogenation (e.g. HCO) and CH$_{3}$/CH$_{2}$. Because the binding energy of CH$_{3}$/CH$_{2}$ is in the right range, its diffusion rate increases significantly with CO coverage. At $T>$14 K and with less influence, enhanced diffusion of HCO also contributes to the increase of the abundances of COM. We find, however, that chemistry is not always enhanced on CO ice and that the temperature and cosmic ray ionization rate of each astronomical object is crucial for this particular chemistry, revealing a highly non-trivial behavior that needs to be addressed on a per-case basis. Our results are highly relevant in the context of interstellar ice observations with JWST.

Ephemeral Fast Radio Bursts (FRBs) must be powered by some of the most energetic processes in the Universe. That makes them highly interesting in their own right and as precise probes for estimating cosmological parameters. This field thus poses a unique challenge: FRBs must be detected promptly and immediately localised and studied based only on that single millisecond-duration flash. The problem is that the burst occurrence is highly unpredictable and that their distance strongly suppresses their brightness. Since the discovery of FRBs in single-dish archival data in 2007, detection software has evolved tremendously. Pipelines now detect bursts in real-time within a matter of seconds, operate on interferometers, buffer high-time and frequency resolution data, and issue real-time alerts to other observatories for rapid multi-wavelength follow-up. In this paper, we review the components that comprise a FRB search software pipeline, we discuss the proven techniques that were adopted from pulsar searches, we highlight newer, more efficient techniques for detecting FRBs, and we conclude by discussing the proposed novel future methodologies that may power the search for FRBs in the era of big data astronomy.

Context: MATISSE, the mid-infrared spectro-imaging instrument of VLTI, was designed to deliver its advertised performance when paired with an external second generation fringe tracker. Science observation started in 2019, demonstrating imaging capabilities and faint science target observations. Now, The GRAVITY fringe tracker stabilizes the MATISSE fringes which allows using all spectroscopic modes and improves sensitivity and data accuracy. Aims: We present how the MATISSE and GRAVITY instruments were adapted to make the GRAVITY fringe tracker work with MATISSE, under the umbrella of the aptly-named GRA4MAT project, led by ESO in collaboration with the two instrument consortia. Methods: We detail the software modifications needed to implement an acquisition and observing sequence specific to GRA4MAT, including simultaneous fringe tracking and chopping and a narrow off-axis capability inspired by the galactic center and exoplanet capability of GRAVITY. We explain the modified data collection and reduction processes. We show how we leveraged the recent fringe tracker upgrade to implement features specific to its use with MATISSE, e.g. fringe jumps mitigation with an improved group delay control and simultaneous fringe tracking and chopping with a new state machine. Results: We successfully demonstrate significant improvements to the MATISSE instrument. Observations can now be performed at higher spectral resolutions of up to $R\sim3300$ and across the full LM bands at once. Long detector integration times, made possible with stabilized fringes, have improved the LM-bands sensitivity by a factor of 10. Low flux biases in coherently-reduced N-band data have been eliminated. The L-band transfer function is now higher and more stable. We finally illustrate the scientific potential of GRA4MAT with a preview of the first exoplanet observation made by MATISSE on $\beta$ Pictoris b.

The Be X-ray binary EXO 2030+375 went through its third recorded giant outburst from June 2021 to early 2022. We present the results of both spectral and timing analysis based on NICER monitoring, covering the 2-10 keV flux range from 20 to 310 mCrab. Dense monitoring with observations carried out about every second day and a total exposure time of 160 ks allowed us to closely track the source evolution over the outburst. Changes in spectral shape and pulse profiles showed a stable luminosity dependence during the rise and decline. The same type of dependence has been seen in past outbursts. The pulse profile is characterized by several distinct peaks and dips. The profiles show a clear dependence on luminosity with a stark transition at a luminosity of 2x10^36 erg/s, indicating a change in the emission pattern. Using relativistic ray-tracing, we demonstrate how anisotropic beaming of emission from an accretion channel with constant geometrical configuration can give rise to the observed pulse profiles over a range of luminosities.

It is generally hard to put robust constraints on progenitor masses of supernovae (SNe) and remnants (SNRs) observationally, while they offer tantalizing clues to understanding explosion mechanisms and mass distribution. Our recent study suggests that ``shell merger'', which is theoretically expected for stellar evolution, can appreciably affect final yields of inter-mediate mass elements (IMEs; such as Ne, Mg, and Si). In light of this, here we report results of X-ray spectral analysis of a Galactic SNR G359.0-0.9, whose abundance pattern may possibly be anomalous according to a previous study. Our spectroscopy using all the available data taken with XMM-Newton reveals that this remnant is classified as Mg-rich SNRs because of its high Mg-to-Ne ratio (Z_Mg/Z_Ne=1.90+0.27-0.19; mass ratio 0.66+0.09-0.07) and conclude that the result cannot be explained without the shell merger. By comparing the observation with theoretical calculations, we prefer the so-called Ne-burning shell intrusion and in this case the progenitor mass M_ZAMS is likely <15M_sun. We confirm the result also by our new molecular line observations with the NRO-45 m telescope: G359.0-0.9 is located in the Scutum-Centaurus arm (2.66--2.94 kpc) and in this case the resultant total ejecta mass ~6.8M_sun is indeed consistent with the above estimate. Our method using mass ratios of IMEs presented in this paper will become useful to distinguish the type of the shell merger, the Ne-burning shell intrusion and the O-burning shell merger, for future SNR studies.

The interstellar medium of galaxies is composed of multiple phases, including molecular, atomic, and ionized gas, as well as dust. Stars are formed within this medium from cold molecular gas clouds, which collapse due to their gravitational attraction. Throughout their life, stars emit strong radiation fields and stellar winds, and they can also explode as supernovae at the end of their life. These processes contribute to stirring the turbulent interstellar medium and regulate star formation by heating up, ionizing, and expelling part of the gas. However, star formation does not proceed uniformly throughout the history of the Universe and decrease by an order of magnitude in the last ten billion years. To understand this winding-down of star formation and assess possible variations in the efficiency of star formation, it is crucial to probe the molecular gas reservoirs from which stars are formed. In this article following my presentation at the 10th International Conference on Frontiers of Plasma Physics and Technology held in Kathmandu from 13-17 March 2023, I review some aspects of the multiphase interstellar medium and star formation, with an emphasis on the interplay between neutral and ionized phases, and present recent and ongoing observations of the molecular gas content in typical star-forming galaxies across cosmic time and in different environments. I also present some of our understanding of star-forming galaxies from theoretical models and simulations.

The Fermi Gamma-Ray Space Telescope reveals two large bubbles in the Galaxy, extending nearly symmetrically $\sim50^{\circ}$ above and below the Galactic center (GC). Previous simulations of bubble formation invoking active galactic nucleus (AGN) jets have assumed that the jets are vertical to the Galactic disk; however, in general, the jet orientation does not necessarily correlate with the rotational axis of the Galactic disk. Using three-dimensional special relativistic hydrodynamic simulations including cosmic rays (CRs) and thermal gas, we show that the dense clumpy gas within the Galactic disk disrupts jet collimation ("failed jets" hereafter), which causes the failed jets to form hot bubbles. Subsequent buoyancy in the stratified atmosphere renders them vertical to form the symmetric Fermi and eROSITA bubbles (collectively, Galactic bubbles). We find that (1) despite the relativistic jets emanated from the GC are at various angles $\le45^{\circ}$ with respect to the rotational axis of the Galaxy, the Galactic bubbles nonetheless appear aligned with the axis; (2) the edge of the eROSITA bubbles corresponds to a forward shock driven by the hot bubbles; (3) followed by the forward shock is a tangling contact discontinuity corresponding to the edge of the Fermi bubbles; (4) assuming a leptonic model we find that the observed gamma-ray bubbles and microwave haze can be reproduced with a best-fit CR power-law spectral index of 2.4; The agreements between the simulated and the observed multi-wavelength features suggest that forming the Galactic bubbles by oblique AGN failed jets is a plausible scenario.

Pair cascades in polar cap regions of neutron stars are considered to be an essential process in various models of coherent radio emissions of pulsars. The cascades produce pair plasma bunch discharges in quasi-periodic spark events. The cascade properties, and therefore also the coherent radiation, depend strongly on the magnetospheric plasma properties and vary significantly across and along the polar cap. It is furthermore still uncertain from where the radio emission emanates in polar cap region. We investigate the generation of electromagnetic waves by pair cascades and their propagation in the polar cap for three representative inclination angles of a magnetic dipole, $0^\circ$, $45^\circ$, and $90^\circ$. 2D particle-in-cell simulations that include quantum-electrodynamic pair cascades are used in a charge limited flow from the star surface. We found that the discharge properties are strongly dependent on the magnetospheric current profile in the polar cap and that transport channels for high intensity Poynting flux are formed along magnetic field lines where the magnetospheric currents approach zero and where the plasma cannot carry the magnetospheric currents. There, the parallel Poynting flux component is efficiently transported away from the star and may eventually escape the magnetosphere as coherent radio waves. The Poynting flux decreases faster with the distance from the star in regions of high magnetospheric currents. Our model shows that no process of energy conversion from particles to waves is necessary for the coherent radio wave emission. Moreover, the pulsar radio beam does not have a cone structure, but rather the radiation generated by the oscillating electric gap fields directly escapes along open magnetic field lines in which no pair creation occurs.

We present an analysis of rest-frame UV colors of 17,243 galaxies at $z\sim2-4$ in the HST UVCANDELS fields: GOODS-N, GOODS-S, COSMOS, and EGS. Here, we study the rest-frame UV spectral slope, $\beta$, measured via model spectra obtained via spectral energy distribution (SED) fitting, $\beta_{SED}$, and explore its correlation with various galaxy parameters (photometric redshift, UV magnitude, stellar mass, dust attenuation, star formation rate [SFR], and specific SFR) obtained via SED fitting with Dense Basis. We also obtain measurements for $\beta$ via photometric power-law fitting and compare them to our SED-fit-based results, finding good agreement on average. While we find little evolution in $\beta$ with redshift from $z=2-4$ for the full population, there are clear correlations between $\beta$ (and related parameters) when binned by stellar mass. For this sample, lower stellar mass galaxies (log[$M_*$] = 7.5-8.5 $M_\odot$) are typically bluer ($\beta_{SED}=-2.0\pm 0.2$ / $\beta_{PL} = -2.1\pm0.4$), fainter ($MUV = -17.8^{+0.7}_{-0.6}$) less dusty ($A{v}=0.4\pm0.1$ mag), exhibit lower rates of star formation (log[SFR]=$0.1\pm0.2 M_\odot/$ yr) and higher specific star formation rates (log[sSFR]=$-8.2\pm0.2 \ \mathrm{yr}^{-1}$) than their high-mass counterparts. Higher-mass galaxies (log[$M_*$] $=10.0-12.0 \ M_\odot$) are on average redder ($\beta_{SED}=-0.9^{+0.8}_{-0.5}$ / $\beta_{PL}=-1.0^{+0.8}_{-0.5}$), brighter ($MUV=-19.6^{+1.0}_{-1.2}$), dustier ($Av = 0.9^{+0.5}_{-0.4}$ mag), have higher SFRs (log[SFR]=$1.2^{+0.6}_{-1.1} M_\odot$ yr), and lower sSFRs (log[sSFR]=$-9.1^{+0.5}_{-1.1} {yr}^{-1}$). This study's substantial sample size provides a benchmark for demonstrating that the rest-frame UV spectral slope correlates with stellar mass-dependent galaxy characteristics at $z\sim2-4$, a relationship less discernible with smaller datasets typically available at higher redshifts.

HESS J1702-420 is a peculiar TeV complex with a morphology changing from a diffuse (HESS J1702-420B source) at $\lesssim 2$ TeV to point-like (HESS J1702-420A) at $\gtrsim 10$ TeV energies. The morphology and the spectral properties of HESS J1702-420 could be understood in terms of a (diffusive) hadronic or leptonic models in which the observed TeV emission arises correpondingly from proton-proton or IC-radiation of relativistic particles present in the region. In this work we perform searches of the X-ray counterpart of HESS J1702-420B source originated from the synchrotron emission of the primary or secondary relativistic electrons produced within leptonic or hadronic models. Such an emission can be largely extent and remain beyond the detection capabilities of a narrow-FoV instruments such as XMM-Newton. We utilise the publicly available first 6-months eROSITA dataset (DR1) fully covering selected for the analysis region of $> 5^\circ$-radius around HESS J1702-420. We discuss biases connected to variable plasma temperature/neutral hydrogen column density in the region and present results based on background modelling approach. The performed analysis does not allow us to detect the extended X-ray counterpart of HESS J1702-420 of $0.07^\circ - 3^\circ$-radii sizes. The derived upper limits are significantly higher than the expected hadronic model flux of the X-ray counterpart. For the leptonic model the derived limits indicate the magnetic field in the region $B\lesssim 2\mu$G. We argue, that the further advances in the diffuse X-ray counterpart searches could be achieved either with next generation missions or Msec-long observational campaigns with currently operating instruments.

To get distances, Type Ia Supernovae magnitudes are corrected for their correlation with lightcurve width and colour. Here we investigate how this standardisation is affected by the SN environment, with the aim to reduce scatter and improve standardisation. We first study the SN Ia stretch distribution, as well as its dependence on environment, as characterised by local and global (g-z) colour and stellar mass. We then look at the standardisation parameter $\alpha$, which accounts for the correlation between residuals and stretch, along with its environment dependence and linearity. We finally compute magnitude offsets between SNe in different astrophysical environments after colour and stretch standardisation, aka steps. This analysis is made possible due to the unprecedented statistics of the ZTF SN Ia DR2 volume-limited sample. The stretch distribution exhibits a bimodal behaviour, as previously found in literature. However, we find the distribution means to decrease with host stellar mass at a 9.0$\sigma$ significance. We demonstrate, at the 14.3$\sigma$ level, that the stretch-magnitude relation is non-linear, challenging the usual linear stretch-residuals relation. Fitting for a broken-$\alpha$ model, we indeed find two different slopes between stretch regimes ($x_1<-0.49\pm0.06$): $\alpha_{low}=0.28\pm0.01$ and $\alpha_{high}=0.09\pm0.01$, a $\Delta_{\alpha}=-0.19\pm0.01$ difference. As the relative proportion of SNe Ia in the high-/low-stretch modes evolves with redshift and environment, this implies that a linear $\alpha$ also evolves with redshift and environment. Concerning the environmental magnitude offset $\gamma$, we find it to be greater than 0.14 mag regardless of the considered environmental tracer used (local or global colour and stellar mass), all measured at the $\geq 6\sigma$ level, increased to $\sim0.18\pm0.01$ mag when accounting for the stretch-non linearity.

We present a detailed investigation of photometric, spectroscopic, and polarimetric observations of the Type II SN 2023ixf. The early detection of highly-ionized flash features, rapid ascent in ultraviolet flux coupled with the blueward shift in near-ultraviolet colors and temperature provides compelling evidence for a delayed shock breakout from a confined dense circumstellar matter (CSM) enveloping the progenitor star. The temporal evolution of polarization in the SN 2023ixf phase revealed three distinct peaks in polarization evolution at 1.4 d, 6.4 d, and 79.2 d, indicating an asymmetric dense CSM, an aspherical shock front and clumpiness in the low-density extended CSM, and an aspherical inner ejecta/He-core. SN 2023ixf displayed two dominant axes, one along the CSM-outer ejecta and the other along the inner ejecta/He-core, showcasing the independent origin of asymmetry in the early and late evolution. The argument for an aspherical shock front is further strengthened by the presence of a high-velocity broad absorption feature in the blue wing of the Balmer features in addition to the P-Cygni absorption post 16 d. Hydrodynamical light curve modeling indicated a progenitor mass of 10 solar mass with a radius of 470 solar radius, explosion energy of 2e51 erg, and 0.06 solar mass of 56Ni. The modeling also indicated a two-zone CSM: a confined dense CSM extending up to 5e14 cm, with a mass-loss rate of 1e-2 solar mass per year, and an extended CSM spanning from 5e14 cm to 1e16 cm with a mass-loss rate of 1e-4 solar mass per year. The early nebular phase observations display an axisymmetric line profile of [OI] and red-ward attenuation of the emission of Halpha post 125 days, marking the onset of dust formation.

The renormalization group for large-scale structure (RG-LSS) describes the evolution of galaxy bias and stochastic parameters as a function of the cutoff $\Lambda$. In this work, we introduce interaction vertices that describe primordial non-Gaussianity into the Wilson-Polchinski framework, thereby extending the free theory to the interacting case. The presence of these interactions forces us to include new operators and bias coefficients to the bias expansion to ensure closure under renormalization. We recover the previously-derived ``scale-dependent bias'' contributions, as well as a new (subdominant) stochastic contribution. We derive the renormalization group equations governing the RG-LSS for a large class of interactions which account for vertices at linear order in $f_{\rm NL}$ that parametrize interacting scalar and massive spinning fields during inflation. Solving the RG equations, we show the evolution of the non-Gaussian contributions to galaxy clustering as a function of scale.

The treatment of convection remains a major weakness in the modelling of stellar evolution with one-dimensional (1D) codes. The ever increasing computing power makes now possible to simulate in 3D part of a star for a fraction of its life, allowing us to study the full complexity of convective zones with hydrodynamics codes. Here, we performed state-of-the-art hydrodynamics simulations of turbulence in a neon-burning convective zone, during the late stage of the life of a massive star. We produced a set of simulations varying the resolution of the computing domain (from 1283 to 10243 cells) and the efficiency of the nuclear reactions (by boosting the energy generation rate from nominal to a factor of 1000). We analysed our results by the mean of Fourier transform of the velocity field, and mean-field decomposition of the various transport equations. Our results are in line with previous studies, showing that the behaviour of the bulk of the convective zone is already well captured at a relatively low resolution (2563), while the details of the convective boundaries require higher resolutions. The different boosting factors used show how various quantities (velocity, buoyancy, abundances, abundance variances) depend on the energy generation rate. We found that for low boosting factors, convective zones are well mixed, validating the approach usually used in 1D stellar evolution codes. However, when nuclear burning and turbulent transport occur on the same timescale, a more sophisticated treatment would be needed. This is typically the case when shell mergers occur.

The Classical T Tauri Star (CTTS) stage is a critical phase of the star and planet formation process. In an effort to better understand the mass accretion process, which can dictate further stellar evolution and planet formation, a multi-epoch, multi-wavelength photometric and spectroscopic monitoring campaign of four CTTSs (TW Hya, RU Lup, BP Tau, and GM Aur) was carried out in 2021 and 2022/2023 as part of the Outflows and Disks Around Young Stars: Synergies for the Exploration of ULYSSES Spectra (ODYSSEUS) program. Here we focus on the HST UV spectra obtained by the HST Director's Discretionary Time UV Legacy Library of Young Stars as Essential Standards (ULLYSES) program. Using accretion shock modeling, we find that all targets exhibit accretion variability, varying from short increases in accretion rate by up to a factor of 3 within 48 hours, to longer decreases in accretion rate by a factor of 2.5 over the course of 1 year. This is despite the generally consistent accretion morphology within each target. Additionally, we test empirical relationships between accretion rate and UV luminosity and find stark differences, showing that these relationships should not be used to estimate the accretion rate for individual target. Our work reinforces that future multi-epoch and simultaneous multi-wavelength studies are critical in our understanding of the accretion process in low-mass star formation.

The mass loss mechanism of red supergiant stars is not well understood, even though it has crucial consequences for their stellar evolution and the appearance of supernovae that occur upon core-collapse. We argue that outgoing shock waves launched near the photosphere can support a dense chromosphere between the star's surface and the dust formation radius at several stellar radii. We derive analytic expressions for the time-averaged density profile of the chromosphere, and we use these to estimate mass loss rates due to winds launched by radiation pressure at the dust formation radius. These mass loss rates are similar to recent observations, possibly explaining the upward kink in mass loss rates of luminous red supergiants. Our models predict that low-mass red supergiants lose less mass than commonly assumed, while high-mass red supergiants lose more. The chromospheric mass of our models is $\sim$0.01 solar masses, most of which lies within a few stellar radii. This may explain the early light curves and spectra of type-II P supernovae without requiring extreme pre-supernova mass loss. We discuss implications for stellar evolution, type II-P supernovae, SN 2023ixf, and Betelgeuse.

With stunning clarity, JWST has revealed the Universe's first billion years. The scientific community is analyzing a wealth of JWST imaging and spectroscopic data from that era, and is in the process of rewriting the astronomy textbooks. Here, 1.5 years into the JWST science mission, we provide a snapshot of the great progress made towards understanding the initial chapters of our cosmic history. We highlight discoveries and breakthroughs, topics and issues that are not yet understood, and questions that will be addressed in the coming years, as JWST continues its revolutionary observations of the Early Universe. While this compendium is written by a small number of authors, invited to ISSI Bern in March 2024 as part of the 2024 ISSI Breakthrough Workshop, we acknowledge the work of a large community that is advancing our collective understanding of the evolution of the Early Universe.

We report results from 13 years of Chandra monitoring of nonthermal X-ray emission from the youngest Galactic supernova remnant G1.9+0.3, the only remnant known to be increasing in brightness. We confirm the spatially-integrated flux increase rate of $(1.2 \pm 0.2)$% yr$^{-1}$ between 1 and 7 keV, but find large spatial variations, from decreases of $-3$% yr$^{-1}$ to increases of 7% yr$^{-1}$, over length scales as small as $10''$ or smaller. We observe relatively little change in spectral slope, though one region shows significant hardening (photon index $\Delta \Gamma \sim 0.4$) as it brightens by 1% yr$^{-1}$. Such rates of change can be accommodated by any of several explanations, including steady evolution of the blast wave, expansion or compression of discrete plasma blobs, strong magnetic turbulence, or variations in magnetic-field aspect angle. Our results do not constrain the mean magnetic-field strength, but a self-consistent picture of the spatially averaged rate of increase can be produced in which the maximum energies of accelerated particles are limited by the remnant age (applying both to electrons and to ions) to about 20 TeV, and the remnant-averaged magnetic field strength is about 30 $\mu$G. The deceleration parameter $m$ (average shock radius varying as $t^m$) is about 0.7, consistent with estimates from overall expansion dynamics, and confirming an explosion date of about 1900 CE. Shock-efficiency factors $\epsilon_e$ and $\epsilon_B$ (fractions of shock energy in relativistic electrons and magnetic field) are 0.003 and 0.0002 in this picture. However, the large range of rates of brightness change indicates that such a global model is oversimplified. Temporal variations of photon index, expected to be small but measurable with longer time baselines, can discriminate among possible models.

Classical T Tauri Stars (CTTSs) are young, low-mass stars which accrete material from their surrounding protoplanetary disk. To better understand accretion variability, we conducted a multi-epoch, multi-wavelength photometric monitoring campaign of four CTTSs: TW Hya, RU Lup, BP Tau, and GM Aur, in 2021 and 2022, contemporaneous with HST UV and optical spectra. We find that all four targets display significant variability in their light curves, generally on days-long timescales (but in some cases year-to-year) often due to periodicity associated with stellar rotation and to stochastic accretion variability. Their is a strong connection between mass accretion and photometric variability in all bands, but the relationship varies per target and epoch. Thus, photometry should be used with caution as a direct measure of accretion in CTTSs.

We derive the central stellar velocity dispersion function for quiescent galaxies in 280 massive clusters with $\log (M_{200} / M_{\odot}) > 14$ in IllustrisTNG300. The velocity dispersion function is an independent tracer of the dark matter mass distribution of subhalos in galaxy clusters. Based on the IllustrisTNG cluster catalog, we select quiescent member subhalos with a specific star formation rate $< 2 \times 10^{-11}$ yr${^-1}$ and stellar mass $\log (M_{*} / M_{\odot}) > 9$. We then simulate fiber spectroscopy to measure the stellar velocity dispersion of the simulated galaxies; we compute the line-of-sight velocity dispersions of star particles within a cylindrical volume that penetrates the core of each subhalo. We construct the velocity dispersion functions for quiescent subhalos within $R_{200}$. The simulated cluster velocity dispersion function exceeds the simulated field velocity dispersion function for $\log \sigma_{*} > 2.2$, indicating the preferential formation of large velocity dispersion galaxies in dense environments. The excess is similar in simulations and in the observations. We also compare the simulated velocity dispersion function for the three most massive clusters with $\log (M_{200} / M_{\odot}) > 15$ with the observed velocity dispersion function for the two most massive clusters in the local universe, Coma and A2029. Intriguingly, the simulated velocity dispersion functions are significantly lower for $\log \sigma_{*} > 2.0$. This discrepancy results from 1) a smaller number of subhalos with $\log (M_{*} / M_{\odot}) > 10$ in TNG300 compared to the observed clusters, and 2) a significant offset between the observed and simulated $M_{*} - \sigma_{*}$ relations. The velocity dispersion function offers a unique window on galaxy and structure formation in simulations.