Papers for Thursday, Jul 28 2022

The quality of astrochemical models is highly dependent on reliable binding energy (BE) values that consider the morphological and energetic variety of binding sites on the surface of ice-grain mantles. Here, we present the Binding Energy Evaluation Platform (BEEP) and database that, using quantum chemical methods, produces full BE distributions of molecules bound to an amorphous solid water (ASW) surface model. BEEP is highly automatized and allows to sample binding sites on set of water clusters and to compute accurate BEs. Using our protocol, we computed 21 BE distributions of interstellar molecules and radicals on an amorphized set of 15-18 water clusters of 22 molecules each. The distributions contain between 225 and 250 unique binding sites. We apply a Gaussian fit and report the mean and standard deviation for each distribution. We compare with existing experimental results and find that the low and high coverage experimental BEs coincide well with the high BE tail and mean value of our distributions, respectively. Previously reported single BE theoretical values are broadly in line with ours, even though in some cases significant differences can be appreciated. We show how the use of different BE values impact a typical problem in astrophysics, such as the computation of snow lines in protoplanetary discs. BEEP will be publicly released so that the database can be expanded to other molecules or ice-models in a community effort.

Identifying active galactic nuclei (AGN) and isolating their contribution to a galaxy's energy budget is crucial for studying the co-evolution of AGN and their host galaxies. Brightness temperature ($T_b$) measurements from high-resolution radio observations at GHz frequencies are widely used to identify AGN. Here we investigate using new sub-arcsecond imaging at 144 MHz with the International LOFAR Telescope to identify AGN using $T_b$ in the Lockman Hole field. We use ancillary data to validate the 940 AGN identifications, finding 83 percent of sources have AGN classifications from SED fitting and/or photometric identifications, yielding 160 new AGN identifications. Considering the multi-wavelength classifications, brightness temperature criteria select over half of radio-excess sources, 32 percent of sources classified as radio-quiet AGN, and 20 percent of sources classified as star-forming galaxies. Infrared colour-colour plots and comparison with what we would expect to detect based on peak brightness in 6 arcsec LOFAR maps, imply that the star-forming galaxies and sources at low flux densities have a mixture of star-formation and AGN activity. We separate the radio emission from star-formation and AGN in unresolved, $T_b$-identified AGN with no significant radio excess and find the AGN comprises $0.49\pm 0.16$ of the radio luminosity. Overall the non-radio excess AGN show evidence for having a variety of different radio emission mechanisms, which can provide different pathways for AGN and galaxy co-evolution. This validation of AGN identification using brightness temperature at low frequencies opens the possibility for securely selecting AGN samples where ancillary data is inadequate.

We have performed a systematic study of the statistical behavior of non-Poissonian template fitting (NPTF), a method designed to analyze and characterize unresolved point sources in general counts datasets. In this paper, we focus on the properties and characteristics of the Fermi-LAT gamma-ray data set. In particular, we have simulated and analyzed gamma-ray sky maps under varying conditions of exposure, angular resolution, pixel size, energy window, event selection, and source brightness. We describe how these conditions affect the sensitivity of NPTF to the presence of point sources, for inner-galaxy studies of point sources within the Galactic Center excess, and for the simplified case of isotropic emission. We do not find opportunities for major gains in sensitivity from varying these choices, within the range available with current Fermi-LAT data. We provide an analytic estimate of the NPTF sensitivity to point sources for the case of isotropic emission and perfect angular resolution, and find good agreement with our numerical results for that case.

We present the Sherwood-Relics simulations, a new suite of large cosmological hydrodynamical simulations aimed at modelling the intergalactic medium (IGM) during and after the cosmic reionization of hydrogen. The suite consists of over 200 simulations that cover a wide range of astrophysical and cosmological parameters. It also includes simulations that use a new lightweight hybrid scheme for treating radiative transfer effects. This scheme follows the spatial variations in the ionizing radiation field, as well as the associated fluctuations in IGM temperature and pressure smoothing. It is computationally much cheaper than full radiation hydrodynamics simulations and circumvents the difficult task of calibrating a galaxy formation model to observational constraints on cosmic reionization. Using this hybrid technique, we study the spatial fluctuations in IGM properties that are seeded by patchy cosmic reionization. We investigate the relevant physical processes and assess their impact on the z > 4 Lyman-alpha forest. Our main findings are: (i) Consistent with previous studies patchy reionization causes large scale temperature fluctuations that persist well after the end of reionization, (ii) these increase the Lyman-alpha forest flux power spectrum on large scales, and (iii) result in a spatially varying pressure smoothing that correlates well with the local reionization redshift. (iv) Structures evaporated or puffed up by photoheating cause notable features in the Lyman-alpha forest, such as flat-bottom or double-dip absorption profiles.

We use Keck/NIRSPEC to survey a sample of of young ($<$1 Gyr), short period mini Neptunes orbiting nearby K dwargs to measure their mass loss via the metastable helium line. We detect helium absorption from all four of the targets in our initial sample. The first detection, around TOI 560b, was announced in a previous paper. We now announce three additional detections around TOI 1430.01, 2076b, and TOI 1683.01. All four planets show an average in-transit excess absorption of 0.7--1.0%. However, the outflows differ in their kinematic properties. TOI 1430b exhibits pre-ingress absorption, while TOI 2076b's outflow is exceptionally optically thick and shows significant post-egress absorption. For all four planets, the width of the measured helium absorption signal is consistent with expectations for a photoevaporative outflow (10--30 km/s, 5000--10,000 K). Unless broadening mechanisms other than thermal velocity and the bulk outflow velocity are significant, our observations disfavor core-powered mass loss models, which predict much slower (1-3 km/s) outflows. We utilize both an isothermal Parker wind model and an order-of-magnitude method to estimate the mass loss timescale, and obtain $\sim$ a few hundred Myr for each planet. We conclude that many, if not all, of these planets will lose their hydrogen-rich envelopes and become super Earths. Our results demonstrate that most mini Neptunes orbiting sun-like stars have primordial atmospheres, and that photoevaporation is an efficient mechanism for stripping these atmospheres and transforming these planets into super Earths.

We have determined bulk proper motions (PMs) for 31 LMC GCs from Gaia eDR3 and Hubble Space Telescope data using multiple independent analysis techniques. Combined with literature values for distances, line-of-sight velocities and existing bulk PMs, we extract full 6D phase-space information for 32 clusters, allowing us to examine the kinematics of the LMC GC system in detail. Except for two GCs (NGC 2159 and NGC 2210) for which high velocities suggest they are not long-term members of the LMC system, the data are consistent with a flattened configuration that rotates like the stellar disk. The one-dimensional velocity dispersions are of order 30 km/s, similar to that of old stellar populations in the LMC disk. Similar to the case for Milky Way disk clusters, the velocity anisotropy is such that the dispersion is smallest in the azimuthal direction; however, alternative anisotropies cannot be ruled out due to distance uncertainties. The data are consistent with a single multi-dimensional Gaussian velocity distribution. Given the non-collisional nature of the LMC disk, this suggests that most, if not all, of the LMC GCs are formed by a single formation mechanism in the stellar disk, despite a significant spread in age and metallicity. Any accreted halo GC population is absent or far smaller in the LMC compared to the Milky Way.

We present the first rest-frame optical size-luminosity relation of galaxies at $z>7$, using the NIRCam imaging data obtained by the GLASS Jame Webb Space Telescope Early Release Science (GLASS-JWST-ERS) program, providing the deepest extragalactic data of the ERS campaign. Our sample consist of 21 photometrically selected bright galaxies with $m_\text{F444W}\leq27.8$ at $7<z<9$ and $m_\text{F444W}<28.2$ at $z\sim9-15$. We measure the size of the galaxies in 5 bands, from the rest-frame optical ($\sim4800\,{\rm \AA}$) to the ultra-violet (UV; $\sim1600\,{\rm \AA}$) based on the S\'ersic model, and analyze the size-luminosity relation as a function of wavelength. Remarkably, the data quality of NIRCam imaging is sufficient to probe the half-light radius $r_e$ down to $\sim 100\,pc$ at $z>7$. Given the limited sample size and magnitude range, we first fix the slope to that observed for larger samples in rest-frame UV using HST samples. The median size $r_0$ at the reference luminosity $M=-21$ decreases slightly from rest-frame optical, $370\pm50\,pc$, to UV, $230\pm50\,pc$. We then re-fit the size-luminosity relation allowing the slope to vary. The slope is consistent with $\beta\sim0.2$ for all bands except F150W, where we find a marginally steeper slope of $\beta=0.51\pm0.13$. The steep UV slope is mainly driven by the smallest and faintest galaxies. If confirmed by larger samples, it implies that the UV size-luminosity relation breaks toward the faint end as suggested by lensing studies.

