Papers for Friday, Mar 04 2022

Most globular clusters (GCs) show evidence for multiple stellar populations, suggesting the occurrence of several distinct star-formation episodes. The large fraction of second population (2P) stars observed requires a very large 2P gaseous mass to have accumulated in the cluster core to form these stars. Hence the first population of stars (1P) in the cluster core has had to become embedded in 2P gas, just prior to the formation of later populations. Here we explore the evolution of binaries in ambient 2P gaseous media of multiple-population GCs. We mostly focus on black hole binaries and follow their evolution as they evolve from wide binaries towards short periods through interaction with ambient gas, followed by gravitational-wave (GW) dominated inspiral and merger. We show this novel GW-merger channel could provide a major contribution to the production of GW-sources. We consider various assumptions and initial conditions and calculate the resulting gas-mediated change in the population of binaries and the expected merger rates due to gas-catalyzed GW-inspirals. For plausible conditions and assumptions, we find an expected GW merger rate observable by aLIGO of the order of up to a few tens of $\rm{Gpc^{-3} yr^{-1}}$, and an overall range for our various models of $0.08-25.51 \ \rm{Gpc^{-3} yr^{-1}}$. Finally, our results suggest that the conditions and binary properties in the early stage of GCs could be critically affected by gas-interactions and may require a major revision in the current modeling of the evolution of GCs.

We analyse how drag forces modify the orbits of objects moving through extended gaseous distributions. We consider how hydrodynamic (surface area) drag forces and dynamical friction (gravitational) drag forces drive the evolution of orbital eccentricity. While hydrodynamic drag forces cause eccentric orbits to become more circular, dynamical friction drag can cause orbits to become more eccentric. We develop a semi-analytic model that accurately predicts these changes by comparing the total work and torque applied to the orbit at periapse and apoapse. We use a toy model of a radial power-law density profile, $\rho \propto r^{-\gamma}$, to determine that there is a critical $\gamma = 3$ power index which separates the eccentricity evolution: orbits become more eccentric for $\gamma < 3$ and circularize for $\gamma > 3$. We apply these findings to the infall of a Jupiter-like planet into the envelope of its host star. The hydrostatic envelopes of stars are defined by steep density gradients near the limb and shallower gradients in the interior. Under the influence of gaseous dynamical friction, an infalling object's orbit will first decrease in eccentricity, then increase. The critical separation that delineates these regimes is predicted by the local density slope and is linearly dependent on polytropic index. More broadly, our findings indicate that binary systems may routinely emerge from common envelope phases with non-zero eccentricities that were excited by the dynamical friction forces that drove their orbital tightening.

The analysis of star cluster ages in tandem with the morphology of their HII regions can provide insight into the processes that clear a cluster's natal gas, as well as the accuracy of cluster ages and dust reddening derived from Spectral Energy Distribution (SED) fitting. We classify 3757 star clusters in 16 nearby galaxies according to their H$\alpha$ morphology (concentrated, partially exposed, no emission), using Hubble Space Telescope (HST) imaging from the Legacy ExtraGalactic Ultraviolet Survey (LEGUS). We find: 1) The mean SED ages of clusters with concentrated (1-2 Myr) and partially exposed HII region morphologies (2-3 Myr) indicate a relatively early onset of gas clearing and a short (1-2 Myr) clearing timescale. 2) The reddening of clusters can be overestimated due to the presence of red supergiants, which is a result of stochastic sampling of the IMF in low mass clusters. 3) The age-reddening degeneracy impacts the results of the SED fitting - out of 1408 clusters with M$_*$ $\geq$ 5000 M$_{\odot}$, we find that at least 46 (3%) have SED ages which appear significantly underestimated or overestimated based on H$\alpha$ and their environment, while the total percentage of poor age estimates is expected to be several times larger. 4) Lastly, we examine the dependence of the morphological classifications on spatial resolution. At HST resolution, our conclusions are robust to the distance range spanned by the sample (3-10 Mpc). However, analysis of ground-based H$\alpha$ images shows that compact and partially exposed morphologies frequently cannot be distinguished from each other.

We present a study of the variation timescales of the chromospheric activity indicator H$\alpha$ on a sample of 13 fully-convective, active mid-to-late M stars with masses between 0.1--0.3 solar masses. Our goal was to determine the dominant variability timescale and, by inference, a possible mechanism responsible for the variation. We gathered 10 or more high-resolution spectra each of 10 stars using the TRES spectrograph at times chosen to span all phases of stellar rotation, as determined from photometric data from the MEarth Observatories. All stars varied in their H$\alpha$ emission. For 9 of these stars, we found no correlation between H$\alpha$ and rotational phase, indicating that constant emission from fixed magnetic structures, such as starspots and plage, are unlikely to be the dominant source of H$\alpha$ emission variability. In contrast, one star, G 7-34, shows a clear relationship between H$\alpha$ and stellar rotational phase. Intriguingly, we found that this star is a member of the AB Doradus moving group and hence has the young age of 149 Myr. High-cadence spectroscopic observations of three additional stars revealed that they are variable on timescales ranging from 20--45 minutes, which we posit may be due to flaring behavior. For one star, GJ 1111, simultaneous TESS photometry and spectroscopic monitoring show an increase in H$\alpha$ emission with increased photometric brightness. We conclude that low-energy flares are able to produce variation in H$\alpha$ on the timescales we observe and thus may be the dominant source of H$\alpha$ variability on active fully-convective M dwarfs.

We report ALMA Band 3 observations of CO(6-5), CO(7-6) and [CI](2-1) in B14-65666 ("Big Three Dragons"), one of the brightest Lyman-Break Galaxies at $z>7$ in the rest-frame ultraviolet continuum, far-infrared continuum, and emission lines of [OIII] 88 $\mu$m and [CII] 158 $\mu$m. None of CO(6-5), CO(7-6) and [CI](2-1) are detected, whose $3\sigma$ upper limits on the luminosities are about 50 times fainter than the [CII] luminosity. Based on three methods of i) [CII] luminosity and a [CII]-to-H$_{2}$ conversion factor reported in local metal-poor dwarf galaxies, ii) a dust mass and metallicity-dependent dust-to-gas mass ratio, and iii) a dynamical mass estimate, we obtain the molecular gas mass ($M_{\rm mol}$) to be $(0.05-11)\times10^{10} M_{\rm \odot}$, which is consistent with its upper limit inferred from the non-detection of mid-$J$ CO and [CI](2-1). Albeit with large uncertainty in $M_{\rm mol}$, we estimate a high molecular gas-to-stellar mass ratio ($\mu_{\rm gas}$) of $0.65-140$ and a short gas depletion time ($t_{\rm dep}$) of $2.5-550$ Myr, which are broadly consistent with extrapolations of $\mu_{\rm gas}$ and $t_{\rm dep}$ as functions of redshift, specific-star formation rate, and stellar mass as reported in previous studies. The short $t_{\rm dep}$ partly reflects the starburst nature of the target, likely to be induced by a major-merger event. B14-65666 can be an ancestor of a passive galaxy at $z\gtrsim4$ if no gas is fueled from outside the galaxy. Overall, our observations highlight difficulty in obtaining robust conclusion from mid-$J$ CO alone particularly in the epoch of reionization, and motivate studies that use other tracers as an alternative gas mass tracer.

Currently time-domain astronomy can scan the entire sky on a daily basis, discovering thousands of interesting transients every night. Classifying the ever-increasing number of new transients is one of the main challenges for the astronomical community. One solution that addresses this issue is the robotically controlled Spectral Energy Distribution Machine (SEDM) which supports the Zwicky Transient Facility (ZTF). SEDM with its pipeline PYSEDM demonstrates that real-time robotic spectroscopic classification is feasible. In an effort to improve the quality of the current SEDM data, we present here two new modules, BYECR and CONTSEP. The first removes contamination from cosmic rays, and the second removes contamination from non-target light. These new modules are part of the automated PYSEDM pipeline and fully integrated with the whole process. Employing BYECR and CONTSEP modules together automatically extracts more spectra than the current PYSEDM pipeline. Using SNID classification results, the new modules show an improvement in the classification rate and accuracy of 2.8% and 1.7%, respectively, while the strength of the cross-correlation remains the same. Improvements to the SEDM astrometry would further boost the improvement of the CONTSEP module. This kind of robotic follow-up with a fully automated pipeline has the potential to provide the spectroscopic classifications for the transients discovered by ZTF and also by the Rubin Observatory's Legacy Survey of Space and Time.

