Papers for Monday, Aug 08 2022

A critical analysis and comparison of different methods for obtaining point spread function (PSF) photometry are carried out. Deep ACS observations of NGC3370 were reduced using four distinct approaches. These reductions explore a number of methodological differences: software packages (DAOPHOT and DOLPHOT), input images (individual and stacked frames), PSF models (synthetic and empirical), and aperture correction methods (automatic and manual). A comparison of the photometry leads to the following results: 1) Photometric incompleteness between individual reductions shows only a minimal difference (<10%). 2) Statistical errors are 20% to 30% smaller for DAOPHOT runs on stacked frames than DOLPHOT runs on individual frames. 3) Statistical errors assigned directly by the photometry codes are 25% to 50% smaller than the errors measured from artificial star tests. 4) Systematic errors are magnitude dependent and become larger at the faint end, at the level of $\sigma_s\sim0.1$ mag. 5) The automatic aperture correction routines in DOLPHOT result in a significant systematic error ($\sigma_s \sim 0.05$ mag). 6) Individual reductions agree well at the 0.02 mag level when the systematic errors are properly corrected through artificial star tests. The reasonable agreement between the reductions leads to important implications that i) the reduction dependent errors can be reduced to a 1% level in the luminosity distance scale, and ii) the stacked frame photometry can be a good means to study non-variable stars in external galaxies.

A wealth of extragalactic populations completely missed at UV-optical wavelengths has been identified in the last decade, combining the deepest HST and Spitzer observations. These dark sources are thought to be very dusty and star-forming systems at 3<z<5, and major contributors to the stellar mass build up. JWST is now promising to detect such objects well beyond the end of the Epoch of Reionization. In this Letter we report an investigation of the deep JWST survey in the SMACS0723 cluster, analysing NIRCam and MIRI images. We search for sources in the F444W band that are undetected in the F200W catalogues. We characterise the main properties of these sources via detailed SED modelling that account for a wide set of parameters and star formation histories, after a careful determination of their photometry. Among a robust sample of 20 candidates, we identify a mixed population of very red sources. We highlight the identification of candidate evolved systems, with stellar masses M*~10^(9-11)Msun at 8<z<13 characterized by unexpectedly important dust content at those epochs (Av up to ~5.8mag), challenging current model predictions. We further identify an extremely red source (F200W-F440W~7mag) that can be reproduced only by the spectrum of a passive, quenched galaxy of M*~10^11.8Msun at z~5, filled of dust (Av~5mag).

X-ray selected samples are known to miss galaxy clusters that are gas poor and have a low surface brightness. This is different for the optically selected samples such as the X-ray Unbiased Selected Sample (XUCS). We characterise the origin of galaxy clusters that are gas poor and have a low surface-brightness by studying covariances between various cluster properties at fixed mass using hydrodynamic cosmological simulations. We extracted approx. 1800 galaxy clusters from a high-resolution Magneticum hydrodynamic cosmological simulation and computed covariances at fixed mass of the following properties: core-excised X-ray luminosity, gas fraction, hot gas temperature, formation redshift, concentration, galaxy richness, fossilness parameter, and stellar mass of the bright central galaxy. We also compared the correlation between concentration and gas fractions in non-radiative simulations, and we followed the trajectories of particles inside galaxy clusters to assess the role of AGN depletion on the gas fraction. In simulations and in observational data, differences in surface brightness are related to differences in gas fraction. Simulations show that the gas fraction strongly correlates with assembly time, in the sense that older clusters are gas poor. Clusters that formed earlier have lower gas fractions because the feedback of the active galactic nucleus ejected a significant amount of gas from the halo. When the X-ray luminosity is corrected for the gas fraction, it shows little or no covariance with other quantities. Older galaxy clusters tend to be gas poor and possess a low X-ray surface brightness because the feedback mechanism removes a significant fraction of gas from these objects. Moreover, we found that most of the $L_X$ covariance with the other quantities is explained by differences in the gas fraction.

Context. Extragalactic surveys are a key tool in better understanding the evolution of galaxies. Both deep and wide surveys are improving the emerging picture of physical processes that take place in and around galaxies and identifying which of these processes are the most important in shaping the properties of galaxies. Aims. The Lockman-SpReSO survey aims to provide one of the most complete optical spectroscopic follow-ups of far-infrared (FIR) sources detected by the $Herschel$ Space Observatory in the Lockman Hole field. Such a large optical spectroscopic sample of FIR-selected galaxies will supply valuable information about the relation between fundamental FIR and optical parameters (including extinction, star formation rate and gas metallicity). In this article, we introduce and provide an in-depth description of the Lockman-SpReSO survey and of its early results. Methods. We have selected FIR sources from the observations of the $Herschel$ telescope over the central 24 arcmin $\times$ 24 arcmin of the Lockman Hole field with an optical counterpart up to 24.5 $R_{\rm C}$(AB). The sample comprises 956 $Herschel$ FIR sources plus 188 interesting additional objects in the field. The faint component of the catalogue ($R_{\rm C}$(AB)$\geq$20) was observed using the OSIRIS instrument on the 10.4 m Gran Telescopio Canarias (GTC) in MOS mode. The bright component was observed using two multifibre spectrographs: the AF2-WYFFOS at the William Herschel Telescope (WHT) and the Hydra instrument at the WYIN telescope.

Compensated isocurvature perturbations (CIPs) are relative density perturbations in which a baryon-density fluctuation is accompanied by a dark-matter-density fluctuation such that the total-matter density is unperturbed. These fluctuations can be produced primordially if multiple fields are present during inflation, and therefore they can be used to differentiate between different models for the early Universe. Kinetic Sunyaev-Zeldovich (kSZ) tomography allows for the reconstruction of the radial-velocity field of matter as a function of redshift. This technique can be used to reconstruct the total-matter-overdensity field, independent of the galaxy-density field obtained from large-scale galaxy surveys. We leverage the ability to measure the galaxy- and matter-overdensity fields independently to construct a minimum-variance estimator for the primordial CIP amplitude, based on a mode-by-mode comparison of the two measurements. We forecast that a configuration corresponding to CMB-S4 and VRO will be able to detect (at $2\sigma$) a CIP amplitude $A$ (for a scale-invariant power spectrum) as small as $A\simeq 5\times 10^{-9}$. Similarly, a configuration corresponding to SO and DESI will be sensitive to a CIP amplitude $A\simeq 1\times 10^{-7}$. These values are to be compared to current constraints $A \leq {\cal O}(0.01)$.

This paper presents the first census of Galactic post-Asymptotic Giant Branch stars in the HR diagram. We combined Gaia DR3 parallax-based distances with extinction corrected integrated fluxes, and derived luminosities for a sample of 185 stars that had been proposed to be post-AGB stars in the literature. The luminosities allow us to create an HR diagram containing the largest number of post-AGB candidate objects to date. A significant fraction of the objects fall outside the typical luminosity range as covered by theoretical evolutionary post-AGB tracks as well as observed for Planetary Nebula central stars. These include massive evolved supergiants and lower luminosity objects. Here we highlight the fact that one third of the post-AGB candidates is underluminous and we identify these with the recently recognised class of post-Red Giant Branch objects thought to be the result of binary evolution.

