Papers for Monday, Sep 05 2022

The Warm-Hot Intergalactic Medium (WHIM) is believed to host a significant fraction of the ``missing baryons'' in the nearby Universe. Its signature has been detected in the X-ray absorption spectra of distant quasars. However, its detection in emission, that would allow us to study the WHIM in a systematic way, is still lacking. Motivated by the possibility to perform these studies with next generation integral field spectrometers, and thanks to the availability of a large suite of state-of-the-art hydrodynamic simulations -- the CAMELS suite -- we study here in detail the emission properties of the WHIM and the possibility to infer its physical properties with upcoming X-ray missions like Athena. We focused on the two most prominent WHIM emission lines, the OVII triplet and the OVIII singlet, and build line surface brightness maps in a lightcone, mimicking a data cube generated through integral field spectroscopy. We confirm that detectable WHIM emission, even with next generation instruments, is largely associated to galaxy-size dark matter halos and that the WHIM properties evolve little from $z\simeq0.5$ to now. Some characteristics of the WHIM, like the line number counts as a function of their brightness, depend on the specific hydrodynamic simulation used, while others, like the WHIM clustering properties, are robust to this aspect. The large number of simulations available in the CAMELS datasets allows us to assess the sensitivity of the WHIM properties to the background cosmology and to the energy feedback mechanisms regulated by AGN and stellar activity. [ABRIDGED]

High-resolution observations of the Perseus B5 "core" have previously revealed that this subsonic region actually consists of several filaments that are likely in the process of forming a quadruple stellar system. Since subsonic filaments are thought to be produced at the $\sim 0.1$ pc sonic scale by turbulent compression, a detailed kinematic study is crucial to test such a scenario in the context of core and star formation. Here we present a detailed kinematic follow-up study of the B5 filaments at a 0.009 pc resolution using the VLA and GBT combined observations fitted with multi-component spectral models. Using precisely identified filament spines, we find a remarkable resemblance between the averaged width profiles of each filament and Plummer-like functions, with filaments possessing FWHM widths of $\sim 0.03$ pc. The velocity dispersion profiles of the filaments also show decreasing trends towards the filament spines. Moreover, the velocity gradient field in B5 appears to be locally well ordered ($\sim 0.04$ pc) but globally complex, with kinematic behaviors suggestive of inhomogeneous turbulent accretion onto filaments and longitudinal flows towards a local overdensity along one of the filaments.

We present a new prediction of the luminosity functions of the [CII] line at 158 $\mu$m, of the CO lines from J=0 to J=24, and of the molecular gas mass density up to z=10, using the Spectro-Photometric Realisations of Infrared-selected Targets at all-z (SPRITZ) simulation (Bisigello et al. 2021). We update the state-of-the-art phenomenological simulation SPRITZ to include both the CO ($J\leq24$) and [CII] line luminosities. This has been performed using different empirical and theoretical relations to convert the total infrared luminosity (or star formation rate) to [CII] or CO luminosity. The resulting line luminosity functions have been compared for validation with a large set of observations available in the literature. We then used the derived CO and [CII] line luminosities to estimate the molecular gas mass density and compare it with available observations. The CO and [CII] luminosity functions presented here are well in agreement with all the available observations. In particular, the best results for [CII] are obtained deriving the [CII] luminosity directly from the star formation rate, but considering a dependence of this relation on the gas metallicity. For all the CO luminosity functions, the estimates favoured by the data are derived considering different relations, depending on the ionisation mechanism dominating each galaxy, i.e. star formation or active galactic nuclei, and, moreover, deriving the $J\geq4$ CO lines directly from the [CII] luminosity. However, further data are necessary to fully discriminate between models. Finally, the best agreement with observations of the molecular gas mass density are derived by converting the [CII] luminosity to H2 mass, using a [CII]-to-H2 conversion ~130 $\rm M_{\odot}/{\rm L}_{\odot}$. All the line luminosity functions, useful for planning and interpreting future observations, are made publicly available.

The most massive stars dominate the chemical enrichment, mechanical and radiative feedback, and energy budget of their host environments. Yet how massive stars initially form and how they evolve throughout their lives is ambiguous. The mass loss of the most massive stars remains a key unknown in stellar physics, with consequences for stellar feedback and populations. In this work, we compare grids of very massive star (VMS) models with masses ranging from 80-1000Msun, for a range of input physics. We include enhanced winds close to the Eddington limit as a comparison to standard O-star winds, with consequences for present-day observations of ~50-100Msun stars. We probe the relevant surface H abundances (Xs) to determine the key traits of VMS evolution compared to O stars. We find fundamental differences in the behaviour of our models with the enhanced-wind prescription, with a convergence on the stellar mass at 1.6 Myr, regardless of the initial mass. It turns out that Xs is an important tool in deciphering the initial mass due to the chemically homogeneous nature of VMS above a mass threshold. We use Xs to break the degeneracy of the initial masses of both components of a detached binary, and a sample of WNh stars in the Tarantula nebula. We find that for some objects, the initial masses are unrestricted and, as such, even initial masses of the order 1000Msun are not excluded. Coupled with the mass turnover at 1.6 Myr, Xs can be used as a 'clock' to determine the upper stellar mass.

Constraints on cosmological parameters are often distilled from sky surveys by fitting templates to summary statistics of the data that are motivated by a fiducial cosmological model. However, recent work has shown how to estimate the distance scale using templates that are more generic: the basis functions used are not explicitly tied to any one cosmological model. We describe a Bayesian framework for (i) determining how many basis functions to use and (ii) comparing one basis set with another. Our formulation provides intuition into how (a) one's degree of belief in different basis sets, (b) the fact that the choice of priors depends on basis set, and (c) the data set itself, together determine the derived constraints. We illustrate our framework using measurements in simulated datasets before applying it to real data.

We started a survey with CHARA/MIRC-X and VLTI/GRAVITY to search for low mass companions orbiting individual components of intermediate mass binary systems. With the incredible precision of these instruments, we can detect astrometric "wobbles" from companions down to a few tens of micro-arcseconds. This allows us to detect any previously unseen triple systems in our list of binaries. We present the orbits of 12 companions around early F to B-type binaries, 9 of which are new detections and 3 of which are first astrometric detections of known RV companions. The masses of these newly detected components range from 0.45-1.3 solar masses. Our orbits constrain these systems to a high astrometric precision, with median residuals to the orbital fit of 20-50 micro-arcseconds in most cases. For 7 of these systems we include newly obtained radial velocity data, which help us to identify the system configuration and to solve for masses of individual components in some cases. Although additional RV measurements are needed to break degeneracy in the mutual inclination, we find that the majority of these inner triples are not well-aligned with the wide binary orbit. This hints that higher mass triples are more misaligned compared to solar and lower mass triples, though a thorough study of survey biases is needed. We show that the ARMADA survey is extremely successful at uncovering previously unseen companions in binaries. This method will be used in upcoming papers to constrain companion demographics in intermediate mass binary systems down to the planetary mass regime.