Some of the radio galaxies show distortion in their jets, forming tailed or winged sources. X-shaped radio galaxies are a subclass of winged sources formation mechanism of which is still unclear. The focus of this work is to understand hydro-dynamical back-flows and their role in dynamics and non-thermal emission signatures (in presence of radiative losses and diffusive shock acceleration) during the initial phase of these galaxies. We have performed relativistic MHD simulations of an under-dense jet travelling in a tri-axial ambient using a hybrid Eulerian-Lagrangian framework to incorporate effects of micro-physical processes. We have demonstrated the dominant role played by pressure gradient in shaping XRGs in thermally dominated cases. We show that the prominence of the formed structure decreases as the jet deviates from the major axis of the ambient. The wing evolution is mainly governed by re-energized particles due to shocks that keep the structure active during the evolution time. The synthetic intensity maps of the radio galaxy show similarities with morphologies that are typically found in observed XRGs. This includes the cases with wider wings than the active lobes. The characteristic emission signatures in terms of its synchrotron spectra and implication of equipartition condition in age estimation are also discussed here. Additionally, we show that discrepancy of age can be attributed to mixing of different aged particle populations. Further, the effect of viewing angle on the difference of spectral index of the active lobes and the wings $(\Delta \alpha)$ shows a large variation and degenerate behaviour. We have demonstrated the role of diffusive shocks in the obtained variation and have concluded that the spread of $(\Delta \alpha)$ is not a dependable characteristic in determining the formation model of XRGs.

We present Hubble Space Telescope Wide Field Camera 3 photometric and grism observations of the candidate ultra-high-redshift (z>7) radio galaxy, GLEAM J0917-0012. This radio source was selected due to the curvature in its 70-230 MHz, low-frequency Murchison Widefield Array radio spectrum and its faintness in K-band. Follow-up spectroscopic observations of this source with the VLA and ALMA were inconclusive as to its redshift. Our F105W and F0986M imaging observations detect the host of GLEAM J0917-0012 and a companion galaxy, ~one arcsec away. The G102 grism observations reveal a single weak line in each of the spectra of the host and the companion. To help identify these lines we utilised several photometric redshift techniques including template fitting to the grism spectra, fitting the UV-to-radio photometry with galaxy templates plus a synchrotron model, fitting of the UV-to-near-infrared photometry with EAZY, and fitting the radio data alone with RAiSERed. For the host of GLEAM J0917-0012 we find a line at 1.12 micron and the UV-to-radio spectral energy distribution fitting favours solutions at z~2 or z~8. While this fitting shows a weak preference for the lower redshift solution, the models from the higher redshift solution are more consistent with the strength of the spectral line. The redshift constraint by RAiSERed of z>6.5 also supports the interpretation that this line could be Lyman-alpha at z=8.21; however EAZY favours the z~2 solution. We discuss the implications of both solutions. For the companion galaxy we find a line at 0.98 micron and the spectral energy distribution fitting favours solutions at z<3 implying that the line could be the [OII]3727 doublet at z=1.63 (although the EAZY solution is z~2.6+/-0.5). Further observations are still required to unambiguously determine the redshift of this intriguing candidate ultra-high-redshift radio galaxy (abridged).

Observations of synchrotron emission at low radio frequencies reveal a labyrinth of polarised Galactic structures. However, the explanation for the wealth of structures remains uncertain due to the complex interactions between the interstellar medium and the magnetic field. A multi-tracer approach to the analysis of large sky areas is needed. This paper aims to use polarimetric images from the LOFAR Two metre Sky Survey (LoTSS) to produce the biggest mosaic of polarised emission in the northern sky at low radio frequencies (150 MHz) to date. The large area this mosaic covers allows for detailed morphological and statistical studies of polarised structures in the high-latitude outer Galaxy, including the well-known Loop III region. We produced a 3100 square degree Faraday tomographic cube using a rotation measure synthesis tool. We calculated the statistical moments of Faraday spectra and compared them with data sets at higher frequencies (1.4 GHz) and with a map of a rotation measure derived from extragalactic sources. The mosaic is dominated by polarised emission connected to Loop III. Additionally, the mosaic reveals an abundance of other morphological structures, mainly {narrow and extended} depolarisation canals, which are found to be ubiquitous. We find a correlation between the map of an extragalactic rotation measure and the LoTSS first Faraday moment image. The ratio of the two deviates from a simple model of a Burn slab (Burn 1966) along the line of sight, which highlights the high level of complexity in the magnetoionic medium that can be studied at these frequencies.

Active galactic nucleus (AGN) feedback is postulated as a key mechanism for regulating star formation within galaxies. Studying the physical properties of the outflowing gas from AGN is thus crucial for understanding the co-evolution of galaxies and supermassive black holes. Here we report 55 pc resolution ALMA neutral atomic carbon [CI] $^3P_1\text{-}^3P_0$ observations toward the central 1 kpc of the nearby type-2 Seyfert galaxy NGC 1068, supplemented by 55 pc resolution CO($J=1\text{-}0$) observations. We find that [CI] emission within the central kpc is strongly enhanced by a factor of $>$5 compared to the typical [CI]/CO intensity ratio of $\sim$0.2 for nearby starburst galaxies (in units of brightness temperature). The most [CI]-enhanced gas (ratio $>$ 1) exhibits a kpc-scale elongated structure centered at the AGN that matches the known biconical ionized gas outflow entraining molecular gas in the disk. A truncated, decelerating bicone model explains well the kinematics of the elongated structure, indicating that the [CI] enhancement is predominantly driven by the interaction between the ISM in the disk and the highly inclined ionized gas outflow (which is likely driven by the radio jet). Our results strongly favor the "CO dissociation scenario" rather than the "in-situ C formation" one which prefers a perfect bicone geometry. We suggest that the high [CI]/CO intensity ratio gas in NGC 1068 directly traces ISM in the disk that is currently dissociated and entrained by the jet and the outflow, i.e., the "negative" effect of the AGN feedback.

The hot ($>10^6$ K) phase of the circumgalactic medium (CGM) contains a large fraction of baryons in galaxies. It also retains signatures of the processes that shaped the galaxies, such as feedback from active galactic nuclei (AGNs) and supernovae, and offers a unique and powerful way to constrain theoretical models of feedback. It is, however, notoriously difficult to detect. By stacking 2643 optically selected galaxies in the eROSITA Final Equatorial Depth Survey (eFEDS), we present spatially resolved properties of the large-scale CGM in both star-forming and quiescent galaxies spanning an order of magnitude in stellar mass. We mask out resolved point sources and galaxy groups/clusters and model the contribution from X-ray binaries and the hot ISM, producing accurate radial profiles of the CGM. We compare these data with mock X-ray observations of galaxy stacks in the IllustrisTNG100 (TNG) and EAGLE cosmological simulations. We detect extended emission from both the high-mass ($10.7<\log(M_*/M_\odot)<11.2$) and low-mass ($10.2<\log(M_*/M_\odot)<10.7$) galaxy stacks. Galaxies have somewhat more luminous CGM between $10-100$~kpc if they are more massive or star-forming. However, the CGM luminosity increases slower with stellar mass than predicted in simulations. In star-forming galaxies, the observed profiles are shallower than in simulations, suggesting that stellar feedback models are not efficient enough at blowing gas out. In quenched, high-mass galaxies the observed profiles are cuspier than observed, suggesting that AGN feedback models are too efficient at blowing gas out. Low-mass, quenched galaxies, where both supernova and AGN feedback models are inefficient, agree with simulations. Our results highlight the need to modify future prescriptions of galaxy formation models, particularly as they pertain to feedback and chemical enrichment of the CGM.

We use high-resolution maps of the molecular interstellar medium (ISM) in the centres of eighty-six nearby galaxies from the millimetre-Wave Interferometric Survey of Dark Object Masses (WISDOM) and Physics at High Angular Resolution in Nearby GalaxieS (PHANGS) surveys to investigate the physical mechanisms setting the morphology of the ISM at molecular cloud scales. We show that early-type galaxies tend to have smooth, regular molecular gas morphologies, while the ISM in spiral galaxy bulges is much more asymmetric and clumpy when observed at the same spatial scales. We quantify these differences using non-parametric morphology measures (Asymmetry, Smoothness and Gini), and compare these measurements with those extracted from idealised galaxy simulations. We show that the morphology of the molecular ISM changes systematically as a function of various large-scale galaxy parameters, including galaxy morphological type, stellar mass, stellar velocity dispersion, effective stellar mass surface density, molecular gas surface density, star formation efficiency and the presence of a bar. We perform a statistical analysis to determine which of these correlated parameters best predicts the morphology of the ISM. We find the effective stellar mass surface (or volume) density to be the strongest predictor of the morphology of the molecular gas, while star formation and bars maybe be important secondary drivers. We find that gas self-gravity is not the dominant process shaping the morphology of the molecular gas in galaxy centres. Instead effects caused by the depth of the potential well such as shear, suppression of stellar spiral density waves and/or inflow affect the ability of the gas to fragment.

