Papers for Thursday, Nov 02 2023

4U 1210-64 is a peculiar X-ray binary that was first discovered in 1978 by the Uhuru satellite. The analysis of the X-ray data revealed a 6.7-day orbital period and an additional long-term modulation that is manifested as low and high flux states. Based on the previous classification of the donor star from the analysis of its optical spectra, the system has been suggested to be a high-mass X-ray binary. We re-visit the optical classification where, based on the spectra from the Southern African Large Telescope (SALT), we conclude that the donor star is of spectral class A8 III-IV, making it a member of the rare intermediate-mass X-ray binaries. We perform radial velocity analysis using the SALT spectra where we consider circular and eccentric orbits. From the mass function derived and the mass constraints of the donor star, we demonstrate that a neutron star is favoured as the compact object in the binary system. We show, for the first time, the folded optical lightcurves, whose shape is interpreted to be due to a combination of ellipsoidal variations, irradiation of the donor star, and mutual eclipses of the star and accretion disk.

Accreting supermassive black holes (SMBHs) produce highly magnetized relativistic jets that tend to collimate gradually as they propagate outward. However, recent radio interferometric observations of the 3C 84 galaxy reveal a stunning, cylindrical jet already at several hundred SMBH gravitational radii, $r\gtrsim350r_{\rm g}$. We explore how such extreme collimation emerges via a suite of 3D general-relativistic magnetohydrodynamic (GRMHD) simulations. We consider an SMBH surrounded by a magnetized torus immersed in a constant-density ambient medium that starts at the edge of the SMBH sphere of influence, chosen to be much larger than the SMBH gravitational radius, $r_{\text{B}}=10^3r_{\text{g}}$. We find that radiatively inefficient accretion flows (e.g., M87) produce winds that collimate the jets into parabolas near the BH. After the disk winds stop collimating the jets at $r\lesssim{}r_\text{B}$, they turn conical. Once outside $r_\text{B}$, the jets run into the ambient medium and form backflows that collimate the jets into cylinders some distance beyond $r_{\text{B}}$. Interestingly, for radiatively-efficient accretion, as in 3C 84, the radiative cooling saps the energy out of the disk winds: at early times, they cannot efficiently collimate the jets, which skip the initial parabolic collimation stage, start out conical near the SMBH, and turn into cylinders already at $r\simeq300r_{\rm g}$, as observed in 3C 84. Over time, jet power remains approximately constant, whereas the mass accretion rate increases: the winds grow in strength and start to collimate the jets, which become quasi-parabolic near the base; the transition point to a nearly cylindrical jet profile moves outward while remaining inside $r_\text{B}$.

We report the detection of cold dust in an apparently quiescent massive galaxy ($\log({M_{\star}/M_{\odot}})\approx11$) at $z\sim2$ (G4). The source is identified as a serendipitous 2 mm continuum source in a deep ALMA observation within the field of Q2343-BX610, a $z=2.21$ massive star-forming disk galaxy. Available multi-band photometry of G4 suggests redshift of $z\sim2$ and a low specific star-formation rate (sSFR), $\log(SFR/M_{\star}) [yr^{-1}] \approx -10.2$, corresponding to $\approx1.2$ dex below the $z=2$ main sequence (MS). G4 appears to be a peculiar dust-rich quiescent galaxy for its stellar mass ($\log({M_{\rm dust}/M_{\star}}) = -2.71 \pm 0.26$), with its estimated mass-weighted age ($\sim$ 1-2 Gyr). We compile $z\gtrsim1$ quiescent galaxies in the literature and discuss their age-$\Delta$MS and $\log({M_{\rm dust}/M_{\star}})$-age relations to investigate passive evolution and dust depletion scale. A long dust depletion time and its morphology suggest morphological quenching along with less efficient feedback that could have acted on G4. The estimated dust yield for G4 further supports this idea, requiring efficient survival of dust and/or grain growth, and rejuvenation (or additional accretion). Follow-up observations probing the stellar light and cold dust peak are necessary to understand the implication of these findings in the broader context of galaxy evolutionary studies and quenching in the early universe.

We present here estimates of the average rates of accretion of neutral gas onto main-sequence galaxies and the conversion of atomic gas to molecular gas in these galaxies at two key epochs in galaxy evolution: (i) $z\approx1.3-1.0$, towards the end of the epoch of peak star-formation activity in the Universe, and (ii) $z\approx1-0$, when the star-formation activity declines by an order of magnitude. We determine the net gas accretion rate $\rm{R_{Acc}}$ and the molecular gas formation rate $\rm{R_{Mol}}$ by combining the relations between the stellar mass and the atomic gas mass, the molecular gas mass, and the star-formation rate (SFR) at three epochs, $z=1.3$, $z=1.0$, and $z=0$, with the assumption that galaxies evolve continuously on the star-forming main-sequence. We find that, for all galaxies, $\rm{R_{Acc}}$ is far lower than the average SFR $\rm{R_{SFR}}$ at $z\approx1.3-1.0$; however, $\rm{R_{Mol}}$ is similar to $\rm{R_{SFR}}$ during this interval. Conversely, both $\rm{R_{Mol}}$ and $\rm{R_{Acc}}$ are significantly lower than $\rm{R_{SFR}}$ over the later interval, $z\approx1-0$. We find that massive main-sequence galaxies had already acquired most of their present-day baryonic mass by $z\approx1.3$. At $z\approx1.3-1.0$, the rapid conversion of the existing atomic gas to molecular gas was sufficient to maintain a high average SFR, despite the low net gas accretion rate. However, at later times, the combination of the lower net gas accretion rate and the lower molecular gas formation rate leads to a decline in the fuel available for star-formation, and results in the observed decrease in the SFR density of the Universe over the last 8 Gyr.

We present a set of six general relativistic, multi-frequency, radiation magnetohydrodynamic simulations of thin accretion disks with different target mass accretion rates around black holes with spins ranging from non-rotating to rapidly spinning. The simulations use the $\mathbf{M}_1$ closure scheme with twelve, independent frequency (or energy) bins ranging logarithmically from $5\times 10^{-3}$ to $5\times 10^3$ keV. The multi-frequency capability allows us to generate crude spectra and energy-dependent light curves directly from the simulations without a need for special post-processing. While we generally find roughly thermal spectra with peaks around 1 to 4 keV, our high-spin cases showed harder than expected tails for the soft or thermally dominant state. This leads to radiative efficiencies that are up to five times higher than expected for a Novikov-Thorne disk at the same spin. We attribute these high efficiencies to the high-energy, coronal emission. These coronae mostly occupy the effectively optically thin regions near the inner edges of the disks and also cover or sandwich the inner $\sim 15 GM/c^2$ of the disks.

We develop a model independent, robust method for determining galaxy rotation velocities across a 2D array of spaxels from an integral field spectrograph. Simulations demonstrate the method is accurate down to lower spectral signal-to-noise than standard methods: 99\% accurate when median $S/N=4$. We apply it to MaNGA data to construct the galaxy velocity map and galaxy rotation curve. We also develop a highly efficient cubic smoothing approach that is $25\times$ faster computationally and only slightly less accurate. Such model independent methods could be useful in studying dark matter properties without assuming a galaxy model.

We study out-of-thermodynamic equilibrium effects in neutron star mergers with 3D general-relativistic neutrino-radiation large-eddy simulations. During merger, the cores of the neutron stars remain cold ($T \sim$ a few MeV) and out of thermodynamic equilibrium with trapped neutrinos originating from the hot collisional interface between the stars. However, within ${\sim}2{-}3$ milliseconds matter and neutrinos reach equilibrium everywhere in the remnant. Our results show that dissipative effects, such as bulk viscosity, if present, are only active for a short window of time after the merger.

Jetted astrophysical phenomena with black hole (BH) engines, including binary mergers, jetted tidal disruption events, and X-ray binaries, require a large-scale vertical magnetic field for efficient jet formation. However, a dynamo mechanism that could generate these crucial large-scale magnetic fields has not been identified and characterized. We have employed 3D global general relativistic magnetohydrodynamical (MHD) simulations of accretion disks to quantify, for the first time, a dynamo mechanism that generates large-scale magnetic fields. This dynamo mechanism primarily arises from the nonlinear evolution of the magnetorotational instability (MRI). In this mechanism, large non-axisymmetric MRI-amplified shearing wave modes, mediated by the axisymmetric azimuthal magnetic field, generate and sustain the large-scale vertical magnetic field through their nonlinear interactions. We identify the advection of magnetic loops as a crucial feature, transporting the large-scale vertical magnetic field from the outer regions to the inner regions of the accretion disk. This leads to a larger characteristic size of the, now advected, magnetic field when compared to the local disk height. We characterize the complete dynamo mechanism with two timescales: one for the local magnetic field generation, $t_{\rm g}$, and one for the large-scale scale advection, $t_{\rm adv}$. Whereas the dynamo we describe is nonlinear, we explore the potential of linear mean field models to replicate its core features. Our findings indicate that traditional $\alpha$-dynamo models, often computed in stratified shearing box simulations, are inadequate and that the effective large-scale dynamics is better described by the shear current effects or stochastic $\alpha$-dynamos.

A key event in the history of the Milky Way is the formation of the bar. This event affects the subsequent structural and dynamical evolution of the entire Galaxy. When the bar formed, gas was likely rapidly funnelled to the centre of the Galaxy settling in a star-forming nuclear disc. The Milky Way bar formation can then be dated by considering the oldest stars in the formed nuclear stellar disc. In this highly obscured and crowded region, reliable age tracers are limited, but bright, high-amplitude Mira variables make useful age indicators as they follow a period--age relation. We fit dynamical models to the proper motions of a sample of Mira variables in the Milky Way's nuclear stellar disc region. Weak evidence for inside-out growth and both radial and vertical dynamical heating with time of the nuclear stellar disc is presented suggesting the nuclear stellar disc is dynamically well-mixed. Furthermore, for Mira variables around a $\sim350$ day period, there is a clear transition from nuclear stellar disc-dominated kinematics to background bar-bulge-dominated kinematics. Using a Mira variable period-age relation calibrated in the solar neighbourhood, this suggests the nuclear stellar disc formed in a significant burst in star formation $(8\pm 1)\,\mathrm{Gyr}$ ago, although the data are also weakly consistent with a more gradual formation of the nuclear stellar disc at even earlier epochs. This implies a relatively early formation time for the Milky Way bar ($\gtrsim8\,\mathrm{Gyr}$), which has implications for the growth and state of the young Milky Way and its subsequent history.

