Papers for Thursday, Aug 04 2022

In this paper we leverage Gaia parallactic distances to deliver revised estimates of planetary parameters and extreme ultraviolet (EUV) fluxes for a distance-limited ($\lesssim 100$ pc) sample of 27 gaseous planets -- from super-Earths to hot Jupiters -- with publicly available X-ray observations. Notably, the revised mass and radius imply a Saturn-like density ($0.86 \pm 0.09$ g cm$^{-3}$) for HD 149026 b -- consistent with the lowest values reported in the literature -- thus removing the need for a high metal fraction. For 20 planets with X-ray detected host stars we also derive updated atmospheric mass outflow rates making use of the 1D photoionization hydrodynamics code ATES. We note, however, how X-ray variability combined with large uncertainties in the conversion between X-rays and EUV fluxes severely affects the inferred instantaneous mass loss rates, and thus the ability to constrain/predict the integrated mass loss over a planet's lifetime.

We present a study of spatially-resolved star formation histories (SFHs) for 60 $z\sim2.3$ main-sequence, star-forming galaxies selected from the MOSDEF spectroscopic survey in the GOODS-N field. Photometry is decomposed into a central and outer spatial component using observed $z_\mathrm{F850LP}-H_\mathrm{F160W}$ colors. The Prospector code is used to model spectral energy distributions for the centers, outskirts, and integrated galaxy using HST/ACS and WFC3, Spitzer/IRAC, and ground-based photometry, with additional constraints on metallicity and spectroscopic redshift from MOSDEF spectroscopy. For the low-resolution bands, spatially-resolved photometry is determined with an iterative approach. The reconstructed SFHs indicate that the majority of galaxies with $\log(M_\star/M_\odot)<10.5$ are observed while their central regions undergo relatively recent ($<100$ Myr) bursts of star formation, while the outskirts have a smooth, quasi-steady SFH. The enhanced star formation activity of the central parts is broadly consistent with the idea that it is produced by highly dissipative gas compaction and accretion. The broad dispersion of central density and size observed in the sample suggests that for the selected galaxies this process has started but is still far from being completed. The implication would be that selecting star-forming galaxies at cosmic noon frequently includes systems in an "evolved" evolutionary phase where the centers have recently started a burst of star formation activity that will likely initiate inside-out quenching in the next several hundred million years.

The Fermi bubbles are large gamma-ray-emitting structures. They are symmetric about the Galactic Centre (GC), and their creation is therefore attributed to intensive energy injection at the GC. In this study, we focus on the X-ray gas structures associated with the bubbles. We show that a combination of the density, temperature, and shock age profiles of the X-ray gas can be used to distinguish the energy injection mechanisms. By comparing the results of numerical simulations with observations, we indicate that the bubbles were created by a fast wind from the GC because it generates a strong reverse shock and reproduces the observed temperature peak there. On the other hand, instantaneous energy injection at the GC cannot reproduce the temperature profile. The wind had a speed of ~1000 km s^-1, and blew for ~10^7 yr. Because the mass flux of the wind is large, the entrainment of interstellar gas by wide-angle outflows from the black hole is required. Thus, the wind may be the same as active galactic nuclei outflows often observed in other galaxies and thought to regulate the growth of galaxies and their central black holes.

Observations of the S-stars, the cluster of young stars in the inner 0.1 pc of the Galactic Center, have been crucial in providing conclusive evidence for a supermassive black hole at the center of our galaxy. Since some of the stars have orbits less than that of a typical human lifetime, it is possible to observe multiple orbits and test the weak-field regime of general relativity. Current calculations of S-star orbits require slow and expensive computations in order to numerically solve geodesic equations for many small time steps. In this paper, we present a computationally efficient, first-order post-Newtonian model for the astrometric and spectroscopic data gathered for the S-stars. We find that future, 30-m class telescopes -- and potentially even current large telescopes with very high spectroscopic resolution -- may be able to detect the Shapiro effect for an S-star in the next decade or so.

Direct imaging studies have mainly used low-resolution spectroscopy ($R\sim20-100$) to study the atmospheres of giant exoplanets and brown dwarf companions, but the presence of clouds has often led to degeneracies in the retrieved atmospheric abundances (e.g. C/O, metallicity). This precludes clear insights into the formation mechanisms of these companions. The Keck Planet Imager and Characterizer (KPIC) uses adaptive optics and single-mode fibers to transport light into NIRSPEC ($R\sim35,000$ in $K$ band), and aims to address these challenges with high-resolution spectroscopy. Using an atmospheric retrieval framework based on petitRADTRANS, we analyze KPIC high-resolution spectrum ($2.29-2.49~\mu$m) and archival low-resolution spectrum ($1-2.2~\mu$m) of the benchmark brown dwarf HD 4747 B ($m=67.2\pm1.8~M_{\rm{Jup}}$, $a=10.0\pm0.2$ au, $T_{\rm eff}\approx1400$ K). We find that our measured C/O and metallicity for the companion from the KPIC high-resolution spectrum agree with that of its host star within $1-2\sigma$. The retrieved parameters from the $K$ band high-resolution spectrum are also independent of our choice of cloud model. In contrast, the retrieved parameters from the low-resolution spectrum are highly sensitive to our chosen cloud model. Finally, we detect CO, H$_2$O, and CH$_4$ (volume mixing ratio of log(CH$_4$)=$-4.82\pm0.23$) in this L/T transition companion with the KPIC data. The relative molecular abundances allow us to constrain the degree of chemical disequilibrium in the atmosphere of HD 4747 B, and infer a vertical diffusion coefficient that is at the upper limit predicted from mixing length theory.

One of the toughest challenges in modern cosmology is to probe the small scales $k \gtrsim 0.5\,{\rm Mpc}^{-1}$ in the matter power spectrum and clustering. Here, information on the nature of dark matter and on the shape of the primordial power spectrum is retained; however, uncertainties are large due to the non-linearity of structure formation and to the contribution of astrophysical processes. We show that such small scales will be accessible via upcoming line-intensity mapping surveys, with carbon monoxide (CO) emission from star-forming galaxies at high redshifts as an example. While these galaxies cannot be individually detected and the two-point correlations of the intensity fluctuation maps is not accessible at these scales, the voxel intensity distribution (VID) of the highly non-Gaussian intensity maps is sensitive to the integrated emission from faint sources. We explore the prospects of planned CO surveys targeting $z\sim 3$ and the epoch of reionization, $z \sim 6$, and show that their VID can probe deviations from $\Lambda$CDM as small as $\lesssim10\%$ at least up to $k\sim 10\,{\rm Mpc}^{-1}$.

We present the results of a stellar-population analysis of Lyman-alpha emitting galaxies (LAES) in GOODS-N at 1.9 < z < 3.5 spectroscopically identified by the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX). We provide a method for connecting emission-line detections from the blind spectroscopic survey to imaging counterparts, a crucial tool needed as HETDEX builds a massive database of ~1 million Lyman-alpha detections. Using photometric data spanning as many as 11 filters covering 0.4-4.5 microns from the Hubble and Spitzer Space Telescopes, we study the objects' global properties and explore which properties impact the strength of Lyman-alpha emission. We measure a median stellar mass of 0.8 (^+2.9_-0.5) x 10^9 Msol and conclude that the physical properties of HETDEX spectroscopically-selected LAEs are comparable to LAEs selected by previous deep narrow band studies. We find that stellar mass and star formation rate correlate strongly with the Lyman-alpha equivalent width. We then use a known sample of z>7 LAEs to perform a proto-study of predicting Lyman-alpha emission from galaxies in the Epoch of Reionization, finding agreement at the 1-sigma level between prediction and observation for the majority of strong emitters.