Two scenarios have been proposed to match the existing observational constraints of the object HR 6819. The system could consist of a close inner B-type giant plus a black hole (BH) binary with an additional Be companion in a wide orbit. Alternatively, it could be a binary composed of a stripped B star and a Be star in a close orbit. Either scenario makes HR 6819 a cornerstone object as the stellar BH closest to Earth, or as an example of an important transitional, non-equilibrium phase for Be stars with solid evidence for its nature. We aimed to distinguish between the two scenarios for HR 6819. Both models predict two luminous stars but with very different angular separations and orbital motions. Therefore, the presence of bright sources in the 1-100 milliarcsec (mas) regime is a key diagnostic for determining the nature of the HR 6819 system. We obtained new high-angular resolution data with VLT/MUSE and VLTI/GRAVITY of HR 6819. The MUSE data are sensitive to bright companions at large scales, whilst the interferometric GRAVITY data are sensitive down to separations on mas scales and large magnitude differences. The MUSE observations reveal no bright companion at large separations and the GRAVITY observations indicate the presence of a stellar companion at an angular separation of ~1.2 mas that moves on the plane of the sky over a timescale compatible with the known spectroscopic 40-day period. We conclude that HR 6819 is a binary system and that no BH is present in the system. The unique nature of HR 6819, and its proximity to Earth make it an ideal system for quantitatively characterising the immediate outcome of binary interaction and probing how Be stars form.

Previous investigations have revealed that eccentric super-Earths represent a class of planets which are particularly effective at transporting minor bodies towards white dwarfs and subsequently polluting their atmospheres with observable chemical signatures. However, the lack of discoveries of these planets beyond a few astronomical units from their host stars prompts a better understanding of their orbital architectures from their nascent birth cluster. Here, we perform stellar cluster simulations of 3-planet and 7-planet systems containing super-Earths on initially circular, coplanar orbits. We adopt the typical stellar masses of main-sequence progenitors of white dwarfs ($1.5\,\mathrm{M}_{\odot}$-$2.5\,\mathrm{M}_{\odot}$) as host stars and include 8,000 main-sequence stars following a Kroupa initial mass function in our clusters. Our results reveal that about 30 per cent of the simulated planets generate eccentricities of at least 0.1 by the time of cluster dissolution, which would aid white dwarf pollution. We provide our output parameters to the community for potential use as initial conditions for subsequent evolution simulations.

M54 is a prototype of a globular cluster embedded in a dark matter halo. Gaia EDR3 photometry and proper motions separate the old, metal-poor stars from the more metal rich and younger dwarf galaxy stars. The metal poor stars dominate the inner 50 pc, with a velocity dispersion profile that declines to a minimum around 30 pc then rises back to nearly the central velocity dispersion, as expected for a globular cluster at the center of a dark matter halo. The Jeans analysis of the three separate stellar populations give consistent masses that rise approximately linearly with radius to 1 kpc, implying a small core or cuspy halo. These data are compatible with an infalling CDM dark matter halo reduced to 3x10^8 M_sun at the 50 kpc apocenter 2.3 Gyr ago, with a central globular cluster surrounded by the remnant of a dwarf galaxy. Tides gradually remove material beyond 1 kpc but have little effect on the stars and dark matter within 300 pc of the center. M54 appears to be a transitional system between globular clusters with and without local dark halos, whose evolution within the galaxy depends sensitively on the time of accretion and orbital pericenter.

Supra-arcade downflows (SADs) are infrequent, wiggly opaque structures observed to descend through the solar corona, mostly in EUV and soft X-ray frequencies. From their physical characteristics, SADs have been interpreted as voided (subdense) bubbles and are related to magnetic reconnection processes during long-term erupting flares. In this work we use numerical MHD simulations to compute flux density maps, which are convolved with telescope beams to synthesise images with the aim to assess the expected SADs emission at radio wavelengths and propose observing strategies, including the instruments that can be used. We assume that the emission is thermal bremsstrahlung from a fully ionised plasma without any appreciable gyroresonance contribution since magnetic fields are of the order of ~10 G. We find that SADs emission should be optically thin in the frequency [10-1000] GHz range, and the spatially integrated flux should be larger than 1 Jy. We conclude, therefore, that observing SADs in radio frequencies between [0.5-1000] GHz is feasible with present instrumentation. Moreover, since the emission is for the most part optically thin, the flux density is proportional to temperature, density and line-of-sight depth, and when combined with EUV and soft X-ray images, may allow a better density and temperature determination of SADs.

All disc-accreting astrophysical objects produce powerful outflows. In binaries containing neutron stars (NS) or black holes, accretion often takes place during violent outbursts. The main disc wind signatures during these eruptions are blue-shifted X-ray absorption lines, which are preferentially seen in disc-dominated "soft states". By contrast,optical wind-formed lines have recently been detected in "hard states", when a hot corona dominates the luminosity. The relationship between these signatures is unknown, and no erupting system has revealed wind-formed lines between the X-ray and optical bands yet, despite the many strong resonance transitions in this ultraviolet (UV) region. Here, we show that the transient NS binary Swift J1858.6-0814 exhibits wind-formed, blue-shifted absorption associated with C IV, N V and He II in time-resolved UV spectroscopy during a luminous hard state. This represents the first evidence for a warm, moderately ionized outflow component in this state. Simultaneously observed optical lines also display transient blue-shifted absorption. Decomposing the UV data into constant and variable components, the blue-shifted absorption is associated with the former. This implies that the outflow is not connect to the luminous flares in the data. The joint presence of UV and optical wind features reveals a multi-phase and/or stratified outflow from the outer disc. This type of persistent mass loss across all accretion states has been predicted by radiation-hydrodynamic simulations and helps to explain the shorter-than-expected outbursts duration.

We present the second part (days 19 to 102) of our high-resolution Echelle atlas of the spectral evolution of RS Oph during the 2021 nova outburst. Some quick conclusions that can be drawn are: (1) the spectra evolution has progressed smoothly, suggesting that both the expansion of the ejecta through the pre-existing RG wind and the photo-ionization from the central source did not encountered sudden discontinuities prior to switch-off of the nuclear burning; (2) the 2021 spectral evolution is similar to that of the 2006 eruption, with a mirror-appearance of the triple-peaked line profiles (blue peak brighter than red in 2006, the opposite in 2021); (3) the evolution of the triple-peaked profiles progressed similarly through different excitation and ionization degrees, suggesting an origin in a bipolar outflow nested to an equatorial torus; (4) the time-behavior of coronal emission lines appears quite smooth, both in terms of radiated flux as well as evolution of the triple-peaked profiles, with a plateau at maximum lasting up to about day 87; (5) there seems to exist within RS Oph two distinct types of ejecta: one fast moving and producing Gaussian-like line profiles (primarely [NII] and [OIII]), that keep expanding at FWHM=1000 km/s with no sign of ongoing deceleration during the two months leading up to Solar conjunction, and slower moving ejecta giving rise to the triple-peaked line profiles (HeI, HeII, coronal lines, and also Balmer lines at later times), for which the separation in velocity of blue and red peaks is still shrinking at the time of Solar conjunction when it reaches an average of 330 km/s; (6) the "symbiotic band" at 6825 Ang, due to Raman scattering of OVI 1032 by neutral hydrogen, appears simultaneously with the coronal lines: its very presence indicates that a sizeable part of the wind of the red giant has remained neutral during the development of coronal lines.

The practical distinctions between asteroids and comets, viewed as products of accretion on either side of the snow line, are less clear-cut than previously understood. In this chapter, we discuss the numerous solar system populations which have physical and dynamical properties that conflict with any simple diagnosis of their nature and origin. Studies of these so-called "continuum" or "transition objects", which include many of the most intriguing bodies in the solar system, have implications for a broad range of scientific topics from the demise of comets and the activation of asteroids to the production of interplanetary debris and the origin of the terrestrial planet volatiles. We present an overview of the current state of knowledge concerning the asteroid-comet continuum and discuss the numerous physical processes behind the activity shown by small bodies in the solar system.

Ahead of the United States' crewed return to the moon in 2024, Intuitive Machines, under a NASA Commercial Lunar Payload Services contract, will land their Nova-C lunar lander in October 2021. At 30 meters altitude during the terminal descent, EagleCam will be deployed, and will capture and transmit the first-ever third-person images of a spacecraft making an extraterrestrial landing. This paper will focus on the structural design, modeling, and impact analysis of a 1.5U CubeSat payload to withstand a ballistic, soft-touch landing on the lunar surface.

We present The MKID Pipeline, a general use science data pipeline for the reduction and analysis of ultraviolet, optical and infrared (UVOIR) Microwave Kinetic Inductance Detector (MKID) data sets. This paper provides an introduction to the nature of MKID data sets, an overview of the calibration steps included in the pipeline, and an introduction to the implementation of the software.