We propose a new approach to improve the precision of astrophysical parameter constraints for the 21cm signal from Epoch of Reionization (EoR). Our method introduces new sets of summary statistics, hereafter evolution compressed statistics, that quantify the spectral evolution of the 2D spatial statistics computed a fixed redshift. We defined such compressed statistics for Power Spectrum (PS), as well as Reduced Wavelet Scattering Transform (RWST) and Wavelet Moments (WM), which also characterise non-Gaussian features. To compare these different statistics with fiducial 3D power spectrum, we estimate their Fisher information on three cosmological parameters from an ensemble of simulations of 21cm EoR data, both in noiseless and noisy scenarios using Square Kilometre Array (SKA) noise levels equivalent to 100 and 1000 hours of observations. For the noiseless case, the compressed wavelet statistics give constraints up to five times higher precision than the 3D isotropic power spectrum, while for 100h SKA noise, for which non-Gaussian features are hard to extract, they still give constraints which are 30% better. From this study, we demonstrate that evolution-compressed statistics extract more information than usual 3D isotropic approaches and that our wavelet-based statistics can consistently outmatch power spectrum-based statistics. When constructing such wavelet-based statistics, we also emphasise the need to choose a set of wavelets with an appropriate spectral resolution concerning the astrophysical process studied.

Despite major progress in our understanding of massive stars, concerning discrepancies still remain between observations and theory. Most notable are the numerous stars observed beyond the theoretical main sequence, an evolutionary phase expected to be short-lived and hence sparsely populated. This is the "Blue Supergiant Problem." Stellar models with abnormal internal structures can provide long-lived solutions for this problem: core hydrogen-burning stars with oversized cores may explain the hotter ones, and core helium-burning stars with undersized cores may explain the cooler ones. Such stars may result from enhanced or suppressed mixing in single stars, or, more likely, as the products of binary interaction and stellar mergers. Here we investigate the potential of asteroseismology to uncover the nature of blue supergiants. We construct stellar models for the above scenarios and show that they predict $g$-mode period spacings that differ by an order of magnitude: $\sim$ 200~min versus $\sim$ 20~min for long-lived core H and He burning stars, respectively. For the classical scenario of H-shell burning stars rapidly crossing the Hertzsprung gap, we furthermore predict changes of the order of $10^{-2}\,\mu\rm{Hz\,yr}^{-1}$ in high-frequency modes; this effect would be in principle observable from $\sim 5$~yr of asteroseismic monitoring if these modes can be identified. This raises the possibility of revealing the internal structure of blue supergiants, and thus determining whether these stars are indeed binary merger products. These asteroseismic diagnostics may be measurable through long timeseries observations from the ongoing TESS mission and upcoming PLATO mission, thereby laying a path toward resolving the blue supergiant problem.

We conduct a statistical study of black hole masses of barred and unbarred galaxies in the IllustrisTNG100 cosmological magneto-hydrodynamical simulations. This work aims to understand the role of the bars in the growth of central supermassive black hole mass and its implications on AGN fueling. Our sample consists of 1191 barred galaxies and 2738 unbarred galaxies in the IllustrisTNG100 simulations. To have an unbiased study, we perform our analysis with an equal number of barred and unbarred galaxies by using various controlled parameters like total galaxy mass, stellar mass, gas mass, dark matter halo mass, etc. Except for the stellar mass controlling, we find that the median of the black hole mass distribution for barred galaxies is higher than that of the unbarred ones indicating that stellar mass is a key parameter influencing the black hole growth. The higher mean accretion rate of the black holes in barred galaxies, averaged since the bar forming epoch (z~2 ), explains the higher mean black hole masses in barred galaxies. Further, we also test that these results are unaffected by other environmental processes like minor/major merger histories and neighboring gas density of black hole. Although the relationship between stellar mass, bar formation, and black hole growth is complex, with various mechanisms involved, our analysis suggests that bars can play a crucial role in feeding black holes, particularly in galaxies with massive stellar disks.

We examine the stellar [$\alpha$/Fe]-[Fe/H] distribution of $\simeq1000$ present-day galaxies in a high-resolution EAGLE simulation. A slight majority of galaxies exhibit the canonical distribution, characterised by a sequence of low-metallicity stars with high [$\alpha$/Fe] that transitions at a "knee" to a sequence of declining [$\alpha$/Fe] with increasing metallicity. This population yields a knee metallicity - galaxy-mass relation similar to that observed in Local Group galaxies, both in slope and scatter. However, many simulated galaxies lack a knee or exhibit more complicated distributions. Knees are found only in galaxies with star formation histories (SFHs) featuring a sustained decline from an early peak ($t\simeq7~{\rm Gyr}$), which enables enrichment by Type Ia supernovae to dominate that due to Type II supernovae (SN II), reducing [$\alpha$/Fe] in the interstellar gas. The simulation thus indicates that, contrary to the common interpretation implied by analytic galactic chemical evolution (GCE) models, knee formation is not a consequence of the onset of enrichment by SN Ia. We use the SFH of a simulated galaxy exhibiting a knee as input to the VICE GCE model, finding it yields an $\alpha$-rich plateau enriched only by SN II, but the plateau comprises little stellar mass and the galaxy forms few metal-poor ([Fe/H]$\lesssim - 1$) stars. This follows from the short, constant gas consumption timescale typically assumed by GCEs, which implies the presence of a readily-enriched, low-mass gas reservoir. When an initially longer, evolving consumption timescale is adopted, VICE reproduces the simulated galaxy's track through the [$\alpha$/Fe]-[Fe/H] plane and its metallicity distribution function.

Because of the electromagnetic radiation produced during the merger, compact binary coalescences with neutron stars may result in multi-messenger observations. In order to follow up on the gravitational-wave signal with electromagnetic telescopes, it is critical to promptly identify the properties of these sources. This identification must rely on the properties of the progenitor source, such as the component masses and spins, as determined by low-latency detection pipelines in real time. The output of these pipelines, however, might be biased, which could decrease the accuracy of parameter recovery. Machine learning algorithms are used to correct this bias. In this work, we revisit this problem and discuss two new implementations of supervised machine learning algorithms, K-Nearest Neighbors and Random Forest, which are able to predict the presence of a neutron star and post-merger matter remnant in low-latency compact binary coalescence searches across different search pipelines and data sets. Additionally, we present a novel approach for calculating the Bayesian probabilities for these two metrics. Instead of metric scores derived from binary machine learning classifiers, our scheme is designed to provide the astronomy community well-defined probabilities. This would deliver a more direct and easily interpretable product to assist electromagnetic telescopes in deciding whether to follow up on gravitational-wave events in real time.

A detailed analysis of ground-based CCD UBV photometry and space-based Gaia Data Release 3 (DR3) data for the open clusters King 6 and NGC 1605 was performed. Using the pyUPMASK algorithm on Gaia astrometric data to estimate cluster membership probabilities, we have identified 112 stars in King 6 and 160 stars in NGC 1605 as the statistically most likely members of each cluster. We calculated reddening and metallicity separately using UBV two-color diagrams to estimate parameter values via independent methods. The color excess $E(B-V)$ and photometric metallicity [Fe/H] for King 6 are $0.515 \pm 0.030$ mag and $0.02 \pm 0.20$ dex, respectively. For NGC 1605, they are $0.840 \pm 0.054$ mag and $0.01 \pm 0.20$ dex. With reddening and metallicity kept constant, we have estimated the distances and cluster ages by fitting PARSEC isochrones to color-magnitude diagrams based on the Gaia and UBV data. Photometric distances are 723 $\pm$ 34 pc for King 6 and 3054 $\pm$ 243 pc for NGC 1605. The cluster ages are $200 \pm 20$ Myr and $400 \pm 50$ Myr for King 6 and NGC 1605, respectively. Mass function slopes were found to be 1.29 $\pm$ 0.18 and 1.63 $\pm$ 0.36 for King 6 and NGC 1605, respectively. These values are in good agreement with the value of Salpeter (1955). The relaxation times were estimated as 5.8 Myr for King 6 and 60 Myr for NGC 1605. This indicates that both clusters are dynamically relaxed since these times are less than the estimated cluster ages. Galactic orbit analysis shows that both clusters formed outside the solar circle and are members of the young thin-disc population.

In every proposed unification scheme for Active Galactic Nuclei (AGN), an integral element is the presence of circumnuclear dust arranged in torus-like structures. A crucial model parameter in this context is the covering factor (CF), defined as the ratio between the infrared luminosity of the dusty torus $L_{\rm IR}$, and the accretion disk bolometric luminosity $L_{\rm agn}$. Our study aims to determine whether CF evolution is genuine or if selection effects significantly influence it. Based on cross-matched multiwavelength photometrical data from the five major surveys (SDSS, GALEX, UKIDSS, WISE, SPITZER), a sample of almost 2,000 quasars was derived. The main parameters of quasars, such as black hole masses and the Eddington ratios, were calculated based on the spectroscopic data. The data were divided into two redshift bins: Low-$z$ (redshift ~1) and High-$z$ (redshift ~2) quasars. We identified an issue with the accuracy of the WISE W4 filter. Whenever feasible, it is recommended to utilize SPITZER MIPS 24 $\mu$m data. The calculated median CF values for the highest quality SPITZER data are comparable within errors $\log$ CF$_{\textrm{low}-z} = -0.19\pm 0.11$ and $\log$ CF$_{\textrm{high}-z}= -0.18\pm 0.11$. The Efron & Petrosian test confirmed the presence of luminosity evolution with redshift for both $L_{\rm IR}$ and $L_{\rm agn}$. Both the Low-$z$ and High-$z$ samples exhibit a similar correlation between $L_{\rm agn}$ and $L_{\rm IR}$. No discernible evolution of the CF was observed in the subsample of quasars with high SMBH mass bin or high luminosities. The relationship between $L_{\rm IR}$ and $L_{\rm agn}$ deviates slightly from the expected 1:1 scaling. However, no statistically significant dependence of CF on luminosities could be claimed across the entire dataset.

