Papers for Monday, May 29 2023

The aim of "search for extraterrestrial intelligence" (SETI) commensal surveys is to scan the sky to find possible technosignatures from the extraterrestrial intelligence (ETI). The mitigation of radio frequency interference (RFI) is the important step for the search, especially for the most sensitive Five-hundred-meter Aperture Spherical radio Telescope (FAST), which can detect more weak RFI. In this paper, we propose our procedure of RFI mitigation using several new methods, and use the procedure to perform a search for ETI signals from the data of FAST's first SETI commensal sky survey. We detect the persistent narrowband RFI by setting a threshold of the signals' sky separation, and detect the drifting RFI (and potentially other types of RFI) using the Hough transform method. We also use the clustering algorithms to remove more RFI and select candidates. The results of our procedure are compared to the earlier work on the same FAST data. We find that our methods, though relatively simpler in computation, remove more RFI, but preserve the simulated ETI signals except those severely affected by the RFI. We also report more interesting candidate signals, about a dozen of which are new candidates that are not previously reported. In addition, we find that the proposed Hough transform method, with suitable parameters, also has the potential to remove the broadband RFI. We conclude that our methods can effectively remove the vast majority of the RFI while preserving and finding the candidate signals that we are interested in.

We derive a new constraint on the expansion history of the Universe by applying the cosmic chronometers method, studying the age evolution of high-redshift galaxies with a full-spectral-fitting approach. We select a sample of 39 massive ($log(M/M_\odot)>10.8$) and passive ($log(sSFR/yr^{-1})<-11$) galaxies from the data release 4 of the VANDELS survey at $1<z<1.5$, combining different selection criteria to minimize the potential contamination by star-forming outliers. We perform full-spectral-fitting jointly on spectra and photometry of our sources with the code BAGPIPES, without any cosmological assumption on the age of the population. The derived physical properties of the selected galaxies are characteristic of a passive population, with short star formation timescales ($<\tau>=0.28\pm0.02$ Gyr), low dust extinction ($<A_{V,dust}>=0.43\pm0.02$ mag), and sub-solar metallicities ($<Z/Z_{\odot}>=0.44\pm0.01$). The ages show a decreasing trend with redshift compatible with a standard cosmological model, even if no cosmological constraint is assumed in the fit, and a clear mass-downsizing pattern. Testing the impact of the star formation history on the results, we find only a maximum 2\% fluctuation in age and metallicity. By fitting the median age-redshift relation with a flat $\Lambda$CDM model and assuming a Gaussian prior on $\Omega_{M,0}= 0.3\pm0.02$ from late-Universe probes, we obtain $H_0=67_{-15}^{+14}\:km\:s^{-1}\:Mpc^{-1}$. In the end, we derive a new estimate of the Hubble parameter with the cosmic chronometers method, $H(z=1.26)=135\pm65\:km\:s^{-1}\:Mpc^{-1}$ including statistical and systematic errors. While the error budget is currently dominated by the scarcity of the sample, this work proves the potential strength of the cosmic chronometers approach up to $z>1$, especially in view of incoming large spectroscopic surveys like Euclid. (abridged)

We present multi-wavelength characterization of 65 high mass X-ray binary (HMXB) candidates in M33. We use the Chandra ACIS survey of M33 (ChASeM33) catalog to select hard X-ray point sources that are spatially coincident with UV-bright point source optical counterparts in the Panchromatic Hubble Andromeda Treasury: Triangulum Extended Region (PHATTER) catalog, which covers the inner disk of M33 at near infrared, optical, and near ultraviolet wavelengths. We perform spectral energy distribution (SED) fitting on multi-band photometry for each point source optical counterpart to measure its physical properties including mass, temperature, luminosity, and radius. We find that the majority of the HMXB companion star candidates are likely B-type main sequence stars, suggesting that the HMXB population of M33 is dominated by Be-XRBs, as is seen in other Local Group galaxies. We use spatially-resolved recent star formation history (SFH) maps of M33 to measure the age distribution of the HMXB candidate sample and the HMXB production rate for M33. We find a bimodal distribution for the HMXB production rate over the last 80 Myr, with a peak at $\sim$10 Myr and $\sim$40 Myr, which match theoretical formation timescales for the most massive HMXBs and Be X-ray binaries (Be-XRBs), respectively. We measure an HMXB production rate of 107$-$136 HMXBs/(M$_{\odot}$ yr$^{-1}$) over the last 50 Myr and 150$-$199 HMXBs/(M$_{\odot}$ yr$^{-1}$) over the last 80 Myr. For sources with compact object classifications from overlapping NuSTAR observations, we find a preference for giant/supergiant companion stars in BH-HMXBs and main sequence companion stars in neutron star HMXBs (NS-HMXBs).

Asteroid diameters are traditionally difficult to estimate. When a direct measurement of the diameter cannot be made through either occultation or direct radar observation, the most common method is to approximate the diameter from infrared observations. Once the diameter is known, a comparison with visible light observations can be used to find the visible geometric albedo of the body. One of the largest datasets of asteroid albedos comes from the NEOWISE mission, which measured asteroid albedos both in the visible and infrared. We model these albedos as a function of proper elements available from the Asteroid Families Portal using an ensemble of neural networks. We find that both the visible and infrared geometric albedos are significantly correlated with asteroid position in the belt and occur in both asteroid families and in the background belt. We find that the ensemble's prediction reduces the average error in albedo by about 37% compared to a model that simply adopts an average albedo, with no regard for the dynamical state of the body. We then use this model to predict albedos for the half million main belt asteroids with proper elements available in the Asteroid Families Portal and provide the results in a catalog. Finally, we show that several presently categorized asteroid families exist within much larger groups of asteroids of similar albedos - this may suggest that further improvements in family identification can be made.

We report on the discovery of cool gas inflows towards three star-forming galaxies at $\left<z\right>\sim$ 2.30. Analysis of Keck Low-Resolution Imaging Spectrometer spectroscopy reveals redshifted low-ionisation interstellar (LIS) metal absorption lines with centroid velocities of 60 - 130 km $\rm{s}^{-1}$. These inflows represent some of the most robust detections of inflowing gas into isolated, star-forming galaxies at high redshift. Our analysis suggests that the inflows are due to recycling metal-enriched gas from previous ejections. Comparisons between the galaxies with inflows and a larger parent sample of 131 objects indicate that galaxies with detected inflows may have higher specific star-formation rates (sSFR) and star-formation-rate surface densities. However, when additional galaxies without robustly detected inflows based on centroid velocity but whose LIS absorption line profiles indicate large red-wing velocities are considered, galaxies with inflows do not show unique properties relative to those lacking inflows. Additionally, we calculate the covering fraction of cool inflowing gas as a function of red-wing inflow velocity, finding an enhancement in high sSFR binned galaxies, likely due to an increase in the amount of recycling gas. Together, these results suggest that the low detection rate of galaxies with cool inflows is primarily related to the viewing angle rather than the physical properties of the galaxies.

