Papers for Friday, Jan 13 2023

Aims. We want to find the distribution of initial conditions that best reproduces disc observations at the population level. Methods. We first ran a parameter study using a 1D model that includes the viscous evolution of a gas disc, dust, and pebbles, coupled with an emission model to compute the millimetre flux observable with ALMA. This was used to train a machine learning surrogate model that can compute the relevant quantity for comparison with observations in seconds. This surrogate model was used to perform parameter studies and synthetic disc populations. Results. Performing a parameter study, we find that internal photoevaporation leads to a lower dependency of disc lifetime on stellar mass than external photoevaporation. This dependence should be investigated in the future. Performing population synthesis, we find that under the combined losses of internal and external photoevaporation, discs are too short lived. Conclusions. To match observational constraints, future models of disc evolution need to include one or a combination of the following processes: infall of material to replenish the discs, shielding of the disc from internal photoevaporation due to magnetically driven disc winds, and extinction of external high-energy radiation. Nevertheless, disc properties in low-external-photoevaporation regions can be reproduced by having more massive and compact discs. Here, the optimum values of the $\alpha$ viscosity parameter lie between $3\times10^{-4}$ and $10^{-3}$ and with internal photoevaporation being the main mode of disc dispersal.

We present the JWST Resolved Stellar Populations Early Release Science (ERS) science program. We obtained 27.5 hours of NIRCam and NIRISS imaging of three targets in the Local Group (Milky Way globular cluster M92, ultra-faint dwarf galaxy Draco II, star-forming dwarf galaxy WLM), which span factors of $\sim10^5$ in luminosity, $\sim10^4$ in distance, and $\sim10^5$ in surface brightness. We describe the survey strategy, scientific and technical goals, implementation details, present select NIRCam color-magnitude diagrams (CMDs), and validate the NIRCam exposure time calculator (ETC). Our CMDs are among the deepest in existence for each class of target. They touch the theoretical hydrogen burning limit in M92 ($<0.08$ $M_{\odot}$; SNR $\sim5$ at $m_{F090W}\sim28.2$; $M_{F090W}\sim+13.6$), include the lowest-mass stars observed outside the Milky Way in Draco II (0.09 $M_{\odot}$; SNR $=10$ at $m_{F090W}\sim29$; $M_{F090W}\sim+12.1$), and reach $\sim1.5$ magnitudes below the oldest main sequence turnoff in WLM (SNR $=10$ at $m_{F090W}\sim29.5$; $M_{F090W}\sim+4.6$). The PARSEC stellar models provide a good qualitative match to the NIRCam CMDs, though are $\sim0.05$ mag too blue compared to M92 F090W$-$F150W data. The NIRCam ETC (v2.0) matches the SNRs based on photon noise from DOLPHOT stellar photometry in uncrowded fields, but the ETC may not be accurate in more crowded fields, similar to what is known for HST. We release beta versions of DOLPHOT NIRCam and NIRISS modules to the community. Results from this ERS program will establish JWST as the premier instrument for resolved stellar populations studies for decades to come.

We present photometric and spectroscopic observations of SN 2020bio, a double-peaked Type IIb supernova (SN) discovered within a day of explosion, primarily obtained by Las Cumbres Observatory and Swift. SN 2020bio displays a rapid and long-lasting initial decline throughout the first week of its light curve, similar to other well-studied Type IIb SNe. This early-time emission is thought to originate from the cooling of the extended outer envelope of the progenitor star that is shock-heated by the SN explosion. We compare SN 2020bio to a sample of other double-peaked Type IIb SNe to investigate its progenitor properties. Analytical model fits to the early-time emission give progenitor radius ($\approx$ 100--1500 $R_\odot$) and H-rich envelope mass ($\approx$ 0.01--0.5 $M_\odot$) estimates that are consistent with other Type IIb SNe. However, SN 2020bio displays several peculiarities, including: 1) weak H spectral features and narrow emission lines indicative of pre-existing circumstellar material; 2) an underluminous secondary light curve peak which implies a small amount of synthesized $^{56}$Ni ($M_{\text{Ni}}$ $\approx$ 0.02 $M_\odot$); and 3) low-luminosity nebular [O I] features. These observations are more consistent with a lower-mass progenitor (M$_{\text{ZAMS}} \approx$ 12 $M_\odot$) that was stripped of most of its H envelope before exploding. This study adds to the growing diversity in the observed properties of Type IIb SNe and their progenitors.

We present the first-ever collection of delta Scuti stars found over the entire area of the Small Magellanic Cloud (SMC). The sample consists of 2810 variables of which over 2600 objects belong to the SMC while the remaining stars are most likely members of the Milky Way's halo. The sample has been divided into 2733 singlemode and 77 multimode pulsators. We provide observational parameters (pulsation periods, mean magnitudes, amplitudes, Fourier coefficients) of all delta Sct stars and the long-term I- and V-band time-series photometric measurements collected during the fourth phase of the Optical Gravitational Lensing Experiment (OGLE-IV).

We investigate the stellar mass-dependence of the galaxy size-dark matter halo radius relation for low redshift galaxies using weak gravitational lensing measurements. Our sample consists of $\sim$38,000 galaxies more massive than $10^{8} {\rm M}_{\odot}h^{-2}$ and within $z<0.3$ drawn from the overlap of GAMA survey DR4 and HSC-SSP PDR2. We divide our sample into a number of stellar mass bins and measure stacked weak lensing signals. We model the signals using a conditional stellar mass function formalism to infer the stellar mass-halo mass relation. We fit a single S\'ersic model to HSC $i$-band images of our sample galaxies and obtain their three-dimensional half-light radii. We use these measurements to construct a median galaxy size-mass relation. We then combine the two relations to derive the relationship between galaxy size and halo radius. We confirm that the galaxy size-halo radius relation is roughly linear over two orders of magnitudes in stellar mass above $\sim 10^{9.35} {\rm M}_{\odot}h^{-2}$. Below this stellar mass, we see evidence of a downward departure from this linearity of up to 25 percent at stellar masses of $10^{8.7} {\rm M}_{\odot}h^{-2}$. If the stellar mass halo mass relation is extrapolated to lower mass scales, our galaxy size measurements imply these deviations can reach even as high as 50 percent at $10^{8} {\rm M}_{\odot}h^{-2}$. The existence of a such trend in dwarf galaxy sectors calls for either modification in models employing a constant fraction of halo angular momentum transferred to explain sizes of dwarfs or else points towards our lack of knowledge about the host dark matter haloes of such low-mass galaxies.

Magnetic white dwarfs with field strengths below 10 MG are easy to recognise since the Zeeman splitting of spectral lines appears proportional to the magnetic field strength. For fields $\geq 100$ MG, however, transition wavelengths become chaotic, requiring quantum-chemical predictions of wavelengths and oscillator strengths with a non-perturbative treatment of the magnetic field. While highly accurate calculations have previously been performed for hydrogen and helium, the variational techniques employed become computationally intractable for systems with more than three to four electrons. Modern computational techniques, such as finite-field coupled-cluster theory, allow the calculation of many-electron systems in arbitrarily strong magnetic fields. Because around 25 percent of white dwarfs have metal lines in their spectra, and some of those are also magnetic, the possibility arises for some metals to be observed in very strong magnetic fields, resulting in unrecognisable spectra. We have identified SDSSJ114333.48+661531.83 as a magnetic DZ white dwarf, with a spectrum exhibiting many unusually shaped lines at unknown wavelengths. Using atomic data calculated from computational finite-field coupled-cluster methods, we have identified some of these lines arising from Na, Mg, and Ca. Surprisingly, we find a relatively low field strength of 30 MG, where the large number of overlapping lines from different elements make the spectrum challenging to interpret at a much lower field strength than for DAs and DBs. Finally we model the field structure of SDSSJ1143+6615 finding the data are consistent with an offset dipole.

We report the detection of 15 GHz radio continuum emission associated with the classical Cepheid variable star delta Cephei based on observations with the Karl G. Jansky Very Large Array. Our results constitute the first probable detection of radio continuum emission from a classical Cepheid. We observed the star at pulsation phase phi~0.43 (corresponding to the phase of maximum radius and minimum temperature) during three pulsation cycles in late 2018 and detected statistically significant emission (>5 sigma) during one of the three epochs. The observed radio emission appears to be variable at a >~10% level on timescales of days to weeks. We also present an upper limit on the 10 GHz flux density at pulsation phase phi=0.31 from an observation in 2014. We discuss possible mechanisms that may produce the observed 15 GHz emission, but cannot make a conclusive identification from the present data. The emission does not appear to be consistent with originating from a close-in, late-type dwarf companion, although this scenario cannot yet be strictly excluded. Previous X-ray observations have shown that delta Cephei undergoes periodic increases in X-ray flux during pulsation phase phi~0.43. The lack of radio detection in two out of three observing epochs at phi~0.43 suggests that either the radio emission is not linked with a particular pulsation phase, or else that the strength of the generated radio emission in each pulsation cycle is variable.

