Papers for Tuesday, Feb 06 2024

We present evidence of variable 1.3 millimeter emission from the 1-3 Myr, SpT G2-G5 class III YSO, HD~283572. HD~283572 was observed on 8 dates with the Submillimeter Array between 2021 December and 2023 May, a total on-source time of 10.2 hours, probing a range of timescales down to 5.2 seconds. Averaging all data obtained on 2022 Jan 17 shows a 4.4 mJy ($8.8\sigma$) point source detection with a negative spectral index ($\alpha{=}{-2.7}{\pm}1.2$), with peak emission rising to 13.8 mJy in one 3 minute span, and 25 mJy in one 29.7 second integration ($L_\nu=4.7\times10^{17}$ erg s$^{-1}$ Hz$^{-1}$). Combining our data for the other 7 dates shows no detection, with an rms noise of 0.24 mJy beam$^{-1}$. The stochastic millimeter enhancements on time frames of seconds--minutes--hours with negative spectral indices are most plausibly explained by synchrotron or gyro-synchrotron radiation from stellar activity. HD 283572's 1.3 mm light-curve has similarities with variable binaries, suggesting HD 283572's activity may have been triggered by interactions with an as-yet undetected companion. We additionally identify variability of HD 283572 at 10 cm, from VLASS data. This study highlights the challenges of interpreting faint mm emission from evolved YSOs that may host tenuous disks, and suggests that a more detailed temporal analysis of spatially unresolved data is generally warranted. The variability of class III stars may open up new ground for understanding the physics of flares in the context of terrestrial planet formation.

This is the second paper in a series presenting the results from a 500 $h^{-1}$Mpc large constrained hydro-dynamical simulation of the local Universe (SLOW). The initial conditions are based on peculiar velocities derived from the CosmicFlows-2 catalogue. The inclusion of galaxy formation treatment, allows to directly predict observable properties of the Intra-Cluster Medium (ICM) within galaxy clusters. Comparing the properties of observed galaxy clusters within the local Universe with the properties of their simulated counterparts, enables us to assess the effectiveness of the initial condition constraints in accurately replicating the non-linear properties of the largest, collapsed objects within the simulation. Based on the combination of several, publicly available surveys, we identified 45 local Universe galaxy clusters in SLOW, including the 13 most massive from the Planck SZ catalog and 70% of those with $M_{500} > 2\times 10^{14}$ M$_{\odot}$. We then derived the probability of the cross identification based on mass, X-ray luminosity, temperature and Compton-y by comparing it to a random selection. In relation to previous constrained simulations of the local volume, we found in SLOW a much larger amount of replicated galaxy clusters, where their simulation based mass prediction falls within the uncertainties of the observational mass estimates. Comparing the median observed and simulated masses of our cross identified sample allows to independently deduce a hydrostatic mass bias of $(1-b)\approx0.87$. The SLOW constrained simulation of the local Universe faithfully reproduces numerous fundamental characteristics of the galaxy clusters within our local neighbourhood, opening a new avenue for studying the formation and evolution of a large set of individual galaxy clusters as well as testing our understanding of physical processes governing the ICM.

Obscuration in high-redshift quasi-stellar objects (QSO) has a profound impact on our understanding of the evolution of supermassive black holes across the cosmic time. An accurate quantification of its relevance is therefore mandatory. We present a study aimed at evaluating the importance of obscuration in high redshift jetted QSO, i.e. those active nuclei characterized by the presence of powerful relativistic jets. We compare the observed number of radio detected QSO at different radio flux density limits with the value predicted by the beaming model on the basis of the number of oriented sources (blazars). Any significant deficit of radio-detected QSO compared to the predictions can be caused by the presence of obscuration along large angles from the jet direction. We apply this method to two sizable samples characterized by the same optical limit (mag=21) but significantly different radio density limits (30 mJy and 1 mJy respectively) and containing a total of 87 independent radio-loud 4<z<6.8 QSO, 31 of which classified as blazars. We find a general good agreement between the numbers predicted by the model and those actually observed, with only a marginal discrepancy at 0.5 mJy that could be caused by the lack of completeness of the sample. We conclude that we have no evidence of obscuration within angles 10-20deg from the relativistic jet direction. We also show how the ongoing deep wide-angle radio surveys will be instrumental to test the presence of obscuration at much larger angles, up to 30-35deg. We finally suggest that, depending on the actual fraction of obscured QSO, relativistic jets could be much more common at high redshifts compared to what is usually observed in the local Universe

We present the X-ray spectral analysis of two complementary sets of intermediate Seyfert galaxies (ISs). Analyzing X-ray data, we estimate the hydrogen abundance $N_H$ and test its connection with the [O III] luminosity acquired from optical observations. The results confirm the conclusions drawn in a previous study concerning the lack of a direct correlation between the obscuration measure ($N_H$) and the intrinsic characteristics of the active nuclei ([O III] luminosity). Instead, we validate the existence of a correlation between the Seyfert type and the $N_H$ parameter, employing a separation threshold of approximately 10$^{22}$ atoms cm$^{-2}$. Simultaneously, our findings align with prior research, corroborating the relationship between X-ray luminosity and the [O III] luminosity.

Protoclusters, the progenitors of galaxy clusters, trace large scale structures in the early Universe and are important to our understanding of structure formation and galaxy evolution. To date, only a handful of protoclusters have been identified in the Epoch of Reionization (EoR). As one of the rarest populations in the early Universe, distant quasars that host active supermassive black holes are thought to reside in the most massive dark matter halos at that cosmic epoch, and could thus potentially pinpoint some of the earliest protoclusters. In this letter, we report the discovery of a massive protocluster around a luminous quasar at $z=6.63$. This protocluster is anchored by the quasar, and includes three [CII] emitters at $z\sim6.63$, 12 spectroscopically confirmed Ly$\alpha$ emitters (LAEs) at $6.54<z\le6.64$, and a large number of narrow-band imaging selected LAE candidates at the same redshift. This structure has an overall overdensity of $\delta=3.3^{+1.1}_{-0.9}$ within $\sim35\times74$ cMpc$^2$ on the sky and an extreme overdensity of $\delta>30$ in its central region (i.e., $R\lesssim2$ cMpc). We estimate that this protocluster will collapse into a galaxy cluster with a mass of $6.9^{+1.2}_{-1.4}\times10^{15}~M_\odot$ at the current epoch, more massive than the most massive clusters known in the local Universe such as Coma. In the quasar vicinity, we discover a double-peaked LAE which implies that the quasar has a UV lifetime greater than 0.8 Myrs and has already ionized its surrounding intergalactic medium.

The inherent stochastic and nonlinear nature of the solar dynamo makes the strength of the solar cycles vary in a wide range, making it difficult to predict the strength of an upcoming solar cycle. Recently, our work has shown that by using the observed correlation of the polar field rise rate with the peak of polar field at cycle minimum and amplitude of following cycle, an early prediction can be made. In a follow-up study, we perform SFT simulations to explore the robustness of this correlation against variation of meridional flow speed, and against stochastic fluctuations of BMR tilt properties that give rise to anti-Joy and anti-Hale type anomalous BMRs. The results suggest that the observed correlation is a robust feature of the solar cycle and can be utilized for a reliable prediction of peak strength of a cycle at least 2 to 3 years earlier than the minimum.

A striking feature of the solar cycle is that at the beginning, sunspots appear around mid-latitudes, and over time the latitudes of emergences migrate towards the equator. The maximum level of activity varies from cycle to cycle. For strong cycles, the activity begins early and at higher latitudes with wider sunspot distributions than for weak cycles. The activity and the width of sunspot belts increase rapidly and begin to decline when the belts are still at high latitudes. However, in the late stages of the cycles, the level of activity, and properties of the butterfly wings all have the same statistical properties independent of the peak strength of the cycles. We have modelled these features using Babcock-Leighton type dynamo model and shown that the toroidal flux loss from the solar interior due to magnetic buoyancy is an essential nonlinearity that leads to all the cycles decline in the same way.

We describe photometry improvements in the La Silla--Quest RR Lyrae star (RRLS) survey that enable it to reach distances from the Sun ($d_{\odot}$) $\sim 140$ kpc. We report the results of surveying $\sim 300~ {\rm deg}^2$ of sky around the large, low-surface-brightness Crater II dwarf spheroidal galaxy. At $d_{\odot}$ $> 80$ kpc, we find a large overdensity of RRLS that extends beyond the traditional isophotal contours used for Crater II. The majority of these RRLS (34) have a linear distribution on the sky, extending over $15^{\circ}$, that runs through Crater II and is oriented along Crater II's proper motion vector. We hypothesize that this unlikely distribution traces extended tidal streams associated with Crater II. To test this, we search for other Crater II stellar populations that should be in the streams. Using Gaia proper motion data, we isolate $\approx$ 17 candidate stars outside of Crater II that are consistent with being luminous stars from the Crater II Red Giant Branch (RGB). Their spatial distribution is consistent with the RRLS one. The inferred streams are long, spanning a distance range $\sim 80 - 135$ kpc from the Galactic Centre. They are oriented at a relatively small angle relative to our line-of-sight ($\sim 25^{\circ}$), which means some stream stars are likely projected onto the main body of the galaxy. Comparing the numbers of RRLS and RGB candidate stars found in the streams to those in the main galaxy, we estimate Crater II has lost $\gtrsim 30\%$ of its stellar mass.

The disparate nature of thermal-nonthermal energy partition during flares, particularly during weak flares, is still an open issue. Following the Neupert effect, quantifying the relative yield of X-ray emission in different energy bands can enable inferring the underlying energy release mechanism. During September 20-25, 2021, the Solar Orbiter mission - being closer to the Sun ($\sim$0.6 AU) and having a moderate separation angle (<40$^{\circ}$) from the Sun-Earth line provided a unique opportunity to analyze multi-wavelength emission from $\sim$200 (mostly weak) flares, commonly observed by the Spectrometer Telescope for Imaging X-rays (STIX), STEREO-A, GOES, and SDO observatories. Associating the quotient (q$_{f}$) of hard X-ray fluence (12-20 keV) and soft X-ray flux (4-10 keV) with the peak SXR flux enabled us to identify strongly non-thermal flares. Multi-wavelength investigation of spectral and imaging mode observations of the 20 strongly non-thermal weak flares reveals an inverse relationship of q$_{f}$ with the emission measure (EM) (and density), and a positive relationship with the flare plasma temperature. This indicates that plasma in tenuous loops attains higher temperatures compared to that in the denser loops, in response to nonthermal energy deposition. This is in agreement with the plasma parameters of the coronal loops, as derived by applying the one-dimensional Palermo Harvard (PH) hydrodynamical code to the coronal loop plasma having different initial coronal loop base pressures when subjected to similar heating input. Our investigation, therefore, indicates that the plasma parameters of the flaring loop in the initial phase have a decisive role in thermal-nonthermal energy partitioning.

We present resolved images of IRAS 23077+6707 in $1.3\,\mathrm{mm}$/$225\,\mathrm{GHz}$ thermal dust and CO gas emission with the Submillimeter Array (SMA) and $grizy$ optical scattered light with the Panoramic Survey Telescope and Rapid Response System (Pan-STARRS). The Pan-STARRS data show a bipolar distribution of optically scattering dust that is characteristic for disks observed at high inclinations. Its radial extent in scattered light emission spans $\sim12''$, with two highly asymmetric filaments extending along the upper bounds of each nebula by $\sim9''$. The SMA data measure $1.3\,\mathrm{mm}$ continuum dust as well as $^{12}$CO, $^{13}$CO and C$^{18}$O $J$=2$-$1 line emission over $11''$-$16''$ extents, with the gas presenting the typical morphology of a disk in Keplerian rotation, in both position-velocity space and in each CO line spectra. IRAS 23077+6707 has no reported distance estimate, but if it is located in the Cepheus star-forming region ($180$-$800\,\mathrm{pc}$), it would have a radius spanning thousands of au. Taken together, we infer IRAS 23077+6707 to be a giant and gas-rich edge-on protoplanetary disk, which to our knowledge is the largest in extent so far discovered.

The equation of state dependence of neutron star's astrophysical features modeling is key to our understanding of dense matter. However, there exists a series of almost equation-of-state independent relations reported in the literature, called quasi-universal relations, that are used to determine neutron star radii and moments of inertia from X-ray and gravitational wave signals. Using sets of equations of state constrained by multi-messenger astronomy measurements and nuclear-physics theory, we discuss quasi-universal relations in the context of future gravitational-wave detectors Cosmic Explorer and Einstein Telescope, and X-ray detector STROBE-X. We focus on relations that involve the moment of inertia $I$, the tidal deformability $\Lambda$ and the compactness $C$: $C(\Lambda)$, $I(\Lambda)$ and $I(C)$. The quasi-universal fits and their associated errors are constructed with three different microphysics approaches which include state of the art nuclear physics theory and astrophysical constraints. Gravitational-wave and X-ray signals are simulated with the sensitivity of the next generation of detectors. Equation of state inference on those simulated signals is performed to assess if quasi-universal relations will offer a better precision on the extraction of neutron star's macroscopic parameters than equation of state dependent relations. We show that detections with the 3rd generation of gravitational wave detectors and the X-ray detector STROBE-X will be sensitive to the fit error marginalization technique. We also find that the sensitivity of those detectors will be sufficient that using full equation of state distributions will offer better precision on extracted parameters than quasi-universal relations.We also note that nuclear physics theory offers a more pronounced equation of state invariance of quasi-universal relations than current astrophysical constraints.