The gravitational wave (GW) antenna LISA will detect the signal from coalescing massive black hole binaries (MBHBs) of $\rm 10^4\,{-}\,10^7\, M_{\odot}$, providing clues on their formation and growth along cosmic history. Some of these events will be localized with a precision of several to less than a deg$^2$, enabling the possible identification of their host galaxy. This work explores the properties of the host galaxies of LISA MBHBs below $z\,{\lesssim}\,3$. We generate a simulated lightcone by using the semi-analytical model $\mathrm{\texttt{L-Galaxies}}$ applied on the merger trees of the high-resolution N-body cosmological simulation $\mathrm{\texttt{Millennium-II}}$. The model shows that LISA MBHBs are expected to be found in optically dim ($r\,{>}\,20$), star-forming ($\rm sSFR\,{>}\,10^{-10}\, \rm yr^{-1}$), gas-rich ($f_{\rm gas}\,{>}\,0.6$) and disc-dominated ($\rm B/T\,{<}\,0.7$) \textit{low-mass galaxies} of stellar masses $10^8\,{-}\,10^9 M_{\odot}$. However, these properties are indistinguishable from those of galaxies harboring single massive black holes with comparable mass, making difficult the selection of LISA hosts among the whole population of low-mass galaxies. Motivated by this, we explore the possibility of using merger signatures to select LISA hosts. We find that 40-80% of the galaxies housing LISA MBHBs display merger features related to the interaction which brought the secondary MBH to the galaxy. Despite this, around 60% of dwarf galaxies placed in the surroundings of the LISA hosts will show such kind of features as well, challenging the unequivocal detection of LISA hosts through the search of merger signatures. Consequently, the detection of an electromagnetic transient associated with the MBHB merger will be vital to pinpoint the star-forming dwarf galaxy where these binary systems evolve and coalesce.

Rapidly rotating magnetars have been associated with gamma-ray bursts (GRBs) and super-luminous supernovae (SLSNe). Using a suite of 2D magnetohydrodynamic simulations at fixed neutrino luminosity and a couple of evolutionary models with evolving neutrino luminosity and magnetar spin period, we show that magnetars are viable central engines for powering GRBs and SLSNe. We also present analytic estimates of the energy outflow rate from the proto-neutron star (PNS) as a function polar magnetic field strength $B_0$, PNS angular velocity $\Omega_{\star}$, PNS radius $R_{\star}$ and mass outflow rate $\dot{M}$. We show that rapidly rotating magnetars with spin periods $P_{\star}\lesssim 4$ ms and polar magnetic field strength $B_0\gtrsim 10^{15}$ G can release $10^{50}-5\times 10^{51}$ ergs of energy during the first $\sim2$ s of the cooling phase. This magnitude of energy release is sufficient to power long duration GRBs. We also show that magnetars with moderate field strengths of $B_0\lesssim 5\times 10^{14}$ G do not release a large fraction of their rotational kinetic energy during the cooling phase and hence, are not likely to power GRBs. Although we cannot simulate to times greater than $\sim 3-5$ s after a supernova, we can hypothesize that moderate field strength magnetars can brighten the supernova light curves by releasing their rotational kinetic energy via magnetic dipole radiation on timescales of days to weeks, since these do not expend most of their rotational kinetic energy during the early cooling phase.

We have analyzed the thermal and turbulent properties of the Low-Latitude Intermediate-Velocity Arch 1 (LLIV1). This was accomplished using archival H I emission and absorption data from two 21,cm line surveys: GHIGLS at $9.^\prime$4 resolution and DHIGLS at $1^\prime$ resolution. The spectral decomposition code $\tt{ROHSA}$ was used to model the column density of different thermal phases and also to analyze an absorption measurement against the radio source 4C~+66.09. From the latter we found spin temperature $T_{\mathrm{s}} \sim 75$K, cold gas mass fraction $f\sim0.5$, and turbulent sonic Mach number $M_t\sim3.4$. Similar to the absorption line modeling against 4C~+66.09, our best emission line decomposition model has no unstable gas across the whole field of view, suggesting that the thermal condensation and phase transition are not on-going but rather have reached an equilibrium state. The cold phase of LLIV1 appears as a collection of elongated filaments that forms a closed structure within the field decomposed. These substructures follow the orientation of the overall large scale cloud, along the diagonal of the GHIGLS field from north-west to south-east (in Galactic coordinates). The angular power spectrum of the cold phase is slightly shallower than that of the warm phase, quantifying that the cold phases have relatively more structure on small scales. Our spatially resolved map of the cold gas mass fraction in LLIV1 from DHIGLS reveals significant variations spanning the possible range of $f$, with mean and standard deviation 0.33 and 0.19, respectively.

The aim of this study is to explore the magnetic and flow properties of fully convective M dwarfs as a function of rotation period Prot and magnetic Reynolds ReM and Prandlt numbers PrM. We performed three-dimensional simulations of fully convective stars using a star-in-a-box setup. This setup allows global dynamo simulations in a sphere embedded in a Cartesian cube. The equations of non-ideal magnetohydrodynamics were solved with the Pencil Code. We used the stellar parameters of an M5 dwarf with 0.21M_odot at three rotation rates corresponding to rotation periods (Prot): 43, 61 and 90 days, and varied the magnetic Prandtl number in the range from 0.1 to 10. We found systematic differences in the behaviour of the large-scale magnetic field as functions of rotation and PrM. For the simulations with Prot = 43 days and PrM <= 2, we found cyclic large-scale magnetic fields. For PrM > 2 the cycles vanish and field shows irregular reversals. In simulations with Prot = 61 days for PrM <= 2 the cycles are less clear and the reversal are less periodic. In the higher-PrM cases, the axisymmetric mean field shows irregular variations. For the slowest rotation case with Prot = 90 days, the field has an important dipolar component for PrM <= 5. For the highest PrM the large-scale magnetic field is predominantly irregular at mid-latitudes, with quasi-stationary fields near the poles. For the simulations with cycles, the cycle period length slightly increases with increasing ReM.

We present the construction of stationary boson-fermion spherically symmetric configurations governed by Newtonian gravity. Bosons are described in the Gross-Pitaevskii regime and fermions are assumed to obey Euler equations for an inviscid fluid with polytropic equation of state. The two components are coupled through the gravitational potential. The families of solutions are parametrized by the central value of the wave function describing the bosons and the central denisty of the fluid. We explore the stability of the solutions using numerical evolutions that solve the time dependent Schr\"odinger-Euler-Poisson system, using the truncation error of the numerical methods as the perturbation. We find that all configurations are stable as long as the polytropic equation of state (EoS) is enforced during the evolution. When the configurations are evolved using the ideal gas EoS they all are unstable that decay into a sort of twin solutions that approach a nearly stationary configuration. We expect these solutions and their evolution serve to test numerical codes that are currently being used in the study of Fuzzy Dark Matter plus baryons.