We present continued analysis of a sample of low-redshift iron low-ionization broad absorption-line quasars (FeLoBALQs). Choi et al. (2022) presented $SimBAL$ spectral analysis of BAL outflows in 50 objects. Leighly et al. (2022) analyzed optical emission lines of 30 of those 50 objects and found that they are characterized by either a high accretion rate ($L_\mathrm{Bol}/L_\mathrm{Edd}>0.3$) or low accretion rate ($0.03<L_\mathrm{Bol}/L_\mathrm{Edd}<0.3$). We report that the outflow velocity is inversely correlated with the BAL location among the high accretion rate objects, with the highest velocities observed in the parsec-scale outflows. In contrast, the low Eddington ratio objects showed the opposite trend. We confirmed the known relationship between outflow velocity and $L_\mathrm{Bol}/L_\mathrm{Edd}$, and found that the scatter plausibly originates in the force multiplier (launch radius) in the low (high) accretion rate objects. A log volume filling factor between $-6$ and $-4$ was found in most outflows, but was as high as $-1$ for low-velocity compact outflows. We investigated the relationship between the observed [O III] emission and that predicted from the BAL gas. We found that these could be reconciled if the emission-line covering fraction depends on Seyfert type and BAL location. The difference between the predicted and observed [O III] luminosity is correlated with the outflow velocity, suggesting that [O III] emission in high Eddington ratio objects may be broad and hidden under Fe II emission. We suggest that the physical differences in the outflow properties as a function of location in the quasar and accretion rate point to different formation, acceleration, and confinement mechanisms for the two FeLoBALQ types.

Most observed stars are part of a multiple star system, but the formation of such systems and the role of environment and various physical processes is still poorly understood. We present a suite of radiation-magnetohydrodynamic simulations of star-forming molecular clouds from the STARFORGE project that include stellar feedback with varied initial surface density, magnetic fields, level of turbulence, metallicity, interstellar radiation field, simulation geometry and turbulent driving. In our fiducial cloud the raw simulation data reproduces the observed multiplicity fractions for Solar-type and higher mass stars, similar to previous works. However, after correcting for observational incompleteness the simulation under-predicts these values. The discrepancy is likely due to the lack of disk fragmentation, as the simulation only resolves multiples that form either through capture or core fragmentation. The raw mass distribution of companions is consistent with randomly drawing from the initial mass function for the companions of $>1\,\mathrm{M_\odot}$ stars, however, accounting for observational incompleteness produces a flatter distribution similar to observations. We show that stellar multiplicity changes as the cloud evolves and anti-correlates with stellar density. This relationship also explains most multiplicity variations between runs, i.e., variations in the initial conditions that increase stellar density (increased surface density, reduced turbulence) decrease multiplicity. While other parameters, such as metallicity, interstellar radiation, and geometry significantly affect the star formation history or the IMF, varying them produces no clear trend in stellar multiplicity properties.

There is a class of binary post-AGB stars (binary system including a post-AGB star) that are surrounded by Keplerian disks and outflows resulting from gas escaping from the disk. To date, there are seven sources that have been studied in detail through interferometric millimeter-wave maps of CO lines (ALMA/NOEMA). For the cases of the Red Rectangle, IW Carinae, IRAS 08544-4431, and AC Herculis, it is found that around greater than 85% of the total nebular mass is located in the disk with Keplerian dynamics. The remainder of the nebular mass is located in an expanding component. This outflow is probably a disk wind consisting of material escaping from the rotating disk. These sources are the disk-dominated nebulae. On the contrary, our maps and modeling of 89 Herculis, IRAS 19125+0343, and R Scuti, which allowed us to study their morphology, kinematics, and mass distribution, suggest that, in these sources, the outflow clearly is the dominant component of the nebula (around 75% of the total nebular mass), resulting in a new subclass of nebulae around binary post-AGB stars: the outflow-dominated sources.Besides CO, the chemistry of this type of source has been practically unknown thus far. We also present a very deep single-dish radio molecular survey in the 1.3, 2, 3, 7, and 13 mm bands (around 600 h of telescope time). Our results and detections allow us to classify our sources as O- or C-rich. We also conclude that the calculated abundances of the detected molecular species other than CO are particularly low, compared with AGB stars. This fact is very significant in those sources where the rotating disk is the dominant component of the nebula.

The Near-Infrared Spectrograph (NIRSpec) is one of the four focal plane instruments on the James Webb Space Telescope which was launched on Dec. 25, 2021. We present an overview of the as-run NIRSpec commissioning campaign, with particular emphasis on the sequence of activities that led to the verification of all hardware components of NIRSpec. We also discuss the mechanical, thermal, and operational performance of NIRSpec, as well as the readiness of all NIRSpec observing modes for use in the upcoming JWST science program.

We study the population of star-forming clumps in A521-sys1, a $\rm z=1.04$ system gravitationally lensed by the foreground ($\rm z=0.25$) cluster Abell 0521. The galaxy presents one complete counter--image with a mean magnification of $\rm \mu\sim4$ and a wide arc containing two partial images of A521-sys1 with magnifications reaching $\rm \mu>20$, allowing the investigations of clumps down to scales of $\rm R_{eff}<50$ pc. We identify 18 unique clumps with a total of 45 multiple images. Intrinsic sizes and UV magnitudes reveal clumps with elevated surface brightnesses, comparable to similar systems at redshifts $\rm z\gtrsim1.0$. Such clumps account for $\sim40\%$ of the galaxy UV luminosity, implying that a significant fraction of the recent star-formation activity is taking place there. Clump masses range from $\rm 10^6\ M_\odot$ to $\rm 10^9\ M_\odot$ and sizes from tens to hundreds of parsec, resulting in mass surface densities from $10$ to $\rm 10^3\ M_\odot\ pc^{-2}$, with a median of $\rm \sim10^2\ M_\odot\ pc^{-2}$. These properties suggest that we detect star formation taking place across a wide range of scale, from cluster aggregates to giant star-forming complexes. We find ages of less than $100$ Myr, consistent with clumps being observed close to their natal region. The lack of galactocentric trends with mass, mass density, or age and the lack of old migrated clumps can be explained either by dissolution of clumps after few $\sim100$ Myr or by stellar evolution making them fall below the detectability limits of our data.

We investigate the dynamics of the accretion disks of young stars with fossil large-scale magnetic field. The author's magnetohydrodynamic (MHD) model of the accretion disks is generalized to consider the dynamical influence of the magnetic field on gas rotation speed and vertical structure of the disks. With the help of the developed MHD model, the structure of an accretion disk of a solar mass T Tauri star is simulated for different accretion rates $\dot{M}$ and dust grain sizes $a_d$. The simulations of the radial structure of the disk show that the magnetic field in the disk is kinematic, and the electromagnetic force does not affect the rotation speed of the gas for typical values $\dot{M}=10^{-8}\,M_\odot\,\mbox{yr}^{-1}$ and $a_d=0.1 \mu$m. In the case of large dust grains, $a_d\geq 1$ mm, the magnetic field is frozen into the gas and a dynamically strong magnetic field is generated at radial distances from the star $r\gtrsim 30$ au, the tensions of which slow down the rotation speed by $\lesssim 1.5$ % of the Keplerian velocity. This effect is comparable to the contribution of the radial gradient of gas pressure and can lead to the increase in the radial drift velocity of dust grains in the accretion disks. In the case of high accretion rate, $\dot{M}\geq 10^{-7}\,M_\odot\,\mbox{yr}^{-1}$, the magnetic field is also dynamically strong in the inner region of the disk, $r<0.2$ au. The simulations of the vertical structure of the disk show that, depending on the conditions on the surface of the disk, the vertical gradient of magnetic pressure can lead to both decrease and increase in the characteristic thickness of the disk as compared to the hydrostatic one by 5-20 %. The change in the thickness of the disk occurs outside the region of low ionization fraction and effective magnetic diffusion (`dead' zone), which extends from $r=0.3$ to $20$ au at typical parameters.