We present best-fit values of porosity -- and the corresponding effective thermal inertiae -- determined from three different depths in Europa's near-subsurface (~1-20 cm). The porosity of the upper ~20 cm of Europa's subsurface varies between 75-50% ($\Gamma_{eff}\approx50-140$ J m$^{-2}$ K$^{-1}$ s$^{-1/2}$) on the leading hemisphere and 50-40% ($\Gamma_{eff}\approx140-180$ J m$^{-2}$ K$^{-1}$ s$^{-1/2}$) on the trailing hemisphere. Residual maps produced by comparison with these models reveal thermally anomalous features that cannot be reproduced by globally homogeneous porosity models. These regions are compared to Europa's surface terrain and known compositional variations. We find that some instances of warm thermal anomalies are co-located with known geographical or compositional features on both the leading and trailing hemisphere; cool temperature anomalies are well correlated with surfaces previously observed to contain pure, crystalline water ice and the expansive rays of Pwyll crater. Anomalous regions correspond to locations with subsurface properties different from those of our best-fit models, such as potentially elevated thermal inertia, decreased emissivity, or more porous regolith. We also find that ALMA observations at ~3 mm sound below the thermal skin depth of Europa (~10-15 cm) for a range of porosity values, and thus do not exhibit features indicative of diurnal variability or residuals similar to other frequency bands. Future observations of Europa at higher angular resolution may reveal additional locations of variable subsurface thermophysical properties, while those at other wavelengths will inform our understanding of the regolith compaction length and the effects of external processes on the shallow subsurface.

The temperature structure of Titan's upper atmosphere exhibits large variability resulting from numerous spatially and temporally irregular external energy sources, seasonal changes, and the influence of molecular species produced via photochemistry. In particular, Titan's relatively abundant HCN is thought to provide substantial cooling to the upper atmosphere through rotational emission, balancing UV/EUV heating and thermal conduction. Here, we present the analysis of ALMA observations of Titan from 2012, 2014, 2015, and 2017, corresponding to planetocentric solar longitudes of ~34-89$^{\circ}$, including vertical HCN and temperature profiles retrieved from the lower mesosphere through the thermosphere (~350-1200 km; $3\times10^{-2}$-$2\times10^{-8}$ mbar). Throughout the atmosphere, temperature profiles differ by 10 to 30 K between observations approximately one Earth year apart, particularly from 600-900 km. We find evidence for a large imbalance in Titan's upper atmospheric energy budget between 2014 and 2015, where the mesospheric thermal structure changes significantly and marks the transition between a mesopause located at ~600 km ($2\times10^{-4}$ mbar) and at ~800 km ($3\times10^{-6}$ mbar). The retrieved HCN abundances vary dramatically during the 2012 to 2017 time period as well, showing close to 2 orders of magnitude difference in abundance at 1000 km. However, the change in HCN abundance does not appear to fully account for the variation in mesospheric temperatures over the $L_S\sim$34-89$^{\circ}$ period. These measurements provide additional insight into the variability of Titan's mesospheric composition and thermal structure following its 2009 vernal equinox, and motivate continued investigation of the origins of such rapid changes in Titan's atmosphere throughout its seasonal cycle.

Hyperactive comet activity typically becomes evident beyond the frost line (3 to 4 au) where it becomes too cold for water-ice to sublimate. If carbon monoxide (CO) and carbon dioxide (CO2) are the species that drive activity at sufficiently large distances, then detailed studies on the production rates of these species are extremely valuable to examine the formation of the solar system because these two species (beyond water) are next culpable for driving cometary activity. The NEOWISE reactivated mission operates at two imaging bandpasses, W1 and W2 at 3.4 and 4.6 microns, respectively, with the W2 channel being fully capable of detecting CO and CO2 at 4.67 and 4.23 microns in the same bandpass. It is extremely difficult to study CO2 from the ground due to contamination in Earth's atmosphere. We present our W1 and W2 photometry, dust measurements, and findings for comets C/2014 B1 (Schwartz), C/2017 K2 (Pan-STARRS), and C/2010 U3 (Boattini), hereafter, B1, K2, and U3, respectively. Our results assess CO and CO2 gas production rates observed by NEOWISE. We have determined: (1) comets B1 and K2 have CO2 and CO gas production rates of 1e27 and 1e29 molecules per second, respectively, if one assumes the excess emission is attributed to either all CO or all CO2; (2) B1 and K2 are considered hyperactive in that their measured AfRho dust production values are on the order of greater than or equal to 1e3 cm; and (3) the CO and CO2 production rates do not always follow the expected convention of increasing with decreased heliocentric distance, while B1 and K2 exhibit noticeable dust activity on their inbound leg orbits.

As next-generation imaging instruments and interferometers search for planets closer to their stars, they must contend with increasing orbital motion and longer integration times. These compounding effects make it difficult to detect faint planets but also present an opportunity. Increased orbital motion makes it possible to move the search for planets into the orbital domain, where direct images can be freely combined with the radial velocity and proper motion anomaly, even without a confirmed detection in any single epoch. In this paper, we present a fast and differentiable multimethod orbit-modeling and planet detection code called Octofitter. This code is designed to be highly modular and allows users to easily adjust priors, change parameterizations, and specify arbitrary function relations between the parameters of one or more planets. Octofitter further supplies tools for examining model outputs including prior and posterior predictive checks and simulation-based calibration. We demonstrate the capabilities of Octofitter on real and simulated data from different instruments and methods, including HD 91312, simulated JWST/NIRISS aperture masking interferometry observations, radial velocity curves, and grids of images from the Gemini Planet Imager. We show that Octofitter can reliably recover faint planets in long sequences of images with arbitrary orbital motion. This publicly available tool will enable the broad application of multiepoch and multimethod exoplanet detection, which could improve how future targeted ground- and space-based surveys are performed. Finally, its rapid convergence makes it a useful addition to the existing ecosystem of tools for modeling the orbits of directly imaged planets.

Exploring physics beyond General Relativity and the Standard Model of Particle Physics involves investigating spacetime variations in natural constants. This study employs an $H_2$-single of QSO 0347-383 observational spectrum to propose a unique approach for detecting potential changes in the proton-to-electron mass ratio. By comparing the ratio from observational and laboratory data in the Lyman-Alpha transition line, we derive a cosmological variation of $\Delta\mu / \mu = (0.120 \pm 0.144) \times 10^{-8}$ at $z_{abs}=3.025$. This approach not only advances fundamental physics understanding but also introduces innovative techniques for analyzing high-redshift QSO systems.

In 2016, an exposure meter was installed on the Lijiang Fiber-fed High-Resolution Spectrograph to monitor the coupling of starlight to the science fiber during observations. Based on it, we investigated a method to estimate the exposure flux of the CCD in real time by using the counts of the photomultiplier tubes (PMT) of the exposure meter, and developed a piece of software to optimize the control of the exposure time. First, by using flat-field lamp observations, we determined that there is a linear and proportional relationship between the total counts of the PMT and the exposure flux of the CCD. Second, using historical observations of different spectral types, the corresponding relational conversion factors were determined and obtained separately. Third, the method was validated using actual observation data, which showed that all values of the coefficient of determination were greater than 0.92. Finally, software was developed to display the counts of the PMT and the estimated exposure flux of the CCD in real-time during the observation, providing a visual reference for optimizing the exposure time control.

X-ray bursts are the thermonuclear runaway of a mixed H/He layer on the surface of a neutron star. Observations suggest that the burning begins locally and spreads across the surface of the star as a flame. Recent multidimensional work has looked in detail at pure He flames spreading across a neutron star. Here we report on progress in multidimensional modeling of mixed H/He flames and discuss the challenges.

We present the fourth data release (DR4) of the SkyMapper Southern Survey (SMSS), the last major step in our hemispheric survey with six optical filters: u, v, g, r, i, z. SMSS DR4 covers 26,000 sq.deg from over 400,000 images acquired by the 1.3m SkyMapper telescope between 2014-03 and 2021-09. The 6-band sky coverage extends from the South Celestial Pole to Dec = +16deg, with some images reaching Dec ~ +28deg. In contrast to previous DRs, we include all good-quality images from the facility taken during that time span, not only those explicitly taken for the public Survey. From the image dataset, we produce a catalogue of nearly 13 billion detections made from ~700 million unique astrophysical objects. The typical 10sigma depths for each field range between 18.5 and 20.5 mag, depending on the filter, but certain sky regions include longer exposures that reach as deep as 22 mag in some filters. As with previous SMSS catalogues, we have cross-matched with a host of other imaging and spectroscopic datasets to facilitate additional science outcomes. SMSS DR4 is now available to the worldwide astronomical community.

To understand the origin of dust in the circum-galactic medium (CGM), we develop a dust enrichment model. We describe each of the central galaxy and its CGM as a single zone, and consider the mass exchange between them through galactic inflows and outflows. We calculate the evolution of the gas, metal, and dust masses in the galaxy and the CGM. In the galaxy, we include stellar dust production and interstellar dust processing following our previous models. The dust in the galaxy is transported to the CGM via galactic outflows, and it is further processed by dust destruction (sputtering) in the CGM. We parameterize the time-scale or efficiency of each process and investigate the effect on the dust abundance in the CGM. We find that the resulting dust mass is sensitive to the dust destruction in the CGM, and the dust supply from galactic outflows, both of which directly regulate the dust abundance in the CGM. The inflow time-scale also affects the dust abundance in the CGM because it determines the gas mass evolution (thus, the star formation history) in the galaxy. The dust abundance in the CGM, however, is insensitive to stellar dust formation in the galaxy at later epochs because the dust production is dominated by dust growth in the interstellar medium. We also find that the resulting dust mass in the CGM is consistent with the value derived from a large sample of SDSS galaxies.

The color index $(J-K)_0$ of tip-red giant branch (TRGB) is used to study the metallicity distribution in the Large and Small Magellanic Cloud. With the most complete and pure sample of red member stars so far, the areas are divided into 154 and 70 bins for the LMC and SMC respectively with similar number of stars by the Voronoi binning. For each bin, the position of TRGB on the near-infrared color-magnitude diagram, specifically $(J-K)_0/K_0$, is determined by the Poison-Noise weighted method. Converting the color index of TRGB into metallicity, the metallicity gradients in the LMC and the SMC are obtained in four major directions. For the LMC, the gradient to the north is $-0.006 \pm 0.004$ dex kpc$^{-1}$, to the south $-0.010 \pm 0.005$ dex kpc$^{-1}$, to the east $-0.006 \pm 0.003$ dex kpc$^{-1}$, and to the west $-0.010 \pm 0.003$ dex kpc$^{-1}$. The farthest distance extends to 16 kpc. For the SMC, the gradients to the north, south, east, and west are $-0.017 \pm 0.031$ dex kpc$^{-1}$, $-0.016 \pm 0.007$ dex kpc$^{-1}$, $-0.003 \pm 0.002$ dex kpc$^{-1}$, and $-0.004 \pm 0.003$ dex kpc$^{-1}$, respectively. The farthest distance for the SMC extends to 27 kpc. The gradient is large from the center to 1 kpc.

Recently, large and homogeneous samples of cataclysmic variables (CVs) identified by the Sloan Digital Sky Survey (SDSS) were published. In these samples, the famous orbital period gap, which is a dearth of systems in the orbital period range ~2-3 hr and the defining feature of most evolutionary models for CVs, has been claimed not to be clearly present. If true, this finding would completely change our picture of CV evolution. In this Letter we focus on potential differences with respect to the orbital period gap between CVs in which the magnetic field of the white dwarf is strong enough to connect with that of the donor star, so-called polars, and non-polar CVs as the white dwarf magnetic field in polars has been predicted to reduce the strength of angular momentum loss through magnetic braking. We separated the SDSS I-IV sample of CVs into polars and non-polar systems and performed statistical tests to evaluate whether the period distributions are bimodal as predicted by the standard model for CV evolution or not. We confirm the existence of a period gap in the SDSS I-IV sample of non-polar CVs with >98 per cent confidence. The boundaries of the orbital period gap are 147 and 191 minutes, with the lower boundary being different to previously published values (129 min). The orbital period distribution of polars from SDSS I-IV is clearly different and does not show a similar period gap. The SDSS samples as well as previous samples of CVs are consistent with the standard theory of CV evolution. Magnetic braking does indeed seem get disrupted around the fully convective boundary, which causes a detached phase during CV evolution. In polars, the white dwarf magnetic field reduces the strength of magnetic braking and consequently the orbital period distribution of polars does not display an equally profound and extended period gap as non-polars.