Warm Jupiters are close-in giant planets with relatively large planet-star separations (i.e., $10< a/R_\star <100$). Given their weak tidal interactions with their host stars, measurements of stellar obliquity may be used to probe the initial obliquity distribution and dynamical history for close-in gas giants. Using spectroscopic observations, we confirm the planetary nature of TOI-1859b and determine the stellar obliquity of TOI-1859 to be $\lambda = 38.9^{+2.8}_{-2.7}\deg$ relative to its planetary companion using the Rossiter-McLaughlin effect. TOI-1859b is a 64-day warm Jupiter orbiting around a late-F dwarf and has an orbital eccentricity of $0.57^{+0.12}_{-0.16}$, inferred purely from transit light curves. The eccentric and misaligned orbit of TOI-1859b is likely an outcome of dynamical interactions, such as planet-planet scattering and planet-disk resonance crossing.

The IceCube Neutrino Observatory is a Cherenkov detector located at the South Pole. Its main component consists of an in-ice array of optical modules instrumenting one cubic kilometer of deep Glacial ice. The DeepCore sub-detector is a denser in-fill array with a lower energy threshold, allowing us to study atmospheric neutrinos oscillations with energy below 100 GeV arriving through the Earth. We present preliminary results of an atmospheric muon neutrino disappearance analysis using data from 2012 to 2021 and employing convolutional neural networks (CNNs) for precise and fast event reconstructions.

COMPLETE is a flagship mission concept combining broadband spectroscopic imaging and comprehensive magnetography from multiple viewpoints around the Sun to enable tomographic reconstruction of 3D coronal magnetic fields and associated dynamic plasma properties, which provide direct diagnostics of energy release. COMPLETE re-imagines the paradigm for solar remote-sensing observations through purposefully co-optimized detectors distributed on multiple spacecraft that operate as a single observatory, linked by a comprehensive data/model assimilation strategy to unify individual observations into a single physical framework. We describe COMPLETE's science goals, instruments, and mission implementation. With targeted investment by NASA, COMPLETE is feasible for launch in 2032 to observe around the maximum of Solar Cycle 26.

Heliophysics image data largely relies on a forty-year-old ecosystem built on the venerable Flexible Image Transport System (FITS) data standard. While many in situ measurements use newer standards, they are difficult to integrate with multiple data streams required to develop global understanding. Additionally, most data users still engage with data in much the same way as they did decades ago. However, contemporary missions and models require much more complex support for 3D multi-parameter data, robust data assimilation strategies, and integration of multiple individual data streams required to derive complete physical characterizations of the Sun and Heliospheric plasma environment. In this white paper we highlight some of the 21$^\mathsf{st}$ century challenges for data frameworks in heliophysics, consider an illustrative case study, and make recommendations for important steps the field can take to modernize its data products and data usage models. Our specific recommendations include: (1) Investing in data assimilation capability to drive advanced data-constrained models, (2) Investing in new strategies for integrating data across multiple instruments to realize measurements that cannot be produced from single observations, (3) Rethinking old data use paradigms to improve user access, develop deep understanding, and decrease barrier to entry for new datasets, and (4) Investing in research on data formats better suited for multi-dimensional data and cloud-based computing.

In this paper, written as a general historical and technical introduction to the various review papers collected in the special issue ``Solar and Stellar Dynamo: A New Era'', we review the evolution and current state of dynamo theory and modelling, with emphasis on the solar dynamo. Starting with a historical survey, we then focus on a set of ``tension points'' that are still left unresolved despite the remarkable progress of the past century. In our discussion of these tension points we touch upon the physical well-posedness of mean-field electrodynamics; constraints imposed by magnetic helicity conservation; the troublesome role of differential rotation; meridional flows and flux transpost dynamos; competing inductive mechanisms and Babcock-Leighton dynamos; the ambiguous precursor properties of the solar dipole; cycle amplitude regulation and fluctuation through nonlinear backreaction and stochastic forcing, including Grand Minima; and the promises and puzzles offered by global magnethydrodynamical numerical simulations of convection and dynamo action. We close by considering the potential bridges to be constructed between solar dynamo theory and modelling, and observations of magnetic activity in late-type stars.

Faraday rotation measures (RMs) have been used for many studies of cosmic magnetism, and in most cases having more RMs is beneficial for those studies. This has lead to development of RM surveys that have produced large catalogs, as well as meta-catalogs collecting RMs from many different publications. However, it has been difficult to take full advantage of all these RMs as the individual catalogs have been published in many different places, and in many different formats. In addition, the polarization spectra used to determine these RMs are rarely published, limiting the ability to re-analyze data as new methods or additional observations become available. We propose a standard convention for RM catalogs, RMTable2023, and a standard for source-integrated polarized spectra of radio sources, PolSpectra2023. These standards are intended to maximize the value and utility of these data for researchers and to make them easier to access. To demonstrate the use of the RMTable2023 standard, we have produced a consolidated catalog of 55 819 RMs collected from 42 published catalogs.

We point out that the recent result of primordial helium-4 ($^4$He) abundance measurement by EMPRESS, which has reported a smaller $^4$He abundance than other measurements, can be well fitted by assuming a time-variation of the fine structure constant $\alpha$ which is slightly smaller than the present value during big bang nucleosynthesis (BBN). We find that the EMPRESS result in combination with deuterium abundance measurement indicates $-2.6\% <\Delta\alpha/\alpha <-1.4 \%$ (68\% C.L.) where $\Delta \alpha$ is the difference between the values of $\alpha$ at the BBN and present epochs, while $-1.2\% <\Delta\alpha/\alpha <0.4 \%$ (68\% C.L.) is obtained from other previous $^4$He abundance data. We also investigate its effects in the framework where the effective number of neutrino species and the lepton asymmetry, which are other typical interpretations of the EMPRESS result, are allowed to vary. Once a smaller $\alpha$ is adopted, the EMPRESS result can be explained without assuming any non-standard values for the effective number of neutrino species and lepton asymmetry.

The relative photometric calibration errors in the DESI Legacy Imaging Surveys (LS), which are used for DESI target selection, can leave imprints on the DESI target densities and bias the resulting cosmological measurements. We characterize the LS calibration systematics by comparing the LS stellar photometry with Gaia DR3 synthetic photometry. We find the stellar photometry of LS DR9 and Gaia has an \textsc{rms} difference of 4.7, 3.7, 4.4 mmag in DECam $grz$ bands, respectively, when averaged over an angular scale of 27 arcmin. There are distinct spatial patterns in the photometric offset resembling the Gaia scan patterns (most notably in the synthesized $g$-band) which indicate systematics in the Gaia spectrophotometry, as well as honeycomb patterns due to LS calibration systematics. We also find large and smoothly varying photometric offsets at $\mathrm{Dec}<-29.25^{\circ}$ in LS DR9 which are fixed in DR10.

We present the observation of a charge-sign dependent solar modulation of galactic cosmic rays (GCRs) with the CALorimetric Electron Telescope onboard the International Space Station over 6 yr, corresponding to the positive polarity of the solar magnetic field. The observed variation of proton count rate is consistent with the neutron monitor count rate, validating our methods for determining the proton count rate. It is observed by the CALorimetric Electron Telescope that both GCR electron and proton count rates at the same average rigidity vary in anticorrelation with the tilt angle of the heliospheric current sheet, while the amplitude of the variation is significantly larger in the electron count rate than in the proton count rate. We show that this observed charge-sign dependence is reproduced by a numerical ``drift model'' of the GCR transport in the heliosphere. This is a clear signature of the drift effect on the long-term solar modulation observed with a single detector.

