Papers for Thursday, Mar 16 2023

The widely known relation between stellar mass and gas metallicity (mass-metallicity relation, MZR) in galaxies is often ascribed to the higher capability of more massive systems to retain metals against the action of galactic outflows. In this scenario the stellar mass would simply be an indirect proxy of the dynamical mass or of the gravitational potential. We test this scenario by using a sample of more than one thousand star-forming galaxies from the MaNGA survey for which dynamical masses have been accurately determined. By using three different methods (average dispersion, Partial Correlation Coefficients, Random Forest) we unambiguously find that the gas metallicity depends primarily and fundamentally on the stellar mass. Once the dependence on stellar mass is taken into account, there is little or no dependence on either dynamical mass or gravitational potential (and, if anything, the metallicity dependence on the latter quantities is inverted). Our result indicates that the MZR is not caused by the retention of metals in more massive galaxies. The direct, fundamental dependence of metallicity on stellar mass suggests the much simpler scenario in which the MZR is just a consequence of the stellar mass being proportional to the integral of metals production in the galaxy.

A joint hadronic model is shown to quantitatively explain the observations of diffuse radio emission from galaxy clusters in the form of minihalos, giant halos, relics, and their hybrid, transitional stages. Cosmic-ray diffusion of order $D\sim 10^{31\text{--}32}\text{ cm}^2\text{ s}^{-1}$, inferred independently from relic energies, the spatial variability of giant-halo spectra, and the spectral evolution of relics, reproduces the observed spatio-spectral distributions, explains the recently discovered mega-halos as enhanced peripheral magnetisation, and quenches electron (re)acceleration by weak shocks or turbulence. For instance, the hard-to-soft evolution along secondary-electron diffusion explains both the soft spectra in most halo peripheries and relic downstreams, and the hard spectra in most halo centres and relic edges, where the photon index can reach $\alpha\simeq -0.5$ regardless of the Mach number $\mathcal{M}$ of the coincident shock. Such spatio-spectral modeling, recent $\gamma$-ray observations, and additional accumulated evidence are thus shown to support a previous claim (Keshet 2010) that the seamless transitions among minihalos, giant halos, and relics, their similar energetics, integrated spectra, and delineating discontinuities, the inconsistent $\mathcal{M}$ inferred from radio vs. X-rays in leptonic models, and additional observations, all indicate that these diffuse radio phenomena are manifestations of the same cosmic-ray ion population, with no need to invoke less natural alternatives.

Strong magnetic fields in the cores of stars are expected to significantly modify the behavior of gravity waves, which is likely the origin of suppressed dipole modes observed in many red giants. However, a detailed understanding of how such fields alter the spectrum and spatial structure of magnetogravity waves has been elusive. For a dipole field, we analytically characterize the horizontal eigenfunctions of magnetogravity modes, assuming that the wavevector is primarily radial. For axisymmetric modes ($m=0$), the magnetogravity wave eigenfunctions become Hough functions, and they have a radial turning point for sufficiently strong magnetic fields. For non-axisymmetric modes ($m\neq0$), the interaction between the discrete $g$ mode spectrum and a continuum of Alfv\'en waves produces nearly discontinuous features in the fluid displacements at critical latitudes associated with a singularity in the fluid equations. We find that magnetogravity modes cannot propagate in regions with sufficiently strong magnetic fields, instead becoming evanescent. When encountering strong magnetic fields, ingoing gravity waves are likely refracted into outgoing slow magnetic waves. These outgoing waves approach infinite radial wavenumbers, which are likely to be damped efficiently. However, it may be possible for a small fraction of the wave power to escape the stellar core as pure Alfv\'en waves or magnetogravity waves confined to a very narrow equatorial band. The artificially sharp features in the WKB-separated solutions suggest the need for global mode solutions which include small terms neglected in our analysis.

We report detections of the [OIII]$\lambda$4364 auroral emission line for 16 galaxies at z=2.1-8.7, measured from JWST/NIRSpec observations obtained as part of the Cosmic Evolution Early Release Science (CEERS) survey program. We combine this CEERS sample with 9 objects from the literature at z=4-9 with auroral-line detections from JWST/NIRSpec and 21 galaxies at z=1.4-3.7 with auroral-line detections from ground-based spectroscopy. We derive electron temperature T_e and direct-method oxygen abundances for the combined sample of 46 star-forming galaxies at z=1.4-8.7. We use these measurements to construct the first high-redshift empirical T_e-based metallicity calibrations for the strong-line ratios [OIII]/H$\beta$, [OII]/H$\beta$, R23=([OIII]+[OII])/H$\beta$, [OIII]/[OII], and [NeIII]/[OII]. These new calibrations are valid over 12+log(O/H)=7.0-8.4 and can be applied to samples of star-forming galaxies at z=2-9, leading to an improvement in the accuracy of metallicity determinations at Cosmic Noon and in the Epoch of Reionization. The high-redshift strong-line relations are offset from calibrations based on typical $z\sim0$ galaxies or HII regions, reflecting the known evolution of ionization conditions between $z\sim0$ and $z\sim2$. Deep spectroscopic programs with JWST/NIRSpec promise to improve statistics at the low and high ends of the metallicity range covered by the current sample, as well as improve the detection rate of [NII]$\lambda$6585 to allow the future assessment of N-based indicators. These new high-redshift calibrations will enable accurate characterizations of metallicity scaling relations at high redshift, improving our understanding of feedback and baryon cycling in the early universe.

Using recent empirical constraints on the dark matter halo--galaxy--supermassive black hole (SMBH) connection from $z=0-7$, we infer how undermassive, typical, and overmassive SMBHs contribute to the quasar luminosity function (QLF) at $z=6$. We find that beyond $L_\mathrm{bol} = 5 \times 10^{46}$ erg/s, the $z=6$ QLF is dominated by SMBHs that are at least 0.3 dex above the $z=6$ median $M_\bullet-M_*$ relation. The QLF is dominated by typical SMBHs (i.e., within $\pm 0.3$ dex around the $M_\bullet-M_*$ relation) at $L_\mathrm{bol} \lesssim 10^{45}$ erg/s. At $z\sim 6$, the intrinsic $M_\bullet-M_*$ relation for all SMBHs is slightly steeper than the $z=0$ scaling, with a similar normalization at $M_* \sim 10^{11} M_\odot$. We also predict the $M_\bullet-M_*$ relation for $z=6$ bright quasars selected by different bolometric luminosity thresholds, finding very good agreement with observations. For quasars with $L_\mathrm{bol} > 3 \times 10^{46}$ ($10^{48}$) erg/s, the scaling relation is shifted upwards by $\sim0.35$ (1.0) dex for $10^{11} M_\odot$ galaxies. To accurately measure the intrinsic $M_\bullet-M_*$ relation, it is essential to include fainter quasars with $L_\mathrm{bol} \lesssim 10^{45}$ erg/s. At high redshifts, low-luminosity quasars are thus the best targets for understanding typical formation paths for SMBHs in galaxies.

Rapidly rotating red giant stars are astrophysically interesting but rare. In this paper we present a catalog of 3217 active red giant candidates in the APOGEE DR16 survey. We use a control sample in the well-studied Kepler fields to demonstrate a strong relationship between rotation and anomalies in the spectroscopic solution relative to typical giants. Stars in the full survey with similar solutions are identified as candidates. We use vsin\textiti measurements to confirm 50+/- 1.2% of our candidates as definite rapid rotators, compared to 4.9+/-0.2% in the Kepler control sample. In both the Kepler control sample and a control sample from DR16, we find that there are 3-4 times as many giants rotating with 5 < vsini < 10 km s$^{-1}$ compared to vsini > 10 km s$^{-1}$, the traditional threshold for anomalous rotation for red giants. The vast majority of intermediate rotators are not spectroscopically anomalous. We use binary diagnostics from APOGEE and \textit{Gaia} to infer a binary fraction of 73+/-2.4%. We identify a significant bias in the reported metallicity for candidates with complete spectroscopic solutions, with median offsets of 0.37 dex in [M/H] from a control sample. As such, up to 10% of stars with reported [M/H]<-1 are not truly metal poor. Finally, we use Gaia data to identify a sub-population of main sequence photometric binaries erroneously classified as giants.