Globular clusters (GCs) play an important role in the formation and evolution of the Milky Way. New candidates are continuously found, particularly in the high-extinction low-latitude regions of the bulge, although their existence and properties have yet to be verified. In order to investigate the new GC candidates, we performed high-resolution NIR spectroscopy of stars toward the bulge using the IGRINS instrument at the Gemini-South telescope. We selected 15 and 10 stars near Camargo 1103 and 1106, respectively, which have recently been reported as metal-poor GC candidates in the bulge. In contrast to the classical approaches used in optical spectroscopy, we determined stellar parameters from a combination of line-depth ratios and the equivalent width of a CO line. The stellar parameters of the stars follow the common trends of nearby APOGEE stars in a similar magnitude range. We also determined the abundances of Fe, Na, Mg, Al, Si, S, K, Ca, Ti, Cr, Ni, and Ce through spectrum synthesis. There is no clear evidence of a grouping in RV-[Fe/H] space that would indicate the characterization of either object as metal-poor GCs. This result emphasizes the necessity of follow-up spectroscopy for new GC candidates toward the bulge, although we cannot completely rule out a low probability that we only observed nonmember stars. We also note discrepancies between the abundances of Al, Ca, and Ti when derived from the H- vs. the K-band spectra. Although the cause of this discrepancy is not clear, the effects of atmosphere parameters or NLTE are discussed. Our approach and results demonstrate that IGRINS spectroscopy is a useful tool for studying the chemical properties of stars toward the Galactic bulge with a statistical uncertainty in [Fe/H] of 0.03 dex, while the systematic error through uncertainties of atmospheric parameter is slightly larger than in measurements from optical spectroscopy.

The He I 1 micron line is a high excitation line which allows us to probe the innermost regions of protostellar disks, and to trace both accreting and outflowing material. We use X-Shooter observations of a sample of 107 young stars in the Lupus (1-3 Myr) and Upper Scorpius (5-10 Myr) star-forming regions to search for correlations between the line properties, as well as the disk inclination and accretion luminosity. We identified eight distinct profile types in the sample. We fitted Gaussian curves to the line features to measure the maximum velocities traced in absorption, the full-width half-maximum (FWHM) of the line features, and the Gaussian area of the features. We compare the proportion of each profile type in our sample to previous studies in Taurus. We find significant variations between Taurus and Lupus in the proportion of P Cygni and inverse P Cygni profiles, and between Lupus and Upper Scorpius in the number of emission-only and combination profile types. We find that the blue-shifted absorption features appear less blue-shifted at disk inclinations close to edge-on, but no such trend with inclination is observed in sources with only red-shifted features. Higher accretion rates were observed in sources with strong blue-shifted features which, along with the changes in the proportions of each profile type observed in the two regions, indicates that younger sources may drive stronger jets or winds. Overall, we observe variations in the proportion of each profile type and in the line properties which indicates and evolution of accretion and ejection signatures over time, and with source properties. These results confirm past works and models of the He I line, but for a larger sample and for multiple star-forming regions. The work highlights the power of the He I line as a probe of the gas in the innermost regions of the disk.

We perturbatively study the effect of non-Gaussianities on the mass fraction of primordial black holes (PBHs) at the time of formation by systematically taking its effect into account in the one-point probability distribution function of the primordial curvature perturbation. We focus on the bispectrum and trispectrum and derive formulas that describe their effects on the skewness and kurtosis of the distribution function. Then considering the case of narrowly peaked spectra, we obtain simple formulas that concisely express the effect of the bi- and trispectra. In particular, together with the $g_{\rm NL}$ and $\tau_{\rm NL}$ parameters of the trispectrum, we find that non-Gaussianity parameters for various types of the bispectrum are linearly combined to give an effective parameter, $f_{\rm NL}^{\rm eff}$, that determines the PBH mass fraction in the narrow spectral shape limit.

Most upcoming CMB experiments are planning to deploy between a few thousand and a few hundred thousand TES bolometers in order to drastically increase sensitivity and unveil the B-mode signal. Differential systematic effects and $1/f$ noise are two of the challenges that need to be overcome in order to achieve this result. In recent years, rotating Half-Wave Plates have become increasingly more popular as a solution to mitigate these effects, especially for those experiments that are targeting the largest angular scales. However, other effects may appear when a rotating HWP is being employed. In this paper we focus on HWP synchronous signals, which are due to intensity to polarization leakage induced by a rotating cryogenic multi-layer sapphire HWP employed as the first optical element of the telescope system. We use LiteBIRD LFT as a case study and we analyze the interaction between these spurious signals and TES bolometers, to determine whether this signal can contaminate the bolometer response. We present the results of simulations for a few different TES model assumptions and different spurious signal amplitudes. Modelling these effects is fundamental to find what leakage level can be tolerated and minimize non-linearity effects of the bolometer response.

Small grains play an essential role in astrophysical processes such as chemistry, radiative transfer, gas/dust dynamics. The population of small grains is mainly maintained by the fragmentation process due to colliding grains. An accurate treatment of dust fragmentation is required in numerical modelling. However, current algorithms for solving fragmentation equation suffer from an over-diffusion in the conditions of 3D simulations. To tackle this challenge, we developed a Discontinuous Galerkin scheme to solve efficiently the non-linear fragmentation equation with a limited number of dust bins.

Changing-look active galactic nucleus NGC~4151, which has attracted a lot of attention, is undergoing the second dramatic outburst stage in its evolutionary history. To investigate the geometry and kinematics of the broad-line region (BLR), and measure the mass of supermassive black hole in NGC~4151, we perform a seven-month photometric and spectroscopic monitoring program in 2020--2021, using the 2.4 m telescope at Lijiang Observatory. We successfully measure the time lags of the responses from broad \ha, \hb, \hg, \hei, and \heii\ emission lines to continuum variation, which are $7.63_{-2.62}^{+1.85}$, $6.21_{-1.13}^{+1.41}$, $5.67_{-1.94}^{+1.65}$, $1.59_{-1.11}^{+0.86}$, and $0.46_{-1.06}^{+1.22}$ days, respectively, following radial stratification. The ratios of time lags among these lines are $1.23 : 1.00 : 0.91 : 0.26 : 0.07$. We find that the continuum lag between the ultraviolet and optical bands can significantly affect the lag measurements of \hei\ and \heii. Virial and infalling gas motions coexist in this campaign, which is different from previous results, implying the evolutionary kinematics of BLR. Based on our measurements and previous ones in the literature, we confirm that the BLR of NGC~4151 is basically virialized. Finally, we compute the black hole mass through multiple lines, and the measurement from \hb\ to be $ 3.94_{-0.72}^{+0.90} \times 10^7 M_{\odot}$, which is consistent with previous results. The corresponding accretion rate is $0.02_{-0.01}^{+0.01} L_{\rm Edd} c^{-2}$, implying a sub-Eddington accretor.