Patchy cosmic reionization resulted in the ionizing UV background asynchronous rise across the Universe. The latter might have left imprints visible in present day observations. Several numerical simulation-based studies show correlations between reionization time and overdensities and object masses today. To remove the mass from the study, as it may not be the sole important parameter, this paper focuses solely on the properties of paired halos within the same mass range as the Milky Way. For this purpose, it uses CoDaII, a fully-coupled radiation hydrodynamics reionization simulation of the local Universe. This simulation holds a halo pair representing the Local Group, in addition to other pairs, sharing similar mass, mass ratio, distance separation and isolation criteria but in other environments, alongside isolated halos within the same mass range. Investigations of the paired halo reionization histories reveal a wide diversity although always inside-out given our reionization model. Within this model, halos in a close pair tend to be reionized at the same time but being in a pair does not bring to an earlier time their mean reionization. The only significant trend is found between the total energy at z = 0 of the pairs and their mean reionization time: pairs with the smallest total energy (bound) are reionized up to 50 Myr earlier than others (unbound). Above all, this study reveals the variety of reionization histories undergone by halo pairs similar to the Local Group, that of the Local Group being far from an average one. In our model, its reionization time is ~625 Myr against 660+/-4 Myr (z~8.25 against 7.87+/-0.02) on average.

Ever since the launch of the NuSTAR mission, the hard X-ray range is being covered to an unprecedented sensitivity. This range encodes the reflection features arising from active galactic nuclei (AGN). Especially, the reflection of the primary radiation off the accretion disk carries the features of the manifestation of General Relativity described by the Kerr metric due to rotating supermassive black holes (SMBHs). We show the results of the broadband analyses of Mrk 876. The spectra exhibit the signature of a Compton hump at energies above 10 keV and a broadened and skewed excess at energies ~6 keV. We establish this spectral excess to be statistically significant at 99.71% (~3sigma) that is the post-trail probability through Monte Carlo simulations. Based on the spectral fit results and the significance of spectral features the relativistic reflection model is favored over the distant reflection scenario. The excess at ~6 keV has a complex shape that we try to recover along with the Compton hump through a self-consistent X-ray reflection model. This allows inferring an upper limit to the black hole spin of a<0.85, while the inclination angle of the accretion disk results in i=32.84 (+12.22\-8.99) degrees, which is in agreement within the errors with a previous independent measurement (i=15.4 (+12.1/-6.8) degrees). While most spin measurements are biased toward high spin values, the black hole mass of Mrk~876 (2.4x10^8 Msun < M < 1.3x10^9 Msun) lies in a range where moderately spinning SMBHs are expected. Moreover, the analyses of twelve Chandra observations reveal for the first time X-ray variability of Mrk876 with an amplitude of 40%.

One of the primary challenges facing upcoming CMB polarization experiments aiming to measure the inflationary B-mode signal is the removal of polarized foregrounds. The thermal dust foreground is often modeled as a single modified blackbody, however overly simplistic foreground models can bias measurements of the tensor-to-scalar ratio r. As CMB polarization experiments become increasingly sensitive, thermal dust emission models must account for greater complexity in the dust foreground while making minimal assumptions about the underlying distribution of dust properties within a beam. We use Planck dust temperature data to estimate the typical variation in dust properties along the line of sight and examine the impact of these variations on the bias in r if a single modified blackbody model is assumed. We then assess the ability of the moment method to capture the effects of spatial averaging and to reduce bias in the tensor-to-scalar ratio for different possible toy models of dust emission. We find that the expected bias due to temperature variations along the line of sight is significant compared to the target sensitivities of future CMB experiments, and that the use of the moment method could reduce bias as well as shed light into the distribution of dust physical parameters.