The role of type Ia supernovae, mainly the Delay Time Distributions (DTDs) determined by the binary systems, and the yields of elements created by different explosion mechanisms, are studied by using the {\sc MulChem} chemical evolution model, applied to our Galaxy. We explored 15 DTDs, and 12 tables of elemental yields produced by different SN Ia explosion mechanisms, doing a total of 180 models. Chemical abundances for $\alpha$-elements (O, Mg, Si, S, Ca) and Fe derived from these models, are compared with recent observational data of $\alpha$-elements over Iron relative abundances, [X/Fe]. These data have been compiled and binned in 13 datasets. By using a $\chi^2$-technique, no model is able to fit simultaneously these datasets. A model computed with the 13 individual best models is good enough to reproduce them. Thus, a power law with a logarithmic slope $\sim -1.1$ and a delay in the range $\Delta \tau=40 --350$ Myr is a possible DTD, but a combination of several channels is more probable. Results of this average model for other disc regions show a high dispersion, as observed, which might be explained by the stellar migration. The dispersion might also come from a combination of DTDs or of explosion channels. The stellar migration joined to a combination of scenarios for SNIa is the probable cause of the observed dispersion.

In 2013 a dearth of close-in planets around fast-rotating host stars was found using statistical tests on Kepler data. The addition of more Kepler and Transiting Exoplanet Survey Satellite (TESS) systems in 2022 filled this region of the diagram of stellar rotation period (Prot) versus the planet orbital period (Porb). We revisited the Prot extraction of Kepler planet-host stars, we classify the stars by their spectral type, and we studied their Prot-Porb relations. We only used confirmed exoplanet systems to minimize biases. In order to learn about the physical processes at work, we used the star-planet evolution code ESPEM (French acronym for Evolution of Planetary Systems and Magnetism) to compute a realistic population synthesis of exoplanet systems and compared them with observations. Because ESPEM works with a single planet orbiting around a single main-sequence star, we limit our study to this population of Kepler observed systems filtering out binaries, evolved stars, and multi-planets. We find in both, observations and simulations, the existence of a dearth in close-in planets orbiting around fast-rotating stars, with a dependence on the stellar spectral type (F, G, and K), which is a proxy of the mass in our sample of stars. There is a change in the edge of the dearth as a function of the spectral type (and mass). It moves towards shorter Prot as temperature (and mass) increases, making the dearth look smaller. Realistic formation hypotheses included in the model and the proper treatment of tidal and magnetic migration are enough to qualitatively explain the dearth of hot planets around fast-rotating stars and the uncovered trend with spectral type.

In 2021, a catalog of 536 fast radio bursts (FRBs) detected with the Canadian Hydrogen Intensity Mapping Experiment (CHIME) radio telescope was released by the CHIME/FRB Collaboration. This large collection of bursts, observed with a single instrument and uniform selection effects, has advanced our understanding of the FRB population. Here we update the results for 140 of these FRBs for which channelized raw voltage (baseband) data are available. With the voltages measured by the telescope's antennas, it is possible to maximize the telescope sensitivity in any direction within the primary beam, an operation called beamforming. This allows us to increase the signal-to-noise ratio (S/N) of the bursts and to localize them to sub-arcminute precision. The improved localization is also used to correct the beam response of the instrument and to measure fluxes and fluences with a ~ 10% uncertainty. Additionally, the time resolution is increased by three orders of magnitude relative to that in the first CHIME/FRB catalog, and, applying coherent dedispersion, burst morphologies can be studied in detail. Polarization information is also available for the full sample of 140 FRBs, providing an unprecedented dataset to study the polarization properties of the population. We release the baseband data beamformed to the most probable position of each FRB. These data are analyzed in detail in a series of accompanying papers.

This paper investigates a hypothesis proposed in previous research relating neutron star (NS) mass and its dark matter (DM) accumulation. As DM accumulates, NS mass decreases, predicting lower NS masses toward the Galactic center. Due to limited NSs data near the galactic center, we examine NSs located within DM clumps. Using the CLUMPY code simulations, we determine the DM clumps distribution, with masses from 10 to $10^{8}$ $M_{\odot}$ and scales from $10^{-3}$ to 10 kpc. These clumps' DM exhibit a peak at the center, tapering toward the outskirts, resembling our Galaxy's DM distribution. We analyse these DM clumps' NS mass variations, considering diverse DM particle masses and galaxy types. We find relatively stable NS mass within 0.01 to 5 kpc from the clump center. This stability supports the initial hypothesis, particularly for NSs located beyond 0.01 kpc from the clump center, where NS mass reaches a plateau around 0.1 kpc. Nevertheless, NS mass near the clump's periphery reveals spatial dependence: NS position within DM clumps influences its mass in Milky Way-type galaxies. Moreover, this dependence varies with the DM model considered. In summary, our study investigates the proposed link between NS mass and DM accumulation by examining NSs within DM clumps. While NS mass remains stable at certain distances from the clump center, spatial dependencies arise near the clump's outer regions, contingent on the specific DM model.

We present the CO(3-2) APEX survey at 6 pc resolution of the bar of the SMC. We aboard the CO analysis in the SMC-Bar comparing the CO(3-2) survey with that of the CO(2-1) of similar resolution. We study the CO(3-2)-to-CO(2-1) ratio (R32) that is very sensitive to the environment properties (e.g., star-forming regions). We analyzed the correlation of this ratio with observational quantities that trace the star formation as the local CO emission, the Spitzer color [70/160], and the total IR surface brightness measured from the Spitzer and Herschel bands. For the identification of the CO(3-2) clouds, we used the CPROPS algorithm, which allowed us to measure the physical properties of the clouds. We analyzed the scaling relationships of such physical properties. We obtained an R32 of 0.65 as a median value for the SMC, with a standard deviation of 0.3. We found that R32 varies from region to region, depending on the star formation activity. In regions dominated by HII and photo-dissociated regions (e.g., N22, N66), R32 tends to be higher than the median values. Meanwhile, lower values were found toward quiescent clouds. We also found that R32 correlates positively with the IR color [70/160] and the total IR surface brightness. This finding indicates that R32 increases with environmental properties like the dust temperature, the total gas density, and the radiation field. We have identified 225 molecular clouds with sizes R > 1.5 pc and signal-to-noise (S/N) ratio > 3 and only 17 well-resolved CO(3-2) clouds increasing the S/N ratio to > 5. These 17 clouds follow consistent scaling relationships to the inner Milky Way clouds but with some departure. The CO(3-2) tends to be less turbulent and less luminous than the inner Milky Way clouds of similar size. Finally, we estimated a median virial-based CO-to-H2 conversion factor of 12.6_{-7}^{+10} Msun/(K km s^{-1} pc^{2}) for the total sample.

During early phases of a protoplanetary disks's life, gravitational instabilities can produce significant mass transport, can dramatically alter disk structure, can mix and shock-process gas and solids, and may be instrumental in planet formation. We present a 3D grid-based radiative hydrodynamics study with varied resolutions of a 0.07 M$_\odot$ disk orbiting a 0.5 M$_\odot$ star as it settles over most of its radial extent into a quasi-steady asymptotic state that maintains approximate balance between heating produced by GIs and radiative cooling governed by realistic dust opacities. We assess disk stability criteria, thermodynamic properties, strengths of GIs, characteristics of density waves and torques produced by GIs, radial mass transport arising from these torques, and the level to which transport can be represented as local or nonlocal processes. Physical and thermal processes display distinct differences between inner optically thick and outer optically thin regions of the disk. In the inner region, gravitational torques are dominated by low-order Fourier components of the azimuthal mass distribution. These torques are strongly variable on the local dynamical time and are subject to rapid flaring presumably driven by recurrent swing amplification. In the outer region, m=1 torques dominate. Ring-like structures exhibiting strong noncircular motions, and vortices develop near the inner edge between 8 and 14 au. We find that GI-induced spiral modes erupt in a chaotic manner over the whole low-Q part of the disk, with many spiral modes appearing and disappearing, producing gravitoturbulence, but dominated by fluctuating large-scale modes, very different from a simple $\alpha$-disk.

The connection between galaxies and dark matter halos is often quantified using the stellar mass-halo mass (SMHM) relation. Optical and near-infrared imaging surveys have led to a broadly consistent picture of the evolving SMHM relation based on measurements of galaxy abundances and angular correlation functions. Spectroscopic surveys at $z \gtrsim 2$ can also constrain the SMHM relation via the galaxy autocorrelation function and through the cross-correlation between galaxies and Ly$\alpha$ absorption measured in transverse sightlines; however, such studies are very few and have produced some unexpected or inconclusive results. We use $\sim$3000 spectra of $z\sim2.5$ galaxies from the Lyman-alpha Tomography IMACS Survey (LATIS) to measure the galaxy-galaxy and galaxy-Ly$\alpha$ correlation functions in four bins of stellar mass spanning $10^{9.2} \lesssim M_* / M_{\odot} \lesssim 10^{10.5}$. Parallel analyses of the MultiDark N-body and ASTRID hydrodynamic cosmological simulations allow us to model the correlation functions, estimate covariance matrices, and infer halo masses. We find that results of the two methods are mutually consistent and are broadly in accord with standard SMHM relations. This consistency demonstrates that we are able to accurately measure and model Ly$\alpha$ transmission fluctuations $\delta_F$ in LATIS. We also show that the galaxy-Ly$\alpha$ cross-correlation, a free byproduct of optical spectroscopic galaxy surveys at these redshifts, can constrain halo masses with similar precision to galaxy-galaxy clustering.