This report presents a down-conversion method involving digital sideband separation for the Yuan Tseh Lee Array to double the processing bandwidth. The receiver consists of a MMIC HEMT LNA frontend operating at a wavelength of 3 mm, and sub-harmonic mixers that output signals at intermediate frequencies of 2 - 18 GHz. The sideband separation scheme involves an analog 90 degree hybrid followed by two mixers that provide down conversion of the IF signal to a pair of in phase (I) and quadrature (Q) signals in baseband. The I and Q baseband signals are digitized using 5 Giga sample per second analog to digital converters. A second hybrid is digitally implemented using field programmable gate arrays to produce two sidebands, each with a bandwidth of 1.6 GHz. The 2 x 1.6 GHz band can be tuned to cover any 3.2 GHz window within the aforementioned IF range of the array. Sideband rejection ratios (SRRs) above 20 dB can be obtained across the 3.2 GHz bandwidth by equalizing the power and delay between the I and Q baseband signals. Furthermore, SRRs above 30 dB can be achieved when calibration is applied.

Binaries with a white dwarf primary and a main sequence secondary can be used to test our understanding of both single and binary star evolution. A small fraction of such systems experienced a common-envelope phase from which they emerged with a relatively short orbital period. Here, we present the characterisation of an eclipsing post-common-envelope binary of this kind, TIC 60040774, based on the light curve provided by the Transiting Exoplanet Survey Satellite (TESS), multi-band photometry collated from the virtual observatory, and spectroscopic data obtained the Southern African Large Telescope (SALT). With an orbital period of $0.404807\pm0.000149$ days, this system consists of a young white dwarf paired with an M6.5 dwarf companion. We estimate the masses of the primary and secondary to be $0.598\pm0.029$ M$_{\odot}$ and $0.107\pm0.020$ M$_{\odot}$, while the effective temperatures are $14050\pm360$ K and $2759\pm50$ K, respectively. The eclipse ingress and egress profile is shallower than expected from a simple geometric model such that more precise high-cadence photometry is required to understand the nature of this system. Given the similarity of TIC 60040774 to systems like GK Vir and NN Ser, it will be worth tracking its eclipse times to check for the presence of one or more circumbinary planets.

We simulate the first minutes of the evolution of a binary-driven hypernova (BdHN) event, with a special focus on the associated accretion processes of supernova (SN) ejecta onto the newborn neutron star ($\nu$NS) and the NS companion. We calculate the rotational evolution of the $\nu$NS and the NS under the torques exerted by the accreted matter and the magnetic field. We take into account general relativistic effects and use realistic hypercritical accretion rates obtained from three-dimensional smoothed-particle-hydrodynamics (SPH) numerical simulations of the BdHN for a variety of orbital periods. We show that the rotation power of the $\nu$NS has a unique double-peak structure while that of the NS has a single peak. These peaks are of comparable intensity and can occur very close in time or even simultaneously depending on the orbital period and the initial angular momentum of the stars. We outline the consequences of the above features in the early emission and their consequent observation in long gamma-ray bursts (GRBs).

Eclipsing, spectroscopic double-lined (SB2) binaries remain to be the prime source of precise and accurate fundamental properties of stars. Furthermore, high-cadence spectroscopic observations of the eclipse phases allow us to resolve the Rossiter-McLaughlin effect whose modelling offers the means to probe spin-orbit misalignment in binaries. In this study, we develop the LSDBinary algorithm that is capable of working with both in-eclipse and out-of-eclipse spectra of SB2 binaries as input and delivers the LSD profiles, LSD-based model spectra, and precise RVs of both binary components as output. We offer an option to account for the Rossiter-McLaughlin effect in the calculation of the initial guess LSD profiles and components' flux ratio such that the effect can be modelled within the algorithm itself. We provide an extensive test of the LSDBinary software package on simulated spectra of artificial binaries. We study the effects of signal-to-noise-ratio of input spectra, resolving power of the instrument, uncertain atmospheric parameters of stars, and orbital properties of the binary system on the resulting LSD profiles and RVs measured from them. We find that atmospheric parameters have negligible effect on the shape of the computed LSD profiles while affecting mostly their global scaling. Our results are barely sensitive to signal-to-noise ratio of the input spectra provided they contain sufficient number of spectral lines, such as in A-type stars and later. Finally, the orbital inclination angle and components' radii ratio are found to have the largest effect on the shapes of the LSD profiles and RV curves extracted from them. The LSDBinary algorithm is specifically developed to perform detailed spectroscopic studies of eclipsing SB2 systems whose orbital configuration and components' atmospheric parameters are estimated by other means.

Multi-mode Cepheids pulsate simultaneously in more than one mode of oscillation. They provide an independent means to test stellar models and pulsation theories. They can also be used to derive metallicities. In recent years, the number of known multi-mode Cepheids has increased dramatically with the discovery of a large number of Galactic double-mode Cepheids. To date, 209 double-mode Cepheids have been detected in the Galactic bulge and disk, mostly based on the Optical Gravitational Lensing Experiment's (OGLE) catalog. In this paper, we conduct a comprehensive search for double-mode Cepheids in the northern sky based on Zwicky Transient Facility Data Release 5. We found 72 such objects in the Milky Way. The periods of the 30 sample objects already included in the OGLE catalog show excellent agreement with the OGLE periods. The period ratios of our new Cepheids are consistent with those of known double-mode Cepheids, as evidenced by their loci in the so-called `Petersen diagram'. Compared with OGLE, the completeness of our double-mode Cepheid sample is around 71\%. The much improved temporal sampling of the Zwicky Transient Facility offers significant scope to find more double-mode Cepheids, especially at the distribution's short-period end.

We report the discovery of the bright reflected emission component in the super-Eddington state of the ULX pulsar Swift J0243.6+6124, based on the NuSTAR observations of the source during its 2017 outburst. The flux of the reflected emission is weakly variable over the pulsar phase while the direct emission shows significantly larger pulsation amplitude. We propose that in this system the neutron star finds itself in the centre of the well formed by the inner edge of the geometrically thick super-Eddington accretion disc truncated by the magnetic field of the pulsar. The aspect ratio of the well is H/R \sim 1. The inner edge of the truncated disc is continuously illuminated by the emission of the accretion column giving rise to the weakly variable reflected emission. As the neutron star rotates, its emission sweeps through the line of sight, giving rise to the pulsating direct emission. From Doppler broadening of the iron line, we measure the truncation radius of the accretion disc \sim 50 R_g. The inferred dipole component of the magnetic field is consistent with previous estimates favouring a not very strong field. The uniqueness of this system is determined by its moderately super-Eddington accretion rate and the moderate magnetic field so that the inner edge of the truncated geometrically thick accretion disc is seen from the neutron star at a large solid angle.

In order to predict the spins of stellar remnants, we need to understand the evolution of the internal rotation of stars, and to identify at which stage the rotation of the contracting cores of evolved stars decouples from their expanding envelopes. The donor stars of mass transferring binaries lose almost their entire envelope and may thus offer a straight view on their core rotation. While after the mass transfer event they contract and fade rapidly, they are well observable when caught in the short-lived B-star phase. The B-type primary of the galactic binary system LB-1, which was originally suggested to contain a massive black hole, is well explained as a stripped star accompanied by a fainter Be star. The narrow absorption lines in the primary's spectrum signify extremely slow rotation, atypical for B-type main-sequence stars. Here, we investigate the evolution of mass donors in generic grids of detailed binary evolution models, where both stars include differential rotation, internal angular momentum transport and spin-orbit coupling. Whereas the mass gainers are typically spun-up during the mass transfer, we find the spins of the stripped donor models to depend sensitively on the employed mechanism for internal angular momentum transport. Purely hydrodynamic transport can not explain the observed slow rotation, while models including magnetic angular momentum transport are able to reproduce the observed rotation of LB-1 and similar stars, independent of the initial rotation rate. In such models, the spin of the white dwarfs which emerge at the end of the evolution is independent of the mass stripping. We find evidence that the mass transfer in LB-1 was moderately non-conservative.