Weakly accreting black hole X-ray binaries launch compact radio jets that persist even in the quiescent spectral state, at X-ray luminosities <1e-5 of the Eddington luminosity. However, radio continuum emission has been detected from only a few of these quiescent systems, and little is known about their radio variability. Jet variability can lead to misclassification of accreting compact objects in quiescence, and affects the detectability of black hole X-ray binaries in next-generation radio surveys. Here we present the results of a radio monitoring campaign of A0620-00, one of the best-studied and least-luminous known quiescent black hole X-ray binaries. We observed A0620-00 at 9.8 GHz using the Karl G. Jansky Very Large Array on 31 epochs from 2017 to 2020, detecting the source ~75% of the time. We see significant variability over all timescales sampled, and the observed flux densities follow a lognormal distribution with a mean of 12.5 uJy and standard deviation of 0.22 dex. In no epoch was A0620-00 as bright as in 2005 (51 +/- 7 uJy), implying either that this original detection was obtained during an unusually bright flare, or that the system is fading in the radio over time. We present tentative evidence that the quiescent radio emission from A0620-00 is less variable than that of V404 Cyg, the only other black hole binary with comparable data. Given that V404 Cyg has a jet radio luminosity ~20 times higher than A0620-00, this comparison could suggest that less luminous jets are less variable in quiescence.