Several studies have detected Ly$\alpha$ from bright ($M_{UV}\lesssim-21.5$) galaxies during the early stages of reionization despite the significantly neutral IGM. To explain these detections, it has been suggested that z>7 Ly$\alpha$ emitters (LAEs) inhabit large, physical Mpc (pMpc)-scale ionized regions. However, systematic searches for the overdensities of faint galaxies expected to be powering these ionized bubbles around LAEs have been challenging. Here, we use CEERS NIRCam imaging to investigate the possibility of galaxy overdensities associated with two very UV-bright LAEs at z=8.7 in the EGS field. We design a color selection to identify objects at z=8.4-9.1, selecting 28 candidates (including the one LAE in the footprint, EGSY8p7). We model the SEDs of these objects and infer that all are moderately faint ($-21.2\lesssim M_{UV}\lesssim-19.1$) with stellar masses of $M_* \approx 10^{7.5 - 8.8}$ $M_\odot$. All are efficient ionizing agents ($\xi_{ion}^*\sim10^{25.5-26.0}$ Hz erg$^{-1}$) and are generally morphologically simple with only one compact ($r_e\lesssim140$ to $\sim650$ pc) star-forming component. Of the 27 objects besides EGSY8p7, 13 lie within 5' of EGSY8p7, leading to a $4\times$ overdensity in projection at separations <5' (1.4 pMpc in projection at z=8.7). Separations of 10'-15' (2.7-4.1 projected pMpc) are consistent with an average field. The spatial distribution of our sample may qualitatively suggest a large ($R\geq2$ pMpc) ionized bubble encompassing both LAEs in the field, which is theoretically unexpected but may be possible for a galaxy population four times more numerous than the average to create given moderate escape fractions ($f_{esc}\gtrsim0.15$) over long times ($\gtrsim200$ Myr). Upcoming spectroscopic follow up will enable characterization of the size of any ionized bubble that may exist and the properties of the galaxies powering such a bubble.

The combined fit of the measured energy spectrum and shower maximum depth distributions of ultra-high-energy cosmic rays is known to constrain the parameters of astrophysical models with homogeneous source distributions. Studies of the distribution of the cosmic-ray arrival directions show a better agreement with models in which a fraction of the flux is non-isotropic and associated with the nearby radio galaxy Centaurus A or with catalogs such as that of starburst galaxies. Here, we present a novel combination of both analyses by a simultaneous fit of arrival directions, energy spectrum, and composition data measured at the Pierre Auger Observatory. We find that a model containing a flux contribution from the starburst galaxy catalog of around 20% at 40 EeV with a magnetic field blurring of around $20^\circ$ for a rigidity of 10 EV provides a fair simultaneous description of all three observables. The starburst galaxy model is favored with a significance of $4.5\sigma$ (considering experimental systematic effects) compared to a reference model with only homogeneously distributed background sources. By investigating a scenario with Centaurus A as a single source in combination with the homogeneous background, we confirm that this region of the sky provides the dominant contribution to the observed anisotropy signal. Models containing a catalog of jetted active galactic nuclei whose flux scales with the $\gamma$-ray emission are, however, disfavored as they cannot adequately describe the measured arrival directions.

A deep survey of the Large Magellanic Cloud at ~0.1-100TeV photon energies with the Cherenkov Telescope Array is planned. We assess the detection prospects based on a model for the emission of the galaxy, comprising the four known TeV emitters, mock populations of sources, and interstellar emission on galactic scales. We also assess the detectability of 30 Doradus and SN 1987A, and the constraints that can be derived on the nature of dark matter. The survey will allow for fine spectral studies of N157B, N132D, LMC P3, and 30 Doradus C, and half a dozen other sources should be revealed, mainly pulsar-powered objects. The remnant from SN 1987A could be detected if it produces cosmic-ray nuclei with a flat power-law spectrum at high energies, or with a steeper index 2.3-2.4 pending a flux increase by a factor >3-4 over ~2015-2035. Large-scale interstellar emission remains mostly out of reach of the survey if its >10GeV spectrum has a soft photon index ~2.7, but degree-scale 0.1-10TeV pion-decay emission could be detected if the cosmic-ray spectrum hardens above >100GeV. The 30 Doradus star-forming region is detectable if acceleration efficiency is on the order of 1-10% of the mechanical luminosity and diffusion is suppressed by two orders of magnitude within <100pc. Finally, the survey could probe the canonical velocity-averaged cross section for self-annihilation of weakly interacting massive particles for cuspy Navarro-Frenk-White profiles.

We have performed a spectral-timing analysis on the black hole X-ray binary MAXI J1535--571 during its 2017 outburst, with the aim of exploring the evolution of the inner accretion flow geometry. X-ray reverberation lags are observed in the hard-intermediate state (HIMS) and soft-intermediate state (SIMS) of the outburst. During the HIMS, the characteristic frequency of the reverberation lags $\nu_0$ (the frequency at which the soft lag turns to zero in the lag-frequency spectra) increases when the spectrum softens. This reflects a reduction of the spatial distance between the corona and accretion disc, when assuming the measured time lags are associated with the light travel time. We also find a strong correlation between $\nu_0$ and type-C Quasi Periodic Oscillation (QPO) centroid frequency $\nu_{QPO}$, which can be well explained by the Lense-Thirring (L-T) precession model under a truncated disk geometry. Despite the degeneracy in the spectral modellings, our results suggest that the accretion disc is largely truncated in the low hard state (LHS), and moves inward as the spectrum softens. Combine the spectral modelling results with the $\nu_0$ - $\nu_{QPO}$ evolution, we are inclined to believe that this source probably have a truncated disk geometry in the hard state.

In the overshoot mixing model with an exponentially decreasing diffusion coefficient, the initial value of the diffusion coefficient plays a crucial role. According to the turbulent convective mixing model, the characteristic length of convection in the convection zone differs from that in the overshoot region, resulting in a rapid decrease of the diffusion coefficient near the convective boundary. To investigate this quick decrease, we conducted an asteroseismic study on the intermediate-mass SPB star KIC 10526294. We generated stellar models with varied input parameters, including the overshoot parameters, and compared the resulting stellar oscillation periods with observations. To mitigate the potential issue arising from large steps in the stellar parameters and stellar age, we employed a comprehensive interpolation scheme for the stellar oscillatory frequencies, considering all stellar parameters and stellar age. Our analysis revealed that the quick decreasing of the diffusion coefficient has discernible effects on the stellar oscillations and a quick decrease with 4 magnitude orders shows the best oscillatory frequencies compared with the observations. This provides weak evidence in support of the prediction made by the turbulent convective mixing model. Furthermore, we examined the residuals of the oscillation periods and discovered a potential association between abundance anomalies in the buoyancy frequency profile and the oscillation-like patterns observed in the residuals.