Large-scale astronomical surveys can capture numerous images of celestial objects, including galaxies and nebulae. Analysing and processing these images can reveal intricate internal structures of these objects, allowing researchers to conduct comprehensive studies on their morphology, evolution, and physical properties. However, varying noise levels and point spread functions can hamper the accuracy and efficiency of information extraction from these images. To mitigate these effects, we propose a novel image restoration algorithm that connects a deep learning-based restoration algorithm with a high-fidelity telescope simulator. During the training stage, the simulator generates images with different levels of blur and noise to train the neural network based on the quality of restored images. After training, the neural network can directly restore images obtained by the telescope, as represented by the simulator. We have tested the algorithm using real and simulated observation data and have found that it effectively enhances fine structures in blurry images and increases the quality of observation images. This algorithm can be applied to large-scale sky survey data, such as data obtained by LSST, Euclid, and CSST, to further improve the accuracy and efficiency of information extraction, promoting advances in the field of astronomical research.

We propose an effective model to describe the backreaction on cosmological observables induced by Laniakea, the gravitational supercluster hosting the Milky Way, which was defined using peculiar velocity data from Cosmicflows-4 (CF4). The structure is well described by an ellipsoidal shape exhibiting triaxial expansion, reasonably approximated by a constant expansion rate along the principal axes. Our best fits suggest that the ellipsoid, after subtracting the background expansion, contracts along the two smaller axes and expands along the longest one, predicting an average expansion of $\sim -1.1 ~\rm{km}/\rm{s}/\rm{Mpc}$. The different expansion rates within the region, relative to the mean cosmological expansion, induce line-of-sight-dependent corrections in the computation of luminosity distances. We apply these corrections to two low-redshift datasets: the Pantheon+ catalog of type Ia Supernovae (SN~Ia), and 63 measurements of Surface Brightness Fluctuations (SBF) of early-type massive galaxies from the MASSIVE survey. We find corrections on the distances of order $\sim 2-3\%$, resulting in a shift in the inferred best-fit values of the Hubble constant $H_0$ of order $\Delta H_0^{\rm{SN~Ia}}\approx 0.5 ~\rm{km}/\rm{s}/\rm{Mpc}$ and $\Delta H_0^{\rm{SBF}}\approx 1.1 ~\rm{km}/\rm{s}/\rm{Mpc}$, seemingly worsening the Hubble tension.

A coherently oscillating ultra-light axion can behave as dark matter. In particular, its coherently oscillating pressure perturbations can source an oscillating scalar metric perturbation, with a characteristic oscillation frequency which is twice the axion Compton frequency. A candidate in the mass range $10^{(-24,-21)}{\rm eV}$ can provide a signal in the frequency range tested by current and future Pulsar Timing Array (PTA) programs. Involving the pressure perturbations in a highly nonlinear environment, such an analysis demands a relativistic and nonlinear treatment. Here, we provide a rigorous derivation of the effect assuming weak gravity and slow-motion limit of Einstein's gravity in zero-shear gauge and show that dark matter's velocity potential determines the oscillation phase and frequency change. A monochromatic PTA signal correlated with the velocity field would confirm the prediction, for example, by cross-correlating the PTA results with the future local velocity flow measurements.

We present a new catalog of Kepler planet candidates that prioritizes accuracy of planetary dispositions and properties over uniformity. This catalog contains 4376 transiting planet candidates, including 1791 residing within 709 multi-planet systems, and provides the best parameters available for a large sample of Kepler planet candidates. We also provide a second set of stellar and planetary properties for transiting candidates that are uniformly-derived for use in occurrence rates studies. Estimates of orbital periods have been improved, but as in previous catalogs, our tabulated values for period uncertainties do not fully account for transit timing variations (TTVs). We show that many planets are likely to have TTVs with long periodicities caused by various processes, including orbital precession, and that such TTVs imply that ephemerides of Kepler planets are not as accurate on multi-decadal timescales as predicted by the small formal errors (typically 1 part in $10^6$ and rarely $ > 10^{-5}$) in the planets' measured mean orbital periods during the Kepler epoch. Analysis of normalized transit durations implies that eccentricities of planets are anti-correlated with the number of companion transiting planets. Our primary catalog lists all known Kepler planet candidates that orbit and transit only one star; for completeness, we also provide an abbreviated listing of the properties of the two dozen non-transiting planets that have been identified around stars that host transiting planets discovered by Kepler.

Inflation with Weyl scaling symmetry provides a viable scenario that can generate both the nearly scaling invariant primordial density fluctuation and a dark matter candidate. Here we point out that, in additional to the primordial gravitational waves (GWs) from quantum fluctuations, the production of high-frequency GWs from preheating in such inflation models can provide an another probe of the inflationary dynamics. We conduct both linear analytical analysis and nonlinear numerical lattice simulation in a typical model. We find that significant stochastic GWs can be produced and the frequency band is located around $10^8$ Hz $\sim$ $10^9$ Hz, which might be probed by future resonance-cavity experiments.