Using the APEX-12m telescope, continuum maps at 350~$\mu$m of eight gas-dust clouds from the southern hemisphere are obtained. Clouds are associated with the regions of massive star and star cluster formation and have dense cores. Core sizes estimated at the half maximum level are $\sim 0.1-0.2$~pc, core masses and mean gas densities lie in the ranges: $\sim 20-1000~M_{\odot}$ and $\sim (0.3-7.3)\times 10^6$~cm$^{-3}$, respectively. Comparison of the 350~$\mu$m data with the data of observations at 1.2~mm has been performed. From the ratios of intensities at two wavelengths convolved to the same angular resolution, spatial distributions of dust temperatures averaged over the line of sight are calculated. Dust temperature maps for most objects correlate with spatial distributions of intensities at 350~$\mu$m. A decrease in dust temperature with a distance from the center is detected in the cores. Dust temperature profiles in most cases are close to linear ones. Using a simple spherically-symmetric cloud model it is shown that temperature profiles similar to the observed ones can be obtained in a model with internal source by varying density profile parameters and setting $\beta$, a power-law index of the dust emissivity dependence on frequency, constant. It is shown that the dust temperatures estimates strongly depend on the accepted value of $\beta$. It is considered how possible variations of $\beta$ in a cloud can affect the results obtained.

The ASTRI Mini-Array is a next-generation system of nine imaging atmospheric Cherenkov telescopes that is going to be built at the Observatorio del Teide site. After a first phase, in which the instrument will be operated as an experiment prioritizing a schedule of primary science cases, an observatory phase is foreseen in which other significant targets will be pointed. We focus on the observational feasibility of extragalactic sources and on astrophysical processes that best complement and expand the ASTRI Mini-Array core science, presenting the most relevant examples that are at reach of detection over long-term time scales and whose observation can provide breakthrough achievements in the very-high energy extragalactic science. Such examples cover a wide range of $\gamma$-ray emitters, including the study of AGN low states in the multi-TeV energy range, the possible detection of Seyfert galaxies with long exposures and the searches of dark matter lines above 10 TeV. Simulations of the presented objects show that the instrument performance will be competitive at multi-TeV energies with respect to current arrays of Cherenkov telescopes.

The ASTRI (Astrofisica con Specchi a Tecnologia Replicante Italiana) Project led by the Italian National Institute for Astrophysics (INAF) is developing and will deploy at the Observatorio del Teide a mini-array (ASTRI Mini-Array) composed of nine telescopes similar to the small-size dual-mirror Schwarzschild-Couder telescope (ASTRI-Horn) currently operating on the slopes of Mt. Etna in Sicily. The ASTRI Mini-Array will surpass the current Cherenkov telescope array differential sensitivity above a few tera-electronvolt (TeV), extending the energy band well above hundreds of TeV. This will allow us to explore a new window of the electromagnetic spectrum, by convolving the sensitivity performance with excellent angular and energy resolution figures. In this paper we describe the Core Science that we will address during the first four years of operation, providing examples of the breakthrough results that we will obtain when dealing with current open questions, such as the acceleration of cosmic rays, cosmology and fundamental physics and the new window, for the TeV energy band, of the time-domain astrophysics.

The ASTRI (Astrofisica con Specchi a Tecnologia Replicante Italiana) Mini-Array will be composed of nine imaging atmospheric Cherenkov telescopes at the Observatorio del Teide site. The array will be best suited for astrophysical observations in the 0.3-200 TeV range with an angular resolution of few arc-minutes and an energy resolution of 10-15\%. A core-science programme in the first four years will be devoted to a limited number of key targets, addressing the most important open scientific questions in the very-high energy domain. At the same time, thanks to a wide field of view of about 10 degrees, ASTRI Mini-Array will observe many additional field sources, which will constitute the basis for the long-term observatory programme that will eventually cover all the accessible sky. In this paper, we review different astrophysical Galactic environments, e.g. pulsar wind nebulae, supernova remnants, and gamma-ray binaries, and show the results from a set of ASTRI Mini-Array simulations of some of these field sources made to highlight the expected performance of the array (even at large offset angles) and the important additional observatory science that will complement the core-science program.

In addition to the large surveys and catalogs of massive young stellar objects and outflows, dedicated studies of particular sources, in which high-angular observations (mainly at near-IR and mm) are analyzed in depth, are needed to shed light on the processes involved in the formation of massive stars. The galactic source G079.1272+02.2782 (G79), a MYSO at about 1.4 kpc, is an ideal source to carry out this kind of studies. Near-IR integral field spectroscopic observations were carried out using NIFS at Gemini-North. The spectral and angular resolutions, allow us to perform a detailed study of the source and its southern jet, resolving structures with sizes between 200 and 300 au. As a complement, millimeter data retrieved from the JCMT and the IRAM 30m telescope databases were analyzed to study the molecular gas at a larger spatial scale. The analysis of a jet extending southwards shows cork-screw like structures at 2.2 um continuum, strongly suggesting that the jet is precessing. The jet velocity is estimated in 30-43 km/s and it is coming to us along the line of sight. We suggest that the precession may be produced by the gravitational tidal effects generated in a probable binary system, and we estimate a jet precession period of about 10^3 yr, indicating a slow-precessing jet, which is in agreement with the observed helical features. An analysis of H2 lines along the jet allows us to investigate in detail a bow-shock produced by this jet. We find that this bow-shock is indeed generated by a C-type shock and it is observed coming to us, with some inclination angle, along the line of sight. This is confirmed by the analysis of molecular outflows at a larger spatial scale. A brief analysis of several molecular species at millimeter wavelengths indicates a complex chemistry developing at the external layers of the molecular clump in which MYSO G79 is embedded.

We regard the Sun-as-a-star magnetic field (i.e. the mean field) as a filter for the spherical harmonic components of the photospheric field, and calculate the transmission coefficients of this filter. The coefficients for each harmonic, $Y_{l}^{m}$, are listed in three tables according to their dependence on $B_{0}$, the observer's latitude in the star's polar coordinate system. These coefficients are used to interpret the 46-yr sequence of daily mean-field measurements at the Wilcox Solar Observatory. We find that the non-axisymmetric part of the field originates in the $Y_{1}^{1}$, $Y_{2}^{2}$, and a combination of the $Y_{3}^{3}$ and $Y_{3}^{1}$ harmonic components. The axisymmetric part of the field originates in $Y_{2}^{0}$ plus a $B_{0}$-dependent combination of the $Y_{1}^{0}$ and $Y_{3}^{0}$ components. The power spectrum of the field has peaks at frequencies corresponding to the ~27-day synodic equatorial rotation period and its second and third harmonics. Each of these peaks has fine structure on its low-frequency side, indicating magnetic patterns that rotate slowly under the influence of differential rotation and meridional flow. The sidebands of the fundamental mode resolve into peaks corresponding to periods of ~28.5 and ~30 days, which tend to occur at the start of sunspot maximum, whereas the ~27-day period tends to occur toward the end of sunspot maximum. We expect similar rotational sidebands to occur in magnetic observations of other Sun-like stars and to be a useful complement to asteroseismology studies of convection and magnetic fields in those stars.

