Papers for Tuesday, Jan 04 2022

Molecular maser lines are signposts of high-mass star formation, probing excitation and kinematics of very compact regions in the close environment of young stellar objects and providing useful targets for trigonometric parallax measurements. Only a few NH$_{3}$ (9,6) masers were known so far, and their origin is still poorly understood. Here we aim to find new NH$_{3}$ (9,6) masers to provide a better observational basis to study their role in high-mass star-forming regions. We carried out NH$_{3}$ (9,6) observations toward Cepheus A and G34.26$+$0.15 with the Effelsberg-100 m telescope and the Karl G. Janksy Very Large Array. We discovered new NH$_{3}$ (9,6) masers in Cep A and G34.26$+$0.25, which increases the number of high-mass star-forming regions hosting NH$_{3}$ (9,6) masers from five to seven. Long term monitoring (20 months) at Effelsberg shows that the intensity of the (9,6) maser in G34.26$+$0.25 is decreasing, while the Cep A maser remains stable. Compared to the Effelsberg data and assuming linear variations between the epochs of observation, the JVLA data indicate no missing flux. This suggests that the NH$_3$ (9,6) emission arises from single compact emission regions that are not resolved by the interferometric measurements. As JVLA imaging shows, the NH$_{3}$ (9,6) emission in Cep A originates from a sub-arcsecond sized region, slightly to the west of the peak position of the 1.36\,cm continuum object, HW2. In G34.26$+$0.25, three NH$_{3}$ (9,6) maser spots are observed: one is close to the head of the cometary ultracompact \h2 region C and the other two are emitted from a compact region to the west of the hypercompact \h2 region A. The newly found (9,6) masers appear to be related to outflows. Higher angular resolution of JVLA and VLBI observations are needed to provide more accurate positions and constraints for pumping scenarios.

Most stellar evolution models predict that black holes (BHs) should not exist above approximately $50-70$ M$_\odot$. However, recent LIGO/Virgo detections indicate the existence of BHs with masses at and above this threshold. We suggest that massive BHs, including intermediate mass black holes (IMBHs), can form in galactic nuclei through collisions between stellar-mass black holes and the surrounding main-sequence stars. Considering dynamical processes such as collisions, mass segregation, and relaxation, we find that this channel can be quite efficient, forming IMBHs as massive as $10^4$ M$_\odot$. Our results suggest that massive black holes and IMBHs may be ubiquitous in galactic centres. This formation channel also has implications for observations. Collisions between stars and BHs can produce electromagnetic signatures, for example, from x-ray binaries and tidal disruption events. Additionally, formed through this channel, both black holes in the mass gap and IMBHs can merge with the supermassive black hole at the center of a galactic nucleus through gravitational waves. These gravitational wave events are extreme and intermediate mass ratio inspirals (EMRIs and IMRIs, respectively).

As is common with the collection of astronomical data, signals are frequently dominated by noise. However, when performing FTs of light curves, re-binning data can improve the signal-to-noise ratio (SNR) at lower frequencies. Using data collected from the Kepler space telescope, we sequentially re-binned data three times to investigate the SNR improvement of lower frequency (< 17 microHz) variability in white dwarf KIC 8626021. We found that the SNR at approximately 5.8 microHz greatly improved through this process, and we postulate that this frequency is linked to the rotation of KIC 8626021.

We present a model of the Cassini state of Mercury that comprises an inner core, a fluid core and a mantle. Our model includes inertial and gravitational torques between interior regions, and viscous and electromagnetic (EM) coupling at the boundaries of the fluid core. We show that the coupling between Mercury's interior regions is sufficiently strong that the obliquity of the mantle spin axis deviates from that of a rigid planet by no more than 0.01 arcmin. The mantle obliquity decreases with increasing inner core size, but the change between a large and no inner core is limited to 0.015 arcmin. EM coupling is stronger than viscous coupling at the inner core boundary and, if the core magnetic field strength is above 0.3 mT, locks the fluid and solid cores into a common precession motion. Because of the strong gravitational coupling between the mantle and inner core, the larger the inner core is, the more this co-precessing core is brought into an alignment with the mantle, and the more the obliquity of the polar moment of inertia approaches that expected for a rigid planet. The misalignment between the polar moment of inertia and mantle spin axis increases with inner core size, but is limited to 0.007 arcmin. Our results imply that the measured obliquities of the mantle spin axis and polar moment of inertia should coincide at the present-day level of measurement errors, and cannot be distinguished from the obliquity of a rigid planet.

We report on a search for persistent radio emission from the one-off Fast Radio Burst (FRB) 20190714A, as well as from two repeating FRBs, 20190711A and 20171019A, using the MeerKAT radio telescope. For FRB 20171019A we also conducted simultaneous observations with the High Energy Stereoscopic System (H.E.S.S.) in very high energy gamma rays and searched for signals in the ultraviolet, optical, and X-ray bands. For this FRB, we obtain a UV flux upper limit of 1.39x10^-16 erg/cm^-2/s/Amstrong, X-ray limit of ~ 6.6x10^-14 erg/cm^-2/s and a limit on the very-high-energy gamma-ray flux (Phi) (E > 120 GeV) < 1.7 x 10^-12 erg/cm^-2/s. We obtain a radio upper limit of ~15 microJy/beam for persistent emission at the locations of both FRBs 20190711A and 20171019A, but detect diffuse radio emission with a peak brightness of ~53 microJy/beam associated with FRB 20190714A at z = 0.2365. This represents the first detection of the radio continuum emission potentially associated with the host (galaxy) of FRB 20190714A, and is only the third known FRB to have such an association. Given the possible association of a faint persistent source, FRB 20190714A may potentially be a repeating FRB whose age lies between that of FRB 20121102A and FRB 20180916A. A parallel search for repeat bursts from these FRBs revealed no new detections down to a fluence of 0.08 Jy ms for a 1 ms duration burst.

Dwarf spheroidal (dSph) galaxies are believed to be strongly dark matter dominated and thus are considered perfect objects to study dark matter distribution and test theories of structure formation. They possess resolved, multiple stellar populations that offer new possibilities for modeling. A promising tool for the dynamical modeling of these objects is the Schwarzschild orbit superposition method. In this work we extend our previous implementation of the scheme to include more than one population of stars and a more general form of the mass-to-light ratio function. We tested the improved approach on a nearly spherical, gas-free galaxy formed in the cosmological context from the Illustris simulation. We modeled the binned velocity moments for stars split into two populations by metallicity and demonstrate that in spite of larger sampling errors the increased number of constraints leads to significantly tighter confidence regions on the recovered density and velocity anisotropy profiles. We then applied the method to the Fornax dSph galaxy with stars similarly divided into two populations. In comparison with our earlier work, we find the anisotropy parameter to be slightly increasing, rather than decreasing, with radius and more strongly constrained. We are also able to infer anisotropy for each stellar population separately and find them to be significantly different.

The technique of normal-mode coupling is a powerful tool with which to seismically image non-axisymmetric phenomena in the Sun. Here we apply mode coupling in the Cartesian approximation to probe steady, near-surface flows in the Sun. Using Doppler cubes obtained from the Helioseismic and Magnetic Imager onboard the Solar Dynamics Observatory, we perform inversions on mode-coupling measurements to show that the resulting divergence and radial vorticity maps at supergranular length scales ($\sim$30 Mm) near the surface compare extremely well with those obtained using the Local Correlation Tracking method. We find that the Pearson correlation coefficient is $\geq$ 0.9 for divergence flows, while $\geq$ 0.8 is obtained for the radial vorticity.

Recent hydrodynamical simulations of convection in a solar-like model suggest that penetrative convective flows at the boundary of the convective envelope modify the thermal background in the overshooting layer. Based on these results, we implement in one-dimensional stellar evolution codes a simple prescription to modify the temperature gradient below the convective boundary of a solar model. This simple prescription qualitatively reproduces the behaviour found in the hydrodynamical simulations, namely a local heating and smoothing of the temperature gradient below the convective boundary. We show that introducing local heating in the overshooting layer can reduce the sound-speed discrepancy usually reported between solar models and the structure of the Sun inferred from helioseismology. It also affects key quantities in the convective envelope, such as the density, the entropy, and the speed of sound. These effects could help reduce the discrepancies between solar models and observed constraints based on seismic inversions of the Ledoux discriminant. Since mixing due to overshooting and local heating are the result of the same convective penetration process, the goal of this work is to invite solar modellers to consider both processes for a more consistent approach.