Fan-spine magnetic structure, as a fundamental three-dimensional topology in magnetic reconnection theory, plays a crucial role in producing solar jets. However, how fan-spine configurations form in the solar atmosphere remains elusive. Using the Chinese H$\alpha$ Solar Explorer (CHASE) and the Solar Dynamics Observatory (SDO), we present a case study on the complete buildup of fan-spine topology driven by flux emergence and the subsequent jet production. Two fan-spine structures and the two associated null points are present. Variations in null-point heights and locations were tracked over time during flux emergence. The north fan-spine structure is found to be created through magnetic reconnection between the newly emerged flux and the background field. Gentle reconnection persistently occurs after formation of the north fan-spine structure, resulting in weak plasma outflows. Subsequently, as flux emergence and magnetic helicity injection continue, the formation and eruption of mini-filaments after reconnection at the quasi-separatrix layer between the two nulls trigger three homologous jets. The CHASE observations reveal that the circular flare ribbon, inner bright patch, and remote brightening all exhibit redshifted signatures during these jet ejections. This work unveils the key role of flux emergence in the formation of fan-spine topology, and highlights the importance of mini-filaments for subsequent jet production.

We have used the Ultraviolet Imaging Telescope (UVIT) aboard AstroSat to study star formation in a sample of nine dual nuclei galaxies with separations ~11 kpc, which is an expected outcome of galaxy mergers. To minimize the contribution of active galactic nuclei (AGN) emission, we have used mid-IR color cut-offs and masked the AGN-dominated nuclei. The UV continuum slope ($\beta$) is used to calculate the internal extinction (A$_V$) which ranges from 0.53 to 4.04 in the FUV band and 0.44 to 3.10 in the NUV band for the sample. We have detected $1-20$ star-forming clumps (SFCs) in our sample galaxies. The extinction-corrected total FUV star-formation rate (SFR) ranges from $\sim$0.35 to $\sim$32 M$_\odot$ yr$^{-1}$. Our analyses of A$_V$, specific SFR (sSFR) show that dual nuclei sources are associated with dusty, star-forming galaxies. The FUV$-$NUV color maps show redder color in the nuclear and disk regions while bluer color is observed in the outskirts of most galaxies which can be due to embedded dust or different stellar populations. We have found some signatures of possible stellar/AGN feedback like a ring of star formation, a redder ring around blue nuclei, etc. However, further observations are required to confirm this.

We present the results from spectral and timing study of the Z source GX 340+0 using AstroSat's SXT and LAXPC data. During the observation the source traced out the complete Z-track, allowing for the spectral evolution study of the Horizontal, Normal and Flaring branches (HB, NB and FB) as well as the hard and soft apexes (HA and SA). The spectra are better and more physically described by a blackbody component and a hot Comptonizing corona with a varying covering fraction, rather than one having a disc component. Along the track, the Comptonized flux (as well as the covering fraction) monotonically decreases. It is the blackbody component (both the temperature and radius) which varies non-monotonically and hence gives rise to the Z-track behaviour. Rapid timing study reveals a prominent Quasi-periodic Oscillation (QPO) at ~ 50 Hz at the HB, HA and upper NB, while a QPO at ~ 6 Hz is seen for the other branches. The fractional r.m.s of the QPOs increase with energy and exhibit soft lags in all branches except SA and FB.

Recently, combining high-contrast imaging and space astrometry we found that Jupiter-like (JL) planets are frequent in the beta Pic moving group (BPMG) around those stars where their orbit can be stable, prompting further analysis and discussion. We broaden our previous analysis to other young nearby associations to determine the frequency, mass, and separation of companions in general and JL in particular and their dependencies on the mass and age of the associations. We collected available data about companions including those revealed by visual observations, eclipses, spectroscopy, and astrometry. We determined search completeness and found that it is very high for stellar companions, while completeness corrections are still large for JL companions. Once these corrections are included, we found a high frequency of companions, both stellar (>0.52+/-0.03) and JL (0.57+/-0.11). The two populations are separated by a gap that corresponds to the brown dwarf desert. Within the population of massive companions, we found trends in frequency, separation, and mass ratios with stellar mass. Planetary companions pile up in the region just outside the ice line and we found them to be frequent once completeness was considered. The frequency of JL planets decreases with the overall mass and possibly the age of the association. We tentatively identify the two populations as due to disk fragmentation and core accretion, respectively. The distributions of stellar companions with a semi-major axis <1000 au is well reproduced by a simple model of formation by disk fragmentation. The observed trends with stellar mass can be explained by a shorter but much more intense phase of accretion onto the disk of massive stars and by a more steady and prolonged accretion on solar-type stars. Possible explanations for the trends in the population of JL planets with association mass and age are briefly discussed.

Using numerical simulations, we study the formation and dynamics of post-flare loops in a local region of the solar atmosphere. The MHD equations rule the post-flare structures' dynamic evolution, including space-dependent magnetic resistivity and highly anisotropic thermal conduction on a 2.5 D slice. We use an initial magnetic configuration consisting of a vertical current sheet, which helps trigger the magnetic reconnection process. Specifically, we study two scenarios, one with only resistivity and the second with resistivity plus thermal conduction. Numerical simulations show differences in the global morphology of the post-flare substructures in both cases. In particular, localized resistivity produces more substructure on the loops related to a Ritchmyer-Meshkov Instability (RMI). On the other hand, in the scenario with resistivity plus thermal conduction, the post-flare loops are smooth, and no apparent substructures develop. Besides, in the $z-$component of the current density for the Res+TC scenario, we observe the development of multiple small magnetic islands along the current sheet.

Recent astrophysical models predict that stellar-mass binary black holes (BBHs) could form and coalesce within a few gravitational radii of a supermassive black hole (SMBH). Detecting the gravitational waves (GWs) from such systems requires numerical tools which can track the dynamics of the binaries while capturing all the essential relativistic effects. This work develops upon our earlier study of a BBH moving along a circular orbit in the equatorial plane of a Kerr SMBH. Here we modify the numerical method to simulate a BBH falling toward the SMBH along a parabolic orbit of arbitrary inclination with respect to the equator. By tracking the evolution in a frame freely falling alongside the binary, we find that the eccentricity of the BBH is more easily excited than it is in the previous equatorial case, and that the cause is the asymmetry of the tidal tensor imposed on the binary when the binary moves out of the equatorial plane. Since the eccentricity reaches maximum around the same time that the BBH becomes the closest to the SMBH, multi-band GW bursts could be produced which are simultaneously detectable by the space- and ground-based detectors. We show that the effective spins of such GW events also undergo significant variation due to the rapid reorientation of the inner BBHs during their interaction with SMBHs. These results demonstrate the richness of three-body dynamics in the region of strong gravity, and highlight the necessity of building new numerical tools to simulate such systems.

A Be/X-ray binary system known as RX J0440.9+4431 (or LS V +44 17) is a potential member of the uncommon gamma-ray binary class. With an orbital period of 150 days, this system consists of a neutron star and a Be star companion. The MAXI observatory discovered an X-ray outburst from the source in December of 2022. Early in January, the outburst reached its peak, which was then followed by a decrease and a subsequent rebrightening. The X-ray flux exceeded 1 Crab in the 15-50 keV range at this second peak. AstroSat observations were conducted close to the peak of the second outburst, from January 11 to January 12, 2023. We report here the results of our search for 3-80 keV X-ray emission in the data of the AstroSat's LAXPC detector. It is found that the pulse period of the source is around 208 seconds. The source is found to be emitting more in the softer part of the X-ray energy range. The spectral characteristics can be described by employing a power-law model with an exponential cutoff, along with a blackbody component, interstellar absorption and an additional 6.4 keV iron fluorescence line.

Study of high-redshift radio galaxies (HzRGs) can shed light on the active galactic nuclei (AGNs) evolution in massive elliptical galaxies. The vast majority of observed high-redshift AGNs are quasars, and there are very few radio galaxies at redshifts $z>3$. We present the radio properties of 173 sources optically identified with radio galaxies at $z\geqslant1$ with flux densities $S_{1.4}\geqslant20$ mJy. Literature data were collected for compilation of broadband radio spectra, estimation of radio variability, radio luminosity, and radio loudness. Almost 60% of the galaxies have steep or ultra-steep radio spectra; 22% have flat, inverted, upturn, and complex spectral shapes, and 18% have peaked spectra (PS). The majority of the PS sources in the sample (20/31) are megahertz-peaked spectrum sources candidates, i.e. possibly very young and compact radio galaxies. The median values of the variability indices at 11 and 5 GHz are $V_{S_{11}}=0.14$ and $V_{S_{5}}=0.13$, which generally indicates a weak or moderate character of the long-term variability of the studied galaxies. The typical radio luminosity and radio loudness are $L_{5}=10^{43}$ - $10^{44}$ erg*s$^{-1}$ and $\log R=3$ - $4$ respectively. We have found less prominent features of the bright compact radio cores for our sample compared to high-redshift quasars at $z\geq3$. The variety of the obtained radio properties shows the different conditions for the formation of radio emission sources in galaxies.

The pulse widths of fast radio bursts are always broadened due to the scattering of the plasma medium through which the electromagnetic wave passes. The recorded pulse width will be further affected by the radio telescopes since the sampling time and the bandwidth cannot be infinitely small. In this study, we focus on the pulse widths of the 3287 bursts detected from FRB 20121102A as of October 2023. Various effects such as the scattering broadening, the redshift broadening and the instrumental broadening are examined. It is found that the instrumental broadening only contributes a fraction of $10^{-3}$--$10^{-1}$ to the observed pulse width. The scattering broadening is even smaller, which constitutes a tiny fraction of $10^{-5}$--$10^{-2}$ in the observed pulse width. After correcting for these broadenings, the intrinsic pulse width is derived for each burst. The maximum and minimum pulse widths at different frequencies are highlighted. Interestingly, both the mean value and the dispersion range of intrinsic pulse width are found to be inversely proportional to the square of the central frequency. The intrinsic widths of most bursts are in a narrow range of 1--10 ms, which leads to a quasi-linear correlation between the fluence and the peak flux.

Context. Propagating streamer blobs could be an important source of disturbances in the solar wind. Direct observations on formation of streamer blobs could be a proxy for understanding the formation of small-scale structures and disturbances in the solar wind. Aims. We aim to investigate how a streamer blob is formed before it is observed in the outer corona. Methods. Usingspecialcoordinated-observationsfromSOHO/LASCO,GOES/SUVIandSDO/AIA, we study the precursors of a streamer blob as seen in the corona below 2.0 solar radii (Rs). Results. We found that the streamer blob was formed due to the gradual merging of three clumps of brightenings initiated from the lower corona at about 1.8Rs, which is likely driven by expansion of the loop system at the base of the streamer. The acceleration of the blob starts from 1.9Rs or lower. It propagates along the south flank of the streamer where an expanding elongated brightening occurs coincidently. Conclusions. Our observations demonstrate that formation of a streamer blob is a complex process. We suggest that the expansion of the loop results in a pinching-off flux-rope-like blob at the loop apex below 2Rs. When the blob moves outward, it can be transferred across the overlying loops through interchange/component magnetic reconnection and then is released into the open field system. When the blob moves toward open field lines, interchange magnetic reconnections might also occur, and that can accelerate the plasma blob intermittently whilst allow it to transfer across the open field lines. Such dynamics in a streamer blob might further trigger small-scale disturbances in the solar wind such as switchbacks in the inner heliosphere.

Young Stellar Objects (YSO) are newly formed stars from molecular clouds. They stay close to where they were born and serve as good tracers to study gas and star formation. During cloud evolution, young massive stars can disrupt the surrounding gas through stellar feedback, changing the gas distribution. We study the distribution of the gas around a sample of YSO associations located at $d \lesssim 3.5 \;\rm kpc$ from the Sun by comparing the location and morphology between $^{12}$CO (J = 1$-$0) emission, Planck 870 $\mu$m maps and YSO associations. Based on the spatial distribution of the gas compared to that of the YSOs, we classify the YSO associations into six types: direct, close, bubble, complex, diffuse, and clean associations. The complex associations are large structures consisting of both gas-rich and gas-poor segments. We study the velocity dispersion-size relation toward different association types. From the ratio between different types, we estimate a feedback time of $\approx$ 1.7 Myr in the solar neighborhood. The sample sets a solid foundation to explore the relationship between interstellar medium evolution, star formation, and Galaxy structure.

In light of the non-perturbative resonance effects that may occur during inflation, we introduce a parametrization for the power spectrum of the stochastic gravitational wave background (SGWB) characterized by narrow-band amplification. We utilize the universal $\Omega_\text{GW}\propto k^3$ infrared limit, applicable to a wide array of gravitational wave sources, to devise a robust yet straightforward parameterization optimized for Markov Chain Monte Carlo (MCMC) analyses. This parameterization is demonstrated through select examples where its application is pertinent, and we discuss the advantages of this approach over traditional parametrizations for narrow-band scenarios. To evaluate the sensitivity of our proposed model parameters, we apply a mock likelihood based on the CMB-Stage4 data. Furthermore, we explicate the computational process for the mapping relationship between the foundational model parameters and our parameterized framework, using a two-field inflation model that resonantly amplifies gravitational waves (GWs) as an example.