Accreted stellar populations are comprised of the remnants of destroyed galaxies, and often dominate the `stellar haloes' of galaxies such as the Milky Way (MW). This ensemble of external contributors is a key indicator of the past assembly history of a galaxy. We introduce a novel statistical method that uses the unbinned metallicity distribution function (MDF) of a stellar population to estimate the mass spectrum of its progenitors. Our model makes use of the well-known mass-metallicity relation of galaxies and assumes Gaussian MDF distributions for individual progenitors: the overall MDF is thus a mixture of MDFs from smaller galaxies. We apply the method to the stellar halo of the MW, as well as the classical MW satellite galaxies. The stellar components of the satellite galaxies have relatively small sample sizes, but we do not find any evidence for accreted populations with L > L_host/100. We find that the MW stellar halo has N~1-3 massive progenitors (L > 10^8 L_Sun) within 10 kpc, and likely several hundred progenitors in total. We also test our method on simulations of MW-mass haloes, and find that our method is able to recover the true accreted population within a factor of two. Future datasets will provide MDFs with orders of magnitude more stars, and this method could be a powerful technique to quantify the accreted populations down to the ultra-faint dwarf mass-scale for both the MW and its satellites.

Prompt cusps are the densest quasi-equilibrium dark matter objects; one forms at the instant of collapse within every isolated peak of the initial cosmological density field. They have power-law density profiles, $\rho \propto r^{-1.5}$ with central phase-space density set by the primordial velocity dispersion of the dark matter. At late times they account for $\sim 1\%$ of the dark matter mass but for $>90\%$ of its annihilation luminosity in all but the densest regions, where they are tidally disrupted. Here we demonstrate that individual stellar encounters, rather than the mean galactic tide, are the dominant disruptors of prompt cusps within galaxies. Their cumulative effect is fully (though stochastically) characterised by an impulsive shock strength $B_* = 2\pi G\int\rho_*({\bf x}(t))\, \mathrm{d}t$ where $\rho_*$, the total mass density in stars, is integrated over a cusp's entire post-formation trajectory. Stellar encounters and mean tides have only a small effect on the halo annihilation luminosity seen by distant observers, but this is not true for the Galactic halo because of the Sun's position. For a 100 GeV WIMP, Earth-mass prompt cusps are predicted, and stellar encounters suppress their mean annihilation luminosity by a factor of two already at 20 kpc, so that their annihilation emission is predicted to appear almost uniform over the sky. The Galactic Center $\gamma$-ray Excess is thus unaffected by cusps. If it is indeed dark matter annihilation radiation, then prompt cusps in the outer Galactic halo and beyond must account for 20-80% of the observed isotropic $\gamma$-ray background in the 1 to 10 GeV range.

We report on the identification of very massive stars (VMS, mass $> 100$\,\msun) possibly segregated in the center of the young massive star cluster at $z$=2.37 hosted in the {\tt Sunburst} lensed galaxy. Such a result is based on two pieces of evidence: (1) the VLT/MUSE spectra of several multiple images of the same star cluster show key spectral signatures of VMS, like the \heii\ broad emission, \nivblue\ emission and \niv\ P-Cygni profile. In particular, \heii\ is broad ($\sim1610\pm300$ \kms) with an equivalent width of 3\AA\ and shows an asymmetric profile. Such features require an extremely young ($\sim2.5$ Myr) stellar population component with masses of the stars exceeding 100~\msun. Assuming a Salpeter IMF and BPASS models for normal massive stars, the observed spectral features require $\sim$400 VMS; (2) the same star cluster is detected at S/N~$\sim100$ in the LyC domain ($\lambda < 900$\AA). The LyC emission emerges from a region with a radius at least 2 times smaller than what is observed at 1700\AA~(independently from magnification) and is located in the center of the cluster. In absolute scales, after de-lensing, the effective radii are R$_{\tt eff}[{\tt LyC}]\sim4.7 \pm 1.5$ pc and R$_{\tt eff}[1700]= 7.8 \pm 1.4$ pc. The LyC radiation is mainly produced by hot and massive stars, implying that their spatial distribution (including VMS) is preferentially more confined in the central parts of the cluster. Approximately 400 VMS hosted by a cluster of $\sim 10^7$ \msun\ are producing $\sim$15\% of the escaping LyC photons, while the rest is produced from other massive early-type stars.

The particle-in-cell (PIC) method is successfully used to study magnetized plasmas. However, this requires large computational costs and limits simulations to short physical run-times and often to setups in less than three spatial dimensions. Traditionally, this is circumvented either via hybrid-PIC methods (adopting massless electrons) or via magneto-hydrodynamic-PIC methods (modelling the background plasma as a single charge-neutral magneto-hydrodynamical fluid). Because both methods preclude modelling important plasma-kinetic effects, we introduce a new fluid-PIC code that couples a fully explicit and charge-conservative multi-fluid solver to the PIC code SHARP through a current-coupling scheme and solve the full set of Maxwell's equations. This avoids simplifications typically adopted for Ohm's Law and enables us to fully resolve the electron temporal and spatial scales while retaining the versatility of initializing any number of ion, electron, or neutral species with arbitrary velocity distributions. The fluid solver includes closures emulating Landau damping so that we can account for this important kinetic process in our fluid species. Our fluid-PIC code is second-order accurate in space and time. The code is successfully validated against several test problems, including the stability and accuracy of shocks and the dispersion relation and damping rates of waves in unmagnetized and magnetized plasmas. It also matches growth rates and saturation levels of the gyro-scale and intermediate-scale instabilities driven by drifting charged particles in magnetized thermal background plasmas in comparison to linear theory and PIC simulations. This new fluid-SHARP code is specially designed for studying high-energy cosmic rays interacting with thermal plasmas over macroscopic timescales.

The formation of super-Earths, the most abundant planets in the Galaxy, remains elusive. These planets have masses that typically exceed that of the Earth by a factor of a few; appear to be predominantly rocky, although often surrounded by H/He atmospheres; and frequently occur in multiples. Moreover, planets that encircle the same star tend to have similar masses and radii, whereas those belonging to different systems exhibit remarkable overall diversity. Here, we advance a theoretical picture for rocky planet formation that satisfies the aforementioned constraints: building upon recent work - which demonstrates that planetesimals can form rapidly at discrete locations in the disk - we propose that super-Earths originate inside rings of silicate-rich planetesimals at approximately ~1 AU. Within the context of this picture, we show that planets grow primarily through pairwise collisions among rocky planetesimals, until they achieve terminal masses that are regulated by isolation and orbital migration. We quantify our model with numerical simulations and demonstrate that our synthetic planetary systems bear a close resemblance to compact, multi-resonant progenitors of the observed population of short-period extrasolar planets. Our results thus indicate that the absence of short-period super-Earths within the solar system can simply be attributed to the comparatively low mass of the primordial planetesimal ring within the protosolar nebula.

Aims. To introduce and develop a Regularized Maximum Likelihood (RML) algorithm designed to address the mathematically ill-posed problem of constructing differential emission measure profiles from a discrete set of EUV intensities in specified wavelength bands, specifically those observed by the Atmospheric Imaging Assembly (AIA) on the NASA Solar Dynamics Observatory. Methods. RML combines features of Maximum Likelihood and regularized approaches used by other authors. It is also guaranteed to produce a positive definite differential emission measure profile. Results. We evaluate and compare the effectiveness of the method against other published algorithms, using both simulated data generated from parametric differential emission profile forms, and AIA data from a solar eruptive event on 2010 November 3. Similarities and differences between the differential emission measure profiles and maps reconstructed by the various algorithms are discussed. Conclusions. The RML inversion method is mathematically rigorous, computationally efficient, and robust in the presence of data noise. As such it shows considerable promise for computing differential emission measure profiles from datasets of discrete spectral lines.

The orientation between a star's spin axis and a planet's orbital plane provides valuable information about the system's formation and dynamical history. For non-transiting planets at wide separations, true stellar obliquities are challenging to measure, but lower limits on spin-orbit orientations can be determined from the difference between the inclination of the star's rotational axis and the companion's orbital plane ($\Delta i$). We present results of a uniform analysis of rotation periods, stellar inclinations, and obliquities of cool stars (SpT $\gtrsim$ F5) hosting directly imaged planets and brown dwarf companions. As part of this effort, we have acquired new $v \sin i_*$ values for 22 host stars with the high-resolution Tull spectrograph at the Harlan J. Smith telescope. Altogether our sample contains 62 host stars with rotation periods, most of which are newly measured using light curves from the Transiting Exoplanet Survey Satellite. Among these, 53 stars have inclinations determined from projected rotational and equatorial velocities, and 21 stars predominantly hosting brown dwarfs have constraints on $\Delta i$. Eleven of these (52$^{+10}_{-11}$% of the sample) are likely misaligned, while the remaining ten host stars are consistent with spin-orbit alignment. As an ensemble, the minimum obliquity distribution between 10-250 AU is more consistent with a mixture of isotropic and aligned systems than either extreme scenario alone--pointing to direct cloud collapse, formation within disks bearing primordial alignments and misalignments, or architectures processed by dynamical evolution. This contrasts with stars hosting directly imaged planets, which show a preference for low obliquities. These results reinforce an emerging distinction between the orbits of long-period brown dwarfs and giant planets in terms of their stellar obliquities and orbital eccentricities.

The spin evolution of main sequence stars has long been of interest for basic stellar evolution, stellar aging, stellar activity, and consequent influence on companion planets. Observations of older than solar late-type main-sequence stars have been interpreted to imply that a change from a dipole-dominated magnetic field to one with more prominent higher multipoles might be necessary to account for the data. The spin-down models that lead to this inference are essentially tuned to the sun. Here we take a different approach which considers individual stars as fixed points rather than just the Sun. We use a time-dependent theoretical model to solve for the spin evolution of low-mass main-sequence stars that includes a Parker-type wind and a time-evolving magnetic field coupled to the spin. Because the wind is exponentially sensitive to the stellar mass over radius and the coronal base temperature, the use of each observed star as a separate fixed point is more appropriate and, in turn, produces a set of solution curves that produces a solution envelope rather than a simple line. This envelope of solution curves, unlike a single line fit, is consistent with the data and does not unambiguously require a modal transition in the magnetic field to explain it. Also, the theoretical envelope does somewhat better track the older star data when thermal conduction is a more dominant player in the corona.

Preservation of the chemical and structural integrity of samples that will be brought back from Mars is paramount to achieving the scientific objectives of MSR. Given our knowledge of the nature of the samples retrieved at Jezero by Perseverance, at least two options need to be tested for opening the sample tubes: (1) One or two radial cuts at the end of the tube to slide the sample out. (2) Two radial cuts at the ends of the tube and two longitudinal cuts to lift the upper half of the tube and access the sample. Strategy 1 will likely minimize contamination but incurs the risk of affecting the physical integrity of weakly consolidated samples. Strategy 2 will be optimal for preserving the physical integrity of the samples but increases the risk of contamination and mishandling of the sample as more manipulations and additional equipment will be needed. A flexible approach to opening the sample tubes is therefore required, and several options need to be available, depending on the nature of the rock samples returned. Both opening strategies 1 and 2 may need to be available when the samples are returned to handle different sample types (e.g., loosely bound sediments vs. indurated magmatic rocks). This question should be revisited after engineering tests are performed on analogue samples. The MSR sample tubes will have to be opened under stringent BSL4 conditions and this aspect needs to be integrated into the planning.

We present an overview and data release of the spectral line component of the SMA Large Program, \textit{CMZoom}. \textit{CMZoom} observed $^{12}$CO(2-1), $^{13}$CO(2-1) and C$^{18}$O(2-1), three transitions of H$_{2}$CO, several transitions of CH$_{3}$OH, two transitions of OCS and single transitions of SiO and SO, within gas above a column density of N(H$_2$)$\ge 10^{23}$\,cm$^{-2}$ in the Central Molecular Zone (CMZ; inner few hundred pc of the Galaxy). We extract spectra from all compact 1.3\,mm \emph{CMZoom} continuum sources and fit line profiles to the spectra. We use the fit results from the H$_{2}$CO 3(0,3)-2(0,2) transition to determine the source kinematic properties. We find $\sim 90$\% of the total mass of \emph{CMZoom} sources have reliable kinematics. Only four compact continuum sources are formally self-gravitating. The remainder are consistent with being in hydrostatic equilibrium assuming that they are confined by the high external pressure in the CMZ. Based on the mass and density of virially bound sources, and assuming star formation occurs within one free-fall time with a star formation efficiency of $10\% - 75\%$, we place a lower limit on the future embedded star-formation rate of $0.008 - 0.06$\,M$_{\odot}$\,yr$^{-1}$. We find only two convincing proto-stellar outflows, ruling out a previously undetected population of very massive, actively accreting YSOs with strong outflows. Finally, despite having sufficient sensitivity and resolution to detect high-velocity compact clouds (HVCCs), which have been claimed as evidence for intermediate mass black holes interacting with molecular gas clouds, we find no such objects across the large survey area.

We fit a new vertically extended corona model to previously measured reverberation time lags observed by \emph{XMM-Newton} in two extremely variable Narrow Line Seyfert 1 Active Galactic Nuclei (AGN), 1H~0707-495 and IRAS~13224-3809, in a variety of similarly observed flux groups and explore the model in all observations over a 16 year period. The model employs two X-ray sources located along the black hole rotational axis at height, $h_1$ and $h_2$ respectively. These sources have their associated photon indices $\Gamma_1$ and $\Gamma_2$ which respond to fluctuations in the disc with a maximum response duration of $t_\text{max}$ and a propagation delay between the response of the two of $t_\text{shift}$. We find that for 1H 0707-495, $h_2$ is significantly correlated with $\Gamma_1$ and anti-correlated with ionisation $\xi$. Whilst the 1H 0707-495 corona extends upwards, the emission appears softer and the disc is less ionised. We find similarities in IRAS 13224-3809, but significant anti-correlation between $\Gamma_2$ and both $t_\text{max}$ and $t_\text{shift}$. This suggests that when the IRAS 13224-3809 corona becomes softer while extending vertically upwards, the overall corona response occurs faster. This may also suggest that the inner disc also becomes more active. In addition, $\Gamma_1$ and $\Gamma_2$ are extreme, relatively less variable, but more separate in IRAS 13224-3809 than in 1H 0707-495. This suggests that the IRAS 13224-3809 corona may be more patchy in the sense that it has two more clear distinct spectral zones of $\Gamma_1$ and $\Gamma_2$ (possibly relating to two distinct zones of coronal temperature) when compared to 1H 0707-495.

Astrometry at centimeter wavelengths using Very Long Baseline Interferometry is approaching accuracies of ~1 uas for the angle between a target and a calibrator source separated by <1 degree on the sky. The BeSSeL Survey and the Japanese VERA project are using this to map the spiral structure of the Milky Way by measuring trigonometric parallaxes of hundreds of maser sources associated with massive, young stars. This paper outlines how micro-arcsecond astrometry is done, including details regarding the scheduling of observations, calibration of data, and measuring positions.

The James Webb Space Telescope (JWST) observed a section of the star forming region NGC 3324 during its Early Release Observations. We make use of the Probabilistic Random Forest machine learning model to identify YSOs within the field of view. We build a matched catalog from photometry data products available on the Mikulski Space Telescope Archive and retrieve 8632 objects, of which Spitzer previously detected 458. We use previously classified data from Spitzer to train on a sample of the Webb data. We retrieve a total of 72 YSO candidates within the data field, 52 of which are only visible with JWST.

In 1181 AD, Chinese and Japanese observers reported a bright `Guest Star' in the constellation Chuanshe, unmoving and visible for 185 days. In 2013, D. Patchick discovered a unique nebula surrounding a unique star, with two groups attributing this structure, named `Pa 30', to be the supernova remnant of SN 1181, as a sub-subclass of supernova, the low-luminosity Type Iax. Here, I provide a wide range of new observational evidence: First, a detailed analysis of the original Chinese and Japanese reports places the `Guest Star' of 1181 into a small region with the only interesting source being Pa 30. Second, the ancient records confidently place the peak magnitude as 0.0 > V_peak > -1.4, and hence peak absolute magnitude of -14.5 > M_V,peak > -16.0 mag. Third, the Pa 30 central star is fading from B=14.90 in 1917, to B=16.20 in 1950, to B=16.58 in 2022. Fourth, recent light curves show typical variability with full-amplitude of 0.24 mag on time-scales of one day and longer, critically with no coherent modulations for periods from 0.00046--10 days to strict limits. Fifth, the spectral energy distribution from the far-infrared to the ultraviolet is a nearly perfect power-law with F_nu proportional to nu^(0.99 +- 0.07), observed luminosity of 128 +- 24 L_Sun, and absolute magnitude M_V = +1.07. I collect my new evidences with literature results to make a confident case to connect the Oriental observations to a supernova, then to Pa 30, then to a low-luminosity Type Iax SN, then to the only possible explosion mechanism as a merger between CO and ONe white dwarfs.

A newly recognized young Galactic SN remnant, Pa 30 (G123.1+4.6), centered on a hot central star with a ~16,000 km/s wind velocity has recently been proposed to be the result of a double-degenerate merger leading to a SN Iax event associated with the guest star of 1181 CE. Here we present deep optical [S II] 6716,6731 images of Pa 30 which reveal an extraordinary and highly structured nebula 170" in diameter with dozens of long (5" - 20") radially aligned filaments with a convergence point near the hot central star. Optical spectra of filaments indicate a peak expansion velocity ~1100 km/s with electron densities of 100 to 700 cm^-3, and a thick shell-like structure resembling its appearance in 22 micron WISE images. No H-alpha emission was seen (6716/H-alpha >8), with the only other line emission detected being faint [Ar III] 7136 suggesting a S, Ar-rich but H-poor remnant. The nebula's angular size, estimated 2.3 kpc distance, and 1100 km/s expansion velocity are consistent with an explosion date around 1181 CE. The remnant's unusual appearance may be due to the photoionization of wind-driven ejecta due to clump-wind interactions caused by the central star's high-luminosity wind.

We present a new optical transmission spectrum of the hot Jupiter HAT-P-32Ab acquired with the Carnegie Observatories Spectrograph and Multiobject Imaging Camera (COSMIC) on the Palomar 200 inch Hale Telescope (P200). The P200/COSMIC transmission spectrum, covering a wavelength range of 3990--9390 \AA, is composed of 25 spectrophotometric bins with widths ranging from 200 to 400 \AA and consistent with previous transit measurements obtained in the common wavelength range. We derive a combined optical transmission spectrum based on measurements from five independent instruments, which, along with the 1.1--1.7 $\mu$m spectrum acquired by the Hubble Space Telescope and two Spitzer measurements, exhibits an enhanced scattering slope blueward of a relatively flat optical continuum, a water absorption feature at 1.4 $\mu$m, and a carbon dioxide feature at 4.4 $\mu$m. We perform Bayesian spectral retrieval analyses on the 0.3--5.1 $\mu$m transmission spectrum and find that it can be well explained by a two-limb approximation of $134^{+45}_{-33}\times$ solar metallicity, with a strongly hazy morning limb of $1134^{+232}_{-194}$ K and a haze-free evening limb of $1516^{+33}_{-44}$~K. This makes HAT-P-32Ab a promising target for James Webb Space Telescope to look for asymmetric signatures directly in the light curves.

In this paper, we report the discovery of an isolated, peculiar compact cloud with a steep velocity gradient at $2\farcm 6$ northwest of Sgr A*. This ``Tadpole'' molecular cloud is unique owing to its characteristic head-tail structure in the position-velocity space. By tracing the CO {\it J}=3--2 intensity peak in each velocity channel, we noticed that the kinematics of the Tadpole can be well reproduced by a Keplerian motion around a point-like object with a mass of $1\!\times\! 10^{5}\,M_{\odot}$. Changes in line intensity ratios along the orbit are consistent with the Keplerian orbit model. The spatial compactness of the Tadpole and absence of bright counterparts in other wavelengths indicate that the object could be an intermediate-mass black hole.

To investigate the dynamical properties of globular clusters, the surface brightness and kinematic data were collected and fitted to a family of lowered isothermal models called LIMEPY models. For 18 studied globular clusters, the amounts of concentration, truncation, and anisotropy were determined. In addition, the cluster mass, half-mass radius, distance, and mass-to-light ratio were also obtained. In general, LIMEPY models could describe these clusters well. Among these 18 clusters, NGC 5139, NGC 6388, and NGC 7078 were claimed to be candidates to host intermediate-mass black holes in literature. The models could not appropriately fit the central proper-motion velocity dispersion of NGC 5139 and the slope of proper-motion velocity-dispersion profile of NGC 6388. Thus, more dedicated models with intermediate-mass black holes or a group of stellar-mass black holes at cluster centers may need to be considered. Considering NGC 7078, our model with some degree of anisotropy can fit the data. Finally, the strong concentration-truncation anti-correlation and truncation-semimajor-axis correlation were revealed, which could be the observational imprint of the dynamical evolution of globular clusters.

We report on observations of the extended environment of the bright Be star $\gamma$-Cas performed using intensity interferometry measurements within its H$\alpha$ emission line. These observations were performed using a modified version of the I2C intensity interferometry instrument installed onto the 1.54 meter M\'{e}O optical metrology telescope and a portable 1-meter telescope (T1M). In order to better constrain the extent of the H$\alpha$ envelope, observations were performed for two different positions of the T1M telescope, corresponding to an intermediate and long baselines in which the extended region was partially and fully resolved. We find that the observed data are consistent with past interferometric observations of $\gamma$-Cas. These observations demonstrate the capability to equip optical telescopes of different optical designs with intensity interferometry capabilities and illustrate the potential to scale a similar system onto many additional telescopes.

Aims. Understanding the activity is vital for deciphering the structure, formation, and evolution of comets. We investigate models of cometary activity by comparing them to the dynamics of 67P/Churyumov-Gerasimenko. Methods. We matched simple thermal models of water activity to the combined Rosetta datasets by fitting to the total outgassing rate and four components of the outgassing induced non-gravitational force and torque, with a final manual adjustment of the model parameters to additionally match the other two torque components. We parametrised the thermal model in terms of a distribution of relative activity over the surface of the comet, and attempted to link this to different terrain types. We also tested a more advanced thermal model based on a pebble structure. Results. We confirm a hemispherical dichotomy and non-linear water outgassing response to insolation. The southern hemisphere of the comet and consolidated terrain show enhanced activity relative to the northern hemisphere and dust-covered, unconsolidated terrain types, especially at perihelion. We further find that the non-gravitational torque is especially sensitive to the activity distribution, and to fit the pole-axis orientation in particular, activity must be concentrated (in excess of the already high activity in the southern hemisphere and consolidated terrain) around the south pole and on the body and neck of the comet over its head. This is the case for both the simple thermal model and the pebble-based model. Overall, our results show that water activity cannot be matched by a simple model of sublimating surface ice driven by the insolation alone, regardless of the surface distribution, and that both local spatial and temporal variations are needed to fit the data.

We present a method to analyze images of the coma of 67P/Churyumov-Gerasimenko obtained using OSIRIS, the main imaging system onboard \textit{Rosetta}, where dust aggregates can be seen as bright tracks because of their relative velocity with respect to the spacecraft. We applied this method to 105 images taken in 2015 July, 2015 December and 2016 January, identifying more than 20000 individual objects. We performed a photometric analysis of them, finding their phase function. This phase function follows the same trend as the one found for the nucleus, consistent with the detected particles having a size larger than $\sim 1$ mm. Additionally, the phase function becomes shallower for increasing heliocentric distances, indicating a decrease in the mean agglomerate size. In order to characterize the agglomerates observed in the image, we developed a simplified model for their ejection and dynamics in the coma, and generated synthetic images based on it. We solved the inverse problem by finding the simulation parameters that give the best fit between synthetic and real images. In doing so, we were able to obtain a mean agglomerate size $\sim$ dm and initial speed $\simeq$ 1 m s$^{-1}$. Both show a decrease with increasing heliocentric distance, sign of the reduction in activity. Also, the sizes obtained by the comparison are not compatible with ejection caused by water activity, so other sources have to be invoked, mainly CO$_2$.

We announce the discovery of pulsations in HD 23642, the only bright eclipsing system in the Pleiades, based on light curves from the Transiting Exoplanet Survey Satellite (TESS). We measure 46 pulsation frequencies and attribute them to delta Scuti pulsations in the secondary component. We find four l=1 doublets, three of which have frequency splittings consistent with the rotation rate of the star. The dipole mode amplitude ratios are consistent with a high stellar inclination angle and the stellar rotation period agrees with the orbital period. Together, these suggest that the spin axis of the secondary is aligned with the orbital axis. We also determine precise effective temperatures and a spectroscopic light ratio, and use the latter to determine the physical properties of the system alongside the TESS data and published radial velocities. We measure a distance to the system in agreement with the Gaia parallax, and an age of 170 +/- 20 Myr based on a comparison to theoretical stellar evolutionary models.

To date, various formation channels of merging events have been heavily explored with the detection of nearly 100 double black hole (BH) merger events reported by the LIGO-Virgo-KAGRA (LVK) Collaboration. We here systematically investigate an alternative formation scenario, i.e., binary BHs (BBHs) formed through double helium stars (hereafter double-core evolution channel). In this scenario, the two helium stars (He-rich stars) could be the outcome of the classical isolated binary evolution scenario involving with and without common-envelope phase (i.e., CE channel and stable mass transfer channel), or alternatively of massive close binaries evolving chemically homogeneously (i.e., CHE channel). We perform detailed stellar structure and binary evolution calculations that take into account internal differential rotation and mass loss of He-rich stars, as well as tidal interactions in binaries. For double He-rich stars with equal masses in binaries, we find that tides start to be at work on the Zero Age Helium Main Sequence (ZAHeMS: the time when a He-rich star starts to burn helium in the core, which is analogous to ZAMS for core hydrogen burning) for initial orbital periods not longer than 1.0 day, depending on the initial metallicities. Besides the stellar mass loss rate and tidal interactions in binaries, we find that the role of the angular momentum transport efficiency in determining the resulting BH spins, becomes stronger when considering BH progenitors originated from a higher metal-metallicity environment. We highlight that double-core evolution scenario does not always produce fast-spinning BBHs and compare the properties of the BBHs reported from the LVK with our modeling.

In the present article, we discuss a numerical method of solution for the so-called "full non-LTE" radiation transfer problem, basic formalism of which was revisited by Paletou & Peymirat (2021; see also Oxenius 1986). More specifically, usual numerical iterative methods for non-LTE radiation transfer are coupled with the above-mentioned formalism. New numerical additions are explained in detail. We benchmark the whole process with the standard non-LTE transfer problem for a two-level atom with Hummer's (1962, 1969) $R_{\rm I-A}$ partial frequency redistribution function. We finally display new quantities such as the spatial distribution of the velocity distribution function of excited atoms, that can only be accessed to by adopting this more general frame for non-LTE radiation transfer.

In a rotation-dominated magnetosphere, there is a region where closed field lines rotate around the planet, and also a region where the open field lines stretch away from the planet, forming the lobes of the magnetotail. This paper shows that there could be a third, significantly different region, where the closed field lines form twisted vortex structures anchored in the magnetotail. Such patterns form when there are significant plasma sources inside the magnetosphere and the time scale of the plasmoid formation process is substantially larger than the planetary rotation period. In the presence of vortices, the Dungey and Vasyliunas cycles act differently. The Dungey flow does not penetrate the central region of the polar cap. Tail reconnection events are rare, thus leaving the plasma time enough to participate in the essentially 3-dimensional vortex-forming plasma motion. The above conditions are fulfilled for Saturn. We discovered vortex-like patterns in the plasma and magnetic field data measured by the Cassini spacecraft in the nightside magnetosphere of Saturn. The plasma whirling around in these vortices never reaches the dayside, instead, it performs a retrograde motion in the high latitude regions of the magnetotail. Low-energy plasma data suggest that the observed patterns correspond to the closed field line vortices.

This study is aimed at investigating the dynamic evolution of the orbits of stellar globular clusters (GCs). To integrate the orbits backward in time, the authors use models of the time-varying potentials derived from cosmological simulations, which are closest to the potential of the Galaxy. This allows for estimating the probability of close passages (collisions) of GCs with respect to each other and the Galactic center (GalC) in the Galaxy undergoing dynamic changes in the past. To reproduce the dynamics of the Galaxy in time, five potentials selected from the IllustrisTNG-100 large-scale cosmological database, which are similar in their characteristics to the current physical parameters of the Milky Way, are used. With these time-varying potentials, we have reproduced the orbital trajectories of 143 GCs 10 Gyr back in time using our original phi-GPU N-body code. Each GC was treated as a single physical particle with the assigned position and velocity of the cluster center from the Gaia DR2 observations. For each of the potentials, 1000 initial conditions were generated with randomized initial velocities of GCs within the errors of the observational data. In this study, we consider close passages to be passages with a relative distance of less than 100 pc and a relative speed of less than 250 km/s. To select clusters that pass at close distances from the GalC, the following criterion is applied based only on the relative distance: it must be less than 100 pc. Applying the above criteria, the authors obtained statistically significant rates of close passages of GCs with respect to each other and to the GalC. It has been determined that GCs during their evolution have approximately 10 intersecting trajectories with each other on the average and approximately 3 to 4 close passages near the GalC in 1 Gyr at a distance of 50 pc for each of the chosen potentials.

Dust grain dynamics in molecular clouds is regulated by its interplay with supersonic turbulent gas motions. The conditions under which dust grains decouple from the dynamics of gas remain poorly constrained. We first aim to investigate the critical dust grain size for dynamical decoupling, using both analytical predictions and numerical experiments. Second, we aim to set the range of validity of two fundamentally different numerical implementations for the evolution of dust and gas mixtures in turbulent molecular clouds. We carried out a suite of numerical experiments using two different schemes. First, we used a monofluid formalism in the terminal velocity approximation (TVA) on a Eulerian grid. Second, we used a two-fluid scheme, in which the dust dynamics is handled with Lagrangian super-particles, and the gas dynamics on a Eulerian grid. The monofluid results are in good agreement with the theoretical critical size for decoupling. We report dust dynamics decoupling for Stokes number St>0.1, that is, dust grains of $s>4~\mu$m in size. We find that the TVA is well suited for grain sizes of 10 $\mu$m in molecular clouds, in particular in the densest regions. However, the maximum dust enrichment measured in the low-density material where St>1 is questionable. In the Lagrangian dust experiments, we show that the results are affected by the numerics for all dust grain sizes. At St<<1, the dust dynamics is largely affected by artificial trapping in the high-density regions, leading to spurious variations of the dust concentration. At St>1, the maximum dust enrichment is regulated by the grid resolution used for the gas dynamics. The results of previous similar numerical work should therefore be revisited with respect to the limitations we highlight in this study. Dust enrichment of submicron dust grains is unlikely to occur in the densest parts of molecular clouds.

The Kepler mission revealed a plethora of stellar variability in the light curves of many stars, some associated with magnetic activity or stellar oscillations. In this work, we analyse the periodic signal in 162 intermediate-mass stars, interpreted as Rossby modes and rotational modulation - the so-called \textit{hump \& spike} feature. We investigate whether the rotational modulation (\textit{spike}) is due to stellar spots caused by magnetic fields or due to Overstable Convective (OsC) modes resonantly exciting g~modes, with frequencies corresponding to the convective core rotation rate. Assuming that the spikes are created by magnetic spots at the stellar surface, we recover the amplitudes of the magnetic fields, which are in good agreement with theoretical predictions. Our data show a clear anti-correlation between the spike amplitudes and stellar mass and possibly a correlation with stellar age, consistent with the dynamo-generated magnetic fields theory in (sub)-surface convective layers. Investigating the harmonic behaviour, we find that for 125 stars neither of the two possible explanations can be excluded. While our results suggest that the dynamo-generated magnetic field scenario is more likely to explain the \textit{spike} feature, we assess further work is needed to distinguish between the two scenarios. One method for ruling out one of the two explanations is to directly observe magnetic fields in \textit{hump \& spike} stars. Another would be to impose additional constraints through detailed modelling of our stars, regarding the rotation requirement in the OsC mode scenario or the presence of a convective-core (stellar age).

PSR J2140$-$2311B is a 13-ms pulsar discovered in 2001 in a 7.8-hour Green Bank Telescope (GBT) observation of the core-collapsed globular cluster M30 and predicted to be in a highly eccentric binary orbit. This pulsar has eluded detection since then, therefore its precise orbital parameters have remained a mystery until now. In this work, we present the confirmation of this pulsar using observations taken with the UHF receivers of the MeerKAT telescope as part of the TRAPUM Large Survey Project. Taking advantage of the beamforming capability of our backends, we have localized it, placing it $1.2(1)^\prime$ from the cluster centre. Our observations have enabled the determination of its orbit: it is highly eccentric ($e = 0.879$) with an orbital period of $6.2$ days. We also measured the rate of periastron advance, $\dot{\omega} = 0.078 \pm 0.002\, \rm deg \, yr^{-1}$. Assuming that this effect is fully relativistic, general relativity provides an estimate of the total mass of the system, $M_{\rm TOT} = 2.53 \pm 0.08$ M$_{\odot}$, consistent with the lightest double neutron star systems known. Combining this with the mass function of the system gives the pulsar and companion masses of $m_p < 1.43 \, \rm M_{\odot}$ and $m_c > 1.10 \, \rm M_{\odot}$ respectively. The massive, undetected companion could either be a massive WD or a NS. M30B likely formed as a result of a secondary exchange encounter. Future timing observations will allow the determination of a phase-coherent timing solution, vastly improving our uncertainty in $\dot{\omega}$ and likely enabling the detection of additional relativistic effects which will determine $m_p$ and $m_c$.

Context: The first tentative detection of a magnetic field on the Hot Jupiter type exoplanet $\tau$ Bo\"otis b was recently reported by Turner et al. (2021). The magnetic field was inferred from observations of circularly-polarized radio emission obtained with the LOFAR telescopes. The radio emission is possibly a consequence of the interaction of the surrounding stellar wind with the planet's magnetic field. Methods: We perform magnetohydrodynamic simulations of the space environment around $\tau$ Bo\"otis b and its interaction with the stellar wind using the PLUTO code. We study the magnetospheric energy fluxes and effects of different magnetic field orientations in order to understand the physical processes which cause energy fluxes leading to the observed radio emission given the proposed magnetic field strength in Turner et al. (2021). Furthermore we study the effect of stellar wind density and pressure on magnetospheric energy fluxes given the uncertainty of extrasolar stellar wind predictions. Results: We find in our simulations that the interaction is most likely super--Alfv\'enic and energy fluxes generated by the stellar wind--planet interaction are consistent with the observed radio powers. Magnetospheric Poynting fluxes are of the order of 1-8 $\times10^{18}$ W for open, semi-open and closed magnetospheres. The Poynting fluxes are energetically consistent with the radio powers in Turner et al. (2021) for a magnetospheric Poynting flux-to-radio efficiency $> 10^{-3}$ when the magnetic fields of the planet and star are aligned. In case of lower efficiency factors the magnetospheric emission scenario is according to the parameter space modeled in this study not powerful enough. In case of a magnetic polarity reversal of the host star towards an anti-aligned field configuration, expected radio powers in the magnetospheric emission scenario fall below the observable threshold.

We report new, extremely precise, photopolarimetry of the rapidly-rotating A0 main-sequence star $\zeta$ Aql, covering the wavelength range $\sim$400--900nm, which reveals a rotationally-induced signal. We model the polarimetry, together with the flux distribution and line profiles, in the framework of Roche geometry with $\omega$-model gravity darkening, to establish the stellar parameters. An additional constraint is provided by TESS photometry, which shows variability with a period, $P_{\rm phot}$, of 11.1 hr. Modelling based on solid-body surface rotation gives rotation periods, $P_{\rm rot}$, that are in only marginal agreement with this value. We compute new ESTER stellar-structure models to predict horizontal surface velocity fields, which depart from solid-body rotation at only the $\sim$2% level (consistent with a reasonably strong empirical upper limit on differential rotation derived from the line-profile analysis). These models bring the equatorial rotation period, $P_{\rm rot,e}$, into agreement with $P_{\rm phot}$, without requiring any 'fine tuning' (for the Gaia parallax). We confirm that surface abundances are significantly subsolar ($\mbox{[M/H]} \simeq -0.5$). The star's basic parameters are established with reasonably good precision: ${M = 2.53\pm0.16\,\mbox{M}_\odot}$, $\log{L/\mbox{L}_\odot} = 1.82\pm0.02$, $R_{\rm p} = 2.21\pm 0.02\,\mbox{R}_\odot$, $T_{\rm eff} = 9693 \pm 50~\mbox{K}$, $i = 85{^{+5}_{-7}}^\circ$, and $\omega/\omega_{\rm c} = 0.95\pm0.02$. Comparison with single-star, solar-abundance stellar-evolution models incorporating rotational effects shows excellent agreement (but somewhat poorer agreement for models at $\mbox{[M/H]} \simeq -0.4$).

Context. Our Galaxy is composed of different stellar populations with varying chemical abundances, which are thought to imprint the composition of planet building blocks (PBBs). As such, the properties of stars should affect the properties of planets and small bodies formed in their systems. In this context, high-resolution spectroscopic surveys open a window into the chemical links between and their host stars. Aims. We aim to determine the PBB composition trends for various stellar populations across the Galaxy by comparing the two large spectroscopic surveys APOGEE and GALAH. We assess the reliability of the PBB composition as determined with these surveys with a propagation error study. Methods. Stellar spectroscopic abundances from the large surveys GALAH-DR3 and APOGEE-DR17 were used as input with a stoichiometric condensation model. We classified stars into different Galactic components and we quantified the PBB composition trends as a function of [Fe/H]. We also analysed the distribution composition patterns in the [$\alpha$/Fe]-[Fe/H] diagram. Results. Our propagation error study suggests that the overall trends with [Fe/H] and [$\alpha$/Fe] are robust, which is supported by the double study of both APOGEE and GALAH. We therefore confirm the existence of a bimodal PBB composition separating the thin disc stars from the thick disc stars. Furthermore, we confirm that the stoichiometric water PBB content is anti-correlated with [Fe/H]. Conclusions. Our results imply that metal-poor stars both in the thin and thick disks are suitable hosts for water-rich PBBs and for ice-rich small bodies. However, for metal-poor stars ([Fe/H]<0), the PBBs around thick disc stars should have a higher water content than that around thin disc stars because of the $\alpha$-contentdependence of the water mass fraction.

The leading theoretical paradigm for the Sun's magnetic cycle is an $\alpha\omega$-dynamo process, in which a combination of differential rotation and turbulent, helical flows produces a large-scale magnetic field that reverses every 11 years. Most $\alpha\omega$ solar dynamo models rely on differential rotation in the solar tachocline to generate a strong toroidal field. The most problematic part of such models is then the production of the large-scale poloidal field, via a process known as the $\alpha$-effect. Whilst this is usually attributed to small-scale convective motions under the influence of rotation, the efficiency of this regenerative process has been called into question by some numerical simulations. Motivated by likely conditions within the tachocline, the aim of this paper is to investigate an alternative mechanism for the poloidal field regeneration, namely the magnetic buoyancy instability in a shear-generated, rotating magnetic layer. We use a local, fully compressible model in which an imposed vertical shear winds up an initially vertical magnetic field. The field ultimately becomes buoyantly unstable, and we measure the resulting mean electromotive force (EMF). For sufficiently rapid rotation, we find that a significant component of the mean EMF is aligned with the direction of the mean magnetic field, which is the characteristic feature of the classical $\alpha\omega$-dynamo model. Our results therefore suggest that magnetic buoyancy could contribute directly to the generation of large-scale poloidal field in the Sun.

Near-Earth Asteroids (NEAs) are a key testbed for investigations into planet formation, asteroid dynamics, and planetary defense initiatives. These studies rely on understanding NEA sizes, albedo distributions, and regolith properties. Simple thermal models are a commonly used method for determining these properties, however they have inherent limitations due to the simplifying assumptions they make about asteroid shapes and properties. With the recent collapse of the Arecibo Telescope and a decrease of direct size measurements, as well as future facilities such as LSST and NEO Surveyor coming online soon, these models will play an increasingly important role in our knowledge of the NEA population. Therefore, it is key to understand the limits of these models. In this work we constrain the limitations of simple thermal models by comparing model results to more complex thermophysical models, radar data, and other existing analyses. Furthermore, we present a method for placing tighter constraints on inferred NEA properties using simple thermal models. These comparisons and constraints are explored using the NEA (285263) 1998 QE2 as a case study. We analyze QE2 with a simple thermal model and data from both the NASA IRTF SpeX instrument and NEOWISE mission. We determine an albedo between 0.05 and 0.10 and thermal inertia between 0 and 425 J m$^{-2}$ s$^{-1/2}$ K$^{-1}$. We find that overall the simple thermal model is able to well constrain the properties of QE2, however we find that model uncertainties can be influenced by topography, viewing geometry, and the wavelength range of data used.

Observations of protoplanetary discs have shown that forming exoplanets leave characteristic imprints on the gas and dust of the disc. In the gas, these forming exoplanets cause deviations from Keplerian motion, which can be detected through molecular line observations. Our previous work has shown that machine learning can correctly determine if a planet is present in these discs. Using our machine learning models, we identify strong, localized non-Keplerian motion within the disc HD 142666. Subsequent hydrodynamics simulations of a system with a 5 Jupiter-mass planet at 75 au recreates the kinematic structure. By currently established standards in the field, we conclude that HD 142666 hosts a planet. This work represents a first step towards using machine learning to identify previously overlooked non-Keplerian features in protoplanetary discs.

We present cosmological-scale 3-dimensional (3D) neutral hydrogen ({\sc Hi}) tomographic maps at $z=2-3$ over a total of 837 deg$^2$ in two blank fields that are developed with Ly$\alpha$ forest absorptions of 14,736 background Sloan Digital Sky Survey (SDSS) quasars at $z$=2.08-3.67. Using the tomographic maps, we investigate the large-scale ($\gtrsim 10$ $h^{-1}$cMpc) average {\sc Hi} radial profiles and two-direction profiles of the line-of-sight (LoS) and transverse (TS) directions around galaxies and AGN at $z=2-3$ identified by the Hobby-Eberly Telescope Dark Energy eXperiment (HETDEX) and SDSS surveys, respectively.The peak of the {\sc Hi} radial profile around galaxies is lower than the one around AGN, suggesting that the dark-matter halos of galaxies are less massive on average than those of AGN. The LoS profile of AGN is narrower than the TS profile, indicating the Kaiser effect. There exist ionized outskirts at $\gtrsim 30$ $h^{-1}$cMpc beyond {\sc Hi} rich structures of galaxies and AGN found in the LoS profiles that can be explained by the influence of high energy photons propagating over a long distance. Our findings indicate that the {\sc Hi} radial profile of AGN has transitions from proximity zones ($\lesssim$ a few $h^{-1}$cMpc) to the {\sc Hi} rich structures ($\sim 1-30$ $h^{-1}$cMpc) and the ionized outskirts ($\gtrsim 30$ $h^{-1}$cMpc). Although there is no significant dependence of AGN types (type-1 vs. type-2) on the {\sc Hi} profiles, the peaks of the radial profiles anti-correlate with AGN luminosities, suggesting that AGN's ionization effects are stronger than the gas mass differences.

We present QUIJOTE intensity and polarization maps in four frequency bands centred around 11, 13, 17 and 19GHz, and covering approximately 29000 deg$^2$, including most of the Northern sky region. These maps result from 9000 h of observations taken between May 2013 and June 2018 with the first QUIJOTE instrument (MFI), and have angular resolutions of around $1^\circ$, and sensitivities in polarization within the range 35-40 $\mu$K per 1-degree beam, being a factor $\sim 2$-$4$ worse in intensity. We discuss the data processing pipeline employed, and the basic characteristics of the maps in terms of real space statistics and angular power spectra. A number of validation tests have been applied to characterise the accuracy of the calibration and the residual level of systematic effects, finding a conservative overall calibration uncertainty of 5%. We also discuss flux densities for four bright celestial sources (Tau A, Cas A, Cyg A and 3C274) which are often used as calibrators at microwave frequencies. The polarization signal in our maps is dominated by synchrotron emission. The distribution of spectral index values between the 11GHz and WMAP 23GHz map peaks at $\beta=-3.09$ with a standard deviation of 0.14. The measured BB/EE ratio at scales of $\ell=80$ is $0.26\pm 0.07$ for a Galactic cut $|b|>10^\circ$. We find a positive TE correlation for 11GHz at large angular scales ($\ell \lesssim 50$), while the EB and TB signals are consistent with zero in the multipole range $30 \lesssim \ell \lesssim 150$. The maps discussed in this paper are publicly available.

We present new intensity and polarisation maps obtained with the QUIJOTE experiment towards the Galactic regions W49, W51 and IC443, covering the frequency range from 10 to 20 GHz at $\sim$ 1 deg angular resolution, with a sensitivity in the range 35-79 ${\mu}$K/beam for total intensity and 13-23 ${\mu}$K/beam for polarisation. For each region, we combine QUIJOTE maps with ancillary data at frequencies ranging from 0.4 to 3000 GHz, reconstruct the spectral energy distribution and model it with a combination of known foregrounds. We detect anomalous microwave emission (AME) in total intensity towards W49 at 4.7${\sigma}$ and W51 at 4.0${\sigma}$ with peak frequencies ${\nu}_{AME}$ = (20.0 $\pm$ 1.4) GHz and ${\nu}_{AME}$ = (17.7 $\pm$ 3.6) GHz respectively; this is the first detection of AME towards W51. The contamination from ultra-compact HII regions to the residual AME flux density is estimated at 10% in W49 and 5% in W51, and does not rule out the AME detection. The polarised SEDs reveal a synchrotron contribution with spectral indices ${\alpha}_s$ = -0.67 $\pm$ 0.10 in W49 and ${\alpha}_s$ = -0.51 $\pm$ 0.07 in W51, ascribed to the diffuse Galactic emission and to the local supernova remnant respectively. Towards IC443 in total intensity we measure a broken power-law synchrotron spectrum with cut-off frequency ${\nu}_{0,s}$ = (114 $\pm$ 73) GHz, in agreement with previous studies; our analysis, however, rules out any AME contribution which had been previously claimed towards IC443. No evidence of polarised AME emission is detected in this study.

The Haze is an excess of microwave intensity emission surrounding the Galactic centre. It is spatially correlated with the $\gamma$-ray Fermi bubbles, and with the S-PASS radio polarization plumes, suggesting a possible common provenance. The models proposed to explain the origin of the Haze, including energetic events at the Galactic centre and dark matter decay in the Galactic halo, do not yet provide a clear physical interpretation. In this paper we present a re-analysis of the Haze including new observations from the Multi-Frequency Instrument (MFI) of the Q-U-I JOint TEnerife (QUIJOTE) experiment, at 11 and 13 GHz. We analyze the Haze in intensity and polarization, characterizing its spectrum. We detect an excess of diffuse intensity signal ascribed to the Haze. The spectrum at frequencies 11$\,\leq\nu\leq\,$70 GHz is a power-law with spectral index $\beta^{\rm H}=-2.79\pm0.08$, which is flatter than the Galactic synchrotron in the same region ($\beta^{\rm S}=-2.98\pm0.04$), but steeper than that obtained from previous works ($\beta^{\rm H}\sim-2.5$ at 23$\,\leq\,\nu\leq\,$70 GHz). We also observe an excess of polarized signal in the QUIJOTE-MFI maps in the Haze area. This is a first hint detection of polarized Haze, or a consequence of curvature of the synchrotron spectrum in that area. Finally, we show that the spectrum of polarized structures associated with Galactic centre activity is steep at low frequencies ($\beta \sim -3.2$ at 2.3 $\leq\nu\leq$ 23 GHz), and becomes flatter above 11 GHz.

The QUIJOTE-MFI Northern Hemisphere Wide-Survey has provided maps of the sky above declinations $-30^\circ$ at 11, 13, 17 and 19$\,$GHz. These data are combined with ancillary data to produce Spectral Energy Distributions in intensity in the frequency range 0.4--3\,000$\,$GHz on a sample of 52 candidate compact sources harbouring anomalous microwave emission (AME). We apply a component separation analysis at 1$^\circ$ scale on the full sample from which we identify 44 sources with high AME significance. We explore correlations between different fitted parameters on this last sample. QUIJOTE-MFI data contribute to notably improve the characterisation of the AME spectrum, and its separation from the other components. In particular, ignoring the 10--20\,GHz data produces on average an underestimation of the AME amplitude, and an overestimation of the free-free component. We find an average AME peak frequency of 23.6 $\pm$ 3.6$\,$GHz, about 4$\,$GHz lower than the value reported in previous studies. The strongest correlation is found between the peak flux density of the thermal dust and of the AME component. A mild correlation is found between the AME emissivity ($A_{\rm AME}/\tau_{250}$) and the interstellar radiation field. On the other hand no correlation is found between the AME emissivity and the free-free radiation Emission Measure. Our statistical results suggest that the interstellar radiation field could still be the main driver of the intensity of the AME as regards spinning dust excitation mechanisms. On the other hand, it is not clear whether spinning dust would be most likely associated with cold phases of the interstellar medium rather than with hot phases dominated by free-free radiation.

We derive linearly polarized astrophysical component maps in the Northern Sky from the QUIJOTE-MFI data at 11 and 13 GHz in combination with the WMAP K and Ka bands (23 and 33 GHz) and all Planck polarized channels (30-353 GHz), using the parametric component separation method B-SeCRET. The addition of QUIJOTE-MFI data significantly improves the parameter estimation of the low-frequency foregrounds, especially the estimation of the synchrotron spectral index, $\beta_s$. We present the first detailed $\beta_s$ map of the Northern Celestial Hemisphere at a smoothing scale of $2^{\circ}$. We find statistically significant spatial variability across the sky. We obtain an average value of $-3.08$ and a dispersion of $0.13$, considering only pixels with reliable goodness-of-fit. The power law model of the synchrotron emission provides a good fit to the data outside the Galactic plane but fails to track the complexity within this region. Moreover, when we assume a synchrotron model with uniform curvature, $c_s$, we find a value of $c_s = -0.0797 \pm 0.0012$. However, there is insufficient statistical significance to determine which model is favoured, either the power law or the power law with uniform curvature. Furthermore, we estimate the thermal dust spectral parameters in polarization. Our CMB, synchrotron, and thermal dust maps are highly correlated with the corresponding products of the PR4 Planck release, although some large-scale differences are observed in the synchrotron emission. Finally, we find that the $\beta_s$ estimation in the high signal-to-noise synchrotron emission areas is prior-independent while, outside these regions, the prior governs the $\beta_s$ estimation.

We present the catalogue of Q-U-I JOint TEnerife (QUIJOTE) Wide Survey radio sources extracted from the maps of the Multi-Frequency Instrument compiled between 2012 and 2018. The catalogue contains 786 sources observed in intensity and polarization, and is divided into two separate sub-catalogues: one containing 47 bright sources previously studied by the \emph{Planck} collaboration and an extended catalogue of 739 sources either selected from the \emph{Planck} Second Catalogue of Compact Sources or found through a blind search carried out with a Mexican Hat 2 wavelet. A significant fraction of the sources in our catalogue (38.7 per cent) are within the $|b| \leq 20^\circ$ region of the Galactic plane. We determine statistical properties for those sources that are likely to be extragalactic. We find that these statistical properties are compatible with currently available models, with a $\sim$1.8 Jy completeness limit at 11 GHz. We provide the polarimetric properties of (38, 33, 31, 23) sources with P detected above the $99.99\%$ significance level at (11, 13, 17, 19) GHz, respectively. Median polarization fractions are in the $2.8$-$4.7$\% range in the 11-19 GHz frequency interval. We do not distinguish between Galactic and extragalactic sources here. The results presented here are consistent with those reported in the literature for flat- and steep-spectrum radio sources.

We use recent microlensing observations toward the central bulge of the Galaxy to probe the overall stellar plus brown dwarf initial mass function (IMF) in these regions well within the brown dwarf domain. We find that the IMF is consistent with the same Chabrier (2005) IMF characteristic of the Galactic disk. In contrast, other IMFs suggested in the literature overpredict the number of short-time events, thus of very-low mass stars and brown dwarfs, compared with observations. This, again, supports the suggestion that brown dwarfs and stars form predominantly via the same mechanism. We show that claims for different IMFs in the stellar and substellar domains rather arise from an incorrect parameterization of the IMF. Furthermore, we show that the IMF in the central regions of the bulge seems to be bottom-heavy, as illustrated by the large number of short-time events compared with the other regions. This recalls our previous analysis of the IMF in massive early type galaxies and suggests the same kind of two-phase formation scenario, with the central bulge initially formed under more violent, burst-like conditions than the rest of the Galaxy.

We examine the binding energy of massive stripped-envelope core collapse supernova (SECCSN) progenitors with the stellar evolution code MESA, and find that only the jittering jets explosion mechanism can account for explosions where carbon-oxygen cores with masses of $>20M_\odot$ collapse to leave a neutron star (NS) remnant. We calculate the binding energy at core collapse under the assumption that the remnant is a NS. Namely, stellar gas above mass coordinate of $~1.5-2.5M_\odot$ is ejected in the explosion. We find that the typical binding energy of the ejecta of stripped-envelope progenitors with carbon-oxygen core masses of $M_{CO}>20M_\odot$ is $E_{bind}>2x10^{51} erg$. Since only jet-driven explosion mechanisms can supply such high energies, we conclude that jets must explode such cores. We apply our results to SN 2020qlb, which is a SECCSN with a claimed core mass of $~30-50M_\odot$, and conclude that the jittering jets explosion mechanism best account for such an explosion that leaves a NS.

The search for life elsewhere in the universe is one of the central aims of science in the 21st century. While most of this work is aimed at planets orbiting other stars, the search for life in our own Solar System is an important part of this endeavour. Venus is often thought to have too harsh an environment for life, but it may have been a more hospitable place in the distant past. If life evolved there in the past then the cloud decks of Venus are the only remaining niche where life as we know it might survive today. The discovery of the molecule phosphine, PH$_3$, in these clouds has reinvigorated research looking into the possibility of life in the clouds. In this review we examine the background to studies of the possibility of life on Venus, discuss the discovery of phosphine, review conflicting and confirming observations and analyses, and then look forward to future observations and space missions that will hopefully provide definitive answers as to the origin of phosphine on Venus and to the question of whether life might exist there.

In this paper we will discuss the application of optimal filtering techniques for the adaptive optics system of the LBT telescope. We have studied the application of both Kalman and H$_\infty$ filters to estimate the temporal evolution of the phase perturbations due to the atmospheric turbulence and the telescope vibrations on tip/tilt modes. We will focus on the H$_\infty$ filter and on its advantages and disadvantages over the Kalman filter.

We present the e-TidalGCs Project which aims at modeling and predicting the extra-tidal features surrounding all Galactic globular clusters for which 6D phase space information, masses and sizes are available (currently 159 globular clusters). We focus the analysis and presentation of the results on the distribution of extra-tidal material on the sky, and on the different structures found at different heliocentric distances. We emphasize the wide variety of morphologies found: beyond the canonical tidal tails, our models reveal that the extra-tidal features generated by globular clusters take a wide variety of shapes, from thin and elongated shapes, to thick, and complex halo-like structures. We also compare some of the most well studied stellar streams found around Galactic globular clusters to our model predictions, namely those associated to the clusters NGC 3201, NGC 4590, NGC 5466 and Pal 5. Additionally, we investigate how the distribution and extension in the sky of the simulated streams vary with the Galactic potential by making use of three different models, containing or not a central spheroid, or a stellar bar. Overall, our models predict that the mass lost by the current globular cluster population in the field from the last 5 Gyrs is between $0.3-2.1\times10^{7}M_{\odot}$, an amount comparable between 7-55 % of current mass. Most of this lost mass is found in the inner Galaxy, with the half-mass radius of this population being between 4-6 kpc. The outputs of the simulations will be publicly available, at a time when the ESA Gaia mission and complementary spectroscopic surveys are delivering exquisite data to which these models can be compared.

There exists a well known relation between the mass of the supermassive black hole (SMBH) in the center of galaxies and their bulge mass or central velocity dispersion. This suggests a co-evolution between SMBH and their galaxy hosts. Our aim is to study this relation specifically for radio loud galaxies, and as a function of redshift $z$. We selected a sample of radio-galaxies and AGN by cross-matching the low radio frequency sources from VLA FIRST with spectroscopically confirmed sources from wide field surveys including SDSS DR14 ugriz and DES DR2 grzY in optical, WISE in infrared, and the Galaxy And Mass Assembly (GAMA) spectroscopic survey. Keeping only high signal to noise (S/N) sources in WISE magnitudes, and those with broad emission lines, we selected a sub sample of 42 radio sources, all with infrared-to-optical counterparts, for which we characterized the stellar, star formation, and black hole properties. We estimated the central SMBH mass, the stellar mass $M_\star$, the Eddington ratio $\eta$ and the jet power, $Q_{\rm jet}$. The relation between SMBH mass, $M_\star$, $\eta$ and $z$ are put into context by comparing them with scaling relations ($M_{\rm BH}$--$M_{\star}$, $M_{\rm BH}/M_\star$--$z$, $M_{\rm BH}$--$Q_{\rm jet}$ and $Q_{\rm jet}$--$\eta$) from the literature. An evolutionary scenario where radio-mode AGN feedback (or the cluster environments) regulate the accretion onto the SMBHs and the stellar mass assembly of the radio sources is discussed, which may explain the observed phenomenology, and in particular the presence of radio sources with high $M_{\rm BH}/M_\star$ ratios. This pilot study represents a benchmark for future ones using wide field surveys such as Euclid and the Vera Rubin telescope.

We present the source structure analysis of 11 calibrator sources in the Southern Hemisphere at 2.3 (S-band) and 8.4 GHz (X-band). We used multi-epoch very long baseline interferometry source maps available in the radio fundamental catalog to analyse jet-structure variability and also used fluxes from the Goddard Space Flight Center database to see whether these two complement each other or not. Also, total fluxes from the maps were plotted with the fluxes from the database. The S/X-band light curve analysis provides a more clear picture of the structural variability at the S/X-band also indicates the possibility of the "core-shift" phenomenon. We found jet-like structures in the majority of the sources in the sample.

Many of the confirmed short period super-Earths and smaller sub-Neptunes are sufficiently irradiated for the surface silicates to be sustained in a long-lasting molten state. While there is no direct evidence of magma ocean influence on exoplanets, theory suggests that due to outgassing and diverse evolution paths, a wide range of resulting atmospheric compositions should be possible. Atmospheric contamination caused by the outgassing of the underlying magma ocean is potentially detectable using low resolution spectroscopy. The James Webb Space Telescope provides the necessary spectral coverage and sensitivity to characterise smaller planets, including lava worlds. In this light, we assess observability of outgassed silicates submerged in volatile atmospheres on the edge of the evaporation valley. By placing a hypothetical 2 R${_\oplus}$ planet around a Sun-like star, we self-consistently model, in 1-D, a wide range of potential atmospheric compositions, including thermal structure and outgassing. We focus on atmospheres rich in H, C and N. We assess diverse chemistry of silicates and volatiles, and what features of outgassed species could be detected via emission spectroscopy using MIRI LRS. Results indicate that even for substantial volatile envelopes, strong in infrared opacity, the presence of silicates causes deep thermal inversions, affecting emission. Similar to pure lava worlds, SiO remains the only outgassed species with major infrared, 5 and 9 \textmu m, bands. However, even a small amount of volatiles, especially of H2O and H-, may hinder its observability. We also find that the C/O ratio plays a large role in determining the abundance of SiO. Detecting SiO on a strongly irradiated planet could indicate an atmosphere with high metallicity and a low C/O ratio, which may be a result of efficient interaction between the atmosphere and the underlying melt.

The next generation of spectroscopic surveys is expected to provide spectra for hundreds of thousands of white dwarf (WD) candidates in the upcoming years. Currently, spectroscopic classification of white dwarfs is mostly done by visual inspection, requiring substantial amounts of expert attention. We propose a data-driven pipeline for fast, automatic selection and spectroscopic classification of WD candidates, trained using spectroscopically confirmed objects with available Gaia astrometry, photometry, and Sloan Digital Sky Survey (SDSS) spectra with signal-to-noise ratios $\geq9$. The pipeline selects WD candidates with improved accuracy and completeness over existing algorithms, classifies their primary spectroscopic type with $\gtrsim 90\%$ accuracy, and spectroscopically detects main sequence companions with similar performance. We apply our pipeline to the Gaia Data Release 3 cross-matched with the SDSS Data Release 17 (DR17), identifying 424 096 high-confidence WD candidates and providing the first catalogue of automated and quantifiable classification for 36 523 WD spectra. Both the catalogue and pipeline are made available online. Such a tool will prove particularly useful for the undergoing SDSS-V survey, allowing for rapid classification of thousands of spectra at every data release.

Luminous Red Galaxies, or LRGs, are representative of the most massive galaxies and were originally selected in the Sloan Digital Sky Survey as good tracers of large scale structure. They are dominated by by uniformly old stellar populations, have low star formation rates, early type morphologies, and little cold gas. Despite having old stellar populations and little in situ star formation, studies have shown that they have grown their stellar mass since z=1, implying that they grow predominantly via the accretion of satellites. Tests of this picture have been limited because of the lack of deep imaging data sets that both covers a large enough area of the sky to contain substantial numbers of LRGs and that also is deep enough to detect faint satellites. We use the 25 square degree Early Data Release (EDR) of the DESI Legacy Imaging Surveys to characterize the satellite galaxy population of LRGs out to z=0.65. The DESI Legacy Imaging Surveys are comprised of grz imaging to 2-2.5 mag deeper than SDSS and with better image quality. We use a new statistical background technique to identify excess populations of putative satellite galaxies around 1823 LRGs at 0.2<z<0.65. In three redshift and luminosity bins we measure the numbers of satellite galaxies and their r- color distribution down to rest-frame $g$-band luminosity limits at least 3.6 times fainter than L*. In addition, we develop a forward modeling technique and apply it to constrain the mean number of satellites in each of our redshift and luminosity bins. Finally, we use these estimates to determine the amount of stellar mass growth in LRGs down to the local Universe.

In this chapter, we explore the origins of cometary material and discuss the clues cometary composition provides in the context of the origin of our solar system. The review focuses on both cometary refractory and volatile materials, which jointly provide crucial information about the processes that shaped the solar system into what it is today. Both areas have significantly advanced over the past decade. We also view comets more broadly and discuss compositions considering laboratory studies of cometary materials, including interplanetary dust particles and meteoritic material that are potential cometary samples, along with meteorites, and in situ/remote studies of cometary comae. In our review, we focus on key areas from elemental/molecular compositions, isotopic ratios, carbonaceous and silicate refractories, short-lived radionuclides, and solar system dynamics that can be used as probes of the solar birth environment. We synthesize this data that points towards the birth of our solar system in a clustered star-forming environment.

The information on Galactic assembly time is imprinted on the chemodynamics of globular clusters. This makes them important probes that help us to understand the formation and evolution of the Milky Way. Discerning between in-situ and ex-situ origin of these objects is difficult when we study the Galactic bulge, which is the most complex and mixed component of the Milky Way. To investigate the early evolution of the Galactic bulge, we analysed the globular cluster NGC 6355. We derived chemical abundances and kinematic and dynamic properties by gathering information from high-resolution spectroscopy with FLAMES-UVES, photometry with the Hubble Space Telescope, and Galactic dynamic calculations applied to the globular cluster NGC 6355. We derive an age of $13.2\pm1.1$ Gyr and a metallicity of [Fe/H]$=-1.39\pm0.08$ for NGC 6355, with $\alpha$-enhancement of [$\alpha$/Fe]$=+0.37\pm0.11$. The abundance pattern of the globular cluster is compatible with bulge field RR Lyrae stars and in-situ well-studied globular clusters. The orbital parameters suggest that the cluster is currently confined within the bulge volume when we consider a heliocentric distance of $8.54\pm0.19$ kpc and an extinction coefficient of $R_V = 2.84\pm0.02$. NGC 6355 is highly likely to come from the main bulge progenitor. {Nevertheless, it still} has a low probability of being formed from an accreted event because its age is uncertain and because of the combined [Mg/Mn] [Al/Fe] abundance. Its relatively low metallicity with respect to old and moderately metal-poor inner Galaxy clusters may suggest a low-metallicity floor for globular clusters that formed in-situ in the early Galactic bulge.

In the FLARES (First Light And Reionisation Epoch Simulations) suite of hydrodynamical simulations, we find the high redshift ($z>5$) intrinsic size-luminosity relation is, surprisingly, negatively sloped. However, after including the effects of dust attenuation we find a positively sloped UV observed size-luminosity relation in good agreement with other simulated and observational studies. In this work, we extend this analysis to probe the underlying physical mechanisms driving the formation and evolution of the compact galaxies driving the negative size-mass/size-luminosity relation. We find the majority of compact galaxies ($R_{1/2, \star}< 1 \mathrm{pkpc}$), which drive the negative slope of the size-mass relation, have transitioned from extended to compact sizes via efficient centralised cooling, resulting in high specific star formation rates in their cores. These compact stellar systems are enshrouded by non-star forming gas distributions as much as $100\times$ larger than their stellar counterparts. By comparing with galaxies from the EAGLE simulation suite, we find that these extended gas distributions `turn on' and begin to form stars between $z=5$ and $z=0$ leading to increasing sizes, and thus the evolution of the size-mass relation from a negative to a positive slope. This explicitly demonstrates the process of inside-out galaxy formation in which compact bulges form earlier than the surrounding discs.

