Papers for Thursday, Apr 20 2023

We report the detection and validation of two planets orbiting TOI-2095 (TIC 235678745). The host star is a 3700K M1V dwarf with a high proper motion. The star lies at a distance of 42 pc in a sparsely populated portion of the sky and is bright in the infrared (K=9). With data from 24 Sectors of observation during TESS's Cycles 2 and 4, TOI-2095 exhibits two sets of transits associated with super-Earth-sized planets. The planets have orbital periods of 17.7 days and 28.2 days and radii of 1.30 and 1.39 Earth radii, respectively. Archival data, preliminary follow-up observations, and vetting analyses support the planetary interpretation of the detected transit signals. The pair of planets have estimated equilibrium temperatures of approximately 400 K, with stellar insolations of 3.23 and 1.73 times that of Earth, placing them in the Venus zone. The planets also lie in a radius regime signaling the transition between rock-dominated and volatile-rich compositions. They are thus prime targets for follow-up mass measurements to better understand the properties of warm, transition radius planets. The relatively long orbital periods of these two planets provide crucial data that can help shed light on the processes that shape the composition of small planets orbiting M dwarfs.

Surveys with a narrow field-of-view can play an important role in probing cosmology, but inferences from these surveys suffer from large sample variance. The standard method for computing the sample variance is based on two key approximations, and we demonstrate that it can lead to a significant underestimate of the sample variance in narrow surveys. We present a new method for accurately computing the sample variance and apply our method to the recent observations of the warm-hot intergalactic medium (WHIM) based on spectroscopic measurements of blazars. We find that the sample variances in these surveys are significantly larger than the quoted measurement errors; for example, the cosmic mean baryon density contained in the WHIM could be lower by $54\%$ at $1\text{-}\sigma$ fluctuation than estimated in one observation. Accurately quantifying the sample variance is essential in deriving correct interpretations of the measurements in surveys with a small field-of-view.

We present a statistical analysis of the large-scale (up to 2 Mpc) environment of an homogeneous and complete sample, both in radio and optical selection, of wide-angle tailed radio galaxies (WATs) in the local Universe (i.e., with redshifts $z\lesssim$ 0.15). The analysis is carried out using the parameters obtained from cosmological neighbors within 2 Mpc of the target source. Results on WATs large-scale environments are then compared with that of Fanaroff-Riley type I (FR Is) and type II (FR IIs) radio galaxies, listed in two others homogeneous and complete catalogs, and selected with the same criterion adopted for the WATs catalog. We obtain indication that at low redshift WATs inhabit environments with a larger number of galaxies than that of FR Is and FR IIs. In the explored redshift range, the physical size of the galaxy group/cluster in which WATs reside appears to be almost constant with respect to FR Is and FR IIs, being around 1 Mpc. From the distribution of the concentration parameter, defined as the ratio between the number of cosmological neighbors lying within 500 kpc and within 1 Mpc, we conclude that WATs tend to inhabit the central region of the group/cluster in which they reside, in agreement with the general paradigm that WATs are the cluster BCG.

We investigate the accuracy requirements for field-level inference of cluster masses and void sizes using data from galaxy surveys. We introduce a two-step framework that takes advantage of the fact that cluster masses are determined by flows on larger scales than the clusters themselves. First, we determine the integration accuracy required to perform field-level inference of cosmic initial conditions on these large scales, by fitting to late-time galaxy counts using the Bayesian Origin Reconstruction from Galaxies (BORG) algorithm. A 20-step COLA integrator is able to accurately describe the density field surrounding the most massive clusters in the Local Super-Volume ($<135\,h^{-1}\mathrm{Mpc}$), but does not by itself lead to converged virial mass estimates. Therefore we carry out `posterior resimulations', using full $N$-body dynamics while sampling from the inferred initial conditions, and thereby obtain estimates of masses for nearby massive clusters. We show that these are in broad agreement with existing estimates, and find that mass functions in the Local Super-Volume are compatible with $\Lambda$CDM.

Jellyfish galaxies are the prototypical examples of satellite galaxies undergoing strong ram pressure stripping (RPS). We analyze the evolution of 535 unique, first-infalling jellyfish galaxies from the TNG50 cosmological hydrodynamical galaxy simulation. These have been visually inspected to be undergoing RPS sometime in the past 5 billion years (since $z=0.5$), have satellite stellar masses $M_{\star}^{\rm sat}\sim10^{8-10.5}{\rm M}_\odot$, and live in hosts with $M_{\rm 200c}\sim10^{12-14.3}{\rm M}_\odot$ at $z=0$. We quantify the cold gas $(\leq10^{4.5}$K) removal using the tracer particles, confirming that for these jellyfish, RPS is the dominant driver of cold gas loss after infall. Half of these jellyfish are completely devoid of cold gas by $z=0$, and these galaxies have earlier infall times and smaller satellite-to-host mass ratios than those that still have some gas today. RPS can act on jellyfish galaxies over long timescales of $\approx1.5-8$Gyr. Jellyfish in more massive hosts are impacted by RPS for a shorter time span and, at a fixed host halo mass, jellyfish with lower stellar masses at $z=0$ have shorter RPS time spans. While RPS may act for long periods of time, the peak RPS period -- where at least 50% of the total RPS occurs -- begins within $\approx1$Gyr of infall and lasts $\lesssim2$Gyr. During this period, the jellyfish are at host-centric distances between $\sim0.2-2R_{\rm 200c}$, illustrating that much of RPS occurs at large distances from the host galaxy. Jellyfish continue forming stars until they have lost $\approx98$% of their cold gas. For groups and clusters in TNG50 $(M_{\rm 200c}\sim10^{13-14.3}{\rm M}_\odot)$, jellyfish galaxies deposit more cold gas ($\sim10^{11-12}{\rm M}_\odot$) into halos than exist in them at $z=0$, demonstrating that jellyfish, and in general satellite galaxies, are a significant source of cold gas accretion.

We develop a framework for self-consistently extracting cosmological information from the clustering of tracers in redshift space, $\textit{without}$ relying on model-dependent templates to describe the baryon acoustic oscillation (BAO) feature. Our approach uses the recently proposed Laguerre reconstruction technique for the BAO feature and its linear point $r_{\rm LP}$, and substantially extends it to simultaneously model the multipoles $\ell=0,2,4$ of the anisotropic galaxy 2-point correlation function (2pcf). The approach is `model-agnostic': it assumes that the non-linear growth of structure smears the BAO feature by an approximately Gaussian kernel with a smearing scale $\sigma_{\rm v}$, but does not assume any fiducial cosmology for describing the shape of the feature itself. Using mock observations for two realistic survey configurations assuming $\Lambda$ cold dark matter ($\Lambda$CDM), combined with Bayesian parameter inference, we show that the linear point $r_{\rm LP}$ and smearing scale $\sigma_{\rm v}$ can be accurately recovered by our method in both existing and upcoming surveys. The precision of the recovery of $r_{\rm LP}$ is always better than $1\%$, while $\sigma_{\rm v}$ can be recovered with $\lesssim10\%$ uncertainty provided the linear galaxy bias $b$ is separately constrained, e.g., using weak lensing observations. Our method is also sensitive to the linear growth rate $f$, albeit with larger uncertainties and systematic errors, especially for upcoming surveys such as DESI. We discuss how our model can be modified to improve the recovery of $f$, such that the resulting constraints on $\{f,\sigma_{\rm v},r_{\rm LP}\}$ can potentially be used as a test of cosmological models including and beyond $\Lambda$CDM.

Due to ram-pressure stripping, jellyfish galaxies are thought to lose large amounts, if not all, of their interstellar medium. Nevertheless, some, but not all, observations suggest that jellyfish galaxies exhibit enhanced star formation compared to control samples, even in their ram pressure-stripped tails. We use the TNG50 cosmological gravity+magnetohydrodynamical simulation, with an average spatial resolution of 50-200 pc in the star-forming regions of galaxies, to quantify the star formation activity and rates (SFRs) of more than 700 jellyfish galaxies at $z=0-1$ with stellar masses $10^{8.3-11}\,\mathrm{M}_\odot$ in hosts with mass $10^{11.5-14.3}\,\mathrm{M}_\odot$. We extract their global SFRs, the SFRs within their main stellar body vs. within the tails, and we follow the evolution of the star formation along their individual evolutionary tracks. We compare the findings for jellyfish galaxies to those of diversely-constructed control samples, including against satellite and field galaxies with matched redshift, stellar mass, gas fraction and host halo mass. According to TNG50, star formation and ram-pressure stripping can indeed occur simultaneously within any given galaxy, and frequently do so. Moreover, star formation can also take place within the ram pressure-stripped tails, even though the latter is typically subdominant. However, TNG50 does not predict elevated population-wide SFRs in jellyfish compared to analog satellite galaxies with the same stellar mass or gas fraction. Simulated jellyfish galaxies do undergo bursts of elevated star formation along their history but, at least according to TNG50, these do not translate into a population-wide enhancement at any given epoch.