The blazar J1924-2914 is a primary Event Horizon Telescope (EHT) calibrator for the Galactic Center's black hole Sagittarius A*. Here we present the first total and linearly polarized intensity images of this source obtained with the unprecedented 20 $\mu$as resolution of the EHT. J1924-2914 is a very compact flat-spectrum radio source with strong optical variability and polarization. In April 2017 the source was observed quasi-simultaneously with the EHT (April 5-11), the Global Millimeter VLBI Array (April 3), and the Very Long Baseline Array (April 28), giving a novel view of the source at four observing frequencies, 230, 86, 8.7, and 2.3 GHz. These observations probe jet properties from the subparsec to 100-parsec scales. We combine the multi-frequency images of J1924-2914 to study the source morphology. We find that the jet exhibits a characteristic bending, with a gradual clockwise rotation of the jet projected position angle of about 90 degrees between 2.3 and 230 GHz. Linearly polarized intensity images of J1924-2914 with the extremely fine resolution of the EHT provide evidence for ordered toroidal magnetic fields in the blazar compact core.

In this paper, an investigation on the relative importance of the thermal gas, radiation, and (minimum-energy) magnetic pressures around $\approx$200 star-forming regions in a sample of nearby normal and luminous infrared galaxies is presented. Given the range of galaxy distances, pressure estimates are made on spatial scales spanning $\sim$0.1$-3$kpc. The ratio of thermal gas-to-radiation pressures does not appear to significantly depend on star formation rate surface density ($\Sigma_{\rm SFR}$), but exhibits a steady decrease with increasing physical size of the aperture over which the quantities are measured. The ratio of magnetic-to-radiation pressures appears to be relatively flat as a function of $\Sigma_{\rm SFR}$ and similar in value for both nuclear and extranuclear regions, but unlike the ratio of thermal gas-to-radiation pressures, exhibits a steady increase with increasing aperture size. Furthermore, it seems that the magnetic pressure is typically weaker than the radiation pressure on sub-kpc scales, and only starts to play a significant role on few-kpc scales. When the internal pressure terms are summed, their ratio to the ($\Sigma_{\rm SFR}$-inferred) kpc-scale dynamical equilibrium pressure estimates is roughly constant. Consequently, it appears that the physical area of the galaxy disk, and not necessarily environment (e.g., nuclear vs. extranuclear regions) or star formation activity, may play the dominant role in determining which pressure term is most active around star-forming regions. These results are consistent with a scenario in which a combination of processes acting primarily on different physical scales work collectively to regulate the star formation process in galaxy disks.

We present the MaNGA Dwarf galaxy, MaNDala, Value-Added-Catalog, VAC, from the final release of the Sloan Digital Sky Survey-IV program. MaNDala consists of 136 randomly selected bright dwarf galaxies with $M_{*}<10^{9.1}\odot$ and $M_{g}>-18.5$ making it the largest Integral Field Spectroscopy homogeneous sample of dwarf galaxies. We release a photometric analysis of the $g,r$ and $z$ broadband imaging based on the DESI Legacy Imaging Surveys as well as the spectroscopic analysis based on the Pipe3D SDSS-IV VAC. Our release includes the surface brightness (SB), geometric parameters and color profiles, S\'ersic fits as well as stellar population properties (such as, stellar ages, metallicities, star formation histories), and emission lines fluxes within the FOV and the effective radii of the galaxies. We find that the majority of the MaNDala galaxies are star forming late-type galaxies with $\langle{}n_{\text{Sersic,r}}\rangle\sim1.6$ that are centrals (central/satellite dichotomy). MaNDala covers a large range of SB values (we find 11 candidates of ultra diffuse galaxies and 3 compact ones), filling the gap between classical dwarfs and low-mass galaxies in the Kormendy Diagram and in the size-mass/luminosity relation, whichseems to flatten at $10^8<M_{*}/\odot<10^{9}$ with $\langle{}R_{e,r}\rangle\sim2.7$ kpc. A large fraction of MaNDala galaxies formed from an early low-metallicity burst of SF but also of late SF events from more metal-enriched gas: half of the MaNDala galaxies assembled $50\%$ of their mass at $\langle{}z\rangle>2$, while the last $20\%$ was at $\langle{}z\rangle<0.3$. Finally, a bending of the sSFR-$M_{*}$ relation at $M_{*}\sim10^{9}\odot$ for the main sequence galaxies seems to be supported by MaNDala.

It is often stated that the removal of gas by ram pressure stripping of a galaxy disk is not a common process in galaxy groups. In this study, with the aid of an observational classification of galaxies and a simple physical model, we show that this may not be true. We examined and identified 45 ram-pressure stripped galaxy candidates from a sample of 1311 galaxy group members within 125 spectroscopically-selected galaxy groups. 13 of these galaxies are the most secure candidates with multiple distinct features. These candidate ram pressure stripped galaxies have similar properties to those found in clusters -- they occur at a range of stellar masses, are largely blue and star-forming and have phase-space distributions consistent with being first infallers into their groups. The only stand-out feature of these candidates is they exist not in clusters, but in groups, with a median halo mass of 10$^{13.5}$ M$_\odot$. Although this may seem surprising, we employ an analytic model of the expected ram pressure stripping force in groups and find that reasonable estimates of the relevant infall speeds and intragroup medium content would result in ram pressure stripped galaxies at these halo masses. Finally, given the considerable uncertainty on the lifetime of the ram-pressure phase, this physical mechanism could be the dominant quenching mechanism in galaxy groups, if our ram pressure stripped candidates can be confirmed.

The sparse interferometric coverage of the Event Horizon Telescope (EHT) poses a significant challenge for both reconstruction and model fitting of black-hole images. PRIMO is a new principal components analysis-based algorithm for image reconstruction that uses the results of high-fidelity general relativistic, magnetohydrodynamic simulations of low-luminosity accretion flows as a training set. This allows the reconstruction of images that are both consistent with the interferometric data and that live in the space of images that is spanned by the simulations. PRIMO follows Monty Carlo Markov Chains to fit a linear combination of principal components derived from an ensemble of simulated images to interferometric data. We show that PRIMO can efficiently and accurately reconstruct synthetic EHT data sets for several simulated images, even when the simulation parameters are significantly different from those of the image ensemble that was used to generate the principal components. The resulting reconstructions achieve resolution that is consistent with the performance of the array and do not introduce significant biases in image features such as the diameter of the ring of emission.

Given that mergers are often invoked to explain many exotic phenomena in massive star evolution, understanding the evolutionary phase directly preceding a merger, the overcontact phase, is of crucial importance. Despite its importance, large uncertainties exist in our understanding of the evolution of massive overcontact binaries. We aim to provide robust observational constraints on the future dynamical evolution of massive overcontact systems by measuring the rate at which the periods change for a sample of six such objects. Furthermore, we aim to investigate whether the periods of unequal mass systems show higher rates of change than their equal mass counterparts as theoretical models predict. Using archival photometric data from various ground- and space-based missions covering up to ~40 years, we measure the periods of each system over several smaller time spans. We then fit a linear regression through the measured periods to determine the rate at which the period is changing over the entire data set. We find that all of the stars in our sample have very small period changes and that there does not seem to be a correlation with the mass ratio. This implies that the orbital periods for these systems are stable on the nuclear timescale, and that the unequal mass systems may not equalize as expected. When comparing our results with population synthesis distributions, we find large discrepancies between the expected mass ratios and period stabilities. We find that these discrepancies can be mitigated to a degree by removing systems with shorter initial periods, suggesting that the observed sample of overcontact systems may originate from binary systems with longer initial orbital periods.

In the era of Extremely Large Telescopes, the current generation of 8-10m facilities are likely to remain competitive at ground-UV wavelengths for the foreseeable future. The Cassegrain U-Band Efficient Spectrograph (CUBES) has been designed to provide high-efficiency (>40%) observations in the near UV (305-400 nm requirement, 300-420 nm goal) at a spectral resolving power of R>20,000 (with a lower-resolution, sky-limited mode of R ~ 7,000). With the design focusing on maximizing the instrument throughput (ensuring a Signal to Noise Ratio (SNR) ~20 per high-resolution element at 313 nm for U ~18.5 mag objects in 1h of observations), it will offer new possibilities in many fields of astrophysics, providing access to key lines of stellar spectra: a tremendous diversity of iron-peak and heavy elements, lighter elements (in particular Beryllium) and light-element molecules (CO, CN, OH), as well as Balmer lines and the Balmer jump (particularly important for young stellar objects). The UV range is also critical in extragalactic studies: the circumgalactic medium of distant galaxies, the contribution of different types of sources to the cosmic UV background, the measurement of H2 and primordial Deuterium in a regime of relatively transparent intergalactic medium, and follow-up of explosive transients. The CUBES project completed a Phase A conceptual design in June 2021 and has now entered the detailed design and construction phase. First science operations are planned for 2028.