Theories beyond general relativity (GR) modify the propagation of gravitational waves (GWs). In some, inhomogeneities (aka. gravitational lenses) allow interactions between the metric and additional fields to cause lens-induced birefringence (LIB): a different speed of the two linear GW polarisations ($+$ and $\times$). Inhomogeneities then act as non-isotropic crystals, splitting the GW signal into two components whose relative time delay depends on the theory and lens parameters. Here we study the observational prospects for GW scrambling, i.e when the time delay between both GW polarisations is smaller than the signal's duration and the waveform recorded by a detector is distorted. We analyze the latest LIGO-Virgo-KAGRA catalog, GWTC-3, and find no conclusive evidence for LIB. The highest log Bayes factor that we find in favour of LIB is $3.21$ for GW$190521$, a particularly loud but short event. However, when accounting for false alarms due to (Gaussian) noise fluctuations, this evidence is below 1-$\sigma$. The tightest constraint on the time delay is $<0.51$ ms (90% C.L.) from GW$200311\_115853$. From the non-observation of GW scrambling, we constrain the optical depth for LIB, accounting for the chance of randomly distributed lenses (eg. galaxies) along the line of sight. Our LIB constraints on a (quartic) scalar-tensor Horndeski theory are more stringent than solar system tests for a wide parameter range and comparable to GW170817 in some limits. Interpreting GW190521 as an AGN binary (i.e. taking an AGN flare as a counterpart) allows even more stringent constraints. Our results demonstrate the potential and high sensitivity achievable by tests of GR, based on GW lensing.