We present the ``Cosmological Jellyfish'' project - a citizen-science classification program to identify jellyfish galaxies within the IllustrisTNG cosmological simulations. Jellyfish (JF) are satellite galaxies that exhibit long trailing gas features -- `tails' -- extending from their stellar body. Their distinctive morphology arises due to ram-pressure stripping (RPS) as they move through the background gaseous medium. Using the TNG50 and TNG100 simulations, we construct a sample of $\sim 80,000$ satellite galaxies spanning an unprecedented range of stellar masses, $10^{8.3-12.3}\,\mathrm{M_\odot}$, and host masses of $M_\mathrm{200,c}=10^{10.4-14.6}\,\mathrm{M_\odot}$ back to $z=2$ \citep[extending the work of][]{yun_jellyfish_2019}. Based on this sample, $\sim 90,000$ galaxy images were presented to volunteers in a citizen-science project on the Zooniverse platform who were asked to determine if each galaxy image resembles a jellyfish. Based on volunteer votes, each galaxy was assigned a score determining if it is a JF or not. This paper describes the project, the inspected satellite sample, the methodology, and the classification process that resulted in a dataset of $5,307$ visually-identified jellyfish galaxies. We find that JF galaxies are common in nearly all group- and cluster-sized systems, with the JF fraction increasing with host mass and decreasing with satellite stellar mass. We highlight JF galaxies in three relatively unexplored regimes: low-mass hosts of $M_\mathrm{200,c}\sim10^{11.5-13}\,\mathrm{M_\odot}$, radial positions within hosts exceeding the virial radius $R_\mathrm{200,c}$, and at high redshift up to $z=2$. The full dataset of our jellyfish scores is publicly available and can be used to select and study JF galaxies in the IllustrisTNG simulations.

The evolution of binaries that become double white dwarf (DWD) can cause the ejection of high amounts of dust and gas. Such material can give rise to circumbinary discs and become the cradle of new planets, yet no studies so far have focused on the formation of circumbinary planets around DWDs. These binaries will be the main sources of gravitational waves (GWs) detectable by the ESA Laser Interferometer Space Antenna (LISA) mission, opening the possibility to detect circumbinary planets around short-period DWDs everywhere in the Milky Way. We investigate the formation of Magrathea planets by simulating multiple planet formation tracks to explore how seeds growing first by pebble accretion, and then by gas accretion, are affected by the disc environments surrounding DWDs. We present both planetary formation tracks taking place in steady-state discs, and formation tracks taking place in discs evolving with time. The time-dependent tracks account for both the disc accretion rate onto the central binary and the disc photoevaporation rate caused by stellar irradiation. Our results show that planetary formation in circumbinary discs around DWDs can be possible. In particular, the extreme planetary formation environment implies three main significant results: (i) the accretion rate and the metallicity of the disc should be high in order to form sub-stellar objects with masses up to 31 M$_J$, this is achieved only if planet formation starts soon after the onset of the disc and if first generation seeds are present in the disc; (ii) seeds formed within 0.1 Myr, or within 1 Myr, from the onset of the disc can only produce sub-Neptune and Neptunian planets, unless the disc accommodates first generation seeds with mass 10 M$_{\oplus}$; (iii) most of the planets are finally located within 1 au from the disc centre, while they are still undergoing the gas accretion phase.

We present Atacama Large Millimeter/Submillimeter Array (ALMA) [CII] and $\sim158$ $\rm\mu m$ continuum observations of REBELS-25, a massive, morphologically complex ultra-luminous infrared galaxy (ULIRG; $L_{\rm IR}=1.5^{+0.8}_{-0.5}\times10^{12}$ L$_\odot$) at $z=7.31$, spectroscopically confirmed by the Reionization Era Bright Emission Line Survey (REBELS) ALMA Large Programme. REBELS-25 has a significant stellar mass of $M_{*}=8^{+4}_{-2}\times10^{9}$ M$_\odot$. From dust-continuum and ultraviolet observations, we determine a total obscured + unobscured star formation rate of SFR $=199^{+101}_{-63}$ M$_\odot$ yr$^{-1}$. This is about four times the SFR estimated from an extrapolated main-sequence. We also infer a [CII]-based molecular gas mass of $M_{\rm H_2}=5.1^{+5.1}_{-2.6}\times10^{10}$ $M_\odot$, implying a molecular gas depletion time of $ t_{\rm depl, H_2}=0.3^{+0.3}_{-0.2}$ Gyr. We observe a [CII] velocity gradient consistent with disc rotation, but given the current resolution we cannot rule out a more complex velocity structure such as a merger. The spectrum exhibits excess [CII] emission at large positive velocities ($\sim500$ km s$^{-1}$), which we interpret as either a merging companion or an outflow. In the outflow scenario, we derive a lower limit of the mass outflow rate of 200 M$_\odot$ yr$^{-1}$, which is consistent with expectations for a star formation-driven outflow. Given its large stellar mass, SFR and molecular gas reservoir $\sim700$ Myr after the Big Bang, we explore the future evolution of REBELS-25. Considering a simple, conservative model assuming an exponentially declining star formation history, constant star formation efficiency, and no additional gas inflow, we find that REBELS-25 has the potential to evolve into a galaxy consistent with the properties of high-mass quiescent galaxies recently observed at $z\sim4$.

Co-orbital systems contain two or more bodies sharing the same orbit around a planet or star. The best-known flavors of co-orbital systems are tadpoles (in which two bodies' angular separations oscillate about the L4/L5 Lagrange points $60^\circ$ apart) and horseshoes (with two bodies periodically exchanging orbital energy to trace out a horseshoe shape in a co-rotating frame). Here, we use N-body simulations to explore the parameter space of many-planet horseshoe systems. We show that up to 24 equal-mass, Earth-mass planets can share the same orbit at 1 au, following a complex pattern in which neighboring planets undergo horseshoe oscillations. We explore the dynamics of horseshoe constellations, and show that they can remain stable for billions of years and even persist through their stars' post-main sequence evolution. With sufficient observations, they can be identified through their large-amplitude, correlated transit timing variations. Given their longevity and exotic orbital architectures, horseshoe constellations may represent potential SETI beacons.

In co-orbital planetary systems, two or more planets share the same orbit around their star. Here we test the dynamical stability of co-orbital rings of planets perturbed by outside forces. We test two setups: i) 'stationary' rings of planets that, when unperturbed, remain equally-spaced along their orbit; and ii) horseshoe constellation systems, in which planets are continually undergoing horseshoe librations with their immediate neighbors. We show that a single rogue planet crossing the planets' orbit more massive than a few lunar masses (0.01-0.04 Earth masses) systematically disrupts a co-orbital ring of 6, 9, 18, or 42 Earth-mass planets located at 1 au. Stationary rings are more resistant to perturbations than horseshoe constellations, yet when perturbed they can transform into stable horseshoe constellation systems. Given sufficient time, any co-orbital ring system will be perturbed into either becoming a horseshoe constellation or complete destabilization.

We present a review of Saturn's interior structure and thermal evolution, with a particular focus on work in the past 5 years. Data from the Cassini mission, including a precise determination of the gravity field from the Grand Finale orbits, and the still ongoing identification of ring wave features in Saturn's C-ring tied to seismic modes in the planet, have led to dramatic advances in our understanding of Saturn's structure. Models that match the gravity field suggest that differential rotation, as seen in the visible atmosphere, extends down to at least a depth of 10,000 km (1/6$^{\rm th}$ the planet's radius). At greater depths, a variety of different investigations all now point to a deep Saturn rotation rate of 10 hours and 33 minutes. There is very compelling evidence for a central heavy element concentration (``core''), that in most recent models is 12-20 Earth masses. Ring seismology strongly suggests that the core is not entirely compact, but is dilute (mixed in with the overlying H/He), and has a substantial radial extent, perhaps out to around one-half of the planet's radius. A wide range of thermal evolution scenarios can match the planet's current luminosity, with progress on better quantifying the helium rain scenario hampered by Saturn's poorly known atmospheric helium abundance. We discuss the relevance of magnetic field data on understanding the planet's current interior structure. We point towards additional future work that combines seismology and gravity within a framework that includes differential rotation, and the utility of a Saturn entry probe.

A blue depression is found in the spectra of M dwarfs from 4000 to 4500A. This depression shows an increase toward lower temperatures though is particularly sensitive to gravity and metallicity. It is the single most sensitive feature in the optical spectra of M dwarfs. The depression appears as centered on the neutral calcium resonance line at 4227A and leads to nearby features being weaker by about two orders of magnitude than predicted. We consider a variety of possible causes for the depression including temperature, gravity, metallicity, dust, damping constants, and atmospheric stratification. We also consider relevant molecular opacities which might be the cause identifying AlH, SiH, and NaH in the spectral region. However, none of these solutions are satisfactory. In the absence of a more accurate determination of the broadening of the calcium line perturbed by molecular hydrogen, we find a promising empirical fit using a modified Lorentzian line profile for the calcium resonance line. Such fits provide a simplistic line-broadening description for this calcium resonance line and potentially other un-modelled resonance lines in cool high-pressure atmospheres. Thus we claim the most plausible cause of the blue depression in the optical spectra of M dwarfs is a lack of appropriate treatment of line broadening for atomic calcium. The broad wings of the calcium resonance line develop at temperatures below about 4000K and are analogous to the neutral sodium and potassium features which dominate the red optical spectra of L dwarfs.

The main scientific goal of TESS is to find planets smaller than Neptune around stars bright enough to allow further characterization studies. Given our current instrumentation and detection biases, M dwarfs are prime targets to search for small planets that are in (or nearby) the habitable zone of their host star. Here we use photometric observations and CARMENES radial velocity measurements to validate a pair of transiting planet candidates found by TESS. The data was fitted simultaneously using a Bayesian MCMC procedure taking into account the stellar variability present in the photometric and spectroscopic time series. We confirm the planetary origin of the two transiting candidates orbiting around TOI-2095 (TIC 235678745). The star is a nearby M dwarf ($d = 41.90 \pm 0.03$ pc, $T_{\rm eff} = 3759 \pm 87$ K, $V = 12.6$ mag) with a stellar mass and radius of $M_\star = 0.44 \pm 0.02 \; M_\odot$ and $R_\star = 0.44 \pm 0.02 \; R_\odot$, respectively. The planetary system is composed of two transiting planets: TOI-2095b with an orbital period of $P_b = 17.66484 \pm (7\times 10^{-5})$ days and TOI-2095c with $P_c = 28.17232 \pm (14\times 10^{-5})$ days. Both planets have similar sizes with $R_b = 1.25 \pm 0.07 \; R_\oplus$ and $R_c = 1.33 \pm 0.08 \; R_\oplus$ for planet b and c, respectively. We put upper limits on the masses of these objects with $M_b < 4.1 \; M_\oplus$ for the inner and $M_c < 7.4 \; M_\oplus$ for the outer planet (95\% confidence level). These two planets present equilibrium temperatures in the range of 300 - 350 K and are close to the inner edge of the habitable zone of their star.

We present a constraint on the abundance of supergiant (SG) stars at redshift z approx. 1, based on recent observations of a strongly lensed arc at this redshift. First we derive a free-form model of MACS J0416.1-2403 using data from the BUFFALO program. The new lens model is based on 72 multiply lensed galaxies that produce 214 multiple images, making it the largest sample of spectroscopically confirmed lensed galaxies on this cluster. The larger coverage in BUFFALO allows us to measure the shear up to the outskirts of the cluster, and extend the range of lensing constraints up to ~ 1 Mpc from the central region, providing a mass estimate up to this radius. As an application, we make predictions for the number of high-redshift multiply-lensed galaxies detected in future observations with JWST. Then we focus on a previously known lensed galaxy at z=1.0054, nicknamed Spock, which contains four previously reported transients. We interpret these transients as microcaustic crossings of SG stars and compute the probability of such events. Based on simplifications regarding the stellar evolution, we find that microlensing (by stars in the intracluster medium) of SG stars at z=1.0054 can fully explain these events. The inferred abundance of SG stars is consistent with either (1) a number density of stars with bolometric luminosities beyond the Humphreys-Davidson (HD) limit (L ~ $6\times10^5 L_{\odot}$) that is below 400 stars per sq. kpc, or (2) the absence of stars beyond the HD limit but with a SG number density of ~ 9000 per sq. kpc for stars with luminosities between $10^5$ and $6\times10^5$. This is equivalent to one SG star per 10x10 pc$^2$. We finally make predictions for future observations with JWST's NIRcam. We find that in observations made with the F200W filter that reach 29 mag AB, if cool red SG stars exist at z~1 beyond the HD limit, they should be easily detected in this arc

Quaoar is a classical Trans-Neptunian Object (TNO) with an area equivalent diameter of 1,100 km and an orbital semi-major axis of 43.3 astronomical units. Based on stellar occultations observed between 2018 and 2021, an inhomogeneous ring (Q1R, Quaoar's first ring) was detected around this body. Aims. A new stellar occultation by Quaoar was observed on August 9th, 2022 aiming to improve Quaoar's shape models and the physical parameters of Q1R while searching for additional material around the body. Methods. The occultation provided nine effective chords across Quaoar, pinning down its size, shape, and astrometric position. Large facilities, such as Gemini North and the Canada-France-Hawaii Telescope (CFHT), were used to obtain high acquisition rates and signal-to-noise ratios. The light curves were also used to characterize the Q1R ring (radial profiles and orbital elements). Results. Quaoar's elliptical fit to the occultation chords yields the limb with an apparent semi-major axis of $579.5\pm4.0$ km, apparent oblateness of $0.12\pm0.01$, and area-equivalent radius of $543\pm2$ km. Quaoar's limb orientation is consistent with Q1R and Weywot orbiting in Quaoar's equatorial plane. The orbital radius of Q1R is refined to a value of $4,057\pm6$ km. The radial opacity profile of the more opaque ring profile follows a Lorentzian shape that extends over 60 km, with a full width at half maximum (FWHM) of $\sim5$ km and a peak normal optical depth of 0.4. Besides the secondary events related to the already reported rings, new secondary events detected during the August 2022 occultation in three different data sets are consistent with another ring around Quaoar with a radius of $2,520\pm20$ km, assuming the ring is circular and co-planar with Q1R. This new ring has a typical width of 10 km and a normal optical depth of $\sim$0.004. Like Q1R, it also lies outside Quaoar's classical Roche limit.

We test the impact of gravitational lensing on the lifetime estimates of seven high-redshift quasars at redshift $z\gtrsim6$. The targeted quasars are identified by their small observed proximity zone sizes, which indicate extremely short quasar lifetimes $(t_Q\lesssim10^5 \text{ yrs})$. However, these estimates of quasar lifetimes rely on the assumption that the observed luminosities of the quasars are intrinsic and not magnified by gravitational lensing, which would bias the lifetime estimates towards younger ages. In order to test possible effects of gravitational lensing, we obtain high-resolution images of the seven quasars with the {\em Hubble Space Telescope (HST)} and look for signs of strong lensing. We do not find any evidence of strong lensing, i.e., all quasars are well-described by point sources, and no foreground lensing galaxy is detected. We estimate that the strong lensing probabilities for these quasars are extremely small $(\sim1.4\times10^{-5})$, and show that weak lensing changes the estimated quasar lifetimes by only $\lesssim0.2$ dex. We thus confirm that the short lifetimes of these quasars are intrinsic. The existence of young quasars indicates a high obscured fraction, radiatively inefficient accretion, and/or flickering light curves for high-redshift quasars. We further discuss the impact of lensing magnification on measurements of black hole masses and Eddington ratios of quasars.

The Atacama Large Millimeter/Submillimeter Array (ALMA) in the sub-millimeter and the James Webb Space Telescope (JWST) in the infrared have achieved robust spectroscopic detections of emission lines from the interstellar medium (ISM) in some of the first galaxies. These unprecedented measurements provide valuable information regarding the ISM properties, stellar populations, galaxy morphologies, and kinematics in these high-redshift galaxies and, in principle, offer powerful tests of state-of-the-art galaxy formation models, as implemented in hydrodynamical simulations. To facilitate direct comparisons between simulations and observations, we develop a fast post-processing pipeline for predicting the line emission from the HII regions around simulated star particles, accounting for spatial variations in the surrounding gas density, metallicity, temperature, and incident radiation spectrum. Our ISM line emission model currently captures H$\alpha$, H$\beta$, and all of the [OIII] and [OII] lines targeted by ALMA and the JWST at $z>6$. We illustrate the power of this approach by applying our line emission model to the publicly available FIRE high-$z$ simulation suite and perform a detailed comparison with current observations. We show that the FIRE mass--metallicity relation is in $1\sigma$ agreement with ALMA/JWST measurements after accounting for the inhomogeneities in ISM properties. We also quantitatively validate the one-zone model description, which is widely used for interpreting [OIII] and H$\beta$ line luminosity measurements. This model is publicly available and can be implemented on top of a broad range of galaxy formation simulations for comparison with JWST and ALMA measurements.

We present a study of the $\rm{[OIII]\lambda\,5007}$ line profile in a sub-sample of 8 active galactic nuclei (AGN) and 6 non-AGN in the optically-selected green valley at $\rm{z\,<\,0.5}$ using long-slit spectroscopic observations with the 11 m Southern African Large Telescope. Gaussian decomposition of the line profile was performed to study its different components. We observe that the AGN profile is more complex than the non-AGN one. In particular, in most AGN (5/8) we detect a blue wing of the line. We derive the FWHM velocities of the wing and systemic component, and find that AGN show higher FWHM velocity than non-AGN in their core component. We also find that the AGN show blue wings with a median velocity width of approximately 600 $\rm{km\,s^{-1}}$, and a velocity offset from the core component in the range -90 to -350 $\rm{km\,s^{-1}}$, in contrast to the non-AGN galaxies, where we do not detect blue wings in any of their $\rm{[OIII]\lambda\,5007}$ line profiles. Using spatial information in our spectra, we show that at least three of the outflow candidate galaxies have centrally driven gas outflows extending across the whole galaxy. Moreover, these are also the galaxies which are located on the main sequence of star formation, raising the possibility that the AGN in our sample are influencing SF of their host galaxies (such as positive feedback). This is in agreement with our previous work where we studied SF, morphology, and stellar population properties of a sample of green valley AGN and non-AGN galaxies.

The high mass X-ray binary system 4U 0114+65 was observed by Nustar in October 2019, and by XMM-Newton in August 2015. Here we performed spectral and timing analysis of the Nustar observation, and carry out timing analysis on the XMM-Newton data. We measured the spin period of the neutron star from both observations and found a spin-up rate $\dot{p} = 1.54 \pm 0.38 \times 10^{-6} s s^{-1}$. During the Nustar observation two flares occured, one occured shortly after the start of the observation and the other near the end separated by a long period of low/quiescent- state. The large and sudden flares mostly resulted from accretion of Corotating Interaction Region (CIR) material. A common spectral model to HMXBs, powerlaw with high energy cutoff and absorption at low energy, gave a good fit to both flaring and quiescent states. A flourescent iron line was not required in fitting any of the states. On the other hand, very tentative evidence of Cyclotron Resonant Scattering Feature (CRSF) at $\sim$ 17 keV was found during fitting using cyclabs model, however fitting improvement was not significant enough to confirm its detection, plus a very narrow width (< 1 keV) was obtained for the line and its first harmonic. Visual inspection of the spectra showed a deficiency of emission near the expected first and second harmonic. Another important feature visually noticed in the spectra is the presence of hard tail above 50 keV. This could be explained by the shocked material bounding the CIR.

We summarize the collective knowledge of physical and surface properties of comet nuclei, focusing on those that are obtained from remote observations. We now have measurements or constraints on effective radius for over 200 comets, rotation periods for over 60, axial ratios and color indices for over 50, geometric albedos for over 25, and nucleus phase coefficients for over 20. The sample has approximately tripled since the publication of Comets II, with IR surveys using Spitzer and NEOWISE responsible for the bulk of the increase in effective radii measurements. Advances in coma morphology studies and long-term studies of a few prominent comets have resulted in meaningful constraints on rotation period changes in nearly a dozen comets, allowing this to be added to the range of nucleus properties studied. The first delay-Doppler radar and visible light polarimetric measurements of comet nuclei have been made since Comets II and are considered alongside the traditional methods of studying nuclei remotely. We use the results from recent in situ missions, notably Rosetta, to put the collective properties obtained by remote observations into context, emphasizing the insights gained into surface properties and the prevalence of highly elongated and/or bilobate shapes. We also explore how nucleus properties evolve, focusing on fragmentation and the likely related phenomena of outbursts and disintegration. Knowledge of these behaviors has been shaped in recent years by diverse sources: high resolution images of nucleus fragmentation and disruption events, the detection of thousands of small comets near the Sun, regular photometric monitoring of large numbers of comets throughout the solar system, and detailed imaging of the surfaces of mission targets. Finally, we explore what advances in the knowledge of the bulk nucleus properties may be enabled in coming years.

The most massive black holes at redshifts z = 6 were already over billion solar masses. In this chapter, we discuss the formation and growth of the first black holes in the Universe. The deaths of massive primordial stars provide potential seeds of supermassive black holes. Theoretical models predict that the seed black hole masses range from 10 to 100,000 solar masses. Their initial fueling may be limited by feedback from its progenitor star, the black hole itself, and nearby star formation. Once the halo and galaxy surpasses a critical mass, black hole growth may accelerate as the central gravitational potential deepens with strong ensuing star formation.

The collapse of massive stars is one of the most-studied paths to black hole formation. In this chapter, we review black hole formation during the collapse of massive stars in the broader context of single and binary stellar evolution and the theory of supernova explosions. We provide a concise overview of the evolutionary channels that may lead to black hole formation -- the classical route of iron core collapse, collapse due to pair instability in very massive stars, and the hypothetical scenario of supermassive star collapse. We then review the current understanding of the parameter space for black hole formation and black hole birth properties that has emerged from theoretical and computational modelling of supernova explosions and transient observations. Finally, we discuss what the intricate interplay between stellar evolution, stellar explosions, and binary interactions implies for the formation of stellar-mass black holes.

Stellar-mass black holes (BHs), with masses comparable to stars, are a major constituent of our Milky Way galaxy. This chapter describes the landscape of challenging, and long-sought efforts to identify these objects in the Galaxy. The first stellar-mass BHs were identified as persistent, but highly variable cosmic X-ray sources. Later, transient BH candidates were detected, and now far outnumber the persistent sources. Decades of effort have also yielded candidate BHs via gravitational microlensing and their orbital effect on binary companions. Populations of BH systems have begun to emerge from these detection strategies, offering insight into the astrophysical context in which BHs exist and driving questions about the formation, assembly, and ongoing evolution of these enigmatic objects.

We identify the first quiescent galaxies in TNG300, the largest volume of the IllustrisTNG cosmological simulation suite, and explore their quenching processes and time evolution to z=0. We find that the first quiescent galaxies with stellar masses M_* > 3 x 10^{10} M_sun and specific star formation rates sSFR < 10^{-11} yr^{-1} emerge at z~4.2 in TNG300. Suppression of star formation in these galaxies begins with a thermal mode of AGN feedback at z~6, and a kinetic feedback mode acts in each galaxy by z~4.7 to complete the quenching process, which occurs on a time-scale of ~0.35 Gyr. Surprisingly, we find that the majority of these galaxies are not the main progenitors of their z=0 descendants; instead, four of the five galaxies fall into more massive galaxies in subsequent mergers at a range of redshifts 2.5 < z < 0.2. By z=0, these descendants are the centres of galaxy clusters with average stellar masses of 8 x 10^{11} M_sun. We make predictions for the first quenched galaxies to be located by the James Webb Space Telescope (JWST).

Analyzing the observations obtained by the LIGO and the Virgo Collaborations, a new era has begun in binary black hole (BBH) merger processes and black hole physics studies. The fact that very massive stars that will become black holes at the end of their evolution are in binary or multiple states adds particular importance to BBH studies. In this study, using the SEOBNRv4$\_opt$ gravitational waveform model developed for compact binary systems, many ($\sim 10^6$) models were produced under different initial conditions, and the pre- and post-merge parameters were compared. In the models, it is assumed that the initial total mass (M$_{\rm{tot}}$) of the binary systems varies between 12-130 $\rm{M}_\odot$ with step interval 1$\rm{M}_\odot$, the mass ratios ($q = \rm{m}_{1i}/\rm{m}_{2i}$) vary between 1 and 2 with step interval 0.004, and the initial spin ($\abs{\rchi_{1i}} = \abs{\rchi_{2i}}$) value varies between $-0.83$ and $+0.83$ with step interval 0.017. Final spin ($\rchi_{f}$), fractional mass loss (M$_{FL}$), and the maximum gravitational wave amplitude (h$_{\rm{max}}$) obtained during the merger were compared with appropriate tables and figures obtained from the results of the relativistic numeric model obtained according to the initial parameters. Our results show that M$_{\rm{FL}}$ in generated BBH coalescences varied about 2.7 to 9.2\%, and $\rchi_{\rm{f}}$ between 0.29 and 0.91. In most of the BBHs we have modeled, we found that M$_{\rm{FL}}$ varies inversely with $q$. However, it has been found that M$_{\rm{FL}}$ values are not always inversely varied to the $q$ parameter in systems of opposite initial spin, where the large mass black hole component is positively oriented. Accordingly, it is understood that the values of M$_{\rm{FL}}$ decrease to a certain point of $q$ and then increase according to the increasing direction of $q$.

We explore properties of stellar kinematics and ionized gas in a sample of 1106 local [U]LIRGs from the AKARI telescope. We combine data from $Wide-field\ Infrared\ Survey\ Explorer$ (WISE) and Sloan Digital Sky Survey (SDSS) Data Release 13 (DR13) to fit the spectral energy distribution (SED) of each source to constrain the contribution of AGN to the total IR luminosity and estimate physical parameters such as stellar mass and star-formation rate (SFR). We split our sample into AGNs and weak/non-AGNs. We find that our sample is considerably above the main sequence. The highest SFRs and stellar masses are associated with ULIRGs. We also fit the H$\beta$ and H$\alpha$ regions to characterize the outflows. We find that the incidence of ionized gas outflows in AGN [U]LIRGs ($\sim$ 72\%) is much higher than that in weak/non-AGN ones ($\sim$ 39\%). The AGN ULIRGs have extreme outflow velocities (up to $\sim$ 2300 km s$^{-1}$) and high mass outflow rates (up to $\sim$ 60 \solarm~yr$^{-1}$). Our results suggest that starbursts are insufficient to produce such powerful outflows. We explore the correlations of SFR and specific SFR (sSFR) with ionized gas outflows. We find that AGN hosts with the highest SFRs exhibit a negative correlation between outflow velocity and sSFR. Therefore, in AGNs containing large amounts of gas, the negative feedback scenario might be suggested.

The hydrogen Balmer decrement (e.g., $\rm H\alpha/H\beta$) is widely adopted as an indicator of the internal reddening of active galactic nuclei (AGNs). This is challenged by some low-luminosity AGNs (LLAGNs) and changing-look AGNs (CLAGNs), which have steep Balmer decrement but without strong evidence for absorption. We compile a sample of normal AGNs and CLAGNs with a wider distribution of bolometric Eddington ratio ($\lambda_{\rm Edd}=L_{\rm bol}/L_{\rm Edd}$) and find a strong negative correlation between $\rm H\alpha/H\beta$ and $\lambda_{\rm Edd}$, which suggests that the Balmer decrement is also accretion-rate dependent. We further explore the Balmer decrement based on the photoionization model using the Cloudy code by considering spectral energy distribution (SED) from the accretion disk with different accretion rates (e.g., disk/corona and truncated disk at high and low Eddington ratios, respectively). Both the standard disk and truncated disk predict a negative correlation of $\rm H\alpha/H\beta-\lambda_{\rm Edd}$, where the relation is steeper in the case of the truncated disk. The negative correlations are also explored in two single CLAGNs. The measured negative correlation of $\rm H\alpha/H\beta$ -- $\lambda_{\rm Edd}$ is mainly caused by the lower responsivity $({\rm dlog}L_{\rm line}/{\rm dlog}L_{\rm cont})$ in $\rm H\alpha$ relative to that in $\rm H\beta$, due to the larger optical depth in the former. We propose that the steep Balmer decrements in low-Eddington-ratio AGNs (e.g., some Seyferts 1.5-1.9 and CLAGNs) are not simply caused by absorption but mainly caused by the relatively low flux of ionizing photons.

The feedback from accretion of central supermassive black holes (SMBHs) is a hot topic in the co-evolution of the SMBHs and their host galaxies. By tracing the large scale outflow by the line profile and bulk velocity shift of $[ \rm O~{\scriptsize III}]~ \lambda 5007$, the evolutionary role of outflow is studied here on a large sample of 221 type 2 quasars (QSO2s) extracted from Reyes et al. By following our previous study on local Seyfert 2 galaxies, the current spectral analysis on the SDSS spectroscopic database enables us to arrive at following results: (1) by using the Lick indices, we confirm that QSO2s are on average associated with younger stellar populations than Seyfert galaxies; (2) QSO2s with a stronger outflow are tend to be associated with a younger stellar population, which implies a coevolution between the feedback from SMBH and the host in QSO2s; (3) although an occupation at the high $L_{\rm bol}/L_{\rm Edd}$ end, the QSO2s follow the $L_{\rm bol}/L_{\rm Edd}$-$D_{n}(4000)$ sequence established from local, less-luminous Seyfert galaxies, which suggests a decrease of accretion activity of SMBH and feedback as the circumnuclear stellar population continuously ages.

In the context of a project aimed at characterizing the properties of star clusters in the Galactic bulge, here we present the determination of the internal kinematics and structure of the massive globular cluster NGC 6569. The kinematics has been studied by means of an unprecedented spectroscopic dataset acquired in the context of the ESO-VLT Multi-Instrument Kinematic Survey (MIKiS) of Galactic globular clusters, combining the observations from four different spectrographs. We measured the line-of-sight velocity of a sample of almost 1300 stars distributed between ~0.8" and 770" from the cluster center. From a sub-sample of high-quality measures, we determined the velocity dispersion profile of the system over its entire radial extension (from ~ 5" to ~ 200" from the center), finding the characteristic behavior usually observed in globular clusters, with a constant inner plateau and a declining trend at larger radii. The projected density profile of the cluster has been obtained from resolved star counts, by combining high-resolution photometric data in the center, and the Gaia EDR3 catalog radially extended out to ~20' for a proper sampling of the Galactic field background. The two profiles are properly reproduced by the same King model, from which we estimated updated values of the central velocity dispersion, main structural parameters (such as the King concentration, the core, half-mass, and tidal radii), total mass, and relaxation times. Our analysis also reveals a hint of ordered rotation in an intermediate region of the cluster (40"<r<90", corresponding to $ 2 r_c<r<4.5 r_c$), but additional data are required to properly assess this possibility.

We argue that observations of the reionization history such as the luminosity function of the Lyman-$\alpha$ emitters can be used as a probe of primordial density fluctuations, particularly on small scales. Although the primordial curvature perturbations are well constrained from measurements of cosmic microwave background (CMB) anisotropies and large-scale structure, these observational data probe the curvature perturbations only on large scales, and hence its information on smaller scales will give us further insight on primordial fluctuations. Since the formation of early galaxies is sensitive to the amplitude of small-scale perturbations, and then, in turn, gives an impact on the reionization history, one can probe the primordial power spectrum on small scales through observations of reionization. In this work, we focus on the running spectral indices of the primordial power spectrum to characterize the small-scale perturbations, and investigate their constraints from observations of the luminosity function of the Lyman-$\alpha$ emitters. We show that the reionization, in combination with large-scale observations such as CMB, would be a useful tool to investigate primordial density fluctuations.

To self-consistently model galactic properties, reionization of the intergalactic medium, and the associated 21-cm signal, we have developed the algorithm polar by integrating the one-dimensional radiative transfer code grizzly with the semi-analytical galaxy formation code L-Galaxies 2020. Our proof-of-concept results are consistent with observations of the star formation rate history, UV luminosity function and the CMB Thomson scattering optical depth. We then investigate how different galaxy formation models affect UV luminosity functions and 21-cm power spectra, and find that while the former are most sensitive to the parameters describing the merger of halos, the latter have a stronger dependence on the supernovae feedback parameters, and both are affected by the escape fraction model.

In this paper we aim to constrain for the first time the dust emission in the mid-to-far infrared domain, in the LMC, with the use of the Spitzer IRS and MIPS SED data, combined with Herschel data. We also consider UV extinction predictions derived from modeling. We selected 10 regions observed as part of the SAGE-Spec program, to probe dust properties in various environments (diffuse, molecular and ionized regions). All data were smoothed to the 40arcsec angular resolution. The SEDs were modeled with DustEM models, using the standard Mathis RF, as well as three additional RFs, with stellar clusters ages ranging from 4 Myr to 600 Myr. Standard dust models used to reproduce the Galactic diffuse medium are clearly not able to reproduce the dust emission in the MIR wavelength domain. This analysis evidences the need of adjusting parameters describing the dust size distribution and shows a clear distinct behavior according to the type of environments. In addition, whereas the small grain emission always seems to be negligible at long wavelengths in our Galaxy, the contribution of this small dust component could be more important than expected, in the submm-mm range, in the LMC averaged SED. Properties of the small dust component of the LMC are clearly different from those of our Galaxy. Its abundance, significantly enhanced, could be the result of large grains shattering due to strong shocks or turbulence. In addition, this grain component in the LMC systematically shows smaller grain size in the ionized regions compared to the diffuse medium. Predictions of extinction curves show significantly distinct behaviors depending on the dust models but also from one region to another. Comparison of model predictions with the LMC mean extinction curve shows that no model gives satisfactory agreement using the Mathis radiation field while using a harder radiation field tends to improve the agreement

Most transit microlensing events due to very low-mass lens objects suffer from extreme finite-source effects. While modeling their light curves, there is a known continuous degeneracy between their relevant lensing parameters, i.e., the source angular radius normalized to the angular Einstein radius $\rho_{\star}$, the Einstein crossing time $t_{\rm E}$, the lens impact parameter $u_{0}$, the blending parameter, and the stellar apparent magnitude. In this work, I numerically study the origin of this degeneracy. I find that these light curves have 5 observational parameters (i.e., the baseline magnitude, the maximum deviation in the magnification factor, the Full Width at Half Maximum $\rm{FWHM}=2 t_{\rm{HM}}$, the deviation from top-hat model, the time of the maximum time-derivative of microlensing light curves $T_{\rm{max}}=t_{\rm E}\sqrt{\rho_{\star}^{2}-u_{0}^{2}}$). For extreme finite-source microlensing events due to uniform source stars we get $t_{\rm{HM}}\simeq T_{\rm{max}}$, and the deviation from the top-hat model tends to zero which both cause the known continuous degeneracy. When either $\rho_{\star}\lesssim10$ or the limb-darkening effect is considerable $t_{\rm{HM}}$, and $T_{\rm{max}}$ are two independent observational parameters. I use a numerical approach, i.e., Random Forests containing $100$-$120$ Decision Trees, to study how these observational parameters are efficient in yielding the lensing parameters. These machine learning models find the mentioned 5 lensing parameters for finite-source microlensing events from uniform, and limb-darkened source stars with the average $R^{2}$-scores of $0.87$, and $0.84$, respectively. $R^{2}$-score for evaluating the lens impact parameter gets worse on adding limb darkening, and for extracting the limb-darkening coefficient itself this score falls as low as $0.67$.

If a hypothetical elementary particle called an axion exists, to solve the strong CP problem, a 57Fe nucleus in the solar core could emit a 14.4-keV monochromatic axion through the M1 transition. If such axions are once more transformed into photons by a 57Fe absorber, a transition edge sensor (TES) X-ray microcalorimeter should be able to detect them efficiently. We have designed and fabricated a dedicated 64-pixel TES array with iron absorbers for the solar axion search. In order to decrease the effect of iron magnetization on spectroscopic performance, the iron absorber is placed next to the TES while maintaining a certain distance. A gold thermal transfer strap connects them. We have accomplished the electroplating of gold straps with high thermal conductivity. The residual resistivity ratio (RRR) was over 23, more than eight times higher than a previous evaporated strap. In addition, we successfully electroplated pure-iron films of more than a few micrometers in thickness for absorbers and a fabricated 64-pixel TES calorimeter structure.

Low mass star formation regions are unlikely to fully populate their initial mass functions, leading to a deficit of massive stars. In binary stellar populations, the full range of binary separations and mass ratios will also be underpopulated. To explore the effects of stochastic sampling in the integrated light of stellar clusters, we calculate models at a broad range of cluster masses, from 10^2 to 10^7 M_sun, using a binary stellar population synthesis code. For clusters with stellar masses less than 10^5 M_sun, observable quantities show substantial scatter and their mean properties reflect the expected deficit of massive stars. In common with previous work, we find that purely stochastic sampling of the initial mass function appears to underestimate the mass of the most massive star in known clusters. However, even with this constraint, the majority of clusters likely inject sufficient kinetic energy to clear their birth clusters of gas. For quantities which directly measure the impact of the most massive stars, such as N_{ion}, xi_{ion} and beta_{UV}, uncertainties due to stochastic sampling dominate over those from the IMF shape or distribution of binary parameters, while stochastic sampling has a negligible effect on the stellar continuum luminosity density.

The Slow Solar Wind Connection Solar Orbiter Observing Plan (Slow Wind SOOP) was developed to utilise the extensive suite of remote sensing and in situ instruments on board the ESA/NASA Solar Orbiter mission to answer significant outstanding questions regarding the origin and formation of the slow solar wind. The Slow Wind SOOP was designed to link remote sensing and in situ measurements of slow wind originating at open-closed field boundaries. The SOOP ran just prior to Solar Orbiter's first close perihelion passage during two remote sensing windows (RSW1 and RSW2) between 2022 March 3-6 and 2022 March 17-22, while Solar Orbiter was at a heliocentric distance of 0.55-0.51 and 0.38-0.34 au from the Sun, respectively. Coordinated observation campaigns were also conducted by Hinode and IRIS. The magnetic connectivity tool was used, along with low latency in situ data, and full-disk remote sensing observations, to guide the target pointing of Solar Orbiter. Solar Orbiter targeted an active region complex during RSW1, the boundary of a coronal hole, and the periphery of a decayed active region during RSW2. Post-observation analysis using the magnetic connectivity tool along with in situ measurements from MAG and SWA/PAS, show that slow solar wind, with velocities between 210 and 600 km/s, arrived at the spacecraft originating from two out of the three of the target regions. The Slow Wind SOOP, despite presenting many challenges, was very successful, providing a blueprint for planning future observation campaigns that rely on the magnetic connectivity of Solar Orbiter.

Strange stars ought to exist in the universe according to the strange quark matter hypothesis, which states that matter made of roughly equal numbers of up, down, and strange quarks could be the true ground state of baryonic matter rather than ordinary atomic nuclei. Theoretical models of strange quark matter, such as the standard MIT bag model, the density-dependent quark mass model, or the quasi-particle model, however, appear to be unable to reproduce some of the properties (masses, radii and tidal deformabilities) of recently observed compact stars. This is different if alternative gravity theory (e.g., non-Newtonian gravity) or dark matter (e.g., mirror dark matter) are considered, which resolve these issues. The possible existence of strange stars could thus provide a clue to new physics, as discussed in this review.

The study of the elusive hot component ($T \gtrsim 10^7$ K) of the Milky Way circumgalactic medium (CGM) is a novel topic to understand Galactic formation and evolution. In this work, we use the stacking technique through 46 lines of sight with Chandra ACIS-S HETG totaling over 10Ms of exposure time and 9 lines of sight with ACIS-S LETG observations totaling over 1Ms of exposure time, to study in absorption the presence of highly ionized metals arising from the super-virial temperature phase of the CGM. Focusing in the spectral range $4 - 8$ $\r{A}$, we were able to confirm the presence of this hot phase with high significance. We detected transitions of Si XIV K$\alpha$ (with total significance of 6.0$\sigma$) and, for the first time, SXVI K (total significance 4.8$\sigma$) in the rest frame of our own Galaxy. For S XVI K$\alpha$ we found a column density of $1.50^{+0.44}_{-0.38} \times 10^{16} \mathrm{cm}^{-2}$. For Si XIV K$\alpha$ we measured a column density of $0.87\pm{0.16} \times 10^{16} \mathrm{cm}^{-2}$. The lines of sight used in this work are spread across the sky, probing widely separated regions of the CGM. Therefore, our results indicate that this newly discovered hot medium extends throughout the halo, and is not related only to the Galactic Bubbles. The hot gas location, distribution, and covering factor, however, remain unknown. This component might contribute significantly to the missing baryons and metals in the Milky Way.

In this paper we present an analysis of near-infrared spectropolarimetric and velocimetric data of the young M dwarf AU Mic, collected with SPIRou at the Canada-France-Hawaii telescope from 2019 to 2022, mostly within the SPIRou Legacy Survey. With these data, we study the large- and small-scale magnetic field of AU Mic, detected through the unpolarized and circularly-polarized Zeeman signatures of spectral lines. We find that both are modulated with the stellar rotation period (4.86 d), and evolve on a timescale of months under differential rotation and intrinsic variability. The small-scale field, estimated from the broadening of spectral lines, reaches $2.61\pm0.05$ kG. The large-scale field, inferred with Zeeman-Doppler imaging from Least-Squares Deconvolved profiles of circularly-polarized and unpolarized spectral lines, is mostly poloidal and axisymmetric, with an average intensity of $550\pm30$ G. We also find that surface differential rotation, as derived from the large-scale field, is $\simeq$30% weaker than that of the Sun. We detect the radial velocity (RV) signatures of transiting planets b and c, although dwarfed by activity, and put an upper limit on that of candidate planet d, putatively causing the transit-timing variations of b and c. We also report the detection of the RV signature of a new candidate planet (e) orbiting further out with a period of $33.39\pm0.10$ d, i.e., near the 4:1 resonance with b. The RV signature of e is detected at 6.5$\sigma$ while those of b and c show up at $\simeq$4$\sigma$, yielding masses of $10.2^{+3.9}_{-2.7}$ and $14.2^{+4.8}_{-3.5}$ Earth masses for b and c, and a minimum mass of $35.2^{+6.7}_{-5.4}$ Earth masses for e.

When a $\gamma$-ray interacts in a semiconductor detector, the resulting electron-hole charge clouds drift towards their respective electrodes for signal collection. These charge clouds will expand over time due to both thermal diffusion and mutual electrostatic repulsion. Solutions to the resulting charge profiles are well understood for the limiting cases accounting for only diffusion and only repulsion, but the general solution including both effects can only be solved numerically. Previous attempts to model these effects have taken into account the broadening of the charge profile due to both effects, but have simplified the shape of the profile by assuming Gaussian distributions. However, the detailed charge profile can have important impacts on charge sharing in multi-electrode strip detectors. In this work, we derive an analytical approximation to the general solution, including both diffusion and repulsion, that closely replicates both the width and the detailed shape of the charge profiles. This analytical solution simplifies the modeling of charge clouds in semiconductor strip detectors.

Ho\v{r}ava-Lifshitz gravity has been proposed as a ghost-free quantum gravity model candidate with an anisotropic UV-scaling between space and time. We present here a cosmological background analysis of two different formulations of the theory, with particular focus on the running of the parameter $\lambda$. Using a large dataset consisting of Cosmic Microwave Background data from {\it Planck}, Pantheon+ supernovae catalogue, SH0ES Cepheid variable stars, Baryon acoustic oscillations (BAO), Cosmic Chronometers, and gamma-ray bursts (GRB), we arrive at new bounds on the cosmological parameters, in particular $\lambda$, which describes deviation from classical general relativity. For the detailed balance scenario we arrive at the bound $\lambda=1.02726\pm0.00012$, and for beyond detailed balance the limit reads $\lambda=0.9949^{+0.0045}_{-0.0046}$. We also study the influence of different data sets and priors, and we find that removing low-redshift data generally moves $\lambda$ closer towards UV values, whilst simultaneously widening the error bars. In the detailed balance scenario, this effect is more noticeable, and $\lambda$ takes on values that are significantly below unity, which corresponds to the infrared limit of the theory.

Stars with an initial mass below ~ 8 Msun evolve through the asymptotic giant branch (AGB) phase, during which they develop strong stellar winds. Recent observations have revealed significant morphological complexities in their outflows, most likely caused by a companion. We study the impact of the radiation force on such companion-perturbed AGB outflows. We present the implementation of a ray tracer for radiative transfer in smoothed particle hydrodynamics (SPH) and compared four different descriptions of radiative transfer: the free-wind, the geometrical, the Lucy, and the attenuation approximation. For both low and high mass-loss rates, the velocity profile of the outflow is modified when going from the free-wind to the geometrical approximation, also resulting in a different morphology. In the case of a low mass-loss rate, the effect of the Lucy and attenuation approximation is negligible due to the low densities but morphological differences appear in the high mass-loss rate regime. By comparing the radiative equilibrium temperature and radiation force to full 3D radiative transfer, we show that the Lucy approximation works best. Although, close to the companion, artificial heating occurs and it fails to simulate the shadow cast by the companion. The attenuation approximation produces a lower equilibrium temperature and weaker radiation force, but it produces the shadow cast by the companion. From the predictions of the 3D radiative transfer, we also conclude that a radially directed radiation force is a reasonable assumption. The radiation force thus plays a critical role in dust-driven AGB winds, impacting the velocity profile and morphological structures. For low mass-loss rates, the geometrical approximation suffices, while high mass-loss rates require a more rigorous method, where the Lucy approximation provides the most accurate results although not accounting for all effects.

The SunPy Project is a community of scientists and software developers creating an ecosystem of Python packages for solar physics. The project includes the sunpy core package as well as a set of affiliated packages. The sunpy core package provides general purpose tools to access data from different providers, read image and time series data, and transform between commonly used coordinate systems. Affiliated packages perform more specialized tasks that do not fall within the more general scope of the sunpy core package. In this article, we give a high-level overview of the SunPy Project, how it is broader than the sunpy core package, and how the project curates and fosters the affiliated package system. We demonstrate how components of the SunPy ecosystem, including sunpy and several affiliated packages, work together to enable multi-instrument data analysis workflows. We also describe members of the SunPy Project and how the project interacts with the wider solar physics and scientific Python communities. Finally, we discuss the future direction and priorities of the SunPy Project.

Combining new HI data from a synergetic survey of ASKAP WALLABY and FAST with the ALFALFA data, we study the effect of ram-pressure and tidal interactions in the NGC 4636 group. We develop two parameters to quantify and disentangle these two effects on gas stripping in HI-bearing galaxies: the strength of external forces at the optical-disk edge, and the outside-in extents of HI-disk stripping. We find that gas stripping is widespread in this group, affecting 80% of HI-detected non-merging galaxies, and that 34% are experiencing both types of stripping. Among the galaxies experiencing both effects, the strengths (and extents) of ram-pressure and tidal stripping are independent of each other. Both strengths are correlated with HI-disk shrinkage. The tidal strength is related to a rather uniform reddening of low-mass galaxies ($M_*<10^9\,\text{M}_\odot$) when tidal stripping is the dominating effect. In contrast, ram pressure is not clearly linked to the color-changing patterns of galaxies in the group. Combining these two stripping extents, we estimate the total stripping extent, and put forward an empirical model that can describe the decrease of HI richness as galaxies fall toward the group center. The stripping timescale we derived decreases with distance to the center, from $\mathord{\sim}1\,\text{Gyr}$ around $R_{200}$ to $\mathord{\lesssim}10\,\text{Myr}$ near the center. Gas-depletion happens $\mathord{\sim}3\,\text{Gyr}$ since crossing $2R_{200}$ for HI-rich galaxies, but much quicker for HI-poor ones. Our results quantify in a physically motivated way the details and processes of environmental-effects-driven galaxy evolution, and might assist in analyzing hydrodynamic simulations in an observational way.

With the help of our previously built MCMC-based parameter estimation package \texttt{CosmoReionMC}, we investigate in detail the potential of 21~cm global signal, when combined with CMB and observations related to the QSO absorption spectra, to constraint the mass of Warm Dark Matter (WDM) particle. For the first time, we simultaneously vary all the free parameters (mass of WDM particle, cosmological parameters, and astrophysical parameters) to address the long-overlooked issue of the possible degeneracies between the Dark Matter particle mass $m_X$ and cosmological/astrophysical parameters. From the existing CMB and QSO absorption spectra data, we can rule out $m_X < 2.8$~keV at 95\% confidence level. Including the mock 21~cm global signal data expected in the future, the forecasted constraint is found to be much tighter $m_X > 7.7$~keV, assuming that the true dark matter model is the usual cold dark matter. In case the mock 21~cm signal is constructed for dark matter particles having $m_X = 7$~keV, our forecasts indicate that $\left(m_X / \text{keV}\right)^{-1}$ is in the range $[0.1, 0.2]$ ($95\%$ confidence level). This implies that the future 21~cm data should allow detection of the WDM particle mass if $m_X \sim 7$~keV

We present results from the first 22 GHz space very-long-baseline interferometric (VLBI) imaging observations of M87 by RadioAstron. As a part of the Nearby AGN Key Science Program, the source was observed in Feb 2014 at 22 GHz with 21 ground stations, reaching projected $(u,v)$-spacings up to $\sim11\,$G$\lambda$. The imaging experiment was complemented by snapshot RadioAstron data of M87 obtained during 2013--2016 from the AGN Survey Key Science Program. Their longest baselines extend up to $\sim25\,$G$\lambda$. For all these measurements, fringes are detected only up to $\sim$2.8 Earth Diameter or $\sim$3 G$\lambda$ baseline lengths, resulting in a new image with angular resolution of $\sim150\,\mu$as or $\sim20$ Schwarzschild radii spatial resolution. The new image not only shows edge-brightened jet and counterjet structures down to submilliarcsecond scales but also clearly resolves the VLBI core region. While the overall size of the core is comparable to those reported in the literature, the ground-space fringe detection and slightly super-resolved RadioAstron image suggest the presence of substructures in the nucleus, whose minimum brightness temperature exceeds $T_{\rm B, min}\sim10^{12}\,$K. It is challenging to explain the origin of this record-high $T_{\rm B, min}$ value for M87 by pure Doppler boosting effect with a simple conical jet geometry and known jet speed. Therefore, this can be evidence for more extreme Doppler boosting due to a blazar-like small jet viewing angle or highly efficient particle acceleration processes occurring already at the base of the outflow.

The light curve of the cataclismic variable ASASSN-18aan is studied using recent observations of the MC589 Observatory, giving an orbital Period and Epoch fully consistent with the data obtained after the discovery flare in 2018. Archival data from ASASSN, ZTF and Gaia were used to check if its flares have a quasi-periodic behaviour. A recurrency time scale of about 11 months is found, confirming a previous tentative result using the historic plate archive of the Asiago Observatory. The next outbursts are expected by April 2023 and March 2024.

We present a new upper limit on the cosmic molecular gas density at $z=2.4-3.4$ obtained using the first year of observations from the CO Mapping Array Project (COMAP). COMAP data cubes are stacked on the 3D positions of 282 quasars selected from the Extended Baryon Oscillation Spectroscopic Survey (eBOSS) catalog, yielding a 95% upper limit for flux from CO(1-0) line emission of 0.210 Jy km/s. Depending on the assumptions made, this value can be interpreted as either an average CO line luminosity $L'_\mathrm{CO}$ of eBOSS quasars of $\leq 7.30\times10^{10}$ K km pc$^2$ s$^{-1}$, or an average molecular gas density $\rho_\mathrm{H_2}$ in regions of the universe containing a quasar of $\leq 2.02\times10^8$ M$_\odot$ cMpc$^{-3}$. The $L'_\mathrm{CO}$ upper limit falls among CO line luminosities obtained from individually-targeted quasars in the COMAP redshift range, and the $\rho_\mathrm{H_2}$ value is comparable to upper limits obtained from other Line Intensity Mapping (LIM) surveys and their joint analyses. Further, we forecast the values obtainable with the COMAP/eBOSS stack after the full 5-year COMAP Pathfinder survey. We predict that a detection is probable with this method, depending on the CO properties of the quasar sample. Based on these achieved sensitivities, we believe that this technique of stacking LIM data on the positions of traditional galaxy or quasar catalogs is extremely promising, both as a technique for investigating large galaxy catalogs efficiently at high redshift and as a technique for bolstering the sensitivity of LIM experiments, even with a fraction of their total expected survey data.

We investigate compact objects in modified teleparallel gravity with realistic equations of state. We propose a modification on Teleparallel Equivalent of General Relativity, then an appropriate tetrad is applied on the field equations. A specific set of relations showing a equivalency between our gravitational model and the New General Relativity is found. The conservation equation implies that our Tolman-Oppenheimer-Volkoff equations are presented with an effective pressure and energy density, where a free parameter \b{eta}3 is used to construct them. Numerical analysis using realistic equations of state is made, the behavior of mass, radius and the relation mass-radius as functions of \b{eta}3 is also investigated.

This work investigates dark matter (DM) effects in compact objects in modified teleparallel gravity (MTG) in which a modification of Teleparallel Equivalent to General Relativity is used. We applied a tetrad to the modified field equations where a set of relations is found. The conservation equation allows us to rewrite our Tolman-Oppenheimer-Volkoff equations with an effective gravitational coupling constant. As input to these new equations, we use a relativistic mean-field (RMF) model with dark matter content included, obtained from a Lagrangian density with both, hadronic and dark particle degrees of freedom, as well as the Higgs boson, used as a mediator in both sectors of the theory. Through numerical calculations, we analyze the mass-radius diagrams obtained from different parametrizations of the RMF-DM model, generated by assuming different values of the dark particle Fermi momentum and running the free parameter coming from the MTG. Our results show that it is possible for the system simultaneously support more DM content, and be compatible with recent astrophysical data provided by LIGO and Virgo Collaboration, as well as by NASA's Neutron star Interior Composition Explorer (NICER).

In the post-LIGO era, there has been a lot of focus on primordial black holes (PBHs) heavier than $\sim 10^{15}$g as potential dark matter (DM) candidates. We point out that the branch of the PBH family that disappeared - PBHs lighter than $\sim 10^9$g that ostensibly Hawking evaporated away in the early Universe - also constitute an interesting frontier for DM physics. Hawking evaporation itself serves as a portal through which such PBHs can illuminate new physics, for example by emitting dark sector particles. Taking a simple DM scalar singlet model as a template, we compute the abundance and mass of PBHs that could have provided, by Hawking evaporation, the correct DM relic density. We consider two classes of such PBHs: those originating from curvature perturbations generated by inflation, and those originating from false vacuum collapse during a first-order phase transition. For PBHs of both origins we compute the gravitational wave (GW) signals emanating from their formation stage: from second-order effects in the case of curvature perturbations, and from sound waves in the case of phase transitions. The GW signals have peak frequencies in the MHz-GHz range typical of such light PBHs. We compute the strength of such GWs compatible with the observed DM relic density, and find that the GW signal morphology can in principle allow one to distinguish between the two PBH formation histories.

The higher order generalisation of the clockwork mechanism to gravitational interactions provides a means to generate an exponentially suppressed coupling to matter from a fundamental theory of multiple interacting gravitons, without introducing large hierarchies in the underlying potential and without the need for a dilaton, suggesting a possible application to the hierarchy problem. We work in the framework of ghost free multi-gravity with "nearest-neighbour" interactions, and present a formalism by which one is able to construct potentials such that the theory will always exhibit this clockwork effect. We also consider cosmological solutions to the general theory, where all metrics are of FRW form, with site-dependent scale factors/lapses. We demonstrate the existence of multiple deSitter vacua where all metrics share the same Hubble parameter, and we solve the modified Einstein equations numerically for an example clockwork model constructed using our formalism, finding that the evolution of the metric that matter couples to is essentially equivalent to that of general relativity at the modified Planck scale. It is important to stress that while we focus on the application to clockwork theories, our work is entirely general and facilitates finding cosmological solutions to any ghost free multi-gravity theory with "nearest-neighbour" interactions. Moreover, we clarify previous work on the continuum limit of the theory, which is generically a scalar-tensor braneworld, using the Randall-Sundrum model as a special case and showing how the discrete-clockwork cosmological results map to the continuum results in the appropriate limit.

We consider a Bardeen-Cooper-Schrieffer (BCS)-like model in the inflationary background. We show that with an axial chemical potential, the attractive quartic fermion self-interaction can lead to a BCS-like condensation. In the de Sitter (dS) limit of inflation, we perform the first computation of the non-perturbative effective potential that includes the full spacetime curvature effects in the presence of the chemical potential. The corresponding BCS phase transition is always first-order, when the varying Hubble is interpreted as an effective Gibbons-Hawking temperature of dS spacetime. In the condensate phase, the theory can be understood from UV and IR sides as fermionic and bosonic, respectively. This leads to distinctive signatures in the primordial non-Gaussianity of curvature perturbations. Namely, the oscillatory cosmological collider signal is smoothly turned off at a finite momentum ratio, since different momentum ratios effectively probe different energy scales. In addition, such BCS phase transitions can also source stochastic gravitational waves, feasible for future experiments.

The emergence of macroscopic coherence in a many-body quantum system is a ubiquitous phenomenon across different physical systems and scales. This Chapter reviews key concepts characterizing such systems (correlation functions, condensation, quasi-condensation) and applies them to the study of emerging non-equilibrium features in the dynamical path towards such a highly-coherent state: particular emphasis is placed on emerging universal features in the dynamics of conservative and open quantum systems, their equilibrium or non-equilibrium nature, and the extent that these can be observed in current experiments with quantum gases. Characteristic examples include symmetry-breaking in the Kibble-Zurek mechanism, coarsening and phase-ordering kinetics, and universal spatiotemporal scalings around non-thermal fixed points and in the context of the Kardar- Parisi-Zhang equation; the Chapter concludes with a brief review of the potential relevance of some of these concepts in modelling the large-scale distribution of dark matter in the universe.

It was recently claimed that black holes can explain the accelerated expansion of the universe. Here we point out that this claim is based on a confusion about the principle of least action, undermining the link between black holes and dark energy.