The possible formation histories of neutron star binaries remain unresolved by current gravitational-wave catalogs. The detection of an eccentric binary system could be vital in constraining compact binary formation models. We present the first search for aligned spin eccentric neutron star-black hole binaries (NSBH) and the most sensitive search for aligned-spin eccentric binary neutron star (BNS) systems using data from the third observing run of the advanced LIGO and advanced Virgo detectors. No new statistically significant candidates are found; we constrain the local merger rate to be less than 150 $\text{Gpc}^{-3}\text{Yr}^{-1}$ for binary neutron stars in the field, and, 50, 100, and 70 $\text{Gpc}^{-3}\text{Yr}^{-1}$ for neutron star-black hole binaries in globular clusters, hierarchical triples and nuclear clusters, respectively, at the 90$\%$ confidence level if we assume that no sources have been observed from these populations. We predict the capabilities of upcoming and next-generation observatory networks; we investigate the ability of three LIGO ($\text{A}^{#}$) detectors and Cosmic Explorer CE (20km) + CE (40km) to use eccentric binary observations for determining the formation history of neutron star binaries. We find that 2 -- 100 years of observation with three $\text{A}^{#}$ observatories are required before we observe clearly eccentric NSBH binaries; this reduces to only 10 days -- 1 year with the CE detector network. CE will observe tens to hundreds of measurably eccentric binaries from each of the formation models we consider.

High-energy neutrino astrophysics is rapidly developing, and in the last two years, new and exciting results have been obtained. Among them are the confirmation of the existence of the diffuse astrophysical neutrino flux by the new independent Baikal-GVD experiment, the discovery of the neutrino emission of our Galaxy, new confirmations of the origin of a part of astrophysical neutrinos in blazars, and much more. This brief review, based on the author's presentation at the session of the RAS Physical Science Division "Gamma quanta and neutrinos from space: what we can see now and what we need to see more", summarizes the results obtained since the publication of the review arXiv:2112.09611, and can be considered as a companion to it.

Be X-ray binaries (Be XRBs) are high-mass X-ray binaries, with a neutron star or black hole orbiting and accreting material from a non-supergiant B-star that is rotating at a near critical rate. These objects are prime targets to understand past binary interactions as the neutron star or black hole progenitor likely experienced Roche lobe overflow to spin up the Be star we observe now. The stellar variability can then allow us to explore the stellar structure of these objects. It was recently demonstrated that the high-mass X-ray binary CPD -29 2176 descended from an ultra-stripped supernova and is a prime target to evolve into an eventual binary neutron star and kilonova. We present the photometric variability from both TESS and ASAS along with the spectral properties and disk variability of the system in this paper. All of the optical lines are contaminated with disk emission except for the He II $\lambda$4686 absorption line. The disk variability time-scales are not the same as the orbital time scale, but could be related to the X-ray outbursts that have been recorded by Swift. We end our study with a discussion comparing CPD -29 2176 to classical Be stars and other Be X-ray binaries, finding the stellar rotation to be near a frequency of 1.5 cycles d$^{-1}$, and exhibiting incoherent variability in three frequency groups.

The unexpected discovery of hot Jupiters challenged the classical theory of planet formation inspired by our solar system. Until now, the origin and evolution of hot Jupiters are still uncertain. Determining their age distribution and temporal evolution can provide more clues into the mechanism of their formation and subsequent evolution. Using a sample of 383 giant planets around Sun-like stars collected from the kinematic catalogs of the Planets Across Space and Time (PAST) project, we find that hot Jupiters are preferentially hosted by relatively younger stars in the Galactic thin disk. We subsequently find that the frequency of hot Jupiters declines with age. In contrast, the frequency of warm/cold Jupiters shows no significant dependence on age. Such a trend is expected from the tidal evolution of hot Jupiters' orbits, and our result offers supporting evidence using a large sample. We also perform a joint analysis on the planet frequencies in the stellar age-metallicity plane. The result suggests that the frequencies of hot Jupiters and warm/cold Jupiters, after removing the age dependence are both correlated with stellar metallicities. Moreover, we show that the above correlations can explain the bulk of the discrepancy in hot Jupiter frequencies inferred from the transit and radial velocity (RV) surveys, given that RV targets tend to be more metal-rich and younger than transits.

The patchiness in the reionization process alters the statistics of Cosmic Microwave Background (CMB), with the kinematic Sunyaev-Zeldovich (kSZ) effect in the CMB temperature power spectrum being a notable consequence. In this work, we aim to explore the potential of future kSZ power spectrum measurements in inferring the details of the reionization process. In this pursuit, we capitalize on the recent developments in foreground mitigation techniques using the Cross-Internal Linear Combination (Cross-ILC) technique, which enables robust detection of the kSZ power spectrum with signal-to-noise ($S/N$) roughly $20-30\sigma$ in this decade by SPT-3G and Simons Observatory (SO); and $\geq 80\sigma$ by CMB-S4, substantially improving on the marginal evidence for kSZ binned at $\ell=3000$ using SPT data (Reichardt et al. 2021). We use a fiducial kSZ power spectrum along with realistic error bars expected from the above technique for SPT-3G, SO, and CMB-S4 to constrain the parameter space for a physical model of reionization. We find that with the improved error bars it will be possible to place stringent constraints on reionization using solely the Cross-ILC recovered SPT-3G kSZ without imposing any prior on $\tau$ in the Bayesian inference. Notably, high-fidelity kSZ measurements from CMB-S4 coupled with $\tau$ measurements through LiteBIRD will enable unprecedented constraint on the midpoint of reionization with an error bar of $\sim 0.25$ and the duration of reionization with an error bar at $\sim 0.21$ exclusively using CMB data. This study highlights the need to capture kSZ power on a broad range of multipoles to gain insights into the inhomogeneous reionization era.

Weak gravitational lensing is an invaluable tool for understanding fundamental cosmological physics. An unresolved issue in weak lensing cosmology is to accurately reconstruct the lensing convergence $\kappa$ maps from discrete shear catalog with survey masks, which the seminal Kaiser-Squire (KS) method is not designed to address. We present the Accurate Kappa Reconstruction Algorithm for masked shear catalog (AKRA) to address the issue of mask. AKRA is built upon the prior-free maximum likelihood mapmaking method (or the unbiased minimum variance linear estimator). It is mathematically robust in dealing with mask, numerically stable to implement, and practically effective in improving the reconstruction accuracy. Using simulated maps with mask fractions ranging from 10\% to 50\% and various mask shapes, we demonstrate that AKRA outperforms KS at both the map level and summary statistics such as the auto power spectrum $C_\kappa$ of the reconstructed map, its cross-correlation coefficient $r_\ell$ with the true $\kappa $ map, the scatter plot and the localization measure. Unlike the Wiener filter method, it adopts no priors on the signal power spectrum, and therefore avoids the Wiener filter related biases at both the map level and cross-correlation statistics. If we only use the reconstructed map in the unmasked regions, the reconstructed $C_\kappa$ is accurate to $1\%$ or better and $1-r_\ell \lesssim 1\%$ (excluding $\ell$ at the smallest scales investigated), even for extreme cases of mask fraction and shape. As the first step, the current version of AKRA only addresses the mask issue and therefore ignores complexities such as curved sky and inhomogeneous shape measurement noise. AKRA is capable of dealing with these issues straightfowrardly, and will be addressed in the next version.

We report a new, rare detection of HI 21-cm absorption associated with a quasar (only six known at $1<z<2$) here towards J2339-5523 at $z_{em}$ = 1.3531, discovered through the MeerKAT Absorption Line Survey (MALS). The absorption profile is broad ($\sim 400$ km/s), and the peak is redshifted by $\sim 200$ km/s, from $z_{em}$. Interestingly, optical/FUV spectra of the quasar from Magellan-MIKE/HST-COS spectrographs do not show any absorption features associated with the 21-cm absorption. This is despite the coincident presence of the optical quasar and the radio `core' inferred from a flat spectrum component of flux density $\sim 65$ mJy at high frequencies ($>5$ GHz). The simplest explanation would be that no large HI column (N(HI)$>10^{17}$ cm$^{-2}$) is present towards the radio `core' and the optical AGN. Based on the joint optical and radio analysis of a heterogeneous sample of 16 quasars ($z_{median}$ = 0.7) and 15 radio galaxies ($z_{median}$ = 0.3) with HI 21-cm absorption detection and matched in 1.4 GHz luminosity (L$_{\rm 1.4\,GHz}$), a consistent picture emerges where quasars are primarily tracing the gas in the inner circumnuclear disk and cocoon created by the jet-ISM interaction. These exhibit L$_{1.4\,\rm GHz}$ - $\Delta V_{\rm null}$ correlation, and frequent mismatch between the radio and optical spectral lines. The radio galaxies show no such correlation and likely trace the gas from the cocoon and the galaxy-wide ISM outside the photoionization cone. The analysis presented here demonstrates the potential of radio spectroscopic observations to reveal the origin of the absorbing gas associated with AGN that may be missed in optical observations.

We present the discovery and timing results of five pulsars discovered in a pilot survey at intermediate Galactic latitudes with the Five-hundred Aperture Spherical Telescope (FAST). Among these pulsars, two belong to the category of millisecond pulsars (MSPs) with spin periods of less than 20\,ms. Two others fall under the classification of `mildly recycled' pulsars, with massive white dwarfs as companions. Remarkably, this small survey, covering an area of 4.7 square degrees, led to the discovery of five pulsars, including four recycled pulsars. Such success underscores the immense potential of future surveys at intermediate Galactic latitudes. In order to assess the potential yield of MSPs, we conducted population simulations and found that both FAST and Parkes new Phased Array Feed surveys, focusing on intermediate Galactic latitudes, have the capacity to uncover several hundred new MSPs.

The accuracy of stellar distances inferred purely from parallaxes degrades rapidly with distance. Proper motion measurements, when combined with some idea of typical velocities, provide independent information on stellar distances. Here I build a direction- and distance-dependent model of the distribution of stellar velocities in the Galaxy, then use this together with parallaxes and proper motions to infer kinegeometric distances and transverse velocities for stars in Gaia DR3. Using noisy simulations I assess the performance of the method and compare its accuracy to purely parallax-based (geometric) distances. Over the whole Gaia catalogue, kinegeometric distances are on average 1.25 times more accurate than geometric ones. This average masks a large variation in the relative performance, however. Kinegeometric distances are considerably better than geometric ones beyond several kpc, for example. On average, kinegeometric distances can be measured to an accuracy of 19% and velocities (sqrt[vra^2 + vdec^2]) to 16 km/s (median absolute deviations). In Gaia DR3, kinegeometric distances are smaller than geometric ones on average for distant stars, but the pattern is more complex in the bulge and disk. With the much more accurate proper motions expected in Gaia DR5, a further improvement in the distance accuracy by a factor of (only) 1.35 on average is predicted (with kinegeometric distances still 1.25 times more accurate than geometric ones). The improvement from proper motions is limited by the width of the velocity prior, in a way that the improvement from better parallaxes is not limited by the width of the distance prior.

We present a multi-wavelength study of the brightest cluster galaxies (BCGs) in a sample of the 95 most massive galaxy clusters selected from South Pole Telescope (SPT) Sunyaev-Zeldovich (SZ) survey. Our sample spans a redshift range of 0.3 < z < 1.7, and is complete with optical spectroscopy from various ground-based observatories, as well as ground and space-based imaging from optical, X-ray and radio wavebands. At z~0, previous studies have shown a strong correlation between the presence of a low-entropy cool core and the presence of star-formation and a radio-loud AGN in the central BCG. We show for the first time that a central entropy threshold for star formation persists out to z~1. The central entropy (measured in this work at a radius of 10 kpc) below which clusters harbor star-forming BCGs is found to be as low as $K_\mathrm{10 ~ kpc} = 35 \pm 4$ keV cm$^2$ at z < 0.15 and as high as $K_\mathrm{10 ~ kpc} = 52 \pm 11$ keV cm$^2$ at z~1. We find only marginal (~1$\sigma$) evidence for evolution in this threshold. In contrast, we do not find a similar high-z analog for an entropy threshold for feedback, but instead measure a strong evolution in the fraction of radio-loud BCGs in high-entropy cores as a function of redshift. This could imply that the cooling-feedback loop was not as tight in the past, or that some other fuel source like mergers are fueling the radio sources more often with increasing redshift, making the radio luminosity an increasingly unreliable proxy for radio jet power. We also find that our SZ-based sample is missing a small (~4%) population of the most luminous radio sources ($\nu L_{\nu} > 10^{42}$ erg/s), likely due to radio contamination suppressing the SZ signal with which these clusters are detected.

Ejecta from Dimorphos following the DART mission impact, significantly increased the brightness of the Didymos-Dimorphos system, allowing us to examine sub-surface material. We report daily near-IR spectroscopic observations of the Didymos system using NASA's IRTF, that follow the evolution of the spectral signature of the ejecta cloud over one week, from one day before the impact. Overall, the spectral features remained fixed (S-type classification) while the ejecta dissipated, confirming both Didymos and Dimorphos are constructed from the same silicate material. This novel result strongly supports binary asteroid formation models that include breaking up of a single body, due to rotational breakup of km-wide bodies. At impact time +14 and +38 hours, the spectral slope decreased, but following nights presented increasing spectral slope that almost returned to the pre-impact slope. However, the parameters of the $1~\mu m$ band remained fixed, and no "fresh" / Q-type-like spectrum was measured. We interpret these as follow: 1. The ejecta cloud is the main contributor ($60-70\%$) to the overall light during the $\sim40$ hours after impact. 2. Coarser debris ($\geq 100~\mu m$) dominated the ejecta cloud, decreasing the spectral slope (after radiation pressure removed the fine grains at $\leq10$ hours after impact); 3. after approximately one week, the ejecta cloud dispersed enough to make the fine grains on Didymos surface the dominating part of the light, increasing the spectral slope to pre-impact level. 4. a negligible amount of non-weathered material was ejected from Dimorphos' sub-surface, suggesting Dimorphos was accumulated from weathered material, ejected from Didymos surface.

Magnetic fields regulate black hole (BH) accretion, governing both the inflow and outflow dynamics. When a BH becomes saturated with large-scale vertical magnetic flux, it enters a magnetically-arrested disk (MAD) state. The dynamically-important BH magnetic flux powers highly efficient relativistic outflows (or jets) and sporadically erupts from the BH into the disk midplane. Here we explore the evolution of MADs when the BH and gas angular momentum are misaligned, which is expected to be more common. Using numerical simulations, we find that jets from rapidly spinning, prograde BHs force the inner accretion flow into alignment with the BH spin via the magneto-spin alignment mechanism for disks initially misaligned at $\mathcal{T}\lesssim 60^{\circ}$. Extremely misaligned MAD disks, on the other hand, exhibit intermittent jets that blow out parts of the disk to $\approx 100$ gravitational radii before collapsing, leaving behind hot cavities and magnetized filaments. These intermittent jet mechanism forms a mini-feedback cycle and could explain some cases of X-ray and radio quasi-periodic eruptions observed in dim AGN. Further, we find that (i) for BHs with low power jets, the BH spin and initial disk tilt angle changes the amount of horizon magnetic flux, and (ii) geometrically-thick, misaligned accretion flows do not undergo sustained Lense-Thirring (LT) precession. Thereby, we suggest that low-luminosity accreting BHs ($\dot{M}\ll 10^{-3} \dot{M}_{\rm Edd}$) are not likely to exhibit quasi-periodic oscillations in lightcurves due to LT precession, in agreement with observations of BH X-ray binaries and AGN in the low-hard/quiescent state. Instead, we suggest that magnetic flux eruptions can mimic precession-like motion, such as observed in the M87 jet, by driving large-scale surface waves in the jets.

Partial eruptions of solar filaments are the typical representative of solar eruptive behavior diversity. Here we investigate a typical filament partial eruption event and present integrated evidence for configuration of the pre-eruption filament and its formation. The CHASE H$\alpha$ observations reveal structured Doppler velocity distribution within the pre-eruption filament, where distinct redshift only appeared in the east narrow part of the south filament region and then disappeared after the partial eruption while the north part dominated by blueshift remained. Combining the SDO, ASO-S observations, and NLFFF modeling results, we verify that there were two independent material flow systems within the pre-flare filament, whose magnetic topology is a special double-decker configuration consisting of two magnetic flux ropes (MFRs) with opposite magnetic twist. During the formation of this filament system, continuous magnetic flux cancellation and footpoint motion were observed around its north end. Therefore, we propose a new double-decker formation scenario that the two MFRs composing such double-decker configuration originated from two magnetic systems with different initial connections and opposite magnetic twist. Subsequent magnetic reconnection with surrounding newly-emerging fields resulted in the motion of footpoint of the upper MFR to the region around footpoint of the lower MFR, thus leading to eventual formation of the double-decker configuration consisting of two MFRs with similar footpoints but opposite signs of magnetic twist. These results provide a potential way to determine unambiguously the progenitor configuration of a partial-eruptive filament and reveal a special type of double-decker MFR configuration and a new double-decker formation scenario.

It is shown that the description of the solar cycle that takes into account the odd zonal harmonic of the solar magnetic field allows us to deepen our knowledge of two important aspects of the solar activity. First, to clarify and expand predictions of the evolution of the cyclic activity of the Sun in the near future. Second, to develop a program for monitoring the spectrophotometric characteristics of radiation of the solar-type stars aimed at obtaining new information about their magnetic fields.

We present a multiwavelength investigation of the hub-filament system RCW 117 (IRAS 17059-4132), which shows intricate filamentary features in the far-infrared, mapped using Herschel images. We obtain the column density and dust temperature maps for the region using the Herschel images, and identify 88 cores and 12 filaments from the column density map of the region ($18'\times18'$). We employ the ThrUMMS $^{13}$CO (J=1-0) data for probing the kinematics in RCW 117, and find velocity gradients ($\sim 0.3-1$ km s$^{-1}$ pc$^{-1}$) with hints of matter inflow along the filamentary structures. Ionised gas emission from the associated HII region is examined using the Giant Metrewave Radio Telescope (GMRT) at 610 and 1280 MHz, and is found to be of extent $5 \times 3$ pc$^2$ with intensity being brightest towards the hub. We estimate the peak electron density towards the hub to be $\sim 750$ cm$^{-3}$. Thirty four Class 0/I young stellar objects (YSOs) have been identified in the region using the Spitzer GLIMPSE colour-colour diagram, with many lying along the filamentary structures. Based on the (i) presence of filamentary structures, (ii) distribution of cores across the region, with $\sim39$% found along the filamentary structures, (iii) massive star-formation tracers in the hub, and (iv) the kinematics, we believe that global hierarchical collapse can plausibly explain the observed features in RCW 117.

Charged particles are constantly accelerated to non-thermal energies by the reconnecting magnetic field in the solar atmosphere. Our understanding of the interactions between the particles and their environment can benefit from three-dimensional atmospheric simulations accounting for non-thermal particle beams. In a previous publication, we presented the first results from such a simulation. However, the original treatment of beam propagation ignores potentially important phenomena. Here, we present a more general beam propagation model incorporating magnetic gradient forces, the return current, acceleration by the ambient electric field, and temperature-dependent collision rates. Neglecting collisional velocity randomisation makes the model sufficiently lightweight to simulate millions of beams. We investigate how each new physical effect changes beam energy transport in a three-dimensional atmosphere. We applied the method of characteristics to the steady-state continuity equation for electron flux to derive ordinary differential equations for the mean evolution of energy, pitch angle, and flux with distance. For each beam, we solved these numerically for a range of initial energies to obtain the evolving flux spectrum, from which we computed the energy deposited into the ambient plasma. Magnetic gradient forces significantly influence the deposition of beam energy. The strong magnetic field convergence leads to a small coronal deposition peak followed by a heavy dip caused by the onset of magnetic mirroring. Mirrored electrons carry away 5 to 10% of the injected beam energy on average. The transition region peak produced by the remaining energetic electrons occurs slightly deeper than in a uniform magnetic field. An initial diverging magnetic field enhances the subsequent impact of mirroring. The other new physical effects are less significant for the studied atmospheric conditions.

Deep learning (DL) has recently been proposed as a novel approach for 21cm foreground removal. Before applying DL to real observations, it is essential to assess its consistency with established methods, its performance across various simulation models and its robustness against instrumental systematics. This study develops a commonly used U-Net and evaluates its performance for post-reionisation foreground removal across three distinct sky simulation models based on pure Gaussian realisations, the Lagrangian perturbation theory, and the Planck sky model. Stable outcomes across the models are achieved provided that training and testing data align with the same model. On average, the residual foreground in the U-Net reconstructed data is $\sim$10% of the signal across angular scales at the considered redshift range. Comparable results are found with traditional approaches. However, blindly using a network trained on one model for data from another model yields inaccurate reconstructions, emphasising the need for consistent training data. The study then introduces frequency-dependent Gaussian beams and gain drifts to the test data. The network struggles to denoise data affected by "unexpected" systematics without prior information. However, after re-training consistently with systematics-contaminated data, the network effectively restores its reconstruction accuracy. This highlights the importance of incorporating prior systematics knowledge during training for successful denoising. Our work provides critical guidelines for using DL for 21cm foreground removal, tailored to specific data attributes. Notably, it is the first time that DL has been applied to the Planck sky model being most realistic foregrounds at present.

The centers of galaxies host a supermassive black hole surrounded by a dense stellar cluster. The cluster is expected to develop mass segregation, in which gravitational scatterings among the stars cause heavier objects to sink closer to the central black hole, while lighter objects will tend to be over concentrated in the outer regions. This work focuses on the implications of mass segregation on the different channels for violent destruction of stars in the cluster: tidal disruptions, gravitational-wave driven inspirals and high-velocity destructive collisions between stars. All such events occur close to the central black hole, where the heavier objects congregate. The analysis is based on a simplified Monte Carlo simulation, which evolves a two-mass population in a cluster surrounding a Milky-Way-like super massive black hole. The simulation is based on the single mass scheme used by Sari and Fragione (2019) and Balberg and Yassur (2023), which has been extended to allow for the dynamical friction effects typical of non-equal mass populations. The effects of mass segregation on the rates of the different destruction channels are analyzed self-consistently in the overall evolution of the cluster. Also considered are stars which are injected into the cluster after being disrupted from a binary system by the supermassive black hole. Such stars are captured in the inner regions of the cluster, and so their orbital evolution, and their destruction rate, are therefore influenced by heavy objects that might be abundant in the vicinity of the supermassive black hole.

Interstellar dust (ISD) particles penetrate the solar system due to the relative motion of the Sun and the local interstellar cloud. Before entering the heliosphere, they pass through the heliospheric interface - the region of the solar wind interaction with the interstellar plasma. The size distribution and number density of dust grains are modified in the interface essentially. The modification depends on the charging of the dust particles along their trajectories. In this paper, we present modeling results of the charging of ISD particles passing through the heliospheric interface. The main physical processes responsible for the charging within the heliospheric conditions are the sticking of primary plasma particles, secondary electron emission, photoemission, and the effects of cosmic ray electrons. We consider two methods to calculate the electric charge of ISD particles based on (1) the classic steady-state assumption that the charge depends only on local plasma and radiation conditions and (2) the dynamical computation of charge along the particle trajectory. We demonstrate that the steady-state assumption is quite justified to model trajectories and number density distributions of relatively big ISD grains (radius of 100 nm and larger) penetrating the heliosphere. The estimates show that ISD grains of these sizes require less than 0.25 years (distance of ~ 1 au) after transition from the LISM into the heliosphere to reach an equilibrium. For small particles (radius of 10 nm), the dynamical computation of charge influences the trajectories and modifies the number density substantially. The dust density accumulations are distributed within a more elongated region along the heliopause in case of dynamically changed charge as compared with the use of a steady-state charge approximation.

We present the results of our Keck/DEIMOS spectroscopic follow-up of candidate galaxies of i-band-dropout protocluster candidate galaxies at $z\sim6$ in the COSMOS field. We securely detect Lyman-$\alpha$ emission lines in 14 of the 30 objects targeted, 10 of them being at $z=6$ with a signal-to-noise ratio of $5-20$, the remaining galaxies are either non-detections or interlopers with redshift too different from $z=6$ to be part of the protocluster. The 10 galaxies at $z\approx6$ make the protocluster one of the riches at $z>5$. The emission lines exhibit asymmetric profiles with high skewness values ranging from 2.87 to 31.75, with a median of 7.37. This asymmetry is consistent with them being Ly$\alpha$, resulting in a redshift range of $z=5.85-6.08$. Using the spectroscopic redshifts, we re-calculate the overdensity map for the COSMOS field and find the galaxies to be in a significant overdensity at the $4\sigma$ level, with a peak overdensity of $\delta=11.8$ (compared to the previous value of $\delta=9.2$). The protocluster galaxies have stellar masses derived from Bagpipes SED fits of $10^{8.29}-10^{10.28} \rm \,M_{\rm \odot}$ and star formation rates of $2-39\,\rm M_{\rm \odot}\rm\,yr^{-1}$, placing them on the main sequence at this epoch. Using a stellar-to-halo-mass relationship, we estimate the dark matter halo mass of the most massive halo in the protocluster to be $\sim 10^{12}\rm M_{\rm \odot}$. By comparison with halo mass evolution tracks from simulations, the protocluster is expected to evolve into a Virgo- or Coma-like cluster in the present day.

The final sizes, composition, and angular momenta of solid planetary bodies depend on the outcomes of collisions between planetary embryos. The most common numerical method for simulating embryo collisions is to combine a gravity solver with a hydrodynamic solver, allowing pressure gradients, shock waves, and gravitational torques to loft material into orbit. Here, we perform the first direct comparison between hydrodynamic methods and a simplified method employing only gravity and a quadratic repulsive force. The formation of Earth's Moon, perhaps the most heavily simulated planetary collision, is used as a test case. Many of the main features of a collision between two planetary embryos, including collisions in which an orbiting disc of material and/or intact moons are formed, are controlled solely by gravitational forces. Comparison of the methods shows that the mass and orbit of the satellite, as well as the extent of physical mixing between the protoearth and impactor, are similar regardless of the inclusion of the inclusion of hydrodynamic effects or the equation of state employed. The study of thermal and chemical effects of the impact, and determining the time scale for lunar accretion, still require a full hydrodynamic calculation. The simplified gravity plus quadratic repulsive force approach allows rapid testing of various initial conditions to identify cases for further detailed study.

Axion-like early dark energy (EDE) as an extension to $\Lambda$CDM has been proposed as a possible solution to the 'Hubble tension'. We revisit this model using a new cosmic microwave background (CMB) temperature and polarization likelihood constructed from the {\it Planck} NPIPE data release. In a Bayesian analysis, we find that the maximum fractional contribution of EDE to the total energy density is $f_{\rm EDE} < 0.061$ (without SH0ES) over the redshift range $z\in[10^3,10^4]$ and that the Hubble constant is constrained to lie within the range $ 66.9 < H_0 < 69.5$ km/s/Mpc (both at 95 \% C.L.). The data therefore favour a model close to $\Lambda$CDM, leaving a residual tension of $3.7\sigma$ with the SH$0$ES Cepheid-based measurement of $H_0$. A comparison with the likelihood profile shows that our conclusions are robust to prior-volume effects. Our new CMB likelihood provides no evidence in favour of a significant EDE component.

To investigate the relative amount of ejecta from high-mass versus intermediate-mass stars and to trace the chemical evolution of the Galaxy, we have performed with the IRAM 30m and the SMT 10m telescopes a systematic study of Galactic interstellar 18O/17O ratios toward a sample of 421 molecular clouds, covering a galactocentric distance range of 1-22 kpc. The results presented in this paper are based on the J=2-1 transition and encompass 364 sources showing both C18O and C17O detections. The previously suggested 18O/17O gradient is confirmed. For the 41 sources detected with both facilities, good agreement is obtained. A correlation of 18O/17O ratios with heliocentric distance is not found, indicating that beam dilution and linear beam sizes are not relevant. For the subsample of IRAM 30 m high-mass star-forming regions with accurate parallax distances, an unweighted fit gives 18O/17O = (0.12+-0.02)R_GC+(2.38+-0.13) with a correlation coefficient of R = 0.67. While the slope is consistent with our J=1-0 measurement, ratios are systematically lower. This should be caused by larger optical depths of C18O 2-1 lines, w.r.t the corresponding 1-0 transitions, which is supported by RADEX calculations and the fact that C18O/C17O is positively correlated with 13CO/C18O. After considering optical depth effects with C18O J=2-1 reaching typically an optical depth of 0.5, corrected 18O/17O ratios from the J=1-0 and J=2-1 lines become consistent. A good numerical fit to the data is provided by the MWG-12 model, including both rotating stars and novae.

Photonics offer numerous functionalities that can be used to realize astrophotonic instruments. The most spectacular example to date is the ESO Gravity instrument at the Very Large Telescope in Chile. Integrated astrophotonic devices stand to offer critical advantages for instrument development, including extreme miniaturization, as well as integration, superior thermal and mechanical stabilization owing to the small footprint, and high replicability offering cost savings. Numerous astrophotonic technologies have been developed to address shortcomings of conventional instruments to date, including for example the development of photonic lanterns, complex aperiodic fiber Bragg gratings, complex beam combiners to enable long baseline interferometry, and laser frequency combs for high precision spectral calibration of spectrometers. Despite these successes, the facility implementation of photonic solutions in astronomical instrumentation is currently limited because of (1) low throughputs from coupling to fibers, coupling fibers to chips, propagation and bend losses, device losses, etc, (2) difficulties with scaling to large channel count devices needed for large bandwidths and high resolutions, and (3) efficient integration of photonics with detectors, to name a few. In this roadmap, we identify 24 areas that need further development. We outline the challenges and advances needed across those areas covering design tools, simulation capabilities, fabrication processes, the need for entirely new components, integration and hybridization and the characterization of devices. To realize these advances the astrophotonics community will have to work cooperatively with industrial partners who have more advanced manufacturing capabilities. With the advances described herein, multi-functional instruments will be realized leading to novel observing capabilities for both ground and space platforms.

We present high angular resolution imaging that reveals the MOA-2008-BLG-379L exoplanet host star using both adaptive optics with the Keck NIRC2 instrument and the Hubble Space Telescope. These observations reveal that the host star and planet masses are $M_{\rm host} = 0.434\pm 0.065 M_\odot$, and $m_p = 2.44 \pm 0.49 M_{\rm Jupiter}$. Their projected separation is $2.70\pm 0.42\,$AU, and they are located at a distance of $D_L = 3.44\pm 0.53\,$kpc toward the Galactic center. These results contribute to our effort to measure the masses of the exoplanet host stars in the \citet{suzuki16} statistical sample. When this program is completed, it will provide a measurement of the dependence of the planet occurrence rate on the mass and Galactocentric distance of the host stars. Our analysis also reveals the importance of including higher order effects, like microlensing parallax and planetary orbital motion even when the light curve data does not provide significant constraints on these effects. The inclusion of these effects is needed to obtain accurate estimates of the uncertainty of other microlensing parameters, such as the Einstein radius crossing time and planetary mass ratio. This will be particularly important for the exoplanet microlensing survey of the Roman Space Telescope, which will employ the same mass measurement method used in our analysis, using its own high angular resolution observations to determine the properties of most of its exoplanet host stars. Roman should also provide a large enough sample of exoplanetary microlens systems so that biases introduced by ignoring these higher order microlensing effects are likely to be significant.

LS V +44 17 is a persistent Be/X-ray binary (BeXRB) that displayed a bright, double-peaked period of X-ray activity in late 2022/early 2023. We present a radio monitoring campaign of this outburst using the Very Large Array. Radio emission was detected, but only during the second, X-ray brightest, peak, where the radio emission followed the rise and decay of the X-ray outburst. LS V +44 17 is therefore the third neutron star BeXRB with a radio counterpart. Similar to the other two systems (Swift J0243.6+6124 and 1A 0535+262), its X-ray and radio luminosity are correlated: we measure a power law slope $\beta = 1.25^{+0.64}_{-0.30}$ and a radio luminosity of $L_R = (1.6\pm0.2)\times10^{26}$ erg/s at a $0.5-10$ keV X-ray luminosity of $2\times10^{36}$ erg/s (i.e. $\sim 1\%$ $L_{\rm Edd}$). This correlation index is slightly steeper than measured for the other two sources, while its radio luminosity is higher. We discuss the origin of the radio emission, specifically in the context of jet launching. The enhanced radio brightness compared to the other two BeXRBs is the first evidence of scatter in the giant BeXRB outburst X-ray - radio correlation, similar to the scatter observed in sub-classes of low-mass X-ray binaries. While a universal explanation for such scatter is not known, we explore several options: we conclude that the three sources do not follow proposed scalings between jet power and neutron star spin or magnetic field, and instead briefly explore the effects that ambient stellar wind density may have on BeXRB jet luminosity.

The galaxy cluster Abell 746 (A746; $z$=0.214), featuring a double radio relic system, two isolated radio relics, a possible radio halo, disturbed V-shaped X-ray emission, and intricate galaxy distributions, is a unique and complex merging system. We present a weak-lensing analysis of A746 based on wide-field imaging data from Subaru/Hyper Suprime-Cam observations. The mass distribution is characterized by a main peak which coincides with the center of the X-ray emission. At this main peak, we detect two extensions toward the north and west, tracing the cluster galaxy and X-ray distributions. Despite the ongoing merger, our estimate of the A746 global mass $M_{500}=4.4\pm1.0\times10^{14}~M_{\odot}$ is consistent with the previous results from SZ and X-ray observations. We conclude that reconciling the distributions of mass, galaxies, and intracluster medium with the double radio relic system and other radio features remains challenging.

V1765 Cyg is a detached eclipsing binary containing a B0.5 supergiant and a B1 main-sequence star, with an orbital period of 13.37 d and an eccentricity of 0.315. The system shows apsidal motion and the supergiant exhibits strong stochastic variability. V1765 Cyg was observed by the Transiting Exoplanet Survey Satellite over four sectors. We analyse these data to obtain the first determinate light curve model for the system. To this we add published spectroscopic orbits to infer masses of 23 +/- 2 Msun and 11.9 +/- 0.7 Msun, and radii of 20.6 +/- 0.8 Rsun and 6.2 +/- 0.3 Rsun. These properties are in good agreement with theoretical predictions for a solar chemical composition and an age around 7 Myr. We also present two epochs of blue-optical spectroscopy that confirm the luminosity classification of the primary star and appear to show absorption lines from the secondary star. Extensive spectroscopy and further analysis of the system is recommended.

NASA's all-sky survey mission, the Transiting Exoplanet Survey Satellite (TESS) is specifically designed to detect transiting exoplanets orbiting bright stars. TESS has already identified about 400 transiting exoplanets, as well as approximately 6000 candidates awaiting validation. In this work, we present the outcomes of the project Validation of Transiting Exoplanets using Statistical Tools (VaTEST), an ongoing endeavor dedicated to validating and characterizing new exoplanets using statistical validation tools. We have validated eight potential super-Earths using a combination of ground-based telescope data, high-resolution imaging, and the statistical validation tool known as \texttt{TRICERATOPS}. These validated planets bear the designations: TOI-238b (1.61$^{+0.09} _{-0.10}$ R$_\oplus$), TOI-771b (1.42$^{+0.11} _{-0.09}$ R$_\oplus$), TOI-871b (1.66$^{+0.11} _{-0.11}$ R$_\oplus$), TOI-1467b (1.83$^{+0.16} _{-0.15}$ R$_\oplus$), TOI-1739b (1.69$^{+0.10} _{-0.08}$ R$_\oplus$), TOI-2068b (1.82$^{+0.16} _{-0.15}$ R$_\oplus$), TOI-4559b (1.42$^{+0.13} _{-0.11}$ R$_\oplus$), and TOI-5799b (1.62$^{+0.19} _{-0.13}$ R$_\oplus$). We also studied the synthetic transmission spectra of all eight validated planets in the HST and JWST band-pass using \texttt{PLATON} and \texttt{PandExo}. Among all these validated planets, six of them fall within the region known as 'keystone planets,' which makes them particularly interesting for study. Based on the location of TOI-771b and TOI-4559b below the radius valley we characterized them as likely super-Earths, though radial velocity mass measurements for these planets will provide more details about their characterization. It is noteworthy that planets within the size range investigated herein are absent from our own solar system, making their study crucial for gaining insights into the evolutionary stages between Earth and Neptune.

Confinement in $SU(N_{\rm DC})$ Yang-Mills theories is known to proceed through first-order phase transition. The wall velocity is bounded by $v_w \lesssim 10^{-6}$ due to the needed time for the substantial latent heat released during the phase transition to dissipate through Hubble expansion. Quarks much heavier than the confinement scale can be introduced without changing the confinement dynamics. After they freeze-out, heavy quarks are squeezed into pockets of the deconfined phase until they completely annihilate with anti-quarks. We calculate the dark baryon abundance surviving annihilation, due to bound-state formation occurring both in the bulk and - for the first time - at the boundary. We find that dark baryons can be dark matter with a mass up to $10^3~\rm TeV$. We study indirect and direct detection, CMB and BBN probes, assuming portals to Higgs and neutrinos.

Thermal relic dark matter below $\sim 10 \ \text{GeV}$ is excluded by cosmic microwave background data if its annihilation to visible particles is unsuppressed near the epoch of recombination. Usual model-building measures to avoid this bound involve kinematically suppressing the annihilation rate in the low-velocity limit, thereby yielding dim prospects for indirect detection signatures at late times. In this work, we investigate a class of cosmologically-viable sub-GeV thermal relics with late-time annihilation rates that are detectable with existing and proposed telescopes across a wide range of parameter space. We study a representative model of inelastic dark matter featuring a stable state $\chi_1$ and a slightly heavier excited state $\chi_2$ whose abundance is thermally depleted before recombination. Since the kinetic energy of dark matter in the Milky Way is much larger than it is during recombination, $\chi_1 \chi_1 \to \chi_2 \chi_2$ upscattering can efficiently regenerate a cosmologically long-lived Galactic population of $\chi_2$, whose subsequent coannihilations with $\chi_1$ give rise to observable gamma-rays in the $\sim 1 \ \text{MeV} - 100 \ \text{MeV}$ energy range. We find that proposed MeV gamma-ray telescopes, such as e-ASTROGAM, AMEGO, and MAST, would be sensitive to much of the thermal relic parameter space in this class of models and thereby enable both discovery and model discrimination in the event of a signal at accelerator or direct detection experiments.

If the dark matter in the Universe is made of $\mu$eV axion-like particles (ALPs), then a rich phenomenology can emerge in connection to their stimulated decay into two photons. We discuss the ALP stimulated decay induced by astrophysical beams of Galactic radio sources. Three signatures, made by two echoes and one collinear emission, are associated with the decay, and can be simultaneously detected, offering a unique opportunity for a clear ALP identification. We derive the formalism associated with such signatures starting from first principles, and providing the relevant equations to be applied to study the ALP phenomenology. We then focus on the case of Galactic pulsars as stimulating sources and derive forecasts for future observations.

Plasmas are beset with instabilities of all types, hydrodynamic, magneto-hydrodynamic, and electromagnetic. These instabilities are complex, occur over a large range of temporal and spatial scales, are most often unmanageable, and have seriously challenged our efforts at applications, even as they have shed light on the understanding of the physics of plasmas in the laboratory and astrophysical environments. A major reason for our limited success in their containment is the lack of direct experimental information on their origins and evolution, both temporal and spatial. In plasmas produced by high-intensity, short, and ultrashort pulse lasers, our knowledge of the instability stems from the (secondary) signals they generate e.g. scattering of electromagnetic waves in the form of Raman or Brillouin scattering. Rarely, if ever, has a direct measurement been made of the instantaneous evolution of the instabilities in plasmas. In this paper, we present direct measurements of the femtosecond evolution of the electromagnetic beam-driven instability that arises from the interaction of forward and return currents in an ultrahigh-intensity laser-produced plasma on a solid target.

We study how the nuclear symmetry energy slope ($L$) can affect the hadron-quark phase transition and neutron star properties. We show that the main physical quantities as the critical chemical potential and pressure are strongly influenced by the symmetry energy slope. In extreme cases, the total amount of deconfined quarks can reach up to 99$\%$ of the hybrid star mass.

We considered the generation of gravitational waves by the binary system associated with a wormhole. In the Newtonian limit, the gravitational potential of a wormhole requires the effective mass of the wormhole taking into account radial tension effects. This definition allows us to derive gravitational wave production in homogeneous and heterogeneous binary systems. Therefore, we studied gravitational waves generation by orbiting wormhole-wormhole and wormhole-black hole binary systems before coalescence. Cases involving negative mass require more careful handling. We also calculated the energy loss to gravitational radiation by a particle orbiting around the wormhole and by a particle moving straight through the wormhole mouth, respectively.

The reheating phase after inflation is one of the least observationally constrained epochs in the evolution of the Universe. The forthcoming gravitational wave observatories will enable us to constrain at least some of the non-standard scenarios. For example, if the radiation bath is produced via the perturbative inflaton decay that oscillates around a flat minimum of the potential of the form $V\propto\phi^{2n}$, with $n>2$. In such scenarios a part of the inflationary gravitational wave spectrum becomes blue tilted, making it observable, depending on the inflation energy scale and the reheating temperature. The degeneracy between the latter two parameters can be broken if dark matter in the Universe is produced by the freeze-in mechanism. The combination of the independent measurement of dark matter mass with gravitational wave observations makes it possible to constrain the reheating temperature and the energy density at the end of inflation.

For the strong gravitational wave model, an explicit transformation is obtained from a privileged coordinate system with a wave variable to a synchronous reference frame with separation of time and space variables. In a synchronous reference frame, a general form of the gravitational wave metric, solutions to the equations of trajectories for test particles in the Hamilton-Jacobi formalism, a solution to the eikonal equation for radiation, and a form of equations for the light cone of an observer in a gravitational wave were found. Using the obtained relations, the form the retarded time of radiation in the gravitational wave was found. The general relations obtained can be applied both in Einstein general theory of relativity and in modified theories of gravity. The obtained relations were applied in the work for an exact model of a gravitational wave in the Bianchi type VI universe based on an exact solution of Einstein vacuum equations.

We address the properties of extreme black holes by considering the Christodoulou-Ruffini/Hawking mass-energy formula. By simple geometrical arguments, we found that the mass/energy formula is satisfied by two meaningful extreme black holes where mass (m), charge (Q), and angular momentum/spin (L) are incommensurable with the black hole's irreducible mass (m_{ir}). These black holes have been studied in the Christodoulou diagram and their topology in E^3 has been investigated by differential geometry. We show that one of the analyzed Kerr-Newman black holes corresponds to the case where the Gaussian curvature becomes zero at the poles. In the second extreme black hole examined, the fundamental quantities m, Q, and L are linked to the irreducible mass by coefficients that depend solely on the golden ratio number -\phi_-. In this case, we show that if this extreme black hole satisfies the Pythagorean fundamental forms relation at the umbilic points, then both the "scale parameter" (corresponding to twice the irreducible mass) and the Gauss curvature of the surface at the poles are equal to the golden ratio numbers. For these two extreme black holes, we calculate the energy extractible by reversible transformations finding that, in percentage, the energy extractable from the latter black hole is higher than the former one.

In this paper we attempt to construct a regular gravastar model using the UV corrected framework of Loop Quantum Cosmology. We find that a stable gravastar model can be constructed with a number of unique features: (i) no thin shell approximation needs to be invoked to obtain solutions in the shell which can be considered to be of a finite thickness, (ii) the central singularity of a self gravitating object can be averted by a bounce mechanism, such that the interior density of the gravastar reaches a maximum critical density and cannot be raised further due to an operative repulsive force, (iii) the inherent isotropy of the effective fluid description does not prevent the formation of a stable gravastar and anisotropic pressures is not an essential requirement.