In-ice radio-detection is a promising technique to discover and characterize ultra-high-energy (UHE) neutrinos, with energies above 100 PeV, adopted by present - ARA, ARIANNA, and RNO-G - and planned - IceCube-Gen2. So far, their ability to measure neutrino flavor has remained unexplored. We show and quantify how the neutrino flavor can be measured with in-ice radio detectors using two complementary detection channels. The first channel, sensitive to $\nu_e$, identifies them via their charged-current interactions, whose radio emission is elongated in time due to the Landau-Pomeranchuk-Migdal effect. The second channel, sensitive to $\nu_\mu$ and $\nu_\tau$, identifies events made up of multiple showers generated by the muons and taus they generate. We show this in state-of-the-art forecasts geared at IceCube-Gen2, for representative choices of the UHE neutrino flux. This newfound sensitivity could allow us to infer the UHE neutrino flavor composition at their sources - and thus the neutrino production mechanism - and to probe UHE neutrino physics.

Context - Around the year 2000, Triton's south pole experienced an extreme summer solstice that occurs every about 650 years, when the subsolar latitude reached about 50{\deg}. Bracketing this epoch, a few occultations probed Triton's atmosphere in 1989, 1995, 1997, 2008 and 2017. A recent ground-based stellar occultation observed on 6 October 2022 provides a new measurement of Triton's atmospheric pressure which is presented here. Aims- The goal is to constrain the Volatile Transport Models (VTMs) of Triton's atmosphere that is basically in vapor pressure equilibrium with the nitrogen ice at its surface. Methods - Fits to the occultation light curves yield Triton's atmospheric pressure at the reference radius 1400 km, from which the surface pressure is induced. Results - The fits provide a pressure p_1400= 1.211 +/- 0.039 microbar at radius 1400 km (47 km altitude), from which a surface pressure of p_surf= 14.54 +/- 0.47 microbar is induced (1-sigma error bars). To within error bars, this is identical to the pressure derived from the previous occultation of 5 October 2017, p_1400 = 1.18 +/- 0.03 microbar and p_surf= 14.1 +/- 0.4 microbar, respectively. Based on recent models of Triton's volatile cycles, the overall evolution over the last 30 years of the surface pressure is consistent with N2 condensation taking place in the northern hemisphere. However, models typically predict a steady decrease in surface pressure for the period 2005-2060, which is not confirmed by this observation. Complex surface-atmosphere interactions, such as ice albedo runaway and formation of local N2 frosts in the equatorial regions of Triton could explain the relatively constant pressure between 2017 and 2022.

Spectropolarimetric observations used to infer the solar magnetic fields are obtained with a limited spatial resolution. The effects of this limited resolution on the inference of the open flux over the observed region have not been extensively studied. We aim to characterize the biases that arise in the inference of the mean flux density by performing an end-to-end study that involves the generation of synthetic data, its interpretation (inversion), and a comparison of the results with the original model. We synthesized polarized spectra of the two magnetically sensitive lines of neutral iron around 630\,nm from a state-of-the-art numerical simulation of the solar photosphere. We then performed data degradation to simulate the effect of the telescope with a limited angular resolution and interpreted (inverted) the data using a Milne-Eddington spectropolarimetric inversion code. We then studied the dependence of the inferred parameters on the telescope resolution. The results show a significant decrease in the mean magnetic flux density -- related to the open flux observed at the disk center -- with decreasing telescope resolution. The original net magnetic field flux is fully resolved by a 1m telescope, but a 20\,cm aperture telescope yields a 30\% smaller value. Even in the fully resolved case, the result is still biased due to the corrugation of the photospheric surface. Even the spatially averaged quantities, such as the open magnetic flux in the observed region, are underestimated when the magnetic structures are unresolved. The reason for this is the presence of nonlinearities in the magnetic field inference process. This effect might have implications for the modeling of large-scale solar magnetic fields; for example, those corresponding to the coronal holes, or the polar magnetic fields, which are relevant to our understanding of the solar cycle.

We report on X-ray polarization measurements of the extra-galactic Crab-like PSR B0540-69 and its Pulsar Wind Nebula (PWN) in the Large Magellanic Cloud (LMC), using a ~850 ks Imaging X-ray Polarimetry Explorer (IXPE) exposure. The PWN is unresolved by IXPE. No statistically significant polarization is detected for the image-averaged data, giving a 99% confidence polarization upper limit (MDP99) of 5.3% in 2-8 keV energy range. However, a phase-resolved analysis detects polarization for both the nebula and pulsar in the 4-6 keV energy range. For the PWN defined as the off-pulse phases, the polarization degree (PD) of (24.5 ${\pm}$ 5.3)% and polarization angle (PA) of (78.1 ${\pm}$ 6.2){\deg} is detected at 4.6${\sigma}$ significance level, consistent with the PA observed in the optical band. In a single on-pulse window, a hint of polarization is measured at 3.8${\sigma}$ with polarization degree of (50.0 ${\pm}$ 13.1)% and polarization angle of (6.2 ${\pm}$ 7.4){\deg}. A 'simultaneous' PSR/PWN analysis finds two bins at the edges of the pulse exceeding 3${\sigma}$ PD significance, with PD of (68 ${\pm}$ 20)% and (62 ${\pm}$ 20)%; intervening bins at 2-3${\sigma}$ significance have lower PD, hinting at additional polarization structure.

The masses and radii of strongly magnetized anisotropic deformed white dwarf stars are investigated using the stellar structure equations in the parameterized $\gamma$-metric formalism. The Equation of State (EoS) of a completely degenerate relativistic electron gas in strong quantizing density-dependent magnetic field is developed. The fluid and field pressure anisotropy among the parallel and perpendicular components to the magnetic field is taken into consideration. This anisotropy in the EoS causes axisymmetric deformation of the star. We found stable solutions of deformed super-Chandrasekhar ultramassive white dwarfs. The masses of anisotropic magnetized white dwarfs at the same central density decrease monotonically with the increase in the strength of the central magnetic field, while the equatorial radii increase monotonically. This is in sharp contrast to the isotropic case where both the mass and radius increase monotonically. High magnetic field increases anisotropy and oblateness. We also see that the maximum mass and its corresponding equatorial radius both decrease as central magnetic field strength increases. We also notice that the maximum mass occurs at higher central density as the magnetic field increases. This shows that increasing magnetic field (hence increasing anisotropy) softens the EoS and makes the star more compact.

In Yukawa cosmology, a recent discovery revealed a relationship between baryonic matter and the dark sector. The relationship is described by the parameter $\alpha$ and the long-range interaction parameter $\lambda$ - an intrinsic property of the graviton. Applying the uncertainty relation to the graviton raises a compelling question: Is there a quantum mechanical limit to the measurement precision of the Hubble constant ($H_0$)? We argue that the uncertainty relation for the graviton wavelength $\lambda$ can be used to explain a running of $H_0$ with redshift. We show that the uncertainty in time has an inverse correlation with the value of the Hubble constant. That means that the measurement of the Hubble constant is intrinsically linked to length scales (redshift) and is connected to the uncertainty in time. On cosmological scales, we found that the uncertainty in time is related to the look-back time quantity. For measurements with a high redshift value, there is more uncertainty in time, which leads to a smaller value for the Hubble constant. Conversely, there is less uncertainty in time for local measurements with a smaller redshift value, resulting in a higher value for the Hubble constant. Therefore, due to the uncertainty relation, the Hubble tension is believed to arise from fundamental limitations inherent in cosmological measurements.

The young stellar source HH 30 is a textbook example of an ionic optical jet originating from a disk in an edge-on system shown by the HST. It has a remnant envelope in $^{12}$CO observed by ALMA. The optical jet is characterized by its narrow appearance, large line width at the base, and high temperature inferred from line diagnostics. Three featured structures can be identified, most evident in the transverse position--velocity diagrams: an extremely--high-velocity (EHV) wide-angle wind component with large spectral widths in the optical, a very--low-velocity (VLV) ambient surrounding medium seen in $^{12}$CO, and a low-velocity (LV) region traced by $^{12}$CO nested both in velocity and location between the primary wind and ambient environment. A layered cavity with multiple shells forms nested morphological and kinematic structures around the optical jet. The atomic gas originating from the innermost region of the disk attains a sufficient temperature and ionization to emit brightly in forbidden lines as an optical jet. The wide-angle portion expands, forming a low-density cavity. The filamentary $^{12}$CO encompassing the wind cavity is mixed and advected inward through the action of the magnetic interplay of the wide-angle wind with the molecular ambient medium. The magnetic interplay results in the layered shells penetrating deeply into the vast cavity of tenuous atomic wind material. The HH 30 system is an ideal manifestation of the unified wind model of \citet{Shang_2020,Shang_2023}, with clearly distinguishable atomic and molecular species mixed through the atomic lightly ionized magnetized wind and the surrounding cold molecular ambient material.

Weak gravitational lensing of distant galaxies provides a powerful probe of dark energy. The aim of this study is to investigate the application of convolutional neural networks (CNNs) to precision shear estimation. In particular, using a shallow CNN, we explore the impact of point spread function (PSF) misestimation and `galaxy population bias' (including `distribution bias' and `morphology bias'), focusing on the accuracy requirements of next generation surveys. We simulate a population of noisy disk and elliptical galaxies and adopt a PSF that is representative of a Euclid-like survey. We quantify the accuracy achieved by the CNN assuming a linear relationship between the estimated and true shears and measure the multiplicative ($m$) and additive ($c$) biases. We make use of an unconventional loss function to mitigate the effects of noise bias and measure $m$ and $c$ when we use either: (i) an incorrect galaxy ellipticity distribution or size-magnitude relation, or the wrong ratio of morphological types, to describe the population of galaxies (distribution bias); (ii) an incorrect galaxy light profile (morphology bias); or (iii) a PSF with size or ellipticity offset from its true value (PSF misestimation). We compare our results to the Euclid requirements on the knowledge of the PSF model shape and size. Finally, we outline further work to build on the promising potential of CNNs in precision shear estimation.

Periodic variables are always of great scientific interest in astrophysics. Thanks to the rapid advancement of modern large-scale time-domain surveys, the number of reported variable stars has experienced substantial growth for several decades, which significantly deepened our comprehension of stellar structure and binary evolution. The Tsinghua University-Ma Huateng Telescopes for Survey (TMTS) has started to monitor the LAMOST sky areas since 2020, with a cadence of 1 minute. During the period from 2020 to 2022, this survey has resulted in densely sampled light curves for ~ 30,000 variables of the maximum powers in the Lomb-Scargle periodogram above the 5sigma threshold. In this paper, we classified 11,638 variable stars into 6 main types using XGBoost and Random Forest classifiers with accuracies of 98.83% and 98.73%, respectively. Among them, 5301 (45.55%) variables are newly discovered, primarily consisting of Delta Scuti stars, demonstrating the capability of TMTS in searching for short-period variables. We cross-matched the catalogue with Gaia's second Data Release (DR2) and LAMOST's seventh Data Release (DR7) to obtain important physical parameters of the variables. We identified 5504 Delta Scuti stars (including 4876 typical Delta Scuti stars and 628 high-amplitude Delta Scuti stars), 5899 eclipsing binaries (including EA-, EB- and EW-type) and 226 candidates of RS Canum Venaticorum. Leveraging the metal abundance data provided by LAMOST and the Galactic latitude, we discovered 8 candidates of SX Phe stars within the class of "Delta Scuti stars". Moreover, with the help of Gaia color-magnitude diagram, we identified 9 ZZ ceti stars.

Lunar sulfides and oxides are a significant source of noble and base metals and will be vital for future human colonies' self-sustainability. Sulfide detection (pyrite and troilite) applies to many technological fields and use cases, for example, as a raw material source (available in situ on the Lunar surface) for new solar panel production methods. Ilmenite is the primary iron and titanium ore on the Moon and can provide helium-3 for nuclear fusion and oxygen for rocket fuel. The most important ore minerals have prominent absorption peaks in a narrow far-infrared (FIR) wavelength range of 20-40 $\mu$m, much stronger than the spectral features of other common minerals, including significant silicates, sulfates, and carbonates. Our simulations based on the linear mixing of pyrite with the silicates mentioned above indicated that areas containing at least 10-20% pyrite could be detected from the orbit in the FIR range. MIRORES, Multiplanetary far-IR ORE Spectrometer, proposed here, would operate with a resolution down to <5 m, enabling the detection of areas covered by 2-3 m$^2$ of pyrite (or ilmenite) on a surface of ~17 m$^2$ from an altitude of 50 km, creating possibilities for detecting large and local smaller orebodies along with their stockworks. The use of the Cassegrain optical system achieves this capability. MIRORES will measure radiation in eight narrow bands (0.3 $\mu$m in width) that can include up to five bands centered on the ore mineral absorption bands, for example, 24.3, 24.9, 27.6, 34.2, and 38.8 $\mu$m for pyrite, marcasite, chalcopyrite, ilmenite, and troilite, respectively. The instrument size is 32 $\times$ 32 $\times$ 42 cm, and the mass is <10 kg, which fits the standard microsatellite requirements.