EMCCDs are efficient imaging devices for low surface brightness UV astronomy from space. The large amplification allows photon counting, the detection of events versus non-events. This paper provides the statistics of the observation process, the photon-counting process, the amplification, process, and the compression. The expression for the signal-to-noise of photon counting is written in terms of the polygamma function. The optimal exposure time is a function of the clock-induced charge. The exact distribution of amplification process is a simple-to-compute powered matrix. The optimal cutoff for comparing to the read noise is close to a strong function of the read noise and a weak function of the electron-multiplying gain and photon rate. A formula gives the expected compression rate.

We use molecular line data from ALMA, SMA, JCMT, and NANTEN2 to study the multi-scale ($\sim$15-0.005 pc) velocity statistics in the massive star formation region NGC 6334. We find that the non-thermal motions revealed by the velocity dispersion function (VDF) stay supersonic over scales of several orders of magnitudes. The multi-scale non-thermal motions revealed by different instruments do not follow the same continuous power-law, which is because the massive star formation activities near central young stellar objects have increased the non-thermal motions in small-scale and high-density regions. The magnitudes of VDFs vary in different gas materials at the same scale, where the infrared dark clump N6334S in an early evolutionary stage shows a lower level of non-thermal motions than other more evolved clumps due to its more quiescent star formation activity. We find possible signs of small-scale-driven (e.g., by gravitational accretion or outflows) supersonic turbulence in clump N6334IV with a three-point VDF analysis. Our results clearly show that the scaling relation of velocity fields in NGC 6334 deviates from a continuous and universal turbulence cascade due to massive star formation activities.

Several planetary formation models have been proposed to explain the observed abundance and variety of compositions of super-Earths and mini-Neptunes. In this context, multitransiting systems orbiting low-mass stars whose planets are close to the radius valley are benchmark systems, which help to elucidate which formation model dominates. We report the discovery, validation, and initial characterization of one such system, TOI-2096, composed of a super-Earth and a mini-Neptune hosted by a mid-type M dwarf located 48 pc away. We first characterized the host star by combining different methods. Then, we derived the planetary properties by modeling the photometric data from TESS and ground-based facilities. We used archival data, high-resolution imaging, and statistical validation to support our planetary interpretation. We found that TOI-2096 corresponds to a dwarf star of spectral type M4. It harbors a super-Earth (R$\sim1.2 R_{\oplus}$) and a mini-Neptune (R$\sim1.90 R_{\oplus}$) in likely slightly eccentric orbits with orbital periods of 3.12 d and 6.39 d, respectively. These orbital periods are close to the first-order 2:1 mean-motion resonance (MMR), which may lead to measurable transit timing variations (TTVs). We computed the expected TTVs amplitude for each planet and found that they might be measurable with high-precision photometry delivering mid-transit times with accuracies of $\lesssim$2 min. Moreover, measuring the planetary masses via radial velocities (RVs) is also possible. Lastly, we found that these planets are among the best in their class to conduct atmospheric studies using the James Webb Space Telescope (JWST). The properties of this system make it a suitable candidate for further studies, particularly for mass determination using RVs and/or TTVs, decreasing the scarcity of systems that can be used to test planetary formation models around low-mass stars.

We report the results of ALMA 1-2 kpc-resolution, three rotational transition line (J=2-1, J=3-2, and J=4-3) observations of multiple dense molecular gas tracers (HCN, HCO$^{+}$, and HNC) for ten nearby (ultra)luminous infrared galaxies ([U]LIRGs). Following the matching of beam sizes to 1-2 kpc for each (U)LIRG, the high-J to low-J transition-line flux ratios of each molecule and the emission line flux ratios of different molecules at each J transition are derived. We conduct RADEX non-LTE model calculations and find that, under a wide range of gas density and kinetic temperature, the observed HCN-to-HCO$^{+}$ flux ratios in the overall (U)LIRGs are naturally reproduced with enhanced HCN abundance compared to HCO$^{+}$. Thereafter, molecular gas properties are constrained primarily through the use of HCN and HCO$^{+}$ data and the adoption of fiducial values for the HCO$^{+}$ column density and HCN-to-HCO$^{+}$ abundance ratio. We quantitatively confirm the following: (1) Molecular gas at the (U)LIRGs' nuclei is dense ($\gtrsim$10$^{3-4}$ cm$^{-3}$) and warm ($\gtrsim$100 K). (2) Molecular gas density and temperature in nine ULIRGs' nuclei are significantly higher than that of one LIRG's nucleus. (3) Molecular gas in starburst-dominated sources tends to be less dense and cooler than ULIRGs with luminous AGN signatures. For six selected sources, we also apply a Bayesian approach by freeing all parameters and support the above main results. Our ALMA 1-2 kpc resolution, multiple transition-line data of multiple molecules are a very powerful tool for scrutinizing the properties of molecular gas concentrated around luminous energy sources in nearby (U)LIRGs' nuclei.

We measure precise orbits and dynamical masses and derive age constraints for six confirmed and one candidate Sirius-like systems, including the Hyades member HD 27483. Our orbital analysis incorporates radial velocities, relative astrometry, and Hipparcos-Gaia astrometric accelerations. We constrain the main-sequence lifetime of a white dwarf's progenitor from the remnant's dynamical mass and semi-empirical initial-final mass relations and infer the cooling age from mass and effective temperature. We present new relative astrometry of HD 27483 B from Keck/NIRC2 observations and archival HST data, and obtain the first dynamical mass of ${0.798}_{-0.041}^{+0.10}$ $M_{\odot}$, and an age of ${450}_{-180}^{+570}$ Myr, consistent with previous age estimates of Hyades. We also measure precise dynamical masses for HD 114174 B ($0.591 \pm 0.011$ $M_{\odot}$) and HD 169889 B (${0.526}_{-0.037}^{+0.039}$ $M_{\odot}$), but their age precisions are limited by their uncertain temperatures. For HD 27786 B, the unusually small mass of $0.443 \pm 0.012$ $M_{\odot}$ suggests a history of rapid mass loss, possibly due to binary interaction in its progenitor's RGB phase. The orbits of HD 118475 and HD 136138 from our RV fitting are overall in good agreement with Gaia DR3 astrometric two-body solutions, despite moderate differences in the eccentricity and period of HD 136138. The mass of ${0.580}_{-0.039}^{+0.052}$ $M_{\odot}$ for HD 118475 B and a speckle imaging non-detection confirms that the companion is a white dwarf. Our analysis shows examples of a rich number of precise WD dynamical mass measurements enabled by Gaia DR3 and later releases, which will improve empirical calibrations of the white dwarf initial-final mass relation.

We analyze a cosmological model featuring an interaction between dark energy and dark matter in light of the measurements of the Cosmic Microwave Background released by three independent experiments: the most recent data by the Planck satellite and the Atacama Cosmology Telescope, and WMAP (9-year data). We show that different combinations of the datasets provide similar results, always favoring an interacting dark sector with a $95\%$~CL significance in the majority of the cases. Remarkably, such a preference remains consistent when cross-checked through independent probes, while always yielding a value of the expansion rate $H_0$ consistent with the local distance ladder measurements. We investigate the source of this preference by scrutinizing the angular power spectra of temperature and polarization anisotropies as measured by different experiments.

Solar surface magneto-convection appears as granulation pattern that impacts spectral lines in terms of both shape and wavelength. Such induced effects also tend to vary over the observed solar disc because of the changing observation angle and, thus, the changing observation height as well. Centre-to-limb observations of the resolved Sun offer an insight into the variable spectral behaviour across different heliocentric observing positions, providing crucial information about limb darkening, convective velocities, and line profile variability relevant to radial velocity (RV) calculations. Thus, RV measurements and exoplanet transit spectroscopy depend on precise reference templates. We want to provide a spectroscopic centre-to-limb solar atlas at high spectral resolution and high-frequency accuracy. The atlas shall help improve the understanding of the solar atmosphere and convection processes. We performed high-resolution observations of the resolved quiet Sun with a Fourier transform spectrograph at the Institut f\"ur Astrophysik und Geophysik in G\"ottingen. Our dataset contains a wavelength range from 4200\r{A} to 8000\r{A}. We obtained 165 spectra in total, with a spectral resolution of $\Delta \nu$ = 0.024cm$^{-1}$, corresponding to a resolving power $R$ of 700,000 at $\sim$6000\r{A}. We present a centre-to-limb solar atlas containing 14 heliocentric positions. To check for consistency, we investigated the FeI ~6175\r{A} absorption line and compared our line profiles with previous centre-to-limb observations and also with simulations. The line profile and also the bisector profiles are generally consistent with previous observations, but we have identified differences to model line profiles, especially close to the solar limb.