We have studied the Be/X-ray binary (BeXRB) pulsar eRASSU J050810.4-660653 recently discovered in the Large Magellanic Cloud (LMC). Timing and spectral features of the source have been discussed in detail using NuSTAR \& XMM-Newton observations. Coherent pulsation of the source was detected at $\sim 40.578\;\pm\;0.001$ s using NuSTAR observation. We analyzed pulse profiles of the source in different energy bands using NuSTAR \& XMM-Newton data. The pulse-profile evolved with time but was generally suggestive of a pencil-beam dominated pattern, which combined with the measured luminosity, indicates that the source may be accreting in the sub-critical regime. The pulse fraction follows a linearly increasing trend with photon energy and is anti-correlated with luminosity. In the 1 year interval between the XMM and NuSTAR observations the pulse period shortened by 0.021 s which could be consistent with spin-up or orbital Doppler effect. The average flux of the source in (3-50) keV energy range is found to be $\sim 5.56\;\times\;10^{-11}\;erg\;cm^{-2}\;s^{-1}$ and the corresponding luminosity is $\sim 1.66\;\times\;10^{37}\;erg\;s^{-1}$. The variation of spectral parameters with pulse phase is studied using phase resolved spectroscopy which reveals that the observed photon index becomes harder with increasing flux.

Starlink is a low-Earth orbit (LEO) satellite constellation operated by Space Exploration Technologies Corp. (SpaceX) which aims to provide global satellite internet access. Thus far, most photometric observations of Starlink satellites have primarily been from citizen scientists' visual observations without using quantitative detectors. This paper aims to characterize Starlink satellites and investigate the impact of mega constellations on ground-based astronomy, considering both the observed magnitude and two-line element (TLE) residuals. We collected 353 observations of 61 different Starlink satellites over a 16-month period and we found an average GAIA G magnitude of 5.5 +/- 0.13 with a standard deviation of 1.12. The average magnitude of V1.0 (pre-VisorSat) Starlinks was 5.1 +/- 0.13 with a standard deviation of 1.13. SpaceX briefly used a low-albedo coating on a Starlink satellite called DarkSat to test light pollution mitigation technologies. The brightness of DarkSat was found to be 7.3 +/- 0.13 with a standard deviation of 0.78, or 7.6 times fainter than V1.0 Starlinks. This concept was later abandoned due to thermal control issues and sun visors were used in future models called VisorSats. The brightness of VisorSats was found to be 6.0 +/- 0.13 with a standard deviation of 0.79, or 2.3 times fainter than V1.0 Starlinks. Over the span of the observations, we found that TLEs were accurate to within an average of 0.12 degrees in right ascension and -0.08 degrees in declination. The error is predominantly along-track, corresponding to a 0.3 second time error between the observed and TLE trajectories. Our observations show that a time difference of 0.3 +/- 0.28 seconds is viable for a proposed 10 second shutter closure time to avoid Starlinks in images.

Using time-sequence vector magnetic field and coronal observations from \textit{Solar Dynamics Observatory}, we report the observations of the magnetic field evolution and coronal activity in four emerging active regions (ARs). The ARs emerge with leading polarity being the same as for the majority of ARs in a hemisphere of solar cycle 24. After emergence, the magnetic polarities separate each other without building a sheared polarity inversion line. In all four ARs, the magnetic fields are driven by foot point motions such that the sign of the helicity injection ($dH/dt$) in the first half of the evolution is changed to the opposite sign in the later part of the observation time. This successive injection of opposite helicity is also consistent with the sign of mean force-free twist parameter ($\alpha_{av}$). Further, the EUV light curves of the ARs in 94\AA~and GOES X-ray flux reveal flaring activity below C-class magnitude. Importantly, the white-light coronagraph images in conjunction with the AR images in AIA 94 \AA~delineate the absence of associated CMEs with the studied ARs. These observations imply that the ARs with successive injection of opposite sign magnetic helicity are not favorable to twisted flux rope formation with excess coronal helicity, and therefore are unable to launch CMEs, according to recent reports. This study provides the characteristics of helicity flux evolution in the ARs referring to the conservative property of magnetic helicity and more such studies would help to quantify the eruptive capability of a given AR.

The phase transition between galaxies and quasars is often identified with the rare population of hyper-luminous, hot dust-obscured galaxies. Galaxy formation models predict these systems to grow via mergers, that can deliver large amounts of gas toward their centers, induce intense bursts of star formation and feed their supermassive black holes. Here we report the detection of 24 galaxies emitting Lyman-alpha emission on projected physical scales of about 400 kpc around the hyper-luminous hot dust-obscured galaxy W0410-0913, at redshift z = 3.631, using Very Large Telescope observations. While this indicates that W0410-0913 evolves in a very dense environment, we do not find clear signs of mergers that could sustain its growth. Data suggest that if mergers occurred, as models expect, these would involve less massive satellites, with only a moderate impact on the internal interstellar medium of W0410-0913, which is sustained by a rotationally-supported fast-rotating molecular disk, as Atacama Large Millimeter Array observations suggest.

The paper develops a method for detecting optical binary stars based on the use of astrometric catalogs in combination with machine learning (ML) methods. A computational experiment was carried out on the example of the HIPPARCOS mission catalog and the Pan-STARRS (PS1) catalog by applying the suggested method. It has shown that the reliability of predicting a stellar binarity reaches 90-95%. We note the prospects and effectiveness of creating a proprietary research platform - Cognotron.

Science archives are cornerstones of modern astronomical facilities. In this paper we describe the version 1.0 milestone of the Atacama Large Millimeter/submillimeter Array Science Archive. This version features a comprehensive query interface with rich metadata and visualisation of the spatial and spectral locations of the observations, a complete set of virtual observatory services for programmatic access, text-based similarity search, display and query for types of astronomical objects in SIMBAD and NED, browser-based remote visualisation, interactive previews with tentative line identification and extensive documentation including video and Jupyter Notebook tutorials. The development is regularly evaluated by means of user surveys and is entirely focused on providing the best possible user experience with the goal of helping to maximise the scientific productivity of the observatory.

The recent discovery of a kilonova associated with an apparent long-duration gamma-ray burst has challenged the typical classification that long gamma-ray bursts originate from the core collapse of massive stars and short gamma-ray bursts are from compact binary coalescence. The kilonova indicates a neutron star merger origin and suggests the viability of gravitational-wave and long gamma-ray burst multimessenger astronomy. Gravitational waves play a crucial role by providing independent information for the source properties. This work revisits the archival 2015-2020 LIGO/Virgo gravitational-wave candidates from the 4-OGC catalog which are consistent with a binary neutron star or neutron star-black hole merger and the long-duration gamma-ray bursts from the Fermi and Swift catalogs. We search for spatial and temporal coincidence with up to 10 s time lag between gravitational-wave candidates and the onset of long-duration GRBs. The most significant candidate association has only a false alarm rate of once every two years; given the LIGO/Virgo observational period, this is consistent with a null result. We report an exclusion distance for each search candidate for a fiducial gravitational-wave signal and conservative viewing angle assumptions.