The dynamics of a binary system moving in the background of a black hole is affected by tidal forces. In this work, for the Kerr black hole, we derive the electric and magnetic tidal moments at quadrupole order, where the latter are computed for the first time in full generality. We make use of these moments in the scenario of a hierarchical triple system made of a Kerr black hole and an extreme-mass ratio binary system consisting of a Schwarzschild black hole and a test particle. We study how the secular dynamics of the test particle in the binary system is distorted by the presence of tidal forces from a much larger Kerr black hole. Our treatment includes strong gravitational effects beyond the post-Newtonian approximation both for the binary system and for the tidal forces since the binary system is allowed to be close to the event horizon of the Kerr black hole. We compute the shifts in the physical quantities for the secular dynamics of the test particle and show that they are gauge-invariant. In particular, we apply our formalism to the innermost stable circular orbit for the test particle and to the case of the photon sphere. Our results are relevant for the astrophysical situation in which the binary system is in the vicinity of a supermassive black hole.

The COSINE-100 experiment has been operating with 106 kg of low-background NaI(Tl) detectors to test the results from the DAMA/LIBRA experiment, which claims to have observed dark matter. However, since the background of the NaI(Tl) crystals used in the COSINE-100 experiment is 2-3 times higher than that in the DAMA detectors, no conclusion regarding the claimed observation from the DAMA/LIBRA experiment could be reached. Therefore, we plan to upgrade the current COSINE-100 experiment to the next phase, COSINE-200, by using ultra-low background NaI(Tl) detectors. The basic principle was already proved with the commercially available Astro-grade NaI powder from Sigma-Aldrich company. However, we have developed a mass production process of ultra-pure NaI powder at the Center for Underground Physics (CUP) of the Institute for Basic Science (IBS), Korea, using the direct purification of the raw NaI powder. We plan to produce more than 1,000 kg of ultra-pure powder for the COSINE200 experiment. With our crystal grower installed at CUP, we have successfully grown a low-background crystal using our purification technique for the NaI powder. We have assembled a low-background NaI(Tl) detector. In this article, we report the performance of this ultra-pure NaI(Tl) crystal detector produced at IBS, Korea.