Detection of the first stars has remained elusive so-far but their presence may soon be unveiled by upcoming JWST observations. Previous studies have not investigated the entire possible range of halo masses and redshifts which may help in their detection. Motivated by the prospects of detecting galaxies up to $z\sim 20$ in JWST early data release, we quantify the contribution of Pop III stars to high-redshift galaxies from $10 \leq z \leq 26$ by employing the semi-analytical model A-SLOTH, which self-consistently models the formation of Pop III and Pop II stars along with their feedback. Our results suggest that the contribution of Pop III stars is the highest in low-mass halos of $\rm 10^7-10^9~M_{\odot}$. While high-mass halos $\rm \geq 10^{10}~M_{\odot}$ contain less than 1\% Pop III stars, they host galaxies with stellar masses of $\rm 10^9~M_{\odot}$ as early as $z \sim 26$. Interestingly, the apparent magnitude of Pop~III populations gets brighter towards higher redshift, but Pop~III-dominated galaxies are too faint to be directly detected with JWST. Our results predict JWST can detect galaxies up to $z\sim 26$, which may help in constraining the IMF of Pop III stars and will guide observers to discern the contribution of Pop~III stars to high-redshift galaxies.

High-mass star formation theories make distinct predictions on the properties of the prestellar seeds of high-mass stars. Observations of the early stages of high-mass star formation can provide crucial constraints, but they are challenging and scarce. We investigate the properties of the prestellar core population embedded in the high-mass clump AGAL014.492-00.139, and we study the kinematics at the clump and the clump-to-core scales. We have analysed an extensive dataset acquired with the ALMA interferometer. Applying a dendrogram analysis to the Band o-$\rm H_2D^+$ data, we identified 22 cores. We have fitted their average spectra in local-thermodinamic-equilibrium conditions, and we analysed their continuum emission at $0.8 \, \rm mm$. The cores have transonic to mildly supersonic turbulence levels and appear mostly low-mass, with $M_\mathrm{core}< 30 \, \rm M_\odot$. Furthermore, we have analysed Band 3 observations of the $\rm N_2H^+$ (1-0) transition, which traces the large scale gas kinematics. Using a friend-of-friend algorithm, we identify four main velocity coherent structures, all of which are associated with prestellar and protostellar cores. One of them presents a filament-like structure, and our observations could be consistent with mass accretion towards one of the protostars. In this case, we estimate a mass accretion rate of $ \dot{M}_\mathrm{acc}\approx 2 \times 10^{-4} \rm \, M_\odot \, yr^{-1}$. Our results support a clump-fed accretion scenario in the targeted source. The cores in prestellar stage are essentially low-mass, and they appear subvirial and gravitationally bound, unless further support is available for instance due to magnetic fields.

The bright supergiant, Betelgeuse (Alpha Orionis, HD 39801), underwent a historic optical dimming during 2020 January 27 $-$ February 13. Many imaging and spectroscopic observations across the electromagnetic spectrum were obtained prior to, during, and subsequent to this dimming event. These observations of Betelgeuse reveal that a substantial surface mass ejection (SME) occurred and moved out through the extended atmosphere of the supergiant. A photospheric shock occurred in 2019 January - March, progressed through the extended atmosphere of the star during the following 11 months and led to dust production in the atmosphere. Resulting from the substantial mass outflow, the stellar photosphere was left with lower temperatures and the chromosphere with a lower density. The mass ejected could represent a significant fraction of the total annual mass loss rate from the star suggesting that episodic mass loss events can contribute an amount comparable to that of the stellar wind. Following the SME, Betelgeuse was left with a cooler average photosphere, an unusual short photometric oscillation, reduced velocity excursions, and the disappearance of the $\sim$400-day pulsation in the optical and radial velocity for more than two years following the Great Dimming.

We introduce the scientific motivations for the development of the Cassegrain U-Band Efficient Spectrograph (CUBES) that is now in construction for the Very Large Telescope. The assembled cases span a broad range of contemporary topics across Solar System, Galactic and extragalactic astronomy, where observations are limited by the performance of current ground-based spectrographs shortwards of 400nm. A brief background to each case is presented and specific technical requirements on the instrument design that flow-down from each case are identified. These were used as inputs to the CUBES design, that will provide a factor of ten gain in efficiency for astronomical spectroscopy over 300-405nm, at resolving powers of R~24,000 and ~7,000. We include performance estimates that demonstrate the ability of CUBES to observe sources that are up to three magnitudes fainter than currently possible at ground-ultraviolet wavelengths, and we place its predicted performance in the context of existing facillities.

In this paper, we report the potential detection of a nonmonotonic radial rotation profile in a low-mass lower-luminosity giant star. For most low- and intermediate-mass stars, the rotation on the main sequence seems to be close to rigid. As these stars evolve into giants, the core contracts and the envelope expands, which should suggest a radial rotation profile with a fast core and a slower envelope and surface. KIC 9267654, however, seems to show a surface rotation rate that is faster than its bulk envelope rotation rate, in conflict with this simple angular momentum conservation argument. We improve the spectroscopic surface constraint, show that the pulsation frequencies are consistent with the previously published core and envelope rotation rates, and demonstrate that the star does not show strong chemical peculiarities. We discuss the evidence against any tidally interacting stellar companion. Finally, we discuss the possible origin of this unusual rotation profile, including the potential ingestion of a giant planet or unusual angular momentum transport by tidal inertial waves triggered by a close substellar companion, and encourage further observational and theoretical efforts.

The identification of AGN in large surveys has been hampered by seemingly discordant classifications arising from differing diagnostic methods, usually tracing distinct processes specific to a particular wavelength regime. However, as shown in Yao et al. (2020), the combination of optical emission line measurements and mid-infrared photometry can be used to optimise the discrimination capability between AGN and star formation activity. In this paper we test our new classification scheme by combining the existing GAMA-WISE data with high-quality MeerKAT radio continuum data covering 8 deg$^2$ of the GAMA G23 region. Using this sample of 1 841 galaxies (z < 0.25), we investigate the total infrared (derived from 12$\mu$m) to radio luminosity ratio, q(TIR), and its relationship to optical-infrared AGN and star-forming (SF) classifications. We find that while q(TIR) is efficient at detecting AGN activity in massive galaxies generally appearing quiescent in the infrared, it becomes less reliable for cases where the emission from star formation in the host galaxy is dominant. However, we find that the q(TIR) can identify up to 70 % more AGNs not discernible at optical and/or infrared wavelengths. The median q(TIR) of our SF sample is 2.57 $\pm$ 0.23 consistent with previous local universe estimates.

Combining constraints from microlensing and Lyman-$\alpha$ forest, we provide a simple argument to show that large spatial clustering of stellar-mass primordial black holes at the time of formation, such as the one induced by the presence of large non-Gaussianities, is ruled out. Therefore, it is not possible to evade existing constraints preventing stellar-mass primordial black holes to be a dominant constituent of the dark matter by boosting their initial clustering.

The corrected dependence of the mean depth of the EAS maximum $X_{max}$ on the energy was obtained from the data of the Tunka-133 array for 7 years and the TAIGA-HiSCORE array for 2 year. The parameter $\langle\ln A\rangle$, characterizing the mean mass compositon was derived from these results. The differential energy spectrum of primary cosmic rays in the energy range of $2\cdot 10^{14}$ - $2\cdot 10^{16}$\,eV was reconstructed using the new parameter $Q_{100}$ the Cherenkov light flux at the core distance 100 m.}

