Papers for Tuesday, Aug 02 2022

We present a catalog of about 25,000 images of massive ($M_{\star} \ge 10^9 M_{\odot}$) galaxies at redshift $3 \leq z \leq 6$ from the TNG50 cosmological simulation, tailored for observations at multiple wavelengths carried out with the James Webb Space Telescope (JWST). The synthetic images were created with the SKIRT radiative transfer code, including the effects of dust attenuation and scattering. The noiseless images were processed with the mirage simulator to mimic the Near Infrared Camera (NIRCam) observational strategy (e.g., noise, dithering pattern, etc.) of the Cosmic Evolution Early Release Science survey (CEERS). In this paper, we analyze the predictions of the TNG50 simulation for the size evolution of galaxies at $3 \leq z \leq 6$ and the expectations for CEERS to probe that evolution. In particular, we investigate how sizes depend on wavelength, redshift, mass, and angular resolution of the images. We find that the effective radius accurately describes the three-dimensional half-mass radius of TNG50 galaxies. Sizes observed at 2 $\mu$m are consistent with those measured at 3.56 $\mu$m at all redshifts and masses. At all masses, the population of higher-$z$ galaxies is more compact than their lower-$z$ counterparts. However, the intrinsic simulation properties and their mock-observed equivalents diverge in describing the most massive galaxies, especially at $z \lesssim 4$. This discrepancy between the mass and light distribution may point to a transition in the galaxy morphology at $z$=4-5, where massive compact systems start to develop more extended stellar structures.

Supermassive stars forming at $z \sim$ 15 - 20 are one of the leading contenders for the origin of the first quasars, over 200 of which have now been discovered at $z >$ 6. These stars likely form in pristine, atomically cooled haloes immersed in strong Lyman-Werner UV backgrounds or in highly supersonic baryon streaming flows. Atomic cooling triggers catastrophic baryon collapse capable of building up stars at rates of up to $\sim$1 M$_{\odot}$ yr$^{-1}$. Here we examine the evolution of supermassive stars with a much larger and finer grid of accretion rates than in previous studies with the MESA stellar evolution code. We find that their final masses range from 3.5 $\times$ 10$^3$ M$_{\odot}$ - 3.7 $\times$ 10$^5$ M$_{\odot}$ at accretion rates of 0.001 M$_{\odot}$ yr$^{-1}$ - 1 M$_{\odot}$ yr$^{-1}$, respectively. We also find that supermassive star evolution diverges at accretion rates of 0.01 M$_{\odot}$ yr$^{-1}$ - 0.02 M$_{\odot}$ yr$^{-1}$, above which they evolve as cool red hypergiants along the Hayashi track and collapse via the general relativistic instability during central hydrogen burning, and below which they evolve as hot blue supergiants and collapse at the end of their nuclear burning lifetimes after exiting the main sequence.

The relative roles of mergers and star formation in regulating galaxy growth are still a matter of intense debate. We here present our DREAM, a new DiscRete statistical sEmi-empiricAl Model specifically designed to predict rapidly and efficiently, in a full cosmological context, galaxy assembly and merger histories for any given input stellar mass-halo mass (SMHM) relation. DREAM generates object-by-object dark matter merger trees (hence discrete) from accurate subhalo mass and infall redshift probability functions (hence statistical) for all subhaloes, including those residing within other subhaloes, with virtually no resolution limits on mass or volume. Merger trees are then converted into galaxy assembly histories via an input, redshift dependent SMHM relation, which is highly sensitive to the significant systematics in the galaxy stellar mass function and on its evolution with cosmic time. DREAM can accurately reproduce the predicted mean galaxy merger rates and assembly histories of hydrodynamic simulations and semi-analytic models, when adopting in input their SMHM relations. In the present work we use DREAM to prove that only specific SMHM relations, namely those implied by stellar mass functions characterized by large abundances of massive galaxies and significant redshift evolution, can simultaneously reproduce the local abundances of satellite galaxies, the galaxy (major merger) pairs since $z \sim 3$, and the growth of Brightest Cluster Galaxies. The same models can also reproduce the local fraction of elliptical galaxies, on the assumption that these are strictly formed by major mergers, but not the full bulge-to-disc ratio distributions, which require additional processes.

Having a massive moon has been considered as a primary mechanism for stabilized planetary obliquity, an example of which being our Earth. This is, however, not always consistent with the exoplanetary cases. This article details the discovery of an alternative mechanism, namely that planets orbiting around binary stars tend to have low spin-axis variations. This is because the large quadrupole potential of the stellar binary could speed up the planetary orbital precession, and detune the system out of secular spin-orbit resonances. Consequently, habitable zone planets around the stellar binaries in low inclination orbits hold higher potential for regular seasonal changes comparing to their single star analogues.

Kepler-289 is a three-planet system containing two sub-Neptunes and one cool giant planet orbiting a young, Sun-like star. All three planets exhibit transit timing variations (TTVs), with both adjacent planet pairs having orbital periods close to the 2:1 orbital resonance. We observe two transits of Kepler-289c with the Wide-field InfraRed Camera (WIRC) on the 200" Hale Telescope at Palomar Observatory, using diffuser-assisted photometry to achieve space-like photometric precision from the ground. These new transit observations extend the original four-year Kepler TTV baseline by an additional 7.5 years. We re-reduce the archival Kepler data with an improved stellar activity correction and carry out a joint fit with the Palomar data to constrain the transit shapes and derive updated transit times. We then model the TTVs to determine the masses of the three planets and constrain their densities and bulk compositions. Our new analysis improves on previous mass and density constraints by a factor of two or more for all three planets, with the innermost planet showing the largest improvement. Our updated atmospheric mass fractions for the inner two planets indicate that they likely have hydrogen-rich envelopes, consistent with their location on the upper side of the radius valley. We also constrain the heavy element composition of the outer saturn-mass planet, Kepler-289c, for the first time, finding that it contains 30.5 $\pm$ 6.9 $M_{\oplus}$ of metals. We use dust evolution models to show that Kepler-289c must have formed beyond 1~au, and likely beyond 3~au, and then migrated inward.

Type-Ia supernovae (SN Ia) are powerful stellar explosions that provide important distance indicators in cosmology. Recently, we proposed a new SN Ia mechanism that involves a nuclear fission chain reaction in an isolated white dwarf (WD) [PRL 126, 1311010]. The first solids that form as a WD starts to freeze are actinide rich and potentially support a fission chain reaction. In this letter we explore thermonuclear ignition from fission heating. We perform thermal diffusion simulations and find at high densities, above about 7x10^8 g/cm^3, that the fission heating can ignite carbon burning. This could produce a SN Ia or another kind of astrophysical transient.

We report discovery of Rhea 468, a 560 Myr-old moving group of stars with 50 members and candidate members at distances 2-50 pc from the Sun using an unsupervised clustering analysis of nearby stars with Gaia DR3 data. This new moving group includes the nearest brown dwarf WISE J104915.57-531906.1 AB (Luhman 16 AB) at a distance of 2 pc, which was previously suspected to be young (600-800 Myr) based on a comparison of its dynamical mass measurements with brown dwarf evolutionary models. Using stellar activity and gyrochronology, we estimate an age of $560 \pm 60$ Myr for this new ensemble of stars. This newly discovered group will be useful to refine the age and chemical composition of Luhman 16 AB, which is already one of the best substellar benchmarks known to date. Furthermore, Rhea 468 is one of the nearest young moving groups identified to date, making it a valuable laboratory for the study of exoplanets and substellar members, with 8 brown dwarf candidate members already identified here.

Elements with weak and blended spectral features in stellar spectra are challenging to measure and require specialized analysis methods to precisely measure their chemical abundances. In this work, we have created a catalog of approximately 120,000 giants with high signal-to-noise APOGEE DR17 spectra, for which we explore weak and blended species to measure Na, P, S, V, Cu, Ce, and Nd abundances and $^{12}$C/$^{13}$C isotopic ratios. We employ an updated version of the BACCHUS (Brussels Automatic Code for Characterizing High accUracy Spectra) code to derive these abundances using the stellar parameters measured by APOGEE's DR17 ASPCAP pipeline, quality flagging to identify suspect spectral lines, and a prescription for upper limits. Combined these allow us to provide our BACCHUS Analysis of Weak Lines in APOGEE Spectra (BAWLAS) catalog of precise chemical abundances for these weak and blended species that agrees well with literature and improves upon APOGEE abundances for these elements, some of which are unable to be measured with APOGEE's current, grid-based approach without computationally expensive expansions. This new catalog can be used alongside APOGEE and provide measurements for many scientific applications ranging from nuclear physics to Galactic chemical evolution and Milky Way population studies. To illustrate this we show some examples of uses for this catalog, such as, showing that we observe stars with enhanced s-process abundances or that we can use the our $^{12}$C/$^{13}$C ratios to explore extra mixing along the red giant branch.

We analyze the physical content of squeezed bispectra involving long-wavelength tensor perturbations, showing that these modes cannot be gauged away, except for the exact (unphysical) limit of infinite wavelength, $k = 0$. This result has a direct implication on the validity of the Maldacena consistency relation, respected by a subclass of inflationary models. Consequently, in the squeezed limit, as in the case of the scalar-scalar-scalar bispectrum, squeezed mixed correlators could be observed by future experiments, remaining a key channel to study Early Universe physics and discriminate among different models of inflation.

High energy gamma-rays propagating in external magnetic fields may convert into axion-like particles (ALPs). The observed gamma-ray spectra are modified by the energy-dependent conversion probability. We use the energy spectra of 20 extra-galactic gamma-ray sources recorded during 10 years of Fermi-LAT observations. We define a test statistics based upon the likelihood ratio to test the hypothesis for a spectral model without vs. a model with photon-ALPs coupling. The conversion probability is calculated for fixed values of the mass and two-photon coupling of the pseudo-scalar particle while the external magnetic field is characterized by the additional free parameters length scale $s$ and average field strength $B$. As a consistency check and in order to extend the analysis to include very high energy gamma-ray data, another test statistics is defined with the $\chi^2$ method. We find for 20 of the 20 sources a favorable fit, with the hypothesis of photon-ALPs coupling in both likelihood and $\chi^2$ analysis. The test statistics of the sources are combined and estimated to correspond to a significance of $2.7\sigma$ (local maxima) and $4.9\sigma$ (global maxima). The estimations of the combined significance are done in such a way where we simulate the individual observations under the null hypothesis to create mock data sets for each source and use bootstrap approach to increase the combinatorial statistics. Particularly, when combining the time-averaged spectrum of PKS~2155-304 from HESS and LAT, the significance increases to about $8\sigma$, whereas no comparable significant improvement is observed in the contemporaneous data set from the same instruments. The locally best-fitting values of $B$ and $s$ fall into the range that is expected for large scale magnetic fields present in the intra-cluster medium of galaxy clusters and in large scale filaments.