In this work, we study the shadow boundary curves of rotating time-dependent black hole solutions which have well-defined Kerr and Vaidya limits. These solutions are constructed by applying the Newman-Janis algorithm to a spherically symmetric seed metric conformal to the Vaidya solution with a mass function that is linear in Eddington-Finkelstein coordinates. Equipped with a conformal Killing vector field, this class of solution exhibits separability of null geodesics, thus allowing one to develop an analytic formula for the boundary curve of its shadow. We find a simple power law describing the dependence of the mean radius and asymmetry factor of the shadow on the accretion rate. Applicability of our model to recent Event Horizon Telescope observations of M87${}^*$ and Sgr A${}^*$ is also discussed.

In this paper we apply the symmetry principle in order to search for an alternative unified explanation of several cosmological puzzles such as the present stage of accelerated expansion of the Universe and the Hubble tension issue, among others. We argue that Weyl gauge symmetry, being a manifest symmetry of gauge invariant theories of gravity operating on Weyl integrable geometry spacetimes, may be an actual (unbroken) symmetry of our present Universe. This symmetry may be at the core of a phenomenologically feasible explanation of modern fundamental issues arising within the framework of general relativity and of its known modifications.

Using the analytic, static and spherically symmetric metric for a Schwarzschild black hole immersed in dark matter (DM) halos with Hernquist type density distribution, we derive analytic formulae for the orbital period and orbital precession, the evolutions of the semi-latus rectum and the eccentricity for eccentric EMRIs with the environment of DM halos. We show how orbital precessions are decreased and even reverse the direction if the density of DM halo is large enough. The presence of local DM halos slows down the decrease of the semi-latus rectum and the eccentricity. Comparing the number of orbital cycles with and without DM halos over one-year evolution before the merger, we find that DM halos with the compactness as small as $10^{-4}$ can be detected. By calculating the mismatch between GW waveforms with and without DM halos, we show that we can use GWs from EMRIs in the environments of galaxies to test the existence of DM halos and detect the compactness as small as $10^{-5}$.