We aim to provide observational signatures of the dust size evolution in the ISM, in particular exploring indicators of polycyclic aromatic hydrocarbon (PAH) mass fraction ($q_{PAH}$) defined as the mass fraction of PAHs relative to total dust grains. Additionally, we validate our dust evolution model by comparing the observational signatures from our simulations to those from observations. We model the evolution of grain size distribution of Milky Way-like and NGC 628-like galaxies representing star-forming galaxies with a hydrodynamic simulation code, GADGET4-OSAKA, which considers dust production and interstellar processing. Furthermore, we perform post-processing dust radiative transfer with SKIRT based on the simulations to predict the observational properties. We find that the intensity ratio between 8 um and 24 um correlates with $q_{PAH}$ and can be used as an indicator of PAH mass fraction. However, this ratio is influenced by the radiation field. As another indicator, we suggest the 8 um-to-total infrared intensity ratio ($\nu I_\nu(8 \mu m)/I$(TIR)), which tightly correlates with $q_{PAH}$. Furthermore, we explore the spatial evolution of $q_{PAH}$ in the simulated Milky Way-like galaxy using $\nu I_\nu(8 \mu m)/I$(TIR). We find that the spatially resolved $q_{PAH}$ increases with metallicity at lower metallicity (Z<0.2 Zsun) due to the interplay between accretion and shattering while it decreases with metallicity at higher metallicity (Z>0.2 Zsun) due to coagulation. Finally, we compare the above indicators in the NGC 628-like simulation with those observed in NGC 628. Consequently, our simulation underestimates the PAH mass fraction throughout the entire galaxy. This is probably because PAH is too efficiently lost by coagulation in the interstellar medium in our model, which suggests that the inhibition of coagulation of the PAHs is key to enhancing PAH formation.

Lithium's susceptibility to burning in stellar interiors makes it an invaluable tracer for delineating the evolutionary pathways of stars, offering insights into the processes governing their development. Observationally, the complex Li production and depletion mechanisms in stars manifest themselves as Li plateaus, and as Li-enhanced and Li-depleted regions of the HR diagram. The Li-dip represents a narrow range in effective temperature close to the main-sequence turn-off, where stars have slightly super-solar masses and strongly depleted Li. To study the modification of Li through stellar evolution, we measure 3D non-local thermodynamic equilibrium (NLTE) Li abundance for 581 149 stars released in GALAH DR3. We describe a novel method that fits the observed spectra using a combination of 3D NLTE Li line profiles with blending metal line strength that are optimized on a star-by-star basis. Furthermore, realistic errors are determined by a Monte Carlo nested sampling algorithm which samples the posterior distribution of the fitted spectral parameters. The method is validated by recovering parameters from a synthetic spectrum and comparing to 26 stars in the Hypatia catalogue. We find 228 613 Li detections, and 352 536 Li upper limits. Our abundance measurements are generally lower than GALAH DR3, with a mean difference of 0.23 dex. For the first time, we trace the evolution of Li-dip stars beyond the main sequence turn-off and up the subgiant branch. This is the first 3D NLTE analysis of Li applied to a large spectroscopic survey, and opens up a new era of precision analysis of abundances for large surveys.

We investigated the multiband emission from the pulsar binaries XSS J12270-4859, PSR J2039-5617, and PSR J2339-0533, which exhibit orbital modulation in the X-ray and gamma-ray bands. We constructed the sources' broadband spectral energy distributions and multiband orbital light curves by supplementing our X-ray measurements with published gamma-ray results, and we modeled the data using intra-binary shock (IBS) scenarios. While the X-ray data were well explained by synchrotron emission from electrons/positrons in the IBS, the gamma-ray data were difficult to explain with the IBS components alone. Therefore, we explored other scenarios that had been suggested for gamma-ray emission from pulsar binaries: (1) inverse-Compton emission in the upstream unshocked wind zone and (2) synchrotron radiation from electrons/positrons interacting with a kilogauss magnetic field of the companion. Scenario (1) requires that the bulk motion of the wind substantially decelerates to ~1000km/s before reaching the IBS for increased residence time, in which case formation of a strong shock is untenable, inconsistent with the X-ray phenomenology. Scenario (2) can explain the data if we assume the presence of electrons/positrons with a Lorentz factor of ~$10^8$ (~0.1 PeV) that pass through the IBS and tap a substantial portion of the pulsar voltage drop. These findings raise the possibility that the orbitally-modulating gamma-ray signals from pulsar binaries can provide insights into the flow structure and energy conversion within pulsar winds and particle acceleration nearing PeV energies in pulsars. These signals may also yield greater understanding of kilogauss magnetic fields potentially hosted by the low-mass stars in these systems.

Resolved stellar morphology of z>1 galaxies was inaccessible before JWST. This limitation, due to the impact of dust on rest-frame UV light, had withheld major observational conclusions required to understand the importance of clumps in galaxy evolution. Essentially independent of this issue, we use the rest-frame near-IR for a stellar-mass dependent clump detection method and determine reliable estimations of selection effects. We exploit publicly available JWST/NIRCam and HST/ACS imaging data from CEERS, to create a stellar-mass based picture of clumps in a mass-complete sample of 418 galaxies within a wide wavelength coverage of 0.5-4.6${\mu}$m and a redshift window of 1<z<2. We find that a near-IR detection gives access to a larger set of clumps within galaxies, with those also detected in UV making up only 28%. Whereas, 85% of the UV clumps are found to have a near-IR counterpart. These near-IR clumps closely follow the UVJ classification of their respective host galaxies, with these hosts mainly populating the star-forming regime besides a fraction of them (16%) that can be considered quiescent. The mass of the detected clumps are found to be within the range of $10^{7.5-9.5}\,\rm M_{\odot}$, therefore expected to drive gas into galaxy cores through tidal torques. However, there is likely contribution from blending of smaller unresolved structures. Furthermore, we observe a radial gradient of increasing clump mass towards the centre of galaxies. This trend could be an indication of clump migration, but accurate star-formation measurements would be required to confirm such a scenario.

DG Tau is a nearby T Tauri star associated with a collimated jet, a circumstellar disk and a streamer a few hundred au long. The streamer connects to the disk at $\sim$50 au from DG Tau. At this location SO emission is observed, likely due to the release of sulphur from dust grains caused by the shock of the impact of the accretion streamer onto the disk. We investigate the possibility that the DG Tau streamer was produced via cloudlet capture on the basis of hydrodynamic simulations, considering a cloudlet initiating infall at 600 au from DG Tau with low angular momentum so that the centrifugal force is smaller than the gravitational force, even at 50 au. The elongation of the cloudlet into a streamer is caused by the tidal force when its initial velocity is much less than the free-fall velocity. The elongated cloudlet reaches the disk and forms a high density gas clump. Our hydrodynamic model reproduces the morphology and line-of-sight velocity of CS ($5-4$) emission from the Northern streamer observed with ALMA. We discuss the conditions for forming a streamer based on the simulations. We also show that the streamer should perturb the disk after impact for several thousands of years.

Slow rotator galaxies are distinct amongst galaxy populations, with simulations suggesting that a mix of minor and major mergers are responsible for their formation. A promising path to resolve outstanding questions on the type of merger responsible, is by investigating deep imaging of massive galaxies for signs of potential merger remnants. We utilise deep imaging from the Subaru-Hyper Suprime Cam Wide data to search for tidal features in massive ($\log_{10}(M_*/M_{\odot}) > 10$) early-type galaxies (ETGs) in the SAMI Galaxy Survey. We perform a visual check for tidal features on images where the galaxy has been subtracted using a Multi-Gauss Expansion (MGE) model. We find that $31\pm 2$ percent of our sample show tidal features. When comparing galaxies with and without features, we find that the distributions in stellar mass, light-weighted mean stellar population age and H$\alpha$ equivalent width are significantly different, whereas spin ($\lambda_{R_e}$), ellipticity and bulge to total ratio have similar distributions. When splitting our sample in age, we find that galaxies below the median age (10.8 Gyr) show a correlation between the presence of shells and lower $\lambda_{R_e}$, as expected from simulations. We also find these younger galaxies which are classified as having "strong" shells have lower $\lambda_{R_e}$. However, simulations suggest that merger features become undetectable within $\sim 2-4$ Gyr post-merger. This implies that the relationship between tidal features and merger history disappears for galaxies with older stellar ages, i.e. those that are more likely to have merged long ago.

We use two-dimensional particle-in-cell simulations to investigate the generation and evolution of the magnetic field associated with the propagation of a jet for various initial conditions. We demonstrate that, in general, the magnetic field is initially grown by the Weibel and Mushroom instabilities. However, the field is saturated by the Alfv'en current limit. For initially non-magnetized plasma, we show that the growth of the magnetic field is delayed when the matter density of the jet environment is lower, which are in agreement with simple analytical predictions. We show that the higher Lorentz factor ($\gtrsim 2$) prevents rapid growth of the magnetic fields. When the initial field is troidal, the position of the magnetic filaments moves away from the jet as the field strength increases. The axial initial field helps the jet maintain its shape more effectively than the troidal initial field.

We calculate the phase space distribution function (DF) and the energy distribution of dark matter particles for a spherical halo in dynamical equilibrium assuming the Navarro-Frenk-White (NFW) density profile. Comparing the results with simulations of a wide range of haloes, we find that with appropriate matching, the energy distribution for a simulated halo can be well described by that derived from the best-fit NFW profile. Deviations occur at low energy when the NFW profile provides a poor fit for $r<0.05R_{vir}$, where $R_{vir}$ is the virial radius. The comparisons of DFs are similar to those of energy distributions, but the DF derived from the best-fit NFW profile has somewhat less accuracy because additional deviations are introduced through the density of energy states. We also compare the NFW fits to the simulated DFs and energy distributions with the DarkEXP fits of Hjorth & Williams (arXiv:1010.0265). We find that these fits have comparable accuracy in the region where both fit well, and that there is an approximate relation between the energy scale of the DarkEXP fits and the parameters of the NFW profile. The DarkEXP fits are better at low energy because they require the central gravitational potential as an input.

We present the results of the single-dish observations using the Korean VLBI Network to search for anomalous microwave emission (AME) in nearby galaxies. The targets were selected from MApping the dense moLecular gAs in the sTrongest stAr-formiNg Galaxies (MALATANG), a legacy survey project of the James Clerk Maxwell Telescope. The MALATANG galaxies are good representatives of local galaxies with enhanced nuclear activity associated with star formation and/or AGN, providing IR-bright galaxy samples; thus, they are good candidates for AME hosts. Combining with the ancillary data, we investigated the radio-IR spectral energy distribution (SED), while searching for the AME signals in five galaxies. The AME in NGC 2903 was well detected at a significant confidence level, whereas that in NGC 2146 and M82 was marginal. NGC 1068 and Arp 299 indicated no significant hints, and we provided the upper limits for the AME. The best-fit SED exhibited local peaks of the AME components at higher frequencies and with stronger peak fluxes than those in the previous studies. This suggested the origin of AME being denser environments such as molecular clouds or photodissociation regions rather than warm neutral/ionized medium as commonly suggested by previous studies. Further, our AME-detected targets were observed to exhibit higher specific star-formation rates than the other extragalactic AME hosts. Furthermore, AME favored starburst galaxies among our sample rather than AGN hosts. Consequently, this might imply that AGNs are excessively harsh environments for tiny dust to survive.

Multi-peaked supernovae with precursors, dramatic light-curve rebrightenings, and spectral transformation are rare, but are being discovered in increasing numbers by modern night-sky transient surveys like the Zwicky Transient Facility (ZTF). Here, we present the observations and analysis of SN 2023aew, which showed a dramatic increase in brightness following an initial luminous (-17.4 mag) and long (~100 days) unusual first peak (possibly precursor). SN 2023aew was classified as a Type IIb supernova during the first peak but changed its type to resemble a stripped-envelope supernova (SESN) after the marked rebrightening. We present comparisons of SN 2023aew's spectral evolution with SESN subtypes and argue that it is similar to SNe Ibc during its main peak. P-Cygni Balmer lines are present during the first peak, but vanish during the second peak's photospheric phase, before H$\alpha$ resurfaces again during the nebular phase. The nebular lines ([O I], [Ca II], Mg I], H$\alpha$) exhibit a double-peaked structure which hints towards a clumpy or non-spherical ejecta. We analyze the second peak in the light curve of SN 2023aew and find it to be broader than normal SESNe as well as requiring a very high $^{56}$Ni mass to power the peak luminosity. We discuss the possible origins of SN 2023aew including an eruption scenario where a part of the envelope is ejected during the first peak which also powers the second peak of the light curve through SN-CSM interaction.

The understanding of neutron star equation of state hinges on a comprehensive analysis of multi-messenger, multi-wavelength data. The recent scrutiny of PSR J0030+0451 data by NICER introduces complexities, unveiling a tension with another X-ray observation of the central compact object in HESS J1731-347, specifically concerning the mass-radius constraint of low-mass neutron stars. This tension persists when integrating NICER's updated data with LIGO/Virgo's gravitational-wave data from the GW170817 binary neutron star merger. Despite attempts to reconcile these disparate observations, the current combined data still can not distinguish different types of neutron stars -- whether they are pure neutron stars or hybrid stars. Bayesian inference indicates only modest changes in the posterior ranges of parameters related to the nuclear matter and deconfinement phase transition. This ongoing exploration underscores the intricate challenges in precisely characterizing neutron stars. It also points out that it is possible to probe the equation of state at different density regimes from future more accurate radii of neutron stars with various masses.