Correlated noise in exoplanet light curves, such as noises from stellar activity, convection noise, and instrumental noises distorts the exoplanet transit light curves, and leads to biases in the best-fit transit parameters. An optimal fitting algorithm is stable against the presence of correlated noises and lead to statistically consistent results, i.e. the actual biases are usually within the error interval. This is not automatically satisfied by most of the algorithms in everyday use, and the testing of the algorithms is necessary. In this paper, we describe a bootstrapping-like test to handle with the general case, and apply this to the wavelet-based TLCM (Transit and Light Curve Modeller) algorithm, testing it for the stability against the correlated noise. We contrast the results to the FITSH algorithm that is based on a white noise assumption. We simulated transit light curves with previously known parameters in the presence of a correlated noise model generated by an ARIMA (Autoregressive Integretad Moving Average) process. Then we solved the simulated observations, and examined the resulting parameters and error intervals. We have found that the assumption of FITSH that only white noise is present led to inconsistencies in the results: the distribution of best-fit parameters is by a factor of 3--6 broader then the determined error intervals. On the other hand, the wavelet-based TLCM algorithm handles the correlated noise properly, leading to properly determined parameter and error intervals which are perfectly consistent with the actual biases.

We present photometric and polarimetric measurements of gamma-ray burst (GRB) optical afterglows observed by the RINGO3 imaging polarimeter over its $\sim$7 year lifetime mounted on the Liverpool Telescope. During this time, RINGO3 responded to 67 GRB alerts. Of these, 28 had optical afterglows and a further ten were sufficiently bright for photometric and polarimetric analysis ($R\lessapprox{17}$). We present high quality multicolour light curves of ten sources: GRB 130606A, GRB 130610A, GRB 130612A, GRB 140430A, GRB 141220A, GRB 151215A, GRB 180325A, GRB 180618A, GRB 190114C, and GRB 191016A and polarimetry for seven of these (excluding GRB 130606A, GRB 130610A, and GRB 130612A, which were observed before the polarimetry mode was fully commissioned). Eight of these ten GRBs are classical long GRBs, one sits at the short-long duration interface with a $T_{90}$ $\sim$ 4 seconds and one is a classical short, hard burst with extended emission. We detect polarization for GRB 190114C and GRB 191016A. While detailed analyses of several of these GRBs have been published previously, here we present a uniform re-reduction and analysis of the whole sample and investigation of the population in a broad context relative to the current literature. We use survival analysis to fully include the polarization upper limits in the comparison with other GRB properties, such as temporal decay rate, isotropic energy and redshift. We find no clear correlation between polarization properties and wider sample properties and conclude that larger samples of early time polarimetry of GRB afterglows are required to fully understand GRB magnetic fields.

We observed the Galactic black hole X-ray binary V4641 Sgr with the high resolution transmission gratings on Chandra during the source's 2020 outburst. Over two epochs of Chandra gratings observations, we see numerous highly ionized metal lines, superimposed on a hot, disc-dominated X-ray continuum. The measured inner disc temperatures and luminosities imply an unfeasibly small inner disc radius, such that we suggest that the central engine of V4641 Sgr is obscured, and we are viewing scattered X-rays. We find that the emission lines in the Chandra spectra cannot be constrained by a single photoionized model, instead finding that two separate photoionized model components are required, one to reproduce the iron lines and a second for the other metals. We compare the observed X-ray spectra of V4641 Sgr to optical studies during previous outbursts of the source, suggesting that the lines originate in an accretion disc wind, potentially with a spherical geometry.

The future Ricochet experiment aims at searching for new physics in the electroweak sector by providing a high precision measurement of the Coherent Elastic Neutrino-Nucleus Scattering (CENNS) process down to the sub-100 eV nuclear recoil energy range. The experiment will deploy a kg-scale low-energy-threshold detector array combining Ge and Zn target crystals 8.8 meters away from the 58 MW research nuclear reactor core of the Institut Laue Langevin (ILL) in Grenoble, France. Currently, the Ricochet collaboration is characterizing the backgrounds at its future experimental site in order to optimize the experiment's shielding design. The most threatening background component, which cannot be actively rejected by particle identification, consists of keV-scale neutron-induced nuclear recoils. These initial fast neutrons are generated by the reactor core and surrounding experiments (reactogenics), and by the cosmic rays producing primary neutrons and muon-induced neutrons in the surrounding materials. In this paper, we present the Ricochet neutron background characterization using $^3$He proportional counters which exhibit a high sensitivity to thermal, epithermal and fast neutrons. We compare these measurements to the Ricochet Geant4 simulations to validate our reactogenic and cosmogenic neutron background estimations. Eventually, we present our estimated neutron background for the future Ricochet experiment and the resulting CENNS detection significance.

Changing-look active galactic nuclei (AGNs) present an important laboratory to understand the origin and physical properties of the broad-line region (BLR). We investigate follow-up optical spectroscopy spanning $\sim 500$ days after the outburst of the changing-look AGN 1ES\,1927+654. The emission lines displayed dramatic, systematic variations in intensity, velocity width, velocity shift, and symmetry. Analysis of optical spectra and multi-band images indicate that the host galaxy contains a pseudobulge and a total stellar mass of $3.56_{-0.35}^{+0.38} \times 10^{9}\, M_\odot$. Enhanced continuum radiation from the outburst produced an accretion disk wind, which condensed into BLR clouds in the region above and below the temporary eccentric disk. Broad Balmer lines emerged $\sim 100$ days after the outburst, together with an unexpected, additional component of narrow-line emission. The newly formed BLR clouds then traveled along a similar eccentric orbit ($e \approx 0.6$). The Balmer decrement of the BLR increased by a factor of $\sim 4-5$ as a result of secular changes in cloud density. The drop in density at late times allowed the production of \hei\ and \heii\ emission. The mass of the black hole cannot be derived from the broad emission lines because the BLR is not virialized. Instead, we use the stellar properties of the host galaxy to estimate $M_\mathrm{BH} = 1.38_{-0.66}^{+1.25} \times 10^{6}\, M_\odot$. The nucleus reached near or above its Eddington limit during the peak of the outburst. We discuss the nature of the changing-look AGN 1ES\,1927+654 in the context of other tidal disruption events.

Speckle Noise is the dominant source of error in high contrast imaging with adaptive optics system. We discuss the potential for wavefront sensing telemetry to calibrate speckle noise with sufficient precision and accuracy so that it can be removed in post-processing of science images acquired by high contrast imaging instruments. In such a self-calibrating system, exoplanet detection would be limited by photon noise and be significantly more robust and efficient than in current systems. We show initial laboratory and on-sky tests, demonstrating over short timescale that residual speckle noise is indeed calibrated to an accuracy exceeding readout and photon noise in the high contrast region. We discuss immplications for the design of space and ground high-contrast imaging systems.

Adaptive Optics projects at Subaru Telescope span a wide field of capabilities ranging from ground-layer adaptive optics (GLAO) providing partial correction over a 20 arcmin FOV to extreme adaptive optics (ExAO) for exoplanet imaging. We describe in this paper current and upcoming narrow field-of-view capabilities provided by the Subaru Extreme Adaptive Optics Adaptive Optics (SCExAO) system and its instrument modules, as well as the upcoming 3000-actuator upgrade of the Nasmyth AO system.

Lyman Break Galaxy (LBG) candidates at $z\gtrsim12$ are rapidly being identified in JWST/NIRCam imaging. Due to the (redshifted) break produced by neutral hydrogen absorption of rest-frame UV photons, these sources are expected to drop out in the bluer filters like F150W and F200W while being well-detected in redder filters (e.g., F277W, F356W, F444W). However, dust-enshrouded star-forming galaxies at lower redshifts ($z\lesssim7$) may also mimic the near-infrared colors of $z>12$ LBGs, representing potential contaminants in LBG candidate samples. Here, we report a galaxy, CEERS_DSFG_1, that drops out in the F115W, F150W, and F200W filters, for which a photometric redshift fit to the JWST data alone predicts a redshift of $z_{\rm phot}\sim18$. However, we show it is a dusty star-forming galaxy (DSFG) at $z\approx5$ based on deep millimeter interferometric observations conducted with NOEMA. We also present a $2.6\sigma$ SCUBA-2 detection at 850$\,\mu\rm m$ around the position of a recently reported $z\approx16.7$ LBG candidate in the same field, CEERS-93316. While we cannot conclusively show this detection is astrophysical or associated with this object, we illustrate that if it is associated, the available photometry are consistent with a $z\sim5$ DSFG with strong nebular emission lines despite its blue NIR colors. Hence, we conclude that robust (sub)millimeter detections in NIRCam dropout galaxies likely imply $z\sim4-6$ redshift solutions, where the observed near-infrared break would be the result of a strong rest-frame optical Balmer break combined with high dust attenuation and strong nebular line emission, rather than the rest-frame UV Lyman break. This provides evidence that DSFGs may contaminate searches for ultra high-redshift LBG candidates from JWST observations.