The comparison of the intrinsic luminosity of extragalactic objects detected at different redshifts requires the evaluation of the observed flux in a common rest-frame emission band. It is usually performed by adding a term to the observed magnitudes denominated K-correction. In this paper, the K-correction transformation has been carefully scrutinized. As a result of this inspection, it has been found that K-corrections systematically downgrades the measured energy flux by a factor of (1+z). C-correction, a more appropriate rest-frame magnitude transformation has been derived. C-correction has been applied to the re-analysis of JLA (Joint Light-curve analysis) supernovae sample. The best cosmological fit to JLA type Ia supernovae (SNIa) sample calibrated with CMB Planck $h=0.675$ for the common luminosity-angular distances relation $L2:D_L=D_A(1+z)^2$ corresponds to a closed Universe $L2(h_0=0.48,\Omega_M=2.25$). Nevertheless, considering the new luminosity-angular distances relation $L1:D_L=D_A(1+z)$ that accounts for expansion lensing, the new SNIa analysis converges to $L1(h_0=0.7558,\Omega_M=0.035,\Omega_\Lambda=0.965)$, providing a coherent explanation of the Hubble constant ($H_0$) tension between Planck and SH0ES ($h=0.7422 \pm 0.0182$) results.

CO is one of the most abundant ice components on interstellar dust grains. When it is mixed with amorphous solid water (ASW) or located on its surface, an absorption band at 2152 cm$^{-1}$ is always present in laboratory measurements. This spectral feature is attributed to the interaction of CO with dangling-OH bonds (dOH) in ASW. However, this band is absent in observational spectra of interstellar ices. This raises the question whether CO forms a relatively pure layer on top of ASW or is in close contact with ASW, but not via dangling bonds. We performed laboratory experiments to simulate the layered structure of the ice mantle, i.e., we grew CO ice on top of: 1) Pure ASW; 2) NH$_3$:H$_2$O=10:100 mixed ice; and 3) CO$_2$:H$_2$O=20:100 mixed ice. We found that annealing the ice reduces the 2152 cm$^{-1}$ band, and that NH$_3$ blocks the dOH on ASW surface and therefore reduces the 2152 cm$^{-1}$ band more effectively than CO$_2$. We suggest that this difference between NH$_3$ and CO$_2$ can be ascribed to the polarity of the guest molecule (NH$_3$ is a polar species, whereas CO$_2$ is apolar). The polarity implies that the formation of an H-bond between the N atom of ammonia and the dOH is a barrier-less reaction. We also determined the pore surface area of the ice mixtures as a function of the annealing temperature, and found that the non-detection of 2152 cm$^{-1}$ band does not necessarily exclude the possibility of a porous ice mantle.

The advent of new all-sky radio surveys such as the VLA Sky Survey (VLASS) and the Rapid ASKAP Continuum Survey (RACS), performed with the latest generation radio telescopes, is opening new possibilities on the classification and study of extragalactic $\gamma$-ray sources, specially the underrepresented ones like radio galaxies. In particular, the enhanced sensitivity (sub-mJy level) and resolution (a few arcsec) provides a better morphological and spectral classification. In this work, we present the reclassification of a Fermi-LAT source as a new FRII radio galaxy from the INTEGRAL sample found to emit at GeV energies. Through a broad-band spectral fitting from radio to $\gamma$-ray, we find that the commonly invoked jet contribution is not sufficient to account for the observed $\gamma$-ray flux. Our modeling suggests that the observed emission could mainly originate in the lobes (rather than in the radio core) by inverse Compton scattering of radio-emitting electrons off the ambient photon fields. In addition, we cross-correlated the latest generation radio surveys with a list of Fermi-LAT objects from the literature considered to be candidate misaligned AGN, finding four new radio galaxies with a double-lobed morphology. Additional four objects could be classified as such thanks to previous studies in the literature, for a total of nine new radio galaxies with GeV emission presented in this work. We foresee that further objects of this class might be found in the near future with the advent of the Square Kilometer Array (SKA), populating the GeV sky.

We consider the bias introduced by a spatially-varying multiplicative shear bias (m-bias) on tomographic cosmic shear angular power spectra. To compute the bias in the power spectra, we estimate the mode-coupling matrix associated with an m-bias map using a computationally-efficient pseudo-Cl method. This allows us to consider the effect of the m-bias to high l. We then conduct a Fisher matrix analysis to forecast resulting biases in cosmological parameters. For a Euclid-like survey with a spatially-varying m-bias, with zero mean and rms of 0.01, we find that parameter biases reach a maximum of ~10% of the expected statistical error, if multipoles up to l_max = 5000 are included. We conclude that the effect of the spatially-varying m-bias may be a sub-dominant but potentially non-negligible contribution to the error budget in forthcoming weak lensing surveys. We also investigate the dependence of parameter biases on the amplitude and angular scale of spatial variations of the m-bias field, and conclude that requirements should be placed on the rms of spatial variations of the m-bias, in addition to any requirement on the mean value. We find that, for a Euclid-like survey, biases generally exceed ~30% of the statistical error for m-bias rms 0.02 - 0.03 and can exceed the statistical error for rms ~0.04 - 0.05. This allows requirements to be set on the permissible amplitude of spatial variations of the m-bias that will arise due to systematics in forthcoming weak lensing measurements.

Tidal disruption events (TDEs) offer a unique way to study dormant black holes. While the number of observed TDEs has grown thanks to the emergence of wide-field surveys in the past few decades, questions regarding the nature of the observed optical, UV, and X-ray emission remain. We present a uniformly selected sample of 30 spectroscopically classified TDEs from the Zwicky Transient Facility Phase I survey operations with follow-up \textit{Swift} UV and X-ray observations. Through our investigation into correlations between light curve properties, we recover a shallow positive correlation between the peak bolometric luminosity and decay timescales. We introduce a new spectroscopic class of TDE, TDE-featureless, which are characterized by featureless optical spectra. We recover a difference between the light curve timescales of the four TDE classes, with TDE-H+He having longer rise times. The new TDE-featureless class shows larger peak bolometric luminosities, peak blackbody temperatures, and peak blackbody radii. We examine the differences between the X-ray bright and X-ray faint populations of TDEs in this sample, finding that X-ray bright TDEs show longer decay timescales than the X-ray faint sub-sample. We present a scheme for classifying the shapes of the optical light curves, finding that TDEs with "structured" light curves have significantly less massive host galaxies than other types of TDEs. This sample of optically selected TDEs is the largest sample of TDEs from a single survey yet, and the systematic discovery, classification, and follow-up of this sample allows for robust characterization of TDE properties, an important stepping stone looking forward toward the Rubin era.

Secondary eclipse observations of hot Jupiters can reveal both their compositions and thermal structures. Previous observations have shown a diversity of hot Jupiter eclipse spectra, including absorption features, emission features, and featureless blackbody-like spectra. We present a secondary eclipse spectrum of the hot Jupiter WASP-77Ab observed between $1-5$ $\mu$m with the Hubble Space Telescope (HST) and the Spitzer Space Telescope. The HST observations show signs of water absorption indicative of a non-inverted thermal structure. We fit the data with both a one-dimensional free retrieval and a grid of one-dimensional self-consistent forward models to confirm this non-inverted structure. The free retrieval places a $3\sigma$ lower limit on the atmospheric water abundance of $\log(n_\mathrm{H_2O})>-4.78$ and can not constrain the CO abundance. The grid fit produces a slightly super-stellar metallicity and constrains the carbon-to-oxygen ratio to less than or equal to the solar value. We also compare our data to recent high-resolution observations of WASP-77Ab taken with the Gemini-South/IGRINS spectrograph and find that our observations are consistent with the best-fit model to the high-resolution data. However, the metallicity derived from the IGRINS data is significantly lower than that derived from our self-consistent model fit. We find that this difference may be due to disequilibrium chemistry, and the varying results between the models applied here demonstrate the difficulty of constraining disequilibrium chemistry with low-resolution, low wavelength coverage data alone. Future work to combine observations from IGRINS, HST, and JWST will improve our estimate of the atmospheric composition of WASP-77Ab.

We report the first-time recovery of a fresh meteorite fall using a drone and a machine learning algorithm. A fireball on the 1st April 2021 was observed over Western Australia by the Desert Fireball Network, for which a fall area was calculated for the predicted surviving mass. A search team arrived on site and surveyed 5.1 km2 area over a 4-day period. A convolutional neural network, trained on previously-recovered meteorites with fusion crusts, processed the images on our field computer after each flight. meteorite candidates identified by the algorithm were sorted by team members using two user interfaces to eliminate false positives. Surviving candidates were revisited with a smaller drone, and imaged in higher resolution, before being eliminated or finally being visited in-person. The 70 g meteorite was recovered within 50 m of the calculated fall line using, demonstrating the effectiveness of this methodology which will facilitate the efficient collection of many more observed meteorite falls.