To understand physical properties of the interstellar medium (ISM) on various scales, we investigate it at parsec resolution on the kiloparsec scale. Here, we report on the sub-kpc scale Gas Density Histogram (GDH) of the Milky Way. The GDH is a density probability distribution function (PDF) of the gas volume density. Using this method, we are free from an identification of individual molecular clouds and their spatial structures. We use survey data of $^{12}$CO and $^{13}$CO ($J$=1-0) emission in the Galactic plane ($l = 10^{\circ}$-$50^{\circ}$) obtained as a part of the FOREST Unbiased Galactic plane Imaging survey with the Nobeyama 45m telescope (FUGIN). We make a GDH for every channel map of $2^{\circ} \times 2^{\circ}$ area, including the blank sky component, and without setting cloud boundaries. This is a different approach from previous works for molecular clouds. The GDH fits well to a single or double log-normal distribution, which we name the low-density log-normal (L-LN) and high-density log-normal (H-LN) components, respectively. The multi-log-normal components suggest that the L-LN and H-LN components originate from two different stages of structure formation in the ISM. Moreover, we find that both the volume ratios of H-LN components to total ($f_{\mathrm{H}}$) and the width of the L-LN along the gas density axis ($\sigma_{\rm{L}}$) show coherent structure in the Galactic-plane longitude-velocity diagram. It is possible that these GDH parameters are related to strong galactic shocks and other weak shocks in the Milky Way.