In Paper I we showed that clumps in high-redshift galaxies, having a high star formation rate density (\Sigma_SFR), produce disks with two tracks in the [Fe/H]-[\alpha/Fe] chemical space, similar to that of the Milky Way's (MW's) thin + thick disks. Here we investigate the effect of clumps on the bulge's chemistry. The chemistry of the MW's bulge is comprised of a single track with two density peaks separated by a trough. We show that the bulge chemistry of an N-body + smoothed particle hydrodynamics clumpy simulation also has a single track. Star formation within the bulge is itself in the high-\Sigma_SFR clumpy mode, which ensures that the bulge's chemical track follows that of the thick disk at low [Fe/H] and then extends to high [Fe/H], where it peaks. The peak at low metallicity instead is comprised of a mixture of in-situ stars and stars accreted via clumps. As a result, the trough between the peaks occurs at the end of the thick disk track. We find that the high-metallicity peak dominates near the mid-plane and declines in relative importance with height, as in the MW. The bulge is already rapidly rotating by the end of the clump epoch, with higher rotation at low [\alpha/Fe]. Thus clumpy star formation is able to simultaneously explain the chemodynamic trends of the MW's bulge, thin + thick disks and the Splash.

Weak lensing provides a direct way of mapping the density distribution in the universe. To reconstruct the density field from the shear catalog, an important step is to build the shear field from the shear catalog, which can be quite nontrivial due to the inhomogeneity of the background galaxy distribution and the shape noise. We propose the PDF-Folding method as a statistically optimal way of reconstructing the shear field. It is an extention of the PDF-SYM method, which is previously designed for optimizing the stacked shear signal as well as the shear-shear correlation. PDF-Folding does not require smoothing kernels as in traditional methods, therefore it suffers less information loss on small scales, and avoids possible biases due to the spatial variation of shear on the scale of the kernel. We show with analytic reasoning as well as numerical examples that the new method can reach the optimal signal-to-noise ratio on the reconstructed shear map under general observing conditions, i.e., with inhomogeneous background densities or masks. We also show the performance of the new method with real shear data around foreground galaxy clusters.

Condensation and sublimation of ices at the surface of the planet is a key part of both the Martian H$_2$O and CO$_2$ cycles, either from a seasonal or diurnal aspect. While most of the ice is located within the polar caps, surface frost is known to be formed during nighttime down to equatorial latitudes. Here, we use data from the Emirates Mars Infrared Spectrometer (EMIRS) onboard the Emirates Mars Mission (EMM) to monitor the diurnal and seasonal evolution of the ices at the surface of Mars over almost one Martian year. The unique local time coverage provided by the instrument allows us to observe the apparition of equatorial CO$_2$ frost in the second half of the Martian night around the equinoxes, to its sublimation at sunrise.

A catalog of more than 43,000 M giant stars has been selected by Li et al. from the ninth data release of LAMOST. Using the data-driven method SLAM, we obtain the stellar parameters (Teff, logg, [M/H], [$\alpha$/M]) for all the M giant stars with uncertainties of 57 K, 0.25 dex, 0.16 dex and 0.06 dex at SNR > 100, respectively. With those stellar parameters, we constrain the absolute magnitude in K-band, which brings distance with relative uncertainties around 25% statistically. Radial velocities are also calculated by applying cross correlation on the spectra between 8000 A $\AA$ and 8950 A $\AA$ with synthetic spectra from ATLAS9, which covers the Ca II triplet. Comparison between our radial velocities and those from APOGEE DR17 and Gaia DR3 shows that our radial velocities have a system offset and dispersion around 1 and 4.6 km s$^{-1}$, respectively. With the distances and radial velocities combining with the astrometric data from Gaia DR3, we calculate the full 6D position and velocity information, which are able to be used for further chemo-dynamic studies on the disk and substructures in the halo, especially the Sagittarius Stream.

We have been conducting a exoplanet search survey using Bohyunsan Observatory Echelle Spectrograph (BOES) for the last 18 years. We present the detection of exoplanet candidate in orbit around HD 18438 from high-precision radial velocity (RV) mesurements. The target was already reported in 2018 (Bang et al. 2018). They conclude that the RV variations with a period of 719 days are likely to be caused by the pulsations because the Lomb-Scargle periodogram of HIPPARCOS photometric and Ha EW variations for HD 18438 show peaks with periods close to that of RV variations and there were no correlations between bisectors and RV measurements. However, the data were not sufficient to reach a firm conclusion. We obtained more RV data for four years. The longer time baseline yields a more accurate determination with a revised period of 803 +/- 5 days and the planetary origin of RV variations with a minimum planetary companion mass of 21 +/- 1 MJup. Our current estimate of the stellar parameters for HD 18438 makes it currently the largest star with a planetary companion.

High-resolution X-ray observations offer a unique tool for probing the still elusive connection between galaxy mergers and active galactic nuclei (AGNs). We present an analysis of nuclear X-ray emission in an optically selected sample of 92 close galaxy pairs (with projected separations $\lesssim 20$ kpc and line-of-sight velocity offsets $<$ 500 km s$^{-1}$) at low redshift ($\bar{z} \sim 0.07$), based on archival Chandra observations. The parent sample of galaxy pairs is constructed without imposing an optical classification of nuclear activity, thus is largely free of selection effect for or against the presence of an AGN. Nor is this sample biased for or against gas-rich mergers. An X-ray source is detected in 70 of the 184 nuclei, giving a detection rate of $38\%^{+5\%}_{-5\%}$, down to a 0.5-8 keV limiting luminosity of $\lesssim 10^{40}\rm~erg~s^{-1}$. The detected and undetected nuclei show no systematic difference in their host galaxy properties such as galaxy morphology, stellar mass and stellar velocity dispersion. When potential contamination from star formation is avoided (i.e., $L_{\rm 2-10~keV} > 10^{41}\rm~erg~s^{-1}$), the detection rate becomes $18\%^{+3\%}_{-3\%}$ (32/184), which shows no excess compared to the X-ray detection rate of a comparison sample of optically classified single AGNs. The fraction of pairs containing dual AGN is only $2\%^{+2\%}_{-2\%}$. Moreover, most nuclei at the smallest projected separations probed by our sample (a few kpc) have an unexpectedly low apparent X-ray luminosity and Eddington ratio, which cannot be solely explained by circumnuclear obscuration. These findings suggest that close galaxy interaction is not a sufficient condition for triggering a high level of AGN activity.

Mid-infrared (mid-IR) variability in young stellar objects (YSOs) is driven by several physical mechanisms, which produce a variety of amplitudes and light curve shapes. One of these mechanisms, variable disk accretion is predicted by models of episodic accretion to drive secular variability, including in the mid-IR. Because the largest accretion bursts are rare, adding new objects to the YSO eruptive variable class aids our understanding of the episodic accretion phenomenon and its possible impact on stellar and planetary formation. A previous analysis of 6.5 yr of NeoWISE light curves (3-5 $\mu$m) of ~7000 nearby YSOs found an increase in the fraction of variability and variability amplitude for objects at younger stages of evolution. To help interpret these light curves, we have obtained low- and high-resolution near-IR spectra of 78 objects from this sample of YSOs. In this work, we present the analysis of nine nearby YSOs (d$<$1 kpc) that show the characteristics of known classes of eruptive variable YSOs. We find one FUor-like source, one EX Lupi-type object, and six YSOs with mixed characteristics, or V1647 Ori-like objects. The varied characteristics observed in our sample are consistent with recent discoveries of eruptive YSOs. We discuss how a wide range in YSO outburst parameters (central mass, maximum accretion rate during outburst, evolutionary stage and/or instability leading to the outburst) may play a significant role in the observed spectro-photometric properties of YSO outbursts.