We present an interstellar medium and stellar population analysis of three spectroscopically confirmed $z>7$ galaxies in the ERO JWST/NIRCam and JWST/NIRSpec data of the SMACS J0723.3-7327 cluster. We use the Bayesian spectral energy distribution (SED) fitting code Prospector with a flexible star-formation history (SFH), a variable dust attenuation law, and a self-consistent model of nebular emission (continuum and emission lines). Importantly, we self-consistently fit both the emission line fluxes from JWST/NIRSpec and the broad-band photometry from JWST/NIRCam, taking into account slit-loss effects. We find that these three $z=7.6-8.5$ galaxies ($M_{\star}\approx10^{8}~M_{\odot}$) are young with rising SFHs and mass-weighted ages of 3-7 Myr, though we find indications for underlying older stellar populations. The inferred gas-phase metallicities broadly agree with the direct metallicity estimates from the auroral lines. The galaxy with the lowest gas-phase metallicity ($\mathrm{Z}_{\rm gas}=0.06~\mathrm{Z}_{\odot}$) has a steeply rising SFH, is very compact ($<0.2~\mathrm{kpc}$) and has a high star-formation rate surface density ($\Sigma_{\rm SFR}\approx38~\mathrm{M}_{\odot}~\mathrm{yr}^{-1}~\mathrm{kpc}^{-2}$), consistent with rapid gas accretion. The two other objects with higher gas-phase metallicity show more complex multi-component morphologies on kpc scales, indicating that their recent increase in star-formation rate is driven by mergers or internal, gravitational instabilities. We discuss effects of assuming different SFH priors or only fitting the photometric data. Our analysis highlights the strength and importance of combining JWST imaging and spectroscopy for fully assessing the nature of galaxies at the earliest epochs.

Neutron stars can capture asymmetric dark matter (ADM), which affects the neutron star's measurable properties and makes compact objects prime targets to search for ADM. In this work, we use Bayesian inference to explore potential neutron star mass-radius measurements, from current and future X-ray telescopes, to constrain the bosonic ADM parameters for the case where bosonic ADM has accumulated in the neutron star interior. We find that the high bosonic ADM particle mass ($m_\chi$) and low effective self-interaction strength ($g_\chi/m_\phi)$ regime is disfavored due to the observationally and theoretically motivated constraint that neutron stars must have at least a mass of $1 \, \mathrm{M_\odot}$. However, within the remaining parameter space, $m_\chi$ and $g_\chi/m_\phi$ are individually unconstrained. On the other hand, the ADM mass-fraction, i.e., the fraction of ADM mass inside the neutron star, can be constrained by such neutron star measurements. The inclusion of bosonic ADM in neutron star cores also relaxes the constraints on the baryonic equation of state space and suggests that ADM should be taken into account when interpreting constraints from mass-radius measurements.

Since the end of 2018, the Transiting Exoplanet Survey Satellite (TESS) has provided stellar photometry to the astronomical community. We have used TESS data to study rotational modulation in the light curves of a sample of chemically peculiar stars with measured large-scale magnetic fields (mCP stars). In general, mCP stars show inhomogeneous distributions of elements in their atmospheres that lead to spectroscopic (line profile) and photometric (light curve) variations commensurate with the rotational period. We analyzed the available TESS data from 50 sectors for eight targets after post-processing them in order to minimize systematic instrumental trends. Analysis of the light curves allowed us to determine rotational periods for all eight of our targets. For each star, we provide a phase diagram calculated using the derived period from the light curves and from the available measurements of the disk-averaged longitudinal magnetic field $\langle B_{\rm z}\rangle$. In most cases, the phased light curve and $\langle B_{\rm z}\rangle$ measurements show consistent variability. Using our rotation periods, and global stellar parameters derived from fitting Balmer line profiles, and from Geneva and Str\"omgren-Crawford photometry, we determined the equatorial rotational velocities and calculated the respective critical rotational fractions $v_{\rm eq}/v_{\rm crit}$. We have shown from our sample that the critical rotational fraction decreases with stellar age, at a rate consistent with the magnetic braking observed in the larger population of mCP stars.

Previous studies on astrophysical dark matter (DM) constraints have all assumed that the Milky Way's (MW) DM halo can be modelled in isolation. However, recent work suggests that the MW's largest dwarf satellite, the Large Magellanic Cloud (LMC), has a mass of 10-20$\%$ that of the MW and is currently merging with our Galaxy. As a result, the DM haloes of the MW and LMC are expected to be strongly deformed. We here address and quantify the impact of the dynamical response caused by the passage of the LMC through the MW on the prospects for indirect DM searches. Utilising a set of state-of-the-art numerical simulations of the evolution of the MW-LMC system, we derive the DM distribution in both galaxies at the present time based on the Basis Function Expansion formalism. Consequently, we build $J$-factor all-sky maps of the MW-LMC system in order to study the impact of the LMC passage on gamma-ray indirect searches for thermally produced DM annihilating in the outer MW halo as well as within the LMC halo standalone. We conduct a detailed analysis of 12 years of Fermi-LAT data that incorporates various large-scale gamma-ray emission components and we quantify the systematic uncertainty associated with the imperfect knowledge of the astrophysical gamma-ray sources. We find that the dynamical response caused by the LMC passage can alter the constraints on the velocity-averaged annihilation cross section for weak scale particle DM at a level comparable to the existing observational uncertainty of the MW halo's density profile and total mass.

The Event Horizon Telescope (EHT) has recently observed the image and shadow of the supermassive compact object Sagittarius A$^*$ (Sgr A$^*$). According to the EHT collaboration, the observed image is consistent with the expected appearance of a Kerr black hole. However, it is well-known that some non-Kerr objects may mimic many of the properties of the Kerr black hole, and hence, their shadows might be consistent with the observed shadow of Sgr A$^*$. In this work, we consider two black hole mimickers and study their shadows. The first mimicker is a rotating generalisation of the recently proposed static, spherically symmetric black-bounce spacetime by Simpson and Visser where the central Schwarzschild singularity is replaced by a minimal surface. The second one is the $\gamma$-metric which is a static, axially symmetric singular solution of the vacuum Einstein's equations without an event horizon. We put constraint on the parameters of these black hole mimickers by comparing their shadows with the one observed for Sgr A$^*$.

Shift-symmetric Horndeski theories admit an interesting class of Schwarzschild-de Sitter black hole solutions exhibiting time-dependent scalar hair. The properties of these solutions may be studied via a bottom-up effective field theory (EFT) based on the background symmetries. This is in part possible by making use of a convenient coordinate choice -- Lema\^itre-type coordinates -- in which the profile of the Horndeski scalar field is linear in the relevant time coordinate. We construct this EFT, and use it to understand the stability of hairy black holes in shift-symmetric Horndeski theories, providing a set of constraints that the otherwise-free functions appearing in the Horndeski Lagrangian must satisfy in order to admit stable black hole solutions. The EFT is analyzed in the decoupling limit to understand potential sources of instability. We also perform a complete analysis of the EFT with odd-parity linear perturbations around general spherically symmetric space-time.