The earliest confirmed interstellar object, `Oumuamua, was discovered in the Solar System by Pan-STARRS in 2017, allowing for a calibration of the abundance of interstellar objects of its size $\sim 100\;$ m. This was followed by the discovery of Borisov, which allowed for a similar calibration of its size $\sim 0.4 - 1 \mathrm{\; km}$. One would expect a much higher abundance of significantly smaller interstellar objects, with some of them colliding with Earth frequently enough to be noticeable. Based on the CNEOS catalog of bolide events, we identified in 2019 the meteor detected at 2014-01-08 17:05:34 UTC as originating from an unbound hyperbolic orbit with 99.999\% confidence. In 2022, the U.S. Department of Defense has since verified that "the velocity estimate reported to NASA is sufficiently accurate to indicate an interstellar trajectory," making the object the first detected interstellar object and the first detected interstellar meteor. Here, we discuss the dynamical and compositional properties of CNEOS 2014-01-08, and describe our plan for an expedition to retrieve meteoritic fragments from the ocean floor.

The properties of the disk/jet coupling in quiescent black hole low mass X-ray binaries (BH LMXBs) are still largely unknown. In this paper we present the first quasi-simultaneous radio and X-ray detection in quiescence of the BH LMXB MAXI J1348$-$630, which is known to display a hybrid disk/jet connection that depends on the accretion rate. We performed deep X-ray and radio observations using the Chandra X-ray Observatory and the Australia Telescope Compact Array. MAXI J1348$-$630 is detected for the first time in quiescence at an X-ray luminosity $L_{\rm X} = (7.5 \pm 2.9) \times 10^{30} (D/2.2 \ {\rm kpc})^2$ erg s$^{-1}$: one of the lowest X-ray luminosities observed for a quiescent BH LMXB, possibly implying a short orbital period for the system. MAXI J1348$-$630 is also detected in radio at $L_{\rm R} = (4.3 \pm 0.9) \times 10^{26} (D/2.2 \ {\rm kpc})^2$ erg s$^{-1}$. These detections allow us to constrain the location of MAXI J1348$-$630 on the radio/X-ray diagram in quiescence, finding that the source belongs to the standard (radio-loud) track in this phase. This provides a strong confirmation that hybrid-correlation sources follow the standard track at low luminosities and down to quiescence, thus improving our knowledge of the disk/jet connection in BH LMXBs.

In addition to the Gnevyshev-Ohl rule (GOR), the relation of the odd cycle with the subsequent even one in the 22-year Hale solar cycle was found. It is shown that 3 years before the 11-year minimum $m$, the value of the relative sunspot number SN in an odd cycle is closely related to the value of the maximum in the next even cycle (correlation coefficient $\rho=0.94$), and the same relation of an odd cycle with the previous even one is weaker. Like GOR, cycles are linked in pairs, but opposite to the Rule. Based on this result, we propose to use SN$_{m-3}$ on the descending phase of the previous odd cycle as a precursor of the subsequent EVEN cycle (Figure 3a) -- a precursor called MI3E. For the prediction of an odd cycle or a prediction without consideration of parity (as in the article by Braj\v{s}a et al., 2022), this method gives less reliable results. To predict the amplitude of an ODD cycle, we propose to use the precursor of the seventh year to its maximum $M$ MA7O -- SN$_{M-7}$ on the descending phase of the previous even cycle (Figure 3b). It turned out that in this case, we can also predict the years near the maximum with a high correlation coefficient ($\rho=0.90{-}0.94$). Thus, the proposed approaches allow us to predict cycles of different parity. According to our prediction, the current solar Cycle 25 in 2023 will reach a maximum of 154 units with a prediction error of $\pm25$ (68% confidence) and $\pm53$ (95% confidence). In 2024, SN will be almost as high as in 2023 -- 147 units, so with smaller time averaging scales, the maximum will fall at the end of 2023.

The landscape of high- and ultra-high-energy astrophysics has changed in the last decade, largely due to the inflow of data collected by large-scale cosmic-ray, gamma-ray, and neutrino observatories. At the dawn of the multimessenger era, the interpretation of these observations within a consistent framework is important to elucidate the open questions in this field. CRPropa 3.2 is a Monte Carlo code for simulating the propagation of high-energy particles in the Universe. This version represents a major leap forward, significantly expanding the simulation framework and opening up the possibility for many more astrophysical applications. This includes, among others: efficient simulation of high-energy particles in diffusion-dominated domains, self-consistent and fast modelling of electromagnetic cascades with an extended set of channels for photon production, and studies of cosmic-ray diffusion tensors based on updated coherent and turbulent magnetic-field models. Furthermore, several technical updates and improvements are introduced with the new version, such as: enhanced interpolation, targeted emission of sources, and a new propagation algorithm (Boris push). The detailed description of all novel features is accompanied by a discussion and a selected number of example applications.

Jets are one of the most common eruptive events in the solar atmosphere, and they are believed to be important in the context of coronal heating and solar wind acceleration. We present an observational study on a sequence of jets with the data acquired with the Solar Dynamics Observatory (SDO) and the Interface Region Imaging Spectrograph (IRIS). This sequence is peculiar in that an EUV jet, $\sim29\arcsec$ long and with a dome-like base, appears to be a consequence of a series of transition region (TR) micro-jets that are a few arcsecs in length.We find that the occurrence of any TR micro-jets is always associated with the change of geometry of micro-loops at the footpoints of the microjets. A bundle of TR flux ropes is seen to link a TR micro-jet to the dome-like structure at the base of the EUV jet. This bundle rises as a response to the TR micro-jets, with the rising motion eventually triggering the EUV jet. We propose a scenario involving a set of magnetic reconnections, in which the series of TR micro-jets are associated with the processes to remove the constraints to the TR flux ropes and thus allow them to rise and trigger the EUV jet. Our study demonstrates that small-scale dynamics in the lower solar atmosphere are crucial in understanding the energy and mass connection between the corona and the solar lower atmosphere, even though many of them might not pump mass and energy to the corona directly.

Gamma-ray bursts (GRBs) are produced during the propagation of ultra-relativistic jets. While our understanding of these jets have improved notably during the last decades, it is currently impossible to study directly the jet close to the central source, due to the high opacity of the medium. In this paper, we present numerical simulations of relativistic jets propagating through a massive, stripped envelope star associated to long GRBs, breaking out of the star and accelerating into the circumstellar medium. We compute the resulting gravitational wave (GW) signal, showing that several key parameters of the jet propagation can be directly determined by the associated GW signal. The signal presents two peaks, the first one corresponding to the jet duration, while the second one corresponding to the end of the acceleration phase. Depending on the observer location (with respect to the jet axis) this peak corresponds to the break-out time for observer located close to the jet axis (which in turn depends on the stellar size), or to much larger times (corresponding to the end of the acceleration phase) for off-axis observers. We also show that the slope of the GW signal before and around the first peak tracks the jet luminosity history and the structure of the progenitor star. The amplitude of the GW signal is $h_+D \sim$ hundreds to several thousands. Although this signal, for extragalactic sources, is outside the range of detectability of current GW detectors, it can be detected by future detectors as BBO, DECIGO and ALIA. Our results illustrate that future detections of GW associated to GRB jets will represent a revolution in our understanding of this phenomenon.

In the first image of the James Webb Space Telescope (JWST) of SMACS J0723.3-7327, one of the most outstanding features is the emergence of a large number of red spiral galaxies, because such red spiral galaxies are only a few percent in the number fraction among nearby spiral galaxies. While these apparently red galaxies were already detected with the Spitzer Space Telescope at $\sim3-4{\rm \mu m}$, the revolutionized view from JWST's unprecedented spatial resolution has unveiled their hidden spiral morphology for the first time. Within the red spiral galaxies, we focus on the three most highly red galaxies that are very faint in the $<0.9\,{\rm \mu m}$ bands and show red colors in the $2-4\,{\rm \mu m}$ bands. Our study finds that the three extremely red spiral galaxies are likely to be in the Cosmic Noon (i.e., $1 < z < 3$) and could be consistent with passive (i.e., $\sim$ zero star-formation rates) galaxies having moderate dust reddening (i.e., $A_{\rm V}\sim1\,{\rm mag}$). These "red spiral" galaxies would be interesting, potentially new population of galaxies, as we start to see their detailed morphology using JWST, for the first time. Finally, we note that the spectral energy distribution of these red $z\sim2.5$ galaxies could mimic $z>10$ Lyman break galaxies and contaminate to $z>10$ galaxy samples, especially when they were faint and small.

Inferring the values and uncertainties of cosmological parameters in a cosmology model is of paramount importance for modern cosmic observations. In this paper, we use the simulation-based inference (SBI) approach to estimate cosmological constraints from a simplified galaxy cluster observation analysis. Using data generated from the Quijote simulation suite and analytical models, we train a machine learning algorithm to learn the probability function between cosmological parameters and the possible galaxy cluster observables. The posterior distribution of the cosmological parameters at a given observation is then obtained by sampling the predictions from the trained algorithm. Our results show that the SBI method can successfully recover the truth values of the cosmological parameters within the 2{\sigma} limit for this simplified galaxy cluster analysis, and acquires similar posterior constraints obtained with a likelihood-based Markov Chain Monte Carlo method, the current state-of the-art method used in similar cosmological studies.

High-energetic gamma rays from astrophysical targets constitute a unique probe for annihilation or decay of heavy particle Dark Matter (DM). After several decades, diverse null detections have resulted in strong constraints for DM particle masses up to the TeV scale. While the gamma-ray signature is expected to be universal from various targets, astrophysical uncertainties strongly vary and weaken the limits. At the same time, spurious signals may originate from non-DM related processes. The many gamma-ray targets in the extragalactic sky being searched for DM play a crucial role to control these uncertainties and to achieve an unambiguous DM detection. Lately, a large progress has been made in combined analyses of TeV DM candidates towards different targets by using data from various instruments and over a wide range of gamma-ray energies. These approaches not only resulted in an optimal exploitation of existing data and a superior sensitivity, but also helped to level out target- and instrument-related uncertainties. This review gathers all searches in the extragalactic sky performed so far with the space-borne Fermi-LAT, and the ground-based Imaging Atmospheric Cherenkov telescopes and the HAWC Water Cherenkov observatory. We discuss the different target classes and provide a complete list of all analyses so far.

We have obtained optical spectroscopy and photometry data during four years after the event. The long-term photometric light-curve and the equivalent widths of the Halpha and He I 6678 lines were used to monitor the state of the Be star disk. The Halpha line profiles show evidence for V/R variability that was accounted for by fitting the Halpha spectral line profile with two Gaussian functions. We divided our data into three phases according to the intensity of the X-ray, optical, and infrared emission. Phase I covers the rise and decay of the giant X-ray outburst that took place in October to November 2017. We interpret phase II as the dissipation of the Be star equatorial disk and phase III as its recovery. The timescale of a complete formation and dissipation process is about 1250 days. The epoch when the dissipation process stopped and the reformation period began is estimated to be around MJD 58530. We find a delay of about 100 to 200 days between the minimum of the optical or infrared intensity and the strength of the Halpha line after the X-ray outburst, which may indicate that the dissipation of the disk begins from the inner parts. The motion of the density perturbation inside the disk is prograde, with a V/R quasi-period of about four years. The source shows a positive correlation in the (B-V) color index versus V-band magnitude diagram, which implies that the system is seen at a small or moderate inclination angle.

The Ultraviolet Transient Astronomy Satellite (ULTRASAT) is a space-borne near UV telescope with an unprecedented large field of view (200 sq. deg.). The mission, led by the Weizmann Institute of Science and the Israel Space Agency in collaboration with DESY (Helmholtz association, Germany) and NASA (USA), is fully funded and expected to be launched to a geostationary transfer orbit in Q2/3 of 2025. With a grasp 300 times larger than GALEX, the most sensitive UV satellite to date, ULTRASAT will revolutionize our understanding of the hot transient universe, as well as of flaring galactic sources. We describe the mission payload, the optical design and the choice of materials allowing us to achieve a point spread function of ~10arcsec across the FoV, and the detector assembly. We detail the mitigation techniques implemented to suppress out-of-band flux and reduce stray light, detector properties including measured quantum efficiency of scout (prototype) detectors, and expected performance (limiting magnitude) for various objects.

We update the constraints on the tensor-to-scalar ratio $r$ and the tensor spectral index $n_t$ using 10 dataset combinations from BICEP/Keck Array 2015/2018, Planck PR3/PR4 and LIGO-Virgo-KAGRA collaboration. This is done by fitting the complete $\Lambda$CDM$+r+n_t$ model with two different approaches for the tensor sector, whose strengths and weaknesses are also analyzed. We find that our most constraining dataset combination yields $r<0.028$ and $-1.37 < n_t < 0.42$ at 95% CL. On a separate fashion, we also explore the consequences on the tensor sector of the common-signal detected by NANOGrav collaboration, finding $r<0.033$ and $0.47 < n_t < 0.85$ at 95% CL. Consistently to what has been already shown in the literature, the latter result confirms that, if the NANOGrav detection would be due to a primordial inflationary gravitational wave background, then a simple power-law cannot reconcile the constraints from different datasets.

The collisional pumping of H$_2$O and CH$_3$OH masers in magnetohydrodynamic nondissociative C-type shocks is considered. A grid of C-type shock models with speeds in the range $5-70$ km s$^{-1}$ and preshock gas densities $n_{\rm H_2,0} = 10^4-10^7$ cm$^{-3}$ is constructed. The large velocity gradient approximation is used to solve the radiative transfer equation in molecular lines. The para-H$_2$O 183.3 GHz and ortho-H$_2$O 380.1 and 448.0 GHz transitions are shown to be inverted and to have an optical depth along the shock velocity $\vert \tau \vert \sim 1$ at relatively low gas densities in the maser zone, $n_{\rm H_2} \gtrsim 10^5-10^6$ cm$^{-3}$. Higher gas densities, $n_{\rm H_2} \gtrsim 10^7$ cm$^{-3}$, are needed for efficient pumping of the remaining H$_2$O masers. Simultaneous generation of H$_2$O and class I CH$_3$OH maser emission in a shock is possible at preshock gas densities $n_{\rm H_2,0} \approx 10^5$ cm$^{-3}$ and shock speeds in the range $u_{\rm s} \approx 17.5-22.5$ km s$^{-1}$. The possibility of detecting class I CH$_3$OH and para-H$_2$O 183.3 GHz masers in star-forming regions and near supernova remnants is investigated.

We present the third data release of the European Space Agency's Gaia mission, GDR3. The GDR3 catalogue is the outcome of the processing of raw data collected with the Gaia instruments during the first 34 months of the mission by the Gaia Data Processing and Analysis Consortium. The GDR3 catalogue contains the same source list, celestial positions, proper motions, parallaxes, and broad band photometry in the G, G$_{BP}$, and G$_{RP}$ pass-bands already present in the Early Third Data Release. GDR3 introduces an impressive wealth of new data products. More than 33 million objects in the ranges $G_{rvs} < 14$ and $3100 <T_{eff} <14500 $, have new determinations of their mean radial velocities based on data collected by Gaia. We provide G$_{rvs}$ magnitudes for most sources with radial velocities, and a line broadening parameter is listed for a subset of these. Mean Gaia spectra are made available to the community. The GDR3 catalogue includes about 1 million mean spectra from the radial velocity spectrometer, and about 220 million low-resolution blue and red prism photometer BPRP mean spectra. The results of the analysis of epoch photometry are provided for some 10 million sources across 24 variability types. GDR3 includes astrophysical parameters and source class probabilities for about 470 million and 1500 million sources, respectively, including stars, galaxies, and quasars. Orbital elements and trend parameters are provided for some $800\,000$ astrometric, spectroscopic and eclipsing binaries. More than $150\,000$ Solar System objects, including new discoveries, with preliminary orbital solutions and individual epoch observations are part of this release. Reflectance spectra derived from the epoch BPRP spectral data are published for about 60\,000 asteroids. Finally, an additional data set is provided, namely the Gaia Andromeda Photometric Survey (abridged)

To understand the origin of the diversity observed in exoplanetary systems, it is crucial to characterize the early stages of their formation, represented by Solar-type protostars. Likely, the gaseous chemical content of these objects directly depends on the composition of the dust grain mantles formed before the collapse. Directly retrieving the ice mantle composition is challenging, but it can be done indirectly by observing the major components, such as NH3 and CH3OH at cm wavelengths, once they are released into the gas-phase during the warm protostellar stage. We observed several CH3OH and NH3 lines toward three Class 0 protostars in NGC1333 (IRAS 4A1, IRAS 4A2, and IRAS 4B), at high angular resolution (1"; ~300 au) with the VLA interferometer at 24-26 GHz. Using a non-LTE LVG analysis, we derived a similar NH3/CH3OH abundance ratio in the three protostars (<0.5, 0.015-0.5, and 0.003-0.3 for IRAS 4A1, 4A2, and 4B, respectively). Hence, we infer they were born from pre-collapse material with similar physical conditions. Comparing the observed abundance ratios with astrochemical model predictions, we constrained the dust temperature at the time of the mantle formation to be ~17 K, which coincides with the average temperature of the southern NGC 1333 diffuse cloud. We suggest that a brutal event started the collapse that eventually formed IRAS 4A1, 4A2 and 4B, which,therefore, did not experience the usual pre-stellar core phase. This event could be the clash of a bubble with NGC 1333 south, that has previously been evoked in the literature.

We present the direct imaging discovery of a low-mass companion to the nearby accelerating F star, HIP 5319, using SCExAO coupled with the CHARIS, VAMPIRES, and MEC instruments in addition to Keck/NIRC2 imaging. CHARIS $JHK$ (1.1-2.4 $\mu$m) spectroscopic data combined with VAMPIRES 750 nm, MEC $Y$, and NIRC2 $L_{\rm p}$ photometry is best matched by an M3--M7 object with an effective temperature of T=3200 K and surface gravity log($g$)=5.5. Using the relative astrometry for HIP 5319 B from CHARIS and NIRC2 and absolute astrometry for the primary from $Gaia$ and $Hipparcos$ and adopting a log-normal prior assumption for the companion mass, we measure a dynamical mass for HIP 5319 B of $31^{+35}_{-11}M_{\rm J}$, a semimajor axis of $18.6^{+10}_{-4.1}$ au, an inclination of $69.4^{+5.6}_{-15}$ degrees, and an eccentricity of $0.42^{+0.39}_{-0.29}$. However, using an alternate prior for our dynamical model yields a much higher mass of 128$^{+127}_{-88}M_{\rm J}$. Using data taken with the LCOGT NRES instrument we also show that the primary HIP 5319 A is a single star in contrast to previous characterizations of the system as a spectroscopic binary. This work underscores the importance of assumed priors in dynamical models for companions detected with imaging and astrometry and the need to have an updated inventory of system measurements.

We present photometric analysis of three bright red nova progenitor contact binary systems: ASAS J082151-0612.6, TYC 7281-269-1 and TYC 7275-1968-1. The primary components in all three systems are solar-type low mass stars with radii somewhat larger than their Zero Age Main Sequence counterparts. The secondaries, as in most contact binary systems, have radii and luminosities well above their main sequence counterparts. All three have extreme low mass ratios ranging from 0.075 to 0.097 and two have high degrees of contact, in excess of 75%. All three have mass ratios and separations below the theoretical values for orbital stability. Chromospheric activity, a hallmark of magnetic activity and magnetic braking, considered important in mediating angular momentum loss, is also explored. All three systems demonstrate the O'Connell effect, and all systems require the introduction of star spots for a better light curve solution. In addition, we show that ASAS J082151-0612.6 and TYC 7281-269-1 have a UV colour excess in the range indicating high chromospheric activity. Another measure of potential significant magnetic activity is X-Ray luminosity; TYC 7275-1968-1, and probably also TYC 7281-269-1, have X-Ray luminosity well above other contact binary systems. We conclude that it is likely that all three are unstable and hence are potential merger candidates.

In this work, we investigated the effect of the cosmic ionization history from Population~III~(Pop.~III) stars in ultracompact minihalos~(UCMHs) using Planck observation data. Although high-redshift astrophysics is not understood yet, UCMHs could host the Pop.~III stars like the halos formed in the standard structure formation scenario. Such Pop.~III stars would emit ionizing photons during their main sequence and facilitate the cosmic reionization in high redshift. To study the effect of the global ionization, we model the cosmic reionization evolution based on the "tanh"-type reionization model which is expressed by $z_{\mathrm{reio}}$ with additional two parameters characterizing the initial mass of UCMHs and the number density of UCMHs. We have implemented the Monte Carlo Markov Chain analysis with the latest Planck observation data for our reionization model. As the result, we found that if the UCMH initial mass is larger than $10^{8.4}\mathrm{M}_{\odot}$, the number density of UCMHs are strictly limited. Then we obtained the constraint on the amplitude of primordial power spectrum through the constraint on the UCMH number density like $\mathcal{A}_{\zeta}\sim 10^{-9}$ in the scales, $k\lesssim 50\mathrm{Mpc}^{-1}$, when we fix $z_{\mathrm{reio}}=6$.

We present here the latest results obtained with the C-RED One camera developed by First Light Imaging for fast ultra-low noise infrared applications. This camera uses the Leonardo Saphira e-APD 320x256 infrared sensor in an autonomous cryogenic environment with a low vibration pulse tube and with embedded readout electronics system. Some recent improvements were made to the camera. The first important one concerns the total noise of the camera. Limited to 1.75 microns wavelength cut-off with proper cold filters, looking at a blackbody at room temperature and f/4 beam aperture, we now measure total noise down to 0.6 e at gain 50 in CDS mode 1720 FPS, dividing previous noise figure by a factor 2. The total camera background of 30-400 e/s is now achieved with a factor 3 of background reduction, the camera also looking at a room temperature blackbody with an F/4 beam aperture. Image bias oscillations, due to electronics grounding scheme, were carefully analyzed and removed. Focal plane detector vibrations transmitted by the pulse tube cooling machine were also analyzed, damped and measured down to 0.3 microns RMS, reducing focal plane vibrations by a factor 3. In addition, a vacuum getter of higher capacity is now used to offer camera operation without camera pumping during months. The camera main characteristics are detailed: pulse tube cooling at 80K with limited vibrations, permanent vacuum solution, ultra-low latency Cameralink full data interface, safety management of the camera by firmware, online firmware update, ambient liquid cooling and reduced weight of 20 kg.

There have been no significant breakthroughs in infrared imagery since the hybridization of III-V or II-VI narrow-bandgap semiconductors on complementary metal-oxide semiconductor (CMOS) read-out integrated circuits (ROICs). The development of third-generation, linear-mode avalanche photodiode arrays (LmAPDs) using mercury cadmium telluride (MCT) has resulted in a significant sensitivity improvement for short-wave infrared (SWIR) imaging. The first dedicated LmAPD device, called SAPHIRA (320x256/24 microns), was designed by Leonardo UK Ltd specifically for SWIR astronomical applications. In the past decade there has been a significant development effort to make larger LmAPD arrays for low-background astronomy. Larger LmAPD formats for ultra-low noise/flux SWIR imaging, currently under development at Leonardo include a 512 x 512 LmAPD array funded by ESO, MPE and NRC Herzberg, a 1k x 1k array funded by NASA and a 2K x 2K device funded by ESA for general scientific imaging applications. The 2048x2048 pixel ROIC has a pitch of 15 microns, 4/8/16 outputs and a maximum frame rate of 10 Hz. The ROIC characterization is scheduled in the third quarter of 2022, while the first arrays will be fabricated by end-2022. The hybridized arrays will be characterized during end-2022. At this time, First Light Imaging will start the development of an autonomous camera integrating this 2Kx2K LmAPD array, based on the unique experience from the C-RED One camera, the only commercial camera integrating the SAPHIRA SWIR LmAPD array.The detector will be embedded in a compact high vacuum cryostat cooled with low vibration pulse at 50-80K which does not require external pumping. Sub-electron readout noise is expected to be achieved with high multiplication gain. Custom cold filters and beam aperture cold baffling will be integrated in the camera.

Abell 3266 is a massive and complex merging galaxy cluster that exhibits significant substructure. We present new, highly sensitive radio continuum observations of Abell 3266 performed with the Australian Square Kilometre Array Pathfinder (0.8$-$1.1 GHz) and the Australia Telescope Compact Array (1.1$-$3.1 GHz). These deep observations provide new insights into recently-reported diffuse non-thermal phenomena associated with the intracluster medium, including a 'wrong-way' relic, a fossil plasma source, and an as-yet unclassified central diffuse ridge, which we reveal comprises the brightest part of a large-scale radio halo detected here for the first time. The 'wrong-way' relic is highly atypical of its kind: it exhibits many classical signatures of a shock-related radio relic, while at the same time exhibiting strong spectral steepening. While radio relics are generally consistent with a quasi-stationary shock scenario, the 'wrong-way' relic is not. We study the spectral properties of the fossil plasma source; it exhibits an ultra-steep and highly curved radio spectrum, indicating an extremely aged electron population. The larger-scale radio halo fills much of the cluster centre, and presents a strong connection between the thermal and non-thermal components of the intracluster medium, along with evidence of substructure. Whether the central diffuse ridge is simply a brighter component of the halo, or a mini-halo, remains an open question. Finally, we study the morphological and spectral properties of the multiple complex radio galaxies in this cluster in unprecedented detail, tracing their evolutionary history.

The axion-nucleon coupling enables the production of axions through the decay of excited ${}^{57}\textrm{Fe}$ isotopes, and axions produced in the Sun through this process are often a target of helioscope searches. We show for the first time that hot, highly magnetic white dwarfs such as ZTF J1901+1458 are a viable target to search for the X-ray signature of axions that were produced by the ${}^{57}\textrm{Fe}$ transition in the core and then converted to photons in the magnetosphere. We calculate that a 100 ks observation of ZTF J1901+1458 with NuSTAR would constrain the coupling of axions to nucleons and photons at a level below the bounds of both current and future planned helioscopes.

The high cosmic abundance and the intermediate volatility and chemical properties of sulfur allow the use of S-bearing species as a tracer of the chemical processes in the atmospheres of hot Jupiters. Nevertheless, despite its properties and relevance as a tracer of the giant planets' formation history, little attention has been paid to this species in the context of hot Jupiter's atmospheres. In this paper, we provide an overview of the abundances of sulfur-bearing species in hot Jupiter atmospheres under different conditions and explore their observability. We use the photochemical kinetics code VULCAN to model hot Jupiter atmospheric disequilibrium chemistry. Transmission spectra for these atmospheres are created using the modelling framework ARCiS. We vary model parameters such as the diffusion coefficient Kzz, and we study the importance of photochemistry on the resulting mixing ratios. Furthermore, we vary the chemical composition of the atmosphere by increasing the metallicity from solar to ~10 times solar. We also explore different C/O ratios. We find that H2S and SO2 are the best candidates for detection between 1 and 10 micron, using a spectral resolution that is representative of the instruments on board the James Webb Space Telescope (JWST). H2S is easiest to detect at an equilibrium temperature of ~1500 K and C/O ratios between 0.7 and 0.9, with the ideal value increasing slightly for increasing metallicity. SO2 is most likely to be detected at an equilibrium temperature of ~1000 K at low C/O ratios and high metallicities. Nevertheless, among these two molecules, we expect SO2 detection to be more common, as is the most favoured scenario from formation models. We conclude that H2S and SO2 will most likely be detected in the coming years with the JWST and that the detection of these species will provide information on atmospheric processes and planet formation scenarios.

We investigate the blue and optical rest-frame sizes (lambda~2300A-4000A) of three compact star-forming regions in a galaxy at z=4 strongly lensed (x30, x45, x100) by the Hubble Frontier Field galaxy cluster A2744 using GLASS-ERS JWST/NIRISS imaging at 1.15um, 1.50mu and 2.0mu with PSF < 0.1". In particular, the Balmer break is probed in detail for all multiply-imaged sources of the system. With ages of a few tens of Myr, stellar masses in the range (0.7-4.0) x 10^6 Msun and optical/ultraviolet effective radii spanning the interval 3 < R_eff < 20 pc, such objects are currently the highest redshift (spectroscopically-confirmed) gravitationally-bound young massive star clusters (YMCs), with stellar mass surface densities resembling those of local globular clusters. Optical (4000A, JWST-based) and ultraviolet (1600A, HST-based) sizes are fully compatible. The contribution to the ultraviolet underlying continuum emission (1600A) is ~30%, which decreases by a factor of two in the optical for two of the YMCs (~4000A rest-frame), reflecting the young ages (<30 Myr) inferred from the SED fitting and supported by the presence of high-ionization lines secured with VLT/MUSE. Such bursty forming regions enhance the sSFR of the galaxy, which is ~10 Gyr^-1. This galaxy would be among the extreme analogs observed in the local Universe having high star formation rate surface density and high occurrence of massive stellar clusters in formation.

The discoveries of two Interstellar Objects (ISOs) in recent years has generated significant interest in constraining their physical properties and the mechanisms behind their formation. However, their ephemeral passages through our Solar System permitted only incomplete characterization. We investigate avenues for identifying craters that may have been produced by ISOs impacting terrestrial Solar System bodies, with particular attention towards the Moon. A distinctive feature of ISOs is their relatively high encounter velocity compared to asteroids and comets. Local stellar kinematics indicate that terrestrial Solar System bodies should have experienced of order unity ISO impacts exceeding 100 km/s. By running hydrodynamical simulations for projectiles of different masses and impact velocities, up to 100 km/s, we show how late-stage equivalence dictates that transient crater dimensions are alone insufficient for inferring the projectile's velocity. On the other hand, the melt volume within craters of a fixed diameter may be a potential route for identifying ISO craters, as faster impacts produce more melt. This method requires that the melt volume scales with the energy of the projectile, while crater diameter scales with the point-source limit (sub-energy). Given that there are probably only a few ISO craters in the Solar System at best, and that transient crater dimensions are not a distinguishing feature for impact velocities at least up to 100 km/s, identification of an ISO crater proves a challenging task. Melt volume and high-pressure petrology may be diagnostic features once large volumes of material can be analyzed in situ.

We present the data release of the Uchuu-SDSS galaxies: a set of 32 high-fidelity galaxy lightcones constructed from the large Uchuu 2.1 trillion particle $N$-body simulation using Planck cosmology. We adopt subhalo abundance matching to populate the Uchuu-box halo catalogues with SDSS galaxy luminosities. These cubic box galaxy catalogues generated at several redshifts are combined to create the set of lightcones with redshift-evolving galaxy properties. The Uchuu-SDSS galaxy lightcones are built to reproduce the footprint and statistical properties of the SDSS main galaxy survey, along with stellar masses and star formation rates. This facilitates direct comparison of the observed SDSS and simulated Uchuu-SDSS data. Our lightcones reproduce a large number of observational results, such as the distribution of galaxy properties, the galaxy clustering, the stellar mass functions, and the halo occupation distributions. Using the simulated and real data we select samples of bright red galaxies at $z_\mathrm{eff}=0.15$ to explore Redshift Space Distortions and Baryon Acoustic Oscillations (BAO) utilizing a full-shape analytical model of the two-point correlation function. We create a set of 5100 galaxy lightcones using GLAM N-body simulations to compute covariance errors. We report a $\sim 30\%$ precision increase on $f\sigma_8$, due to our better estimate of the covariance matrix. From our BAO-inferred $\alpha_{\parallel}$ and $\alpha_{\perp}$ parameters, we obtain the first SDSS measurements of the Hubble and angular diameter distances $D_\mathrm{H}(z=0.15) / r_d = 27.9^{+3.1}_{-2.7}$, $D_\mathrm{M}(z=0.15) / r_d = 5.1^{+0.4}_{-0.4}$. Overall, we conclude that the Planck LCDM cosmology nicely explains the observed large-scale structure statistics of SDSS. All data sets are made publicly available.

Understanding the faint end of quasar luminosity function at a high redshift is important since the number density of faint quasars is a critical element in constraining ultraviolet (UV) photon budgets for ionizing the intergalactic medium (IGM) in the early universe. Here, we present quasar LF reaching $M_{1450} \sim -22.0$ AB mag at $z\sim5$, about one magnitude deeper than previous UV LFs. We select quasars at $z\sim5$ with a deep learning technique from deep data taken by the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP), covering a 15.5 deg$^2$ area. Beyond the traditional color selection method, we improved the quasar selection by training an artificial neural network for distinguishing $z\sim5$ quasars from non-quasar sources based on their colors and adopting the Bayesian information criterion that can further remove high-redshift galaxies from the quasar sample. When applied to a small sample of spectroscopically identified quasars and galaxies, our method is successful in selecting quasars at $\sim83 \%$ efficiency ($5/6$) while minimizing the contamination rate of high-redshift galaxies ($1/8$) by up to three times compared to the selection using color selection alone ($3/8$). The number of our final quasar candidates with $M_{1450} < -22.0$ mag is 35. Our quasar UV LF down to $M_{1450} = -22$ mag or even fainter ($M_{1450} = -21$ mag) suggests a rather low number density of faint quasars and the faint-end slope of $-1.6^{+0.21}_{-0.19}$, favoring a scenario where quasars play a minor role in ionizing the IGM at high redshift.

Ground-based, high-resolution bolometric (sub)millimeter continuum mapping observations on spatially extended target sources are often subject to significant missing fluxes. This hampers accurate quantitative analyses. Missing flux can be recovered by fusing high-resolution images with observations that preserve extended structures. However, the commonly adopted image fusion approaches do not maintain the simplicity of the beam response function and do not try to elaborate the details of the yielded beam response functions. These make the comparison of the observations at multiple wavelengths not straightforward. We present a new algorithm, J-comb, which combines the high and low-resolution images linearly. By applying a taper function to the low-pass filtered image and combining it with a high-pass filtered image using proper weights, the beam response functions of our combined images are guaranteed to have near-Gaussian shapes. This makes it easy to convolve the observations at multiple wavelengths to share the same beam response functions. Moreover, we introduce a strategy to tackle the specific problem that the imaging at 850 um from the present-date ground-based bolometric instrument and that taken with the Planck satellite do not overlap in the Fourier domain. We benchmarked our method against two other widely-used image combination algorithms, CASA-feather and MIRIAD-immerge, with mock observations of star-forming molecular clouds. We demonstrate that the performance of the J-comb algorithm is superior to those of the other two algorithms. We applied the J-comb algorithm to real observational data of the Orion A star-forming region. We successfully produced dust temperature and column density maps with ~10" angular resolution, unveiling much greater details than the previous results.

Small, medium and large arrays for the study of cosmic rays of ultra-high energies existing are aimed at obtaining information about our Galaxy and metagalactic space. Concretely search and study of astronomical objects, that forms flux of relativistic particles that fill outer space. The drift and interaction of such particles with magnetic field and shock waves taking place in interstellar space causes the same interest. The shape of the energy spectrum of cosmic rays in the energy range $10^{15}-10^{18}$ eV, where "knee" and "second knee" is observed, can be formed as a superposition of the partial spectra of various chemical elements. Verification of galactic models, using recent experimental spectral data, makes it possible to study the nature of the galactic and metagalactic components of cosmic rays. The paper presents the result of the energy spectrum of cosmic rays in the range $10^{16}-10^{18}$ eV measurements obtained at the Small Cherenkov array -- a part of the Yakutsk array.

This paper investigates the impact of dust size distribution on magnetic resistivity. In particular, we focus on its impact when the maximum dust size significantly increases from sub-micron. The first half of the paper describes our calculation method for magnetic resistivity based on the model of \citet{1987ApJ...320..803D} and shows that the method reproduces the results of a more realistic chemical reaction network calculations reasonably well. Then, we describe the results of the resistivity calculations for dust distributions with large maximum dust grains. Our results show that resistivity tends to decrease with dust growth, which is particularly true when the dust size power exponent $q$ is $q=2.5$. On the other hand, the decrease is less pronounced when the dust size power exponent $q$ is $q=3.5$, i.e., when the small dust is also responsible for the dust cross-section. Our results suggest that detailed dust coagulation and fragmentation processes play a vital role in the magnetic resistivities in protostar formation.

The paper presents data on the muon component with a threshold \(\varepsilon_{thr} \geq\) 1 GeV. Air showers were registered at the Yakutsk array during almost 50 years of continuous air shower observations. The characteristics of muons are compared with calculations of QGSjetII-04 and EPOS LHC models for a proton and an iron nucleus. There is a muon deficit in the models, at energies greater than 5 EeV. To make an agreement between experimental data and simulations on muons, further tuning of the models is required.

We report the discovery of four quasars with $M_{1450} \gtrsim -25.0$ mag at $z\sim5$ and supermassive black hole mass measurement for one of the quasars. They were selected as promising high-redshift quasar candidates via deep learning and Bayesian information criterion, which are expected to be effective in discriminating quasars from the late-type stars and high-redshift galaxies. The candidates were observed by the Double Spectrograph on the Palomar 200-inch Hale Telescope. They show clear Ly$\alpha$ breaks at about 7000-8000 \r{A}, indicating they are quasars at $4.7 < z < 5.6$. For HSC J233107-001014, we measure the mass of its supermassive black hole (SMBH) using its C\Romannum{4}$\lambda 1549$ emission line. The SMBH mass and Eddington ratio of the quasar are found to be $\sim 10^8 M_{\odot}$ and $\sim 0.6$, respectively. This suggests that this quasar possibly harbors a fast growing SMBH near the Eddington limit despite its faintness ($L_{\rm Bol} < 10^{46}$ erg s$^{-1}$). Our 100 $\%$ quasar identification rate supports high efficiency of our deep learning and Bayesian information criterion selection method, which can be applied to future surveys to increase high-redshift quasar sample.

The study presents a simplified method to identify potential bright red nova progenitors based on the amplitude of the light curve and infrared (J-H) colour of a contact binary system. We employ published criteria for contact binary orbital instability to show that the amplitude of the light curve for a given contact system with a low mass (< 1.4Msun) primary must be less than a specified value for it to be potentially unstable. Using this we search the photometric data of a large survey to identify about 50 potential bright red nova progenitors. We analyse the survey photometry of each to determine the mass ratio and from the estimated mass of the primary other physical parameters of the systems. We show that each system has physical characteristics indicating potential orbital instability. Using the absolute parameters from our sample we model the expected instability separation and period for low mass contact binary systems

Can atmospheric waves in planet-hosting solar-like stars substantially resonate to tidal forcing? Substantially at a level of impacting the space weather or even of being dynamo-relevant? In particular, low-frequency Rossby waves, which have been detected in the solar near-surface layers, are predestined at responding to sunspot cycle-scale perturbations. In this paper, we seek to address these questions as we formulate a forced wave model for the tachocline layer, which is widely considered as the birthplace of several magnetohydrodynamic planetary waves, i.e., Rossby, inertia-gravity (Poincar\'{e}), Kelvin, Alfv\'{e}n and gravity waves. The tachocline is modeled as a shallow plasma atmosphere with an effective free surface on top that we describe within the Cartesian $\beta$-plane approximation. As a novelty to former studies, we equip the governing equations with a conservative tidal potential and a linear friction law to account for dissipation. We combine the linearized governing equations to one decoupled wave equation, which facilitates an easily approachable analysis. Analytical results are presented and discussed within several interesting free, damped and forced wave limits for both mid-latitude and equatorially trapped waves. For the idealized case of a single tide generating body following a circular orbit, we derive an explicit analytic solution that we apply to our Sun for estimating leading-order responses to Jupiter. Our analysis reveals that Rossby waves resonating to low-frequency perturbations can potentially reach considerable velocity amplitudes in the order of $10^1 - 10^2\, {\rm cm}\, {\rm s}^{-1}$, which, however, strongly rely on the yet unknown total dissipation.

The logarithmic extinction coefficient, c(H$\beta$), is usually derived using the H$\alpha$/H$\beta$ ratio for case B recombination and assuming standard values of electron density and temperature. However, the use of strong Balmer lines can lead to selection biases when studying regions with different surface brightness, such as extended nebulae, with the use of single integral field spectroscopy observations, since, in some cases, the H$\alpha$ line can be saturated in moderate to long exposures. In this work, we present a method to derive extinction corrections based only on the weaker lines of HeI, taking into account the presence of triplet states in these atoms and its influence on recombination lines. We have applied this procedure to calculate the extinction of different regions of the 30 Doradus nebula from MUSE integral-field spectroscopy data. The comparison between helium and hydrogen c(H$\beta$) determinations has been found to yield results fully compatible within the errors and the use of both sets of lines simultaneously reduces considerably the error in the derivation.

We consider the largest sample of 561 edge-on galaxies observed with integral field units by the MaNGA survey and find 300 galaxies where the ionised gas shows a negative vertical gradient (lag) in its rotational speed. We introduce the stop altitude as the distance to the galactic midplane at which the gas rotation should stop in the linear approximation. We find correlations between the lags, stop altitude and galactic mass, stellar velocity dispersion and overall Sersic index. We do not find any correlation of the lags or stop altitude with the star formation activity in the galaxies. We conclude that low mass galaxies (log(M*/Mo) < 10) with low Sersic index and with low stellar velocity dispersion posses a wider "zone of influence" in the extragalactic gas surrounding them with respect to higher mass galaxies that have a significant spherical component. We estimated the trend of the vertical rotational gradient with radius and find it flat for most of the galaxies in our sample. A small subsample of galaxies with negative radial gradients of lag has an enhanced fraction of objects with aged low surface brightness structures around them (e.g. faint shells), which indicates that noticeable accretion events in the past affected the extraplanar gas kinematics and might have contributed to negative radial lag gradients. We conclude that an isotropic accretion of gas from the circumgalactic medium plays a significant role in the formation of rotation velocity lags.

The images of the supermassive black holes Sgr A* and M87* by the Event Horizon Telescope (EHT) collaboration mark a special milestone in the history of the subject. For the first time we are able to see the shadow of black holes, testing basic predictions of the theory of general relativity. We are also now learning more about the fundamental astrophysical processes close to the event horizon that help to shape entire galaxies and even parts of our cosmos. The ultimate result was only possible due to a large collaborative effort of scientists and institutions around the world. The road towards these images was the result of a long sociological and scientific process. It started with early pathfinder experiments and a few simple ideas that were remarkably successful in predicting the basic observational signatures to look for. This was based on the premise that black holes are inherently simple objects. Here I describe this journey and some lessons learned from a personal perspective.

We put constraints on the normalized energy density in gravitaional waves from cosmic domain walls (DWs) by searching for the stochastic gravitational-wave background (SGWB) in the data of Advanced LIGO and Virgo's first three observing runs. By adopting a phenomenological broken power-law model, we obtain the upper limit of normalized energy density of SGWB generated by DWs in the peak frequency band $10\sim200$ Hz, and get the most stringent limitation at the peak frequency $f_*=35$ Hz, namely $\Omega_\text{DW}(f_*=35\,\text{Hz})<1.4\times10^{-8}$ at $95\%$ confidence level (CL). Subsequently, we work out the constraints on the parameter space in the appealing realization of DW structure -- the heavy axion model which can avoid the so-called quality problem.

Lately a specific kind of blazars drew the attention of the gamma-ray astronomy community: the extreme TeV BL Lacs, blazars that present an extremely energetic and hard emission at very high-energy. Explaining their features is still an open challenge, in fact the most used phenomenological models have difficulties to satisfactorily reproduce their SED. Based on a scenario we have recently proposed, we suppose that the non-thermal particles are firstly accelerated by a jet recollimation shock, which induces turbulence in the rest of the jet. Non-thermal particles are further accelerated by the turbulence, which hardens the particle spectra and accordingly the radiative emission. Given the physical properties of the plasma, as inferred by emission models, we expect a strong impact of the accelerating particles on the turbulence. Assuming isotropy and homogeneity, the interaction between non-thermal particles and turbulence and their spectra is modeled solving a system of two non-linear, coupled Fokker-Planck equations, while the radiative emission is calculated through the Synchrotron Self Compton model. The emission predicted by our model is then compared with the prototype extreme TeV BL Lac object 1ES 0229+200 and the parameters obtained to reproduce its SED are in line with the expectations.

The earliest JWST observations have revealed an unexpected abundance of super-early ($z>10$), massive ($M_*\approx 10^9\, M_\odot$) galaxies at the bright-end ($M_{\rm UV}\approx -21$) of the ultraviolet luminosity function (UV LF). We present a minimal physical model that explains the observed galaxy abundance at $z=10-14$. The model primarily combines (a) the halo mass function, with (b) an obscured star formation fraction prescription that is consistent with findings of the ALMA REBELS dusty galaxy survey. It has been successfully tested on well-known UV LFs up to $z=7$. The weak evolution from $z=7$ to $z\approx 14$ of the LF bright-end arises from a conspiracy between a decreasing dust attenuation, making galaxies brighter, that almost exactly compensates for the increasing shortage of their host halos. The model also predicts that galaxies at $z > 11$ should contain negligible amounts of dust. We speculate that dust could have been efficiently ejected during the very first phases of galaxy build-up.

Galaxy quenching is a critical step in galaxy evolution. In this work, we present a statistical study of galaxy quenching in 17 cluster candidates at 0.5<z<1.0 in the COSMOS field. We selected cluster members with a wide range of stellar mass and environment to study their mass and environment dependence. Member galaxies are classified into star-forming, quiescent and recently-quenched galaxies (RQG) using the rest-frame UVJ diagram. We further separated fast and slow quenching RQGs by model evolutionary tracks on the UVJ diagram. We defined the quenching efficiency as the ratio of RQGs over star-forming galaxies and the quenching stage as the ratio of RQGs over quiescent galaxies to quantify the quenching processes. We found quenching efficiency is enhanced by both higher stellar mass and denser environment. Massive or dense environment galaxies quench earlier. Slow quenching is more dominant for massive galaxies and at lower redshifts, but no clear dependence on the environment is found. Our results suggest that low-mass galaxies in dense environments are likely quenched through a short-timescale process such as ram pressure stripping, while massive galaxies in a sparse environment are mostly quenched by a longer-timescale process. Using the line strength of H$\delta$ and [OII], we confirmed that our UVJ method to select RQGs agrees with high S/N DEIMOS spectra. However, we caution that the visibility time (duration of a galaxy's stay in the RQG region on the UVJ diagram) may also depend on mass or environment. The method introduced in this work can be applied to RQG candidates for future statistical RQG spectroscopic surveys. The systematic spectroscopic RQG study will disentangle the degeneracy between visibility time and quenching properties.

We present a star formation rate function (SFRF) at $z\sim5.8$ based on star formation rate (SFR) derived by spectral energy distribution (SED) fitting on data from rest-frame UV to optical wavelength of galaxies in the CANDELS GOODS-South and North fields. The resulting SFRF shows an excess to the previous estimations by using rest-frame UV luminosity functions (LFs) corrected for the dust extinction with the IRX-$\beta$ relation, and is comparable to that estimated from a far-infrared LF. This suggests that the contribution of dust-obscured intensively star-forming galaxies to the total star formation activities at $z\sim6$ is underestimated with the previous approach based only on rest-frame UV observations. We parameterize the SFRF with using the Schechter function and obtain the best-fit parameter of the characteristic SFR (${\rm SFR}^*$) when the faint-end slope and characteristic number density are fixed. The best-fit ${\rm SFR}^*$ is comparable to that at $z\sim2$, when the cosmic star formation activity reaches its peak. Together with SFRF estimations with similar approach using rest-frame UV to optical data, the ${\rm SFR}^*$ is roughly constant from $z\sim2$ to $z\sim6$ and may decrease above $z\sim6$. Since the ${\rm SFR}^*$ is sensitive to the high-SFR end of the SFRF, this evolution of ${\rm SFR}^*$ suggests that the high-SFR end of the SFRF grows in the epoch of reionization to a similar level as at $z\sim2$.

As part of the Advanced LIGO+ (A+) project we have developed, produced, and characterised sensors and electronics to interrogate new optical suspensions. The central element is a displacement sensor with an integrated electromagnetic actuator known as a BOSEM and its readout and drive electronics required to integrate them into LIGO's control and data system. In this paper we report on improvements to the sensors and testing procedures undertaken to meet enhanced performance requirements set out by the A+ upgrade to the detectors. The best devices reach a noise level of $4.5\times 10^{-11}{\rm m}/\sqrt{\rm Hz}$ at a measurement frequency of 1 Hz.

Fast Radio Bursts (FRBs) are a class of highly energetic, mostly extragalactic radio transients lasting for a few milliseconds. While over 600 FRBs have been published so far, their origins are presently unclear, with some theories for extragalactic FRBs predicting accompanying high-energy emission. In this work, we use the Canadian Hydrogen Intensity Mapping Experiment (CHIME) Fast Radio Burst (CHIME/FRB) Project to explore whether any FRB-like radio emission coincides in space and time with 81 gamma-ray bursts (GRBs) detected between 2018 July 17 and 2019 July 8 by Swift/BAT and Fermi/GBM. We do not find any statistically significant, coincident pairs within 3sigma of each other's spatial localization regions and within a time difference of up to one week. In addition to searching for spatial matches between known FRBs and known GRBs, we use CHIME/FRB to constrain FRB-like radio emission before, at the time of, or after the reported high-energy emission at the position of 39 GRBs. Our most constraining radio flux limits in the 400- to 800-MHz band for short gamma-ray bursts (SGRBs) are <50 Jy at 18.6 ks pre-high-energy emission, and <5 Jy at 28.4 ks post-high-energy emission, assuming a 10-ms radio burst width with each limit valid for 60 seconds. We use these limits to constrain models that predict FRB-like prompt radio emission before and after SGRBs. We also place limits as low as 2 Jy for long gamma-ray bursts (LGRBs), but there are no strong theoretical predictions for coincident FRB-like radio emission for LGRBs.

We constrain the value of the Hubble constant $H_0$ by using a detailed model of Fast Radio Burst (FRB) observations from the Australian Square Kilometre Array Pathfinder (ASKAP) and Murriyang (Parkes) radio telescopes. We make use of the redshift-dispersion measure ('Macquart') relationship of FRB populations after accounting for the intrinsic luminosity function, cosmological gas distribution, population evolution, host galaxy contributions to the dispersion measure (DM$_{host}$), and observational biases due to burst duration, dispersion measure and telescope beamshape. Using an updated sample of 16 ASKAP FRBs detected by the Commensal Real-time ASKAP Fast Transients (CRAFT) Survey and localised to their host galaxies, and 60 unlocalised FRBs from Parkes and ASKAP, our best-fitting value of $H_0$ is calculated to be $73_{-8}^{+12}$ km s$^{-1}$ Mpc$^{-1}$. The larger uncertainty than previous FRB works stems from our comprehensive treatment of FRB energetics and DM$_{\rm host}$. Using a prior on $H_0$ covering the 67--74 km s$^{-1}$ Mpc$^{-1}$ range, we estimate a median DM$_{\rm host}$ = $186_{-48}^{+59}$ pc cm$^{-3}$ which exceeds previous estimates from galaxy formation simulations and suggests a greater contribution from the progenitor environment of FRBs. Our large sample of FRB measurements confirms that the FRB population evolves similarly to the star-formation rate. We perform a forecast using a sample of 100 mock FRBs, demonstrating the accuracy of our model, and the potential for high-precision measurements ($\pm 2.5$ km s$^{-1}$ Mpc$^{-1}$) with the coherent FRB search upgrade to ASKAP. This may clarify the current Hubble tension in the near future as the sample of FRBs grows exponentially. Last, we explore a range of sample and selection biases that affect FRB analyses, in both current and future surveys.

We measure the optical variability in $\sim$ 16500 low-redshift (z $\sim$ 0.1) galaxies to map the relations between AGN activity and galaxy stellar mass, specific star-formation rate, half-light radius and bulge-to-total ratio. To do this, we use a reduced $\chi ^2$ variability measure on > 10 epoch lightcurves from the Zwicky Transient Facility and combine with spectroscopic data and derived galaxy parameters from the Sloan Digital Sky Survey. We find that below stellar mass of $10^{11} M_\odot$, galaxies classed as star-forming via the BPT diagram have higher mean variabilities than AGN or composite galaxies. Revealingly, the highest mean variabilities occur in star-forming galaxies in a narrow range of specific star-formation, $-11<\log($sSFR/yr$^{-1})<-10$. In very actively star-forming galaxies $(\log($sSFR/yr$^{-1})>-10)$, the reduced variability implies a lack of instantaneous correlation with star-formation rate. Our results may indicate that a high level of variability, and thus black hole growth, acts as a precursor for reduced star-formation, bulge growth, and revealed AGN-like emission lines. These results add to the mounting evidence that optical variability can act as a viable tracer for low-mass AGNs and that such AGNs can strongly affect their host galaxy.

The MAGIC collaboration has recently analyzed data from a long-term multiwavelength campaign of the $\gamma$-ray blazar TXS 0506+056. In December 2018, it was flaring in the very-high-energy (VHE; $E>100$ GeV) $\gamma$-ray band, but no simultaneous neutrino event was detected. We explore prospects for detecting $\gamma$-rays and neutrinos of hadronic origin, produced both inside and outside the jet of TXS 0506+056, while coherently modeling the observed spectral energy distribution (SED) and neutrino flux upper limits. We constrain the neutrino flux through the restriction from observed X-ray flux on the secondary radiation due to hadronic cascade. We propagate the escaping ultra-high-energy cosmic rays (UHECRs; $E\gtrsim0.1$ EeV) in a random, turbulent extragalactic magnetic field (EGMF). The leptonic emission from the jet dominates the GeV range, whereas the cascade emission from CR interactions in the jet contributes substantially to the X-ray and VHE range. The line-of-sight cosmogenic $\gamma$ rays from UHECRs produce a hardening in the VHE range of the SED. Neutrino signal from the jet showed little or no variability during the MAGIC campaign. Therefore, we infer that the correlation between VHE $\gamma$-rays and neutrino flare is minimal. The luminosity in CRs limits the cosmogenic $\gamma$-ray flux, which, in turn, bounds the RMS value of the EGMF to $\gtrsim 10^{-5}$ nG. The cosmogenic neutrino flux is lower than the IceCube-Gen2 detection potential for 10 yrs of observation. VHE $\gamma$-ray variability should arise from an increased activity inside the jet. Upcoming $\gamma$-ray imaging telescopes, such as the CTA, will be able to constrain the cosmogenic $\gamma$-ray component in the SED of TXS 0506+056. Detecting a steady flux at multi-TeV energies will validate blazars as unambiguous sources of UHECRs.

We demonstrate that VPMS J170850.95+433223.7 is a weak line quasar (WLQ) which is remarkable in several respects. It was already classified as a probable quasar two decades ago, but with considerable uncertainty. The non-significant proper motion and parallax from the Gaia early data release 3 have solidified this assumption. Based on previously unpublished spectra, we show that VPMS J170850.95+433223.7 is a WLQ at z = 2.345 with immeasurably faint broad emission lines in the rest-frame ultraviolet. A preliminary estimate suggests that it hosts a supermassive black hole of ~10^9 M_sun accreting close to the Eddington limit, perhaps at the super-Eddington level. We identify two absorber systems with blueward velocity offsets of 0.05c and 0.1c, which could represent high-velocity outflows, which are perhaps related to the high accretion state of the quasar.

Two fibre integral field units (IFU) are being built in the SAAO fibre-lab for the Robert Stobie Spectrograph's visible arm and the future red arm. Each IFU sits in its own slit-mask cassette and is referred to as a slit-mask IFU (SMI). They will be available some time in 2022. The smaller, 200 micron fibre IFU has 309 X 0.9 arcsec diameter spatial elements covering an elongated hexagonal footprint of 18 X 23 arcsec. The larger, 400 micron fibre IFU has 178 X 1.8 arcsec diameter spatial elements covering an on-sky area of 21 X 44 arcsec. In both cases there are two groups of 13 fibres offset by roughly 50 arcsec on either side of the primary array to sample sky. The 1.8 and 0.9 arcsec spatial resolution SMIs provide median spectral resolution of 1200 and 2400 respectively at H alpha wavelengths in the low resolution mode covering 320 to 740 nm bandpass. At a higher grating angle the SMI will deliver spectral resolution up to 5000 and 10000 with 400 and 200 micron core fibre respectively. A future red-arm will extend the simultaneous wavelength coverage up to 900 nm at a median resolution of 3000/6000 for the same flavors of IFUs. SMIs are inserted in the same fashion as the existing longslit cassettes at the SALT focal plane. Prismatic fold mirrors direct the focal plane into the fibre IFU and then back into the RSS collimator after the fibres are routed 180 deg within the cassette and formatted into a pseudo-slit. Fold prisms ensure that the spectrograph collimator continues to see the same focal plane. In this paper we describe the design, fabrication, assembly and characterization of Slit Mask IFUs.

The nature of galaxy spin is still not fully known. Iye et al (2021) applied a 3D analysis to a dataset of bright SDSS galaxies that was used in the past for photometric analysis. They showed that the distribution of spin directions of spiral galaxies is random, providing a dipole axis with low statistical significance of 0.29$\sigma$. However, to show random distribution, two decisions were made, each can lead to random distribution regardless of the real distribution of the spin direction of galaxies. The first decision was to limit the dataset arbitrarily to z$<$0.1, which is a redshift range in which previous literature already showed that random distribution is expected. More importantly, while the 3D analysis requires the redshift of each galaxy, the analysis was done with the photometric redshift. If the asymmetry existed, its signal is expected to be an order of magnitude weaker than the error of the photometric redshift, and therefore the low statistical signal under these conditions is expected. When using the exact same data without limiting to $z_{phot}<0.1$ and without using the photometric redshift, the distribution of the spin directions in that dataset shows a statistical signal of $>2\sigma$. Code and data for reproducing the analysis are publicly available. These results are in agreement with other experiments with SDSS, Pan-STARRS, HST, and the DESI Legacy Survey. The paper also examines other previous studies that showed random distribution in galaxy spin directions. While further research will be required, the current evidence suggest that large-scale asymmetry between the number of clockwise and counterclockwise galaxies cannot be ruled out.

In this work, we use angular diameter distances of 38 galaxy clusters with joint X-ray/SZE observation to circumvent the circularity problem in the Amati relation for Gamma-ray Bursts (GRBs). Assuming the validity of cosmic-distance duality relation, we obtain the luminosity distance from the cluster angular diameter distance and use that to calculate the isotropic equivalent energy of two different GRB datasets, after restricting the GRB redshift range to $z<0.9$. We then check the validity of the Amati relation for both these datasets. The best-fit Amati relation parameters using galaxy cluster distances as low-redshift anchors are consistent with a previous estimate for the same dataset. The intrinsic scatter which we obtain for the two datasets is about 45% and 15%, and is comparable with that found by other distance anchors used to vet the Amati relation.

The paper shows an analysis of the large-scale distribution of galaxy spin directions of 739,286 galaxies imaged by DES. The distribution of the spin directions of the galaxies exhibits a large-scale dipole axis. Comparison of the location of the dipole axis to a similar analysis with data from SDSS, Pan-STARRS, and DESI Legacy Survey shows that all sky surveys exhibit dipole axes within 52$^o$ or less from each other, well within 1$\sigma$ error. While non-random distribution is unexpected, the findings are consistent across all sky surveys, regardless of the telescope or whether the data were annotated manually or automatically. Possible errors that can lead to the observation are discussed. The paper also discusses previous studies showing opposite conclusions,and analyzes the decisions that led to these results. Although the observation is provocative, and further research will be required, the existing evidence justifies to consider the contention that galaxy spin directions as observed from Earth are not necessarily randomly distributed. Possible explanations can be related to mature cosmological theories, but also to the internal structure of galaxies.

We perform 3D hydrodynamical simulations of new-born neutron stars (NSs) colliding with main-sequence binary companions after supernova explosions. Based on those hydrodynamical models, we construct a semi-analytical formula that describes the drag force inside stars with steep density gradients. We then compute the outcome of NS--companion collisions over a wide range of parameters using the semi-analytical formula. Depending on the direction and magnitude of the natal kick, we find that the collision may lead to various outcomes. For relatively fast kicks and high impact parameters, the NS may penetrate the companion star envelope without merging. By allowing the NS to plunge through companions, the companion can be accelerated to have runaway velocities up to $\sim10$ per cent above the theoretical upper limit considered in classical binary disruption scenarios. The NS can capture and carry away up to a few per cent of the companion envelope as it escapes, which may form pulsar planets or cause outflows through accretion to heat the ejecta from inside and power the supernova light curve. For lower impact parameters, the neutron star will directly merge with the companion and form a Thorne-\.Zytkow object. In intermediate cases, the NS penetrates the companion envelope several times before merging, possibly causing multiple bumps in the supernova light curve like in SN2015bn and SN2019stc.

Several new radio facilities have a field of view and sensitivity well suited for transient searches. This makes it more important than ever to accurately determine transient rates in radio surveys. The work presented here seeks to do this task by using Monte-Carlo simulations. In particular, the user inputs either a real or simulated observational setup, and the simulations code calculates transient rate as a function of transient duration and peak flux. These simulations allow for simulating a wide variety of scenarios including observations with varying sensitivities and durations, multiple overlapping telescope pointings, and a wide variety of light curve shapes with the user having the ability to easily add more. While the current scientific focus is on the radio regime, with examples given here from the MeerKAT telescope in South Africa, the simulations code can be easily adapted to other wavelength regimes.

Measurement of chemical and kinematic structures in prestellar cores is essential for better understanding the star formation process. Here, we study the three prestellar cores (TMC-1C, L1544, and TMC-1) of the Taurus molecular cloud by means of the thioxoethenylidene (CCS) radical and ammonia (NH3) molecule observed with Karl G. Jansky Very Large Array telescope in D, C, and CNB configurations. Our main results are based on the CCS observation of the TMC-1C core, showing complex structures are present. Spatial offset relative to dust emission is observed in the CCS radical. Across a wide region around the dust peak, inward motion is found through the CCS radical. We have calculated the infall velocity and measured the turbulence inside the core. The turbulence is found to be subsonic. We obtain the virial parameter {\alpha} is < 1. Thus, thermal and non-thermal motions cannot prevent the collapse. Spatial incoherence of the CCS and NH3 is observed from the integrated intensity maps in these cores, suggesting that these molecules trace different environments of the cores.We compare the integrated flux densities of CCS with previous single-dish data and find that a small amount of flux is recovered in the interferometric observations, indicating the presence of significant diffuse emission in favorable conditions for producing CCS.

Black-hole neutron-star binary mergers, whose existence has been confirmed by gravitational-wave detectors, can lead to an electromagnetic counterpart called a kilonova if the neutron star is disrupted prior to merger. The observability of a kilonova depends crucially on the amount of neutron star ejecta, which is sensitive to the aligned component of the black hole spin. These binaries likely originate from the evolution of isolated stellar binaries. We explore the dependence of the ejected mass on two main mechanisms that provide high black hole spin. When the black hole inherits a high spin from a Wolf-Rayet star that was born with least $\sim 10\%$ of its breakup spin under weak stellar core-envelope coupling, which is relevant for all formation pathways, the median of the ejected mass is $\gtrsim 10^{-2}$ M$_{\odot}$. Though only possible for certain formation pathways, similarly large ejected mass results when the BH accretes $\gtrsim 20\%$ of its companion's envelope to gain a high spin, and a more massive stellar progenitor provides smaller ejected mass compared to when the black hole inherits high spin. Together, these signatures suggest that a population analysis of black hole masses and spins in black-hole neutron-star binary mergers may help distinguish between mechanisms for spin and possible formation pathways. Using a novel kilonova light curve model we show that current capabilities are unlikely to observe a counterpart, however future facilities such as the Vera Rubin Observatory will likely detect counterparts even if the aligned dimensionless spin of the disrupting black hole is as low as $\sim 0.2$. Our model predicts kilonovae as bright as $M_i \sim -14.5$ for an aligned black hole spin of $\sim 0.9$.

The star formation and metal enrichment histories of galaxies - at any epoch - constitute one of the key properties of galaxies, and their measurement is a core aim of observational extragalactic astronomy. The lack of deep rest-frame optical coverage at high-redshift has made robust constraints elusive, but this is now changing thanks to the \emph{James Webb Space Telescope (JWST)}. In preparation for the constraints provided by \emph{JWST} we explore the star formation and metal enrichment histories of galaxies at $z=5-13$ using the First Light And Reionisation Epoch Simulations (FLARES) suite. Built on the EAGLE model, the unique strategy of FLARES allows us to simulate a wide range of stellar masses (and luminosities) and environments. While we predict significant redshift evolution of average ages and specific star formation rates our core result is a mostly flat relationship of age and specific star formation rate with stellar mass. We also find that galaxies in this epoch predominantly have strongly rising star formation histories, albeit with the magnitude dropping with redshift and stellar mass. In terms of chemical enrichment we predict a strong stellar mass - metallicity relation present at $z=10$ and beyond alongside significant $\alpha$-enhancement. Finally, we find no environmental dependence of the relationship between age, specific star formation rate, or metallicity with stellar mass.

Ground-based high contrast exoplanet imaging requires state-of-the-art adaptive optics (AO) systems in order to detect extremely faint planets next to their brighter host stars. For such extreme AO systems (with high actuator count deformable mirrors over a small field of view), the lag time of the correction (which can impact our system by the amount the wavefront has changed by the time the system is able to apply the correction) which can be anywhere from ~1-5 milliseconds, can cause wavefront errors on spatial scales that lead to speckles at small angular separations from the central star in the final science image. One avenue for correcting these aberrations is predictive control, wherein previous wavefront information is used to predict the future state of the wavefront in one-system-lag's time, and this predicted state is applied as a correction with a deformable mirror. Here, we consider two methods for predictive control: data-driven prediction using empirical orthogonal functions and the physically-motivated predictive Fourier control. The performance and robustness of these methods have not previously been compared side-by-side. In this paper, we compare these predictors by applying them as post-facto methods to simulated atmospheres and on-sky telemetry, to investigate the circumstances in which their performance differs, including testing them under different wind speeds, C_n^2 profiles, and time lags. We also discuss future plans for testing both algorithms on the Santa Cruz Extreme AO Laboratory (SEAL) testbed.

We report a robust sample of 9 massive quiescent galaxies at redshift, $z > 3$, selected using the first data from the JWST CEERS programme. Three of these galaxies are at $4 < z < 5$, constituting the best evidence to date for quiescent galaxies significantly before $z=4$. These extreme galaxies have stellar masses in the range log$_{10}(M_*/$M$_\odot) = 10.5-11.3$, and formed the bulk of their mass at $6 < z < 9$, with two objects having star-formation histories that suggest they had already reached log$_{10}(M_*/$M$_\odot) > 10$ by $z\simeq8$. We report number densities for our sample, demonstrating that previous work underestimated the number of quiescent galaxies at $3 < z < 4$ by at least a factor of $3-6$, due to a lack of ultra-deep imaging data at $\lambda>2\,\mu$m. This result deepens the existing tension between observations and theoretical models, which already struggle to reproduce previous estimates of $z>3$ quiescent galaxy number densities. Upcoming wider-area JWST imaging surveys will provide larger samples of such galaxies, as well as providing opportunities to search for quiescent galaxies at $z>5$. The galaxies we report are excellent potential targets for JWST NIRSpec spectroscopy, which will be required to understand in detail their physical properties, providing deeper insights into the processes responsible for quenching star formation during the first billion years.

Studying the spatial distribution of extragalactic source populations is vital in understanding the matter distribution in the Universe. It also enables understanding the cosmological evolution of dark matter density fields and the relationship between dark matter and luminous matter. Clustering studies are also required for EoR foreground studies since it affects the relevant angular scales. This paper investigates the angular and spatial clustering properties and the bias parameter of radio-selected sources in the Lockman Hole field at 325 MHz. The data probes sources with fluxes $\gtrsim$0.3 mJy within a radius of 1.8$^\circ$ around the phase center of a $6^\circ \times 6^\circ$ mosaic. Based on their radio luminosity, the sources are classified into Active Galactic Nuclei (AGNs) and Star-Forming Galaxies (SFGs). Clustering and bias parameters are determined for the combined populations and the classified sources. The spatial correlation length and the bias of AGNs are greater than SFGs -- indicating that more massive haloes host the former. This study is the first reported estimate of the clustering property of sources at 325 MHz, intermediate between the preexisting studies at high and low-frequency bands. It also probes a well-studied deep field at an unexplored frequency with moderate depth and area. Clustering studies require such observations along different lines of sight, with various fields and data sets across frequencies to avoid cosmic variance and systematics. Thus, an extragalactic deep field has been studied in this work to contribute to this knowledge.

Much effort is devoted to measuring the nuclear symmetry energy through neutron star (NS) and nuclear observables. Since matter in the NS core may be non-hadronic, observables like radii and tidal deformability may not provide reliable constraints on properties of nucleonic matter. We demonstrate that coincident timing of a resonant shattering flare (RSF) and gravitational wave signal during binary NS inspiral probes the crust-core transition region and provides constraints on the symmetry energy comparable to terrestrial nuclear experiments. We show that nuclear masses, RSFs and measurements of NS radii and tidal deformabilities constrain different density ranges of the EOS, providing complementary probes.

The James Webb Space Telescope (JWST) will revolutionize the field of high-contrast imaging and enable both the direct detection of Saturn-mass planets and the characterization of substellar companions in the mid-infrared. While JWST will feature unprecedented sensitivity, angular resolution will be the key factor when competing with ground-based telescopes. Here, we aim to characterize the performance of several extreme angular resolution imaging techniques available with JWST in the 3-5 micron regime based on data taken during commissioning. Firstly, we introduce custom tools to simulate, reduce, and analyze NIRCam and MIRI coronagraphy data and use these tools to extract companion detection limits from on-sky NIRCam round and bar mask coronagraphy observations. Secondly, we present on-sky NIRISS aperture masking interferometry (AMI) and kernel phase imaging (KPI) observations from which we extract companion detection limits using the publicly available fouriever tool. Scaled to a total integration time of one hour and a target of the brightness of AB Dor, we find that NIRISS AMI and KPI reach contrasts of $\sim$7-8 mag at $\sim$70 mas and $\sim$9 mag at $\sim$200 mas. Beyond $\sim$250 mas, NIRCam coronagraphy reaches deeper contrasts of $\sim$13 mag at $\sim$500 mas and $\sim$15 mag at $\sim$2 arcsec. While the bar mask performs $\sim$1 mag better than the round mask at small angular separations $\lesssim$0.75 arcsec, it is the other way around at large angular separations $\gtrsim$1.5 arcsec. Moreover, the round mask gives access to the full 360 deg field-of-view which is beneficial for the search of new companions. We conclude that already during the instrument commissioning, JWST high-contrast imaging in the L- and M-bands performs close to its predicted limits.

In a cold and stable space environment, the James Webb Space Telescope (JWST or "Webb") reaches unprecedented sensitivities at wavelengths beyond 2 microns, serving most fields of astrophysics. It also extends the parameter space of high-contrast imaging in the near and mid-infrared. Launched in late 2021, JWST underwent a six month commissioning period. In this contribution we focus on the NIRCam Coronagraphy mode which was declared "science ready" on July 10 2022, the last of the 17 JWST observing modes. Essentially, this mode will allow to detect fainter/redder/colder (less massive for a given age) self-luminous exoplanets as well as other faint astrophysical signal in the vicinity of any bright object (stars or galaxies). Here we describe some of the steps and hurdles the commissioning team went through to achieve excellent performances. Specifically, we focus on the Coronagraphic Suppression Verification activity. We were able to produce firm detections at 3.35$\mu$m of the white dwarf companion HD 114174 B which is at a separation of $\simeq$ 0.5" and a contrast of $\simeq$ 10 magnitudes ($10^{4}$ fainter than the K$\sim$5.3 host star). We compare these first on-sky images with our latest, most informed and realistic end-to-end simulations through the same pipeline. Additionally we provide information on how we succeeded with the target acquisition with all five NIRCam focal plane masks and their four corresponding wedged Lyot stops.

Context. In the past few years, there has been a rise in the detection of streamers, asymmetric flows of material directed toward the protostellar disk with material from outside the star's natal core. It is unclear how they affect the process of mass accretion, in particular beyond the Class 0 phase. Aims. We investigate the gas kinematics around Per-emb-50, a Class I source in the crowded star-forming region NGC 1333. Our goal is to study how the mass infall proceeds from envelope to disk scales in this source. Results. We discover a streamer delivering material toward Per-emb-50 in H$_2$CO and C$^{18}$O emission. The streamer's emission can be well described by the analytic solutions for an infalling parcel of gas along a streamline with conserved angular momentum, both in the image plane and along the line of sight velocities. The streamer has a mean infall rate of $1.3 \times 10^{ -6}$ M$_{ \odot}$ yr$^{ -1}$, $5 -10$ times higher than the current accretion rate of the protostar. SO and SO$_2$ emission reveal asymmetric infall motions in the inner envelope, additional to the streamer around Per-emb-50. Furthermore, the presence of SO$_2$ could mark the impact zone of the infalling material. Conclusions. The streamer delivers sufficient mass to sustain the protostellar accretion rate and might produce an accretion burst, which would explain the protostar's high luminosity with respect to other Class I sources. Our results highlight the importance of late infall for protostellar evolution: streamers might provide a significant amount of mass for stellar accretion after the Class 0 phase.

Combining different observational probes, such as galaxy clustering and weak lensing, is a promising technique for unveiling the physics of the Universe with upcoming dark energy experiments. The galaxy redshift sample from the Dark Energy Spectroscopic Instrument (DESI) will have a significant overlap with major ongoing imaging surveys specifically designed for weak lensing measurements: the Kilo-Degree Survey (KiDS), the Dark Energy Survey (DES) and the Hyper Suprime-Cam (HSC) survey. In this work we analyse simulated redshift and lensing catalogues to establish a new strategy for combining high-quality cosmological imaging and spectroscopic data, in view of the first-year data assembly analysis of DESI. In a test case fitting for a reduced parameter set, we employ an optimal data compression scheme able to identify those aspects of the data that are most sensitive to the cosmological information, and amplify them with respect to other aspects of the data. We find this optimal compression approach is able to preserve all the information related to the growth of structure; we also extend this scheme to derive weights to be applied to individual galaxies, and show that these produce near-optimal results.

Searches for empirical clues beyond Einstein's general relativity (GR) are crucial to understand gravitation and spacetime. Radio pulsars have been playing an important role in testing gravity theories since 1970s. Because radio timing of binary pulsars is very sensitive to changes in the orbital dynamics, small deviations from what GR predicts can be captured or constrained. In this sense, the gravity sector in the standard-model extension was constrained tightly with a set of pulsar systems. Moreover, compact objects like pulsars are possible to develop nonperturbative deviations from GR in some specific alternative gravity theories, thus radio pulsars also provide rather unique testbeds in the strong-gravity regime.

The subject of this thesis is cosmological implications of string compactifications understood in a broad sense. In the first half of the thesis, we will begin by reviewing the four-dimensional description of the tree-level perturbative type IIB action. We will then review a number of open questions in cosmology and their relevance with regards to the remainder of the thesis. We will first explore some of these cosmological questions from the perspective of effective field theories motivated by supergravity. From the naturalness of dark energy and how to obtain a naturally light dark energy field in terms of the clockwork mechanism and the Dvali-Kaloper-Sorbo four-form mixing. We also discuss the coincidence problem for dynamical models of dark energy consistent with a quintessence field slowly rolling down a potential slope, of the type one would expect from the asymptotics of moduli space. In the second half of the thesis, we introduce the effects of perturbative and non-perturbative corrections to the tree-level type IIB action. We then focus on obtaining a viable model of quintessence from the type IIB effective field theory. However, we are able to show that such a model must have a non-supersymmetric Minkowski vacuum at leading order. When we consider the effects of quantum fluctuations during the early Universe, we see that such models must have extremely fine-tuned initial conditions to describe a slow-rolling scalar field at present times. We conclude that quintessence faces more challenges than a true cosmological constant, to the point that quintessence is very unattractive for model building modulo a ruling out of the cosmological constant by observations. Following this line of reasoning, we consider whether other perturbative corrections can generate de Sitter solutions in an appropriate setting.

Standard detection and analysis techniques for transient gravitational waves make the assumption that detector data contains, at most, one signal at any time. As detectors improve in sensitivity, this assumption will no longer be valid. In this paper we examine how current search techniques for transient gravitational waves will behave under the presence of more than one signal. We perform searches on data sets containing time-overlapping compact binary coalescences. This includes a modelled, matched filter search (PyCBC), and an unmodelled coherent search, coherent WaveBurst (cWB). Both of these searches are used by the LIGO-Virgo-KAGRA collaboration. We find that both searches are capable of identifying both signals correctly when the signals are dissimilar in merger time, $|\Delta t_c| \geq 1$ second, with PyCBC losing only $1\%$ of signals for overlapping binary black hole mergers. Both pipelines can find signal pairings within the region $|\Delta t_c| < 1$ second. However, clustering routines in the pipelines will cause only one of the two signals to be recovered, as such the efficiency is reduced. Within this region, we find that cWB can identify both signals. We also find that matched filter searches can be modified to provide estimates of the correct parameters for each signal.

Using a meta model for nuclear Equation of State (EOS) with its parameters constrained by astrophysical observations and terrestrial nuclear experiments, we examine effects of nuclear EOS especially its symmetry energy \esym term on the speed of sound squared $C^2_s(\rho)$ and the hadron-quark transition density $\rho_t$ where $C^2_s(\rho_t)$ vanishes via the spinodal decomposition mechanism in both dense neutron-rich nucleonic matter relevant for relativistic heavy-ion collisions and the $n+p+e+\mu$ matter in neutron stars at $\beta$-equilibrium. Unlike in nucleonic matter with fixed values of the isospin asymmetry $\delta$, in neutron stars with a density dependent isospin profile $\delta(\rho)$ determined consistently by the $\beta$ equilibrium and charge neutrality conditions, the $C^2_s(\rho)$ almost always show a peak and then vanishes at $\rho_t$ strongly depending on the high-density behavior of \esym if the skewness parameter $J_0$ characterizing the stiffness of high-density symmetric nuclear matter (SNM) EOS is not too far above its currently known most probable value of about $-190$ MeV inferred from recent Bayesian analyses of neutron star observables. Moreover, in the case of having a super-soft \esym that is decreasing with increasing density above about twice the saturation density of nuclear matter, the $\rho_t$ is significantly lower than the density where the \esym vanishes (indicating the onset of pure neutron matter formation) in neutron star cores.

The properties of hadron-quark hybrid stars are studied when the quark phase is described in terms of a local SU(3) Nambu--Jona-Lasinio (NJL) model taking into account the contribution of the vector and axial-vector interaction between the quarks, and the hadronic phase, in the relativistic mean field (RMF) model. For different values of the vector coupling constant $G_V$, the equations of state of the quark matter are calculated and the parameters of the hadron-quark phase transition are determined under the assumption that the phase transition takes place in accordance with Maxwell's construction. It is shown that for a larger vector coupling constant, the equation of state of the quark matter will be "stiffer" and the coexistence pressure $P_0$ of the phases will be greater. Using the resulting hybrid equations of state, the TOV equations are integrated numerically and the mass and radius of the compact star are determined for different values of the central pressure $P_c$. It is shown that when $G_V$ is larger, the maximum mass of the compact star will be larger and thereby, the radius of the configuration with maximum mass will be smaller. Questions of the stability of hybrid stars are also discussed. It is shown that in terms of the model examined here, for all values of the vector coupling constant, a hybrid star with an infinitely small quark core is stable. These results are compared with recent measurements of the mass and radius of the pulsars PSR J0030+0451 and PSR J0740+6620, carried out at the International Space Station with the NASA's Neutron star Interior Composition Explorer (NICER) X-ray telescope. A comparison of the theoretical results with observational data does not exclude the possibility of quark deconfinement in the interiors of compact stars.

A particle moving on curved magnetic field lines in the wormhole metrics is considered. The dynamics of a single particle is studied on a co-rotating parameterized magnetic field line. Being interested in the force-free dynamics we choose the rotating Archimedes' spiral embedded in the generalized Ellis-Bronnikov metrics. By analysing the phase space, we explore how the physical parameters affect the motion of the particle and it is shown that the particles may reach the force free regime, eventually escaping the wormhole.

Electron ring velocity space distributions have previously been seen in numerical simulations of magnetic reconnection exhausts and have been suggested to be caused by the magnetization of the electron outflow jet by the compressed reconnected magnetic fields [Shuster et al., ${\it Geophys.~Res.~Lett.}, {\bf 41}$, 5389 (2014)]. We present a theory of the dependence of the major and minor radii of the ring distributions solely in terms of upstream (lobe) plasma conditions, thereby allowing a prediction of the associated temperature and temperature anisotropy of the rings in terms of upstream parameters. We test the validity of the prediction using 2.5-dimensional particle-in-cell (PIC) simulations with varying upstream plasma density and temperature, finding excellent agreement between the predicted and simulated values. We confirm the Shuster et al. suggestion for the cause of the ring distributions, and also find that the ring distributions are located in a region marked by a plateau, or shoulder, in the reconnected magnetic field profile. The predictions of the temperature are consistent with observed electron temperatures in dipolarization fronts, and may provide an explanation for the generation of plasma with temperatures in the 10s of MK in super-hot solar flares. A possible extension of the model to dayside reconnection is discussed. Since ring distributions are known to excite whistler waves, the present results should be useful for quantifying the generation of whistler waves in reconnection exhausts.

Black holes hold a tremendous discovery potential. In this paper the extent to which the Event Horizon Telescope and its next generation upgrade can resolve their structure is quantified. Black holes are characterized by a perfectly absorptive boundary, with a specific area determined by intrinsic parameters of the black hole. We use a general parametrization of spherically symmetric spacetimes describing deviations from this behavior, with parameters controlling the size of the central object and its interaction with light, in particular through a specular reflection coefficient $\Gamma$ and an intrinsic luminosity measured as a fraction $\eta$ of that of the accretion disc. This enables us to study exotic compact objects and compare them with black holes in a model-independent manner. We determine the image features associated with the existence of a surface in the presence of a geometrically thin and optically thick accretion disc, identifying requirements for VLBI observations to be able to cast meaningful constraints on these parameters, in particular regarding angular resolution and image dynamic range. For face-on observations, constraints of order $\eta\lesssim 10^{-4}, \Gamma\lesssim 10^{-1}$ are possible with an enhanced EHT array, imposing strong constraints on the nature of the central object.

The traditional description of the orbital evolution of compact-object binaries, like double white dwarfs (DWDs), assumes that the system is driven only by gravitational wave (GW) radiation. However, the high magnetic fields with intensities of up to gigagauss measured in WDs alert a potential role of the electromagnetic (EM) emission in the evolution of DWDs. We evaluate the orbital dynamics of DWDs under the effects of GW radiation, tidal synchronization, and EM emission by a unipolar inductor generated by the magnetic primary and the relative motion of the non-magnetic secondary. We show that the EM emission can affect the orbital dynamics for magnetic fields larger than megagauss. We applied the model to two known DWDs, SDSS J0651+2844 and ZTF J1539+5027, for which the GW radiation alone does not fully account for the measured orbital decay rate. We obtain upper limits to the primary's magnetic field strength, over which the EM emission causes an orbital decay faster than observed. The contribution of tidal locking and the EM emission is comparable, and together they can contribute up to $20\%$ to the measured orbital decay rate. We show that the gravitational waveform for a DWD modeled as purely driven by GWs and including tidal interactions and EM emission can have large relative dephasing detectable in the mHz regime of frequencies relevant for space-based detectors like LISA. Therefore, including physics besides GW radiation in the waveform templates is essential to calibrate the GW detectors using known sources, e.g., ZTF J1539+5027, and to infer binary parameters.

The impact of exposure to astrophysical ionizing radiation on astronaut health is one of the main concerns in planning crewed missions to Mars. Astronauts will be exposed to energetic charged particles from Galactic and Solar origin for a prolonged period with little protection from a thin spacecraft shield in transit and from the rarefied Martian atmosphere when on the surface. Adverse impacts on astronaut health include, for example, Acute Radiation Syndrome, damage to the nervous system, and increased cancer risk. Using a combination of radiation measurements and numerical modeling with the GEANT4 package, we calculate the distribution of radiation dose in various organs of the human body for various expected scenarios simulated with a model human phantom. We rely on medical studies to assess the impact of enhanced levels of radiation dose on various physiological systems and on the overall health of astronauts. We suggest mitigation strategies such as improved ways of shielding and dietary supplements, and make recommendations for the safety of astronauts in future crewed missions to Mars.