We present an analysis of the theoretical and observed light curve parameters of the fundamental mode (FU) classical Cepheids in the Magellanic Clouds in $V$- and $I$- photometric bands. The state-of-the-art 1D non-linear radial stellar pulsation (RSP) code in MESA (\textsc{mesa-rsp}) has been utilized to generate the theoretical light curves using four sets of convection parameters. Theoretical light curves with two chemical compositions: $Z=0.008$ and $Z=0.004$ appropriate for the Large Magellanic Cloud (LMC) and Small Magellanic Cloud (SMC), respectively, covered a wide range of periods ($3<P (\rm{d})<32$). The observed light curves are taken from the OGLE-IV database. We compare theoretical and observed Fourier parameters (FPs), and investigate the period-luminosity (PL), period-colour (PC), and amplitude-colour (AC) relations as a function of pulsation phase for short ($\log{P}<1$), long ($\log{P}>1$) and all periods. The multiphase relations obtained from theoretical and observed light curves in the PL/PC/AC plane are found to be dynamic in nature, with the effect more pronounced at $\Phi \sim 0.75-0.85$. Furthermore, a contrasting behaviour of the theoretical/observed multiphase PL and PC relations between the short and long periods has been found for both LMC and SMC. The analysis shows that multiphase PL relations are more stringent to test the models with observations over the FPs. Distances to the LMC/SMC determined using long period Cepheids are found to be in good agreement with the literature values when the term $R_{21}$ is added to the PL relation.

Very light pseudoscalar fields, often referred to as axions, are compelling dark matter candidates and can potentially be detected through their coupling to the electromagnetic field. Recently a novel detection technique using the cosmic microwave background (CMB) was proposed, which relies on the fact that the axion field oscillates at a frequency equal to its mass in appropriate units, leading to a time-dependent birefringence. For appropriate oscillation periods this allows the axion field at the telescope to be detected via the induced sinusoidal oscillation of the CMB linear polarization. We search for this effect in two years of POLARBEAR data. We do not detect a signal, and place a median $95 \%$ upper limit of $0.65 ^\circ$ on the sinusoid amplitude for oscillation frequencies between $0.02\,\text{days}^{-1}$ and $0.45\,\text{days}^{-1}$, which corresponds to axion masses between $9.6 \times 10^{-22} \, \text{eV}$ and $2.2\times 10^{-20} \,\text{eV}$. Under the assumptions that 1) the axion constitutes all the dark matter and 2) the axion field amplitude is a Rayleigh-distributed stochastic variable, this translates to a limit on the axion-photon coupling $g_{\phi \gamma} < 2.4 \times 10^{-11} \,\text{GeV}^{-1} \times ({m_\phi}/{10^{-21} \, \text{eV}})$.

We present optical photometric and spectroscopic studies of three supernovae (SNe) SN 2013bz, PSN J0910+5003 and ASASSN-16ex. UV-optical photometric data of ASASSN-16ex obtained with Swift-UVOT are also analyzed. These objects were initially classified as 09dc-like type Ia SNe. The decline rate parameters ($\Delta m_{15}(B)_{true}$) are derived as 0.92 $\pm$ 0.04 (SN 2013bz), 0.70 $\pm$ 0.05 (PSN J0910+5003) and 0.73 $\pm$ 0.03 (ASASSN-16ex). The estimated $B$ band absolute magnitudes at maximum: $-$19.61 $\pm$ 0.20 mag for SN 2013bz, $-$19.44 $\pm$ 0.20 mag for PSN J0910+5003 and $-$19.78 $\pm$ 0.20 mag for ASASSN-16ex indicate that all the three objects are relatively bright. The peak bolometric luminosities for these objects are derived as $\log L_\text{bol}^\text{max} =$ 43.38 $\pm$ 0.07 erg s$^{-1}$, 43.26 $\pm$ 0.07 erg s$^{-1}$ and 43.40 $\pm$ 0.06 erg s$^{-1}$, respectively. The spectral and velocity evolution of SN 2013bz is similar to a normal SN Ia, hence it appears to be a luminous, normal type Ia supernova. On the other hand, the light curves of PSN J0910+5003 and ASASSN-16ex are broad and exhibit properties similar to 09dc-like SNe Ia. Their spectroscopic evolution shows similarity with 09dc-like SNe, strong CII lines are seen in the pre-maximum spectra of these two events. Their photospheric velocity evolution is similar to SN 2006gz. Further, in the UV bands, ASASSN-16ex is very blue like other 09dc-like SNe Ia.

In recent years, scientific CMOS (sCMOS) sensors have been vigorously developed and have outperformed CCDs in several aspects: higher readout frame rate, higher radiation tolerance, and higher working temperature. For silicon image sensors, image lag will occur when the charges of an event are not fully transferred inside pixels. It can degrade the image quality for optical imaging, and deteriorate the energy resolution for X-ray spectroscopy. In this work, the image lag of a sCMOS sensor is studied. To measure the image lag under low-light illumination, we constructed a new method to extract the image lag from X-ray photons. The image lag of a customized X-ray sCMOS sensor GSENSE1516BSI is measured, and its influence on X-ray performance is evaluated. The result shows that the image lag of this sensor exists only in the immediately subsequent frame and is always less than 0.05% for different incident photon energies and under different experimental conditions. The residual charge is smaller than 0.5 e- with the highest incident photon charge around 8 ke-. Compared to the readout noise level around 3 e-, the image lag of this sensor is too small to have a significant impact on the imaging quality and the energy resolution. The image lag shows a positive correlation with the incident photon energy and a negative correlation with the temperature. However, it has no dependence on the gain setting and the integration time. These relations can be explained qualitatively by the non-ideal potential structure inside the pixels. This method can also be applied to the study of image lag for other kinds of imaging sensors.

Radio mini-halos are clouds of diffuse, low surface brightness synchrotron emission that surround the Brightest Cluster Galaxy (BCG) in massive cool-core galaxy clusters. In this paper, we use third generation calibration (3GC), also called direction-dependent (DD) calibration, and point source subtraction on MeerKAT extragalactic continuum data. We calibrate and image archival MeerKAT L-band observations of a sample of five galaxy clusters (ACO 1413, ACO 1795, ACO 3444, MACS J1115.8+0129, MACS J2140.2-2339). We use the CARACal pipeline for direction-independent (DI) calibration, DDFacet and killMS for 3GC, followed by visibility-plane point source subtraction to image the underlying mini-halo without bias from any embedded sources. Our 3GC process shows a drastic improvement in artefact removal, to the extent that the local noise around severely affected sources was halved and ultimately resulted in a 7\% improvement in global image noise. Thereafter, using these spectrally deconvolved Stokes I continuum images, we directly measure for four mini-halos the flux density, radio power, size and in-band integrated spectra. Further to that, we show the in-band spectral index maps of the mini-halo (with point sources). We present a new mini-halo detection hosted by MACS J2140.2-2339, having flux density $S_{\rm 1.28\,GHz} = 2.61 \pm 0.31$ mJy, average diameter 296 kpc and $\alpha^{\rm 1.5\,GHz}_{\rm 1\,GHz} = 1.21 \pm 0.36$. We also found a $\sim$100 kpc southern extension to the ACO 3444 mini-halo which was not detected in previous VLA L-band observations. Our description of MeerKAT wide-field, wide-band data reduction will be instructive for conducting further mini-halo science.

The Gas Pixel Detector is a gas detector, sensitive to the polarization of X-rays, currently flying on-board IXPE - the first observatory dedicated to X-ray polarimetry. It detects X-rays and their polarization by imaging the ionization tracks generated by photoelectrons absorbed in the sensitive volume, and then reconstructing the initial direction of the photoelectrons. The primary ionization charge is multiplied and ultimately collected on a finely-pixellated ASIC specifically developed for X-ray polarimetry. The signal of individual pixels is processed independently and gain variations can be substantial, of the order of 20%. Such variations need to be equalized to correctly reconstruct the track shape, and therefore its polarization direction. The method to do such equalization is presented here and is based on the comparison between the mean charge of a pixel with respect to the other pixels for equivalent events. The method is shown to finely equalize the response of the detectors on board IXPE, allowing a better track reconstruction and energy resolution, and can in principle be applied to any imaging detector based on tracks.

Membership studies characterising open clusters with Gaia data, most using DR2, are so far limited at magnitude G = 18 due to astrometric uncertainties at the faint end. Our goal is to extend current open cluster membership lists with faint members and to characterise the low-mass end, which members are important for many applications, in particular for ground-based spectroscopic surveys. We use a deep neural network architecture to learn the distribution of highly reliable open cluster member stars around known clusters. After that, we use the trained network to estimate new open cluster members based on their similarities in a high-dimensional space, five-dimensional astrometry plus the three photometric bands. Due to the improved astrometric precisions of Gaia DR3 with respect to DR2, we are able to homogeneously detect new faint member stars (G > 18) for the known open cluster population. Our methodology can provide extended membership lists for open clusters down to the limiting magnitude of Gaia, which will enable further studies to characterise the open cluster population, e.g. estimation of their masses, or their dynamics. These extended membership lists are also ideal target lists for forthcoming ground-based spectroscopic surveys.