Empirical evidences suggest the existence of a period-age relation for long-period variables (LPVs). Yet, this property has so far been little studied on theoretical grounds. We aim to examine the period-age relation using the results from recent nonlinear pulsation calculations. We combine isochrone models with theoretical periods to simulate the distribution of fundamental mode LPV pulsators, that include Miras, in the period-age plane, and compare it with observations of LPVs in Galactic and Magellanic Clouds clusters. In agreement with observations, models predict that the fundamental mode period decreases with increasing age because of the dominant role of mass in shaping stellar structure and evolution. At a given age, the period distribution shows a non-negligible width and is skewed towards short periods, except for young C-rich stars. As a result, the period-age relations of O-rich and C-rich models are predicted to have different slopes. We derive best-fit relations describing age and initial mass as a function of the fundamental mode period for both O- and C-rich models. The study confirms the power of the period-age relations to study populations of LPVs of specific types, either O-rich or C-rich, on statistical grounds. In doing so, it is recommended not to limit a study to Miras, which would make it prone to selection biases, but rather to include semi-regular variables that pulsate predominantly in the fundamental mode. The use of the relations to study individual LPVs, on the other hand, requires more care given the scatter in the period distribution predicted at any given age.

Here, we study the temperature structure of flaring and non-flaring coronal loops, using extracted loops from images taken in six extreme ultraviolet (EUV) channels recorded by Atmospheric Imaging Assembly (AIA)/ Solar Dynamic Observatory (SDO). We use data for loops of X2.1-class-flaring active region (AR11283) during 22:10UT till 23:00UT, on 2011, September 6; and non-flaring active region (AR12194) during 08:00:00UT till 09:00:00UT on 2014, October 26. By using spatially-synthesized Gaussian DEM forward-fitting method, we calculate the peak temperatures for each strip of the loops. We apply the Lomb-Scargle method to compute the oscillations periods for the temperature series of each strip. The periods of the temperature oscillations for the flaring loops are ranged from 7 min to 28.4 min. These temperature oscillations show very close behavior to the slow-mode oscillation. We observe that the temperature oscillations in the flaring loops are started at least around 10 minutes before the transverse oscillations and continue for a long time duration even after the transverse oscillations are ended. The temperature amplitudes are increased at the flaring time (during 20 min) in the flaring loops. The periods of the temperatures obtained for the non-flaring loops are ranged from 8.5 min to 30 min, but their significances are less (below 0.5) in comparison with the flaring ones (near to one). Hence the detected temperature periods for the non-flaring loops' strips are less probable in comparison with the flaring ones, and maybe they are just fluctuations. Based on our confined observations, it seems that the flaring loops' periods show more diversity and their temperatures have wider ranges of variation than the non-flaring ones. More accurate commentary in this respect requires more extensive statistical research and broader observations.

We present a spectroscopic investigation with VLT/X-shooter of seven candidate extremely strong damped Lyman-$\alpha$ absorption systems (ESDLAs, $N(\text{HI})\ge 5\times 10^{21}$ cm$^{-2}$) observed along quasar sightlines. We confirm the extremely high column densities, albeit slightly (0.1~dex) lower than the original ESDLA definition for four systems. We measured low-ionisation metal abundances and dust extinction for all systems. For two systems we also found strong associated H$_2$ absorption $\log N(\text{H$_2$)[cm$^{-2}$]}=18.16\pm0.03$ and $19.28\pm0.06$ at $z=3.26$ and $2.25$ towards J2205+1021 and J2359+1354, respectively), while for the remaining five we measured conservative upper limits on the H$_2$ column densities of typically $\log N(\text{H$_2$)[cm$^{-2}$]}<17.3$. The increased H$_2$ detection rate ($10-55$% at 68% confidence level) at high HI column density compared to the overall damped Lyman-$\alpha$ population ($\sim 5-10$%) confirms previous works. We find that these seven ESDLAs have similar observed properties as those previously studied towards quasars and gamma-ray burst afterglows, suggesting they probe inner regions of galaxies. We use the abundance of ionised carbon in excited fine-structure level to calculate the cooling rates through the CII $\lambda$158$\mu$m emission, and compare them with the cooling rates from damped Lyman-$\alpha$ systems in the literature. We find that the cooling rates distribution of ESDLAs also presents the same bimodality as previously observed for the general (mostly lower HI column density) damped Lyman-$\alpha$ population.

Several conjectures have been put forward to explain the RRATs, the newest subclass of neutron stars, and their connections to other radio pulsars. This work discusses these conjectures in the context of the characteristic properties of the RRAT population. Contrary to expectations, it is seen that - a) the RRAT population is statistically un-correlated with the nulling pulsars, and b) the RRAT phenomenon is unlikely to be related to old age or death-line proximity. It is more likely that the special emission property of RRATs is a result of certain restructuring of their magnetic fields, resulting from processes like accretion or glitch.

In frame of an implementation of the cooperative program studying of asteroids between the Kharadze Abastumani Astrophysical Observatory and the Astronomical Institute of V.N. Karazin Kharkiv National University the observations for five main-belt asteroids were performed to obtain their magnitude-phase relations and other physical characteristics. Preliminary results of the photometrical observations for the large dark asteroid (1390) Abastumani are presented.

The accreting compact objects in most of ultraluminous X-ray sources (ULXs) are likely to be neutron stars rather than black holes as suggested by the recent detection of periodic pulsations from some of these sources located in neighboring galaxies and one ULX that has hitherto been discovered in our own galaxy. As a member of the ULX family, NGC 300 ULX1 is a new pulsating ULX (PULX) spinning up at substantially high rates compared with other PULXs. In this paper, the strength of the magnetic field on the surface of the neutron star is inferred from the energy of the cyclotron absorption line detected in the pulsed X-ray spectrum of NGC 300 ULX1 and the plausible ranges for the neutron-star mass and beaming fraction are estimated using the observed spin period and period derivative of the pulsar and the measured X-ray flux of the source. Our analysis favors proton cyclotron resonance scattering as a viable mechanism to account for both the observed cyclotron energy and high spin-up rates provided that the absorption line is generated close to the surface of the neutron star.

Primordial black holes of planetary masses captured by compact stars are widely studied to constrain their composition fraction of dark matter. Such a capture may lead to an inspiral process and be detected through gravitational wave signals. In this Letter, we study the post-capture inspiral process by considering two different kinds of compact stars, i.e., strange stars and neutron stars. The dynamical equations are numerically solved and the gravitational wave emission is calculated. It is found that the next generation gravitational wave detectors can detect the inspiraling of a $10^{-5}\,M_\odot$ primordial black hole at a distance of 1 kpc. A Jovian-mass case can even be detected at megaparsecs. Moreover, the kilohertz gravitational wave signal shows significant differences for strange stars and neutron stars, potentially making it a novel probe to the dense matter equation of state.

In this research, the simulation of lateral distribution function (LDF) of Cherenkov radiation was performed using CORSIKA software for two hadronic models QGSJET and GHEISHA. This simulation was performed for several elementary particles such as protons, iron nuclei, electrons and gamma quanta, in the range of energies 1-20 PeV for three zenith angles 0, 20 and 30 degrees. A parameterization of Cherenkov light LDF was performed for that simulated curves using Lorentzian function. The comparison between the obtained results for LDF of Cherenkov light with that measured with Yakutsk EAS array gave a good agreement within the distances of 100-1000 m from the shower axis.

The simulation of the extensive air showers was performed by investigating the longitudinal development parameters (N and Xmax) by using AIRES system version 19.04.0. The simulation was performed at the energy range (10^18-10^20 eV) for different primary particles (such as primary proton and iron nuclei) and different zenith angles. The longitudinal development curves of EAS are fitted using Gaussian function that gave a new parameters for different primary particles and different zenith angles at the energy range (10^18-10^20 eV).