The accumulation and circulation of carbon-hydrogen dictate the chemical evolution of ice giant planets. Species separation and diamond precipitation have been reported in carbon-hydrogen systems, verified by static and shock-compression experiments. Nevertheless, the dynamic formation processes for the above-mentioned phenomena are still insufficiently understood. Here, combing deep learning model, we demonstrate that diamonds form through a three-step process involving decomposition, species separation and nucleation procedures. Under shock condition of 125 GPa and 4590 K, hydrocarbons are decomposed to give hydrogen and low-molecular-weight alkanes (CH4 and C2H6), which escape from the carbon chains resulting in C/H species separation. The remaining carbon atoms without C-H bonds accumulate and nucleate to form diamond crystals. The process of diamond growth is found to associated with a critical nucleus size where dynamic energy barrier plays a key role. These dynamic processes for diamonds formation are insightful in establishing the model for ice giant planet evolution.

We carry out a photometric search for new cataclysmic variable stars (CVs), with the goal of identification for candidates of AR~Scorpii-type binary systems. We select GAIA sources that are likely associated with unidentified X-ray sources, and analyze the light curves taken by the Zwicky Transient Facility, Transiting Exoplanet Survey Satellite, and Lulin One-meter Telescope in Taiwan. We investigate eight sources as candidates for CVs, among which six sources are new identifications. Another two sources have been recognized as CVs in previous studies, but no detailed investigations have been done. We identify two eclipsing systems that are associated with an unidentified XMM-Newton or Swift source, and one promising candidate for polar associated with an unidentified ASKA source. Two polar candidates may locate in the so-called period gap of a CV, and the other six candidates have an orbital period shorter than that of the period gap. Although we do not identify a promising candidate for AR~Scorpii-type binary systems, our study suggests that CV systems that have X-ray emission and do not show frequent outbursts may have been missed in previous surveys.

One--dimensional stellar structure and evolution programs are built using different physical prescriptions and algorithms, which means there can be variations between models' predictions even when using identical input physics. This leads to questions about whether such deviations are physical or numerical; code validation studies are important and necessary tools for studying these questions. We provide the first direct comparison between the Monash stellar evolution program and MESA for a $2M_{\odot}$ model evolved from the zero-age main sequence to the tip of the thermally pulsing asymptotic giant branch. We compare the internal structure of the two models at six critical evolutionary points and find that they are in excellent agreement in characteristics like central temperature, central density and the temperature at the base of the convective envelope during the thermally pulsing asymptotic giant branch. The hydrogen-exhausted core mass between the models differs by less than 4.2% throughout the entire evolution, the final values vary only by 1.5%. Surface quantities such as luminosity and radius vary by less than 0.2% prior to the asymptotic giant branch. During thermal pulses, the difference extends to 3.4%, largely due to uncertainties in mixing and the treatment of atmospheric boundary conditions. Given that the veteran Monash code is closed source, the present work provides the first fully open-source computational analog. This increases accessibility to precision modeling on the asymptotic giant branch and lays the groundwork for higher-mass calculations that are performed with MESA but preserve the standards of the Monash code during the AGB.

Supermassive black holes in active galactic nuclei launch relativistic jets, as indicated by observed superluminal radio blobs. The energy source of these jets is widely discussed in the theoretical framework of Blandford-Znajek process, the electromagnetic energy extraction from rotating black holes (BHs), while formation mechanism of the radio blobs in the electromagnetically-dominated jets has been a long-standing problem. Recent high-resolution magnetohydrodynamic simulations of magnetically arrested disks exhibited magnetic reconnection in a transient magnetically-dominated part of the equatorial disk near the BH horizon, which led to a promising scenario of efficient MeV gamma-ray production and subsequent electron-positron pair loading into BH magnetosphere. We develop this scenario to build a theoretical framework on energetics, timescales and particle number density of the superluminal radio blobs and discuss observable signatures in other wavebands. We analytically show that the non-thermal electrons emit broadband photons from optical to multi-MeV bands. The electron-positron pairs produced in the magnetosphere are optically thick for synchrotron-self absorption, so that the injected energy is stored in the plasma. The stored energy is enough to power the superluminal radio blobs observed in M87. This scenario predicts rather dim radio blobs around Sgr A*, which are consistent with no clear detection by current facilities. In addition, this scenario inevitably produces strong X-ray flares in a short timescale, which will be detectable by future X-ray satellites.

Molecular gas disks are generally Toomre stable ($Q_T>$1) and yet clearly gravitationally unstable to structure formation as evidenced by the existence of molecular clouds and ongoing star formation. This paper adopts a 3D perspective to obtain a general picture of instabilities in flattened rotating disks, using the 3D dispersion relation to describe how disks evolve when perturbed over their vertical extents. By explicitly adding a vertical perturbation to an unperturbed equilibrium disk, stability is shown to vary with height above the mid-plane. Near to $z$=0 where the equilibrium density is roughly constant, instability takes on a Jeans-like quality, occurring on scales larger than the Jeans length and subject to a threshold $Q_M=\kappa^2/(4\pi G\rho)=1$ or roughly $Q_T\approx 2$. Far from the mid-plane, on the other hand, stability is pervasive, and the threshold for the total disk (out to $z=\pm\infty$) to be stabilized is lowered to $Q_T=1$ as a consequence. In this new framework, gas disks are able to fragment through partial 3D instability even where total 2D instability is suppressed. The growth rates of the fragments formed via 3D instability are comparable to, or faster than, Toomre instabilities. The rich structure in molecular disks on the scale of 10s of pc can thus be viewed as a natural consequence of their 3D nature and their exposure to a variety of vertical perturbations acting on roughly a disk scale height, i.e. due to their situation within the more extended galaxy potential, participation in the disk-halo flow, and exposure to star formation feedback.

Study of pre-stellar cloud evolution requires observations with high sensitivity and resolution, and regions of high-mass star formation are particularly challenging. We wish to quantify, to what accuracy the physical conditions within a massive star-forming cloud can be determined from observations. We are particularly interested in the possibilities offered by the Next Generation VLA (ngVLA) interferometer. We use data from a magnetohydrodynamic simulation of star-formation and concentrate on a filamentary structure that has physical properties similar to an infrared-dark cloud. We produce synthetic ngVLA observations of spectral lines and analyse the column density, gas temperature, and kinematics. The results are compared to ideal observations and the actual 3D model. For a distance of 4 kpc, ngVLA provides a resolution of 0.01 pc even in its most compact configuration. For abundant molecules, such as HCO+, NH3, N2H+, and CO isotopomers, cloud kinematics and structure can be mapped down to sub-arcsec scales. For NH3, a reliable column density map could be obtained for the entire 15 * 40 arcsec cloud, even without additional single-dish data, and kinetic temperatures are recovered to a precision of 1 K. At higher frequencies, the loss of large-scale emission is noticeable. The line observations accurately trace the cloud kinematics, except for the largest scales. The line-of-sight confusion complicates the interpretation of the kinematics, and limits the usefulness of collapse indicators based on the blue asymmetry of optically thick lines. The ngVLA will provide accurate data on the small-scale structure and the physical and chemical state of clouds, even in high-mass star-forming regions at kiloparsec distances. Complementary single-dish data are still essential for estimates of the total column density and the large-scale kinematics.