Understanding effects driven by rotation in the solar convection zone is essential for many problems related to solar activity, such as the formation of differential rotation, meridional circulation, and others. We present realistic 3D radiative hydrodynamics simulations of solar subsurface dynamics in the presence of rotation in a local domain 80 Mm-wide and 25 Mm deep, located at 30 degrees latitude. The simulation results reveal the development of a shallow 10-Mm deep near-surface shear layer ("leptocline"), characterized by a strong radial rotational gradient and self-organized meridional flows. This shear layer is located in the hydrogen ionization zone associated with enhanced anisotropic convective flows overshooting into a relatively stable zone between the H and HeII ionization zones. The radial variations of the differential rotation and meridional circulation profiles obtained from the simulations agree with helioseismic observations, indicating that a major role in forming the leptocline and subsurface meridional flows is played by the local Reynolds stresses.

The mobility of lighter species on the surface of interstellar dust grains plays a crucial role in forming simple through complex molecules. Carbon monoxide is one of the most abundant molecules, its surface diffusion on the grain surface is essential to forming many molecules. Recent laboratory experiments found a diverse range of diffusion barriers for CO on the grain surface, their use can significantly impact the abundance of several molecules. The impact of different diffusion barriers of CO, in the astrochemical models, is studied to understand its effect on the abundance of solid CO and the species for which it is a reactant partner. A gas-grain network is used for three different physical conditions; cold-core and warm-up models with slow and fast heating rates. Two different ratios (0.3 and 0.5) between diffusion and desorption barrier are utilized for all the species. For each physical condition and ratio, six different models are run by varying diffusion barriers of CO. Solid CO abundance for the models with the lowest diffusion barrier yields less than 0.1% of water ice for cold clouds and a maximum of 0.4% for slow and fast warm-up models. Also, solid CO$_2$ in dense clouds is significantly overproduced (140 % of water). The abundance of H$_2$CO and CH$_3$OH showed an opposite trend, and HCOOH, CH$_3$CHO, NH$_2$CO, and CH$_3$COCH$_3$ are produced in lower quantities for models with low diffusion barriers for CO. Considerable variation in abundance is observed between models with the high and low diffusion barrier. Models with higher diffusion barriers provide a relatively better agreement with the observed abundances when compared with the models having lower diffusion barriers.

The interaction of turbulence, magnetic fields, self-gravity, and stellar feedback within molecular clouds is crucial for understanding star formation. We study the effects of self-gravity and outflow feedback on the properties of the turbulent velocity via the structure function over length scales from $\sim$ 0.01 pc to 2 pc. We analyze a series of three-dimensional, magnetohydrodynamical (MHD) simulations of star cluster formation, including self-gravity, turbulence, magnetic fields, stellar radiative heating, and outflow feedback. We observe that self-gravity and protostellar outflows increase the velocity fluctuations over all length scales. In particular, outflows can amplify the velocity fluctuations by up to a factor of $\sim$7 on scales $\sim$ 0.01 - 0.2 pc and drive turbulence up to a scale of $\sim$ 1 pc. The amplified velocity fluctuations provide more support against gravity and enhance fragmentation on small scales. The role of self-gravity is more significant on smaller dense clumps and it increases the fraction of the compressive velocity component up to a scale of $\sim$ 0.2 pc. However, outflow feedback drives both solenoidal and compressive modes, but it induces a higher fraction of solenoidal modes relative to compressive modes. Thus, with outflows, the dense core ends up with a slightly higher fraction of solenoidal modes. We find that the compressible fraction is fairly constant with about 1/3 on scales $\sim$ 0.1 - 0.2 pc. The combined effect of enhanced velocity dispersion and reduced compressive fraction contributes to a reduction in the star formation rate compared to when outflow feedback is not included.

The predicted size of dark matter substructures in kilo-parsec scales is model-dependent. Therefore, if the correlations between dark matter mass densities, as a function of the distances between them, are measured via observations, we can scrutinize dark matter scenarios. In this paper, we use Gaia's data to infer the single-body phase-space density of the stars in the Fornax dwarf spheroidal galaxy. The latter together with the Jeans equation, after eliminating the gravitational potential using the Poisson equation, reveals the mass density of dark matter as a function of its position in the galaxy. We derive the correlations between dark matter mass densities as a function of distances between them. For distances beyond 100 parsec, no statistically significant correlation is observed. We use the results to limit the parameter space of the Ginzburg-Landau statistical field theory of dark matter mass densities, and subsequently shrink the parameter space of (i) a classic gas dark matter and (ii) a superfluid dark matter. Our results can be used to leave bounds on high-resolution N-body simulations by comparing the correlations computed using their outputs and the measured ones via observations.

The gamma-ray halo around Geminga is formed owing to the slow diffusion of the electrons released by the Geminga PWN. The latest HAWC and HESS observations exhibit complex features in the TeV gamma-ray spectrum of the Geminga halo. We first show that the new results cannot be interpreted by the commonly used simple model, where a single power-law injection spectrum and an energy index of $\delta=1/3$ for the diffusion coefficient are assumed. We then propose a two-population electron injection model based on the x-ray observations of the Geminga PWN, which consists of a population of freshly accelerated electrons escaping from the PWN through rapid outflows and a population trapped longer inside the PWN before escaping. The two-population model interprets the HAWC and HESS data well, and the goodness of fit improves significantly compared with the single power-law injection model. It also predicts a different energy dependency of the gamma-ray profile from the single power-law model, which could be tested by LHAASO in the coming future. We note that a $\delta$ slightly larger than 1 is needed to fit the HAWC and HESS data consistently. We also discuss the possible improvements by adopting the two-zone diffusion model.

Earth's magnetosphere extension is controlled by solar activity level via solar wind properties. Understanding such a relation in the Solar System is useful to predict the condition of exoplanetary magnetosphere near Sun-like stars. We use measurements of a chromospheric proxy, the Ca II K index, and solar wind OMNI parameters to connect the solar activity variations on the decennial time scales to solar wind properties. The dataset span over the time interval 1965-2021, which almost entirely covers the last 5 solar cycles. Using both cross-correlation and mutual information analysis, a 3.2-year lag of the solar wind speed with respect to the Ca II K index is found. Analogously, a 3.6-year lag is found with respect to the dynamic pressure. A correlation between the solar wind dynamic pressure and the solar UV emission is therefore found and used to derive the Earth's magnetopause standoff distance. Moreover, the advantage of using a chromospheric proxy, such as the Ca II K index, opens the possibility to extend the relation found for the Sun to Sun-like stars, by linking stellar variability to stellar wind properties. The model is applied to a sample of Sun-like stars as a case study, where we assume the presence of an Earth-like exoplanet at 1 AU. Finally, we compare our results with previous estimates of the magnetosphere extension for the same set of sun-like stars.

We used the INPOP19a planetary ephemerides to perform the orbital adjustment of 14099 asteroids based on Gaia-DR2 observations, and compare for 23 of them the resulting orbits to radar data. As Gaia-DR2 has been processed using the planetary ephemeris INPOP10e, the primary goal of this paper is to confirm the portability of the data when using an updated version of the solar system model. In particular, we point out the fact that the Gaia satellite positions -- provided with respect to the INPOP10e solar system barycenter -- must be corrected when using another planetary ephemeris . We also present a convenient least square formalism that only handles small matrices and allows the adjustment of global parameters, such as masses. In order to check the consistency of the Gaia observations with other types of observations, we perform an orbital adjustment in combining Gaia and radar range observations for 23 objects, together with a careful post-fit analysis including an estimation of the Gaia systematic errors. Finally, we show that to ensure the combined use of Gaia angular DR2 observations and radar ranging, a more developed than firstly proposed dynamical modeling is required together with the addition of the systematic Gaia bias in the fit procedure.These results give promising directions for the next Gaia delivery, Gaia-DR3.

In this paper we revisit potential biases in cosmic shear power spectra caused by bias terms that multiply quadratic powers of the shear. We identify spin-2 ($m_2$), spin-(-2) ($m_{-2}$) and spin-6 ($m_6$) multiplicative biases and propagate these biases into shape measurement statistics and the cosmic shear power spectrum. We find that such biases can be measured by performing regression on calibration data. We find that that the impact on the power spectrum is an additional bispectrum dependency. Ignoring quadratic terms can lead to biases in cosmological parameters of up to $2(m_2+m_{-2}-m_6)0.4\sigma$ for Stage-IV dark energy experiments, and a reasonable requirement on the amplitude of these terms is at least $|m_2+m_{-2}-m_6|\leq 0.031$. If one ignores quadratic terms entirely and instead fits only a linear dependence to calibration data then this can in addition cause spurious multiplicative and additive biases which leads to a requirement on the uncertainty of these terms of $\sigma(m_s)\simeq 5\times 10^{-3}$, if ignored in regression. In future, Stage-IV dark energy experiments should seek to measure and minimise quadratic bias terms, and/or decrease susceptibility to them by using methods to remove small-scale sensitivity.