The data from the Gaia satellite led us to revise our conception of the Galaxy structure and history. Hitherto unknown components have been discovered and a deep re-thinking of what the Galactic halo is is in progress. We selected from the Gaia catalogue stars with extreme transverse velocities with respect to the Sun ($|V_T| > 500 $ and observed them with FORS2 at the ESO VLT, to classify them using both their chemical and dynamical properties. Two apparently young stars, identified in paper\,I, were observed with UVES. We derived abundances for Na, Mg, Ca, Ti, Mn, and Fe, analysing the spectra with while for Ba we used line profile fitting. We computed actions from parallaxes and kinematical data. The stars span the metallicity range $ Fe/H -0.5$ with $ Fe/H = -1.6$. Star GHS143 has a total speed of about 1440 which is almost three times faster than the local escape velocity of 522 strongly implying this star is unbound to the Galaxy. Remarkably, this star is not escaping from the Galaxy, but it is falling into it. Ten stars are apparently young with masses in excess of 1.3M. Their interpretation as evolved blue stragglers is doubtful. The existence of a young metal-poor population is possible. The two stars observed with UVES show no lithium, suggesting they are blue stragglers. We detected a metal-poor population, confined to the bulge, that we call SpiteF, and argue that it is the result of a recent accretion event. We detect 102 candidates of the Aurora population that should have formed prior to the formation of the disc. Our sample is non-homogeneous and mainly retrograde. The stars are metal poor, and 23<!PCT!> have Fe/H -2.0$. Our selection is efficient at finding very metal-poor stars, but it selects peculiar populations.

The Gaia satellite has already provided the astronomical community with three data releases, and the Radial Velocity Spectrometer (RVS) on board Gaia has provided the radial velocity for 33 million stars. When deriving the radial velocity from the RVS spectra, several stars are measured to have large values. To verify the credibility of these measurements, we selected some bright stars with the modulus of radial velocity in excess of 500\ to be observed with SOPHIE at OHP and UVES at VLT. This paper is devoted to investigating the chemical composition of the stars observed with UVES. We derived atmospheric parameters using Gaia photometry and parallaxes, and we performed a chemical analysis using the code. We find that the sample consists of metal-poor stars, although none have extremely low metallicities. The abundance patterns match what has been found in other samples of metal-poor stars selected irrespective of their radial velocities. We highlight the presence of three stars with low Cu and Zn abundances that are likely descendants of pair-instability supernovae. Two stars are apparently younger than 1\,Ga, and their masses exceed twice the turn-off mass of metal-poor populations. This makes it unlikely that they are blue stragglers because it would imply they formed from triple or multiple systems. We suggest instead that they are young metal-poor stars accreted from a dwarf galaxy. Finally, we find that the star RVS721 is associated with the Gjoll stream, which itself is associated with the Globular Cluster NGC\,3201.

The previous, high angular resolution 225 GHz ($\sim$1.3 mm) continuum observations on the transitional disk DM Tau have resolved an outer ring at 20-120 au radii that is weakly azimuthally asymmetric. We aimed to examine dust growth and filtration in the outer ring. We performed the $\sim$0$''$.06 ($\sim$8.7 au) resolution Karl G. Jansky Very Large Array (JVLA) 40-48 GHz ($\sim$7 mm; Q band) continuum observations and the complementary observations at lower frequencies. In addition, we analyzed the archival JVLA observations that were taken since 2010. Intriguingly, the Q band image resolved the azimuthally highly asymmetric, knotty dust emission sources close to the inner edge of the outer ring. Fitting the 8-700 GHz spectral energy distribution (SED) with two dust components indicates that the maximum grain size in these knotty dust emission sources is likely $\gtrsim$300 $\mu$m while it is $\lesssim$50 $\mu$m in the rest of the ring. These results may be explained by trapping of inward migrating grown dust close to the ring inner edge. The exact mechanism for developing the azimuthal asymmetry has not yet been identified, which may be due to planet-disk interaction that might also be responsible for the creation of the dust cavity and pressure bump, or the fluid instabilities and vortex formation due to shear motions. Finally, we remark that the asymmetries in DM Tau are hard to diagnose from the $\gtrsim$225 GHz observations owing to a high optical depth at the ring. In other words, the apparent symmetric or asymmetric morphology of the transitional disks may be related to the optical depths of those disks at the observing frequency.

We present photometric and spectroscopic data sets for SN 2020pvb, a Type IIn-P supernova (SN) similar to SNe 1994W, 2005cl, 2009kn and 2011ht, with a precursor outburst detected (PS1 w-band ~ -13.8 mag) around four months before the B-band maximum light. SN 2020pvb presents a relatively bright light curve peaking at M_B = -17.95 +- 0.30 mag and a plateau lasting at least 40 days before it went in solar conjunction. After this, the object is no longer visible at phases > 150 days above -12.5 mag in the B-band, suggesting that the SN 2020pvb ejecta interacts with a dense spatially confined circumstellar envelope. SN 2020pvb shows in its spectra strong Balmer lines and a forest of FeII lines with narrow P Cygni profiles. Using archival images from the Hubble Space Telescope, we constrain the progenitor of SN 2020pvb to have a luminosity of log(L/L_sun) <= 5.4, ruling out any single star progenitor over 50 M_sun. All in all, SN 2020pvb is a Type IIn-P whose progenitor star had an outburst ~ 0.5 yr before the final explosion, the material lost during this outburst is probably playing a role in shaping the physical properties of the supernova.

The HERMES Technologic and Scientific Pathfinder project is a constellation of six CubeSats aiming to observe transient high-energy events such as the Gamma Ray Bursts (GRBs). HERMES will be the first space telescope to include a siswich detector, able to perform spectroscopy in the 2 keV to 2 MeV energy band. The particular siswich architecture, which combines a solid-state Silicon Drift Detector and a scintillator crystal, requires specific calibration procedures that have not been yet standardized in a pipeline. We present in this paper the HERMES calibration pipeline, mescal, intended for raw HERMES data energy calibration and formatting. The software is designed to deal with the particularities of the siswich architecture and to minimize user interaction, including also an automated calibration line identification procedure, and an independent calibration of each detector pixel, in its two different operating modes. The mescal pipeline can set the basis for similar applications in future siswich telescopes.

The emergence of the disc in our Galaxy and the relation of the thick and thin disc formation and evolution is still a matter of debate. The chemo-dynamical characterization of disc stars is key to resolve this question, in particular at parameter regimes where both disc components overlap, such as the region around [Fe/H] $\sim$ $-0.7$ corresponding to the thin disc metal-poor end. In this paper we re-assess the recent detection of a metal-poor extension of stars moving with thin-disc-like rotational velocities between -2 < [Fe/H] < -0.7 that was made based on metallicity estimates obtained from photometric data and their rotational velocity distribution. We explore the chemo-dynamical properties of metal-poor stars within the recent Gaia third data release (DR3), which includes the first catalogue of metallicity estimates from the Radial Velocity Spectrometer (RVS) experiment. We complement them with the two largest high-resolution ($\lambda/d\lambda$ > 20,000) spectroscopic surveys available, the GALAH DR3 and the APOGEE DR17. We confirm that there are high angular-momentum stars moving in thin-disc-like orbits, i.e., with high angular momentum $\rm L_{z}/J_{tot}$ > 0.95, and close to the Galactic plane, $\rm |Z_{max}|$ < 750 pc, reaching metallicity values down to [Fe/H] $\sim-1.5$. We also find tentative evidence of stars moving on such orbits at lower metallicities, down to [Fe/H] $\sim -2.5$, although in smaller numbers. Based on their chemical trends the fast rotators with [Fe/H] < -1 would have formed in a medium less chemically evolved than the bulk of the thick disc. Fast rotators with chemical abundances typical of the thin disc appear at metallicities between -1 < [Fe/H] < -0.7.

We investigated the ionised and atomic gas kinematics and excitation state in the central region of ongoing star formation of the nearby low-metallicity dwarf galaxy Sextans B. The analysis is based on the new observations performed in Ha emission line with high resolution ($R \sim 16000$) scanning Fabry-Perot interferometer at the 6-m BTA SAO RAS telescope, and on the long-slit spectral observations at the 9.2-m SALT and 2.5-m CMO SAI MSU telescopes. Strong non-circular gas motions detected in the studied regions probably resulted from the off-plane gas motions and impact of stellar feedback. We identified six regions of elevated Ha velocity dispersion, five of which exhibit asymmetric or two-component Ha line profiles. Three of these regions are young ($<1.1$ Myr) expanding ($V_{\rm exp} \sim 25-50\ {\rm km\ s^{-1}}$) superbubbles. We argue that at least three regions in the galaxy could be supernova remnants. We conclude that supernovae feedback is the dominant source of energy for superbubbles in Sextans B, which is expected for such a low metallicity, although we cannot rule out a strong impact of pre-supernova feedback for one superbubble.

In Chen et al., 2023 (C23; arXiv:2304.09898), we introduced a new method to directly measure temperature fluctuations and applied it to a nearby dwarf galaxy, Mrk 71, finding a temperature fluctuation parameter $t^2 = 0.008\pm 0.043$. This result is lower by $\sim 2\sigma$ than the value required to explain the abundance discrepancy (AD) in this object. In the Matters Arising article submitted by Mendez-Delgado et al. (arXiv:2310.01197), the authors claim that using the same data presented in C23 in a different way, it is possible to conclude that the measurements are consistent with a larger $t^2 \simeq 0.1$ inferred indirectly from recombination lines (RLs). However, this requires a higher density such that the infrared [O III] 52 $\mu$m and [O III] 88 $\mu$m lines -- which form the basis of the direct measurement method -- are mutually inconsistent. Moreover, to reach agreement between the direct $t^2$ measurement and the larger $t^2$ value inferred from RLs requires systematically varying four parameters by $\sim 1\sigma$ from their best-determined values, which collectively amount to a $\sim2\sigma$ difference, consistent with the significance ($\sim 2 \sigma$) originally reported in C23. Therefore, we conclude that the results of C23 hold, and that the combined optical and infrared [O III] data disfavour $t^2 \simeq 0.1$ at the $\approx2\sigma$ level in Mrk 71. Future work is nonetheless warranted to better understand the AD associated with both optical and infrared emission line analysis.

Relativistic pair beams produced in the cosmic voids by TeV gamma rays from blazars are expected to produce a detectable GeV-scale cascade that is missing in the observations. The suppression of this secondary cascade implies either the deflection of the pair beam by intergalactic magnetic fields or, alternatively, an energy loss of the beam due to the beam-plasma instability. Here, we study how the beam-plasma instability feeds back on the beam, using a realistic two-dimensional beam distribution. We find that the instability broadens the beam opening angles significantly without any significant energy loss, thus confirming a recent feedback study on a simplified one-dimensional beam distribution. However, narrowing diffusion feedback of the beam particles with Lorentz factors less than $10^6$ might become relevant even though initially it is negligible. Finally, when considering the continuous creation of TeV pairs, we find that the beam distribution and the wave spectrum reach a new quasi-steady state, in which the scattering of beam particles persists and the beam opening angle may increase by a factor of hundreds. Understanding the implications on the GeV cascade emission requires accounting for inverse-Compton cooling.

The formation of binary stars is highly influenced by magnetic fields, which play a crucial role in transporting angular momentum. We conducted three-dimensional numerical simulations of binary star accretion via a circumbinary disk, taking into account a magnetic field perpendicular to the disk and an infalling envelope. Our simulations reproduce the following phenomena: (1) spiral arms associated with circumstellar disks, (2) turbulence in the circumbinary disk, induced by magneto-rotational instability (MRI), (3) a fast outflow launched from each circumstellar disk, and (4) a slow outflow from the circumbinary disk. The binary models exhibit a higher $\alpha$-parameter than the corresponding single star models, indicating that the binary stars enhance MRI turbulence. Moreover, an infalling envelope also enhance the turbulence, leading to a high $\alpha$-parameter. While the spiral arms promotes radial flow, causing transfer of mass and angular momentum within the circumbinary disk, the MRI turbulence and outflows are main drivers of angular momentum transfer to reduce the specific angular momentum of the system.

Observations of neutron star mergers have the potential to unveil detailed physics of matter and gravity in regimes inaccessible by other experiments. Quantitative comparisons to theory and parameter estimation require nonlinear numerical simulations. However, the detailed physics of energy and momentum transfer between different scales, and the formation and interaction of small scale structures, which can be probed by detectors, are not captured by current simulations. This is where turbulence enters neutron star modelling. This review will outline the theory and current status of turbulence modelling for relativistic neutron star merger simulations.

This paper presents a comprehensive spectroscopic analysis of the nova RS Ophiuchi during its quiescent stage, spanning a duration of approximately 13 years. The spectra exhibit prominent low-ionization emission features, including hydrogen, helium, iron, and TiO absorption features originating from the cool secondary component. The CLOUDY photoionization code is employed to model these spectra, allowing us to estimate various physical parameters such as temperature, luminosity, and hydrogen density, along with elemental abundances and accretion rate. The central ionizing sources exhibit temperatures in the range of $1.05 - 1.8~\times 10^4$ K and luminosities between $0.1 - 7.9~\times 10^{30}$ \ergs. Notably, \ion{He}{} displays an overabundance from 2008 to 2016, returning to solar values by 2020, while \ion{Fe}{} appears subsolar from 2008 to 2014 but becomes overabundant from 2006 onward. The mean accretion rate, as calculated from the model, is approximately $1.254 \times 10^{-8} M_{\odot}$ yr$^{-1}$. About 47\% of the critical mass was accreted after April, 2020 ($\sim$15 months before the 2021 outburst), and approximately 88\% of the critical mass was accreted after July 20, 2018. This non-uniform accretion rate suggests a more rapid approach towards reaching the critical mass in the final years, possibly attributed to the heightened gravitational pull resulting from previously accreted matter, influencing the accretion dynamics as the system approaches the critical mass limit.

Scenarios such as the QCD axion with the Peccei-Quinn symmetry broken after inflation predict an enhanced matter power spectrum on sub-parsec scales. These theories lead to the formation of dense dark matter structures known as minihalos, which provide insights into early Universe dynamics and have implications for direct detection experiments. We examine the mass loss of minihalos during stellar encounters, building on previous studies that derived formulas for mass loss and performed N-body simulations. We propose a new formula for the mass loss that accounts for changes in the minihalo profile after disruption by a passing star. We also investigate the mass loss for multiple stellar encounters. We demonstrate that accurately assessing the mass loss in minihalos due to multiple stellar encounters necessitates considering the alterations in the minihalo's binding energy after each encounter, as overlooking this aspect results in a substantial underestimation of the mass loss.

We present geomagnetic storms (GSs) selected from three solar cycles, spanning the years 1995 to 2022. We studied the development of the main phase of storms within disturbance storm time (Dst) amplitudes ranging from Dst =-64 nT to Dst=- 422 nT. In order to determine the solar wind (SW) parameters that mainly influence the main phase development of a GS, which can best describe the SW-magnetosphere coupling, we divided our selected GSs into four groups based on main phase duration. Superposed epoch analysis was performed on the selected geomagnetic indices, SW plasma and field parameters, and their derivatives separately for each group. To that end, the dynamics of GS main phase development is mainly guided by interplanetary driver magnetic field southward component, (-Bz). It has been determined that there is a temporal difference between the peak values of Bz and Dst. As a result, Dst is delayed from Bz by 1-4 hours, which is crucial for space weather forecasting. The peak of Dst has a direct relationship with the amplitude of storm sudden commencement (SSC) and an inverse relationship with the duration of SSC. The inter-relationship between the peaks of the three indices (Dst, AE, and ap) during GS, is also obtained. Dst is found to be more closely related to ap than AE. To determine the best fit SW parameter to the geomagnetic activity indices, we used a linear correlation between the peak values of individual geomagnetic indices and SW plasma and field parameters and their derivatives. An electric field related function involving speed and IMF (v$^{4/3}$Bz) when coupled with a viscous term ($\rho^{1/2}$) correlates very well with the intensity of the GS (Dst$ _{min} $ or $ \Delta$Dst) and the magnitude of (ap$_{max}$) and (AE$_{max}$) during storms. However, a related function (v$^{4/3}$B$\rho^{1/2}$) represents slightly better the peak of AE$_{max}$ during the storms.

Kinematically-persistent planes of satellites (KPPs) are fixed sets of satellites co-orbiting around their host galaxy, whose orbital poles are conserved and clustered across long cosmic time intervals. They play the role of 'skeletons', ensuring the long-term durability of positional planes. We explore the physical processes behind their formation in terms of the dynamics of the local Cosmic Web (CW), characterized via the so-called Lagrangian Volumes (LVs) built up around two zoom-in, cosmological hydro-simulations of MW-mass disk galaxy + satellites systems, where three KPPs have been identified. By analyzing the LVs deformations in terms of the reduced Tensor of Inertia (TOI), we find an outstanding alignment between the LV principal directions and KPP satellites' orbital poles. The most compressive local mass flows (along the $\hat{e}_3$ eigenvector) are strong at early times, feeding the so-called $\hat{e}_3$-structure, while the smallest TOI axis rapidly decreases. The $\hat{e}_3$-structure collapse marks the end of this regime and is the timescale for the establishment of satellite orbital pole clustering when the Universe is $\lesssim$ 4 Gyr old. KPP proto-satellites aligned with $\hat{e}_3$ are those whose orbital poles are either aligned from early times, or have been successfully bent at $\hat{e}_3$-structure collapse. KPP satellites associated to $\hat{e}_1$ tend to have early trajectories already parallel to $\hat{e}_3$. We show that KPPs can arise as a result of the $\Lambda$CDM-predicted large-scale dynamics acting on particular sets of proto-satellites, the same dynamics that shape the local CW environment.

We propose a novel cosmological framework within the $f(R,T)$ type modified gravity theory, incorporating a non-minimally coupled with the higher order of the Ricci scalar ($R$) as well as the trace of the energy-momentum tensor ($T$). Therefore, our well-motivated chosen $f(R,T)$ expression is $ R + R^m + 2 \lambda T^n$, where $\lambda$, $m$, and $n$ are arbitrary constants. Taking a constant jerk parameter ($j$), we derive expressions for the deceleration parameter ($q$) and the Hubble parameter ($H$) as functions of the redshift $z$. We constrained our model with the recent Observational Hubble Dataset (OHD), $Pantheon$, and $ Pantheon $ + OHD datasets by using the analysis of Markov Chain Monte Carlo (MCMC). Our model shows early deceleration followed by late-time acceleration, with the transition occurring in the redshift range $1.10 \leq z_{tr} \leq 1.15$. Our findings suggest that this higher-order model of $f(R,T)$ gravity theory can efficiently provide a dark energy model for addressing the current scenario of cosmic acceleration.

The LISA gravitational wave observatory is due to launch in 2035. Researchers are currently working out the data analysis picture, including determining requirements of, and prospects for, parameterised gravitational wave models used in Bayesian inference, of which early stage inspirals of low and intermediate mass black hole binaries constitute a significant and important part. With datasets consisting of up to hundreds of millions of datapoints, their likelihood functions can be extremely expensive to compute. We present an approximation procedure for accurately reproducing the likelihood function (from time-domain data) for such signals in simulated environments, significantly reducing this cost. The method is simply to discard most of the datapoints and define a new inner product operator in terms of the remaining data, which closely resembles the original inner product on the model parameter space. Highly accurate reproductions of the likelihood can be achieved using just a few hundred to a few thousand datapoints, which can reduce parameter estimation (PE) time by factors of up to tens of thousands (compared to estimated frequency domain PE time), however, the end user must be responsible for ensuring convergence. The time-domain implementation is particularly useful for modelling waveform modifications arising from non-trivial astrophysical environments. The contents of this article provide the theoretical basis of the software package dolfen.

Primordial black holes (PBHs), if captured by neutron stars (NSs), would emit a characteristic gravitational wave (GW) signal as they orbit inside the host star. We identify a specific and qualitatively new feature of these signals, namely quasi-periodic beats caused by the precession of noncircular PBH orbits. We demonstrate numerically and analytically that the beat frequency depends rather sensitively on the NS structure, so that hypothetical future observations with next-generation GW detectors would provide valuable constraints on the nuclear equation of state.

We model and constrain the unique indirect detection signature produced by dark matter particles that annihilate through a $U(1)$ gauge symmetry into dark photons that subsequently decay into three-photon final states. We focus on scenarios where the dark photon is long-lived, and show that $\gamma$-ray probes of celestial objects can set strong constraints on the dark matter/baryon scattering cross section that in many cases surpass the power of current direct detection constraints, and in some cases even peer into the neutrino fog.

Clement and Noiucer submitted a note {\tt arXiv:2401.16008 [gr-qc]} replying to our criticism {\tt arXiv:2401.10479 [gr-qc]} of their previous submission. We reply to the contents of this note and remark that these authors have not addressed our arguments. This will be our last response to them. Readers are advised to look at all material and judge by themselves

In this work we shall prove that the tensor spectral index of the primordial tensor perturbations for GW170817-compatible Einstein-Gauss-Bonnet theories, takes the approximate simplified form $n_{\mathcal{T}}\simeq 2\left(-1+\frac{1}{\lambda(\phi)} \right)\epsilon_1$ at leading order, with $\lambda (\phi)$ being a function of the scalar field which depends on the scalar field potential and the second derivative of the scalar-Gauss-Bonnet coupling $\xi''(\phi)$. With our analysis we aim to provide a definitive criterion for selecting Einstein-Gauss-Bonnet models that can provide a blue-tilted inflationary phenomenology, by simply looking at the scalar potential and the scalar-Gauss-Bonnet coupling. We shall prove this using two distinct approaches and as we show the tilt of the tensor spectral index is determined by the values of the potential $V(\phi)$ and of scalar-Gauss-Bonnet coupling at first horizon crossing. Specifically the blue-tilted tensor spectral index can occur when $\xi''(\phi_*)V(\phi_*)>0$ at first horizon crossing.

In this work we aim to revive the interest for non-minimal derivative coupling theories of gravity, in light of the GW170817 event. These theories include a string motivated non-minimal kinetic term for the scalar field of the form $\sim \xi(\phi) G^{\mu\nu}\partial_\mu\phi\partial_\nu\phi$ and predict that the primordial tensor perturbations have a speed that it is distinct from the speed of light. Due to the fact that the Universe is classical during and in the post-inflationary epoch, there is no fundamental reason for the graviton to change its mass, the GW170817 event severely constrained these theories. We analyze and formalize the inflationary phenomenology of these theories, using both the latest Planck data and the GW170817 event data to constrain these theories. Due to the fact that there is no constraint in the choice of the scalar potential and the non-minimal coupling function, we provide several classes of viable models, using convenient forms of the minimal coupling function in terms of the scalar potential, aiming for analyticity, and we discuss the advantages and disadvantages of each viable model.

Recent pulsar timing array collaborations have reported evidence of the stochastic gravitational wave background. The gravitational waves induced by primordial curvature perturbations, referred to as scalar-induced gravitational waves (SIGWs), could potentially be the physical origins of this gravitational wave background. In addition to the statistical properties of SIGWs derived from the primordial fluctuations, SIGWs also have intrinsic non-Gaussianity originating from the non-linear interactions of Einstein's gravity, which, however, is rarely explored. In this paper, we study the intrinsic non-Gaussianity of SIGWs with the bispectrum and the skewness under the assumption of the Gaussian primordial curvature perturbations, phenomenologically modeled as a lognormal spectrum. The bispectrum is shown to be vanishing in the collinear limit, which is independent of the initial conditions and the dynamics of SIGWs. For the SIGWs generated during the radiation-dominated era, the bispectrum is of flatten-type non-Gaussianity, with only four polarization components left to be non-vanishing. Additionally, in order to obtain the correct bispectrum, we also propose a time-oscillation average scheme and a regularization scheme. Utilizing the skewness for quantifying the degree of non-Gaussianity, it is found that the curvature power spectrum with a narrow width can result in an enhancement of the third-order non-Gaussianity. The conclusion holds for both the SIGWs generated in the radiation-dominated era and the matter-dominated era.

TianQin is a proposed mission for space-based gravitational-wave detection that features a triangular constellation in circular high Earth orbits. The mission entails three drag-free controlled satellites and long-range laser interferometry with stringent beam pointing requirements at remote satellites. For the payload architecture and pointing control strategies, having two test masses per satellite, one for each laser arm, and rotating entire opto-mechanical assemblies (each consisting of a telescope, an optical bench, an inertial sensor, etc.) for constellation breathing angle compensation represent an important option for TianQin. In this paper, we examine its applicability from the perspectives of test mass and satellite control in the science mode, taking into account of perturbed orbits and orbital gravity gradients. First, based on the orbit-attitude coupling relationship, the required electrostatic forces and torques for the test mass suspension control are estimated and found to be sufficiently small for the acceleration noise budget. Further optimization favors configuring the centers of masses of the two test masses collinear and equidistant with the center of mass of the satellite, and slightly offsetting the assembly pivots from the electrode housing centers forward along the sensitive axes. Second, the required control forces and torques on the satellites are calculated, and thrust allocation solutions are found under the constraint of having a flat-top sunshield on the satellite with varying solar angles. The findings give a green light to adopting the two test masses and telescope pointing scheme for TianQin.

Pointing-related displacement noises are crucial in space-based gravitational wave detectors, where point-ahead angle control of transmitted laser beams may contribute significantly. For TianQin that features a geocentric concept, the circular high orbit design with a nearly fixed constellation plane gives rise to small variations of the point-ahead angles within $\pm 25$ nrad in-plane and $\pm 10$ nrad off-plane, in addition to a static bias of 23 $\mu$rad predominantly within the constellation plane. Accordingly, TianQin may adopt fixed-value compensation for the point-ahead angles and absorb the small and slow variations into the pointing biases. To assess the in-principle feasibility, the far-field tilt-to-length (TTL) coupling effect is discussed, and preliminary requirements on far-field wavefront quality are derived, which have taken into account of TTL noise subtraction capability in post processing. The proposed strategy has benefits in simplifying the interferometry design, payload operation, and TTL noise mitigation for TianQin.

We generalize, imposing the field equations only at dominant order, the Isaacson formula for the gravitational wave (GW) energy-momentum tensor (EMT) to the class of Horndeski theories in which the tensor modes travel at the speed of light (reduced Horndeski theories) and scalar waves are present. We discuss important particular cases such as: theories where scalar waves are also luminal and theories in which the transverse-traceless gauge can be achieved in an arbitrary open set. The vanishing of the trace of the gravitational wave energy-momentum tensor is obtained for theories in which all wave perturbations propagate at the speed of light. The trace is shown not to vanish trivially in other cases. We obtain, as a particular case of our general result, the GW EMTs, in a Brans-Dicke theory, both in the Einstein frame, recovering previous results in the literature, and in the Jordan frame, thereby showing the GW EMT is not conformally invariant. We further prove that there exists a subclass of reduced Horndeski theories where, in contrast to general relativity, the divergence of the GW EMT does not vanish even after the imposition of the full equations of motion, assuming an eikonal solution.

The ringdown portion of a binary black hole merger consists of a sum of modes, each containing an infinite number of tones that are exponentially damped sinusoids. In principle, these can be measured as gravitational-waves with observatories like LIGO/Virgo/KAGRA, however in practice it is unclear how many tones can be meaningfully resolved. We investigate the consistency and resolvability of the overtones of the quadrupolar l = m = 2 mode by starting at late times when the gravitational waveform is expected to be well-approximated by the l m n = 220 tone alone. We present a Bayesian inference framework to measure the tones in numerical relativity data. We measure tones at different start times, checking for consistency: we classify a tone as stably recovered if and only if the 95% credible intervals for amplitude and phase at time t overlap with the credible intervals at all subsequent times. We test the first four overtones of the fundamental mode and find that the 220 and 221 tones can be measured consistently with the inclusion of additional overtones. The 222 tone measurements can be stabilised when we include the 223 tone, but only in a narrow time window, after which it is too weak to measure. The 223 tone recovery appears to be unstable, and does not become stable with the introduction of the 224 tone. Within our framework, the ringdown of the fundamental mode of binary black hole waveforms can be self-consistently described by three to four tones, which are stable from 10 M after the peak strain. However, additional tones are not obviously required because the fit amplitudes are consistent with zero. We conclude that recent claims of overtone detection are not necessarily an exercise in over-fitting; the observed tones can be self-consistently modelled, although, not until ~ 10 M from peak strain with our four-tone model.

Gravitational wave (GW) astronomy has opened new frontiers in understanding the cosmos, while the integration of artificial intelligence (AI) in science promises to revolutionize data analysis methodologies. However, a significant gap exists, as there is currently no dedicated platform that enables scientists to develop, test, and evaluate AI algorithms efficiently. To address this gap, we introduce GWAI, a pioneering AI-centered software platform designed for gravitational wave data analysis. GWAI contains a three-layered architecture that emphasizes simplicity, modularity, and flexibility, covering the entire analysis pipeline. GWAI aims to accelerate scientific discoveries, bridging the gap between advanced AI techniques and astrophysical research.

Alfv\'en waves (AWs) excited by the cosmic-ray (CR) streaming instability (CRSI) are a fundamental ingredient for CR confinement. The effectiveness of self-confinement relies on a balance between CRSI growth rate and damping mechanisms acting on quasi-parallel AWs excited by CRs. One relevant mechanism is the so-called turbulent damping, in which an AW packet injected in pre-existing turbulence undergoes a cascade process due to its nonlinear interaction with fluctuations of the background. The turbulent damping of an AW packet in pre-existing magnetohydrodynamic turbulence is re-examined, revised, and extended to include most-recent theories of MHD turbulence that account for dynamic alignment and reconnection-mediated regime. The case in which the role of feedback of CR-driven AWs on pre-existing turbulence is important will also be discussed. Particular attention is given to the nonlinearity parameter $\chi^w$ that estimates the strength of nonlinear interaction between CR-driven AWs and background fluctuations. We point out the difference between $\chi^w$ and $\chi^z$ that instead describes the strength of nonlinear interactions between pre-existing fluctuations. When $\chi^w$ is properly taken into account, one finds that (i) the turbulent damping rate of quasi-parallel AWs in anisotropic turbulence depends on the background-fluctuations' amplitude to the third power, hence is strongly suppressed, and (ii) the dependence on the AW's wavelength (and thus on the CR gyro-radius from which it is excited) is different from what has been previously obtained. Finally, (iii) when dynamic alignment of cascading fluctuations and the possibility of a reconnection-mediated range is included in the picture, the turbulent damping rate exhibits novel regimes and breaks. Finally, a criterion for CR-feedback is derived and simple phenomenological models of CR-modified turbulent scaling are provided.

Hydrodynamic and magnetohydrodynamic convective attractors in three-dimensional rotating Rayleigh-B\'enard convection are studied numerically by varying the Taylor and Rayleigh numbers as control parameters. First, an analysis of hydrodynamic attractors and their bifurcations is conducted, where routes to chaos via quasiperiodicity are identified. Second, the behaviour of the magnetohydrodynamic system is investigated by introducing a seed magnetic field and measuring its growth or decay as a function of the Taylor number, while keeping the Rayleigh number fixed. Analysis of the attractors shows that rotation has a significant impact on magnetic field generation in Rayleigh-B\'enard convection, with the critical magnetic Prandtl number changing nonmonotonically with the rotation rate. It is argued that a nonhysteretic blowout bifurcation with on-off intermittency is responsible for the transitions to dynamo.

We study the accretion flows towards a central Kerr super-spinning attractor, discussing the formation of the flow inversion points, defined by condition $u^{\phi}=0$ on the particles flow axial velocity. We locate two closed surfaces, defining \emph{inversion coronas} (spherical shells), surrounding the central attractor. The coronas analysis highlights observational aspects distinguishing the central attractors and providing indications on their spin and the orbiting fluids. The inversion corona is a closed region, generally of small extension and thickness, which is for the counter-rotating flows of the order of $\lesssim 1.4 M$ (central attractor mass) on the vertical rotational axis. There are no co-rotating inversion points (from co-rotating flows). The results point to strong signatures of the Kerr super-spinars, provided in both accretion and jet flows. With very narrow thickness, and varying little with the fluid initial conditions and the emission process details, inversion coronas can have remarkable observational significance for primordial Kerr super-spinars predicted by string theory. The corona region closest to the central attractor is the most observably recognizable and active part, distinguishing black holes solutions from super-spinars. Our analysis expounds the Lense--Thirring effects and repulsive gravity effects in the super-spinning ergoregions.

A first order electroweak phase transition probes physics beyond the Standard Model on multiple frontiers and therefore is of immense interest for theoretical exploration. We conduct a model-independent study of the effects of relevant dimension 6 operators, of the Standard Model Effective Field Theory, on electroweak phase transition. We use a thermally corrected and renormalization group improved potential and study its impact on nucleation temperature. We then outline bubble dynamics that lead to ultra-relativistic bubble wall velocities which are mainly motivated from the viewpoint of gravitational wave detection. We highlight the ranges of the Wilson coefficients that give rise to such bubble wall velocities and predict gravitational wave spectra generated by such transitions which can be tested in future experiments.

A modified gravity model of Starobinsky inflation and primordial black hole production was proposed in good (within $1\sigma$) agreement with current measurements of the cosmic microwave background radiation. The model is an extension of the singularity-free Appleby-Battye-Starobinsky model by the $R^4$-term with different values of the parameters whose fine-tuning leads to efficient production of primordial black holes on smaller scales with the asteroid-size masses between $10^{16}$ g and $10^{20}$ g. Those primordial black holes may be part (or the whole) of the current dark matter, while the proposed model can be confirmed or falsified by detection or absence of the induced gravitational waves with the frequencies about $10^{-2}$ Hz. The main purpose of this paper was to estimate the size of quantum (loop) corrections in the model. It was found their relative contribution to the power spectrum of scalar perturbations is about $10^{-5}$ or less, so that the model is not ruled out by the quantum corrections.

We investigate the intricate relationships between the non-radial \(f\) mode oscillation frequencies of neutron stars (NS)s and the corresponding nuclear matter equation of state (EOS) using a machine learning (ML) approach within the ambit of the relativistic mean field (RMF) framework for nuclear matter. With two distinct parameterizations of the Walecka model, namely, (1) with non-linear self interactions of the scalar field (NL) and, (2) a density dependent Bayesian model (DDB), we perform a thorough examination of the \(f\) mode frequency in relation to various nuclear saturation properties. The correlations between the \(f\) mode frequencies and nuclear saturation properties reveal, through various analytical and ML methods, the complex nature of NSs and their potential as the cosmic laboratory for studying extreme states of matter. A principal component analysis (PCA) has been performed using mixed datasets from DDB and NL models to discriminate the relative importance of the different components of the EOS on the $f$ mode frequencies. Additionally, a {\it Random forest feature importance} analysis also elucidates the distinct roles of these properties in determining the \(f\) mode frequency across a spectrum of NS masses. Our findings are further supported by symbolic regression searches, yielding high-accuracy relations with strong Pearson coefficients and minimal errors. These relations suggest new methodologies for probing NS core characteristics, such as energy density, pressure, and speed of sound from observations of non-radial \(f\) mode oscillations of NSs.

Advancements in cosmology through next-generation ground-based gravitational wave observatories will bring in a paradigm shift. We explore the pivotal role that gravitational-wave standard sirens will play in inferring cosmological parameters with next-generation observatories, not only achieving exquisite precision but also opening up unprecedented redshifts. We examine the merits and the systematic biases involved in gravitational-wave standard sirens utilizing binary black holes, binary neutron stars, and neutron star-black hole mergers. Further, we estimate the precision of bright sirens, golden dark sirens, and spectral sirens for these binary coalescences and compare the abilities of various next-generation observatories (A^sharp, Cosmic Explorer, Einstein Telescope, and their possible networks). When combining different sirens, we find sub-percent precision over more than 10 billion years of cosmic evolution for the Hubble expansion rate $H(z)$. This work presents a broad view of opportunities to precisely measure the cosmic expansion rate, decipher the elusive dark energy and dark matter, and potentially discover new physics in the uncharted Universe with next-generation gravitational-wave detectors.

We present a novel mechanism for producing topologically stable monopoles (TSMs) from the quantum mechanical decay of metastable cosmic strings in the early universe. In an $SO(10)$ model this mechanism yields TSMs that carry two units ($4\pi/e$) of Dirac magnetic charge as well as some color magnetic charge which is screened. For a dimensionless string tension parameter $G\mu \approx 10^{-9} - 10^{-5}$, the monopoles are superheavy with masses of order $10^{15} - 10^{17}$ GeV. Monopoles with masses of order $10^8 - 10^{14}$ GeV arise from metastable strings for $G\mu$ values from $\sim 10^{-22}$ to $10^{-10}$. We identify the parameter space for producing these monopoles at an observable level with detectors such as IceCube and KM3NeT. For lower $G\mu$ values the ultra-relativistic monopoles should be detectable at Pierre Auger and ANITA. The stochastic gravitational wave emission arises from metastable strings with $G\mu\sim 10^{-9}-10^{-5}$ and should be accessible at HLVK and future detectors including the Einstein Telescope and Cosmic Explorer. An $E_6$ extension based on this framework would yield TSMs from the quantum mechanical decay of metastable strings that carry three units ($6\pi/e$) of Dirac magnetic charge.