$\epsilon$ Eridani is a highly active young K2 star with an activity cycle of about three years established using Ca II H & K line index measurements (S_MWO). This relatively short cycle has been demonstrated to be consistent with X-ray and magnetic flux measurements. Recent work suggested a change in the cyclic behaviour. Here we report new X-ray flux and S_MWO measurements and also include S_MWO measurements from the historical Mount Wilson program. This results in an observational time baseline of over 50 years for the S_MWO data and of over 7 years in X-rays. Moreover, we include Ca II infrared triplet (IRT) index measurements (S_IRT) from 2013-2022 in our study. With the extended X-ray data set, we can now detect the short cycle for the first time using a periodogram analysis. Near-simultaneous S_MWO data and X-ray fluxes, which are offset by 20 days at most, are moderately strongly correlated when only the lowest activity state (concerning short-term variability) is considered in both diagnostics. In the S_MWO data, we find strong evidence for a much longer cycle of about 34 years and an 11-year cycle instead of the formerly proposed $12$-year cycle in addition to the known 3-year cycle. The superposition of the three periods naturally explains the recent drop in S_MWO measurements. The two shorter cycles are also detected in the S_IRT data, although the activity cycles exhibit lower amplitudes in the S_IRT than in the S_MWO data. Finally, the rotation period of $\epsilon$Eri can be found more frequently in the S_MWO as well as in the S_IRT data for times near the minimum of the long cycle. This may be explained by a scenario in which the filling factor for magnetically active regions near cycle maximum is too high to allow for notable short-term variations.

Basaltic V-type asteroids are leftovers from the formation and evolution of differentiated planetesimals. They are thought to originate from mantles and crusts of multiple different parent bodies. Identifying the links between individual V-type asteroids and multiple planetesimals is challenging, especially in the inner part of the main asteroid belt, where the majority of V-type asteroids are expected to have originated from a single planetesimal, namely, (4) Vesta. In this work, we aim to trace the origin of a number of individual V-type asteroids from the inner part of the main asteroid belt. The main goal is to identify asteroids that may not be traced back to (4) Vesta and may therefore originate from other differentiated planetesimals. We performed a 2 Gy backward numerical integration of the orbits of the selected V-type asteroids. For each asteroid, we used 1001 clones to map the effect of orbital uncertainties. In the integration, we use information on physical properties of the considered V-type asteroids such as pole orientation, rotational period, and thermal parameters. The majority of V-types in the inner main belt outside the Vesta family are clearly Vesta fugitives. Two objects, namely, (3307) Athabasca and (17028) 1999 FJ$_{5}$, show no clear dynamical link to (4) Vesta. Together with (809) Lundia (from our previous work), these objects could represent the parent bodies of anomalous HED meteorites such as the Banbura Rockhole. Furthermore, some objects of the low-inclination population cannot be traced back to (4) Vesta within the 2 Gy integration.

An $\alpha\Omega$ dynamo, combining shear and cyclonic convection in the tachocline, is believed to generate the solar cycle. However, this model cannot explain cycles in fast rotators (with minimal shear) or in fully convective stars (no tachocline); analysis of such stars could therefore provide key insights into how these cycles work. We reexamine ASAS data for 15 M dwarfs, 11 of which are presumed fully convective; the addition of newer ASAS-SN data confirms cycles in roughly a dozen of them, while presenting new or revised rotation periods for five. The amplitudes and periods of these cycles follow $A_{\rm cyc} \propto P_{\rm cyc}^{0.94 \pm 0.11}$, with $P_{\rm cyc}/P_{\rm rot} \propto {\rm Ro}^{-1.02 \pm 0.06}$ (where Ro is the Rossby number), very similar to $P_{\rm cyc}/P_{\rm rot} \propto {\rm Ro}^{-0.81 \pm 0.17}$ that we find for 40 previously studied FGK stars, although $P_{\rm cyc}/P_{\rm rot}$ and $\alpha$ are a factor of $\sim$20 smaller in the M stars. The very different $P_{\rm cyc}/P_{\rm rot}$-Ro relationship seen here compared to previous work suggests that two types of dynamo, with opposite Ro dependences, operate in cool stars. Initially, a (likely $\alpha^2$ or $\alpha^2\Omega$) dynamo operates throughout the convective zone in mid-late M and fast rotating FGK stars, but once magnetic breaking decouples the core and convective envelope, a tachocline $\alpha\Omega$ dynamo begins and eventually dominates in older FGK stars. A change in $\alpha$ in the tachocline dynamo generates the fundamentally different $P_{\rm cyc}/P_{\rm rot}$-Ro relationship.

We investigate superfluid dark matter (SFDM), a model that promises to reproduce the successes of both particle dark matter on cosmological scales and those of Modified Newtonian Dynamics (MOND) on galactic scales. But SFDM reproduces MOND only up to a certain distance from the galactic center and only for kinematic observables. Most importantly, it does not affect trajectories of light. We test whether or not this is in conflict with a recent analysis of weak gravitational lensing which has probed accelerations around galaxies at unprecedentedly large radii. This analysis found the data to be close to the prediction of MOND, suggesting they might be difficult to fit with SFDM. To investigate this matter, we solved the equations of motion of the model and compared the result to observational data. Our results show that the SFDM model is incompatible with the weak-lensing observations, at least in its current form.

Cool giant and supergiant stars are among the largest and most luminous stars in the Universe and, therefore, dominate the integrated light of their host galaxies. These stars were extensively studied during last few decades, however their photometric variability and mass loss are still poorly constrained. The atmospheres of evolved stars are characterized by complex dynamics due to different interacting processes, such as convection, pulsation, formation of molecules and dust, and the development of mass loss. These dynamical processes impact the formation of spectral lines producing their asymmetries and Doppler shifts. Thus, by studying the line-profile variations on spatial and temporal scales it is possible to reconstruct atmospheric motions in stars and link them to the photometric variability and mass loss. The tomographic method, which is based on the cross-sectioning through the stellar atmosphere and recovering the velocity field for each atmospheric slice, is an ideal technique for this purpose. We present the tomographic method and its application to spectroscopic and spectro-interferometric observations of giant and supergiant stars as well as to state-of-the-art three-dimensional numerical simulations to constrain their atmospheric motions on spatial and temporal scales and better understand respective mechanisms responsible for their photometric variability and mass loss.

We report photometry of the intranight variability of the dwarf nova RX And in two bands (B and V). The observations are carried out during three nights in November-December 2022 at the 50/70~cm Schmidt telescope of the Rozhen National Astronomical Observatory. The observations indicate that the amplitude of the flickering is about 0.5 mag in B band when the star is in faint state (m_V=13.5), but it is considerably lower (less than 0.1 mag) when the star is bright (m_V=10.9). The mass accretion rate in high state is estimated to be $1.2 \times 10^{-9}$ M$_\odot$ yr$^{-1}$. Combining our data and GAIA distances we find for the mass donor in RX~And spectral type K6V-K7V. The data are available upon request from the authors.

We present a novel approach for the dimensionality reduction of galaxy images by leveraging a combination of variational auto-encoders (VAE) and domain adaptation (DA). We demonstrate the effectiveness of this approach using a sample of low redshift galaxies with detailed morphological type labels from the Galaxy-Zoo DECaLS project. We show that 40-dimensional latent variables can effectively reproduce most morphological features in galaxy images. To further validate the effectiveness of our approach, we utilised a classical random forest (RF) classifier on the 40-dimensional latent variables to make detailed morphology feature classifications. This approach performs similarly to a direct neural network application on galaxy images. We further enhance our model by tuning the VAE network via DA using galaxies in the overlapping footprint of DECaLS and BASS+MzLS, enabling the unbiased application of our model to galaxy images in both surveys. We observed that noise suppression during DA led to even better morphological feature extraction and classification performance. Overall, this combination of VAE and DA can be applied to achieve image dimensionality reduction, defect image identification, and morphology classification in large optical surveys.