The polarimetry of gamma rays converting to an $e^+e^-$ pair would open a new window on the high-energy gamma-ray sky by, among other things, providing insight into the radiation mechanism in pulsars (curvature or synchrotron) or deciphering the composition of the gamma-ray emitting jets in blazars (leptonic or lepto-hadronic). The performance of polarimeters based on homogeneous active targets (gas detectors (MeV, HARPO) or emulsions (GeV, GRAINE) has been studied both with simulation and by the analysis of data collected with telescope prototypes on linearly-polarised gamma-ray beams, and found to be excellent. The present (Fermi LAT), AGILE and future project (AMEGO, ASTROGAM) gamma-ray missions, though, are using active targets based on silicon strip detectors (SSD). No demonstration of a non-zero effective polarisation asymmetry with SSDs has been published to date, be it only with simulated data, and sensitivity estimations were obtained from an assumed value of the effective polarisation asymmetry. I present a characterisation of the potential of SSD-based active targets for polarimetry with gamma-ray conversions to pairs and the development of various methods to improve on the sensitivity. This work could pave the way to providing the polarimetry of the brightest gamma-ray sources of the sky from the decade of data collected by the Fermi LAT and by AGILE, and to guiding the design of future missions.

Advances in computing power and numerical methodologies over the past several decades sparked a prolific output of dynamical investigations of the late stages of terrestrial planet formation. Among other peculiar inner solar system qualities, the ability of simulations to reproduce the small mass of Mars within the planets' geochemically inferred accretion timescale of <10 Myr after the appearance of calcium aluminum-rich inclusions (CAIs) is arguably considered the gold standard for judging evolutionary hypotheses. At present, a number of independent models are capable of consistently generating Mars-like planets and simultaneously satisfying various important observational and geochemical constraints. However, all models must still account for the effects of the epoch of giant planet migration and orbital instability; an event which dynamical and cosmochemical constraints indicate occurred within the first 100 Myr after nebular gas dispersal. If the instability occurred in the first few Myr of this window, the disturbance might have affected the bulk of Mars' growth. In this manuscript, we turn our attention to a scenario where the instability took place after t=50 Myr. Specifically, we simulate the instability's effects on three nearly-assembled terrestrial systems that were generated via previous embryo accretion models and contain three large proto-planets with orbits interior to a collection of ~Mars-mass embryos and debris. While the instability consistently triggers a Moon-forming impact and efficiently removes excessive material from the Mars-region in our models, we find that our final systems are too dynamically excited and devoid of Mars and Mercury analogs. Thus, we conclude that, while possible, our scenario is far more improbable than one where the instability either occurred earlier, or at a time where Earth and Venus' orbits were far less dynamically excited.

Close binary systems consisting of two neutron stars (BNS) emit gravitational waves, that allow them to merge on timescales shorter than Hubble time. It is widely believed, that NS-NS mergers in such systems power short gamma-ray bursts (GRB). Several mechanisms which could lead to electromagnetic energy release prior to a merger have been proposed. We estimate the ability to observe the possible pre-burst emission with telescopes of Spectrum-Roentgen-Gamma. We also investigate first such event, GRB210919A, which fell into the field of view of the SRG telescopes less than two days before the burst.

Astrophysical gases such as the interstellar-, circumgalactic- or intracluster-medium are commonly multiphase, which poses the question of the structure of these systems. While there are many known processes leading to fragmentation of cold gas embedded in a (turbulent) hot medium, in this work, we focus on the reverse process: coagulation. This is often seen in wind-tunnel and shearing layer simulations, where cold gas fragments spontaneously coalesce. Using 2D and 3D hydrodynamical simulations, we find that sufficiently large ($\gg c_{\rm s} t_{\rm cool}$), perturbed cold gas clouds develop overstable sound waves which ensure cold gas mass growth over an extended period of time ($\gg r / c_{\rm s}$). This mass growth efficiently accelerates hot gas which in turn can entrain cold droplets, leading to coagulation. The attractive inverse square force between cold gas droplets has interesting parallels with gravity; the `monopole' is surface area rather than mass. We develop a simple analytic model which reproduces our numerical findings.

A statistical model for the polarization of pulsar radio emission is enhanced to account for the heavy modulation of the emission, the possible covariance of the Stokes parameters, and the observed asymmetries in the distributions of total intensity, polarization, and fractional polarization by treating the intensities of the orthogonal polarization modes as exponential random variables. The model is used to derive theoretical distributions to compare with what is observed. The resulting distributions are unimodal and generally asymmetric. The unimodality arises from the model's fundamental assumption that the orthogonal modes are superposed. The asymmetry originates primarily from different fluctuations in mode intensities. The distributions of fractional polarization are truncated at the degree of linear and circular polarization intrinsic to the modes. A number of observable parameters that quantify the statistical properties of the emission and its polarization are derived and are shown to be functions only of the ratio of the modes' mean intensities, M, suggesting their spectra coevolve according to the frequency dependence of M. This particular implementation of the model requires the modes to fluctuate differently in order to replicate the observations. Since a single underlying emission mechanism seems unlikely to selectively modulate the mode intensities, the different fluctuations are attributed either to different emission mechanisms for the modes or to mode-dependent propagation or scattering effects in the pulsar magnetosphere.

Presently large systematic uncertainties remain in the description of hadronic interactions at ultra-high energies and a fully consistent description of air-shower experimental data is yet to be reached. The amount of data collected by the Pierre Auger Observatory using simultaneously the fluorescence and surface detectors in the energy range $10^{18.5}-10^{19.0}$ eV has provided opportunity to perform a multi-parameter test of model predictions. We apply a global method to simultaneously fit the mass composition of cosmic rays and adjustments to the simulated depth of shower maximum ($X_\text{max}$), and hadronic signals at ground level ($R_\text{Had}$). The best description of hybrid data is obtained for a deeper scale of simulated $X_\text{max}$ than predicted by hadronic interaction models tuned to the LHC data. Consequently, the deficit of the simulated hadronic signal at ground level, dominated by muons, is alleviated with respect to the unmodified hadronic interaction models. Because of the size of the adjustments $\Delta X_\text{max}$ and $R_\text{Had}$ and the large number of events in the sample, the statistical significance of these assumed adjustments is large, greater than 5$\sigma_\text{stat}$, even for the combination of the systematic experimental shifts within 1$\sigma_\text{sys}$ that are the most favorable for the models.

Context. The space weathering timescale of near-Earth S-type asteroids has been investigated by several approaches (i.e., experiments, sample analyses, and theoretical approaches), yet there are orders of magnitude differences. Aims. We aim to examine the space weathering timescale on a near-Earth S-type asteroid, Itokawa using Hayabusa-AMICA images and further investigate the evolutional process of the asteroid. Methods. We focused on bright mottles on the boulder surfaces generated via impacts with interplanetary dust particles (IDPs). We compared the bright mottle size distribution with an IDP flux model to determine the space weathering timescale. Results. As a result, we found that the space weathering timescale on Itokawa's boulder surfaces is 10$^3$ years (in the range of 10$^2$-10$^4$ years), which is consistent with the timescale of space weathering by light ions from the solar wind. Conclusions. From this result, we conclude that Itokawa's surface has been weathered shortly in 10$^3$ years but portions of the surface are exposed via seismic shaking triggered by a recent impact that created the Kamoi crater.

We report the serendipitous discovery of an elliptical shell of CO associated with the faint stellar object SSTc2d J163134.1-24006 as part of the "Ophiuchus Disk Survey Employing ALMA" (ODISEA), a project aiming to study the entire population of protoplanetary disks in the Ophiuchus Molecular Cloud from 230 GHz continuum emission and $^{12}$CO (J=2-1), $^{13}$CO (J=2-1) and C$^{18}$CO (J=2-1) lines readable in Band-6. Remarkably, we detect a bright $^{12}$CO elliptical shape emission of $\sim$ 3$^{"}$ $\times$ 4$^{"}$ towards SSTc2d J163134.1-24006 without a 230 GHz continuum detection. Based on the observed near-IR spectrum taken with the Very Large Telescope (KMOS), the brightness of the source, its 3-dimensional motion, and Galactic dynamic arguments, we conclude that the source is not a giant star in the distant background ($>$5 - 10 kpc) and is most likely to be a young brown dwarf in the Ophiuchus cloud, at a distance of just $\sim$139 pc. This is the first report of quasi-spherical mass loss in a young brown dwarf. We suggest that the observed shell could be associated with a thermal pulse produced by the fusion of deuterium, which is not yet well understood, but for a sub-stellar object is expected to occur during a short period of time at an age of a few Myr, in agreement with the ages of the objects in the region. Other more exotic scenarios, such as a merger with planetary companions, cannot be ruled out from the current observations.

Distance-redshift diagrams probe expansion history of the Universe. We show that the stellar mass-binding energy (massE) relation of galaxies proposed in our previous study offers a new distance ruler at cosmic scales. By using elliptical galaxies in the main galaxy sample of the Sloan Digital Sky Survey Data Release 7, we construct a distance-redshift diagram over the redshift range from 0.05 to 0.2 with the massE ruler. The best-fit dark energy density is 0.675+-0.079 for flat Lambda-CDM, consistent with those by other probes. At the median redshift of 0.11, the median distance is estimated to have a fractional error of 0.34%, much lower than those by supernova (SN) Ia and baryonic acoustic oscillation (BAO) and even exceeding their future capability at this redshift. The above low-z measurement is useful for probing dark energy that dominates at the late Universe. For a flat dark energy equation of state model (flat wCDM), the massE alone constrains w to an error that is only a factor of 2.2, 1.7 and 1.3 times larger than those by BAO, SN Ia, and cosmic microwave background (CMB), respectively.

We study the impact of local environment on the transformation of spiral galaxies in three nearby ($z < 0.08$) Abell clusters: A85/A496/A2670. These systems were observed in HI with the Very Large Array, covering a volume extending beyond the virial radius and detecting 10, 58, 38 galaxies, respectively. High fractions (0.40--0.86) of bright spirals [log$(M_{*}/M_{\odot})=9-10$] are not detected in HI. We provide further evidence of environmental effects consisting in significant fractions (0.10--0.33) of abnormal objects and a number of red (passive) spirals, suggesting an ongoing process of quenching. Ram-pressure profiles, and the sample of the brightest spirals used as test particles for environmental effects, indicate that ram-pressure plays an important role in stripping and transforming late-types. Phase-space diagrams and our search for substructures helped to trace the dynamical stage of the three systems. This was used to compare the global cluster effects $vs.$ pre-processing, finding that the former is the dominating mechanism in the studied clusters. By contrasting the global distribution of HI normal $vs.$ HI disturbed spirals in the combined three clusters, we confirm the expected correlation of disturbed objects located, on average, at shorter projected radii. However, individual clusters do not necessarily follow this trend and we show that A496 and A2670 present an atypical behavior. In general we provide conclusive evidence about the dependence of the transformation of infalling spirals on the ensemble of cluster properties like mass, ICM density, dynamical stage and surrounding large-scale structure.

Time-domain (TD) spectroscopic data from the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) can provide accurate and high-cadence radial velocities (RVs). In this work, we search for binaries with compact components with RV monitoring method by using the LAMOST TD survey of four $K$2 plates. Three binary systems including an unseen white dwarf or neutron star are found. For each binary system, we estimate the stellar parameters of the visible star and orbital parameters, and finally calculate the binary mass function and the minimum mass of the unseen star. No obvious double-lined feature is seen from the LAMOST medium-resolution spectra of the three sources. In addition, we found no X-ray counterpart for all these sources but UV companions for two of them. Spectral disentangling also shows no additional component with optical absorption spectra, supporting that these systems contain compact objects.

The Sun's proximity offers us a unique opportunity to study in detail the physical processes on a star's surface; however, the highly dynamic nature of the stellar surface -- in particular, energetic eruptions such as flares and coronal mass ejections -- presents tremendous observational challenges. Spectroscopy probes the physical state of the solar atmosphere, but conventional scanning spectrographs and spectrometers are unable to capture the full evolutionary history of these dynamic events with a sufficiently wide field of view and high spatial, spectral, and temporal resolution. Resolving the physics of the dynamic sun requires gathering simultaneous spectra across a contiguous area over the full duration of these events, a goal now tantalizingly close to achievable with continued investment in developing powerful new Integral Field Spectrographs to serve as the foundation of both future ground- and space-based missions. This technology promises to revolutionize our ability to study solar flares and CMEs, addressing NASA's strategic objective to "understand the Sun, solar system, and universe." Since such events generate electromagnetic radiation and high-energy particles that disrupt terrestrial electric infrastructure, this investment not only advances humanity's scientific endeavors but also enhances our space weather forecasting capability to protect against threats to our technology-reliant civilization.

It is commonly assumed that the stochastic background of gravitational waves on cosmological scales follows an almost scale-independent power spectrum, as generically predicted by the inflationary paradigm. However, it is not inconceivable that the spectrum could have strongly scale-dependent features, generated, e.g., via transient dynamics of spectator axion-gauge fields during inflation. Using the temperature and polarisation maps from the \textit{Planck} and BICEP/Keck datasets, we search for such features, taking the example of a log-normal bump in the primordial tensor spectrum at CMB scales. We do not find any evidence for the existence of bump-like tensor features at present, but demonstrate that future CMB experiments such as LiteBIRD and CMB-S4 will greatly improve our prospects of determining the amplitude, location and width of such a bump. We also highlight the role of delensing in constraining these features at angular scales $\ell\gtrsim 100$.

In preparation for the release of the astrometric orbits of Gaia, Shahaf et al. (2019) proposed a triage technique to identify astrometric binaries with compact companions based on their astrometric semi-major axis, parallax, and primary mass. The technique requires the knowledge of the appropriate mass-luminosity relation to rule out single or close-binary main-sequence companions. The recent publication of the Gaia DR3 astrometric orbits used a schematic version of this approach, identifying 735 astrometric binaries that might have compact companions. In this communication, we return to the triage of the DR3 astrometric binaries with more careful analysis, estimating the probability for its astrometric secondary to be a compact object or a main-sequence close binary. We compile a sample of 177 systems with highly-probable non-luminous massive companions, which is smaller but cleaner than the sample reported in Gaia DR3. The new sample includes 8 candidates to be black-hole systems with compact-object masses larger than 2.4 $M_\odot$. The orbital-eccentricity$-$secondary-mass diagram of the other 169 systems suggests a tentative separation between the white-dwarf and the neutron-star binaries. Most white-dwarf binaries are characterized by small eccentricities of about 0.1 and masses of 0.6 $M_\odot$, while the neutron star binaries display typical eccentricities of 0.4 and masses of 1.3 $M_\odot$.

We investigate the impact of radiation from primordial black holes (PBHs), in the mass range of $10^{15} \lesssim M_{\rm PBH} \lesssim 10^{17}~\rm g$ and $10^{2} \lesssim M_{\rm PBH} \lesssim 10^{4}~M_{\odot}$, on the 21-cm angular-power spectrum in the dark ages. PBHs in the former mass range effect the 21-cm angular-power spectrum through the evaporation known as Hawking radiation, while the radiation from the accretion process in the latter mass range. In the dark ages, radiation from PBHs can increase the ionization fraction and temperature of the intergalactic medium, change the global 21-cm differential brightness temperature and then effect the 21-cm angular-power spectrum. Taking into account the effects of PBHs, we find that in the dark ages, $30 \lesssim z \lesssim 100$, the amplitude of the 21-cm angular-power spectrum is decreased depending on the mass and mass fraction of PBHs. We also investigate the potential constraints on the mass fraction of PBHs in the form of dark matter for the future radio telescope in lunar obit or on the farside surface of the Moon.

Context. After landing on C-type asteroid Ryugu, MASCOT imaged brightly colored, submillimeter-sized inclusions in a small rock. Hayabusa2 successfully returned a sample of small particles from the surface of Ryugu, but none of these appear to harbor such inclusions. The samples are considered representative of Ryugu. Aims. To understand the apparent discrepancy between MASCOT observations and Ryugu samples, we assess whether the MASCOT landing site, and the rock by implication, is perhaps atypical for Ryugu. Methods. We analyzed observations of the MASCOT landing area acquired by three instruments on board Hayabusa2: a camera (ONC), a near-infrared spectrometer (NIRS3), and a thermal infrared imager (TIR). We compared the landing area properties thus retrieved with those of the average Ryugu surface. Results. We selected several areas and landforms in the landing area for analysis: a small crater, a collection of smooth rocks, and the landing site itself. The crater is relatively blue and the rocks are relatively red. The spectral and thermophysical properties of the landing site are very close to those of the average Ryugu surface. The spectral properties of the MASCOT rock are probably close to average, but its thermal inertia may be somewhat higher. Conclusions. The MASCOT rock can also be considered representative of Ryugu. Some of the submillimeter-sized particles in the returned samples stand out because of their atypical spectral properties. Such particles may be present as inclusions in the MASCOT rock.

We study an intermediate-age open cluster NGC 2506 using the \textit{ASTROSAT}/UVIT data and other archival data. We identified 2175 cluster members using a machine learning-based algorithm, ML--MOC, on Gaia EDR3 data. Among the cluster members detected in UVIT filters, F148W, F154W, and F169M, we detect 9 blue straggler stars (BSS), 3 yellow straggler stars (YSS) and 3 red clump (RC) stars. We construct multi-wavelength spectral energy distributions (SEDs) of these objects to characterize them and to estimate their parameters. We discovered hot companions to 3 BSS, 2 YSS and 3 RC candidates and estimated their properties. The hot companions with estimated temperatures, T$\mathrm{_{eff}}$ $\sim$ 13250--31000 K, are WDs of extremely low-mass ($\sim$ 0.20 M$_\odot$), low-mass ($\sim$ 0.20--0.40 M$_\odot$), normal mass ($\sim$ 0.40--0.60 M$_\odot$), and high-mass ($\sim$ 0.8 M$_\odot$). We suggest that systems with extremely low mass and low mass WDs as companions are formed via Case-A/Case-B mass transfer mechanism. A BSS is the likely progenitor of the high mass WD, as a star with more than the turn-off mass of the cluster is needed to form a high mass WD. Thus, systems with high mass WD are likely to be formed through merger in triple systems. We conclude that mass transfer as well as merger pathways of BSS formation are present in this cluster.

Recently, there has been an increasing demand for deep imaging surveys to investigate the history of the mass assembly of galaxies in detail by examining the remnants of mergers and accretions, both of which have very low surface brightness (LSB). In addition, the nature of star formation in LSB regions, such as galaxy outer disks, is also an intriguing topic in terms of understanding the physical mechanisms of disk evolution. To address these issues, this study conducts a survey project, called the Korea Microlensing Telescope Network (KMTNet) Nearby Galaxy Survey to construct a deep imaging data set of nearby galaxies in the southern hemisphere using KMTNet. It provides deep and wide-field images with a field-of-view of $\sim$12 deg$^2$ for 13 nearby galaxies drawn from the Carnegie-Irvine Galaxy Survey catalog, in optical broadbands ($BRI$) and an H$\alpha$ narrowband. Through a dedicated data reduction, the surface brightness limit in 10$^{\prime\prime}\times10^{\prime\prime}$ boxes was found to reach as deep as $\mu_{1\sigma}\sim29$-31 mag arcsec$^{-2}$ in the optical broadbands and $f_{1\sigma}\sim1$-$2\times 10^{-18}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$ in the H$\alpha$ narrowband. To conclude the paper, several possible scientific applications for this data set are described.

Multi-conjugate adaptive optics (MCAO) will assist a new era of ground-based astronomical observations with the extremely large telescopes and the Very Large Telescope. High precision relative astrometry is among the main science drivers of these systems and challenging requirements have been set for the astrometric measurements. A clear understanding of the astrometric error budget is needed and the impact of the MCAO correction has to be taken into account. In this context, we propose an analytical formulation to estimate the residual phase produced by an MCAO correction in any direction of the scientific field of view. The residual phase, computed in the temporal frequency domain, allows to consider the temporal filtering of the turbulent phase from the MCAO loop and to extract the temporal spectrum of the residuals, as well as to include other temporal effects such as the scientific integration time. The formulation is kept general and allows to consider specific frameworks by setting the telescope diameter, the turbulence profile, the guide stars constellation, the deformable mirrors configuration, the modes sensed and corrected and the tomographic reconstruction algorithm. The formalism is presented for both a closed loop and a pseudo-open loop control. We use our results to investigate the effect of tip-tilt residuals on MCAO-assisted astrometric observations. We derive an expression for the differential tilt jitter power spectrum that also includes the dependence on the scientific exposure time. Finally, we investigate the contribution of the differential tilt jitter error on the future astrometric observations with MAVIS and MAORY.

The millisecond pulsar MAXI J1816--195 was recently discovered by MAXI in 2022 May. We have studied different properties of the pulsar using data from NuSTAR and NICER observations. The position of the source is measured by NuSTAR as RA = $18^h 16^m 52^s.40$, Dec = $-19^o37^{'} 58^{''}.35$. The unstable burning of accreted material on the surface of neutron stars induces thermonuclear (Type-I) bursts. Several thermonuclear bursts have been detected from the source during the outburst. We study the evolution of burst profile with flux and energy using NuSTAR and NICER observations. During the NuSTAR observation, a total of four bursts were detected from the source. The duration of each burst was around $\sim$ 30 s and the ratio of peak to persistent count rate is $\sim$ 26 as seen from the NuSTAR data. The thermonuclear bursts are modeled to determine the burst timing parameters using a sharp linear rise and exponential decay function. The burst profiles show a relatively long tail in lower energies. The hardness ratio during the thermonuclear bursts shows significant variation as observed by NuSTAR. We successfully model the broadband burst-resolved spectra with a combination of an absorbed blackbody along with a non-thermal component to account for the persistent emission. The burst-resolved spectral parameters show significant evolution during the burst. During the peak of the burst, the Eddington luminosity is found to be $\sim 3.7 \times 10^{38}$ erg s$^{-1}$. The burst-resolved spectral parameters provide a source distance of $8.5\pm1.2$ kpc for isotropic burst emission.

A new era of ground-based observations, either in the infrared with the next-generation of 25-40m extremely large telescopes or in the visible with the 8m Very Large Telescope, is going to be assisted by multi-conjugate adaptive optics (MCAO) to restore the unprecedented resolutions potentially available for these systems in absence of atmospheric turbulence. Astrometry is one of the main science drivers, as MCAO can provide good quality and uniform correction over wide field of views ($\sim$ 1 arcmin) and offer a large number of reference sources with high image quality. The requirements have been set to very high precisions on the differential astrometry (e.g. 50$\mu$as for MICADO/MORFEO - formerly known as MAORY - at the Extremely Large Telescope) and an accurate analysis of the astrometric error budget is needed. In this context, we present an analysis of the impact of MCAO atmospheric tip-tilt residuals on relative astrometry. We focus on the effects of the scientific integration time on tip-tilt residuals, that we model through the temporal transfer function of the exposure. We define intra- and inter-exposure tip-tilt residuals that we use in the estimation of the centroiding error and the differential tilt jitter error within the astrometric error budget. As a case study, we apply our results in the context of the MORFEO astrometric error budget.

In the context of adaptive optics for astronomy, one can rely on the statistics of the turbulent phase to assess a part of the system's performance. Temporal statistics with one source and spatial statistics with two sources are well-known and are widely used for classical adaptive optics systems. A more general framework, including both spatial and temporal statistics, can be useful for the analysis of the existing systems and to support the design of the future ones. In this paper, we propose an expression of the temporal cross power spectral densities of the turbulent phases in two distinct beams, that is from two different sources to two different apertures. We either consider the phase as it is, without piston, or as its decomposition on Zernike modes. The general formulas allow to cover a wide variety of configurations, from single-aperture to interferometric telescopes equipped with adaptive optics, with the possibility to consider apertures of different sizes and/or sources at a finite distance. The presented approach should lead to similar results with respect to existing methods in the Fourier domain, but it is focused on temporal frequencies rather than spatial ones, which might be convenient for some aspects such as control optimization. To illustrate this framework with a simple application, we demonstrate that the wavefront residual due to the anisoplanatism error in a single-conjugated adaptive optics system is overestimated when it is computed from covariances without taking into account the temporal filtering of the adaptive optics loop. We also show this overestimation in the case of a small-baseline interferometer, for which the two beams are significantly correlated.

Various theories have been proposed to explain the Moon's current inclined orbit. We test the viability of these theories by reconstructing the thermal-orbital history of the Moon. We build on past thermal-orbital models and incorporate the evolution of the lunar figure including a fossil figure component. Obliquity tidal heating in the lunar magma ocean would have produced rapid inclination damping, making it difficult for an early inclination to survive to the present-day. An early inclination is preserved only if the solid-body of the early Moon were less dissipative than at present. If instabilities at the Laplace plane transition were the source of the inclination, then the Moon had to recede slowly, which is consistent with previous findings of a weakly dissipative early Earth. If collisionless encounters with planetesimals up to 140 Myr after Moon formation excited the inclination, then the Moon had to migrate quickly to pass through the Cassini state transition at 33 Earth radii and reach a period of limited inclination damping. The fossil figure was likely established before 16 Earth radii to match the present-day degree-2 gravity field observations.

For the past decade, it has been suggested that intermediate polars (IPs), a subclass of magnetic cataclysmic variables (CVs), are one of the main contributors to the hard diffuse X-ray emission from the Galactic center (GC) and Galactic ridge. In our ongoing \emph{XMM-Newton} survey of the central region of the Galactic disk ($20^\circ\times2^\circ$), we detected a persistent IP candidate, $1.7^\circ$ away from the GC. In this work, we better characterize the behavior of this source by looking at the new and archival XMM-Newton data. We performed a detailed X-ray spectral modeling of the source. Furthermore, we searched for X-ray pulsations in the light curve as well as its counterpart at other wavelengths. The XMM-Newton spectrum (0.8--10 keV) of the source is described by a partial covering collisionally ionized diffuse gas with plasma temperature $kT=15.7^{+20.9}_{-3.6}$ keV. In addition, the spectrum shows the presence of iron lines at $E=6.44$, 6.65, and 6.92 keV with equivalent widths of $194^{+89}_{-70}$, $115^{+79}_{-75}$, and $98^{+93}_{-74}$ eV, respectively. The X-ray light curve shows a coherent modulation with a period of $P=432.44\pm0.36$ s, which we infer is the spin period of the white dwarf. The white dwarf mass estimated from fitting a physical model to the spectrum results in $M_{\rm WD}=1.05^{+0.16}_{-0.21}\ M_{\odot}$. We were able to find a likely optical counterpart in the Gaia catalog with a G magnitude of 19.26, and the distance to the source derived from the measured Gaia parallax is $\sim$4.3 kpc. We provide an improved source localization with subarcsec accuracy. The spectral modeling of the source indicates the presence of intervening circumstellar gas, which absorbs the soft X-ray photons. The measured equivalent width of the iron lines and the detection of the spin period in the light curve are consistent with those from IPs.

We present the physical and dynamical properties of the recently discovered active asteroid (248370) 2005 QN173 (aka 433P). From our observations, we derived two possible rotation period solutions of 2.7 and 4.1 hours. The corresponding light curve amplitudes computed after correcting for the effect of coma are 0.28 and 0.58 mag, respectively. Both period solutions are shorter than the critical rotation limit computed for a strengthless triaxial ellipsoid, suggesting that rotation mass shedding should at least partly be responsible for the observed activity. We confirm that the activity level is fading further, but at a very modest rate of only 0.006 mag/day, still also compatible with sublimation-driven activity. We found that 248370 likely belongs to the Themis asteroid family, making it a fourth main-belt comet associated with this group. Orbital characteristics of 248370 are also consistent with its origin in the young 288P cluster of asteroids. The 288P cluster is associated with its namesake main-belt comet, providing an exciting possibility for a comparative analysis of intriguing main-belt comets 248370 and 288P.

Upcoming large-scale spectroscopic surveys such as WEAVE and 4MOST will provide thousands of spectra of massive stars, which need to be analysed in an efficient and homogeneous way. Studies on massive stars are usually based on samples of a few hundred objects which pushes current spectroscopic analysis tools to their limits because visual inspection is necessary to verify the spectroscopic fit. The novel spectroscopic analysis pipeline takes advantage of the statistics that large samples provide, and determines the model error to account for imperfections in stellar atmosphere codes due to simplified, wrong or missing physics. Considering observational plus model uncertainties improve spectroscopic fits. The pipeline utilises the entire spectrum rather than selected diagnostic lines allowing a wider range of temperature from B to early O stars to be analysed. A small fraction of stars like peculiar, contaminated or spectroscopic binaries require visual inspection, which are identified through their larger uncertainties.

Mass loss through stellar winds plays a dominant role in the evolution of massive stars. Very massive stars (VMSs, $> 100 M_{\odot}$) display Wolf-Rayet spectral morphologies (WNh) whilst on the main-sequence. Bestenlehner (2020) extended the elegant and widely used stellar wind theory by Castor, Abbott & Klein (1975) from the optically thin (O star) to the optically thick main-sequence (WNh) wind regime. The new mass-loss description is able to explain the empirical mass-loss dependence on the Eddington parameter and is suitable for incorporation into stellar evolution models for massive and very massive stars. The prescription can be calibrated with the transition mass-loss rate defined in Vink & Gr\"afener (2012). Based on the stellar sample presented in Bestenlehner et al. (2014) we derive a mass-loss recipe for the Large Magellanic Cloud using the new theoretical mass-loss prescription of Bestenlehner (2020).

The importance of photometric galaxy redshift estimation is rapidly increasing with the development of specialised powerful observational facilities. We develop a new photometric redshift estimation workflow TOPz to provide reliable and efficient redshift estimations for the upcoming large-scale survey J-PAS which will observe 8500 deg2 of the northern sky through 54 narrow-band filters. TOPz relies on template-based photo-z estimation with some added J-PAS specific features and possibilities. We present TOPz performance on data from the miniJPAS survey, a precursor to the J-PAS survey with an identical filter system. First, we generated spectral templates based on the miniJPAS sources using the synthetic galaxy spectrum generation software CIGALE. Then we applied corrections to the input photometry by minimising systematic offsets from the template flux in each filter. To assess the accuracy of the redshift estimation, we used spectroscopic redshifts from the DEEP2, DEEP3, and SDSS surveys, available for 1989 miniJPAS galaxies with r < 22 magAB. We also tested how the choice and number of input templates, photo-z priors, and photometric corrections affect the TOPz redshift accuracy. The general performance of the combination of miniJPAS data and the TOPz workflow fulfills the expectations for J-PAS redshift accuracy. Similarly to previous estimates, we find that 38.6% of galaxies with r < 22 mag reach the J-PAS redshift accuracy goal of dz/(1 + z) < 0.003. Limiting the number of spectra in the template set improves the redshift accuracy up to 5%, especially for fainter, noise-dominated sources. Further improvements will be possible once the actual J-PAS data become available.

The magnetic fields observed in Ap-stars, white dwarfs, and neutron stars are known to be stable for long times. However, the physical conditions inside the stellar interiors that allow these states are still a matter of research. It has been formally demonstrated that both purely toroidal and purely poloidal magnetic fields develop instabilities at some point in the star. On the other hand, numerical simulations have proved the stability of roughly axisymmetric magnetic field configurations inside stably stratified stars. These configurations consist of mutually stabilizing toroidal and poloidal components in a twisted torus shape. Previous studies have proposed rough upper and lower bounds on the ratio of the magnetic energy in the toroidal and poloidal components of the magnetic field. With the purpose of mapping out the parameter space under which such configurations remain stable, we used the Pencil Code to perform 3D magnetohydrodynamic simulations of the evolution of the magnetic field in non-rotating, non-degenerate stars in which viscosity is the only dissipation mechanism, both for stars with a uniform (barotropic) and radially increasing (stably stratified) specific entropy. Furthermore, we considered different conditions regarding the degree of stable stratification and the magnetic energy in each component, roughly confirming the previously suggested stability boundaries for the magnetic field.

We propose a new scenario of early dark energy (EDE) with a dark Higgs trapped at the origin. To keep this dark Higgs trapped until around the matter-radiation equality, we use dark photons produced non-thermally by coherent oscillations of axions, which have a much stronger trapping effect than thermal mass. When the trapping ends, the dark Higgs quickly decays into dark photons, which are then red-shifted as radiation. The dark Higgs EDE scenario works well for an ordinary Mexican-hat potential, and the dark Higgs naturally sits at the origin from the beginning, since it is the symmetry-enhanced point. Thus, unlike the axion EDE, there is no need for elaborate potentials or fine-tuning with respect to the initial condition. Interestingly, the axion not only produces dark photons, but also explains dark matter. We find the viable parameter region of the axion decay constant and the axion mass where dark matter and the $H_0$ tension can be simultaneously explained. We also discuss the detectability of the axion in the presence of axion-photon coupling, and show that the axion can be the QCD axion.

In this work, we explore the possibility of using artificial neural networks to impose constraints on teleparallel gravity and its $f(T)$ extensions. We use the available Hubble parameter observations from cosmic chronometers and baryon acoustic oscillations from different galaxy surveys. We discuss the procedure for training a network model to reconstruct the Hubble diagram. Further, we describe the procedure to obtain $H'(z)$, the first order derivative of $H(z)$, in a novel way. These analyses are complemented with two presently debated values of $H_0$, namely, the local measurements by the SH0ES team ($H_0^{\text{R20}} = 73.2 \pm 1.3$~km~Mpc$^{-1}$~s$^{-1}$) and the updated TRGB calibration from the Carnegie Supernova Project ($H_0^{\text{TRGB}} = 69.8 \pm 1.9$~km~Mpc$^{-1}$~s$^{-1}$), respectively. Additionally, we investigate the validity of the concordance model, through some cosmological null tests with these reconstructed data sets. Finally, we reconstruct the allowed $f(T)$ functions for different combinations of the observational Hubble data sets. Results show that the $\Lambda$CDM model lies comfortably included at the 1$\sigma$ confidence level for all the examined cases.

We introduce our Radial Velocity Survey for Planets around Young stars (RVSPY), characterise our target stars, and search for substellar companions at orbital separations smaller than a few au from the host star. We use the FEROS spectrograph to obtain high signal-to-noise spectra and time series of precise radial velocities (RVs) of 111 stars most of which are surrounded by debris discs. Our target stars have spectral types between early F and late K, a median age of 400 Myr, and a median distance of 45 pc. We determine for all target stars their basic stellar parameters and present the results of the high-cadence RV survey and activity characterization. We achieve a median single-measurement RV precision of 6 m/s and derive the short-term intrinsic RV scatter of our targets (median 22 m/s), which is mostly caused by stellar activity and decays with age from >100 m/s at <20 Myr to <20 m/s at >500 Myr. We discover six previously unknown close companions with orbital periods between 10 and 100 days, three of which are low-mass stars, and three are in the brown dwarf mass regime. We detect no hot companion with an orbital period <10 days down to a median mass limit of ~1 M_Jup for stars younger than 500 Myr, which is still compatible with the established occurrence rate of such companions around main-sequence stars. We find significant RV periodicities between 1.3 and 4.5 days for 14 stars, which are, however, all caused by rotational modulation due to starspots. We also analyse the TESS photometric time series data and find significant periodicities for most of the stars. For 11 stars, the photometric periods are also clearly detected in the RV data. We also derive stellar rotation periods ranging from 1 to 10 days for 91 stars, mostly from TESS data. From the intrinsic activity-related short-term RV jitter, we derive the expected mass-detection thresholds for longer-period companions.

Identifying long-period eclipsing binaries with space-based photometry is still a challenge even in the century of space telescopes due to the relatively short observation sequences and short lifetime of these missions. The Transiting Exoplanet Survey Satellite (TESS) space telescope is an appropriate tool to supplement previous space-based observations. In this paper we report the first results of the eclipse timing variation (ETV) analyses of eclipsing binaries (EBs) measured by CoRoT and TESS space telescopes. Among the 1428 EB candidates we found 4 new potential triple candidates, for which ETV was analysed and fitted by the well-known light-travel-time effect (LTTE). One of them shows significant phase shift in its folded light curve which required extra care. In this paper we also present some other systems showing significant ETV signals that could be explained by mass transfer or apsidal motion.

WEAVE is the new wide-field spectroscopic facility for the prime focus of the William Herschel Telescope in La Palma, Spain. Its fibre positioner is essential for the accurate placement of the spectrograph's ~960-fibre multiplex. To maximise the assignment of its optical fibres, WEAVE uses a simulated annealing algorithm called Configure, which allocates the fibres to targets in the field of view. We have conducted an analysis of the algorithm's behaviour using a subset of mid-tier WEAVE-LOFAR fields, and adjusted the priority assignment algorithm to optimise the total fibres assigned per field, and the assignment of fibres to the higher priority science targets. The output distributions have been examined, to investigate the implications for the WEAVE science teams.

WEAVE is the new wide-field spectroscopy facility for the prime focus of the William Herschel Telescope in La Palma, Spain. Its fibre positioner is essential for the accurate placement of the spectrograph's 960 fibre multiplex. We provide an overview of the recent maintenance, flexure modifications, and calibration measurements conducted at the observatory prior to the final top-end assembly. This work ensures that we have a complete understanding of the positioner's behaviour as it changes orientation during observations. All fibre systems have been inspected and repaired, and the tumbler structure contains new clamps to stiffen both the internal beam and the retractor support disk onto which the field plates attach. We present the updated metrology procedures and results that will be verified on-sky.

The nearby, luminous infrared galaxy (LIRG) NGC 7469 hosts a Seyfert nucleus with a circumnuclear star-forming ring and is thus the ideal local laboratory for investigating the starburst--AGN connection in detail. We present integral-field observations of the central 1.3 kpc region in NGC 7469 obtained with the JWST Mid-InfraRed Instrument. Molecular and ionized gas distributions and kinematics at a resolution of {\sim}100 pc over the 4.9 - 7.6{\mu}m region are examined to study gas dynamics influenced by the central AGN. The low-ionization [Fe II] {\lambda}5.34{\mu}m and [Ar II] {\lambda}6.99{\mu}m lines are bright on the nucleus and in the starburst ring, as opposed to H2 S(5) {\lambda}6.91{\mu}m which is strongly peaked at the center and surrounding ISM. The high-ionization [Mg V] line is resolved and shows a broad, blueshifted component associated with the outflow. It has a nearly face-on geometry that is strongly peaked on the nucleus, where it reaches a maximum velocity of -650 km/s, and extends about 400 pc to the East. Regions of enhanced velocity dispersion in H2 and [Fe II] {\sim}180 pc from the AGN that also show high L(H2)/L(PAH) and L([Fe II])/L(Pf{\alpha}) ratios to the W and N of the nucleus pinpoint regions where the ionized outflow is depositing energy, via shocks, into the dense interstellar medium between the nucleus and the starburst ring. These resolved mid-infrared observations of the nuclear gas dynamics demonstrate the power of JWST and its high-sensitivity integral-field spectroscopic capability to resolve feedback processes around supermassive black holes in the dusty cores of nearby LIRGs.

Understanding the long-term evolution of hierarchical triple systems is challenging due to its inherent chaotic nature, and it requires computationally expensive simulations. Here we propose a convolutional neural network model to predict the stability of hierarchical triples by looking at their evolution during the first $5 \times 10^5$ inner binary orbits. We employ the regularized few-body code TSUNAMI to simulate $5\times 10^6$ hierarchical triples, from which we generate a large training and test dataset. We develop twelve different network configurations that use different combinations of the triples' orbital elements and compare their performances. Our best model uses 6 time-series, namely, the semimajor axes ratio, the inner and outer eccentricities, the mutual inclination and the arguments of pericenter. This model achieves an area under the curve of over $95\%$ and informs of the relevant parameters to study triple systems stability. All trained models are made publicly available, allowing to predict the stability of hierarchical triple systems $200$ times faster than pure $N$-body methods.

From October 2004 to May 2022, the concentration of radon in the air was measured at a depth of 700 m in the Yangyang underground laboratory. The average rates in the two experimental areas, called A6 and A5, were measured as 53.4\pm0.2 Bq/m^3 and 33.5\pm0.1 Bq/m^3, respectively. The lower rate in the A5 area was caused by the improved temperature control and ventilation. In particular, these radon rates are correlated to the local temperature of the area, with a correlation coefficient r = 0.22. Therefore, the radon rates displayed a seasonal variation, because the local temperature driven by the overground season influences air ventilation in the experimental areas. A cosine fit on the annual residual rates exhibited the maximum amplitude on August 31 \pm 6 d every year.

The equation of state (EoS) for cold dense matter inside neutron stars is investigated by using holographic QCD models in the framework of the Einstein-Maxwell-dilaton (EMD) system and the improved Karch-Katz-Son-Stephanov (KKSS) action for matter part. This method of describing holographic nuclear matter in the EMD$+$KKSS framework is different from that by using the Dirac-Born-Infeld (DBI) action and the Chern-Simons (CS) terms. Combining with the Hebeler-Lattimer-Pethick-Schwenk (HLPS) intermediate equation of state (EoS), the hybrid EoS inside the neutron stars is constructed. The obtained hybrid EoS is located in the range that is defined by the low-density chiral effective theory, the high-density perturbative QCD, and the polytropic interpolations between them, and is constrained by the astrophysics observations. The square of the sound velocity reaches a maximum value larger than $0.8$ in the region of $2-5$ times the saturation baryon number density and approaches the conformal limit at the high baryon density range. The mass-radius relation and the tidal deformability of the neutron stars are in agreement with astrophysical measurements. The possible maximum mass for the neutron star is about $2.5 M_{\odot}$ and the radius is about $12 \mathrm{km}$ then. It is noticed that the holographic quark matter branch in the mass-radius relation is always unstable and the holographic nuclear matter can produce a stable branch. These results indicate that even in the core of the NS, the matter is still in the confinement phase and the quark matter is not favored.

Gravitational waves travel through the distributions of matter and dark energy during propagation. For this reason, gravitational waves emitted from binary compact objects serve as a useful tool especially to probe the nature of dark energy. The geometrical optics approximation is a conventional way of investigating wave propagation. However, the approximation becomes less accurate as the wavelength approaches the curvature radius of the background, which can occur in generic situations. In this paper, we suggest a formulation for higher-order corrections of the geometrical optics expansion, applied to Horndeski theory which accommodates many dark energy models. At the level of the background, assuming that the derivative of the scalar field is non-vanishing and timelike, we choose the time slices to coincide with the contours of the scalar field. This choice of the background time slices is advantageous as the sound cones of both scalar and tensor gravitational waves are upright with respect to the background time slices whenever the scalar field behaves as a perfect fluid. We then analyze the equations of motion for scalar and tensor components of gravitational waves at the leading and next-to-leading order in the geometrical optics expansion, deriving the evolution equations for their amplitudes under certain conditions. In particular, for Generalized Brans-Dicke theories, we find a simple description of equations for gravitational waves in terms of an effective metric.

The crossed-sine wavefront sensor (WFS) is a pupil plane wavefront sensor that measures the first derivatives of the wavefront. It is made by three main components: a gradient transmission filter (GTF) built from a product of sine functions rotated by 45 degrees around the optical axis, a 2x2 mini-lens array (MLA) at the focus of the tested optical system and a detector array located on a plane conjugated to the pupil. The basic principle consists in acquiring four pupil images simultaneously, each image being observed from different points located behind the GTF. After the simulation work which demonstrated the wavefront reconstruction capability, we are now in the phase of implementation of the prototype in the lab. The crossed-sine WFS could achieve a simultaneous high spatial resolution at the pupil of the tested optics and absolute measurement accuracy comparable to that attained by laser-interferometers. In this paper we introduce seven customized phase masks and make measurements of them.First tests and resultsare demonstrated, based on which we explore the performance of our crossed-sine WFS and make comparisons with that of the laser-interferomete

The long-standing model-independent annual modulation effect measured by DAMA deep underground at Gran Sasso Laboratory with different experimental configurations is summarized and perspectives will be highlighted. DAMA/LIBRA-phase2 set-up, $\simeq$ 250 kg highly radio-pure NaI(Tl) confirms the evidence of a signal that meets all the requirements of the model independent Dark Matter annual modulation signature at high C.L.; the full exposure is 2.86 ton $\times$ yr over 22 annual cycles. The experiment is currently collecting data in the DAMA/LIBRA-phase2 empowered configuration with an even lower software energy threshold. Other recent claims are shortly commented.

We construct a holographic model for dark energy in the Brans-Dicke cosmology by using the holographic principle considering the Barrow entropy instead of the standard Bekenstein-Hawking one. The former arises from the effort to account for quantum-gravitational effects in black-hole physics and, according to the gravity-thermodynamic conjecture, in the cosmological framework. In order to explore the cosmological consequences of our model, we consider the Hubble horizon as the IR cutoff. We investigate both the non-interacting and interacting cases with the sign-changeable and linear interactions, showing that they can explain the present accelerated phase of the Universe expansion, in contrast to the standard Holographic Dark Energy model. We then perform the classical stability analysis using the squared sound speed. We find that, whilst the non-interacting model is unstable against the small perturbations, the sign-changeable interacting one can be stable only for suitable values of the model parameters. On the other hand, the linear interacting model always predicts a stable Universe. The consistency of our model with cosmological observations is discussed.

We considered a vacuum polarization inside a galaxy in the eikonal approximation and found that two possible types of polarization exist. The first type is described by the equation of state $p=\rho/3$, similar to radiation. Using the conformally-unimodular metric allows constructing a nonsingular solution for this vacuum ``substance'', if a compact astrophysical object exists in the galaxy's center. As a result, a ``dark'' galactical halo appears that increases the rotation velocity of a test particle as a function of the distance from a galactic center. The second type of vacuum polarization has a more complicated equation of state. As a static physical effect, it produces renormalization of the gravitational constant, thus, causing no static halo. However, a nonstationary polarization of the second type, resulting from an exponential increase (or decrease) of the galactic nuclei mass with time in some hypothetical time-dependent process, produces a gravitational potential looking like a dark matter halo.

Recently, a new nonlinear mechanism for black hole scalarization, different from the standard spontaneous scalarization, was demonstrated to exist for scalar Gauss-Bonnet theories in which no tachyonic instabilities can occur. Thus Schwarzschild black hole is linearly stable but instead nonlinear instability can kick-in. In the present paper we extend on this idea in the case of multi-scalar Gauss-Bonnet gravity with exponential coupling functions of third and fourth leading order in the scalar field. The main motivation comes from the fact that these theories admit hairy compact objects with zero scalar charge, thus zero scalar-dipole radiation, that automatically evades the binary pulsar constraints on the theory parameters. We demonstrate numerically the existence of scalarized black holes for both coupling functions and for all possible maximally symmetric scalar field target spaces. The thermodynamics and the stability of the obtained solution branches is also discussed.

We study the gravitational radiation emission efficiency $\Gamma$ of superconducting cosmic strings. We demonstrate, by using a solvable model of transonic strings, that the presence of a current leads to a suppression of the gravitational emission of cusps, kinks and different types of loops. We also show that, when a current is present, the spectrum of emission of loops with cusps is exponentially suppressed as the harmonic mode increases, thus being significantly different from the power law spectrum of currentless loops. Furthermore, we establish a phenomenological relationship between $\Gamma$ and the value of the current on cosmic strings. We conjecture that this relation should be valid for an arbitrary type of current-carrying string. We use this result to study the potential impact of current on the stochastic gravitational wave background generated by cosmic strings with additional degrees of freedom and show that both the amplitude and shape of the spectrum may be significantly affected.