Timing observations are crucial for determining the physical properties of newly discovered pulsars. Using the L-band 19-beam receiver of the Five-hundred-meter Aperture Spherical radio Telescope (FAST), the FAST Galactic Plane Pulsar Snapshot (GPPS) survey has discovered many faint and weak pulsars, which are hardly detected using other radio telescopes in limited observation time. To obtain accurate position, spin parameters, dispersion measure, and to calculate derived parameters such as characteristic age and surface magnetic fields, we collect available FAST pulsar data obtained either through targeted following-up observations or coincidental survey observations with one of the 19 beams of the L-band 19-beam receiver. From these data we get the time of arrival (TOA) measurements for 30 newly discovered pulsars as well as 13 known pulsars. We demonstrate that TOA measurements from any beams of the L-band 19-beam receiver acquired through any FAST observation mode (e.g., the tracking mode or the snapshot mode) can be combined together for getting timing solutions. We update the ephemeris of 13 previously known pulsars and obtained the first phase-coherent timing results for 30 isolated pulsars discovered in the FAST GPPS survey. Notably, PSR J1904+0853 is an isolated millisecond pulsar, PSR J1906+0757 is a disrupted recycled pulsar, and PSR J1856+0211 is a long-period pulsar that can constrain pulsar death lines. Based on these timing solutions, all available FAST data can be added together to get the best pulse profiles.

We present a new algorithm for the identification and physical characterization of current sheets and reconnection sites in 2D and 3D large scale relativisticmagnetohydrodynamic numerical simulations. This has been implemented in the PLUTO code and tested in the cases of a single current sheet, a 2D jet and a 3D unstable plasma column. Its main features are: a) a computational cost which allows its use in large scale simulations; b) the capability to deal with complex 2D and 3D structures of the reconnection sites. In the performed simulations, we identify the computational cells that are part of a current sheet by a measure of the gradient of the magnetic field along different directions. Lagrangian particles, which follow the fluid, are used to sample plasma parameters before entering the reconnection sites that form during the evolution of the different configurations considered. Specifically, we track the distributions of the magnetization parameter $\sigma$ and the thermal to magnetic pressure ratio $\beta$ that - according to particle-in-cell simulation results - control the properties of particle acceleration in magnetic reconnection regions. Despite the initial conditions of the simulations were not chosen "ad hoc", the 3D simulation returns results suitable for efficient particle acceleration and realistic non-thermal particle distributions.

We present the results from deep 21 cm H I mapping of two nearby blue compact dwarf galaxies (BCDGs), W1016+3754 and W2326+0608, using the Giant Metrewave Radio Telescope (GMRT). These BCDGs are bright in mid-infrared (MIR) data and undergoing active star formation. With the GMRT observations, we investigate the role of cold neutral gas as the fuel resource of the current intensive star formation activity. Star formation in these galaxies is likely to be due to the infall of H I gas triggered by gravitational perturbation from nearby galaxies. The BCDG W2326+0608 and nearby galaxy SDSS J232603.86+060835.8 share a common H I envelope. We find star formation takes place in the high H I column density gas ($\gtrsim 10^{21}$\,cm$^{-2}$) regions for both BCDGs. The recent starburst and infall of metal-free gas have kept the metallicity low for the BCDG W1016+3754. The metallicity for W2326+0608 is higher, possibly due to tidal interaction with the nearby galaxy SDSS J232603.86+060835.8.

Many of the short-lived radioactive nuclei that were present in the early Solar System can be produced in massive stars. In the first paper in this series (Brinkman et al. 2019), we focused on the production of $^{26}$Al in massive binaries. In our second paper (Brinkman et al. 2021), we considered rotating single stars, two more short-lived radioactive nuclei, $^{36}$Cl and $^{41}$Ca, and the comparison to the early Solar System data. In this work, we update our previous conclusions by further considering the impact of binary interactions. We used the MESA stellar evolution code with an extended nuclear network to compute massive (10-80 M$ _{\odot} $), binary stars at various initial periods and solar metallicity (Z=0.014), up to the onset of core collapse. The early Solar System abundances of $^{26}$Al and $^{41}$Ca can be matched self-consistently by models with initial masses $\geq$25 M$_{\odot}$, while models with initial primary masses $\geq$35 M$_{\odot}$ can also match $^{36}$Cl. Almost none of the models provide positive net yields for $^{19}$F, while for $^{22}$Ne the net yields are positive from 30 M$_{\odot}$ and higher. This leads to an increase by a factor of approximately 4 in the amount of $^{22}$Ne produced by a stellar population of binary stars, relative to single stars. Also, besides the impact on the stellar yields, our 10 M$_{\odot}$ primary star undergoing Case A mass-transfer ends its life as a white dwarf instead of as a core-collapse supernova. This demonstrates that binary interactions can also strongly impact the evolution of stars close to the supernova boundary.

We review the state of the art of three dimensional numerical simulations of solar and stellar dynamos. We summarize fundamental constraints of numerical modelling and the techniques to alleviate these restrictions. Brief summary of the relevant observations that the simulations seek to capture is given. We survey the current progress of simulations of solar convection and the resulting large-scale dynamo. We continue to studies that model the Sun at different ages and to studies of stars of different masses and evolutionary stages. Both simulations and observations indicate that rotation, measured by the Rossby number which is the ratio of rotation period and convective turnover time, is a key ingredient in setting the overall level and characteristics of magnetic activity. Finally, efforts to understand global 3D simulations in terms of mean-field dynamo theory are discussed.

It is known that the single-field inflation with a transient ultra-slow-roll phase can produce a large curvature perturbation at small scales for the formation of primordial black holes. In our previous work, we have considered quantum loop corrections to the curvature perturbation and found that the growth of these small-scale modes would affect the curvature perturbation at large scales probed by cosmic microwave background observation. In this work, we will further derive the constraints on the growing modes in the transition between the slow-roll and the ultra-slow-roll phases under the effect of the loop corrections. Our results would help clarify the recent controversy on whether or not the primordial-black-hole formation from the single-field inflation is ruled out at one-loop level.

Context: Chondrules originate from the reprocessing of dust grains. They are key building blocks of telluric planets, yet their formation, which must happen in strongly localized regions of high temperature, remains poorly understood. Aims: We examine the dust spatial distribution near regions of strong local heating produced by current sheets, as a step toward exploring a potential path for chondrule formation. We further aim to investigate current sheet formation under various conditions in protoplanetary disks in the presence of ambipolar diffusion and Ohmic resistivity and the effect of current sheet morphology on dust dynamics in their vicinity. Methods: We use the RAMSES code including modules for non-ideal magnetohydrodynamics and solution of the dynamics of multiple sizes of dust grains to compute unstratified shearing box simulations of current sheet formation. We investigate, through seven models the effect of the ambipolar diffusion and Ohmic resistivity strength, the initial density, and magnetic field, as well as the resolution and box size. Results: We find that current sheets form in all our models, with typical widths of 0.001-0.01 AU and that strong dust fraction variations occur for millimeter-sized grains. These variations are typically of an order of magnitude and up to two orders of magnitude for the most favorable cases. We also show that the box size and resolution has a strong impact on the current sheet distribution and intensity. Conclusions: The formation of current sheets that can intensely heat their surroundings near strong dynamical dust fraction variations could have important implications for chondrule formation, as it appears likely to happen in regions of large dust fraction.