The Kuiper belt objects, the Centaurs, and the Jupiter-family comets form an evolutionary continuum of small outer Solar System objects, and their study allows us to gain insight into the history and evolution of the Solar System. Broadband photometry can be used to measure their phase curves, allowing a first-order probe into the surface properties of these objects, though limited telescope time makes measuring accurate phase curves difficult. We make use of serendipitous broadband photometry from the long-baseline, high-cadence ATLAS survey to measure the phase curves for a sample of 18 Kuiper belt objects, Centaurs, and Jupiter-family comets with unprecedentedly large datasets. We find phase curves with previously reported negative slopes become positive with increased data and are thus due to insufficient sampling of the phase curve profile, and not a real physical effect. We search for correlations between phase curve parameters, finding no strong correlations between any parameter pair, consistent with the findings of previous studies. We search for instances of cometary activity in our sample, finding a previously reported outburst by Echeclus and a new epoch of increased activity by Chiron. Applying the main belt asteroid HG1G2 phase curve model to three Jupiter-family comets in our sample with large phase angle spans, we find their slope parameters imply surfaces more consistent with those of carbonaceous main belt asteroids than silicaceous ones.

The Legacy e-MERLIN Multi-band Imaging of Nearby Galaxies survey (LeMMINGs) is a statistically-complete census of nuclear accretion and star formation processes in the local Universe. The LeMMINGs observations at 1.5 and 5 GHz yield angular resolutions on 10s milliarcsecond-scales, with sensitivities of 10s $\mu$Jy. Awarded 810 hours of observing time, the full statistical sample (at 1.5 GHz) plus several studies of individual objects have now been published. Combined with multi-wavelength follow up observations, this survey will provide a unique legacy data set of our Galactic back yard. We present an overview of the LeMMINGs results so far, including the 1.5 GHz sample results and associated Chandra X-ray data. We describe the next steps for LeMMINGs to analyse the 5 GHz survey and produce widefield images to categories all radio sources in the LeMMINGs galaxies.

Despite centuries of rigorous theoretical and observational research, the origin and acceleration mechanism of Galactic Cosmic Rays (GCRs) remain a mystery. In 1949, Fermi proposed a diffusive shock acceleration model that includes a prominent mechanism for GCR acceleration. However, observational evidence, on the other hand, remains elusive. Here, we provided the first apparent verification of GCR acceleration at 1 AU using measurements from the CRIS instrument onboard the ACE spacecraft.

Coronal and interplanetary shock waves produced by coronal mass ejections (CMEs) are major drivers of space-weather phenomena, inducing major changes in the heliospheric radiation environment and directly perturbing the near-Earth environment, including its magnetosphere. A better understanding of how these shock waves evolve from the corona to the interplanetary medium can therefore contribute to improving nowcasting and forecasting of space weather. Early warnings from these shock waves can come from radio measurements as well as coronagraphic observations that can be exploited to characterise the dynamical evolution of these structures. Our aim is to analyse the geometrical and kinematic properties of 32 CME shock waves derived from multi-point white-light and ultraviolet imagery taken by the Solar Dynamics Observatory (SDO), Solar and Heliospheric Observatory (SoHO), and Solar-Terrestrial Relations Observatory (STEREO) to improve our understanding of how shock waves evolve in 3D during the eruption of a CME. We use our catalogue to search for relations between the shock wave's kinematic properties and the flaring activity associated with the underlying genesis of the CME piston. Past studies have shown that shock waves observed from multiple vantage points can be aptly reproduced geometrically by simple ellipsoids. The catalogue of reconstructed shock waves provides the time-dependent evolution of these ellipsoidal parameters. From these parameters, we deduced the lateral and radial expansion speeds of the shocks evolving over time. We compared these kinematic properties with those obtained from a single viewpoint by SoHO in order to evaluate projection effects. Finally, we examined the relationships between the shock wave and the associated flare when the latter was observed on the disc by considering the measurements of soft and hard X-rays.

Thermal emission of neutron stars in soft X-ray transients (SXTs) in a quiescent state is believed to be powered by the heat deposited in the stellar crust due to nuclear reactions during accretion (deep crustal heating paradigm). Confronting observations of SXTs with simulations helps to verify theoretical models of the dense matter in the neutron stars. Usually, such simulations were carried out assuming that the free neutrons and nuclei in the inner crust move together. A recently proposed thermodynamically consistent approach allows for independent motion of the free neutrons. We simulate the thermal evolution of the SXTs within the thermodynamically consistent approach and compare the results with the traditional approach and with observations. For the latter, we consider a collection of quasi-equilibrium thermal luminosities of the SXTs in quiescence and the observed neutron star crust cooling in SXT MXB 1659$-$29. We test different models of the equation of state and baryon superfluidity and take into account additional heat sources in the shallow layers of neutron-star crust (the shallow heating). We find that the observed quasi-stationary thermal luminosities of the SXTs can be equally well fitted using the traditional and thermodynamically consistent models, provided that the shallow heat diffusion into the core is taken into account. The observed crust cooling in MXB 1659$-$29 can also be fitted in the frames of both models, but the choice of the model affects the derived parameters responsible for the thermal conductivity in the crust and for the shallow heating.

Neutron stars are one of the most extreme objects in the universe, with densities that can exceed those of atomic nuclei and gravitational fields that are among the strongest known. Theoretical and observational research on neutron stars has revealed a wealth of information about their structural characteristics and physical properties. The structural characteristics of neutron stars are determined by the equations of state that describe the relationship between their density, pressure, and energy. These equations of state are still not well understood, and ongoing theoretical research aims to refine our understanding of the behavior of matter under these extreme conditions. Observational research on neutron stars, such as measurements of their masses and radii, can provide valuable constraints on the properties of the equation of state. The physical properties of neutron stars are also of great interest to researchers. Neutron stars have strong magnetic fields, which can produce observable effects such as pulsations and emission of X-rays and gamma rays. The surface temperature of neutron stars can also provide insight into their thermal properties, while observations of their gravitational fields can test predictions of Einstein's theory of general relativity. Observational research on neutron stars is carried out using a variety of techniques, including radio and X-ray telescopes, gravitational wave detectors, and optical telescopes. These observations are often combined with theoretical models to gain a more complete understanding of the properties of neutron stars.

The upcoming Nancy Grace Roman Space Telescope will carry out a wide-area survey in the near infrared. A key science objective is the measurement of cosmic structure via weak gravitational lensing. Roman data will be undersampled, which introduces new challenges in the measurement of source galaxy shapes; a potential solution is to use linear algebra-based coaddition techniques such as Imcom that combine multiple undersampled images to produce a single oversampled output mosaic with a desired "target" point spread function (PSF). We present here an initial application of Imcom to 0.64 square degrees of simulated Roman data, based on the Roman branch of the Legacy Survey of Space and Time (LSST) Dark Energy Science Collaboration (DESC) Data Challenge 2 (DC2) simulation. We show that Imcom runs successfully on simulated data that includes features such as plate scale distortions, chip gaps, detector defects, and cosmic ray masks. We simultaneously propagate grids of injected sources and simulated noise fields as well as the full simulation. We quantify the residual deviations of the PSF from the target (the "fidelity"), as well as noise properties of the output images; we discuss how the overall tiling pattern as well as Moir\'e patterns appear in the final fidelity and noise maps. We include appendices on interpolation algorithms and the interaction of undersampling with image processing operations that may be of broader applicability. The companion paper ("Paper II") explores the implications for weak lensing analyses.

One challenge for applying current weak lensing analysis tools to the Nancy Grace Roman Space Telescope is that individual images will be undersampled. Our companion paper presented an initial application of Imcom - an algorithm that builds an optimal mapping from input to output pixels to reconstruct a fully sampled combined image - on the Roman image simulations. In this paper, we measure the output noise power spectra, identify the sources of the major features in the power spectra, and show that simple analytic models that ignore sampling effects underestimate the power spectra of the coadded noise images. We compute the moments of both idealized injected stars and fully simulated stars in the coadded images, and their 1- and 2-point statistics. We show that the idealized injected stars have root-mean-square ellipticity errors (1 - 6) x 10-4 per component depending on the band; the correlation functions are >= 2 orders of magnitude below requirements, indicating that the image combination step itself is using a small fraction of the overall Roman 2nd moment error budget, although the 4th moments are larger and warrant further investigation. The stars in the simulated sky images, which include blending and chromaticity effects, have correlation functions near the requirement level (and below the requirement level in a wide-band image constructed by stacking all 4 filters). We evaluate the noise-induced biases in the ellipticities of injected stars, and explain the resulting trends with an analytical model. We conclude by enumerating the next steps in developing an image coaddition pipeline for Roman.