The Orion molecular cloud complex harbours the nearest GMCs and site of high-mass star formation. Its YSO populations are thoroughly characterized. The region is therefore a prime target for the study of star formation. Here, we verify the performance of the SuperCAM 64 pixel heterodyne array on APEX. We give a descriptive overview of a set of wide-field CO(3-2) spectral cubes obtained towards the Orion GMC complex, aimed at characterizing the dynamics and structure of the extended molecular gas in diverse regions of the clouds, ranging from very active sites of clustered star formation in Orion B to comparatively quiet regions in southern Orion A. We present a 2.7 square degree (130pc$^2$) mapping survey in the CO(3-2) transition, obtained using SuperCAM on APEX at an angular resolution of 19'' (7600AU or 0.037pc at a distance of 400pc), covering L1622, NGC2071, NGC2068, OriB9, NGC2024, and NGC2023 in Orion B, and the southern part of the L1641 cloud in Orion A. We describe CO integrated emission and line moment maps and position-velocity diagrams and discuss a few sub-regions in some detail. Evidence for expanding bubbles is seen with lines splitting into double components, most prominently in NGC2024, where we argue that the bulk of the molecular gas is in the foreground of the HII region. High CO(3-2)/CO(1-0) line ratios reveal warm CO along the western edge of Orion B in the NGC2023/NGC2024 region facing the IC434 HII region. Multiple, well separated radial velocity components seen in L1641-S suggest that it consists of a sequence of clouds at increasingly larger distances. We find a small, spherical cloud - the 'Cow Nebula' globule - north of NGC2071. We trace high velocity line wings for the NGC2071-IR outflow and the NGC2024 CO jet. The protostellar dust core FIR4 (rather than FIR5) is the true driving source of the NGC2024 monopolar outflow.

The supermassive black holes in most galaxies in the universe are powered by hot accretion flows. Both theoretical analysis and numerical simulations have indicated that, depending on the degree of magnetization, black hole hot accretion flow is divided into two modes, namely SANE (standard and normal evolution) and MAD (magnetically arrested disk). It has been an important question which mode the hot accretion flows in individual sources should belong to in reality, SANE or MAD. This issue has been investigated in some previous works but they all suffer from various uncertainties. By using the measured rotation measure values in the prototype low-luminosity active galactic nuclei in {M87} at 2, 5, and 8 GHz along the jet at various distances from the black hole, combined with three dimensional general relativity magnetohydrodynamical numerical simulations of SANE and MAD, we show in this paper that the predicted rotation measure values by MAD are well consistent with observations, while the SANE model overestimates the rotation measure by over two orders of magnitude thus is ruled out.

We study the effect of cosmic ray (CR) acceleration in the massive compact star cluster Westerlund 1 in light of its recent detection in $\gamma$-rays. Recent observations reveal a $1/r$ radial distribution of the CR energy density. Here we theoretically investigate whether or not this profile can help to distinguish between (1) continuous CR acceleration in the star cluster stellar wind-driven shocks and (2) discrete CR acceleration in multiple supernovae shocks -- which are often debated in the literature. Using idealized two-fluid simulations and exploring different acceleration sites and diffusion coefficients, we obtain the CR energy density profile and luminosity to find the best match for the $\gamma$-ray observations. We find that the inferred CR energy density profiles from observations of $\gamma$-ray luminosity and mass can be much different from the true radial profile. CR acceleration at either the cluster core region or the wind termination shock can explain the observations, if the diffusion coefficient is $\kappa_{\rm cr}\sim 10^{27}$ cm$^2$ s$^{-1}$ and a fraction of $\approx 10\%-20\%$ of the shock power/post-shock pressure is deposited into the CR component. We also study the possibility of discrete supernovae (SN) explosions being responsible for CR acceleration and find that with an injection rate of 1 SN in every $\sim 0.03$ Myr, one can explain the observed $\gamma$-ray profile. This multiple SN scenario is consistent with X-ray observations only if the thermal conductivity is close to the Spitzer value.

Earlier studies have suggested that deviations from the local thermodynamic equilibrium (LTE) play a minor role in the formation of Be II 313 nm resonance lines in solar and stellar atmospheres. Recent improvements in the atomic data allow a more complete model atom of Be to be constructed and the validity of these claims to be reassessed using more up-to-date atomic physics. The main goal of this study therefore is to refocus on the role of non-local thermodynamic equilibrium (NLTE) effects in the formation of Be II 313.04 and 313.11 nm resonance lines in solar and stellar atmospheres. For this, we constructed a model atom of Be using new atomic data that recently became available. The model atom contains 98 levels and 383 radiative transitions of Be I and Be II and uses the most up-to-date collision rates with electrons and hydrogen. This makes it the most complete model atom used to determine 1D NLTE solar Be abundance and to study the role of NLTE effects in the formation of Be II 313 nm resonance lines. We find that deviations from LTE have a significant influence on the strength of the Be II 313 nm line in solar and stellar atmospheres. For the Sun, we obtained the 1D NLTE Be abundance of A(Be)_NLTE = 1.32+/-0.05, which is in excellent agreement with the meteoritic value of A(Be)=1.31+/-0.04. Importantly, we find that NLTE effects become significant in FGK stars. Moreover, there is a pronounced variation in 1D NLTE-LTE abundance corrections with the effective temperature and metallicity. Therefore, contrary to our previous understanding, the obtained results indicate that NLTE effects play an important role in Be line formation in stellar atmospheres and have to be properly taken into account in Be abundance studies, especially in metal-poor stars.