New large observational surveys such as Gaia are leading us into an era of data abundance, offering unprecedented opportunities to discover new physical laws through the power of machine learning. Here we present an end-to-end strategy for recovering a free-form analytical potential from a mere snapshot of stellar positions and velocities. First we show how auto-differentiation can be used to capture an agnostic map of the gravitational potential and its underlying dark matter distribution in the form of a neural network. However, in the context of physics, neural networks are both a plague and a blessing as they are extremely flexible for modeling physical systems but largely consist in non-interpretable black boxes. Therefore, in addition, we show how a complementary symbolic regression approach can be used to open up this neural network into a physically meaningful expression. We demonstrate our strategy by recovering the potential of a toy isochrone system.

Investigations of the phase diagram of quantum chromodynamics (QCD) have revealed that exotic new phases, the so called {\it color superconducting phases}, may arise at very high baryon densities. It is speculated that these exotic phases may arise in the cores of neutron stars. Focus on neutrons stars has tremendously intensified in recent years with the direct detection of gravitational waves (GW) by LIGO/Virgo from BNS merger events which has allowed the possibility of directly probing the properties of the interior of a neutron star. A remarkable phenomenon manifested by rapidly rotating neutron stars is in their {\it avatar} as {\it Pulsars}. The accuracy of pulsar timing can reach the level of one part in 10$^{15}$, comparable to that of atomic clocks. This suggests that even a tiny deformation of the pulsar can leave its imprints on the pulses by inducing tiny perturbations in the entire moment of inertia (MI) tensor affecting the pulse timings, as well as the pulse profile (from wobbling induced by off-diagonal MI components). This may allow a new probe of various phase transitions occurring inside a pulsar core through induced density fluctuations affecting the MI tensor. Such perturbations also naturally induce a rapidly changing quadrupole moment of the star, thereby providing a new source of gravitational wave emission. Another remarkable possibility arises when we consider the effect of an external GW on neutron star. With the possibility of detecting any minute changes in its configuration through pulse observations, the neutron star has the potential of performing as a Weber detector of gravitational wave. This brief review will focus on these specific aspects of a pulsar. Specifically, the focus will be on the type of physics which can be probed by utilizing the effect of changes in the MI tensor of the pulsar on pulse properties.

In this study we investigate the growth index $\gamma_L$, which characterizes the growth of linear matter perturbations, while analysing different cosmological datasets. We compare the approaches implemented by two different patches of the cosmological solver CAMB: MGCAMB and CAMB_GammaPrime_Growth. In our analysis we uncover a deviation of the growth index from its expected $\Lambda$CDM value of $0.55$ when utilizing the Planck dataset, both in the MGCAMB case and in the CAMB_GammaPrime_Growth case, but in opposite directions. This deviation is accompanied by a change in the direction of correlations with derived cosmological parameters. However, the incorporation of CMB lensing data helps reconcile $\gamma_L$ with its $\Lambda$CDM value in both cases. Conversely, the alternative ground-based telescopes ACT and SPT consistently yield growth index values in agreement with $\gamma_L=0.55$. We conclude that the presence of the A$_{\mathrm{lens}}$ problem in the Planck dataset contributes to the observed deviations, underscoring the importance of additional datasets in resolving these discrepancies.

The newly operational JWST offers the potential to study the atmospheres of distant worlds with precision that has not been achieved before. One of the first exoplanets observed by JWST in the summer of 2022 was WASP-96 b, a hot-Saturn orbiting a G8 star. As part of the Early Release Observations program, one transit of WASP-96 b was observed with NIRISS/SOSS to capture its transmission spectrum from 0.6-2.85 microns. In this work, we utilise four retrieval frameworks to report precise and robust measurements of WASP-96 b's atmospheric composition. We constrain the logarithmic volume mixing ratios of multiple chemical species in its atmosphere, including: H$_2$O = $-3.59 ^{+ 0.35 }_{- 0.35 }$, CO$_2$ = $-4.38 ^{+ 0.47 }_{- 0.57 }$ and K = $-8.04 ^{+ 1.22 }_{- 1.71 }$. Notably, our results offer a first abundance constraint on potassium in WASP-96 b's atmosphere, and important inferences on carbon-bearing species such as CO$_2$ and CO. Our short wavelength NIRISS/SOSS data are best explained by the presence of an enhanced Rayleigh scattering slope, despite previous inferences of a clear atmosphere - although we find no evidence for a grey cloud deck. Finally, we explore the data resolution required to appropriately interpret observations using NIRISS/SOSS. We find that our inferences are robust against different binning schemes. That is, from low $R = 125$ to the native resolution of the instrument, the bulk atmospheric properties of the planet are consistent. Our systematic analysis of these exquisite observations demonstrates the power of NIRISS/SOSS to detect and constrain multiple molecular and atomic species in the atmospheres of hot giant planets.

The standard $\Lambda$CDM model of cosmology is largely successful in describing many observations, including precise measurements of the Cosmic Microwave Background (CMB) radiation. However, some intriguing anomalies remain currently unexplained within this theoretical framework. Such discrepancies can be broadly categorized into two groups: those involving CMB-independent probes and those within different CMB experiments. Examples of the former category include the $H_0$-tension between the value of the present-day expansion rate of the Universe inferred by CMB observations and local distance ladder measurements. The latter category involves anomalies between the values of cosmological parameters obtained by different CMB experiments, as well as their consistency with the predictions of $\Lambda$CDM. In this chapter, we primarily focus on this second category and study the agreement among the most recent CMB measurements to acquire a better understanding of the limitations and uncertainties underlying both the current data and the cosmological model. Finally, we discuss the implications for the $H_0$-tension.

Context. The effects of planetesimal fragmentation on planet formation has been studied by various models on single embryos therefore neglecting concurrent effects mostly in the outer disk. They show that planetesimal fragmentation can both hinder or aid planet formation due to the introduction of competing effects, namely speeding up accretion and depleting the feeding zone of forming planets. Aims. We investigate the influence of the collisional fragmentation of planetesimals on the planet formation process using a population synthesis approach. We aim to investigate its effects for a large set of initial conditions and also explore the consequences on the formation of multiple embryos in the same disk. Methods. We run global planet formation simulations including fragmentation, drift and an improved ice line description. To do this we use a fragmentation model in our code. The initial conditions for the simulations that are informed by observations are varied to generate synthetic exoplanet populations. Results. Our synthetic populations show that depending on the typical size of solids generated in collisions, fragmentation in tandem with the radial drift can either enhance or hinder planet formation. For larger fragments we see increased accretion throughout the populations especially beyond the ice line. However, the shorter drift timescale of smaller fragments, due to their stronger coupling to the gas, can hinder the formation process. Furthermore, beyond the ice line fragmentation promotes late growth when the damping by gas drag fades Conclusions. Fragmentation significantly affects the planet formation process in various ways for all types of planet and warrants further investigation.