We present a pedagogical review of the halo model, a flexible framework that can describe the distribution of matter and its tracers on non-linear scales for both conventional and exotic cosmological models. We start with the premise that the complex structure of the cosmic web can be described by the sum of its individual components: dark matter, gas, and galaxies, all distributed within spherical haloes with a range of masses. The halo properties are specified through a series of simulation-calibrated ingredients including the halo mass function, non-linear halo bias and a dark matter density profile that can additionally account for the impact of baryon feedback. By incorporating a model of the galaxy halo occupation distribution, the properties of central and satellite galaxies, their non-linear bias and intrinsic alignment can be predicted. Through analytical calculations of spherical collapse in exotic cosmologies, the halo model also provides predictions for non-linear clustering in beyond-$\Lambda$CDM models. The halo model has been widely used to model observations of a variety of large-scale structure probes, most notably as the primary technique to model the underlying non-linear matter power spectrum. By documenting these varied and often distinct use cases, we seek to further coherent halo model analyses of future multi-tracer observables. This review is accompanied by the release of pyhalomodel: https://github.com/alexander-mead/pyhalomodel , flexible software to conduct a wide range of halo-model calculations.

The direction and magnitude of the dipole anisotropy of ultra-high-energy cosmic rays with energies above 8 EeV observed by the Pierre Auger Observatory indicate their extragalactic origin. The observed dipole on Earth does not necessarily need to correspond to the anisotropy of the extragalactic cosmic-ray flux due to the effects of propagation in the Galactic magnetic field. We estimate the size of these effects via numerical simulations using the CRPropa 3 package. The Jansson-Farrar and Terral-Ferri\`ere models of the Galactic magnetic field are used to propagate particles from the edge of the Galaxy to an observer on Earth. We identify allowed directions and amplitudes of the dipole outside the Galaxy that are compatible with the measured features of the dipole on Earth for various mass composition scenarios at the 68% and 95% confidence level.

Several MHD works and, in particular, the recent one by Medina-Torrejon et al. (2021) based on three-dimensional MHD simulations of relativistic jets, have evidenced that particle acceleration by magnetic reconnection driven by the turbulence in the flow occurs from the resistive up to the large injection scale of the turbulence. Particles experience Fermi-type acceleration up to ultra-high-energies, predominantly of the parallel velocity component to the local magnetic field, in the reconnection layers in all scales due to the ideal electric fields of the background fluctuations ($V\times B$, where $V$ and $B$ are the velocity and magnetic field of the fluctuations, respectively). In this work, we show MHD-particle-in-cell (MHD-PIC) simulations following the early stages of the particle acceleration in the relativistic jet which confirm these previous results, demonstrating the strong potential of magnetic reconnection driven by turbulence to accelerate relativistic particles to extreme energies in magnetically dominated flows. Our results also show that the dynamical time variations of the background magnetic fields do not influence the acceleration of the particles in this process.

There is no source for cosmic vorticity within the cold dark matter cosmology. However, vorticity has been observed in the universe, especially on the scales of clusters, filaments, galaxies, etc. Recent results from high-resolution general relativistic N-body simulation show that the vorticity power spectrum dominates over the power spectrum of the divergence of the peculiar velocity field on scales where the effective field theory of large-scale structure breaks down. Incidentally, this scale also corresponds to the scale where shell-crossing occurs. Several studies have suggested a link between shell crossing in the dark matter fluid and the vorticity generation in the universe, however, no clear proof of how it works exists yet. We describe for the first time how vorticity is generated in a universe such as ours with expanding and collapsing regions. We show how vorticity is generated at the boundary of the expanding and collapsing regions. Our result indicates that the amplitude of the generated vorticity is determined by the difference in gradients of the gravitational potential, pressure and the expansion rate of a one-parameter family of geodesics across the boundary. Additionally, we argue that the presence of vorticity in the matter fields implies a non-vanishing magnetic part of the Weyl tensor. This has implications for the generation of Maxwell's magnetic field and the dynamics of clusters. The impact on accelerated expansion of the universe and the existence of causal limit for massive particles are discussed

We present detailed modeling of periodic flaring events in the 6.7 GHz and 12.2 GHz methanol lines as well as the OH 1665 MHz and 1667 MHz transitions observed in the G9.62+0.20E star-forming region. Our analysis is performed within the framework of the one-dimensional Maxwell-Bloch equations, which intrinsically cover the complementary quasi-steady state maser and transient superradiance regimes. We find that the variations in flaring time-scales measured for the different species/transitions, and sometimes even for a single spectral line, are manifestations of and are best modeled with Dicke's superradiance, which naturally accounts for a modulation in the duration of flares through corresponding changes in the inversion pump. In particular, it can explain the peculiar behaviour observed for some features, such as the previously published result for the OH 1667 MHz transition at $v_\mathrm{lsr}=+1.7$ km s$^{-1}$ as well as the methanol 6.7 GHz line at $v_\mathrm{lsr}=-1.8$ km s$^{-1}$, through a partial quenching of the population inversion during flaring events.

The Galactic nova rate is intimately linked to our understanding of its chemical enrichment and progenitor channels of Type Ia supernovae. Yet past estimates have varied by more than an order of magnitude ($\approx10-300$ yr$^{-1}$) owing to limitations in both discovery methods as well as assumptions regarding the Galactic dust distribution and extragalactic stellar populations. Recent estimates utilizing synoptic near-infrared surveys have begun to provide a glimpse of a consensus ($\approx25-50$ yr$^{-1}$); however, a consistent estimate remains lacking. Here, we present the first all-sky search for Galactic novae using 8 years of data from the NEOWISE mid-infrared (MIR) survey. Operating at $3.4$ and $4.6$ $\mu$m where interstellar extinction is negligible, the 6-month cadence NEOWISE dataset offers unique sensitivity to discover slowly evolving novae across the entire Galaxy. Using a novel image subtraction pipeline together with systematic selection criteria, we identify a sample of 49 rapidly evolving MIR outbursts as candidate Galactic novae. While 27 of these sources are known novae, the remaining are previously missed nova candidates discovered in this work. The unknown events are spatially clustered along the densest and most heavily obscured regions of the Galaxy where previous novae are severely underrepresented. We use simulations of the NEOWISE survey strategy, the pipeline detection efficiency and our criteria to derive a Galactic nova rate of $47.9^{+3.1}_{-8.3}$ yr$^{-1}$. The discovery of these exceptionally bright (yet overlooked) nova candidates confirm emerging suggestions that optical surveys have been highly incomplete in searches for Galactic novae, highlighting the potential for MIR searches in revealing the demographics of Galactic stellar outbursts.

We present new estimates on the fraction of heavily X-ray obscured, Compton-thick (CT) active galactic nuclei (AGNs) out to a redshift of $z \leq$ 0.8. From a sample of 540 AGNs selected by mid-IR (MIR) properties in observed X-ray survey fields, we forward model the observed-to-intrinsic X-ray luminosity ratio ($R_{L_{\text{X}}}$) with a Markov chain Monte Carlo (MCMC) simulation to estimate the total fraction of CT AGNs ($f_{\text{CT}}$), many of which are missed in typical X-ray observations. We create model $N_{\text{H}}$ distributions and convert these to $R_{L_{\text{X}}}$ using a set of X-ray spectral models. We probe the posterior distribution of our models to infer the population of X-ray non-detected sources. From our simulation we estimate a CT fraction of $f_{\text{CT}}$ = $\text{0.555}^{+\text{0.037}}_{-\text{0.032}}$. We perform an X-ray stacking analysis for sources in Chandra X-ray Observatory fields and find that the expected soft (0.5-2 keV) and hard (2-7 keV) observed fluxes drawn from our model to be within 0.48 and 0.12 dex of our stacked fluxes, respectively. Our results suggests at least 50% of all MIR-selected AGNs, possibly more, are Compton-thick ($N_{\text{H}} \gtrsim$ 10$^{\text{24}}$ cm$^{-\text{2}}$), which is in excellent agreement with other recent work using independent methods. This work indicates that the total number of AGNs is higher than can be identified using X-ray observations alone, highlighting the importance of a multiwavelength approach. A high $f_{\text{CT}}$ also has implications for black hole (BH) accretion physics and supports models of BH and galaxy co-evolution that include periods of heavy obscuration.

In this work, we perform a phenomenological derivation of the first- and second-order relativistic hydrodynamics of dissipative fluids. To set the stage, we start with a review of the ideal relativistic hydrodynamics from energy-momentum and particle number conservation equations. We then go on to discuss the matching conditions to local thermodynamical equilibrium, symmetries of the energy-momentum tensor, decomposition of dissipative processes according to their Lorentz structure, and finally, the definition of the fluid velocity in the Landau and Eckart frames. With this preparatory work, we first formulate the first-order (Navier-Stokes) relativistic hydrodynamics from the entropy flow equation, keeping only the first-order gradients of thermodynamical forces. A generalized form of diffusion terms is found with a matrix of diffusion coefficients describing the relative diffusion between various flavors. The procedure of finding the dissipative terms is then extended to the second order to obtain the most general form of dissipative function for multiflavor systems up to the second order in dissipative fluxes. The dissipative function now includes in addition to the usual second-order transport coefficients of Israel-Stewart theory also second-order diffusion between different flavors. The relaxation-type equations of second-order hydrodynamics are found from the requirement of positivity of the dissipation function, which features the finite relaxation times of various dissipative processes that guarantee the causality and stability of the fluid dynamics. These equations contain a complete set of nonlinear terms in the thermodynamic gradients and dissipative fluxes arising from the entropy current, which are not present in the conventional Israel-Stewart theory.

We review applications of string theory to cosmology, from primordial times to the present-day accelerated expansion. Starting with a brief overview of cosmology and string compactifications, we discuss in detail moduli stabilisation, inflation in string theory, the impact of string theory on post-inflationary dynamics (reheating, moduli domination, kination), dark energy (the cosmological constant from a string landscape and models of quintessence) and various alternative scenarios (string/brane gases, the pre big-bang scenario, rolling tachyons, ekpyrotic/cyclic cosmologies, bubbles of nothing, S-brane and holographic cosmologies). The state of the art in string constructions is described in each topic and, where relevant, connections to swampland conjectures are made. The possibilities for novel particles and excitations (axions, moduli, cosmic strings, branes, solitons, oscillons and boson stars) are emphasised. Implications for the physics of the CMB, gravitational waves, dark matter and dark radiation are discussed along with potential observational signatures.

Without extensive maintenance, Dyson spheres inevitably disintegrate by asteroid impacts over billions of years. The resulting fragments would appear as anomalous interstellar objects, potentially sharing the unusual shape and motion of 1I/Oumuamua or the unusual material strength of the first two interstellar meteors, IM1 and IM2. If the Dyson sphere tiles are light sails, the number of fragments could exceed that of interstellar asteroids because of their reduced escape speed from the star and the increase in stellar luminosity during the red giant phase.

This work outlines a fast, high-precision time-domain solver for scalar, electromagnetic and gravitational perturbations on hyperboloidal foliations of Kerr space-times. Time-domain Teukolsky equation solvers have typically used explicit methods, which numerically violate Noether symmetries and are Courant-limited. These restrictions can limit the performance of explicit schemes when simulating long-time extreme mass ratio inspirals, expected to appear in LISA band for 2-5 years. We thus explore symmetric (exponential, Pad\'e or Hermite) integrators, which are unconditionally stable and known to preserve certain Noether symmetries and phase-space volume. For linear hyperbolic equations, these implicit integrators can be cast in explicit form, making them well-suited for long-time evolution of black hole perturbations. The 1+1 modal Teukolsky equation is discretized in space using polynomial collocation methods and reduced to a linear system of ordinary differential equations, coupled via mode-coupling arrays and discretized (matrix) differential operators. We use a matricization technique to cast the mode-coupled system in a form amenable to a method-of-lines framework, which simplifies numerical implementation and enables efficient parallelization on CPU and GPU architectures. We test our numerical code by studying late-time tails of Kerr spacetime perturbations in the sub-extremal and extremal cases.

The IBEX-Lo instrument on the Interstellar Boundary Explorer (IBEX) mission measures interstellar neutral (ISN) helium atoms. The detection of helium atoms is made through negative hydrogen (H$^-$) ions sputtered by the helium atoms from the IBEX-Lo conversion surface. The energy spectrum of ions sputtered by ISN helium atoms is broad and overlaps the four lowest IBEX-Lo electrostatic analyzer (ESA) steps. Consequently, the energy response function for helium atoms does not correspond to the nominal energy step transmission. Moreover, laboratory calibration is incomplete because it is difficult to produce narrow-energy neutral atom beams that are expected for ISN helium atoms. Here, we analyze the ISN helium observations in ESA steps 1-4 to derive the relative in-flight response of IBEX-Lo to helium atoms. We compare the ratios of the observed count rates as a function of the mean ISN helium atom energy estimated using the Warsaw Test Particle Model (WTPM). The WTPM uses a global heliosphere model to calculate charge exchange gains and losses to estimate the secondary ISN helium population. We find that the modeled mean energies of ISN helium atoms, unlike their modeled fluxes, are not very sensitive to the very local interstellar medium parameters. The obtained relative responses supplement the laboratory calibration and enable more detailed quantitative studies of the ISN helium signal. A similar procedure that we applied to the IBEX-Lo observations may be used to complement laboratory calibration of the next-generation IMAP-Lo instrument on the Interstellar Mapping and Acceleration Probe (IMAP) mission.

Magnetic reconnection and plasma turbulence are ubiquitous processes important for laboratory, space and astrophysical plasmas. Reconnection has been suggested to play an important role in the energetics and dynamics of turbulence by observations, simulations and theory for two decades. The fundamental properties of reconnection at kinetic scales, essential to understanding reconnection in turbulence, remain largely unknown at present. Here we present an application of the magnetic flux transport method that can accurately identify reconnection in turbulence to a three-dimensional simulation. Contrary to ideas that reconnection in turbulence would be patchy and unpredictable, highly extended reconnection X-lines, on the same order of magnitude as the system size, form at kinetic scales. Extended X-lines develop through bi-directional reconnection spreading, and satisfy critical balance characteristic of turbulence. This presents a new picture of fundamentally extended reconnection in kinetic-scale turbulence.

Neutron capture cross-section measurements are fundamental in the study of astrophysical phenomena, such as the slow neutron capture (s-) process of nucleosynthesis operating in red-giant and massive stars. However, neutron capture measurements via the time-of-flight (TOF) technique on key $s$-process nuclei are often challenging. Difficulties arise from the limited mass ($\sim$mg) available and the high sample-related background in the case of the unstable $s$-process branching points. Measurements on neutron magic nuclei, that act as $s$-process bottlenecks, are affected by low (n,$\gamma$) cross sections and a dominant neutron scattering background. Overcoming these experimental challenges requires the combination of facilities with high instantaneous flux, such as n\_TOF-EAR2, with detection systems with an enhanced detection sensitivity and high counting rate capabilities. This contribution reviews some of the latest detector developments in detection systems for (n,$\gamma$) measurements at n\_TOF, such as i-TED, an innovative detection system which exploits the Compton imaging technique to reduce the dominant neutron scattering background and s-TED, a highly segmented total energy detector intended for high flux facilities. The discussion will be illustrated with results of the first measurement of key the $s$-process branching-point reaction $^{79}$Se(n,$\gamma$).

Convolutional neural networks (CNNs) have seen extensive applications in scientific data analysis, including in neutrino telescopes. However, the data from these experiments present numerous challenges to CNNs, such as non-regular geometry, sparsity, and high dimensionality. Consequently, CNNs are highly inefficient on neutrino telescope data, and require significant pre-processing that results in information loss. We propose sparse submanifold convolutions (SSCNNs) as a solution to these issues and show that the SSCNN event reconstruction performance is comparable to or better than traditional and machine learning algorithms. Additionally, our SSCNN runs approximately 16 times faster than a traditional CNN on a GPU. As a result of this speedup, it is expected to be capable of handling the trigger-level event rate of IceCube-scale neutrino telescopes. These networks could be used to improve the first estimation of the neutrino energy and direction to seed more advanced reconstructions, or to provide this information to an alert-sending system to quickly follow-up interesting events.