Three-dimensional (3D) maps of Galactic interstellar dust are a tool for a wide range of uses. We aim to construct 3D maps of dust extinction in the Local Arm and surrounding regions. Gaia EDR3 photometric data were combined with 2MASS measurements to derive extinction towards stars with accurate photometry and relative uncertainties on parallaxes of less than 20%. We applied our hierarchical inversion algorithm adapted to inhomogeneous spatial distributions of target stars to this catalogue of extinctions. We present the updated 3D dust extinction distribution and provide an estimate of the error on integrated extinctions from the Sun to each area in the 3D map. The computational area is similar to the one of the previous DR2 map, a 6 kpc x 6 kpc x 0.8 kpcAstrophysics volume around the Sun. Due to the addition of fainter target stars, the volume in which the clouds can be reconstructed has increased. Due to the improved accuracy of the parallaxes and photometric data in EDR3, extinctions among neighbouring targets are more consistent, allowing one to reach an increased contrast in the dense areas, while cavity contours are more regular. We show several comparisons with recent results on dust and star distributions. The wavy pattern around the Plane of the dust concentrations is better seen and exists over large regions. Its mean vertical peak-to-peak amplitude is of the order of 300 pc; interestingly, it is similar to the vertical period of the spectacular snail-shaped stellar kinematical pattern discovered in Gaia data. The Gaia EDR3 catalogue allows for a significant improvement of the extinction maps to be made and the hierarchical technique confirms its efficiency for massive datasets. Future comparisons between 3D maps of interstellar matter and stellar distributions may help to understand which mergers or internal perturbations have shaped the Galaxy within the first 3 kpc.

Euclid will be the first space mission to survey most of the extragalactic sky in the 0.95-2.02 $\mu$m range, to a 5$\sigma$ point-source median depth of 24.4 AB mag. This unique photometric data set will find wide use beyond Euclid's core science. In this paper, we present accurate computations of the Euclid Y_E, J_E and H_E passbands used by the Near-Infrared Spectrometer and Photometer (NISP), and the associated photometric system. We pay particular attention to passband variations in the field of view, accounting among others for spatially variable filter transmission, and variations of the angle of incidence on the filter substrate using optical ray tracing. The response curves' cut-on and cut-off wavelengths - and their variation in the field of view - are determined with 0.8 nm accuracy, essential for the photometric redshift accuracy required by Euclid. After computing the photometric zeropoints in the AB mag system, we present linear transformations from and to common ground-based near-infrared photometric systems, for normal stars, red and brown dwarfs, and galaxies separately. A Python tool to compute accurate magnitudes for arbitrary passbands and spectral energy distributions is provided. We discuss various factors from space weathering to material outgassing that may slowly alter Euclid's spectral response. At the absolute flux scale, the Euclid in-flight calibration program connects the NISP photometric system to Hubble Space Telescope spectrophotometric white dwarf standards; at the relative flux scale, the chromatic evolution of the response is tracked at the milli-mag level. In this way, we establish an accurate photometric system that is fully controlled throughout Euclid's lifetime.

In this review we discuss all the relevant solar/stellar radiation and plasma parameters and processes that act together in the formation and modification of atmospheres and exospheres that consist of surface-related minerals. Magma ocean degassed silicate atmospheres or thin gaseous envelopes from planetary building blocks, airless bodies in the inner Solar System, and close-in magmatic rocky exoplanets such as CoRot-7b, HD219134b and 55 Cnc e are addressed. The depletion and fractionation of elements from planetary embryos, which act as the building blocks for protoplanets are also discussed. In this context the formation processes of the Moon and Mercury are briefly reviewed. The Lunar surface modification since its origin by micrometeoroids, plasma sputtering, plasma impingement as well as chemical surface alteration and the search of particles from the early Earth's atmosphere that were collected by the Moon on its surface are also discussed. Finally, we address important questions on what can be learned from the study of Mercury's environment and its solar wind interaction by MESSENGER and BepiColombo in comparison with the expected observations at exo-Mercurys by future space-observatories such as the JWST or ARIEL and ground-based telescopes and instruments like SPHERE and ESPRESSO on the VLT, and vice versa.

Our knowledge of the heating mechanisms that are at work in the chromosphere of plage regions remains highly unconstrained from observational studies. The purpose of our study is to estimate the chromospheric heating terms from a plage dataset, characterize their spatio-temporal distribution and set constraints on the heating processes that are at work. We make use of NLTE inversions to infer a model of the photosphere and chromosphere of a plage dataset acquired with the Swedish 1-m Solar Telescope. We use this model atmosphere to calculate the chromospheric radiative losses from H i, Ca ii and Mg ii atoms. We approximate the chromospheric heating terms by the net radiative losses predicted by the inverted model. In order to make the analysis of time-series over a large field-of-view computationally tractable, we make use of a neural network. In the lower chromosphere, the contribution from the Ca ii lines is dominant and located in the surroundings of the photospheric footpoints. In the upper chromosphere, the H i contribution is dominant. Radiative losses in the upper chromosphere form an homogeneous patch that covers the plage region. The net radiative losses can be split in a periodic component with an average amplitude of ampQ = 7.6 kW m^{-2} and a static (or very slowly evolving) component with a mean value of -26.1 kW m^{-2}. Our interpretation is that in the lower chromosphere, the radiative losses are tracing the sharp lower edge of the hot magnetic canopy, where the electric current is expected to be large. In the upper chromosphere, both the magnetic field and the distribution of net radiative losses are room-filling, whereas the amplitude of the periodic component is largest. Our results suggest that acoustic wave heating may be responsible for one third of the energy deposition in the upper chromosphere, whereas other heating mechanisms must be responsible for the rest.

We report observations of five dwarf galaxies in the vicinity of the luminous S0 galaxy NGC 3115 performed with the Advanced Camera for Surveys on the Hubble Space Telescope. Their distances determined via the Tip of the Red Giant Branch are: 10.05 Mpc (UGCA 193), 9.95 Mpc (KKSG 17), 10.13 Mpc (2MASX-J0957-0915), 10.42 Mpc (2dFGRS-TGN218Z179) and 11.01 Mpc (KKSG 19). With their typical distance error of about 0.75 Mpc all the five dwarfs are consistent to be true satellites of the host galaxy NGC 3115 (10.2$\pm$0.2 Mpc). Using the DESI Legacy Imaging Surveys we also found 5 new probable dwarf satellites of NGC3115, as well as 4 new probable members of the neighboring group around NGC 3521 situated 3 Mpc away from the NGC 3115 group. Based on the radial velocities and projected separations of 10 dwarf companions, we derived the total (orbital) mass of NGC 3115 to be (4.89$\pm$1.48) $10^{12}$ $M_{\odot}$. The ratio of the total mass-to-K-luminosity of NGC 3115 is (50$\pm$15) $M_\odot/L_\odot$, which is typical for the early-type luminous galaxies.

Establishing a large sample of post common envelope binaries (PCEBs) that consist of a white dwarf plus an intermediate mass companion star of spectral type AFGK, offers the potential to provide new constraints on theoretical models of white dwarf binary formation and evolution. Here we present a detailed analysis of two new systems, TYC 110-755-1 and TYC 3858-1215-1. Based on radial velocity measurements we find the orbital periods of the two systems to be $\sim$ 0.85 and $\sim$ 1.64 days, respectively. In addition, HST spectroscopy of TYC 110-755-1 allowed us to measure the mass of the white dwarf in this system (0.78 M$_\odot$). We furthermore analysed TESS high time resolution photometry and find both secondary stars to be magnetically extremely active. Differences in the photometric and spectroscopic periods of TYC 110-755-1 indicate that the secondary in this system is differentially rotating. Finally, studying the past and future evolution of both systems, we conclude that the common envelope efficiency is likely similar in close white dwarf plus AFGK binaries and PCEBs with M-dwarf companions and find a wide range of possible evolutionary histories for both systems. While TYC 3858-1215-1 will run into dynamically unstable mass transfer that will cause the two stars to merge and evolve into a single white dwarf, TYC 110-755-1 is a progenitor of a cataclysmic variable system with an evolved donor star.

The system of two transiting Neptune-sized planets around the bright, young M-dwarf AU Mic provides a unique opportunity to test models of planet formation, early evolution, and star-planet interaction. However, the intense magnetic activity of the host star makes measuring the masses of the planets via the radial velocity (RV) method very challenging. We report on a 1-year, intensive monitoring campaign of the system using 91 observations with the HARPS spectrograph, allowing for detailed modelling of the $\sim 600\, \rm{m\,s^{-1}}$ peak-to-peak activity-induced RV variations. We used a multidimensional Gaussian Process framework to model these and the planetary signals simultaneously. We detect the latter with semi-amplitudes of $\rm{K_b} = 5.8 \pm 2.5\, \rm{m\,s^{-1}}$ and $\rm{K_c} = 8.5 \pm 2.5\, \rm{m\, s^{-1}}$, respectively. The resulting mass estimates, $\rm{M_b} = 11.7 \pm 5.0\, \rm{M_{\rm \oplus}}$ and $\rm{M_c} = 22.2 \pm 6.7\, \rm{M_{\rm \oplus}}$, suggest that planet b might be less dense, and planet c considerably denser than previously thought. These results are in tension with the current standard models of core-accretion. They suggest that both planets accreted a H/He envelope that is smaller than expected, and the trend between the two planets' envelope fractions is the opposite of what is predicted by theory.