Using a quasi-spherical approximation of an affine-null metric adapted to an asymptotic Bondi inertial frame, we present high order approximations of the metric functions in terms of the specific angular momentum for a slowly rotating stationary and axi-symmetric vacuum spacetime. The metric is obtained by following the procedure of integrating the hierarchy of Einstein equations in a characteristic formulation utilizing master functions for the perturbations. It is further verified its equivalence with the Kerr metric in the slowly rotation approximation by carrying out an explicit transformation between the Boyer-Lindquist coordinates to the employed affine-null coordinates.

Phenomenological work in the last few years has provided significant support to the idea that the vacuum energy density (VED) is a running quantity with the cosmological evolution and that this running helps to alleviate the cosmological tensions afflicting the $\Lambda$CDM. On the theoretical side, recent devoted studies have shown that the properly renormalized $\rho_{\rm vac}$ in FLRW spacetime adopts the "running vacuum model" (RVM) form. While in three previous studies by two of us (CMP and JSP) such computations focused solely on scalar fields non-minimally coupled to gravity, in the present work we compute the spin-$1/2$ fermionic contributions and combine them both. The calculation is performed using a new version of the adiabatic renormalization procedure based on subtracting the UV divergences at an off-shell renormalization point $M$. The quantum scaling of $\rho_{\rm vac}$ with $M$ turns into cosmic evolution with the Hubble rate, $H$. As a result the "cosmological constant" $\Lambda$ appears in our framework as the nearly sustained value of $8\pi G(H)\rho_{\rm vac}(H)$ around (any) given epoch $H$, where $G(H)$ is the gravitational coupling, which is also running, although very mildly (logarithmically). We find that the VED evolution at present reads $\delta \rho_{\rm vac}(H)\sim \nu_{\rm eff} m_{\rm Pl}^2 \left(H^2-H_0^2 \right)\ (|\nu_{\rm eff}|\ll 1)$. The coefficient $\nu_{\rm eff}$ receives contributions from all the quantized fields, bosons and fermions. Remarkably, there also exist higher powers ${\cal O}(H^{6})$ which can trigger inflation in the early universe. Finally, the equation of state (EoS) of the vacuum receives also quantum corrections from bosons and fermion fields, shifting its value from -1. The remarkable consequence is that the EoS of the quantum vacuum may nowadays effectively look like quintessence.