We explore the growth of planetary embryos by planetesimal accretion up to and beyond the point where pebble accretion becomes efficient at the so-called Hill-transition mass. Both the transition mass and the characteristic mass of planetesimals formed by the streaming instability increase with increasing distance from the star. We developed a model for the growth of a large planetesimal (embryo) embedded in a population of smaller planetesimals formed in a filament by the streaming instability. The model includes in a self-consistent way the collisional mass growth of the embryo, the fragmentation of the planetesimals, the velocity evolution of all involved bodies, as well as the viscous spreading of the filament. We find that the embryo accretes all available material in the filament during the lifetime of the protoplanetary disc only in the inner regions of the disc. In contrast, we find little or no growth in the outer parts of the disc beyond 5--10 AU. Overall, our results demonstrate very long timescales for collisional growth of planetesimals in the regions of the protoplanetary disc where giant planets form. As such, in order to form giant planets in cold orbits, pebble accretion must act directly on the largest bodies present in the initial mass-function of planetesimals with little or no help from mutual collisions.

IceCube is a cubic-kilometer Cherenkov detector in the deep ice at the geographic South Pole. The dominant event yield in the deep ice detector consists of penetrating atmospheric muons with energies above approximately 300 GeV, produced in cosmic ray air showers. In addition, the surface array, IceTop, measures the electromagnetic component and GeV muons of air showers. Hence, IceCube and IceTop yield unique opportunities to study cosmic rays with unprecedented statistics in great detail. We will present recent results of comic ray measurements from IceCube and IceTop. In this overview, we will highlight measurements of the energy spectrum of cosmic rays from 250 TeV up to the EeV range and their mass composition above 3 PeV. We will also report recent results from measurements of the muon content in air showers and discuss their consistency with predictions from current hadronic interaction models.

It is generally assumed that galaxies are a bimodal population in both star formation and structure: star-forming galaxies are disks, while passive galaxies host large bulges or are entirely spheroidal. Here, we test this scenario by presenting a full census of the kinematic morphologies of a volume-limited sample of galaxies in the local Universe extracted from the MaNGA galaxy survey. We measure the integrated stellar line-of-sight velocity to velocity dispersion ratio ($V/\sigma$) for 4574 galaxies in the stellar mass range $9.75 < \log M_{\star}[\rm{M}_{\odot}] < 11.75$. We show that at fixed stellar mass, the distribution of $V/\sigma$ is not bimodal, and that a simple separation between fast and slow rotators is over-simplistic. Fast rotators are a mixture of at least two populations, referred to here as dynamically-cold disks and intermediate systems, with disks dominating in both total stellar mass and number. When considering star-forming and passive galaxies separately, the star-forming population is almost entirely made up of disks, while the passive population is mixed, implying an array of quenching mechanisms. Passive disks represent $\sim$30% (both in number and mass) of passive galaxies, nearly a factor of two higher than that of slow rotators, reiterating that these are an important population for understanding galaxy quenching. These results paint a picture of a local Universe dominated by disky galaxies, most of which become somewhat less rotation-supported upon or after quenching. While spheroids are present to a degree, they are certainly not the evolutionary end-point for the majority of galaxies.

The Mars Microphone is one of the five measurement techniques of SuperCam, an improved version of the ChemCam instrument that has been functioning aboard the Curiosity rover for several years. SuperCam is located on the Rover's Mast Unit, to take advantage of the unique pointing capabilities of the rover's head. In addition to being the first instrument to record sounds on Mars, the SuperCam Microphone can address several original scientific objectives: the study of sound associated with laser impacts on Martian rocks to better understand their mechanical properties, the improvement of our knowledge of atmospheric phenomena at the surface of Mars: atmospheric turbulence, convective vortices, dust lifting processes and wind interactions with the rover itself. The microphone will also help our understanding of the sound signature of the different movements of the rover: operations of the robotic arm and the mast, driving on the rough floor of Mars, monitoring of the pumps, etc ... The SuperCam Microphone was delivered to the SuperCam team in early 2019 and integrated at the Jet Propulsion Laboratory (JPL, Pasadena, CA) with the complete SuperCam instrument. The Mars 2020 Mission launched in July 2020 and landed on Mars on February 18th 2021. The mission operations are expected to last until at least August 2023. The microphone is operating perfectly.

Using high angular resolution ALMA observations ($0.02^{\prime\prime}\approx0.34$ pc), we study the thermal structure and kinematics of the proto super star cluster $13$ in the central region of NGC253 through their continuum and vibrationally excited HC$_3$N emission from $J=24-23$ and $J=26-25$ lines arising from vibrational states up to $v_4=1$. We have carried 2D-LTE and non-local radiative transfer modelling of the radial profile of the HC$_3$N and continuum emission in concentric rings of $0.1$ pc width. From the 2D-LTE analysis, we found a Super Hot Core (SHC) of $1.5$ pc with very high vibrational temperatures ($>500$ K), and a jump in the radial velocity ($21$ km s$^{-1}$) in the SE-NW direction. From the non-local models, we derive the HC$_3$N column density, H$_2$ density and dust temperature ($T_\text{dust}$) profiles. Our results show that the thermal structure of the SHC is dominated by the greenhouse effect due to the high dust opacity in the IR, leading to an overestimation of the LTE $T_\text{dust}$ and its derived luminosity. The kinematics and $T_\text{dust}$ profile of the SHC suggest that star formation was likely triggered by a cloud-cloud collision. We compare proto-SSC $13$ to other deeply embedded star-forming regions, and discuss the origin of the $L_\text{IR}/M_{\text{H}_2}$ excess above $\sim100$ L$_\odot$ M$_\odot^{-1}$ observed in (U)LIRGs.

The earthward journey of ultra high energy electrons ($\sim 600$ TeV) produced in the Pulsar atmosphere by Landau damping of magneto-centrifugally excited Langmuir waves (drawing energy form the rotational slowdown) on primary electrons, is charted. It is shown, that just as they escape the light cylinder zone, the ultra-high energy particles, interacting with the medium of the Crab nebula, rapidly loose their energy via the quantum synchrotron process, producing highly energetic gamma rays ~ $\sim 0.6$PeV. Interacting with the cosmic background radiation in the interstellar medium, only a tiny fraction of these ultra high energy photons (via the $\gamma\gamma$ channel) are, then transformed into electron-positron pairs. Detected flux of these photons imposes an upper limit on the fraction ($4\times 10^{-7}$) of the magnetospheric particles involved in the process of generation of ultra-high energy photons (up to $600$ TeV).

Context. The study of the oldest and most metal-poor stars in our Galaxy promotes our understanding of the Galactic chemical evolution and the beginning of Galaxy and star formation. However, they are notoriously difficult to find, with only five stars at $\mathrm{[Fe/H]<-5.0}$ having been detected to date. Thus, the spectrophotometric data of 219 million sources which became available in the third Gaia Data Release comprise a very promising dataset for the identification of metal-poor stars. Aims. We want to use the low-resolution Gaia Blue Photometer / Red Photometer (BP/RP) spectra to identify metal-poor stars. Our primary aspiration is to help populate the poorly constrained tail of the metallicity distribution function of the stellar halo of the Galaxy. Methods. We developed a metal-poor candidate selection method based on flux ratios from the BP/RP Gaia spectra, using simulated synthetic spectra. Results. We found a relation between the relative iron abundance and the flux ratio of the Ca H \& K region to that of the $\mathrm{H\beta}$ line. This relation is temperature and surface gravity dependent, and it holds for stars with $\mathrm{4800\,K \leq T_{eff}\leq6300\,K}$. We applied it to noisy simulated synthetic spectra and inferred $\mathrm{[Fe/H]}$ with an uncertainty of $\sigma_{\mathrm{[Fe/H]}}\lessapprox0.65$ dex for $\mathrm{-3\leq[Fe/H]}\leq 0.5$ and G=15-17mag, which is sufficient to identify stars at $\mathrm{[Fe/H]<-2.0 }$ reliably. We predict that by selecting stars with inferred $\mathrm{[Fe/H]}\leq-2.5$ dex, we can retrieve 80% of the stars with $\mathrm{[Fe/H]}\leq-3$ and have a success rate of about 50%, that is one in two stars we select would have $\mathrm{[Fe/H]}\leq-3$. We do not take into account the effect of reddening, so our method should only be applied to stars which are located in regions of low extinction.