A number of modern millimeter, sub-millimeter, and far-infrared detectors are read out using superconducting microwave (1-10GHz) resonators. The main detector technologies are Transition Edge Sensors, read out using Microwave SQUID Multiplexers ($\mu$mux) and Microwave Kinetic Inductance Detectors. In these readout schemes, sky signal is encoded as resonance frequency changes. One way to interrogate these superconducting resonators is to calibrate the probe tone phase such that any sky signal induced frequency shifts from the resonators show up primarily as voltage changes in only one of the two quadratures of the interrogation tone. However, temperature variations in the operating environment produce phase drifts that degrade the phase calibration and can source low frequency noise in the final detector time ordered data if left to drift too far from optimal calibration. We present a method for active software monitoring of the time delay through the system which could be used to feedback on the resonator probe tone calibration angle or to apply an offline cleaning. We implement and demonstrate this monitoring method using the SLAC Microresonator RF Electronics on a 65 channel $\mu$mux chip from NIST.

We find a significant destructive interference among Kerr overtones in the early ringdown induced by an extreme mass-ratio merger of a massive black hole and a compact object, and that the ringdown spectrum apparently follows the Fermi-Dirac distribution. We numerically compute the spectral amplitude of gravitational waves induced by a particle plunging into a Kerr black hole and study the excitation of multiple quasi-normal (QN) modes. We find that the start time of ringdown is before the strain peak of the signal and corresponds to the time when the particle passes the photon sphere. When the black hole has the near-extremal rotation, the Kerr QN frequencies are close to the fermionic Matsubara frequencies with the Hawking temperature and a chemical potential of the superradiant frequency. We indeed find that the absolute square of the spectral amplitude apparently follows the Fermi-Dirac distribution with the chemical potential of around the real QN frequency of the fundamental mode. Fitting the Boltzmann distribution to the data in higher frequencies, the best-fit temperature is found out to be close to the Hawking temperature, especially for rapid rotations. In the near-extremal limit, the gravitational-wave spectrum exhibits a would-be Fermi degeneracy with the Fermi surface at the superradiant frequency $\omega = \mu_{\rm H}$. This opens a new possibility that we can test the holographic nature of Kerr black holes, such as the Kerr/CFT correspondence, by observationally searching for the Boltzmann distribution in frequencies higher than $\mu_{\rm H}$ without extracting overtones out of ringdown.

Information on the phase structure of strongly interacting matter at high baryon densities can be gained from observations of neutron stars and their detailed analysis. In the present work Bayesian inference methods are used to set constraints on the speed of sound in the interior of neutron stars, based on recent multi-messenger data in combination with limiting conditions from nuclear physics at low densities. Two general parametric representations are introduced for the sound speed $c_s$ in order to examine the independence with respect to choices for the parametrisation of Priors. Credible regions for neutron star properties are analysed, in particular with reference to the quest for possible phase transitions in cold dense matter. The evaluation of Bayes factors implies extreme evidence for a violation of the conformal bound, $c_s^2 \leq 1/3$, inside neutron stars. Given the presently existing data base, it can be concluded that the occurrence of a first-order phase transition in the core of even a two-solar-mass neutron star is unlikely, while a continuous crossover cannot be ruled out. At the same time it is pointed out that the discovery of a superheavy neutron star with a mass $ M \sim 2.3 - 2.4\, M_\odot$ would strengthen evidence for a phase change in the deep interior of the star.

We argue that cubic order interactions between two scalar gravitons and one tensor graviton are ubiquitous in models of dark energy where the strong coupling scale is $\Lambda_3$. These interactions can potentially provide efficient decay channels for gravitational waves. They can also lead to gradient instabilities of the scalar perturbations in the presence of large amplitude gravitational waves, e.g. those detected by LIGO. In contrast with models in scalar-tensor theories, there is an infinite number of higher order interactions in generic $\Lambda_3$ models, which make it difficult to predict the fate of these instabilities inferred from cubic order interactions.

I discuss why state-of-the art perturbative QCD calculations of the equation of state at large chemical potential that are reliable at asymptotically high densities constrain the same equation of state at neutron-star densities. I describe how these theoretical calculations affect the EOS at lower density. I argue that the ab-initio calculations in QCD offer significant information about the equation of state of the neutron-star matter, which is complementary to the current astrophysical observations.

A quantitative description of the properties of hot nuclear matter will be needed for the interpretation of the available and forthcoming astrophysical data, providing information on the post merger phase of a neutron star coalescence. We have employed a recently developed theoretical model, based on a phenomenological nuclear Hamiltonian including two- and three-nucleon potentials, to study the temperature dependence of average and single-particle properties of nuclear matter relevant to astrophysical applications. The possibility to represent the results of microscopic calculations using simple and yet physically motivated parametrisations of thermal effects, suitable for use in numerical simulations of astrophysical processes, is also discussed.

Wavelike, bosonic dark matter candidates like axion and dark photons can be detected using microwave cavities commonly referred to as haloscopes. Traditionally, haloscopes consist of tunable copper cavities operating in the TM$_{010}$ mode, but the performance of these cavities has been limited by ohmic losses. In contrast, superconducting radio frequency (SRF) cavities can achieve quality factors of $\sim 10^{10}$, perhaps five orders of magnitude better than copper cavities, which would lead to more sensitive dark matter detectors. In this paper, we first derive that the scan rate of a haloscope experiment is proportional to the loaded quality factor $Q_L$, even if the cavity bandwidth is much narrower than the dark matter halo lineshape. We then present a proof-of-concept search for dark photon dark matter using a non-tunable ultra-high quality SRF cavity. We exclude dark photon dark matter with kinetic mixing strengths of $\chi > 2\times 10^{-16}$ for a dark photon mass of $m_{A^{\prime}} = 5.37 \mu\text{eV}$, achieving the deepest exclusion to wavelike dark photons by almost an order of magnitude.

Using a new sample of extremely metal poor systems, the EMPRESS survey has recently reported a primordial helium abundance that is $3\sigma$ smaller than the prediction from the Standard BBN scenario. This measurement could be interpreted as a hint for a primordial lepton asymmetry in the electron neutrino flavor. Motivated by the EMPRESS results, we present a comprehensive analysis of the lepton asymmetry using measurements of the abundances of primordial elements, along with CMB data from Planck. Assuming that there is no dark radiation in our Universe, we find an electron neutrino chemical potential $\xi_{\nu_e} = 0.037 \pm 0.013$, which deviates from zero by $2.8\sigma$. If no assumption is made on the abundance of dark radiation in the Universe, the chemical potential is $\xi_{\nu_e} = 0.037 \pm 0.020$, which deviates from zero by $1.9\sigma$. We also find that this result is rather insensitive to the choice of nuclear reaction rates. If the true helium abundance corresponds to the EMPRESS central value, future CMB observations from the Simons Observatory and CMB-S4 will increase the significance for a non-zero lepton asymmetry to $4\sigma$ and $5\sigma$ respectively, assuming no dark radiation, or to $3\sigma$ when no assumption is made on the abundance of dark radiation.

A chiral chemical potential present in the early universe can source helical hypermagnetic fields through the chiral plasma instability. If these hypermagnetic fields survive until the electroweak phase transition, they source a contribution to the baryon asymmetry of the universe. In this letter, we demonstrate that lepton flavour asymmetries above $|\mu|/T \sim 4 \times 10^{-3}$ trigger this mechanism even for vanishing total lepton number. This excludes the possibility of such large lepton flavour asymmetries present at temperatures above $10^6$ GeV, setting a constraint which is about two orders of magnitude stronger than the current CMB and BBN limits.