Gyrochronology allows the derivation of ages for cool main sequence stars from their observed rotation periods and masses, or a suitable proxy of the latter. It is increasingly well explored for FGK stars, but requires further measurements for older ages and K-M-type stars. Recent work has shown that the behavior of stellar spindown differs significantly from prior expectations for late-type stars. We study the 4Gyr-old benchmark open cluster M67 to explore this behavior further. We combined a Gaia DR3 sample with the Kepler K2 superstamp of Campaign 5 around M67 and created new light curves from aperture photometry. The light curves are subjected to an extensive correction process to remove instrumental systematics and trending, followed by period analysis to measure stellar rotation. We identify periodic signals in 136 light curves, 47 of which are from the rotation of effectively single main-sequence stars that span from early-G to mid-M type. These results connect well to prior work on M67 and extend it to much later spectral types. We find that the rotation periods of single stars of age 4Gyr define a tight relationship with color, ranging from spectral types F through M. The corresponding surface of rotation period against age and mass is therefore well-defined to an older age than was previously known. However, the deviations from prior expectations of the stellar spindown behavior are even more pronounced at 4Gyr. The binary cluster members do not follow the single star relationship. The majority are widely scattered below the single star sequence. Consequently, they do not seem to be suitable for gyrochronology at present.

The future is now - after its long-awaited launch in December 2021, JWST began science operations in July 2022 and is already revolutionizing exoplanet astronomy. The Early Release Observations (ERO) program was designed to provide the first images and spectra from JWST, covering a multitude of science cases and using multiple modes of each on-board instrument. Here, we present transmission spectroscopy observations of the hot-Saturn WASP-96b with the Single Object Slitless Spectroscopy (SOSS) mode of the Near Infrared Imager and Slitless Spectrograph, observed as part of the ERO program. As the SOSS mode presents some unique data reduction challenges, we provide an in-depth walk-through of the major steps necessary for the reduction of SOSS data: including background subtraction, correction of 1/f noise, and treatment of the trace order overlap. We furthermore offer potential routes to correct for field star contamination, which can occur due to the SOSS mode's slitless nature. By comparing our extracted transmission spectrum with grids of atmosphere models, we find an atmosphere metallicity between 1x and 5x solar, and a solar carbon-to-oxygen ratio. Moreover, our models indicate that no grey cloud deck is required to fit WASP-96b's transmission spectrum, but find evidence for a slope shortward of 0.9$\mu$m, which could either be caused by enhanced Rayleigh scattering or the red wing of a pressure-broadened Na feature. Our work demonstrates the unique capabilities of the SOSS mode for exoplanet transmission spectroscopy and presents a step-by-step reduction guide for this new and exciting instrument.

We present the first catalog of very-high energy and ultra-high energy $\gamma$-ray sources detected by the Large High Altitude Air Shower Observatory (LHAASO), using 508 days of data collected by the Water Cherenkov Detector Array (WCDA) from March 2021 to September 2022 and 933 days of data recorded by the Kilometer Squared Array (KM2A) from January 2020 to September 2022. This catalog represents the most sensitive $E > 1$ TeV gamma-ray survey of the sky covering declination from $-$20$^{\circ}$ to 80$^{\circ}$. In total, the catalog contains 90 sources with extended size smaller than $2^\circ$ and with significance of detection at $> 5\sigma$. For each source, we provide its position, extension and spectral characteristics. Furthermore, based on our source association criteria, 32 new TeV sources are proposed in this study. Additionally, 43 sources are detected with ultra-high energy ($E > 100$ TeV) emission at $> 4\sigma$ significance level.

We present second observations by K2 of OJ~287 and 7 other $\gamma$-ray AGNs obtained in 2017-2018, second and third observations of the lobe-dominated, steep spectrum quasar 3C~207, and observations of 9 additional blazars not previously observed with K2. The AGN were observed simultaneously with K2 and the Fermi Large Area Telescope for 51-81 days. Our full sample, observed in 2014-2018, contained 16 BL Lac objects (BL Lacs), 9 Flat Spectrum Radio Quasars (FSRQs), and 4 other $\gamma$-ray AGNs. Twelve BL Lacs and 7 FSRQs exhibited fast, jagged light curves while 4 BL Lacs and 2 FSRQs had slow, smooth light curves. Some objects changed their red-noise character significantly between repeated K2 observations. The optical characteristics of OJ~287 derived from the short-cadence K2 light curves changed between observations made before and after the predicted passage of the suspected secondary supermassive black hole through the accretion disk of the primary supermassive black hole. The average slopes of the periodogram power spectral densities of the BL Lacs' and FSRQs' light curves differed significantly, by $\approx 12$\%, with the BL Lac slopes being steeper, and a KS test with a $p$-value of 0.039 indicates that these samples probably come from different populations; however, this result is not as strongly supported by PSRESP analyses. Differences in the origin of the jets from the ergosphere or accretion disk in these two classes could produce such a disparity, as could different sizes or locations of emission regions within the jets.

Kernel phase interferometry (KPI) is a data processing technique that allows for the detection of asymmetries (such as companions or disks) in high-Strehl images, close to and within the classical diffraction limit. We show that KPI can successfully be applied to hyperspectral image cubes generated from integral field spectrographs (IFSs). We demonstrate this technique of spectrally-dispersed kernel phase by recovering a known binary with the SCExAO/CHARIS IFS in high-resolution K-band mode. We also explore a spectral differential imaging (SDI) calibration strategy that takes advantage of the information available in images from multiple wavelength bins. Such calibrations have the potential to mitigate high-order, residual systematic kernel phase errors, which currently limit the achievable contrast of KPI. The SDI calibration presented here is applicable to searches for line emission or sharp absorption features, and is a promising avenue toward achieving photon-noise-limited kernel phase observations. The high angular resolution and spectral coverage provided by dispersed kernel phase offers novel opportunities for science observations which would have been challenging to achieve otherwise.

Very recently, HAWC observatory discovered the high-energy gamma ray emission from the solar disk during the quiescent stage of the sun, extending the Fermi-LAT detection of intense, hard emission between 0.1 - 200 GeV to TeV energies. The flux of these observed gamma-rays is significantly higher than that theoretically expected from hadronic interactions of galactic cosmic rays with the solar atmosphere. More importantly, spectral slope of Fermi and HAWC observed gamma ray energy spectra differ significantly from that of galactic cosmic rays casting doubt on the prevailing galactic cosmic ray ancestry model of solar disk gamma rays. In this letter, we argue that the quiet sun can accelerate cosmic rays to TeV energies with an appropriate flux level in the solar chromosphere, as the solar chromosphere in its quiet state probably possesses the required characteristics to accelerate cosmic rays to TeV energies. Consequently, the mystery of the origin of observed gamma rays from the solar disc can be resolved consistently through the hadronic interaction of these cosmic rays with solar matter above the photosphere in a quiet state. The upcoming IceCube-Gen2 detector should be able to validate the proposed model in future through observation of TeV muon neutrino flux from the solar disk. The proposed idea should have major implications on the origin of galactic cosmic rays.