Missions to small celestial bodies rely heavily on optical feature tracking for characterization of and relative navigation around the target body. While deep learning has led to great advancements in feature detection and description, training and validating data-driven models for space applications is challenging due to the limited availability of large-scale, annotated datasets. This paper introduces AstroVision, a large-scale dataset comprised of 115,970 densely annotated, real images of 16 different small bodies captured during past and ongoing missions. We leverage AstroVision to develop a set of standardized benchmarks and conduct an exhaustive evaluation of both handcrafted and data-driven feature detection and description methods. Next, we employ AstroVision for end-to-end training of a state-of-the-art, deep feature detection and description network and demonstrate improved performance on multiple benchmarks. The full benchmarking pipeline and the dataset will be made publicly available to facilitate the advancement of computer vision algorithms for space applications.

The Magellanic Clouds, two dwarf galaxy companions to the Milky Way, are among the Fermi Large Area Telescope (LAT) brightest gamma-ray sources. Aiming at a comprehensive modeling of the non-thermal electromagnetic and neutrino emission in both Clouds, we self-consistently model the radio and gamma-ray spectral energy distribution from their disks based on recently published Murchison Widefield Array and Fermi/LAT data. All relevant radiative processes involving relativistic and thermal electrons (synchrotron, Compton scattering, and bremsstrahlung) and relativistic protons (neutral-pion decay following interaction with thermal protons) are considered, using exact emission formulae. Our joint spectral analyses indicate that radio emission in the Clouds has both primary and secondary electron synchrotron and thermal bremsstrahlung origin, whereas gamma rays originate mostly from neutral-pion decay with some contributions from relativistic bremsstrahlung and Compton scattering off starlight. The proton spectra in both galaxies are modeled as power laws in energy with similar spectral indices, ~2.4, and energy densities, ~1 eV/cm3. The predicted 0.1-10 GeV neutrino flux is too low for detection by current and upcoming experiments. Our analyses confirm earlier suggestions of a largely hadronic origin of the gamma-ray emission in both Magellanic Clouds.

Ground-based optical photometry of the counterparts of High-Mass X-ray Binaries (HMXBs) has revealed the presence of periodic modulations on timescales of ~0.3-0.5 d. More recent space-based observations Corot and TESS of OB and Be stars have shown that pulsations caused by p and g modes are common in early type stars. We have therefore undertaken a systematic search for variability in the optical counterparts of 23 HMXBs (mostly neutron star systems, but including one black hole, Cyg X-1) using TESS data primarily in 2 min cadence mode. After removing the orbital period modulation in four systems, we find that all 23 sources show evidence for quasi-periodic variability on periods shorter than ~1 d. We compare their power spectra with those from observations of other OB and Be type stars. In two systems, V725 Tau and HD 249179 (which may not be a HMXB), we find evidence for an outburst, the former being simultaneous with an X-ray flare. We search for changes in the power spectra over the outburst duration, and compare them with outbursts seen in other Be systems.

We investigate the origin in the early Solar System of the short-lived radionuclide 244Pu (with a half life of 80 Myr) produced by the rapid (r) neutron-capture process. We consider two large sets of r-process nucleosynthesis models and analyse if the origin of 244Pu in the ESS is consistent with that of the other r and slow (s) neutron-capture process radioactive nuclei. Uncertainties on the r-process models come from both the nuclear physics input and the astrophysical site. The former strongly affects the ratios of isotopes of close mass (129I/127I, 244Pu/238U, and 247Pu/235U). The 129I/247Cm ratio, instead, which involves isotopes of a very different mass, is much more variable than those listed above and is more affected by the physics of the astrophysical site. We consider possible scenarios for the evolution of the abundances of these radioactive nuclei in the galactic interstellar medium and verify under which scenarios and conditions solutions can be found for the origin of 244Pu that are consistent with the origin of the other isotopes. Solutions are generally found for all the possible different regimes controlled by the interval ($\delta$) between additions from the source to the parcel of interstellar medium gas that ended up in the Solar System, relative to decay timescales. If r-process ejecta in interstellar medium are mixed within a relatively small area (leading to a long $\delta$), we derive that the last event that explains the 129I and 247Cm abundances in the early Solar System can also account for the abundance of 244Pu. Due to its longer half life, however, 244Pu may have originated from a few events instead of one only. If r-process ejecta in interstellar medium are mixed within a relatively large area (leading to a short $\delta$), we derive that the time elapsed from the formation of the molecular cloud to the formation of the Sun was 9-16 Myr.

Giant radio sources (GRSs) defined to be > 0.7 Mpc are the largest single objects in the Universe and can be associated with both galaxies (GRGs) and quasars (GRQs). They are important for understanding the evolution of radio galaxies and quasars whose sizes range from pc to Mpc scales and are also valuable probes of their environment. These radio-loud active galactic nuclei (RLAGN) interact with the interstellar medium of the host galaxy on small scales and the large-scale intracluster or intergalactic medium for the GRSs. With several new and sensitive surveys over the last few years, the number of known GRSs has increased many fold which has led a resurgence of interest in the field. This review article summarises our current understanding of these sources based on nearly five decades of research, and discusses the importance of the Square Kilometer Array (SKA) in addressing some of the outstanding questions.

Context. Weak lensing and clustering statistics beyond two-point functions can capture non-Gaussian information about the matter density field, thereby improving the constraints on cosmological parameters relative to the mainstream methods based on correlation functions and power spectra. Aims. This paper presents a cosmological analysis of the fourth data release of the Kilo Degree Survey based on the density split statistics, which measures the mean shear profiles around regions classified according to foreground densities. The latter is constructed from a bright galaxy sample, which we further split into red and blue samples, allowing us to probe their respective connection to the underlying dark matter density. Methods. We use the state-of-the-art model of the density splitting statistics and validate its robustness against mock data infused with known systematic effects such as intrinsic galaxy alignment and baryonic feedback. Results. After marginalising over the photometric redshift uncertainty and the residual shear calibration bias, we measure for the full KiDS-bright sample a structure growth parameter of $S_8 = \sigma_8 \sqrt{\Omega_\mathrm{m}/0.3} = 0.74^{+0.03}_{-0.02}$ that is competitive to and consistent with two-point cosmic shear results, a matter density of $\Omega_\mathrm{m} = 0.28 \pm 0.02$, and a constant galaxy bias of $b = 1.32^{+0.12}_{-0.10}$.

We present the details of a simulation suite for modeling the effects of readout with SLAC Microresonator RF (SMuRF) electronics. The SMuRF electronics are a warm readout and control system for use with superconducting microwave resonator-based detector systems. The system has been used with the BICEP/Keck program and will be used on the upcoming Simons Observatory and BICEP Array experiments. This simulation suite is a software implementation of the main SMuRF algorithms for offline analysis, modeling, and study. The firmware-implemented algorithms for calibration, resonator frequency estimation, and tone tracking present sources of potential bias or errors if not modeled properly. The simulator takes as input true detector signal, realistic resonator properties, and SMuRF-related user-controlled readout settings. It returns the final flux ramp-demodulated output of a detector timestream as would be passed to the experiment data acquisition system, enabling the analysis of the impact of readout-related parameters on the final science data. It is publicly available in Python with accompanying Jupyter notebooks for user tutorials.

Using data from JWST, we analyze the compact sources ("sparkles") located around a remarkable $z_{\rm spec}=1.378$ galaxy (the "Sparkler") that is strongly gravitationally lensed by the $z=0.39$ galaxy cluster SMACS J0723.3-7327. Several of these compact sources can be cross-identified in multiple images, making it clear that they are associated with the host galaxy. Combining data from JWST's {\em Near-Infrared Camera} (NIRCam) with archival data from the {\em Hubble Space Telescope} (HST), we perform 0.4-4.4$\mu$m photometry on these objects, finding several of them to be very red and consistent with the colors of quenched, old stellar systems. Morphological fits confirm that these red sources are spatially unresolved even in strongly magnified JWST/NIRCam images, while JWST/NIRISS spectra show [OIII]5007 emission in the body of the Sparkler but no indication of star formation in the red compact sparkles. The most natural interpretation of these compact red companions to the Sparkler is that they are evolved globular clusters seen at $z=1.378$. Applying \textsc{Dense Basis} SED-fitting to the sample, we infer formation redshifts of $z_{form} \sim 7-11$ for these globular cluster candidates, corresponding to ages of $\sim 3.9-4.1$ Gyr at the epoch of observation and a formation time just $\sim$0.5~Gyr after the Big Bang. If confirmed with additional spectroscopy, these red, compact "sparkles" represent the first evolved globular clusters found at high redshift, could be amongst the earliest observed objects to have quenched their star formation in the Universe, and may open a new window into understanding globular cluster formation. Data and code to reproduce our results will be made available at \faGithub\href{https://niriss.github.io/sparkler.html}{this http URL}.