We employed the {\it VISTA near-infrared $YJK_\mathrm{s}$ survey of the Magellanic System} (VMC), to analyse the $Y,\,J,\,K_\mathrm{s}$ light curves of $\delta$ Cepheid stars (DCEPs) in the Large Magellanic Cloud (LMC). Our sample consists of 4408 objects accounting for 97 per cent of the combined list of OGLE\,IV and {\it Gaia}\,DR2 DCEPs. We determined a variety of period-luminosity ($PL$) and period-Wesenheit $PW$ relationships for Fundamental (F) and First Overtone (1O) pulsators. We discovered for the first time a break in these relationships for 1O DCEPs at $P$=0.58 d. We derived relative individual distances for DCEPs in the LMC with a precision of $\sim$1 kpc, calculating the position angle of the line of nodes and inclination of the galaxy: $\theta$=145.6$\pm$1.0 deg and $i$=25.7$\pm$0.4 deg. The bar and the disc are seen under different viewing angles. We calculated the ages of the pulsators, finding two main episodes of DCEP formation lasting $\sim$40 Myr which happened 93 and 159 Myr ago. Likely as a result of its past interactions with the SMC, the LMC shows a non-planar distribution, with considerable structuring: the bar is divided into two distinct portions, the eastern and the western displaced by more than 1 kpc from each other. Similar behaviour is shown by the spiral arms. The LMC disc appears "flared" and thick, with a disc scale height of $h\sim 0.97$ kpc. This feature can be explained by strong tidal interactions with the Milky Way and/or the Small Magellanic Cloud or past merging events with now disrupted LMC satellites.

Several diffuse interstellar bands (DIBs) have profiles with resolved sub-peaks that resemble rotational bands of large molecules. Analysis of these profiles can constrain the sizes and geometries of the DIB carriers, especially if the profiles exhibit clear variations along lines of sight probing different physical conditions. Using the extensive data set from the EDIBLES survey we searched for systematic variations in the peak-to-peak separation of these sub-peaks for the $\lambda\lambda$5797, 6379, and 6614 DIBs in lines of sight with a single dominant interstellar cloud. We used the spectra of twelve single-cloud sight lines to measure the peak-to-peak separation in the band profile substructures for these DIBs. We adopted the rotational contour formalism to infer the rotational constant for each DIB carrier and the rotational excitation temperature in the sight lines. We compared these to rotational constants for linear and spherical molecules to estimate the DIB carrier sizes. All three DIBs have peak separations that vary systematically between lines of sight, indicating correlated changes in the rotational excitation temperatures. We derived $B_{6614}$=$(22.2\pm8.9)\times 10^{-3}$ cm$^{-1}$, consistent with previous estimates. Assuming a similar rotational temperature for the $\lambda$6614 DIB carrier and assuming a linear carrier, we found B$_{5797}^{\rm linear}=(5.1\pm2.0)\times10^{-3}~{\rm cm}^{-1}$ and B$_{6379}^{\rm linear} =(2.3\pm0.9)\times10^{-3}~{\rm cm}^{-1}$. If the carriers of those DIBs however are spherical species, their rotational constants are half that value, $B_{5797}^{\rm spherical} = (2.6\pm1.0)\times10^{-3}~{\rm cm}^{-1}$ and $B_{6379}^{\rm spherical} = (1.1\pm0.4)\times10^{-3}~{\rm cm}^{-1}$. We estimate molecule sizes that range from 7--9 carbon atoms ($\lambda$6614 carrier, linear) to 77--114 carbon atoms ($\lambda$6379, spherical).

I show that the quadrupole of the Cosmic Microwave Background (CMB) evolves more rapidly than previously expected, as a result of the acceleration of the Sun towards the Galactic center. The acceleration, measured most recently by Gaia EDR3, implies a fractional change in the quadrupole of ~10^{-9} per year, an order of magnitude larger than expected from the evolution in the last scattering surface of the CMB.

The Gaussian linear model provides a unique way to obtain the posterior probability distribution as well as the Bayesian evidence analytically. Considering the expansion rate data, the Gaussian linear model can be applied for $\Lambda$CDM, wCDM, and a non-flat $\Lambda$CDM. In this paper, we simulate the expansion data with various precision and obtain the Bayesian evidence, then it has been used to discriminate the models. The data uncertainty is in the range $\sigma\in(0.5,10)\%$ and two different sampling rates have been considered. Our results indicate that it is possible to discriminate $w=-1.02$ (or $w=-0.98$) model from the $\Lambda$CDM $(w=-1)$ with $\sigma=0.5\%$ uncertainty in expansion rate data. Finally, we perform a parameters inference in both the MCMC and Gaussian linear model, using currently available expansion rate data, and compare the results.

Inverse techniques are used to extract information about an exoplanet's atmosphere. These techniques are prone to biased results if the appropriate forward model is not used. One assumption used in a forward model is to assume that the radius of the planet is constant with wavelength, however a more realistic assumption is that the photospheric radius varies with each wavelength. We explore the bias induced when attempting to extract the molecular abundance from an emission spectrum which was generated with a variable radius. We find that for low gravity planets, the retrieval model is not able to fit the data if a constant radius model is used. We find that biased results are obtained when studying a typical hot Jupiter in the MIRI LRS wavelength range. Finally, we show that high gravity planets do not suffer a bias. We recommend that future spectral retrievals that interpret exoplanet emission spectra should take into account a variable radius.

GRO J1750-27, discovered during an outburst in 1995 with CGRO-BATSE, is one of the farthest known Be-X-ray binary systems. This relatively poorly studied system recently went into an outburst. From the measurements of the spin-up rate of this pulsar, the X-ray fluxes during the two previous outbursts in 2008 and 2014-2015, and using the standard theory of accretion torque on magnetized neutron stars, the magnetic field strength was estimated to be $2 \times 10^{12}$ G and $3.5-4.5 \times 10^{12}$ G respectively. The uncertainty in the distance causes large uncertainty in the estimated magnetic field. The source was observed during the latest outburst using the NuSTAR telescope during the rising phase of the outburst. We estimate the spin period of the source to be ~ 4.45s using which we produced energy-resolved pulse profiles between 3 and 79 keV. We find that the profiles appear to be double-peaked at low energies (<20 keV) while evolving into a single peak at higher energies (>20 keV). The broadband spectrum of this source was described by a power-law modified by an exponential cut off and we report the discovery of a cyclotron resonant scattering feature (CRSF) in this hard X-ray spectrum of this source at ~43 keV indicating a magnetic field strength of 3.7 $\times 10^{12}$ G. Our estimate of the magnetic field strength using the cyclotron line is consistent with the estimate made using the accretion torque model.

While Jupiter's massive gas envelope consists mainly of hydrogen and helium, the key to understanding Jupiter's formation and evolution lies in the distribution of the remaining (heavy) elements. Before the Juno mission, the lack of high-precision gravity harmonics precluded the use of statistical analyses in a robust determination of the heavy-elements distribution in Jupiter's envelope. In this paper, we assemble the most comprehensive and diverse collection of Jupiter interior models to date and use it to study the distribution of heavy elements in the planet's envelope. We apply a Bayesian statistical approach to our interior model calculations, reproducing the Juno gravitational and atmospheric measurements and constraints from the deep zonal flows. Our results show that the gravity constraints lead to a deep entropy of Jupiter corresponding to a 1 bar temperature 5-15 K higher than traditionally assumed. We also find that uncertainties in the equation of state are crucial when determining the amount of heavy elements in Jupiter's interior. Our models put an upper limit to the inner compact core of Jupiter of 7 Earth masses, independently on the structure model (with or without dilute core) and the equation of state considered. Furthermore, we robustly demonstrate that Jupiter's envelope is inhomogenous, with a heavy-element enrichment in the interior relative to the outer envelope. This implies that heavy element enrichment continued through the gas accretion phase, with important implications for the formation of giant planets in our solar system and beyond.

A model-independent, weighted semi-discrete, fast optimal transport algorithm to reconstruct the Lagrangian positions of proto-halos from their evolved Eulerian positions is presented. Tests with state-of-art cosmological simulations show that the positions of proto-halos are reconstructed accurately, without having to assume a background cosmology. The algorithm, which makes use of a mass estimate of the biased tracers, but is robust to errors in this estimation, recovers the shape and amplitude of the initial pair correlation function of the tracers, enabling sub-percent precision in the BAO distance scale that is not tied to a cosmological model. In principle, our algorithm also allows direct and independent determinations of the bias factor and the smearing scale, potentially providing new methods for breaking the degeneracy between the bias factor $b$ and $\sigma_8$.