We present the integrated 3-point correlation functions (3PCF) involving both the cosmic shear and the galaxy density fields. These are a set of higher-order statistics that describe the modulation of local 2-point correlation functions (2PCF) by large-scale features in the fields, and which are easy to measure from galaxy imaging surveys. Based on previous works on the shear-only integrated 3PCF, we develop the theoretical framework for modelling 5 new statistics involving the galaxy field and its cross-correlations with cosmic shear. Using realistic galaxy and cosmic shear mocks from simulations, we determine the regime of validity of our models based on leading-order standard perturbation theory with an MCMC analysis that recovers unbiased constraints of the amplitude of fluctuations parameter $A_s$ and the linear and quadratic galaxy bias parameters $b_1$ and $b_2$. Using Fisher matrix forecasts for a DES-Y3-like survey, relative to baseline analyses with conventional 3$\times$2PCFs, we find that the addition of the shear-only integrated 3PCF can improve cosmological parameter constraints by $20-40\%$. The subsequent addition of the new statistics introduced in this paper can lead to further improvements of $10-20\%$, even when utilizing only conservatively large scales where the tree-level models are valid. Our results motivate future work on the galaxy and shear integrated 3PCFs, which offer a practical way to extend standard analyses based on 3$\times$2PCFs to systematically probe the non-Gaussian information content of cosmic density fields.

M64, often called the "Evil Eye" galaxy, is unique among local galaxies. Beyond its dramatic, dusty nucleus, it also hosts an outer gas disk that counter-rotates relative to its stars. The mass of this outer disk is comparable to the gas content of the Small Magellanic Cloud (SMC), prompting the idea that it was likely accreted in a recent minor merger. Yet, detailed follow-up studies of M64's outer disk have shown no evidence of such an event, leading to other interpretations, such as a "flyby" interaction with the distant diffuse satellite Coma P. We present Subaru Hyper Suprime-Cam observations of M64's stellar halo, which resolve its stellar populations and reveal a spectacular radial shell feature, oriented $\sim$30$^{\circ}$ relative to the major axis and along the rotation axis of the outer gas disk. The shell is $\sim$45 kpc southeast of M64, while a similar but more diffuse plume to the northwest extends to $>$100 kpc. We estimate a stellar mass and metallicity for the southern shell of $M_{\star} {=} 1.80~{\pm}~0.54{\times}10^8~M_{\odot}$ and [M/H] $=$ $-$1.0, respectively, and a similar mass of $1.42~{\pm}~0.71{\times}10^8 M_{\odot}$ for the northern plume. Taking into account the accreted material in M64's inner disk, we estimate a total stellar mass for the progenitor satellite of $M_{\rm \star,prog}~{\simeq}~5{\times}10^8~M_{\odot}$. These results suggest that M64 is in the final stages of a minor merger with a gas-rich satellite strikingly similar to the SMC, in which M64's accreted counter-rotating gas originated, and which is responsible for the formation of its dusty inner star-forming disk.

We discuss the production of primordial gravitational waves (GW) from radiative inflaton decay during the period of reheating, assuming perturbative decay of the inflaton either into a pair of bosons or fermions, leading to successful reheating satisfying constraint from Big Bang nucleosynthesis. Assuming that the inflaton $\phi$ oscillates in a general monomial potential $V(\phi)\propto \phi^n$, which results in a time-dependent inflaton decay width, we show that the resulting stochastic GW background can have optimistic detection prospects, especially in detectors that search for a high-frequency GW spectrum, depending on the choice of $n$ that determines the shape of the potential during reheating. We also discuss how this GW energy density may affect the measurement of $\Delta N_{\text{eff}}$ for bosonic and fermionic reheating scenarios.

We consider the solutions of the Guderley problem, consisting of a converging and diverging hydrodynamic shock wave in an ideal gas with a power law initial density profile. The self-similar solutions, and specifically the reflected shock coefficient, which determines the path of the reflected shock, are studied in detail, for cylindrical and spherical symmetries and for a wide range of values of the adiabatic index and the spatial density exponent. Finally, we perform a comprehensive comparison between the analytic solutions and Lagrangian hydrodynamic simulations, by setting proper initial and boundary conditions. A very good agreement between the analytical solutions and the numerical simulations is obtained. This demonstrates the usefulness of the analytic solutions as a code verification test problem.

Buoyancy-driven turbulent convection leads to a fully compressible flow with a strong top-down asymmetry of first- and second-order statistics when the adiabatic equilibrium profiles of temperature, density and pressure decay very strongly across the convection layer. The growth of this asymmetry and the formation of an increasingly thicker stabilized sublayer with a negative mean convective heat flux at the top of the convection zone is reported here by a series highly resolved three-dimensional direct numerical simulations beyond the Oberbeck-Boussinesq and anelastic limits for dimensionless dissipation numbers $0.1 \le D\le 0.8$ at fixed Rayleigh number $Ra=10^6$ and superadiabaticity. The highly stratified compressible convection regime appears for $D > D_{\rm crit}\approx 0.65$, when density fluctuations collapse to those of pressure; it is characterized by an up to nearly 50\% reduced global turbulent heat transfer and a sparse network of focused thin thermal plumes falling through the top sublayer deep into the bulk.

By adding a matter-coupled dark energy field to Einstein's General Relativity (GR), this paper proves that the dynamical dark energy field can change the frequency of photons from distant galaxies as well as from background radiation of remote Universe. Therefore, when the observed frequency-shift of the photons is entirely attributed to the temporal variation of the cosmic scale factor, the calculated expansion rate of the Universe will be slightly greater than its actual value. The predicted values of the temperature of the cosmic blackbody radiation in the past (future) of the Universe are slightly larger (gradually smaller and smaller) than those in the standard cosmology. Since the blackbody radiation becomes the present cosmic microwave background (CMB) and its present-day temperature is directly estimated according to the Planck's law of blackbody radiation, the measured value of the CMB temperature is independent of whether to consider the scalar field or not.

We investigate how the resonant conversion at a temperature $\bar{T}=25$-$65$ keV of a fraction of the CMB photons into an axion-like majoron affects BBN. The scenario, that assumes the presence of a primordial magnetic field and the subsequent decay of the majorons into neutrinos at $T\approx 1$ eV, has been proposed to solve the $H_0$ tension. We find two main effects. First, since we lose photons to majorons at $\bar{T}$, the baryon to photon ratio is smaller at the beginning of BBN $(T>\bar{T})$ than during decoupling and structure formation ($T\ll \bar{T}$). This relaxes the $2\sigma$ mismatch between the observed deuterium abundance and the one predicted by the standard $\Lambda$CDM model. Second, since the conversion implies a sudden drop in the temperature of the CMB during the final phase of BBN, it interrupts the synthesis of lithium and beryllium and reduces their final abundance, possibly alleviating the lithium problem.