Calcium-rich (SN 2005E-like) explosions are very faint (typical -15.5, type I supernovae (SNe) showing strong Ca-lines, mostly observed in old stellar environments. Several models for such SNe had been explored and debated, but non were able to consistently reproduce the observed properties of Ca-rich SNe, nor their rates and host-galaxy distributions. Here we show that the disruptions of low-mass carbon-oxygen (CO) white dwarfs (WDs) by hybrid Helium-CO (HeCO) WDs during their merger could explain the origin and properties of such SNe. We make use of detailed multi-dimensional hydrodynamical-thermonuclear (FLASH) simulations models of HeCO-CO WD-WD mergers to characterize such explosions. We find that the accretion of CO material onto a HeCO-WD heats its He-shell and eventually leads to its "weak" detonation and ejection and the production of a sub-energetic $\sim10^{49}$erg Ca-rich SN, while leaving the CO core of the HeCO-WD intact as a hot remnant white dwarf, possibly giving rise to X-ray emission as it cools down. We model the detailed light-curve and spectra of such explosions to find an excellent agreement with observations of Ia/c Ca-rich SNe. In these models most of the He is burned, producing Ia/c Ca-rich SNe but higher He-enrichment levels might potentially also explain the Ib Ca-rich SNe. We thereby provide a viable, consistent model for the origins of Ca-rich. These findings can shed new light on the role of Ca-rich in the chemical evolution of galaxies and the intra-cluster medium, and their contribution to the observed 511 kev signal in the Galaxy originating from positrons produced from ${}^{44}\textrm{Ti}$ decay. Finally, the origins of such SNe points to the key-role of HeCO WDs as SN progenitors, and their potential role as progenitors of other thermonuclear SNe including normal type Ia SNe.

We devise a method to constrain self-interacting dark matter (SIDM) from observations of quadruply-imaged quasars, and apply it to five self-interaction potentials with a long-range dark force. We consider several SIDM models with an attractive potential that allows for the formation of quasi-bound states, giving rise to resonant features in the cross section localized at particular velocities below $50 \ \rm{km} \ \rm{s^{-1}}$. We propose these resonances, which amplify or suppress the cross section amplitude by over an order of magnitude, accelerate or delay the onset of core collapse in low-mass dark matter halos, and derive constraints on the timescale for core collapse for the five interaction potentials we consider. Our data strongly disfavors scenarios in which a majority of halos core collapse, with the strongest constraints obtained for cross section strengths exceeding $100 \ \rm{cm^2} \rm{g}^{-1}$ at relative velocities below $30 \ \rm{km} \ \rm{s^{-1}}$. This work opens a new avenue to explore the vast landscape of possible SIDM theories.

We report the first gas-phase metallicity map of a distant galaxy measured with the James Webb Space Telescope (JWST). We use the NIRISS slitless spectroscopy acquired by the GLASS Early Release Science program to spatially resolve the rest-frame optical nebular emission lines in a gravitationally lensed galaxy at $z=3.06$ behind the Abell 2744 galaxy cluster. This galaxy (dubbed GLASS-Zgrad1) has stellar mass $\sim10^9 M_\odot$, instantaneous star formation rate $\sim30$ $M_\odot$/yr (both corrected for lensing magnification), and global metallicity one-third solar. From its emission line maps ([O III], H$\beta$, H$\gamma$, [O II]) we derive its spatial distribution of gas-phase metallicity using a well-established forward-modeling Bayesian inference method. The exquisite resolution and sensitivity of JWST/NIRISS, combined with lensing magnification, enable us to resolve this $z\sim3$ dwarf galaxy in $\gtrsim$50 resolution elements with sufficient signal, an analysis hitherto not possible. We find that the radial metallicity gradient of GLASS-Zgrad1 is strongly inverted (i.e. positive): $\Delta\log({\rm O/H})/\Delta r$ = $0.165\pm0.023$ $\mathrm{dex~kpc^{-1}}$, in stark contrast to recent cosmological hydrodynamic simulations which generally predict flat or negative gradients at these redshifts. These first results showcase the power of JWST wide-field slitless spectroscopic modes to resolve the mass assembly and chemical enrichment of low-mass galaxies in and beyond the peak epoch of cosmic star formation ($z\gtrsim2$). Reaching masses $\lesssim 10^9~M_\odot$ at these redshifts is especially valuable to constrain the effects of galactic feedback, and is possible only with JWST's new capabilities.

The photon counting imaging paradigm in the visible and the infrared comes from the very small energy carried by a single photon at these wavelengths. Usually to detect photons the photoelectric effect is used. It converts a photon to a single electron making it very difficult to detect because of the readout noise of the electronics. To overcome this there are two strategies, either to amplify the signal to make it larger than the readout noise (used in the so called gain or amplified detectors), or to lower the readout noise in a standard image sensor. For a long time, only amplified detectors were able to do some photon counting. Since the first photon counting systems in the visible, developed by Boksenberg and his collaborators in 1972, many groups around the world improved photon counting techniques. In the 2000's in the visible, EMCCDs (electron multiplying charge coupled devices) allowed to replace the classical image intensifier photon counting systems by solid state devices and improved a lot the QE. But EMCCDs suffer from several issues, and the most important of them is the excess noise factor which prevents to know the exact incoming number of photons in the case of multiple photons per pixel. In the infrared there was no equivalent to EMCCDs up to the development of e-APD sensors and cameras made with HgCdTe material (electron initiated avalanche photo diode). With an excess noise factor near 1 at low temperatures, photon counting is possible with these devices but only in the infrared. We will show that having excess noise factor prevents from being able to do multiple photon counting (quanta imaging) and the only solution is to lower the readout noise.

We present new Big Bang Nucleosynthesis (BBN) limits on the cosmic expansion rate or relativistic energy density, quantified via the number $N_\nu$ of equivalent neutrino species. We use the latest light element observations, neutron mean lifetime, and update our evaluation for the nuclear rates $d+d \rightarrow He3 + n$ and $d+d \rightarrow H3 + p$. Combining this result with the independent constraints from the cosmic microwave background (CMB) yields tight limits on new physics that perturbs $N_\nu$ and $\eta$ prior to cosmic nucleosynthesis: a joint BBN+CMB analysis gives $N_\nu = 2.898 \pm 0.141$, resulting in $N_\nu < 3.180$ at $2\sigma$. We apply these limits to a wide variety of new physics scenarios including right-handed neutrinos, dark radiation, and a stochastic gravitational wave background. We also search for limits on potential {\em changes} in $N_\nu$ and/or the baryon-to-photon ratio $\eta$ between the two epochs. The present data place strong constraints on the allowed changes in $N_\nu$ between BBN and CMB decoupling; for example, we find $-0.708 < N_\nu^{\rm CMB}-N_\nu^{\rm BBN} < 0.328$ in the case where $\eta$ and the primordial helium mass fraction $Y_p$ are unchanged between the two epochs; we also give limits on the allowed variations in $\eta$ or in $(\eta,N_\nu)$ jointly. Looking to the future, we forecast the tightened precision for $N_\nu$ arising from both CMB Stage 4 measurements as well as improvements in astronomical \he4 measurements. We find that CMB-S4 combined with present BBN and light element observation precision can give $\sigma(N_\nu) \simeq 0.03$. Such future precision would reveal the expected effect of neutrino heating ($N_{\rm eff}-3=0.044$) of the CMB during BBN, and would be near the level to reveal any particle species ever in thermal equilibrium with the standard model.

We study 147 star clusters in the Large Magellanic Cloud (LMC) in order to determine their mean metallicities and ages, as well as the mean metallicities of 80 surrounding fields. We construct an age-metallicity relation (AMR) for the clusters in the LMC. For this purpose, we used Str\"omgren photometry obtained with the SOI camera on the 4.1 m SOAR telescope. We found the mean metallicity and age for 110 star clusters. For the remaining 37, we provide an age estimation only. To the best of our knowledge, for 29 clusters from our sample, we provide both the metallicity and age for the first time, whereas for 66 clusters, we provide a first determination of the metallicity, and for 43 clusters, the first estimation of the age. We also calculated the mean metallicities for stars from 80 fields around the clusters. The old, metal-poor star clusters occur both in and out of the LMC bar region, while intermediate-age clusters are located mostly outside of the bar. The majority of star clusters younger than 1 Gyr are located in the bar region. We find a good agreement between our AMR and theoretical models of the LMC chemical enrichment, as well as with AMRs for clusters from the literature. Next, we took advantage of 26 stellar clusters from our sample which host Cepheid variables and used them as an independent check of the correctness of our age determination procedure. We used period-age relations for Cepheids to calculate the mean age of a given cluster and compared it with the age obtained from isochrone fitting. We find good agreement between these ages, especially for models taking into account additional physical processes (e.g., rotation). We also compared the AMR of the LMC and Small Magellanic Cloud (SMC) derived in a uniform way and we note that they indicate possible former interaction between these two galaxies.

[Abridged] Planets and their atmospheres are built from gas and solid material in protoplanetary disks. This solid material grows from smaller, micron-sized grains to larger sizes in the disks, during the process of planet formation. Our goal is to model the compositional evolution of volatile ices on grains of different sizes, assuming both time-dependent grain growth and several constant grain sizes. The state-of-the-art Walsh chemical kinetics code is utilised for modeling chemical evolution. This code has been upgraded to account for the time-evolving sizes of solids. Chemical evolution is modelled locally at four different radii in a protoplanetary disk midplane for up to 10Myr. The evolution is modelled for five different constant grain sizes, and one model where the grain size changes with time according to a grain growth model appropriate for the disk midplane. Local grain growth, with conservation of total grain mass and the assumption of spherical grains, acts to reduced the total grain-surface area that is available for ice-phase reactions. This reduces these reactions efficiency compared to a chemical scenario with a conventional grain-size choice of 0.1$\mu$m. The modelled chemical evolution with grain growth leads to increased abundances of H$_{2}$O ice. For carbon in the inner disk, grain growth leads CO gas to overtake CO$_{2}$ ice as dominant carrier, and in the outer disk, CH$_{4}$ ice to become the dominant carrier. Overall, a constant grain size adopted from a grain evolution model leads to almost identical chemical evolution, when compared with chemical evolution with evolving grain sizes. A constant grain size choice, albeit larger than 0.1$\mu$m, may therefore be an appropriate simplification when approximating the impact of grain growth on chemical evolution.

The evolution of the dust grain size distribution has been studied in recent years with great detail in cosmological hydrodynamical simulations taking into account all the channels under which dust evolves in the interstellar medium. We present a systematic analysis of the observed spectral energy distribution of a large sample of galaxies in the local universe in order to derive not only the total dust masses but also the relative mass fraction between small and large dust grains (DS/DL). Simulations reproduce fairly well the observations except for the high stellar mass regime where dust masses tend to be overestimated. We find that ~45% of galaxies exhibit DS/DL consistent with the expectations of simulations, while there is a sub-sample of massive galaxies presenting high DS/DL (log(DS/DL)~-0.5), and deviating from the prediction in simulations. For these galaxies, which also have high molecular gas mass fractions and metallicities, coagulation is not an important mechanism affecting the dust evolution. Including diffusion, transporting large grains from dense regions to a more diffuse medium where they can be easily shattered, would explain the observed high DS/DL values in these galaxies. With this study we reinforce the use of the small-to-large grain mass ratio to study the relative importance of the different mechanisms in the dust life cycle. Multi-phase hydrodynamical simulations with detailed feedback prescriptions and more realistic subgrid models for the dense phase could help to reproduce the evolution of the dust grain size distribution traced by observations.

CMB-S4 will map the cosmic microwave background to unprecedented precision, while simultaneously surveying the millimeter-wave time-domain sky, in order to advance our understanding of cosmology and the universe. CMB-S4 will observe from two sites, the South Pole and the Atacama Desert of Chile. A combination of small- and large-aperture telescopes with hundreds of thousands of polarization-sensitive detectors will observe in several frequency bands from 20-300 GHz, surveying more than 50 percent of the sky to arcminute resolution with unprecedented sensitivity. CMB-S4 seeks to make a dramatic leap in sensitivity while observing across a broad range of largely unprotected spectrum which is increasingly being utilized for terrestrial and satellite transmissions. Fundamental aspects of CMB instrument technology leave them vulnerable to radio frequency interference (RFI) across a wide range of frequencies, including frequencies outside of their observing bands. Ground-based CMB instruments achieve their extraordinary sensitivities by deploying large focal planes of superconducting bolometers to extremely dry, high-altitude sites, with large fractional bandwidths, wide fields of view, and years of integration time. Suitable observing sites have historically offered significant protection from RFI, both naturally through their extremely remote locations as well as through restrictions on local emissions. Since the coupling mechanisms are complex, safe levels or frequencies of emission that would not interfere with CMB measurements cannot always be determined through straightforward calculations. We discuss models of interference for various types of RFI relevant to CMB-S4, mitigation strategies, and the potential impacts on survey sensitivity.

Planetary engulfment events can occur while host stars are on the main sequence. The addition of rocky planetary material during engulfment will lead to refractory abundance enhancements in the host star photosphere, but the level of enrichment and its duration will depend on mixing processes that occur within the stellar interior, such as convection, diffusion, and thermohaline mixing. We examine engulfment signatures by modeling the evolution of photospheric lithium abundances. Because lithium can be burned before or after the engulfment event, it produces unique signatures that vary with time and host star type. Using MESA stellar models, we quantify the strength and duration of these signatures following the engulfment of a 1, 10, or 100 $M_{\oplus}$ planetary companion with bulk Earth composition, for solar-metallicity host stars with masses ranging from 0.5$-$1.4 $M_{\odot}$. We find that lithium is quickly depleted via burning in low-mass host stars ($\lesssim 0.7 \, M_\odot$) on a time scale of a few hundred Myrs, but significant lithium enrichment signatures can last for Gyrs in G-type stars ($\sim \! 0.9 \, M_{\odot}$). For more massive stars (1.3$-$1.4 $M_{\odot}$), engulfment can enhance internal mixing and diffusion processes, potentially decreasing the surface lithium abundance. Our predicted signatures from exoplanet engulfment are consistent with observed lithium-rich solar-type stars and abundance enhancements in chemically inhomogeneous binary stars.

Repeating and apparently non-repeating Fast Radio Bursts are distinct classes of events produced by distinct classes of sources. As for conference proceedings, this paper reviews the evidence for that division, and then discusses the statistics of and possible models of each class of source: black hole accretion discs for repeating FRB and hypermagnetized neutron stars (SGR) for apparently non-repeating FRB.

We hereby report a low-speed (about~21~km$\cdot$~s$^{-1}$ with respect to the Sun) intruder member in the Hyades cluster based on the data in the literature. The results show that the star is a non-native member star for the Hyades, with its radial velocity being smaller than the radial velocity of the Hyades cluster, even exceeding the standard deviation of the radial velocity of the cluster by a factor of 9. Furthermore, by analyzing and comparing the orbits of this star and its host, it may have intruded into its host in the past 2~Myr. If the star's current motion orbit remains unchanged, it may leave its host in the next 2~Myr. This implies that the intruder star may be temporarily residing in the cluster. This study presents the first observational evidence of a star intrusion into a cluster, which suggests that more evidence may be found.

We present a study of scintillation induced by the mid-latitude ionosphere. By implementing methods currently used in Interplanetary Scintillation studies to measure amplitude scintillation at low frequencies, we have proven it is possible to use the Murchison Widefield Array to study ionospheric scintillation in the weak regime, which is sensitive to structures on scales $\sim$300 m at our observing frequency of 154 MHz, where the phase variance on this scale was 0.06 rad$^{2}$ in the most extreme case observed. Analysing over 1000 individual 2-minute observations, we compared the ionospheric phase variance with that inferred with previous measurements of refractive shifts, which are most sensitive to scales almost an order of magnitude larger. The two measurements were found to be highly correlated (Pearson correlation coefficient 0.71). We observed that for an active ionosphere, the relationship between these two metrics is in line with what would be expected if the ionosphere's structure is described by Kolmogorov turbulence between the relevant scales of 300m and 2000m. In the most extreme ionospheric conditions, the refractive shifts were sometimes found to underestimate the small-scale variance by a factor of four or more, and it is these ionospheric conditions that could have significant effects on radio astronomy observations.

Warps are vertical distortions of the stellar or gaseous disks of galaxies. One of the proposed scenarios for the formation of warps involves tidal interactions among galaxies. A recent study identified a stellar warp in the outer regions of the south-western (SW) disk of the Large Magellanic Cloud (LMC) and suggested that it might have originated due to the tidal interaction between the LMC and the Small Magellanic Cloud (SMC). Due to the limited spatial coverage of the data, the authors could not investigate the counterpart of this warp in the north-eastern (NE) region, which is essential to understanding the global shape, nature, and origin of the outer LMC warp. In this work, we study the structure of the LMC disk using data on red clump stars from the Gaia Early Data Release 3 (EDR3), which cover the entire Magellanic system. We detected a warp in the NE outer LMC disk which is deviated from the disk plane in the same direction as that of the SW outer warp, but with a lower amplitude. This suggests that the outer LMC disk has an asymmetric stellar warp, which is likely to be a U-shaped warp. Our result provides an observational constraint to the theoretical models of the Magellanic system aimed at improving the understanding the LMC-SMC interaction history.

On June 14, 2022, the EHT collaboration (hereafter EHTC) made the web page (https://eventhorizontelescope.org/blog/imaging-reanalyses-eht-data) with the title "Imaging Reanalyses of EHT Data," in which they made comments on our recent Miyoshi et al .2022 published in the Astrophysical Journal. We investigated the EHTC comments and found that all of the five points raised by the EHTC are subjective and unsubstantiated claims. Thus they do not prove the correctness of the result of EHTC. Sincerely we hope that the EHTC will publish, not a collection of unsubstantiated claims, but a discussion based on scientific arguments.

Comets are a rich reservoir of complex organic molecules, such as alcohols, aldehydes, acids, carbon chains, diols, the N-bearing organics formamide, acetonitrile, cyanoacetylene, imines and the simple amino acid glycine. It has not yet been proven that all of the molecules identified in the coma of comets originate from the frozen ices inside the cometary nucleus. Comets showing moderate to high activity can reach sufficient coma densities for molecules to form by active gas-phase coma chemistry. Using a multifluid chemical-hydrodynamical model and an updated chemical network, we studied the effect of gas-phase chemistry that occurs in the cometary coma of four Oort cloud comets, namely C/1996 B2 (Hyakutake), C/2012 F6 (Lemmon), C/2013 R1 (Lovejoy), and C/2014 Q2 (Lovejoy). These comets are significantly enriched in complex organics and show moderate to high activity. We studied the formation of a large number of CHO molecules, N-bearing molecules and the simplest amino acid glycine. We found that by incorporating new pathways, the production rates for HCOOH, HCOOCH$_3$, CH$_3$CHO, CH$_3$CN, CH$_3$COOH, HC$_3$N, HC$_5$N, and NH$_2$CHO can be increased, which can account at least partially towards the total production rates. However, the production rates for C$_2$H$_5$OH, (CH$_2$OH)$_2$ , CH$_2$OHCHO and the simplest amino acid glycine are low for most comets; therefore, their formation requires surface chemistry almost solely. We also found that factors such as initial cometary abundance, relative abundances of the reactants and temperature of the reacting species significantly affect the formation of molecular species in different regions of the coma.

The two most common components of several upcoming CMB experiments are large arrays of superconductive TES (Transition-Edge Sensor) detectors and polarization modulator units, e.g. continuously-rotating Half-Wave Plates (HWP). A high detector count is necessary to increase the instrument raw sensitivity, however past experiments have shown that systematic effects are becoming one of the main limiting factors to reach the sensitivity required to detect primordial $B$-modes. Therefore, polarization modulators have become popular in recent years to mitigate several systematic effects. Polarization modulators based on HWP technologies require a rotating mechanism to spin the plate and modulate the incoming polarized signal. In order to minimize heat dissipation from the rotating mechanism, which is a stringent requirement particularly for a space mission like $LiteBIRD$, we can employ a superconductive magnetic bearing to levitate the rotor and achieve contactless rotation. A disadvantage of this technique is the associated magnetic fields generated by those systems. In this paper we investigate the effects on a TES detector prototype and find no detectable $T_c$ variations due to an applied constant (DC) magnetic field, and a non-zero TES response to varying (AC) magnetic fields. We quantify a worst-case TES responsivity to the applied AC magnetic field of $\sim10^5$ pA/G, and give a preliminary interpretation of the pick-up mechanism.

FU Ori type objects (FUors) are decades-long outbursts of accretion onto young stars that are strong enough to viscously heat disks so that the disk outshines the central star. We construct models for FUor objects by calculating emission components from a steady-state viscous accretion disk, a passively-heated dusty disk, magnetospheric accretion columns, and the stellar photosphere. We explore the parameter space of the accretion rate and stellar mass to investigate implications on the optical and near-infrared spectral energy distribution and spectral lines. The models are validated by fitting to multi-wavelength photometry of three confirmed FUor objects, FU Ori, V883 Ori and HBC 722 and then comparing the predicted spectrum to observed optical and infrared spectra. The brightness ratio between the viscous disk and the stellar photosphere, $\eta$, provides an important guide for identifying viscous accretion disks, with $\eta=1$ ("transition line") and $\eta=5$ ("sufficient dominance line") marking turning points in diagnostics, evaluated here in the near-infrared. These turning points indicate the emergence and complete development of FUor-characteristic strong CO absorption, weak metallic absorption, the triangular spectral continuum shape in the $H$-band, and location in color-magnitude diagrams. Lower stellar mass $M_*$ and higher accretion rate $\dot{M}$ lead to larger $\eta$; for $M_*=0.3~{\rm M_\odot}$, $\eta=1$ corresponds to $\dot{M}=2\times10^{-7}~{\rm M_\odot}/$yr and $\eta=5$ to $\dot{M}=6\times10^{-7}~{\rm M_\odot}/$yr. The sufficient dominance line also coincides with the expected accretion rate where accreting material directly reaches the star. We discuss implications of the models on extinction diagnostics, FUor brightening timescales, viscous disks during initial protostellar growth, and eruptive young stellar objects (YSO) associated with FUors.

The presence of an early-formed giant planet in the protoplanetary disk has mixed influence on the growth of other planetary embryos. Gravitational perturbation from the planet can increase the relative velocities of planetesimals at the mean motion resonances to very high values and impede accretion at those locations. However, gas drag can also align the orbital pericenters of equal-size planetesimals in certain disk locations and make them dynamically quiet and "accretion-friendly" locations for planetesimals of similar sizes. Following the previous paper, where we investigated the effect of a Jupiter-like planet on an external planetesimal disk, we generalize our findings to extrasolar planetary systems by varying the planet parameters. In particular, we focus on the dependence of the planetesimal relative velocities on the mass and eccentricity of the existing planet. We found that the velocity dispersion of identical-mass particles increases monotonically with increasing planet mass. Meanwhile, the dependence of the relative velocity between different-mass planetesimals on their mass ratio becomes weaker as the planet mass increases. While the relative velocities generally increases with increasing planet eccentricity, the velocity dispersion of smaller-mass particles ($m \lesssim 10^{18}~\rm{g}$) is almost independent of planet eccentricity owing to their strong coupling to gas. We find that the erosion limits are met for a wider range of parameters (planet mass/eccentricity, planetesimal mass ratio) when the planetesimal size decreases. Our results could provide some clues for the formation of Saturn's core as well as the architecture of some exoplanetary systems with multiple cold giant planets.

We analyze the cooling and feedback properties of 48 galaxy clusters at redshifts $0.4 < z < 1.3$ selected from the South Pole Telescope (SPT) catalogs to evolve like the progenitors of massive and well-studied systems at $z{\sim}0$. We estimate the radio power at the brightest cluster galaxy (BCG) location of each cluster from an analysis of Australia Telescope Compact Array (ATCA) data. Assuming that the scaling relation between radio power and active galactic nucleus (AGN) cavity power $P_{\mathrm{cav}}$ observed at low redshift does not evolve with redshift, we use these measurements in order to estimate the expected AGN cavity power in the core of each system. We estimate the X-ray luminosity within the cooling radius $L_{\mathrm{cool}}$ of each cluster from a joint analysis of the available $Chandra$ X-ray and SPT Sunyaev-Zel'dovich (SZ) data. This allows us to characterize the redshift evolution of the $P_{\mathrm{cav}} / L_{\mathrm{cool}}$ ratio. When combined with low-redshift results, these constraints enable investigations of the properties of the feedback/cooling cycle across 9~Gyr of cluster growth. We model the redshift evolution of this ratio measured for cool core clusters by a log-normal distribution $\mathrm{Log}$-$\mathcal{N}(\alpha + \beta z, \sigma^2)$ and constrain the slope of the mean evolution $\beta = -0.05\pm 0.47$. This analysis improves the constraints on the slope of this relation by a factor of two. We find no evidence of redshift evolution of the feedback/cooling equilibrium in these clusters which suggests that the onset of radio-mode feedback took place at an early stage of cluster formation. High values of $P_{\mathrm{cav}} / L_{\mathrm{cool}}$ are found at the BCG location of non-cool core clusters which might suggest that the timescales of the AGN feedback cycle and the cool core / non-cool core transition are different.

We present a multiwavelength study of RXC J2014.8-2430, the most extreme cool-core cluster in the Representative $XMM-Newton$ Cluster Structure Survey (REXCESS), using $Chandra$ X-ray, Southern Astrophysical Research (SOAR) Telescope, Atacama Large Millimeter/submillimeter Array (ALMA), Very Large Array (VLA), and Giant Metrewave Radio Telescope (GMRT) observations. While feedback from an active galactic nucleus (AGN) is thought to be the dominant mechanism by which a cooling flow is suppressed, the $Chandra$ imaging observations surprisingly do not reveal the bi-lateral X-ray cavities expected in the intracluster medium (ICM) of an extreme cool core hosting a powerful radio source. We discuss the limits on the presence of any radio bubbles associated with any undetected X-ray cavities. We place upper limits on any significant X-ray AGN in the brightest cluster galaxy, and show that the X-ray peak is offset from the central radio source, which exhibits a steep low frequency radio spectrum indicative of electron ageing. The SOAR data reveal an extended, luminous emission line source. From our narrowband H$\alpha$ imaging of the BCG, the central H$\alpha$ peak is coincident with the radio observations, yet offset from the X-ray peak, consistent with sloshing found previously in this cluster. ALMA observations reveal a large reservoir of molecular gas that traces the extended H$\alpha$ emission. We conclude either that the radio source and its cavities in the X-ray gas are nearly aligned along the line of sight, or that ram pressure induced by sloshing has significantly displaced the cool molecular gas feeding it, perhaps preempting the AGN feedback cycle. We argue that the sloshing near the core is likely subsonic, as expected, given the co-location of the H$\alpha$, CO(1-0), radio continuum, and stellar emission peaks and their proximity to the intact cool core seen in X-ray.

Identifying the young optically visible population in a star-forming region is essential for fully understanding the star formation event. In this paper, We identify 211 candidate members of the Perseus molecular cloud based on Gaia astronomy. We use LAMOST spectra to confirm that 51 of these candidates are new members, bringing the total census of known members to 856. The newly confirmed members are less extincted than previously known members. Two new stellar aggregates are identified in our updated census. With the updated member list, we obtain a statistically significant distance gradient of $\rm 4.84\;pc\;deg^{-1}$ from west to east. Distances and extinction corrected color-magnitude diagrams indicate that NGC 1333 is significantly younger than IC 348 and the remaining cloud regions. The disk fraction in NGC 1333 is higher than elsewhere, consistent with its youngest age. The star formation scenario in the Perseus molecular cloud is investigated and the bulk motion of the distributed population is consistent with the cloud being swept away by the Per-Tau Shell.

We present the first results from our ongoing project to study extremely low mass (ELM) white dwarfs (WDs) ($M$ $\leq$ 0.3$M_{\sun}$) with the Large Sky Area Multi-Object Fibre Spectroscopic Telescope (LAMOST) spectra. Based on the LAMOST DR8 spectral database, we analyzed 136 ELM WD candidates selected from $\it Gaia$ DR2 data and 12 known objects previously identified by the ELM Survey. The atmospheric parameters and radial velocities of these stars were obtained by fitting the LAMOST low-resolution spectra. After comparing the atmospheric parameters of the 12 known objects from this work to the results reported by the ELM Survey, we demonstrated the potential of LAMOST spectra in probing into the nature of ELM WDs. Based on the atmospheric parameters and $\it Gaia$ EDR3 data, we identified 21 new high-probability ELM WDs with masses $M$ $\leq$ 0.3$M_{\sun}$ and parallax estimates that agree to within a factor of 3. Two of them, J0338+4134 and J1129+4715, show significant radial velocity variability and are very likely to be binary systems containing at least one ELM WD.

More than 20 years after the end of NASA's Compton Gamma-Ray Observatory mission, the data collected by its Imaging Compton Telescope (COMPTEL) still provide the most comprehensive and deepest view of our Universe in MeV gamma rays. While most of the COMPTEL data are archived at NASA's High Energy Astrophysics Science Archive Research Center (HEASARC), the absence of any publicly available software for their analysis means the data cannot benefit from the scientific advances made in the field of gamma-ray astronomy at higher energies. To make this unique treasure again accessible for science, we developed open source software that enables a comprehensive and modern analysis of the archived COMPTEL telescope data. Our software is based on a dedicated plug-in to the GammaLib library, a community-developed toolbox for the analysis of astronomical gamma-ray data. We implemented high-level scripts for building science analysis workflows in ctools, a community-developed gamma-ray astronomy science analysis software framework. We describe the implementation of our software and provide the underlying algorithms. Using data from the HEASARC archive, we demonstrate that our software reproduces derived data products that were obtained in the past using the proprietary COMPTEL software. We furthermore demonstrate that our software reproduces COMPTEL science results published in the literature. This brings the COMPTEL telescope data back into life, allowing them to benefit from recent advances in gamma-ray astronomy, and gives the community a means to unveil its still hidden treasures.

Close pairs of galaxies have been broadly studied in the literature as a way to understand galaxy interactions and mergers. In observations they are usually defined by setting a maximum separation in the sky and in velocity along the line of sight, and finding galaxies within these ranges. However, when observing the sky, projection effects can affect the results, by creating spurious pairs that are not close in physical distance. In this work we mimic these observational techniques to find pairs in The Three Hundred simulations of clusters of galaxies. The galaxies' 3D coordinates are projected into 2D, with Hubble flow included for their line-of-sight velocities. The pairs found are classified into "good" or "bad" depending on whether their 3D separations are within the 2D spatial limit or not. We find that the fraction of good pairs can be between 30 and 60 per cent depending on the thresholds used in observations. Studying the ratios of observable properties between the pair member galaxies, we find that the likelihood of a pair being "good" can be increased by around 40, 20 and 30 per cent if the given pair has, respectively, a mass ratio below 0.2, metallicity ratio above 0.8, or colour ratio below 0.8. Moreover, shape and stellar-to-halo mass ratios respectively below 0.4 and 0.2 can increase the likelihood by 50 to 100 per cent. These results suggest that these properties can be used to increase the chance of finding good pairs in observations of galaxy clusters and their environment.

Wide Field Slitless Spectroscopy (WFSS) provides a powerful tool for detecting strong line emission in star forming galaxies (SFGs) without the need for target pre-selection. As part of the GLASS-JWST-ERS program, we leverage the near-infrared wavelength capabilities of NIRISS ($1-2.2\mu$m) to observe rest-optical emission lines out to $\rm{z}\sim 3.4$, to a depth and with a spatial resolution higher than ever before (H$\alpha$ to z<2.4; [OIII]+H$\beta$ to z<3.4). In this letter we constrain the rest-frame [OIII]$\lambda5007$ equivalent width (EW) distribution for a sample of 76 $1<\rm{z}<3.4$ SFGs in the Abell 2744 Hubble Frontier Field and determine an abundance fraction of extreme emission line galaxies with EW$>750$A in our sample to be $12\%$. We determine a strong correlation between the measured H$\beta$ and [OIII]$\lambda5007$ EWs, supporting that the high [OIII]$\lambda5007$ EW objects require massive stars in young stellar populations to generate the high energy photons needed to doubly ionise oxygen. We extracted spectra for objects up to 2 mag fainter in the near-infrared than previous WFSS studies with the Hubble Space Telescope. Thus, this work clearly highlights the potential of JWST/NIRISS to provide high quality WFSS datasets in crowded cluster environments.

We report the discovery and characterization of high-speed (>100 km/s) horizontal flows in solar active regions, making use of the Sun-as-a-star spectroscopy in the range 5-105 nm provided by the EVE (Extreme Ultraviolet Variability Experiment) spectrometers on the Solar Dynamics Observatory. These apparent flows are persistent on time scales of days, and are well observed in lines of Mg X, Si XII and Fe XVI for example. They are prograde, as evidenced directly by blueshifts/redshifts peaking at the east/west limb passages of isolated active regions. The high-speed flow behavior does not depend upon active-region latitude or solar cycle, with similar behavior in Cycles 24 and 25.

In the Gaia era it is increasingly apparent that traditional static, parameterized models are insufficient to describe the mass distribution of our complex, dynamically evolving Milky Way (MW). In this work, we compare different time-evolving and time-independent representations of the gravitational potentials of simulated MW-mass galaxies from the FIRE-2 suite of cosmological baryonic simulations. Using these potentials, we calculate actions for star particles in tidal streams around three galaxies with varying merger histories at each snapshot from 7 Gyr ago to the present day. We determine the action-space coherence preserved by each model using the Kullback-Leibler Divergence to gauge the degree of clustering in actions and the relative stability of the clusters over time. We find that all models produce a clustered action space for simulations with no significant mergers. However, a massive (mass ratio prior to infall more similar than 1:8) interacting galaxy not present in the model will result in mischaracterized orbits for stars most affected by the interaction. The locations of the action space clusters (i.e. the orbits of the stream stars) are only preserved by the time-evolving model, while the time-independent models can lose significant amounts of information as soon as 0.5--1 Gyr ago, even if the system does not undergo a significant merger. Our results imply that reverse-integration of stream orbits in the MW using a fixed potential is likely to give incorrect results if integrated longer than 0.5 Gyr into the past.

Neutron stars provide a unique physical laboratory to study the properties of matter at high density. We study a diagnostic of the composition of high-density matter, namely, g-mode oscillations, which are driven by buoyancy forces. These oscillations can be excited by tidal forces and couple to gravitational waves. We extend prior results for the g-mode spectrum of cold neutron star matter to temperatures that are expected to be achieved in neutron star mergers using a parameterization for finite-temperature effects recently proposed by Raithel, \"Ozel and Psaltis. We find that the g-modes of canonical mass neutron stars ($\approx$1.4$M_{\odot}$) are suppressed at high temperature, and core $g$-modes are supported only in the most massive ($\geq $2$M_{\odot}$) of hot neutron stars.

The past decades have witnessed a lot of progress in gravitational lensing with two main targets: stars and galaxies (with active galactic nuclei). The success is partially attributed to the continuous luminescence of these sources making the detection and monitoring relatively easy. With the running of ongoing and upcoming large facilities/surveys in various electromagnetic and gravitational-wave bands, the era of time-domain surveys would guarantee constant detection of strongly lensed explosive transient events, for example, supernovae in all types, gamma ray bursts with afterglows in all bands, fast radio bursts and even gravitational waves. Lensed transients have many advantages over the traditional targets in studying the Universe and magnification effect helps to understand the transients themselves at high redshifts. In this review article, basing on the recent achievements in the literature, we summarize the methods of searching for different kinds of lensed transient signals, the latest results on detection and their applications in fundamental physics, astrophysics and cosmology. At the same time, we give supplementary comments as well as prospects of this emerging research direction that may help readers who are interested in entering this field.

We use the EAGLE cosmological simulations to perform a comprehensive and systematic analysis of the $z=0.1$ Fundamental Plane (FP), the tight relation between galaxy size, mass and velocity dispersion. We first measure the total mass and velocity dispersion (including both random and rotational motions) within the effective radius to show that simulated galaxies obey a total mass FP that is very close to the virial relation ($<10\%$ deviation), indicating that the effects of non-homology are weak. When we instead use the stellar mass, we find a strong deviation from the virial plane, which is driven by variations in the dark matter content. The dark matter fraction is a smooth function of the size and stellar mass, and thereby sets the coefficients of the stellar mass FP without substantially increasing the scatter. Hence, both star-forming and quiescent galaxies obey the same FP, with equally low scatter (0.02 dex). We employ simulations with a variable stellar initial mass function (IMF) to show that IMF variations have a modest additional effect on this FP. Moreover, when we use luminosity-weighted mock observations of the size and spatially-integrated velocity dispersion, the inferred FP changes only slightly. However, the scatter increases significantly, due to the luminosity-weighting and line-of-sight projection of the velocity dispersions, and measurement uncertainties on the half-light radii. Importantly, we find significant differences between the simulated FP and observations, which likely reflects a systematic difference in the stellar mass distributions. Therefore, we suggest the stellar mass FP offers a simple test for cosmological simulations, requiring minimal post-processing of simulation data.

We performed a unique Venus observation campaign to measure the disk brightness of Venus over a broad range of wavelengths in August and September 2020. The primary goal of the campaign is to investigate the absorption properties of the unknown absorber in the clouds. The secondary goal is to extract a disk mean SO$_2$ gas abundance, whose absorption spectral feature is entangled with that of the unknown absorber at the ultraviolet (UV) wavelengths. A total of 3 spacecraft and 6 ground-based telescopes participated in this campaign, covering the 52 to 1700~nm wavelength range. After careful evaluation of the observational data, we focused on the data sets acquired by 4 facilities. We accomplished our primary goal by analyzing the reflectivity spectrum of the Venus disk over the 283-800 nm wavelengths. Considerable absorption is present in the 350-450 nm range, for which we retrieved the corresponding optical depth by the unknown absorber. The result shows a consistent wavelength dependence of the relative optical depth with that at low latitudes during the Venus flyby by MESSENGER in 2007 (P\'erez-Hoyos et al. 2018), which was expected because the overall disk reflectivity is dominated by low latitudes. Last, we summarize the experience obtained during this first campaign that should enable us to accomplish our second goal in future campaigns.

In this review, we discuss various open-source software for modelling the broad-band emission of extragalactic sources from radio up to the highest gamma ray energies. As we provide an overview of the different tools available, we discuss the physical processes such tools implement and detail the computations they can perform. We also examine their conformity with modern good software practices. After considering the currently available software as a first generation of open-source modelling tools, we outline some desirable characteristics for the next generation.

Star-forming galaxies can exhibit strong morphological differences between the rest-frame far-UV and optical, reflecting inhomogeneities in star-formation and dust attenuation. We exploit deep, high resolution NIRCAM 7-band observations to take a first look at the morphology of galaxies in the epoch of reionization ($z>7$), and its variation in the rest-frame wavelength range between Lyman $\alpha$ and 6000-4000\AA, at $z=7-12$. We find no dramatic variations in morphology with wavelength -- of the kind that would have overturned anything we have learned from the Hubble Space Telescope. No significant trends between morphology and wavelengths are detected using standard quantitative morphology statistics. We detect signatures of mergers/interactions in 4/21 galaxies. Our results are consistent with a scenario in which Lyman Break galaxies -- observed when the universe is only 400-800 Myrs old - are growing via a combination of rapid galaxy-scale star formation supplemented by accretion of star forming clumps and interactions.

Gamma-ray pulsar halos are ideal indicators of cosmic-ray propagation in localized regions of the Galaxy and electron injection from pulsar wind nebulae. HESS~J1831$-$098 is a candidate pulsar halo observed by both H.E.S.S. and HAWC experiments. We adopt the flux map of the H.E.S.S. Galactic plane survey and the spectrum measurements of H.E.S.S. and \textit{Fermi}-LAT to study HESS~J1831$-$098. We find that HESS~J1831$-$098 meets all the criteria for a pulsar halo. The diffusion coefficient inside the halo and the conversion efficiency from the pulsar spin-down energy to the electron energy are both similar to the Geminga halo, a canonical pulsar halo. The injection spectrum can be well described by an exponentially-cutoff power law. However, the needed power-law term is very hard with $p\lesssim1$ if the diffusion coefficient is spatially and temporally independent. Considering the possible origins of the slow-diffusion environment, we adopt the two-zone diffusion model and the time-delayed slow-diffusion model. Both the models can interpret the H.E.S.S. and \textit{Fermi}-LAT results with a milder $p$. A modified injection time profile may have a similar effect.

I present an analytic model for the early post-collapse evolution of a spherical density peak on the coherence scale of the initial fluctuations in a universe filled with collisionless and pressure-free "dust". On a time-scale which is short compared to the peak's collapse time $t_0$, its inner regions settle into an equilibrium cusp with a power-law density profile, $\rho\propto r^{-12/7}$. Within this cusp, the circular orbit period $P$ at each radius is related to the enclosed mass $M$ by $P = t_0 (M/M_c)^{2/3}$ where $M_c$ is a suitably defined characteristic mass for the initial peak. The relaxation mechanism which produces this cusp gives insight into those which are active in high-resolution simulations of first halo formation in Cold or Warm Dark Matter universes, and, indeed, a simple argument suggests that the same power-law index $\gamma=-12/7$ should describe the prompt cusps formed during the collapse of generic peaks, independent of any symmetry assumption. Further work is needed to investigate whether additional factors are required to explain the slightly flatter exponent, $\gamma\approx -1.5$, found in high-resolution numerical simulations of peak collapse.

Whether there is a quark matter core in the neutron star (NS) is a fundamental question. The increasing multi-messenger data set of NSs provide a valuable chance to examine such an attractive possibility. Here we carry out the Bayesian nonparametric inference of the NS equation of state (EOS) via a single-layer feed-forward neural network, taking into account the data of GW170817, PSR J0030+0451, and PSR J0740+6620, and incorporating the latest constraints from the chiral effective theory ($\chi$EFT) and perturbative quantum chromodynamics (pQCD) at low and very high energy densities, respectively. It is found out that a sizable quark matter core ($\geq 10^{-3}M_\odot$) is plausible ($\geq 90\%$ probability) for the very massive NS with a gravitational mass above about $0.97M_{\rm TOV}$, where $M_{\rm TOV}$, the maximum gravitational mass of a non-rotating cold NS, is simultaneously constrained to be $2.18^{+0.27}_{-0.13}M_\odot$ ($90\%$ credibility). The average density of the quark matter core is found to be $\sim 2.2$ times that of the host NS. A few percent of the posterior EOSs, which do not predict quark matter cores even in the heaviest NSs, are characterized by a quicker rising of sound speed at relatively low densities. We also find that sound speed may reach close to zero near the center density of NS with $M \approx M_{\rm TOV}$ and hence only allows the presence of the strong first-order phase transition in the center of the most massive NSs.

Measurements of short-timescale star formation variations (i.e., "burstiness") are integral to our understanding of star formation feedback mechanisms and the assembly of stellar populations in galaxies. We expand upon the work of Broussard et al. (2019) by introducing a new analysis of galaxy star formation burstiness that accounts for variations in the $Q_{sg}=\mathrm{E(B-V)_{stars}}/\mathrm{E(B-V)_{gas}}$ distribution, a major confounding factor. We use Balmer decrements from the MOSFIRE Deep Evolution Field (MOSDEF) survey to measure $Q_{sg}$, which we use to construct mock catalogs from the Santa Cruz Semi-Analytic Models and Mufasa cosmological hydrodynamical simulation based on 3D-HST, Fiber Multi-Object Spectrograph (FMOS)-COSMOS, and MOSDEF galaxies with H$\alpha$ detections. The results of the mock catalogs are compared against observations using the burst indicator $\eta = \log_{10}(\mathrm{SFR_{H\alpha}/SFR_{NUV}})$, with the standard deviation of the $\eta$ distribution indicating burstiness. We find decent agreement between mock and observed $\eta$ distribution shapes; however, the FMOS-COSMOS and MOSDEF mocks show a systematically low median and scatter in $\eta$ in comparison to the observations. This work also presents the novel approach of analytically deriving the relationship between the intrinsic scatter in $\eta$, scatter added by measurement uncertainties, and observed scatter, resulting in an intrinsic burstiness measurement of $0.06-0.16$ dex.

How passive galaxies form, and the physical mechanisms which prevent star formation over long timescales, are some of the most outstanding questions in understanding galaxy evolution. The properties of quiescent galaxies over cosmic time provide crucial information to identify the quenching mechanisms. Passive galaxies have been confirmed and studied out to $z\sim4$, but all of these studies have been limited to massive systems (mostly with $\log{(M_{\rm star}/M_{\odot})}>10.8$). Using James Webb Space Telescope (JWST) NIRISS grism slitless spectroscopic data from the GLASS JWST ERS program, we present spectroscopic confirmation of two quiescent galaxies at $z_{\rm spec}=2.650^{+0.004}_{-0.006}$ and $z_{\rm spec}=2.433^{+0.032}_{-0.016}$ (3$\sigma$ errors) with stellar masses of $\log{(M_{\rm star}/M_{\odot})}=10.53^{+0.18}_{-0.06}$ and $\log{(M_{\rm star}/M_{\odot})}=9.93^{+0.06}_{-0.07}$ (corrected for magnification factors of $\mu=2.0$ and $\mu=2.1$, respectively). The latter represents the first spectroscopic confirmation of the existence of low-mass quiescent galaxies at cosmic noon, showcasing the power of JWST to identify and characterize this enigmatic population.

Using a data-driven machine learning tool we report $T_{\text{eff}}$, $\log{(g)}$, $v\sin{(i)}$, and elemental abundances for 15 elements (C, N, O, Na, Mg, Al, Si, Ca, Ti, V, Cr, Mn, Fe, Ni, Y) for a sample of 25 exoplanet host stars targeted by JWST's first year of observations. The chemical diversity of these stars show that, while a number of their companion planets may have formed in a disk with chemistry similar to Solar, some JWST targets likely experienced different disk compositions. This sample is part of a larger forthcoming catalog that will report homogeneous abundances of $\sim$4,500 FGK stars derived from Keck/HIRES spectra.

FRB 20180916B is a repeating fast radio burst (FRB) with an activity period of 16.33 days. In previous observations ranging from $\sim 150-1400$ MHz, the activity window was found to be frequency dependent, with lower frequency bursts occurring later. In this work, we present the highest-frequency detections of bursts from this FRB, using the 100-m Effelsberg Radio Telescope at 4$-$8 GHz. We present the results from two observing campaigns. We performed the first campaign over an entire activity period which resulted in no detections. The second campaign was in an active window at 4$-$8 GHz which we predicted from our modelling of chromaticity, resulting in eight burst detections. The bursts were detected in a window of 1.35 days, 3.6 days preceding the activity peak seen by CHIME, suggesting the chromaticity extends to higher frequency. The detected bursts have narrower temporal widths and larger spectral widths compared to lower frequencies. All of them have flat polarization position angle sweeps and high polarization fractions. The bursts also exhibit diffractive scintillation due to the Milky Way, following a $f^{3.90\pm0.05}$ scaling, and vary significantly over time. We find that burst rate across frequency scales as $f^{-2.6\pm0.2}$. Lastly, we examine implications of the frequency dependency on the source models.

We provide an accurate comparison, against large cosmological $N$-body simulations, of different prescriptions for modelling nonlinear matter power spectra in the presence of massive neutrinos and dynamical dark energy. We test the current most widely used approaches: fitting functions (HALOFIT and HMcode), the halo-model reaction (ReACT) and emulators (baccoemu and EuclidEmulator2). Focussing on redshifts $z\leq2$ and scales $k\lesssim 1 \ h/$Mpc (where the simulation mass resolution provides $\sim 1\%$ accuracy), we find that HMcode and ReACT considerably improve over the HALOFIT prescriptions of Smith and Takahashi (both combined with the Bird correction), with an overall agreement of 2\% for all the cosmological scenarios considered. Concerning emulators, we find that, especially at low redshifts, EuclidEmulator2 remarkably agrees with the simulated spectra at $\lesssim 1\%$ level in scenarios with dynamical dark energy and massless neutrinos, reaching a maximum difference of $\sim 2\%$ at $z=2$. baccoemu has a similar behaviour as EuclidEmulator2, except for a couple of dark energy models. In cosmologies with massive neutrinos, at $z=0$ all the nonlinear prescriptions improve their agreement with respect to the massless neutrino case, except for the Bird and TakaBird models which, however, are not tailored to $w_0$--$w_a$ models. At $z>0$ we do not find a similar improvement when including massive neutrinos, probably due to the lower impact of neutrino free-streaming at higher redshifts; rather at $z=2$ EuclidEmulator2 exceeds $2\%$ agreement for some dark energy EoS. When considering ratios between the matter power spectrum computed in a given cosmological model and its $\Lambda$CDM counterpart, all the tested prescriptions agree with simulated data, at sub-percent or percent level, depending on $z$. [ABRIDGED]

We investigate an interacting dark energy model which allows for the kinetic term of the scalar field to couple to dark matter via a power-law interaction. The model is characterised by scaling solutions at early times, which are of high interest to alleviate the coincidence problem, followed by a period of accelerated expansion. We discuss the phenomenology of the background evolution and of the linear scalar perturbations and we identify measurable signatures of the coupling in the dark sector on the cosmic microwave background, the lensing potential auto-correlation and the matter power spectra. We also perform a parameter estimation analysis using data of cosmic microwave background temperature, polarisation and lensing, baryonic acoustic oscillations and supernovae. We find that the strength of the coupling between the dark sectors, regulated by the parameter $\alpha$, is constrained to be of order $10^{-4}$. A model selection analysis does not reveal a statistical preference between $\Lambda$CDM and the Kinetic model.

We introduce a modification of the Press-Schechter formalism aimed to derive general mass functions for primordial black holes (PBHs). In this case, we start from primordial power spectra (PPS) which include a monochromatic spike, typical of ultra slow-roll inflation models. We consider the PBH formation as being associated to the amplitude of the spike on top of the linear energy density fluctuations, coming from a PPS with a blue index. By modelling the spike with a log-normal function, we study the properties of the resulting mass function spikes, and compare these to the underlying extended mass distributions. When the spike is at PBH masses which are much lower than the exponential cutoff of the extended distribution, very little mass density is held by the PBHs within the spike, and there is little statistical connection between the power spectrum and the resulting mass function. Instead, the connection is strong when the spike mass is similar to, or larger than the cutoff mass, and it can hold a similar mass density as the extended part. Such particular mass functions also contain large numbers of small PBHs, especially if stable PBH relics are considered. The constraints on the fraction of dark matter in PBHs for monochromatic mass functions are somewhat relaxed when there is an additional underlying extended distribution of masses.

Galaxy mergers provide a mechanism for galaxies to effectively funnel gas and materials toward their nuclei and fuel the central starbursts and accretion of supermassive black holes. In turn, the active nuclei drive galactic-scale outflows that subsequently impact the evolution of the host galaxies. The details of this transformative process as they pertain to the supermassive black holes remain ambiguous, partially due to the central obscuration commonly found in the dust-reddened merger hosts, and also because there are relatively few laboratories in the nearby universe where the process can be studied in depth. This review highlights the current state of the literature on the role of accreting supermassive black holes in local luminous infrared galaxies as seen from various windows within the electromagnetic spectrum. Specifically, we discuss the multiwavelength signatures of the active nucleus, its associated feeding and feedback processes, and the implications of multiple supermassive black holes found in nearby interacting galaxy systems for galaxy evolution from the observational perspective. We conclude with a future outlook on how the topic of active nuclei in low- and high-redshift galaxy mergers will benefit from the advent of next-generation observing facilities with unparalleled resolving power and sensitivity in the coming decade.

We present a detailed analysis of the rest-frame optical emission line ratios for three spectroscopically confirmed galaxies at $z>7.5$. The galaxies were identified in the \emph{James Webb Space Telescope} (\emph{JWST}) Early Release Observations field SMACS J0723.3$-$7327. By quantitatively comparing Balmer and oxygen line ratios of these galaxies with various low-redshift "analogue" populations (e.g. Green Peas, Blueberries, etc.), we show that no single analogue population captures the diversity of line ratios of all three galaxies observed at $z>7.5$. We find that S06355 at $z=7.67$ and S10612 at $z=7.66$ are similar to local Green Peas and Blueberries. In contrast, S04590 at $z=8.50$ appears to be significantly different from the other two galaxies, most resembling extremely low-metallicity systems in the local Universe. Perhaps the most striking spectral feature in S04590 is the curiously high [O {\small III}]\,$\lambda4363$/[O {\small III}]\,$\lambda5007$ ratio (RO3) of $0.047$ (or $0.059$ when dust-corrected), implying either extremely high electron temperatures, $>3\times10^4$~K, or gas densities $>10^4\ {\rm cm^{-3}}$. Observed line ratios indicate that this galaxy is unlikely to host an AGN. Using photoionization modelling, we show that the inclusion of high-mass X-ray binaries or a high cosmic ray background in addition to a young, low-metallicity stellar population can provide the additional heating necessary to explain the observed high RO3 while remaining consistent with other observed line ratios. Our models represent a first step at accurately characterising the dominant sources of photoionization and heating at very high redshifts, demonstrating that non-thermal processes may become important as we probe deeper into the Epoch of Reionization.

In many hypothesis testing applications, we have mixed priors, with well-motivated informative priors for some parameters but not for others. The Bayesian methodology uses the Bayes factor and is helpful for the informative priors, as it incorporates Occam's razor via multiplicity or trials factor in the Look Elsewhere Effect. However, if the prior is not known completely, the frequentist hypothesis test via the false positive rate is a better approach, as it is less sensitive to the prior choice. We argue that when only partial prior information is available, it is best to combine the two methodologies by using the Bayes factor as a test statistic in the frequentist analysis. We show that the standard frequentist likelihood-ratio test statistic corresponds to the Bayes factor with a non-informative Jeffrey's prior. We also show that mixed priors increase the statistical power in frequentist analyses over the likelihood ratio test statistic. We develop an analytic formalism that does not require expensive simulations using a statistical mechanics approach to hypothesis testing in Bayesian and frequentist statistics. We introduce the counting of states in a continuous parameter space using the uncertainty volume as the quantum of the state. We show that both the p-value and Bayes factor can be expressed as energy versus entropy competition. We present analytic expressions that generalize Wilks' theorem beyond its usual regime of validity and work in a non-asymptotic regime. In specific limits, the formalism reproduces existing expressions, such as the p-value of linear models and periodograms. We apply the formalism to an example of exoplanet transits, where multiplicity can be more than $10^7$. We show that our analytic expressions reproduce the p-values derived from numerical simulations.

We analyse the dynamical properties of disformally transformed theories of gravity. We show that disformal transformation typically introduces novel degrees of freedom, equivalent to the mimetic dark matter, which possesses a Weyl-invariant formulation. We demonstrate that this phenomenon occurs in a wider variety of disformal transformations than previously thought.

Resonant excitations of $f$-modes in binary neutron star coalescences influence the gravitational waves (GWs) emission in both quasicircular and highly eccentric mergers and can deliver information on the star interior. Most models of resonant tides are built using approximate, perturbative approaches and thus require to be carefully validated against numerical relativity (NR) simulations in the high-frequency regime. We perform detailed comparisons between a set of high-resolution NR simulations and the state of the art effective one body (EOB) model ${\tt TEOBResumS}$ with various tidal potentials and including a model for resonant tides. For circular mergers, we find that $f$-mode resonances can improve the agreement between EOB and NR, but there is no clear evidence that the tidal enhancement after contact is due to a resonant mechanism. Tidal models with $f$-mode resonances do not consistently reproduce, at the same time, the NR waveforms and the energetics within the errors, and their performances is comparable to resummed tidal models without resonances. For highly eccentric mergers, we show for the first time that our EOB model reproduces the bursty NR waveform to a high degree of accuracy. However, the considered resonant model does not capture the $f$-mode oscillations excited during the encounters and present in the NR waveform. Finally, we analyze GW170817 with both adiabatic and dynamical tides models and find that the data shows no evidence in favor of models including dynamical tides. This is in agreement with the fact that resonant tides are measured at very high frequencies, which are not available for GW170817 but might be tested with next generation detectors.

While most searches for cosmic axions so far focused on their cold relics as (a component of) dark matter, various well-motivated cosmological sources can produce "boosted" axions that remain relativistic today. We demonstrate that existing/upcoming neutrino experiments such as Super-Kamiokande, Hyper-Kamiokande, DUNE, JUNO, and IceCube can probe such energetic axion relics. The characteristic signature is the mono-energetic single photon signal from axion absorption induced by the axion-photon coupling. This proposal offers to cover parameter ranges that are complementary to existing axion searches and provides new opportunities for discovery with neutrino facilities.

We propose a unimodular version of the Brans-Dicke theory designed with a constrained Lagrangian formulation. The resulting field equations are traceless. The vacuum solutions in the cosmological background reproduce the corresponding solutions of the usual Brans-Dicke theory but with a cosmological constant term. A perturbative analysis of the scalar modes is performed and stable and unstable configurations appear in contrast with the Brans-Dicke case for which only stable configurations occur. On the other hand, tensorial modes in this theory remains the same as in the traditional Brans-Dicke theory.

We present results of a performance study of an astrophysical radiation hydrodynamics code, V2D, on the Arm-based A64FX processor developed by Fujitsu. The code solves sparse linear systems, a task for which the A64FX architecture should be well suited. We performed the performance analysis study on Ookami, an Apollo 80 platform utilizing the A64FX processor. We explored several compilers and performance analysis packages and found the code did not perform as expected under scalable vector extension optimization, suggesting that a "deeper dive" into analyzing the code is worthwhile. However, a simple driver program that exercised basic sparse linear algebra routines used by V2D did show significant speedup with the use of the scalable vector extension optimization. We present the initial results from the study which used V2D on a relatively simple test problem that emphasized the repeated solution of sparse linear systems.

In this paper, we construct the first asymmetric strongly interacting massive particles (SIMP) dark matter (DM) model, where a new vector-like fermion and a new complex scalar both having nonzero chemical potentials can be asymmetric DM particles. After the spontaneous breaking of a U(1)$^{}_\textsf{D}$ dark gauge symmetry, these two particles can have accidental $\mathbb{Z}^{}_4$ charges making them stable. By adding one more complex scalar as a mediator between the SIMP DM, the relic density of DM is determined by $3 \to 2$ and two-loop induced $2 \to 2$ annihilations in this model. On the other hand, the SIMP DM can maintain kinetic equilibrium with the thermal bath until the DM freeze-out temperature via the new gauge interaction. Interestingly, this model can have a bouncing effect on DM, whereby the DM number density rises after the chemical freeze-out of DM. With this effect, the prediction of the DM self-interacting cross section in this model can be consistent with astrophysical observations, and the ratio of the DM energy density to the baryonic matter energy density can be explained by primordial asymmetries. We also predict the DM-electron elastic scattering cross section that can be used to test this model in future projected experiments.

In order to model neutron stars, an equation of state (EoS) correlating total energy density with pressure is essential for providing physically allowable mass-radius relationship of neutron stars and calculating the neutron star masses. The mass, radius and crustal fraction of moment of inertia in neutron stars have been determined using $\beta$-equilibrated dense $npe\mu$ neutron star matter obtained using the Skyrme effective interaction with the NRAPR set. The maximum mass of neutron star calculated from this set is able to reach $\sim$2$M_\odot$, the mass of highly massive compact stars. The study of pulsar glitches facilitates extraction of crustal fraction of the moment of inertia. This fraction is greatly dependent on the pressure and corresponding density at core-crust transition. The core-crust transition density and pressure together with the observed minimum crustal fraction of the total moment of inertia provide a limit for the radius of the Vela pulsar, $R \geq 3.65 + 3.67 M/M_\odot$ kms. Present calculations imply that due to crustal entrainment coming from the Bragg reflection of unbound neutrons by lattice ions, the crustal fraction of the total moment of inertia is $\sim$6.3$\%$.

We present efforts at improving the performance of FLASH, a multi-scale, multi-physics simulation code principally for astrophysical applications, by using huge pages on Ookami, an HPE Apollo 80 A64FX platform. FLASH is written principally in modern Fortran and makes use of the PARAMESH library to manage a block-structured adaptive mesh. We explored options for enabling the use of huge pages with several compilers, but we were only able to successfully use huge pages when compiling with the Fujitsu compiler. The use of huge pages substantially reduced the number of translation lookaside buffer misses, but overall performance gains were marginal.