Significant absorption of radiation is usually accompanied by refraction. This is not the case for $\gamma$ rays travelling cosmic distances. We show that the real and imaginary parts of the refraction index are indeed commensurable, as they are related by dispersion relations, but when turning to physical observables, the (finite) optical depth is way larger than the (infinitesimal) time delay of the gamma rays relative to gravitational radiation. The numerically large factor solving the apparent contradiction is $E_\gamma/H_0$ arising from basic wave properties (Bouguer-Beer-Lambert law) and the standard cosmological model, respectively.

The most distant known trans-Neptunian objects (perihelion distance above 38 au and semimajor axis above 150 au) are of interest for their potential to reveal past, external, or present but unseen perturbers. Realizing this potential requires understanding how the known planets influence their orbital dynamics. We use a recently-developed Poincare mapping approach for orbital phase space studies of the circular planar restricted three body problem, which we have extended to the case of the three-dimensional restricted problem with $N$ planetary perturbers. With this approach, we explore the dynamical landscape of the 23 most distant TNOs under the perturbations of the known giant planets. We find that, counter to common expectations, almost none of these TNOs are far removed from Neptune's resonances. Nearly half (11) of these TNOs have orbits consistent with stable libration in Neptune's resonances; in particular, the orbits of TNOs 148209 and 474640 overlap with Neptune's 20:1 and 36:1 resonances, respectively. Five objects can be ruled currently non-resonant, despite their large orbital uncertainties, because our mapping approach determines the resonance boundaries in angular phase space in addition to semimajor axis. Only three objects are in orbital regions not appreciably affected by resonances: Sedna, 2012 VP113 and 2015 KG163. Our analysis also demonstrates that Neptune's resonances impart a modest (few percent) non-uniformity in the longitude of perihelion distribution of the currently observable distant TNOs. While not large enough to explain the observed clustering, this small dynamical sculpting of the perihelion longitudes could become relevant for future, larger TNO datasets.

Building upon several recent advances in the development of effective-one-body models for spin-aligned eccentric binaries with individual masses $(m_1,m_2)$ we introduce a new EOB waveform model that aims at describing inspiralling binaries in the large mass-ratio regime, $m_1\gg m_2$. The model exploits the current state-of-the-art TEOBResumS-Dali model for eccentric binaries, but the standard EOB potentials $(A,\bar{D},Q)$, informed by Numerical Relativity (NR) simulations, are replaced with the corresponding functions that are linear in the symmetric mass ratio $\nu\equiv m_1 m_2/(m_1+m_2)^2$ taken at 8.5PN accuracy. To improve their strong-field behavior, these functions are: (i) suitably factorized and resummed using Pad\'e approximants and (ii) additionally effectively informed to state-of-the-art numerical results obtained by gravitational self-force theory (GSF). For simplicity, the spin-sector of the model is taken to be the one of TEOBResumS-Dali, though removing the NR-informed spin-orbit effective corrections. We propose the current GSF-informed EOB framework as a conceptually complete analytical tool to generate waveforms for eccentric Extreme (and Intermediate) Mass Ratio Inspirals for future gravitational wave detectors.

Gravitational waves might help resolve the tension between early and late Universe measurements of the Hubble constant, and this possibility can be enhanced with a gravitational wave detector in the decihertz band as we will demonstrate in this study. Such a detector is particularly suitable for the multiband observation of stellar-mass black hole binaries between space and ground, which would significantly improve the source localization accuracy thanks to a long baseline for timing triangulation, hence promoting the "dark siren" cosmology. Proposed decihertz concepts include DECIGO/B-DECIGO, TianGO, and others. We consider here the prospects of multiband observation of dark siren binaries with a variety of network configurations. We find that a multiband observation can uniquely identify a black hole binary to a single galaxy to a cosmological distance, and thus a dark siren behaves as if it had an electromagnetic counterpart. Considering only fully localized dark sirens, we use a Fisher matrix approach to estimate the error in the Hubble constant and matter density parameter. We find that a decihertz detector substantially improves our ability to measure cosmological parameters because it enables host galaxies to be identified out to a larger distance without the systematics from statistical techniques based on comparing the population distribution.

After presenting the main characteristics of the Sky of Salamanca, we analyse whether what is represented in it is motivated by astrological or astronomical considerations. We will see that the astrological explanation based on the system of planetary houses described in Claudius Ptolemy's Tetrabiblos does not correspond to the planetary configuration we see in the salmantine vault. Finally, we will conclude that the astronomical interpretation of the Sky of Salamanca is justified by the representation of the Universe known in the 15th century according to the Almagest, by the placement of the stars in the constellations according to that work and by the dating of the painted planetary configuration, which circumstantial evidence places in August 1475.

The detection of eccentricity from a gravitational wave signal is expected to help distinguish between formation channels for a given binary. In this study, we reassess all previously-reported binary black holes with previous claims of possible eccentricity as well as a few binaries with more interesting source parameters, for the first time using a model (TEOBResumSGeneral) which accounts for the full eccentricity range possible and incorporates higher-order gravitational emission critical to model emission from highly eccentric orbits. We estimate the eccentricity of these five events. For the first time, we present marginal evidence of eccentricity for one of the events: GW190929. Contrary to previous work with different settings, we do not find evidence supporting eccentric orbits for the same systems. We find the incorporation of eccentricity in our analyses dramatically shifts the posterior in multiple parameters for several events, features could negatively impact other analyses.

The focal plane of Athena's WFI consists of spectroscopic single photon X-ray detectors that contain arrays of DEPFETs (DEpleted P-channel Field-Effect Transistor) as well as ASICs that are used for steering, readout and analog signal shaping. These components have to be examined regarding the effect of ionizing radiation. A Total Ionizing Dose (TID) test was done with prototype detector modules with 64x64 DEPFETs and one SWITCHER and VERITAS ASIC each. The current design of the WFI detector head features a proton shield equivalent to 4 cm of aluminum in order to prevent a strong increase of leakage current in the fully depleted 450 $\mu$m thick bulk of the sensor. This keeps the expected doses and dose rates during the nominal mission relatively low ($\sim$5 Gy). It is nevertheless important to study the current system in a dedicated TID test in order to exclude unforeseen effects and to study any radiation related changes that can have an effect on the very sensitive readout chain and the detector performance. The combination of low doses, low dose rates, low operating temperature (<-60{\deg}C) but high sensitivity on small changes of the threshold voltages represent somehow unusual boundary conditions in comparison to TID tests for standard radiation hard electronic components. Under these circumstances it was found beneficial to do the test in our own laboratory with an X-ray source in order to realize irradiation during nominal operation conditions. Furthermore, it facilitated to take annealing effects into account. Reasonably accurate dosimetry is achieved by measuring the X-ray spectrum and intensity with the device under test. After irradiation to a total dose of 14 Gy and subsequent annealing the threshold voltage of the DEPFETs were shifted by a mean value of 80 mV, the performance remained unchanged apart from a slight increase in readout noise by 10%.

The non-relativistic cross section from Rayleigh scattering by hydrogen atoms in the ground state is calculated over a wide range of photon energies ($< 0.8$ keV). Evaluations are performed in terms of the real and imaginary components of the atomic polarizability. The sum over intermediate states that characterizes this second-order radiative process is performed using exact analytic expressions for oscillator strengths of bound and continuum states. Damping terms associated with the finite lifetimes of excited states and their splitting into two fine-structure levels ($p_{1/2}$ and $p_{3/2}$) are taken into account in resonance cross sections. Fitting formulas required for cross-section evaluation are presented for incident photon energy i) redward of the first resonance (Lyman-$\alpha_{1/2}$), ii) in the spectral region corresponding to resonances (for an arbitrary number of them), and iii) above the ionization threshold.