High-contrast imaging in the visible and near-infrared (VIS/NIR) has revealed the presence of a plethora of substructures in circumstellar disks (CSDs). One of the most commonly observed substructures are concentric gaps that may hint at the presence of embedded forming planets. However, direct detections of them are rare, and thus ambiguity regarding the origin of most gap features remains. The aim of this study is to investigate the capabilities of high-contrast VIS/NIR imaging of directly detecting and characterizing low-mass giant planets in gaps. To this end, a grid of models of protoplanetary disks was generated. The models include a central T Tauri star surrounded by a face-on CSD harboring an accreting planet, which itself is surrounded by a circumplanetary disk (CPD) and carves a gap. We use the Monte Carlo radiative transfer code Mol3D to generate temperature distributions and synthetic observations. Based on these simulations, we measured the impact the planet and its CPD have on contrast curves and quantified the impact of the observing wavelength and of five key parameters on the determined signal strength. Then, we applied a detection criterion on our results and assess the capabilities of SPHERE/VLT of detecting the embedded planets. We find that a part of the investigated parameter space includes detectable planets, and we elaborate on the implications of a non-detection. Furthermore, we analyze the potential loss of valuable information by the use of a too small coronagraphic mask. However, we find this outcome to be very unlikely. Finally, within the VIS/NIR we identify for each of the investigated basic properties of the planets and the disks the most promising observing wavelengths that enable us to distinguish between different underlying parameter values. We find, that the detectability and the characterization often benefit from different observing wavelengths.

This paper focuses on the observations of sunspots made by Ru{\dj}er Bo\v{s}kovi\'c in 1777. We derived the expressions needed to calculate the elements of the Sun's rotation and period from observations. We used modern ephemeris data in the processing of the observation results. Obtained results are very similar to Bo\v{s}kovi\'c's original calculations. In addition to historical significance, they also provide scientifically valuable data on the Sun's differential rotation, which plays a significant role in generating and maintaining solar magnetic activity.

We report the first detection of the S(1) pure rotational line of ortho-H2 at 17.04 um in an asymptotic giant branch star, using observations of IRC+10216 with the Echelon-cross-echelle Spectrograph (EXES) mounted on the Stratospheric Observatory for Infrared Astronomy (SOFIA). This line, which was observed in a very high sensitivity spectrum (RMS noise ~0.04% of the continuum), was detected in the wing of a strong telluric line and displayed a P Cygni profile. The spectral ranges around the frequencies of the S(5) and S(7) ortho-H2 transitions were observed as well but no feature was detected in spectra with sensitivities of 0.12% and 0.09% regarding the continuum emission, respectively. We used a radiation transfer code to model these three lines and derived a mass-loss rate of 2.43(0.21)E-05 M_sun/yr without using the CO abundance. The comparison of this rate with previous estimates derived from CO observations suggests that the CO abundance relative to H2 is 6.7(1.4)E-04. From this quantity and previously reported molecular abundances, we estimate the O/H and C/H ratios to be 3.3(0.7)E-04 and >5.2(0.9)E-04, respectively. The C/O ratio is >1.5(0.4). The absence of the S(5) and S(7) lines of ortho-H2 in our observations can be explained by the opacity of hot dust within 5R* from the center of the star. We estimate the intensity of the S(0) and S(2) lines of para-H2 to be ~0.1% and 0.2% of the continuum, respectively, which are below the detection limit of EXES.

We consider the effective field theory formulation of torsional gravity in a cosmological framework to alter the background evolution. Then we use the latest $H_0$ measurement from the SH0ES Team as well as observational Hubble data from cosmic chronometer (CC) and radial baryon acoustic oscillations (BAO) and we reconstruct the $f(T)$ form in a model-independent way by applying Gaussian processes. Since the special square-root term does not affect the evolution at the background level, we finally result to a family of functions that can produce the background evolution required by the data. Finally, performing a fitting using polynomial functions we find an analytic expression that may describe the cosmological evolution in reat agreement with observations.

Heavy axion-like particles (ALPs), with masses up to a few 100 keV and coupled with photons can be efficiently produced in stellar plasmas, contributing to a significant energy-loss. This argument has been applied to helium burning stars in Globular Clusters (GCs) to obtain stringent bounds on the ALP-photon coupling $g_{a\gamma}$. However, for sufficiently large values of the ALP mass and coupling to photons, one should expect a significant fraction of ALPs to decay inside the star. These ALPs do not contribute to the energy loss but rather lead to an efficient energy transfer inside the star. We present a new ballistic recipe that covers both the energy-loss and energy-transfer regimes and we perform the first dedicated simulation of GC stars including the ALP energy transfer. This argument allows us to constrain ALPs with $m_a \lesssim 0.4$ MeV and $g_{a\gamma} \simeq 10^{-5}$ GeV$^{-1}$, probing a section of the ALP parameter space informally known as "cosmological triangle". This region is particularly interesting since it has been excluded only using standard cosmological arguments that can be evaded in nonstandard scenarios.

The type-II seesaw model is a possible candidate for simultaneously explaining non-vanishing neutrino masses and the observed baryon asymmetry of the Universe. In this work, we study in detail the pattern of phase transition and the gravitational wave production of this model. We find a strong first-order electroweak phase transition generically prefers positive Higgs portal couplings and a light triplet below $\sim550$ GeV. In addition, we find the gravitational wave yield generated during the phase transition would be at the edge of BBO sensitivity and could be further examined by Ultimate-DECIGO.

We consider a cosmological model consisting of two scalar fields defined in the hyperbolic plane known as hyperbolic inflation. For the background space we consider a homogeneous and isotropic spacetime with nonzero curvature. We study the asymptotic behaviour of the dynamics and we search for attractors in the expanding regime. We prove that the hyperobolic inflation period is an unstable solution and the flatness problem is solved when the universe goes through a second inflationary era or when exits the inflation the universe reaches a matter era.

For the minimal wino and Higgsino benchmark models we provide accurate energy spectra of high-energy photons from TeV scale dark-matter annihilation $\chi\chi\to \gamma+X$ by merging electroweak Sudakov resummation near maximal energy with the electroweak parton-shower PPPC4DM, and accounting for the Sommerfeld effect. Electroweak resummation significantly changes the shape of the photon-energy spectrum in the wide range $E_\gamma \sim (0.6\ldots 1)\, m_\chi$ and hence the form of the so-called "line-signal".

Self-gravitating bodies can have an arbitrarily complex shape, which implies a much richer multipolar structure than that of a black hole in General Relativity. With this motivation, we study the corrections to the dynamics of a binary system due to generic, nonaxisymmetric mass quadrupole moments to leading post-Newtonian (PN) order. Utilizing the method of osculating orbits and a multiple scale analysis, we find analytic solutions to the precession and orbital dynamics of a (generically eccentric) binary in terms of the dimensionless modulus parameters $\epsilon_{m}$, corresponding to axial $m=1$ and polar $m=2$ corrections from oblateness/prolateness. The solutions to the precession dynamics are exact for $0 \le \epsilon_{2} < 1$, and perturbative in $\epsilon_{1} \ll 1$. We further compute the leading order corrections to the gravitational wave amplitude and phase for a quasi-circular binary due to mass quadrupole effects. Making use of the stationary phase approximation and shifted uniform asymptotics (SUA), the corrections to the phase enter at relative 2PN order, while the amplitude modulations enter at -0.5PN order with a SUA amplitude correction at 5.5PN order, relative 2PN order to the leading order SUA correction. By investigating the dephasing due to generic quadrupole moments, we find that a phase difference $\gtrsim 0.1$~radians is achievable for $\epsilon_{m} \gtrsim 10^{-3}$, which suggests that constraints with current and future ground-based gravitational wave detectors are possible. Our results can be implemented in parameter estimation studies to constrain generic multipolar deformations of the Kerr geometry and of neutron stars.

Continuous gravitational waves are nearly monochromatic signals emitted by asymmetries in rotating neutron stars. These signals have not yet been detected. Deep all-sky searches for continuous gravitational waves from isolated neutron stars require significant computational expense. Deep searches for neutron stars in binary systems are even more expensive, but potentially these targets are more promising emitters, especially in the hundreds-Hz region, where ground-based gravitational wave detectors are most sensitive. We present here an all-sky search for continuous signals with frequency between 300 and 500 Hz, from neutron stars in binary systems with orbital period between 15 and 60 days, and projected semi-major axis between 10 and 40 light-seconds. This is the only binary search on Advanced-LIGO data that probes this frequency range. Compared to previous results, our search is over an order of magnitude more sensitive. We do not detect any signals, but our results exclude plausible and unexplored neutron star configurations, for example, neutron stars with relative deformations greater than 3e-6 within 1 kpc from Earth and r-mode emission at the amplitude level of ~ 1e-4 within the same distance.

We present several novel ionization-electron and scintillation-photon recording concepts in noble-liquid detectors, for future applications in particle and astroparticle physics and in other fields. These involve both single- and dual-phase detector configurations with combined electroluminescence and small charge multiplication in gas and liquid media.