We investigate the origin of the extraplanar diffuse ionized gas (eDIG) and its predominant ionization mechanisms in five nearby (17-46 Mpc) low-mass ($10^9\text{-}10^{10}$ $M_{\odot}$) edge-on disk galaxies: ESO 157-49, ESO 469-15, ESO 544-27, IC 217, and IC 1553. We acquired Multi Unit Spectroscopic Explorer (MUSE) integral field spectroscopy and deep narrowband H$\alpha$ imaging of our sample galaxies. To investigate the connection between in-plane star formation and eDIG, we perform a photometric analysis of our narrowband H$\alpha$ imaging. We measure eDIG scale heights of $h_{z\text{eDIG}} = 0.59 \text{-} 1.39$ kpc and find a positive correlation between them and specific star formation rates. In all galaxies, we also find a strong correlation between extraplanar and midplane radial H$\alpha$ profiles. Using our MUSE data, we investigate the origin of eDIG via kinematics. We find ionized gas rotation velocity lags above the midplane with values between 10 and 27 km s$^{-1}$ kpc$^{-1}$. While we do find hints of an accretion origin for the ionized gas in ESO 157-49, IC 217, and IC 1553, overall the ionized gas kinematics of our galaxies do not match a steady galaxy model or any simplistic model of accretion or internal origin for the gas. We also construct standard diagnostic diagrams and emission-line maps (EW(H$\alpha$), [NII]/H$\alpha$, [SII]//H$\alpha$, [OIII]/H$\beta$) and find regions consistent with mixed OB star and hot low-mass evolved stars (HOLMES) ionization, and mixed OB-shock ionization. Our results suggest that OB stars are the primary driver of eDIG ionization, while both HOLMES and shocks may locally contribute to the ionization of eDIG to a significant degree. Despite our galaxies' similar structures and masses, we find a surprisingly composite image of ionization mechanisms and a multifarious origin for the eDIG.

Evidence for high-energy astrophysical PeV neutrinos has been found in the IceCube experiment from an analysis with a 7.5 year (2010 - 2017) data. Active galactic nuclei (AGN) are among the most prominent objects in the universe, and are widely speculated to be emitters of ultra-high-energy (UHE) cosmic rays with proton domination. Based on the standard two-step LLCD mechanism of particle acceleration, a transformation of energy occurs from AGN's central super-massive black hole (SMBH) rotation to high-energy protons. Protons can be accelerated up to $\sim 0.1$ EeV energies and above, and might generate PeV neutrinos in the energy range $1$--$10$~ PeV through plausible hadronic interactions. The theoretically estimated revised extragalactic diffuse muon neutrino flux employing the "luminosity-dependent density evolution (LDDE)" model for the AGN luminosity function (LF) is found consistent with the IceCube level if only a fraction, $6.56\%$ of the total bolometric luminosity (BL) of AGN is being realizable to power the PeV neutrinos. In the $\Lambda$~CDM cosmological framework with the LDDE modeled LF and photon index distribution, about $5.18\%$ of the total BL is enough to power the IceCube neutrinos.

We present a phenomenological analysis of current observational constraints on classes of FLRW cosmological models in which the matter side of Einstein's equations includes, in addition to the canonical term, a term proportional to some function of the energy-momentum tensor ($T^2=T_{\alpha\beta}T^{\alpha\beta}=\rho^2+3p^2$), or of its trace ($T=\rho-3p$). Qualitatively, one may think of these models as extensions of general relativity with a nonlinear matter Lagrangian. As such they are somewhat different from the usual dynamical dark energy or modified gravity models: in the former class of models one adds further dynamical degrees of freedom to the Lagrangian (often in the form of scalar fields), while in the latter the gravitational part of the Lagrangian is changed. We study both of these models under two different scenarios: (1) as phenomenological two-parameter or three-parameter extensions of the standard $\Lambda$CDM, in which case the model still has a cosmological constant but the nonlinear matter Lagrangian leads to additional terms in Einstein's equations, which cosmological observations tightly constrain, and (2) as alternatives to $\Lambda$CDM, where there is no cosmological constant, and the nonlinear matter term would have to provide the acceleration (which would be somewhat closer in spirit to the usual modified gravity models). A comparative analysis of the observational constraints obtained in the various cases provides some insight on the level of robustness of the $\Lambda$ model and on the parameter space still available for phenomenological alternatives.

Matter distribution in the environment of galaxy clusters, from their cores to their connected cosmic filaments, must be in principle related to the underlying cluster physics and it evolutionary state. We aim to investigate how radial and azimuthal distribution of gas is affected by cluster environments, and how it can be related to cluster mass assembly history. Radial physical properties of gas (velocity, temperature, and density) is first analysed around 415 galaxy cluster environments from IllustrisTNG simulation at z=0. Whereas hot plasma is virialised inside clusters (< R200), the dynamics of warm hot inter-galactic medium (WHIM) can be separated in two regimes: accumulating and slowly infalling gas at cluster peripheries (~ R200) and fast infalling motions outside clusters (> 1.5 R200). The azimuthal distribution of dark matter (DM), hot and warm gas phases is secondly statistically probed by decomposing their 2-D distribution in harmonic space. Inside clusters, the azimuthal symmetries of DM and hot gas are well tracing cluster structural properties, such as their center offsets, substructure fractions and elliptical shapes. Beyond cluster virialised regions, we found that WHIM gas follows the azimuthal distribution of DM by tracing cosmic filament patterns. Azimuthal symmetries of hot and warm gas distribution is finally shown to be imprints of cluster mass assembly history, by strongly correlating with the formation time, mass accretion rate, and dynamical state of clusters. Azimuthal mode decomposition of 2-D gas distribution is a promising probe to assess the 3-D physical and dynamical cluster properties up to their connected cosmic filaments.

We present herein a systematic X-ray analysis of blue galaxy-clusters at $z=0.84$ discovered by the Subaru telescope. The sample consisted of 43 clusters identified by combining red-sequence and blue-cloud surveys, covering a wide range of emitter fractions (i.e., 0.3--0.8). The spatial extent of the over-density region of emitter galaxies was approximately 1~Mpc in radius. The average cluster mass was estimated as $0.6 (<1.5)\times10^{14}~{\rm M_\odot}$ from the stacked weak-lensing measurement. We analyzed the XMM-Newton archival data, and measured the X-ray luminosity of the hot intracluster medium. As a result, diffuse X-ray emission was marginally detected in 14 clusters, yielding an average luminosity of $5\times 10^{42}~{\rm erg\,s^{-1}}$. On the contrary, it was not significant in 29 clusters. The blue clusters were significantly fainter than the red-dominated clusters, and the X-ray luminosity did not show any meaningful correlation either with emitter fraction or richness. The X-ray surface brightness was low, but the amount of gas mass was estimated to be comparable to that observed in the $10^{13-14}~{\rm M_{\odot}}$ cluster. Based on the results, we suggest that the blue clusters are at the early formation stage, and the gas is yet to be compressed and heated up to produce appreciable X-rays. Follow-up spectroscopic measurements are essential to clarify the dynamical status and co-evolution of galaxies and hot gas in the blue clusters.

Nanohertz gravitational waves (GWs) generated by the individual inspiraling supermassive black hole binaries (SMBHBs) in the centers of galaxies may be detected by pulsar timing arrays (PTAs) in the future. The GW signals from individual SMBHBs can be employed as standard sirens to measure absolute cosmic distances, and further provide constraints on cosmological parameters via the distance--redshift relation. Namely, these SMBHBs with known redshifts can serve as bright sirens. In this paper, we analyze the ability of SKA-era PTAs to detect existing SMBHB candidates by means of simulating the timing residuals of pulsar signals, and use the mock data of SMBHB bright sirens to perform cosmological parameter estimations. We find that once the root mean square of timing residuals ($\sigma_t$) could be reduced to 20 ns, about only 100 millisecond pulsars are needed to provide a tight constraint on the Hubble constant ($H_0$) with the precision close to $1\%$, which meets the criterion of precision cosmology. We further show that the SMBHB bright sirens can effectively break the cosmological parameter degeneracies inherent in the CMB data, and thus improve the constraint precision of the equation of state of dark energy ($w$) to $3.5\%$, which is comparable with the result of Planck 2018 TT,TE,EE+lowE+lensing+SNe+BAO. We conclude that the bright sirens from SKA-era PTAs will play an important role in breaking the cosmological parameter degeneracies and exploring the nature of dark energy.

Platelets corrugated horn is a promising technology for their scalability to a large corrugated horn array. In this paper, we present the design, fabrication, measurement and uncertainty analysis of a wideband 170-320 GHz platelet corrugated horn that features with low sidelobe across the band (<-30 dB). We also propose an accurate and universal method to analyze the axial misalignment of the platelets for the first time. It is based on the mode matching (MM) method with a closed-form solution to off-axis circular waveguide discontinuities obtained by using Graf addition theorem for the Bessel functions. The uncertainties introduced in the fabrication have been quantitatively analyzed using the Monte Carlo method. The analysis shows the cross-polarization of the corrugated horn degrades significantly with the axial misalignment. It well explains the discrepancy between the designed and the measured cross-polarization of platelets corrugated horn fabricated in THz band. The method can be used to determine the fabrication tolerance needed for other THz corrugated horns and evaluate the impact of the corrugated horn for astronomical observations.

Single-pulse studies are important to understand the pulsar emission mechanism and the noise floor in precision timing. We study total intensity and polarimetry properties of three bright millisecond pulsars - PSRs J1022+1001, J1713+0747, and B1855+09 - that have detectable single pulses at multiple frequencies. We report for the first time the detection of single pulses from PSR J1022+1001 and PSR J1713+0747 at 4.5 GHz. In addition, for those two pulsars the fraction of linear polarization in the average profile is significantly reduced at 4.5 GHz, compared to 1.38 GHz, which could support the expected deviation from a dipolar field closer to the pulsar surface. While the fraction of linear polarization in single pulses of PSR J1022+1001 also follows this trend, the few single pulses in our data set of J1713+0747 seem to show a higher fraction of linear polarization at 4.5 GHz than at 1.38 GHz. We do not find evidence for orthogonal modes in single pulses for any of the pulsars. More sensitive multi-frequency observations may be useful to confirm these findings. The jitter noise contributions at 1.38 GHz, scaled to one hour, for PSR J1022+1001 and PSR J1713+0747 are ~220 ns and ~41 ns respectively and are consistent with previous studies, but for B1855+09 the jitter noise is ~86 ns, which is lower than previously measured. We also show that selective bright-pulse timing of PSR J1022+1001 yields improved root-mean-square residuals of ~23 mus, which is a factor of ~4 better than timing using single pulses alone.

Only recently, OD, the deuterated isotopolog of hydroxyl, OH, has become accessible in the interstellar medium; spectral lines from both species have been observed in the supra-Terahertz and far infrared regime. Here we study rotational lines of OD and OH towards 13 Galactic high-mass star forming regions, with the aim to constrain the OD abundance and infer the deuterium fractionation of OH in their molecular envelopes. We used the Stratospheric Observatory for Infrared Astronomy (SOFIA) to observe the $^2\Pi_{3/2}$ $J=5/2-3/2$ ground-state transition of OD at 1.3 THz ($215~\mu$m) and the rotationally excited OH line at 1.84 THz ($163~\mu$m). We also used published high-spectral-resolution SOFIA data of the OH ground-state transition at 2.51 THz ($119.3~\mu$m). Our results show that absorption from the $^2\Pi_{3/2}$ OD $J=5/2-3/2$ ground-state transition is prevalent in the dense clumps surrounding active sites of high-mass star formation. We performed detailed radiative transfer modelling to investigate the OD abundance profile in the inner envelope for a large fraction of our sample. Our modelling suggests that part of the absorption arises from the denser inner parts, while the bulk of it as seen with SOFIA originates in the outer, cold layers of the envelope for which our constraints on the molecular abundance suggest a strong enhancement in deuterium fractionation. We find a weak negative correlation between the OD abundance and the bolometric luminosity to mass ratio, an evolutionary indicator, suggesting a slow decrease of OD abundance with time. A comparison with HDO shows a similarly high deuterium fractionation for the two species in the cold envelopes, which is of the order of 0.48% for the best studied source, G34.26+0.15. Our results are consistent with chemical models that favour rapid exchange reactions to form OD in the dense cold gas.

We analyze spectral properties of solar convection in the range of depths from 0 to 19~Mm using subsurface flow maps obtained by the time-distance heiioseismology analysis of solar-oscillation data from the Helioseismic and Magnetic Imager (HMI) onboard Solar Dynamics Observatory (SDO) from May 2010 to September 2020. The results reveal a rapid increase of the horizontal flow scales with the depth, from supergranulation to giant-cell scales, and support the evidence of large-scale convection, previously detected by tracking the motion of supergranular cells on the surface. The total power of convective flows correlates with the solar activity cycle. During the solar maximum, the total power decreases in shallow subsurface layers and increases in the deeper layers.

We report the detection of strong thermal spectral component in the short Gamma-Ray Burst 170206A with three intensive pulses in its lightcurves, throughout which the fluxes of this thermal component exhibit fast temporal variability same as that of the accompanied non-thermal component. The values of the time-resolved low-energy photon index in the non-thermal component are between about -0.78 and -0.17, most of which are harder than the line-of-death (-2/3) of the synchrotron emission process. In addition, we found the plausible common evolution between the thermal component and the non-thermal component, indicating a positive correlation between the fluxes as well as peak energies of both components. Based on the observations, we explore the possible origin of the thermal component and the non-thermal component based on the models of one common radiation region and two radiation regions.

The evaluation of sky characteristics plays a fundamental role for many astrophysical experiments and ground-based observations. In solar physics, the main requirement for such observations is a very low sky brightness value, that is, less than one-millionth of the solar disk brightness. Few places match such requirement for ground-based, out-of-eclipse coronagraphic measurements. A candidate coronagraphic site is the Dome C plateau in Antarctica. In this paper, we show the first results of the sky brightness measurements at Dome C with the Extreme Solar Coronagraphy Antarctic Program Experiment (ESCAPE) at Italian-French Concordia Station, on Dome C, Antarctica (3300 m) during the summer XXXIV and XXXV Expeditions of the "Italian Piano Nazionale Ricerche Antartiche" (PNRA). The sky brightness measurements were carried out with the internally-occulted Antarctic coronagraph - AntarctiCor. In optimal atmospheric conditions, the sky brightness of Dome C has reached values of the order of 1.0 to 0.7 millionth of of the solar disk brightness.

During its commissioning phase in 2020, the Spectrometer/Telescope for Imaging X-rays (STIX) on board the Solar Orbiter spacecraft observed 69 microflares. The two most significant events from this set (of GOES class B2 and B6) were observed on-disk from the spacecraft as well as from Earth and analysed in terms of the spatial, temporal, and spectral characteristics. We complement the observations from the STIX instrument with EUV imagery from SDO/AIA and GOES soft X-ray data by adding imaging and plasma diagnostics over different temperature ranges for a detailed microflare case study that is focussed on energy release and transport. Spectral fitting of the STIX data shows clear nonthermal emission for both microflares studied here. The deduced plasma parameters from DEM reconstruction as well as spectral fitting roughly agree with the values in the literature for microflares as do the nonthermal fit parameters from STIX. The observed Neupert effects and impulsive and gradual phases indicate that both events covered in this study are consistent with the standard chromospheric evaporation flare scenario. For the B6 event on 7 June 2020, this interpretation is further supported by the temporal evolution seen in the DEM profiles of the flare ribbons and loops. For this event, we also find that accelerated electrons can roughly account for the required thermal energy. The 6 June 2020 event demonstrates that STIX can detect nonthermal emission for GOES class B2 events that is nonetheless smaller than the background rate level. We demonstrate for the first time how detailed multi-instrument studies of solar flares can be performed with STIX.

Using parallaxes from Gaia Early Data Release 3 (EDR3), we determine multi-wavelength BVI, JHK, and [3.6] & [4.5] micron absolute magnitudes for 37 nearby Milky Way Cepheids, covering the period range between 5 and 60 days. We apply these period-luminosity relations to Cepheids in the Large and Small Magellanic Clouds, and find that the derived distances are significantly discrepant with the geometric distances according to detached eclipsing binaries (DEBs). We explore several potential causes of these issues, including reddening, metallicity, and the existence of an additional zero-point offset, but none provide a sufficient reconciliation with both DEB distances. We conclude that the combination of the systematic uncertainties on the EDR3 parallaxes with the uncertainties on the effect of metallicity on the Cepheid distance scale leads to a systematic error floor of approximately 3%. We therefore find that the EDR3 data is not sufficiently accurate in the regime of these bright Cepheids to determine extragalactic distances precise to the 1% level at this time, in agreement with a number of contemporary studies.

Kinematic lensing (KL) is a new cosmological measurement technique that combines traditional weak lensing (WL) shape measurements of disc galaxies with their kinematic information. Using the Tully-Fisher relation KL breaks the degeneracy between intrinsic and observed ellipticity and significantly reduces the impact of multiple systematics that are present in traditional WL. We explore the performance of KL given the instrument capabilities of the $\textit{Roman Space Telescope}$, assuming overlap of the High Latitude Imaging Survey (HLIS), the High Latitude Spectroscopy Survey (HLSS) over 2,000 deg$^2$. Our KL suitable galaxy sample has a number density of $n_{\mathrm{gal}}=4~\mathrm{arcmin}^{-1}$ with an estimated shape noise level of $\sigma_{\epsilon}=0.035$. We quantify the cosmological constraining power on $\Omega_{\mathrm{m}}$-$S_8$, $w_p$-$w_a$ by running simulated likelihood analyses that account for redshift and shear calibration uncertainties, intrinsic alignment and baryonic feedback. Compared to a traditional WL survey we find that KL significantly improves the constraining power on $\Omega_{\mathrm{m}}$-$S_8$ (FoM$_{\mathrm{KL}}$=1.70FoM$_{\mathrm{WL}}$) and $w_p$-$w_a$ (FoM$_{\mathrm{KL}}$=3.65FoM$_{\mathrm{WL}}$). We also explore a "narrow tomography KL survey" using 30 instead of the default 10 tomographic bins, however we find no meaningful enhancement to the FoM even when assuming a significant time-dependence in our fiducial dark energy input scenarios.

A good knowledge of the angular diameters of stars used to calibrate the observables in stellar interferometry is fundamental. As the available precision for giant stars is worse than the required per cent level, we aim to improve the knowledge of many diameters using MATISSE (Multiple AperTure mid-Infrared SpectroScopic Experiment) data in its different instrumental configurations. Using the squared visibility MATISSE observable, we compute the angular diameter value, which ensures the best-fitting curves, assuming an intensity distribution of a uniform disc. We take into account that the transfer function varies over the wavelength and is different from one instrumental configuration to another. The uncertainties on the diameters are estimated using the residual bootstrap method. Using the low spectral resolution mode in the L band, we observed a set of 35 potential calibrators selected in the Mid-infrared stellar Diameter and Flux Compilation Catalogue with diameters ranging from about 1 to 3 mas. We reach a precision on the diameter estimates in the range 0.6 per cent to 4.1 per cent. The study of the stability of the transfer function in visibility over two nights makes us confident in our results. In addition, we identify one star, 75 Vir initially present in the calibrator lists, for which our method does not converge, and prove to be a binary star. This leads us to the conclusion that our method is actually necessary to improve the quality of the astrophysical results obtained with MATISSE, and that it can be used as a useful tool for 'bad calibrator' detection.

Cosmic microwave background radiation can supply us some most significant parts of the information on the universe. Some researchers believe that The gravitational system cannot be decribed by the standard statistical mechanics. In this article we apply Tsallis nonextensive statistical mechanics to investigate CMB spectrum and related cosmological processes. Based on recent observational data we find that the nonextensive statistical mechanics can modify the values of related physical quantites. Since the value of physical quantites have changed, some processes, such as recombination, can be affected. We have investigated the anisotropy of the CMB for two effects: the dipole anisotropy of CMB and the Sunyaev-Zel'dovich effect. We find that the dipole anisotropy of CMB cannot be modified by the nonextensive statistical mechanics. However, the standard result of the Sunyaev-Zel'dovich effect should be modified by nonextensive statistical mechanics. In principle, future work can distinguish these effects.

Archival [3.6] and [4.5] images are used to identify and characterize variable stars in the Magellanic-type galaxies Holmberg II, NGC 2366, and IC 2574. Using parametric and non-parametric detection methods, 74 confirmed or suspected long period variables (LPVs) are found. The period distributions of the LPVs in NGC 2366 and IC 2574 are similar. While the period distribution of LPVs in Ho II is uncertain due to small number statistics there appears to be a deficiency of LPVs with periods between 550 and 650 days when compared with NGC 2366 and IC 2574. The LPVs are diffusely distributed on the sky, and do not follow the underlying light from unresolved stars, as expected if episodes of star formation within the past few hundred Myr have occurred throughout the galaxies, including their outer regions. Distances computed for Ho II and NGC 2366 from the period-luminosity relations (PLRs) agree to within ~0.1 magnitudes with those based on the tip of the red giant branch (RGB). Efforts to estimate an LPV-based distance modulus for IC 2574 are complicated by the presence of first overtone pulsators among LPVs with periods ~ 600 days, although the PLR at the long period end is consistent with the distance estimated from the RGB-tip. In addition to the LPVs, 10 candidate sgB[e] or luminous blue variables and 2 candidate red supergiant variables are also identified. Nine candidate sgB[e] stars that do not show evidence of variability are also identified based on their locations in the color-magnitude diagram.

The PTAL (Planetary Terrestrial Analogues Library) project aims at building and exploiting a database involving several analytical techniques, to help characterizing the mineralogical evolution of terrestrial bodies, starting with Mars. Around 100 natural Earth rock samples have been collected from selected locations to gather a variety of analogues for Martian geology, from volcanic to sedimentary origin with different levels of alteration. All samples are to be characterized within the PTAL project with different mineralogical and elemental analysis techniques, including techniques brought on actual and future instruments at the surface of Mars (Near InfraRed spectroscopy, Raman spectroscopy and Laser Induced Breakdown Spectroscopy). This paper presents the NIR measurements and interpretations acquired with the ExoMars MicrOmega spare instrument. MicrOmega is a NIR hyperspectral microscope, mounted in the analytical laboratory of the ExoMars rover Rosalind Franklin. All PTAL samples have been observed at least once with MicrOmega using a dedicated setup. For all PTAL samples data description and interpretation are presented. For some chosen examples, RGB images and spectra are presented a well. A comparison with characterizations by NIR and Raman spectrometry is discussed for some of the samples. In particular, the spectral imaging capacity of MicrOmega allows detections of mineral components and potential organic molecules that were not possible with other one-spot techniques. Additionally, it enables to estimate heterogeneities in the spatial distribution of various mineral species. The MicrOmega/PTAL data shall support the future observations and analyses performed by MicrOmega/Rosalind Franklin instrument.

The spin axes of the mantle, fluid core and solid inner core of the Moon precess at frequency $\Omega_p=2\pi/18.6$ yr$^{-1}$ though with different orientations, leading to viscous friction at the core-mantle boundary (CMB) and inner core boundary (ICB). Here, we use a rotational model of the Moon with a range of inner core and outer core radii to investigate the relative importance of viscous dissipation at the CMB and ICB, and to show how this dissipation is connected to the phase lead angle ($\phi_p$) of the mantle ahead of its Cassini state. We show that when the inner core radius is $>80$ km and the free inner core nutation frequency $\Omega_{ficn}$ approaches $\Omega_p$, viscous dissipation at the ICB can be comparable to that at the CMB, and in the most extreme cases exceed it by as much as a factor 10. If so, the viscous dissipation in the lunar core projected back in time depends on how $\Omega_{ficn}$ has evolved relative to $\Omega_p$. We further show that constraints on the CMB and ICB radii of the lunar core can in principle be extracted by matching the observed phase lead of $\phi_p=0.27$ arcsec; this requires an improved estimate of tidal dissipation and an accurate model of the turbulent viscous torque. Lastly, when our rotational model is constrained to match $\phi_p=0.27$ arcsec, our results suggest that the viscous dissipation at the ICB is likely insufficient to have ever been above the threshold to power a thermally driven dynamo.

Mass accretion from the circumstellar disk onto the protostar is a fundamental process during star formation. Measuring the mass accretion rate is particularly challenging for stars belonging to binary systems, because it is often difficult to discriminate which component is accreting. DQ~Tau is an almost equal-mass spectroscopic binary system where the components orbit each other every 15.8~days. The system is known to display pulsed accretion, i.e., the periodic modulation of the accretion by the components on eccentric orbit. We present multiepoch ESO/VLT X-Shooter observations of DQ~Tau, with the aim to determine which component of this system is the main accreting source. We use the absorption lines in the spectra to determine the radial velocity of the two components, and measure the continuum veiling as a function of wavelength and time. We fit the observed spectra with non-accreting templates to correct for the photospheric and chromospheric contribution. In the corrected spectra we study in details the profiles of the emission lines and calculate mass accretion rates for the system as a function of orbital phase. In accordance with previous findings, we detect elevated accretion close to periastron. We measure that the accretion rate varies between $10^{-8.5}$ and $10^{-7.3}$ Msun/yr. The emission line profiles suggest that both stars are actively accreting and the dominant accretor is not always the same component, varying in few orbits.

Analyses of Lunar Laser Ranging data show that the spin-symmetry axis of the Moon is ahead of its expected Cassini state by an angle of $\phi_p$ = 0.27 arcsec. This indicates the presence of one or more dissipation mechanisms acting on the lunar rotation. A combination of solid-body tides and viscous core-mantle coupling have been proposed in previous studies. Here, we investigate whether viscoelastic deformation within a solid inner core at the centre of the Moon can also account for a part of the observed phase lead angle $\phi_p$. We build a rotational dynamic model of the Cassini state of the Moon that comprises an inner core, a fluid core and a mantle, and where solid regions are allowed to deform viscoelastically in response to an applied forcing. We show that the presence of an inner core does not change the global monthly Q of the Moon and hence, that the contribution from solid-body tides to $\phi_p$ is largely unaffected by an inner core. However, we also show that viscoelastic deformation within the inner core, acting to realign its figure axis with that of the mantle, can contribute significantly to $\phi_p$ through inner core-mantle gravitational coupling. We show that the contribution to $\phi_p$ is largest when the inner core viscosity is in the range of $10^{13}$ to $10^{14}$ Pa s, when the inner core radius is large and when the free inner core nutation frequency approaches a resonance with the precession frequency of $2\pi/18.6$ yr$^{-1}$.

The Extreme Universe Space Observatory on a Super Pressure Balloon 2 (EUSO-SPB2), in preparation, aims to make the first observations of Ultra-High Energy Cosmic Rays (UHECRs) from near space using optical techniques. EUSO-SPB2 will prototype instrumentation for future satellite-based missions, including the Probe of Extreme Multi-Messenger Astrophysics (POEMMA) and K-EUSO. The payload will consist of two telescopes. The first is a Cherenkov telescope (CT) being developed to quantify the background for future below-the-limb very high energy (E>10 PeV) astrophysical neutrino observations, and the second is a fluorescence telescope (FT) being developed for detection of UHECRs. The FT will consist of a Schmidt telescope, and a 6192 pixel ultraviolet camera with an integration time of 1.05 microseconds. The first step in the data acquisition process for the FT is a hardware level trigger in order to decide which data to record. In order to maximize the number of UHECR induced extensive air showers (EASs) which can be detected, a novel trigger algorithm has been developed based on the intricacies and limitations of the detector. The expected performance of the trigger has been characterized by simulations and, pending hardware verification, shows that EUSO-SPB2 is well positioned to attempt the first near-space observation of UHECRs via optical techniques.

The Cassini state equilibrium associated with the precession of the Moon predicts that the mantle, fluid core and solid inner core precess at different angles. We present estimates of the dissipation from viscous friction associated with the differential precession at the core-mantle boundary (CMB), $Q_{cmb}$, and at the inner core boundary (ICB), $Q_{icb}$, as a function of the evolving lunar orbit. We focus on the latter and show that, provided the inner core was larger than 100 km, $Q_{icb}$ may have been as high as $10^{10}-10^{11}$ W for most of the lunar history for a broad range of core density models. This is larger than the power required to maintain the fluid core in an adiabatic state, therefore the heat released by the differential precession at the ICB can drive a past lunar dynamo by thermal convection. This dynamo can outlive the dynamo from precession at the CMB and may have shutoff only relatively recently. Estimates of the magnetic field strength at the lunar surface are of the order of a few $\mu$T, compatible with the lunar paleomagnetic intensities recorded after 3 Ga. We further show that it is possible that a transition of the Cassini state associated with the inner core may have occurred as a result of the evolution of the lunar orbit. The heat flux associated with $Q_{icb}$ can be of the order of a few mW m$^{-2}$, which should slow down inner core growth and be included in thermal evolution models of the lunar core.

We present a model of the precession dynamics of the Moon that comprises a fluid outer core and a solid inner core. We show that three Cassini states associated with the inner core exist. The tilt angle of the inner core in each of these states is determined by the ratio between the free inner core nutation frequency ($\omega_{ficn}$) and the precession frequency $\Omega_p = 2\pi/18.6$ yr $^{-1}$. All three Cassini states are possible if $|\omega_{ficn}| > 2\pi/16.4$ yr $^{-1}$, but only one is possible otherwise. Assuming that the lowest energy state is favoured, this transition marks a discontinuity in the tilt angle of the inner core, transiting from $-33^\circ$ to $17^\circ$ as measured with respect to the mantle figure axis, where negative angles indicate a tilt towards the orbit normal. Possible Lunar interior density structures cover a range of $\omega_{ficn}$, from approximately half to twice as large as $\Omega_p$, so the precise tilt angle of the inner core remains unknown, though it is likely large because $\Omega_p$ is within the resonant band of $\omega_{ficn}$. Adopting one specific density model, we suggest an inner core tilt of approximately $-17^\circ$. Viscoelastic deformations within the inner core and melt and growth at the surface of a tilted inner core, both neglected in our model, should reduce this amplitude. If the inner core is larger than approximately 200 km, it may contribute by as much as a few thousandths of a degree on the observed mantle precession angle of $1.543^\circ$.

Measuring the quasi-normal mode~(QNM) spectrum emitted by a perturbed black-hole~(BH) --~also known as BH spectroscopy~-- provides an excellent opportunity to test the predictions of general relativity in the strong-gravity regime. We investigate the prospects and precision of BH spectroscopy in massive binary black hole ringdowns, one of the primary science objectives of the future Laser Interferometric Space Antenna~(LISA) mission. We simulate various massive binary BH population models, featuring competing prescriptions for the Delays between galaxy and BH mergers, for the impact of supernova feedback on massive BH growth, and for the initial population of high redshift BH seeds (light versus heavy seeds). For each of these scenarios, we compute the average number of expected events for precision BH spectroscopy using a Fisher-matrix analysis. We find that, for any heavy seed scenario, LISA will measure the dominant mode frequency within ${\cal O}(0.1) \%$ relative uncertainty and will estimate at least 3 independent QNM parameters within $1 \%$ error. The most optimistic heavy seed scenarios produce $\mathcal{O}(100)$ events with $1 \%$ measurability for 3 or more QNM quantities during LISA's operational time. On the other hand, light seed scenarios produce lighter merger remnants, which ring at frequencies higher than LISA's sensitivity. Interestingly, the light seed models give rise to a fraction of mergers in the band of Einstein Telescope, allowing for the measurement of 3 QNM parameters with $\sim 10 \%$ relative errors in approximately a few to ten events/yr. More precise BH spectroscopy in the light seed scenarios would require instruments operating in the deciHertz band.

The existence of black hole horizons has not been strictly proven observationally, and indeed it may not be possible to do so. However, alternatives may be established by the observation of gravitational wave echoes that probe possible near-horizon structure. These echoes are proposed to be generated in exotic compact objects that are horizonless and feature a partially reflecting "wall" inside their light rings, creating a cavity in which gravitational perturbations may echo, while leaking out through the angular momentum barrier with each pass. The characteristic signature of echoes is a comb of nearly evenly spaced spectral resonances. While approximately true, deviations from this simple picture can lead to severe observational signal losses. In this paper, we explore such subtleties with the latest results for echo sourcing and geometry. A physically motivated echo model is then developed as a sum over Lorentzian spectral lines, parametrized by functions of the horizon frame frequency and the size of the cavity. Our final spectrum is a function of only the mass and spin of the black hole, as well as the UV scale of the near-horizon physics. We then apply this model in a search for echoes in the gravitational wave event with the loudest ringdown signal in LIGO/Virgo, i.e. GW190521. We interpret our findings as a measurement of the fractional energy in post-merger echoes equal to $E_{echoes} / E_{GR} = 8.9 \pm 4.5\%$, where the uncertainty range represents the 90% credible region. The robustness of this result is tested against noise backgrounds and simulated injections, and we find that a signal persists through modifications to the model and changes in the data search.

Being the most massive binary black hole merger event observed to date, GW190521 is in a class of its own. The exceptionally loud ringdown of this merger makes it an ideal candidate to search for gravitational wave echoes, a proposed smoking gun for the quantum structure of black hole horizons. We perform an unprecedented multi-pronged search for echoes via two well-established and independent pipelines: a template-based search for stimulated emission of Hawking radiation, or Boltzmann echoes, and the model-agnostic coherent WaveBurst (cWB) search. Stimulated Hawking radiation from the merger is expected to lead to post-merger echoes at horizon mode frequency of $\sim 50$ Hz (for quadrupolar gravitational radiation), repeating at intervals of $\sim 1$ second, due to partial reflection off Planckian quantum structure of the horizon. A careful analysis using dynamic nested sampling yields a Bayesian evidence of $ 7\pm 2$ (90% confidence level) for this signal following GW190521, carrying an excess of $10^{+9}_{-7}\%$ in gravitational wave energy, relative to the main event. Similarly, the reconstructed waveform of the first echo in cWB carries an energy excess of $13^{+16}_{-7}\%$. Accounting for the "look-elsewhere" effects, we estimate a p-value for false detection probability of $5.1 \times 10^{-3}$ (or 2.6$\sigma$) using cWB pipeline, although the verdict on the co-localization of the post-merger echo and the main event in the sky is inconclusive. While the current evidence for stimulated Hawking radiation does not reach the gold standard of $5\sigma$, our findings are in line with expectations for stimulated Hawking radiation at current detector sensitivities. The next generation of gravitational wave observatories can thus draw a definitive conclusion on the quantum nature of black hole horizons.

Intraplate volcanic islands are often considered as stable relief with constant vertical motion and used for relative sea-level curves reconstruction. This study shows that large landslides cause non-negligible isostatic adjustment. The vertical motion that occurred after landslide is quantified using a modelling approach. We show that a giant landslide caused a coastline uplift of 80-110 m for an elastic thickness of 15 km < $T_e$ < 20 km in Tahiti. Theoretical cases also reveal that a coastal motion of 1 m occurred for a landslide involving a displaced volume of 0.2 $km^3$ and influence relative sea-level reconstruction. In Tahiti, a change in the subsidence rate of 0.1 mm/yr (from 0.25 mm/yr to 0.15 mm/yr) occurred during the last 6 kyr and could be explained by a landslide involving a minimum volume of 0.2 $km^3$, $6 \pm 1$ kyr ago.

Cosmic ray muons have been considered as a non-conventional radiation probe in various applications. To utilize cosmic ray muons in engineering applications, two important quantities, trajectory and momentum, must be known. The muon trajectories are easily reconstructed using two-fold detector arrays with a high spatial resolution. However, precise measurement of muon momentum is difficult to be achieved without deploying large and expensive spectrometers such as solenoid magnets. Here, we propose a new method to estimate muon momentum using multi-layer pressurized gas Cherenkov radiators. This is accurate, portable, compact (< 1m3), and easily coupled with existing muon detectors without the need of neither bulky magnetic nor time-of-flight spectrometers. The results show that not only our new muon spectrometer can measure muon momentum with a resolution of +-0.5 GeV/c in a momentum range of 0.1 to 10.0 GeV/c, but also we can reconstruct cosmic muon spectrum with high accuracy (~90%).

We study a family of equations of state (EoS) for hybrid neutron star (NS) matter. Our hybrid EoS are based on an instantaneous nonlocal version of NJL model in $SU(2)_f$ with vector interactions and color superconductivity describing the quark matter (QM) with a Maxwell construction phase transition to hadronic matter phase described by the "DD2" EoS with excluded volume and a crust at low baryonic densities. The form factor in the QM nonlocal model was fitted to lattice QCD (LQCD) results in the Coulomb gauge. To simultaneously fulfill constraints from NICER and tidal deformability from GW170817 it is necessary to consider a $\mu$ dependent bag constant that mimics confinement. Our results show an asymptotic constant behaviour for the speed of sound that reproduces the conjectured value of $1/3$ from QCD in the free case, and larger constant values in the range 0.4 - 0.6 when interactions are turned on.

In a recent work [Bret, EPL \textbf{135} (2021) 35001], quantum electrodynamic (QED) effects were evaluated for the two-stream instability. It pertains to the growth of perturbations with a wave vector oriented along the flow in a collisionless counter-streaming system. Here, the analysis is extended to every possible orientation of the wave vector. The previous result for the two-stream instability is recovered, and it is proved that even within the framework of a 3D analysis, this instability remains fundamentally 1D even when accounting for QED effects. The filamentation instability, found for wave vectors normal to the flow, is weakly affected by QED corrections. As in the classical case, its growth rate saturates at large $k_\perp$. The saturation value is found independent of QED corrections. Also, the smallest unstable $k_\perp$ is independent of QED corrections. Surprisingly, unstable modes found for oblique wave vectors do \emph{not} follow the same pattern. For some, QED corrections do reduce the growth rate. But for others, the same corrections increase the growth rate instead. The possibility for QED effects to play a role in un-magnetized systems is evaluated. Pair production resulting from gamma emission by particles oscillating in the exponentially growing fields, is not accounting for.

First-order phase transitions (FOPTs) are ubiquitous in physics beyond the Standard Model (SM). Recently, models with no dimensionful parameters in the tree-level action have been attracting much attention because they can predict a very strong FOPT with ultra-supercooling. In this paper, we study the cosmological signatures of such a supercooling model. As a concrete model, we consider the SM with two additional real scalars $\phi$ and $S$, which can realize the electroweak symmetry breaking via Coleman-Weinberg mechanism. One of the additional scalars $S$ can naturally become a Dark Matter (DM) candidate due to the $Z_2^{}$ symmetry of the action. We study the FOPT of this model and calculate the Gravitational Wave (GW) signals and the thermal relic abundance of $S$. As a result, we find that (i) the GW peak amplitude can reach $\sim 10^{-10}$ around the frequency $f\sim 10^{-3}~$Hz for some model parameters and (ii) the scalar mixing coupling $\lambda_{\phi S}^{}$ is constrained to be $0.8\lesssim \lambda_{\phi S}^{}\lesssim 1$ in order to achieve both FOPT and DM relic abundance.

We examine the gravitational wave frequencies from neutron stars during thermal evolution, adopting the relativistic Cowling approximation. We particularly focus on the neutron star models, in which the direct Urca (rapid cooling process) does not work, without the superfluidity and superconductivity. For such models, the cooling curve hardly depends on the equation of state (EOS) as well as the mass of neutron star, while we show that the gravitational wave frequencies strongly depend on the both properties. Then, we find that the frequencies of the fundamental and the 1st pressure mode multiplied with the stellar mass are well expressed as a function of the stellar compactness almost independently of the EOS. We also find that the frequency of the 1st gravity mode in later phase of the thermal evolution is strongly correlated with the stellar compactness. In addition, we derive the empirical formula estimating the threshold mass for the onset of the direct Urca inside the neutron star as a function of the nuclear saturation parameter. This formula will give us a constraint on the neutron star properties, if it would be observationally found that the direct Urca occurs (or does not work) inside the neutron star.

Solitons in space-time capable of transporting time-like observers at superluminal speeds have long been tied to violations of the weak, strong, and dominant energy conditions of general relativity. This trend was recently broken by a new approach that identified soliton solutions capable of superluminal travel while being sourced by purely positive energy densities. This is the first example of hyper-fast solitons satisfying the weak energy condition, reopening the discussion of superluminal mechanisms rooted in conventional physics. This article summarizes the recent finding and its context in the literature. Remaining challenges to autonomous superluminal travel, such as the dominant energy condition, horizons, and the identification of a creation mechanism are also discussed.

We study cosmological effects of explicit Peccei-Quinn breaking on the QCD axion dark matter. We find that the axion abundance decreases or increases significantly depending on the initial position, even for a tiny Peccei-Quinn breaking that satisfies the experimental bound of the neutron electric dipole measurements. If the axion first starts to oscillate around a wrong vacuum and if it gets trapped there until the false vacuum disappears due to non-perturbative QCD effects, its abundance increases significantly and is independent of the decay constant $f_a$, as first pointed out in [JHEP 06 (2016) 150]. Thus, the axion produced by the trapping mechanism can explain dark matter even when the decay constant is close to the lower limit due to stellar cooling arguments. On the other hand, if the axion starts to oscillate about a potential minimum close to the low-energy vacuum, its abundance is significantly reduced because of the adiabatic suppression mechanism. This relaxes the upper limit of the axion window to large values of $f_a$. We also discuss how the axionic isocurvature perturbation is affected by the Peccei-Quinn breaking term, and show that it can be suppressed in both regimes. In particular, the isocurvature bound on the inflation scale is relaxed by many orders of magnitudes for $f_a \lesssim 10^{11}{\rm GeV}$ compared to the conventional scenario.

We present results of an all-sky search for continuous gravitational waves which can be produced by spinning neutron stars with an asymmetry around their rotation axis, using data from the third observing run of the Advanced LIGO and Advanced Virgo detectors. Four different analysis methods are used to search in a gravitational-wave frequency band from 10 to 2048 Hz and a first frequency derivative from $-10^{-8}$ to $10^{-9}$ Hz/s. No statistically-significant periodic gravitational-wave signal is observed by any of the four searches. As a result, upper limits on the gravitational-wave strain amplitude $h_0$ are calculated. The best upper limits are obtained in the frequency range of 100 to 200 Hz and they are ${\sim}1.1\times10^{-25}$ at 95\% confidence-level. The minimum upper limit of $1.10\times10^{-25}$ is achieved at a frequency 111.5 Hz. We also place constraints on the rates and abundances of nearby planetary- and asteroid-mass primordial black holes that could give rise to continuous gravitational-wave signals.