Papers for Thursday, Apr 13 2023

X-ray binaries, as bright local sources with short variability timescales for a wide range of accretion processes, represent ideal targets for high-resolution X-ray spectroscopy. In this chapter, we present a high-resolution X-ray spectral perspective on X-ray binaries, focusing on black holes and neutron stars. The majority of the chapter is devoted to observational and theoretical signatures of mass ejection via accretion disk winds: we discuss their appearance (including an overview of photoionization and thermodynamic processes that determine their visibility in X-ray spectra) and their life cycles (including efforts to constrain their time-dependent mass loss rates), and we provide a broad overview of the primary accretion disk wind driving mechanisms that have been considered in the literature: (1) radiation pressure, where radiation accelerates a wind by scattering off electrons or atoms in the disk or its atmosphere; (2) thermal driving, where Compton heating of the outer accretion disk causes gas thermal velocities to exceed the local escape speed; and (3) magnetohydrodynamic processes, where gas may be ejected from the disk via magnetic pressure gradients or magnetocentrifugal effects. We then turn to spectroscopic constraints on the geometry of accreting systems, from relativistically blurred emission lines to dipping sources, clumpy, structured stellar winds, and baryonic jets. We conclude with discussions of measurements of the interstellar medium and the potential of next-generation high-resolution X-ray spectroscopy for X-ray binaries.

We re-visit the calculation of mode oscillations in the ocean of a rotating neutron star, which may be excited during thermonuclear X-ray bursts. Our present theoretical understanding of ocean modes relies heavily on the traditional approximation, commonly employed in geophysics. The approximation elegantly decouples the radial and angular sectors of the perturbation problem by neglecting the vertical contribution from the Coriolis force. However, as the implicit assumptions underlying it are not as well understood as they ought to be, we examine the traditional approximation and discuss the associated mode solutions. The results demonstrate that, while the approximation may be appropriate in certain contexts, it may not be accurate for rapidly rotating neutron stars. In addition, using the shallow-water approximation, we show analytically how the solutions that resemble r-modes change their nature in neutron-star oceans to behave like gravity waves. We also outline a simple prescription for lifting Newtonian results in a shallow ocean to general relativity, making the result more realistic.

Cosmological simulations of structure formation predict that galaxy clusters continue to grow and evolve through ongoing mergers with group-scale systems. During these merging events, the ram pressure applied by the intracluster medium acts to strip the gas from the infalling groups, forming large tails of stripped gas, which eventually become part of the main cluster. In this work, we present a detailed analysis of our new deep Chandra observations of the NGC 4839 group falling into the nearby Coma cluster, providing a unique opportunity to explore the way galaxy clusters in the local universe continue to grow. Our analysis reveals a cold front feature at the leading head of the group, preceded by a bow shock of hot gas in front with a Mach number of $\sim\! 1.5$. The power spectrum of surface brightness fluctuations in the tail shows that the slope gets less steep as the distance from the leading head increases, changing from $-2.35_{-0.06}^{+0.07}$ at the inner part of the tail to $-1.37_{-0.07}^{+0.09}$ at the outermost part of the tail. These values are shallower than the slope of the Kolmogorov 2D power spectrum, indicating that thermal conduction is being suppressed throughout the tail, enabling long-lived small-scale turbulence, which would typically be washed out if thermal conduction was not inhibited. The characteristic amplitude of surface brightness fluctuations in the tail suggests a mild level of turbulence with a Mach number in the range of 0.1-0.5, agreeing with that found for the infalling group in Abell 2142.

We report the discovery of new lensed quasar candidates in the imaging data of HSC-SSP DR4. In addition to two catalogs of MILLIQUAS and AllWISEAGN, we preselected quasar sources using color cuts from the HSC ($grizy$) and unWISE ($W1+W2$) photometric data based on SDSS spectroscopic catalogs. We search for the presence of multiple point sources with similar color through the convolution of the Laplacian of the preselected quasar image cutouts with the Laplacian of the point spread function, resulting in a reduction of lens candidates from 1,652,329 to 234,424 (14.2%). After conducting a visual binary classification, we grade 6,199 (0.4%) potential lenses on a scale of 0 to 3, with 3 indicating a lens and 0 indicating a non-lens. Finally we obtain 162 lens candidates with a grade of $G_{\rm av}$, recovering 18 known lenses. By fitting the light distribution and removing the known contaminants, we discover that 73 new systems contain at least two point sources and a galaxy in between, including 10 possible quadruply lensed quasars. Spectroscopic or high-resolution imaging follow up on these newly discovered lensed quasar candidates will further allow their natures to be confirmed.

SOFIA/HAWC+ 154 $\mu$m Far-Infrared polarimetry observations of the well-studied edge-on galaxy NGC 891 are analyzed and compared to simple disk models with ordered (planar) and turbulent magnetic fields. The overall low magnitude and the narrow dispersion of fractional polarization observed in the disk require significant turbulence and a large number of turbulent decorrelation cells along the line-of-sight through the plane. Higher surface brightness regions along the major axis to either side of the nucleus show a further reduction in polarization and are consistent with a view tangent to a spiral feature in our disk models. The nucleus also has a similar low polarization, and this is inconsistent with our model spiral galaxy where the ordered magnetic field component would be nearly perpendicular to the line-of-sight through the nucleus on an edge-on view. A model with a barred spiral morphology with a magnetic field geometry derived from radio synchrotron observations of face-on barred spirals fits the data much better. There is clear evidence for a vertical field extending into the halo from one location in the disk coincident with a polarization null point seen in near-infrared polarimetry, probably due to a blowout caused by star formation. Although our observations were capable of detecting a vertical magnetic field geometry elsewhere in the halo, no clear signature was found. A reduced polarization due to a mix of planar and vertical fields in the dusty regions of the halo best explains our observations, but unusually significant turbulence cannot be ruled out.

Supernova neutrinos are invaluable signals that offer information about the interior of supernovae. Because a nearby supernova can occur at any time, preparing for future supernova neutrino observation is an urgent task. For the prompt analysis of supernova neutrinos, we have developed a new analysis code, ``Supernova Parameter Estimation Code based on Insight on Analytic Late-time Burst Light curve at Earth Neutrino Detector (SPECIAL BLEND)''. This code estimates the parameters of supernova based on an analytic model of supernova neutrinos from the proto-neutron star cooling phase. For easy availability to the community, this code is public and easily runs on web environments. SPECIAL BLEND can estimate the parameters better than the analysis pipeline we developed in the previous paper. By using SPECIAL BLEND, we can estimate the supernova parameters within $10\%$ precision up to $\sim 20\,{\rm kpc}$ and $\sim 60\,{\rm kpc}$ (Large Magellanic Cloud contained) with Super Kamiokande and Hyper Kamiokande, respectively.

The study of jet-inflated X-ray cavities provides a powerful insight into the energetics of hot galactic atmospheres and radio-mechanical AGN feedback. By estimating the volumes of X-ray cavities, the total energy and thus also the corresponding mechanical jet power required for their inflation can be derived. Properly estimating their total extent is, however, non-trivial, prone to biases, nearly impossible for poor-quality data, and so far has been done manually by scientists. We present a novel and automated machine-learning pipeline called Cavity Detection Tool (CADET), developed to detect and estimate the sizes of X-ray cavities from raw Chandra images. The pipeline consists of a convolutional neural network trained for producing pixel-wise cavity predictions and a DBSCAN clustering algorithm, which decomposes the predictions into individual cavities. The convolutional network was trained using mock observations of early-type galaxies simulated to resemble real noisy Chandra-like images. The network's performance has been tested on simulated data obtaining an average cavity volume error of 14 % at an 89 % true-positive rate. For simulated images without any X-ray cavities inserted, we obtain a 5 % false-positive rate. When applied to real Chandra images, the pipeline recovered 91 out of 100 previously known X-ray cavities in nearby early-type galaxies and all 14 cavities in chosen galaxy clusters. Besides that, the CADET pipeline discovered 8 new cavity pairs in atmospheres of early-type galaxies and galaxy clusters (IC4765, NGC533, NGC2300, NGC3091, NGC4073, NGC4125, NGC4472, NGC5129) and a number of potential cavity candidates.

The accretion of matter onto black holes and neutron stars often leads to the launching of outflows that can greatly affect the environments surrounding the compact object. In supermassive black holes, these outflows can even be powerful enough to dictate the evolution of the entire host galaxy, and yet, to date, we do not understand how these so-called accretion disk winds are launched - whether by radiation pressure, magnetic forces, thermal irradiation, or a combination thereof. An important means of studying disk winds produced near the central compact object is through X-ray absorption line spectroscopy, which allows us to probe outflow properties along a single line of sight, but usually provides little information about the global 3D disk wind structure that is vital for understanding the launching mechanism and total wind energy budget. Here, we study Hercules X-1, a unique, nearly edge-on X-ray binary with a warped accretion disk precessing with a period of about 35 days. This disk precession results in changing sightlines towards the neutron star, through the ionized outflow. We perform time-resolved X-ray spectroscopy over the precession phase and detect a strong decrease in the wind column density by three orders of magnitude as our sightline progressively samples the wind at greater heights above the accretion disk. The wind becomes clumpier as it rises upwards and expands away from the neutron star. Modelling the warped disk shape, we create a 2D map of wind properties. This unique measurement of the vertical structure of an accretion disk wind allows direct comparisons to 3D global simulations to reveal the outflow launching mechanism.

A Ribbon of enhanced energetic neutral atom (ENA) emissions was discovered by the Interstellar Boundary Explorer (IBEX) in 2009, redefining our understanding of the heliosphere boundaries and the physical processes occurring at the interstellar interface. The Ribbon signal is intertwined with that of a globally distributed flux (GDF) that spans the entire sky. To a certain extent, Ribbon separation methods enabled examining its evolution independent of the underlying GDF. Observations over a full solar cycle revealed the Ribbon's evolving nature, with intensity variations closely tracking those of the solar wind (SW) structure after a few years delay accounting for the SW-ENA recycling process. In this work, we examine the Ribbon structure, namely, its ENA fluxes, angular extent, width, and circularity properties for two years, 2009 and 2019, representative of the declining phases of two adjacent solar cycles. We find that, (i) the Ribbon ENA fluxes have recovered in the nose direction and south of it down to ~ 25{\deg} (for energies below 1.7 keV) and not at mid and high ecliptic latitudes; (ii) The Ribbon width exhibits significant variability as a function of azimuthal angle; (iii) Circularity analysis suggests that the 2019 Ribbon exhibits a statistically consistent radius with that in 2009. The Ribbon's partial recovery is aligned with the consensus of a heliosphere with its closest point being southward of the nose region. The large variability of the Ribbon width as a function of Azimuth in 2019 compared to 2009 is likely indicative of small-scale processes within the Ribbon.

A new determination of the temperature of the intergalactic medium over $3.9 \leq z \leq 4.3$ is presented. We applied the curvature method on a sample of 10 high resolution quasar spectra from the Ultraviolet and Visual Echelle Spectrograph on the VLT/ESO. We measured the temperature at mean density by determining the temperature at the characteristic overdensity, which is tight function of the absolute curvature irrespective of $\gamma$. Under the assumption of fiducial value of $\gamma = 1.4$, we determined the values of temperatures at mean density $T_{0} = 7893^{+1417}_{-1226}$ K and $T_{0} = 8153^{+1224}_{-993}$ K for redshift range of $3.9 \leq z \leq 4.1$ and $4.1 \leq z \leq 4.3$, respectively. Even though the results show no strong temperature evolution over the studied redshift range, our measurements are consistent with an intergalactic medium thermal history that includes a contribution from He II reionization.

The current expansion of the Universe has been observed to be accelerating, and the widely accepted spatially-flat concordance model of general relativistic cosmology attributes this phenomenon to a constant dark energy, a cosmological constant, which is measured to comprise about 70% of the total energy budget of the current Universe. However, observational discrepancies and theoretical puzzles have raised questions about this model, suggesting that alternative cosmological models with non-zero spatial curvature and/or dark energy dynamics might provide better explanations. To explore these possibilities, we have conducted a series of studies using standardized, lower-redshift observations to constrain six different cosmological models with varying degrees of flatness and dark energy dynamics. Through comparing these observations with theoretical predictions, we aim to deepen our understanding of the evolution of the Universe and shed new light on its mysteries. Our data provide consistent cosmological constraints across all six models, with some suggesting the possibility of mild dark energy dynamics and slight spatial curvature. However, these joint constraints do not rule out the possibility of dark energy being a cosmological constant and the spatial hypersurfaces being flat. Overall, our findings contribute to the ongoing efforts to refine our understanding of the Universe and its properties, and suggest that multiple cosmological models remain viable.

We present a probabilistic cross-correlation approach to estimate time delays in the context of reverberation mapping (RM) of Active Galactic Nuclei (AGN). We reformulate the traditional interpolated cross-correlation method as a statistically principled model that delivers a posterior distribution for the delay. The method employs Gaussian processes as a model for observed AGN light curves. We describe the mathematical formalism and demonstrate the new approach using both simulated light curves and available RM observations. The proposed method delivers a posterior distribution for the delay that accounts for observational noise and the non-uniform sampling of the light curves. This feature allow us to fully quantify its uncertainty and propagate it to subsequent calculations of dependent physical quantities, e.g., black hole masses. It delivers out-of-sample predictions, which enables us to subject it to model selection and it can calculate the joint posterior delay for more than two light curves. Because of the numerous advantages of our reformulation and the simplicity of its application, we anticipate that our method will find favour not only in the specialised community of RM, but in all fields where cross-correlation analysis is performed. We provide the algorithms and examples of their application as part of our Julia GPCC package.

We study the radio galaxies with known redshift detected by the Fermi satellite after 10 years of data (4FGL-DR2). We use a one-zone leptonic model to fit the quasi-simultaneous multiwavelength data of these radio galaxies and study the distributions of the derived physical parameter as a function of black hole mass and accretion disk luminosity. The main results are as follows. (1) We find that the jet kinetic power of most radio galaxies can be explained by the hybrid jet model based on ADAFs surrounding Kerr black holes. (2) After excluding the redshift, there is a significant correlation between the radiation jet power and the accretion disk luminosity, while the jet kinetic power is weakly correlated with the accretion disk luminosity. (3) We also find a significant correlation between inverse Compton luminosity and synchrotron luminosity. The slope of the correlation for radio galaxies is consistent with the synchrotron self-Compton (SSC) process. The result may suggest that the high-energy component of radio galaxies is dominated by the SSC process.

Despite advances in our understanding of low luminosity active galactic nuclei (LLAGNs), the fundamental details about the mechanisms of radiation and flare/outburst in hot accretion flow are still largely missing. We have systematically analyzed the archival Chandra and NuSTAR X-ray data of the nearby LLAGN M81*, whose $L_{\rm bol}\sim 10^{-5} L_{\rm Edd}$. Through a detailed study of X-ray light curve and spectral properties, we find that the X-ray continuum emission of the power-law shape more likely originates from inverse Compton scattering within the hot accretion flow. In contrast to Sgr A*, flares are rare in M81*. Low-amplitude variation can only be observed in soft X-ray band (amplitude usually $\lesssim 2$). Several simple models are tested, including sinusoidal-like and quasi-periodical. Based on a comparison of the dramatic differences of flare properties among Sgr A*, M31* and M81*, we find that, when the differences in both the accretion rate and the black hole mass are considered, the flares in LLAGNs can be understood universally in a magneto-hydrodynamical model.

About 25-50% of white dwarfs (WDs) show metal lines in their spectra. Among the widely accepted explanations for this effect is that the these WDs are accreting asteroids that are perhaps flung onto the WDs by a planet via resonance, for instance. A number of theoretical works have looked into the accretion of asteroids onto WDs and obtained a fair agreement with the observed accretion rate. However, it is solely a very recent study (referenced in this work) that has taken asteroid binarity into consideration, examining the scattering between an asteroid binary and planets and showing that a dissociation and ejection of the former might result and the effect on WD metal accretion is likely to be weak. Here, we investigate the close encounter between an asteroid binary and the central WD and consider how the binary's dissociation may affect the WD's accretion. We find that depending on the orbital and physical properties, the components may acquire orbits that are significantly different (even on the order of unity) from that of the parent binary. We assumed all the inner main belt asteroids are binaries and we studied their accretion onto the solar WD under the perturbation of the giant planets. We find that compared to the case without binaries, the components' accretion may be postponed (or put forward) by millions of years or more, as the objects may be taken out of (or driven deeper into) the resonance due to the sudden orbital change upon dissociation. However, the overall influence of binary dissociation on the accretion rate is not very significant.

Astronomy and Astrophysics is an observational science dealing with celestial objects. Aryabhatta Research Institute of Observational Sciences (ARIES) is one of the premier institutions in astronomy and astrophysics and has contributed significantly in this field. No doubt, India is a part of several mega-science projects in the domain of Astronomy and Astrophysics, such as the Thirty Meter Telescope (TMT); Square Kilometer Array (SKA) and Laser Interferometer Gravitational-wave Observatory (LIGO) projects. Growing engagement of India with mega-science projects has brought a positive impact on its science and technology landscape. A few such collaborations are mentioned to demonstrate that international cooperation are necessary in the field of Astrophysical sciences.

Using a sample of 361 nearby brown dwarfs, we have searched for 4.6$\mu$m variability indicative of large-scale rotational modulations or large-scale long-term changes on timescales of over 10 years. Our findings show no statistically significant variability in \textit{Spitzer} ch2 or \textit{WISE} W2 photometry. For \textit{Spitzer} the ch2 1$\sigma$ limits are $\sim$8 mmag for objects at 11.5 mag and $\sim$22 mmag for objects at 16 mag. This corresponds to no variability above 4.5$\%$ at 11.5 mag and 12.5$\%$ at 16 mag. We conclude that highly variable brown dwarfs, at least two previously published examples of which have been shown to have 4.6$\mu$m variability above 80 mmag, are very rare. While analyzing the data, we also developed a new technique for identifying brown dwarfs binary candidates in \textit{Spitzer} data. We find that known binaries have IRAC ch2 PRF (point response function) flux measurements that are consistently dimmer than aperture flux measurements. We have identified 59 objects that exhibit such PRF versus apertures flux differences and are thus excellent binary brown dwarf candidates.

Thermal inflation was proposed as a mechanism to dilute the density of cosmological moduli. Thermal inflation is driven by a complex scalar field possessing a large vacuum expectation value and a very flat potential, called a `flaton'. Such a model admits cosmic string solutions, and a network of such strings will inevitably form in the symmetry breaking phase transition at the end of the period of thermal inflation. We discuss the differences of these strings compared to the strings which form in the Abelian Higgs model. Specifically, we find that the upper bound on the symmetry breaking scale is parametrically lower than in the case of Abelian Higgs strings, and that the lower cutoff on the string loop distribution is determined by cusp annihilation rather than by gravitational radiation (for the value of the transition temperature proposed in the original work on thermal inflation).

We present light curves and flares from a seven day, multi-wavelength observational campaign of AU Mic, a young and active dM1e star with exoplanets and a debris disk. We report on 73 unique flares between the X-ray to optical data. We use high-time resolution NUV photometry and soft X-ray (SXR) data from XMM-Newton to study the empirical Neupert effect, which correlates the gradual and impulsive phase flaring emissions. We find that 65% (30 of 46) flares do not follow the Neupert effect, which is three times more excursions than seen in solar flares, and propose a four part Neupert effect classification (Neupert, Quasi-Neupert, Non-Neupert I & II) to explain the multi-wavelength responses. While the SXR emission generally lags behind the NUV as expected from the chromospheric evaporation flare models, the Neupert effect is more prevalent in larger, more impulsive flares. Preliminary flaring rate analysis with X-ray and U-band data suggests that previously estimated energy ratios hold for a collection of flares observed over the same time period, but not necessarily for an individual, multi-wavelength flare. These results imply that one model cannot explain all stellar flares and care should be taken when extrapolating between wavelength regimes. Future work will expand wavelength coverage using radio data to constrain the nonthermal empirical and theoretical Neupert effects to better refine models and bridge the gap between stellar and solar flare physics.

I build a toy model in the frame of the jittering jets explosion mechanism (JJEM) of core collapse supernovae (CCSNe) that incorporates both the stochastically varying angular momentum component of the material that the newly born neutron star (NS) accretes and the constant angular momentum component and show that the JJEM can account for the ~2.5-5Mo mass gap between NSs and black holes (BHs). The random component of the angular momentum results from pre-collapse core convection fluctuations that are amplified by post-collapse instabilities. The fixed angular momentum component results from pre-collapse core rotation. For slowly rotating pre-collapse cores the stochastic angular momentum fluctuations form intermittent accretion disks (or belts) around the NS with varying angular momentum axes in all directions. The intermittent accretion disk/belt launches jets in all directions that expel the core material in all directions early on, hence leaving a NS remnant. Rapidly rotating pre-collapse cores form an accretion disk with angular momentum axis that is about the same as the pre-collapse core rotation. The NS launches jets along this axis and hence the jets avoid the equatorial plane region. In-flowing core material continues to feed the central object from the equatorial plane increasing the NS mass to form a BH. The narrow transition from slow to rapid pre-collapse core rotation, i.e., from an efficient to inefficient jet feedback mechanism, accounts for the sparsely populated mass gap.

Given its large plate scale of 21" / pixel, analyses of data from the TESS space telescope must be wary of source confusion from blended light curves, which creates the potential to attribute observed photometric variability to the wrong astrophysical source. We explore the impact of light curve contamination on the detection of fast yellow pulsating supergiant (FYPS) stars as a case study to demonstrate the importance of confirming the source of detected signals in the TESS pixel data. While some of the FYPS signals have already been attributed to contamination from nearby eclipsing binaries, others are suggested to be intrinsic to the supergiant stars. In this work, we carry out a detailed analysis of the TESS pixel data to fit the source locations of the dominant signals reported for 17 FYPS stars with the Python package TESS_localize. We are able to reproduce the detections of these signals for 14 of these sources, obtaining consistent source locations for four. Three of these originate from contaminants, while the signal reported for BZ Tuc is likely a spurious frequency introduced to the light curve of this 127-day Cepheid by the data processing pipeline. Other signals are not significant enough to be localized with our methods, or have long periods that are difficult to analyze given other TESS systematics. Since no localizable signals hold up as intrinsic pulsation frequencies of the supergiant targets, we argue that unambiguous detection of pulsational variability should be obtained before FYPS are considered a new class of pulsator.

There is a lack of exoplanets with sizes similar to Neptune orbiting their host stars with periods $\lesssim 3$ days -- hence the name ``sub-Jovian/Neptune desert''. Recently, several exoplanets have been confirmed to reside in the desert transforming it into a ``savanna'' with several ``giraffe'' planets (such as LTT 9779 b and TOI-674 b). The most prominent scenarios put forward for the explanation of the formation of the desert are related to the stellar irradiation destroying the primary atmosphere of certain specific exoplanets. We aim to present three targets (LTT 9779 b, TOI-674 b and WASP-156 b) which, when observed at wide wavelength ranges in infrared (IR), could prove the presence of these processes, and therefore improve the theories of planetary formation/evolution. We simulate and analyse realistic light curves of the selected exoplanets with PLATO/NCAM and the three narrow-band filters of Ariel (VISPhot, FGS1 and FGS2) based on TESS observations of these targets. We improved the precision of the transit parameters of the three considered planets from the TESS data. We find that the combination of the three narrow-band filters of Ariel can yield inner precision of $\lesssim 1.1\%$ for the planetary radii. Data from the three telescopes together will span decades, allowing the monitoring of changes in the planetary atmosphere through radius measurements. The three selected ``giraffe'' planets can be golden targets for Ariel, whereby the loss of planetary mass due to stellar irradiation could be studied with high precision, multi-wavelength (spectro-)photometry.

Space debris, also known as "space junk," presents a significant challenge for all space exploration activities, including those involving human-onboard spacecraft such as SpaceX's Crew Dragon and the International Space Station. The amount of debris in space is rapidly increasing and poses a significant environmental concern. Various studies and research have been conducted on space debris capture mechanisms, including contact and contact-less capturing methods, in Earth's orbits. While advancements in technology, such as telecommunications, weather forecasting, high-speed internet, and GPS, have benefited society, their improper and unplanned usage has led to the creation of debris. The growing amount of debris poses a threat of collision with the International Space Station, shuttle, and high-value satellites, and is present in different parts of Earth's orbit, varying in size, shape, speed, and mass. As a result, capturing and removing space debris is a challenging task. This review article provides an overview of space debris statistics and specifications, and focuses on ongoing mitigation strategies, preventive measures, and statutory guidelines for removing and preventing debris creation, emphasizing the serious issue of space debris damage to space agencies and relevant companies.

In situ observations of interplanetary (IP) coronal mass ejections (ICMEs) and IP shocks are important to study as they are the main components of the solar activity. Hundreds of IP shocks have been detected by various space missions at different times and heliocentric distances. Some of these are followed by clearly identified drivers, while some others are not. In this study, we carry out a statistical analysis of the distributions of plasma and magnetic parameters of the IP shocks recorded at various distances to the Sun. We classify the shocks according to the heliocentric distance, namely from 0.29 to 0.99 AU (Helios-1/2); near 1 AU (Wind, ACE and STEREO-A/B); and from 1.35 to 5.4 AU (Ulysses). We also differentiate the IP shocks into two populations, those with a detected ICME and those without one. We find, as expected, that there are no significant differences in the results from spacecraft positioned at 1 AU. Moreover, the distributions of shock parameters, as well as the shock normal have no significant variations with the heliocentric distance. Additionally, we investigate how the number of shocks associated to stream-interaction regions (SIRs) increases with distance in proportion of ICME/shocks. From 1 to 5 AU, SIRs/ shock occurrence increases slightly from 21% to 34%, in contrast ICME/shocks occurrence decreases from 47% to 17%. We find also indication of an asymmetry induced by the Parker spiral for SIRs and none for ICMEs.

We present photometric and kinematic analysis of an intermediate age open cluster NGC 1027 using $UBV(RI)_c$ and Gaia Early Data Release 3 (EDR3) data. Structural and fundamental parameters such as cluster center, cluster extent, reddening, age and distance are estimated in this study. Cluster center is found about 2 arcmin away from the center reported earlier. Radius has been estimated to be about 8.00 arcmin(2.65pc). Using proper motion Gaia EDR3 data, membership probabilities has been derived for the stars in the region of cluster radius. We find mean proper motion of the cluster to be $\sim$(-0.84, 2.04) mas yr$^{-1}$ in (RA, DEC). We find 217 most probable (P$_\mu>$ 70\%) cluster members with mean parallax 0.892 $\pm$ 0.088 mas. Out of these, 160 members have counterparts in our optical observations. Few stars having P$_\mu>$ 70\% are found out of the cluster radius showing imprints of dynamical evolution. The color-color and color-magnitude diagrams for the cluster members found within 8.00 arcmin have been constructed using $UBV(RI)_c$ photometry and Gaia EDR3 data. This yields a reddening E($B$-$V$) $\sim$ 0.36 mag, age $\sim$ 130 Myr and distance $\sim$ 1.14 kpc. The mass function slope in the cluster region is $\Gamma$ $\sim$ -1.46 $\pm$ 0.15, which is similar to other Galactic open clusters. The dynamical study shows lack of faint stars in its inner region leading to mass segregation effect. A comparison of dynamical age with cluster age indicates that NGC 1027 is a dynamically relaxed cluster suggesting that mass segregation may be imprint of its dynamical relaxation.

The exoplanet population with orbital periods $P<100$ d around solar-type stars is dominated by super-Earths and sub-Neptunes. These planets are, however, missing in our Solar System, and the reason for that is unknown. Two theoretical scenarios invoke the role of Jupiter as the possible culprit: Jupiter may have acted as a dynamical barrier to the inward migration of sub-Neptunes from beyond the water iceline or, alternatively, may have reduced considerably the inward flux of material (pebbles) required to form super-Earths inside that iceline. Both scenarios predict an anti-correlation between the presence of small planets (SPs) and that of cold Jupiters (CJs) in exoplanetary systems. To test that prediction, we homogeneously analyzed the radial-velocity (RV) measurements of 38 Kepler and K2 transiting SP systems gathered over almost 10 years with the HARPS-N spectrograph, as well as publicly available RVs collected with other facilities. We detected five CJs in three systems, two in Kepler-68, two in Kepler-454, and a very eccentric one in K2-312. We derived an occurrence rate of $9.3^{+7.7}_{-2.9}\%$ for CJs with $0.3-13~M_{Jup}$ and 1-10 au, which is lower but still compatible at $1.3\sigma$ with that measured from RV surveys for solar-type stars, regardless of SP presence. This does not allow us to draw a firm conclusion about the predicted anti-correlation between SPs and CJs, which would require a considerably larger sample. Nevertheless, we found no evidence of previous claims of an excess of CJs in SP systems. As an important by-product of our analyses, we homogeneously determined the masses of 64 Kepler and K2 small planets, reaching a precision better than 5, 7.5 and 10$\sigma$ for 25, 13 and 8 planets, respectively. Finally, we release to the scientific community the 3661 HARPS-N radial velocities used in this work. [Abridged]

External FUV irradiation of protoplanetary disks has an important impact on their evolution and ability to form planets. However, nearby (<300 pc) star-forming regions lack sufficiently massive young stars, while the Trapezium Cluster and NGC 2024 have complicated star-formation histories and their O-type stars' intense radiation fields ($>10^4\,G_0$) destroy disks too quickly to study this process in detail. We study disk mass loss driven by intermediate (10 - 1000 $G_0$) FUV radiation fields in L1641 and L1647, where it is driven by more common A0 and B-type stars. Using the large (N=873) sample size offered by the Survey of Orion Disks with ALMA (SODA), we search for trends in the median disk dust mass with FUV field strength across the region as a whole and in two separate regions containing a large number of irradiated disks. For radiation fields between 1 - 100 $G_0$, the median disk mass in the most irradiated disks drops by a factor $\sim 2$ over the lifetime of the region, while the 95th percentile of disk masses drops by a factor 4 over this range. This effect is present in multiple populations of stars, and localized in space, to within 2 pc of ionizing stars. We fit an empirical irradiation - disk mass relation for the first time: $M_{\rm{dust,median}} = -1.3^{+0.14}_{-0.13} \log_{10}(F_{\rm{FUV}} / G_0) + 5.2^{+0.18}_{-0.19}$. This work demonstrates that even intermediate FUV radiation fields have a significant impact on the evolution of protoplanetary disks.

Kilonovae produced by mergers of binary neutron stars (BNSs) are important transient events to be detected by time domain surveys with the alerts from the ground-based gravitational wave detectors. The observational properties of these kilonovae depend on the physical processes involved in the merging processes and the equation of state (EOS) of neutron stars (NSs). In this paper, we investigate the dependence of kilonova luminosities on the parameters of BNS mergers, and estimate the distribution functions of kilonova peak luminosities (KLFs) at the u-, g-, r-, i-, y-, and z-bands as well as its dependence on the NS EOS, by adopting a comprehensive semi-analytical model for kilonovae (calibrated by the observations of GW170817), a population synthesis model for the cosmic BNSs, and the ejecta properties of BNS mergers predicted by numerical simulations. We find that the kilonova light curves depend on both the BNS properties and the NS EOS, and the KLFs at the considered bands are bimodal with the bright components mostly contributed by BNS mergers with total mass $\lesssim 3.2M_\odot$/$2.8M_\odot$ and fainter components mostly contributed by BNS mergers with total mass $\gtrsim 3.2M_\odot$/$2.8M_\odot$ by assuming a stiff/soft (DD2/SLy) EOS. The emission of the kilonovae in the KLF bright components are mostly due to the radiation from the wind ejecta by the remnant disks of BNS mergers, while the emission of the kilonovae in the KLF faint components are mostly due to the radiation from the dynamical ejecta by the BNS mergers.

We investigate phenomenologically the viability of fuzzy dark matter (FDM). We do this by confronting the predictions of the model, in particular the formation of a solitonic core at the centre of dark matter haloes, with a homogeneous and robust sample of high-resolution rotation curves from the ``LITTLE THINGS in 3D'' catalog. This comprises a collection of isolated, dark matter dominated dwarf-irregular galaxies that provides an optimal benchmark for cosmological studies. We use a statistical framework based on Markov-Chain Monte Carlo techniques that allows us to extract relevant parameters such as the axion mass, the mass of the solitonic core, the mass of the dark matter halo and its concentration parameter with a rather loose set of priors except for the implementation of a core-halo relation that is predicted by simulations. The results of the fits are used to perform various diagnostics on the predictions of the model. FDM provides an excellent fit to the rotation curves of the ``LITTLE THINGS in 3D'' catalog, with axion masses determined from different galaxies clustering around $m_a\approx2\times10^{-23}$ eV. However we find two major problems in our analysis. First, the data follow scaling relations of the properties of the core which are not consistent with the predictions of the soliton. This problem is particularly acute in the core radius - mass relation with a tension that, at face value, has a significance $\gtrsim5\sigma$. The second problem is related to the strong suppression of the linear power spectrum that is predicted by FDM for the axion mass preferred by the data. This can be constrained very conservatively by the galaxy counts in our sample, which leads to a tension exceeding again $5\sigma$. We estimate the effects of baryons in our analysis and discuss whether they could alleviate the tensions of the model with observations.

We investigate the early (t < 1 day) kilonova from the neutron star merger by deriving atomic opacities for all the elements from La to Ra (Z = 57 - 88) ionized to the states V - XI. The opacities at high temperatures for the elements with open f-shells (e.g., lanthanides) are exceptionally high, reaching kappa_{exp} ~ 10^4 cm2/g at lambda < 1000 A at T ~ 70,000 K, whereas, the opacities at the same temperature and wavelengths for the elements with the open d-, p-, and s-shells reach kappa_{exp} ~ 1 cm2/g, 0.1 cm2/g, and 0.01 cm2/g, respectively. Using the new opacity dataset, we derive the early kilonovae for various compositions and density structures expected for neutron star merger ejecta. The bolometric luminosity for the lanthanide-rich ejecta shows distinct signatures and is fainter than that for the lanthanide-free ejecta. The early luminosity is suppressed by the presence of a thin outer layer, agreeing with the results of Kasen et al. (2017) and Banerjee et al. (2020). The early brightness in Swift UVOT filters and in the optical g-, r-, i-, z-filters for a source at 100 Mpc are ~ 22 - 20 mag and ~ 21 - 19 mag, respectively, at t ~ 0.1 days. Such kilonovae are ideal targets for the upcoming UV satellites, such as ULTRASAT, UVEX, and DORADO, and the upcoming surveys, e.g., Vera Rubin Observatory. We suggest the gray opacities to reproduce the bolometric light curves with and without lanthanides are ~ 1 - 20 cm2/g and ~ 0.8 - 1 cm2/g.

Remnant radio galaxies represent the dying phase of radio-loud active galactic nuclei (AGN). Large samples of remnant radio galaxies are important for quantifying the radio galaxy life cycle. The remnants of radio-loud AGN can be identified in radio sky surveys based on their spectral index, or, complementary, through visual inspection based on their radio morphology. However, this is extremely time-consuming when applied to the new large and sensitive radio surveys. Here we aim to reduce the amount of visual inspection required to find AGN remnants based on their morphology, through supervised machine learning trained on an existing sample of remnant candidates. For a dataset of 4107 radio sources, with angular sizes larger than 60 arcsec, from the LOw Frequency ARray (LOFAR) Two-Metre Sky Survey second data release (LoTSS-DR2), we started with 151 radio sources that were visually classified as 'AGN remnant candidate'. We derived a wide range of morphological features for all radio sources from their corresponding Stokes-I images: from simple source catalogue-derived properties, to clustered Haralick-features, and self-organising map (SOM) derived morphological features. We trained a random forest classifier to separate the 'AGN remnant candidates' from the not yet inspected sources. The SOM-derived features and the total to peak flux ratio of a source are shown to be most salient to the classifier. We estimate that $31\pm5\%$ of sources with positive predictions from our classifier will be labelled 'AGN remnant candidates' upon visual inspection, while we estimate the upper bound of the $95\%$ confidence interval for 'AGN remnant candidates' in the negative predictions at $8\%$. Visual inspection of just the positive predictions reduces the number of radio sources requiring visual inspection by $73\%$.

Like the solar cycle, stellar activity cycles are also irregular. Observations reveal that rapidly rotating (young) Sun-like stars exhibit a high level of activity with no Maunder-like grand minima and rarely display smooth regular activity cycles. On the other hand, slowly rotating old stars like the Sun have low activity levels and smooth cycles with occasional grand minima. We, for the first time, try to model these observational trends using flux transport dynamo models. Following previous works, we build kinematic dynamo models of one solar mass star with different rotation rates. Differential rotation and meridional circulation are specified with a mean-field hydrodynamic model. We include stochastic fluctuations in the Babcock-Leighton source of the poloidal field to capture the inherent fluctuations in the stellar convection. Based on extensive simulations, we find that rapidly rotating stars produce highly irregular cycles with strong magnetic fields and rarely produce Maunder-like grand minima, whereas the slowly-rotating stars (with a rotation period of 10 days and longer) produce smooth cycles of weaker strength, long-term modulation in the amplitude, and occasional extended grand minima. The average duration and the frequency of grand minima increase with decreasing rotation rate. These results can be understood as the tendency of less supercritical dynamo in slower rotating stars to be more prone to produce extended grand minima

Context. Massive stars are generally believed to form in supersonic turbulent environment. However, recent observations have challenged this traditional view. High spatial and spectral resolution observations of the Orion Molecular Cloud and an infrared dark cloud G35.39 show a resolution-dependent turbulence, and that high-mass stars are forming exclusively in subsonic to transonic cores in those clouds. These studies demand a re-evaluation of the role of the turbulence in massive star formation. Aims. We aim to study the turbulence in a typical massive star-forming region G35.20-0.74 N with a sufficient spatial resolution to resolve the thermal Jeans length, and a spectral resolution to resolve the thermal linewidth. Methods. We use the ALMA dust continuum emission to resolve fragmentation, JVLA 1.2 cm continuum to trace ionized gas, and JVLA NH3 (1,1) to (7,7) inversion transition lines to trace linewidth, temperature, and dynamics. We fit those lines and remove line broadening due to channel width, thermal pressure, and velocity gradient to obtain a clean map of intrinsic turbulence. Results. We find that (1) the turbulence in G35.20 is overall supersonic, with mean and median Mach numbers 3.7 and 2.8, respectively. (2) Mach number decreases from 6-7 at 0.1 pc scale to <3 towards the central cores at 0.01 pc scale. (3) The central ALMA cores appear to be decoupled form the host filament, evident by an opposite velocity gradient and significantly reduced turbulence. Because of intense star formation activities in G35.20 (as compared to the relatively young and quiescent IRDC G35.39), the supersonic turbulence is likely replenished by protostellar outflows. G35.20 is, thus, representative of an evolved form of IRDC G35.39. More observations of a sample of IRDCs are highly demanded to further investigate the role of turbulence in initial conditions for the massive star formation.

Recently, the $\mu$-deformation-based approach to modeling dark matter, which exploits $\mu$-deformed thermodynamics, was extended to the study of galaxy halo density profile and of the rotation curves of a number of (dwarf or low brightness) galaxies. For that goal, $\mu$-deformed analogs of the Lane--Emden equation (LEE) have been proposed, and their solutions describing density profiles obtained. There are two seemingly different versions of $\mu$-deformed LEE which possess the same solution, and so we deal with their equivalence. From the latter property we derive new, rather unusual, $\mu$-deformed Heisenberg algebra (HA) for the position and momentum operators, and present the $\mu$-HA in few possible forms (each one at $\mu\to0$ recovers usual HA). The generalized uncertainty relation linked with the new $\mu$-HA is studied, along with its interesting implications including the appearance of the quadruple of both maximal and minimal lengths and momenta.

Quantifying changes in galaxies' interstellar medium (ISM) abundance after quenching star formation is an important aspect of galaxy evolution, but it is poorly constrained beyond the local universe. We characterise the dust-related properties in 548 quiescent galaxies observed at $0.1<z<0.6$ as part of the hCOSMOS spectroscopic survey. This is the largest sample of quiescent galaxies at intermediate redshifts, for which the co-evolution of dust, metals and stars have been estimated. We reveal the complex relations between the key markers of galaxies' dust life-cycles, such as specific dust mass ($M_{\rm dust}$/$M_{\rm \star}$), with gas-metallicity ($Z_{\rm gas}$), time since quenching ($t_{\rm quench}$), stellar age and size. We find morphology to be important factor of a large scatter ($\sim2$ orders of magnitude) in $M_{\rm dust}/M_{\rm \star}$. Through modelling the star formation histories of our objects, we derive a broad dynamical range of post-quenching timescales ($60\:\rm Myr<t_{\rm quench}<3.2\:\rm Gyr$). We find that $M_{\rm dust}/M_{\rm \star}$ is the highest in recently quenched systems ($t_{\rm quench}<500$ Myr), but its further evolution is non-monotonic as a consequence of diverse pathways for prolonged dust formation, or removal on various timescales. Our data are well reproduced by the SIMBA cosmological simulation and chemical models that include dust growth in the ISM. While this process is prevalent in dusty quiescent galaxies, $\sim15\%$ of objects show signs of external dust acquisition, most likely via minor mergers. Our results strongly suggest that prolonged dust production on a timescale $0.5-1\:\rm Gyr$ since quenching may be common in dusty quiescent galaxies at intermediate redshifts, even if their gas reservoirs are heavily exhausted (i.e., cold gas fraction $<1-5\%$).

Accurate quasar classifications and redshift measurements are increasingly important to precision cosmology experiments. Broad absorption line (BAL) features are present in 15-20\% of all quasars, and these features can introduce systematic redshift errors, and in extreme cases produce misclassifications. We quantitatively investigate the impact of BAL features on quasar classifications and redshift measurements with synthetic spectra that were designed to match observations by the Dark Energy Spectroscopic Instrument (DESI) survey. Over the course of five years, DESI aims to measure spectra for 40 million galaxies and quasars, including nearly three million quasars. Our synthetic quasar spectra match the signal-to-noise ratio and redshift distributions of the first year of DESI observations, and include the same synthetic quasar spectra both with and without BAL features. We demonstrate that masking the locations of the BAL features decreases the redshift errors by about 1\% and reduces the number of catastrophic redshift errors by about 80\%. We conclude that identifying and masking BAL troughs should be a standard part of the redshift determination step for DESI and other large-scale spectroscopic surveys of quasars.

Generalized symmetrons are models that have qualitatively similar features to the archetypal symmetron, but have barely been studied. In this article, we investigate for what parameter values the fifth forces induced by disformally coupling generalized symmetrons can provide an explanation for the difference between baryonic and lens masses of galaxies. While it is known that the standard symmetron struggles with providing an alternative source for the lensing otherwise attributed to particle dark matter, we show that some generalized symmetron models are more suitable for complying with existing constraints on disformal couplings. This motivates future studies of these only little explored models.

Mergers of Primordial Black Holes (PBHs) may contribute to the gravitational wave mergers detected by the LIGO-Virgo-KAGRA (LVK) Collaboration. We study the dynamics of PBH binaries dressed with dark matter (DM) spikes, for PBHs with extended mass functions. We analyze the impact of DM spikes on the orbital parameters of the PBH binaries formed in the early Universe and calculate their merger rates at the age of the Universe today. We consider two possible scenarios for the dynamics of the dressed binaries: assuming that either the DM spikes are completely evaporated from the binaries before merger or they remain static until the merger. Contrary to previous studies, we find that the presence of spikes may increase or decrease the present-day PBH merger rates, in some cases dramatically. Comparing with merger rates reported by the LVK Collaboration in the third Gravitational Wave Transient Catalog (GWTC-3), we derive approximate constraints on the fraction of Solar-mass PBHs in cold dark matter as $f_\mathrm{pbh}\leq \mathcal{O}(10^{-5} - 10^{-3})$, depending on the mass function. Our calculations are valid only for the idealized scenarios in which the DM spikes are either evaporated or static. However, they suggest that the impact of DM spikes on PBH merger rates may be more complicated than previously thought and motivate the development of a more general description of the merger dynamics, including feedback of the DM spikes in highly eccentric PBH binaries.

We present the peculiar case of PSZ2G091.83+26.11 at z=0.822. This cluster hosts a Mpc-scale radio halo and an elongated radio source, whose location with the respect to the intracluster medium (ICM) distribution and to the cluster centre is not consistent with a simple merger scenario. We use VLA data at 1-4 GHz to investigate the spectral and polarisation properties of the diffuse radio emission. We combine them with previously published data from LOFAR n the 120-168 MHz band, and from the uGMRT at 250-500 and 550-900 MHz. We also complement the radio data with Chandra X-ray observations to compare the thermal and non-thermal emission of the cluster. The elongated radio emission is visible up to 3.0 GHz and has an integrated spectral index of $-1.24\pm0.03$, with a steepening from $-0.89\pm0.03$ to $-1.39\pm0.03$. These values correspond to Mach numbers $\mathcal{M}_{\rm radio,int}=3.0\pm0.19$ and $\mathcal{M}_{\rm radio,inj}=2.48\pm0.15$. Chandra data reveals a surface brightness discontinuity at the location of the radio source, with a compression factor of $\mathcal{C}=2.22^{+0.39}_{-0.30}$ (i.e. $\mathcal{M}_{\rm Xray}=1.93^{+0.42}_{-0.32}$). We also find that the source is polarised at GHz frequencies. We estimate an intrinsic polarisation fraction of $\sim0.2$, a Rotation Measure of $\sim50~{\rm rad~m^{-2}}$ (including the Galactic contribution) and an external depolarisation of $\sim60~{\rm rad~m^{-2}}$. The $B$-vectors are aligned with the major axis of the source, suggesting magnetic field compression. Hence, we classify this source as a radio relic. We also find a linear/super-linear correlation between the non-thermal and thermal emission. We propose an off-axis merger and/or multiple merger events to explain the position and orientation of the relic. Given the properties of the radio relic, we speculate that PSZ2G091.83+26.11 is in a fairly young merger state.

Observations have shown that the effective temperature of hydrogen-free Wolf-Rayet (WR) stars is considerably lower than that of the standard model, which means that the radius of the observed H-free WR stars is several times larger than that estimated by the standard model. The envelope inflation structure (EIS) caused by the radiation luminosity being close to the Eddington luminosity in the iron opacity peak region of H-free WR stars may be the key to resolve the radius problem of H-free WR stars. We try to explain the H-free WR stars observed in the Milk Way (MW) and the Large Magellanic Cloud (LMC) by the He stars. Using the Modules for Experiments in Stellar Astrophysics code, we compute the evolution of He stars with and without MLT++ prescriptions and discuss their effects on the EIS. We have calculated the evolution of He stars using a new mass-loss rate formula and three different relative rotational velocity and compared our results with observations on Hertzsprung-Russell diagrams. The low luminosity (log$(L/L_{{\odot}})\leq5.2$) H-free WR stars in the MW and the LMC can be explained by the helium giant phase in low-mass He stars, the high $X_{C}$ and $X_{O}$ in WC stars can only evolve through low-mass He stars with a rapid rotation. High-mass He stars with the EIS can explain H-free WR stars with a luminosity exceeding $10^{5.7} L_{{\odot}}$ and an effective temperature above $10^{4.7}$ K in the MW. They can also explain H-free WR stars on the right-hand side of the He zero-age main sequence in the LMC. High-mass stars with the EIS evolve into WO stars at the final evolution stage, and the shorter lifetime fraction is consistent with the small number of observed WO stars.

Asymptotic Giant Branch (AGB) stars shed a significant amount of their mass in the form of a stellar wind, creating a vast circumstellar envelope (CSE). Owing to the ideal combination of relatively high densities and cool temperatures, CSEs serve as rich astrochemical laboratories. {While the chemical structure of AGB outflows has been modelled and analysed in detail for specific physical setups, there is a lack of understanding regarding the impact of changes in the physical environment on chemical abundances. A systematic sensitivity study is necessary to comprehend the nuances in the physical parameter space, given the complexity of the chemistry. This is crucial for estimating uncertainties associated with simulations and observations. In this work, we present the first sensitivity study of the impact of varying outflow densities and temperature profiles on the chemistry. With the use of a chemical kinetics model, we report on the uncertainty in abundances, given a specific uncertainty on the physical parameters. }Additionally, we analyse the molecular envelope extent of parent species and compare our findings to observational studies. Mapping the impact of differences in physical parameters throughout the CSE on the chemistry is a strong aid to observational studies.

Inertial modes have been observed on the Sun at low longitudinal wavenumbers. These modes probe the dynamics and structure of the solar convection zone down to the tachocline. While linear analysis allows the complex eigenfrequencies and eigenfunctions of these modes to be computed, it gives no information about their excitation nor about their amplitudes. We tested the hypothesis that solar inertial modes are stochastically excited by the turbulent motions entailed by convection. We have developed a theoretical formalism where the turbulent velocity fluctuations provide the mechanical work necessary to excite the modes. The modes are described by means of a 2D linear wave equation, relevant for the quasi-toroidal modes observed on the Sun, with a source term, under the beta plane approximation. Latitudinal differential rotation is included in the form of a parabolic profile that approximates the solar differential rotation at low and mid latitudes. We obtain synthetic power spectra for the wave's latitudinal velocity, longitudinal velocity, and radial vorticity, with azimuthal orders between 1 and 20. The synthetic power spectra contain the classical equatorial Rossby modes, as well as a rich spectrum of additional modes. The mode amplitudes are found to be of the same order of magnitude as observed on the Sun (~ 1 m/s). There is a qualitative transition between low and high azimuthal orders: the power spectra for m < 5 show modes that are clearly resolved in frequency space, while the power spectra for m > 5 display regions of excess power that consist of many overlapping modes. The general agreement between the predicted and observed inertial mode amplitudes supports the assumption of stochastic excitation by turbulent convection. Our work shows that the power spectra are not easily separable into individual modes, thus complicated the interpretation of the observations.

We utilize the probability distribution function (PDF) of convergence maps reconstructed from the Subaru Hyper Suprime-Cam (HSC) Y1 shear catalogue, in combination with the power spectrum, to measure the matter clustering amplitude $S_8=\sigma_8\sqrt{\Omega_m/0.3}$. The large-scale structure's statistical properties are incompletely described by the traditional two-point statistics, motivating our investigation of the PDF -- a complementary higher-order statistic. By defining the PDF over the standard deviation-normalized convergence map we are able to isolate the non-Gaussian information. We use tailored simulations to compress the data vector and construct a likelihood approximation. We mitigate the impact of survey and astrophysical systematics with cuts on smoothing scales, redshift bins, and data vectors. We find $S_8=0.852^{+0.086}_{-0.094}$ from the PDF alone and $S_8=0.827^{+0.033}_{-0.044}$ from the combination of PDF and power spectrum (68% CL). The PDF improves the power spectrum-only constraint by about 10%.

A massive naked singularity would be cloaked by accreted matter, and thus may appear to a distant observer as an opaque \mbox{(quasi-)}spherical surface of a fluid, not unlike that of a star or planet. We present here analytical solutions for levitating atmospheres around a wide class of spherically symmetric naked singularities. Such an atmosphere can be constructed in every spacetime which possesses a zero-gravity radius and which is a solution of a (modified-)gravity theory possessing the usual conservation laws for matter. Its density peaks at the zero-gravity radius and the atmospheric fluid is supported against infall onto the singularity by gravity alone. In an astrophysical context, an opaque atmosphere would be formed in a very short time by accretion of ambient matter onto the singularity -- in a millisecond for an X-ray binary, in a thousand seconds for a singularity traversing interstellar space, and a thousand years for a singularity that is the central engine of an AGN.

Primordial non-Gaussianity arising from inflationary models is a unique probe of non-trivial dynamics of the inflaton field and its interactions with other fields. Often when examining and constraining the scalar non-Gaussianity arising from inflation, certain templates are adopted for the scalar non-Gaussianity parameter $f_{_{\rm NL}}$, in classifying their behaviors in terms of wavenumbers. The current constraints from cosmic microwave background (CMB) on such templates of $f_{_{\rm NL}}$ are weak and provide rather large bounds on their amplitudes. In this work, we explore a different method of constraining $f_{_{\rm NL}}$ through their effect on the scalar power. We compute the correction to the scalar power due to $f_{_{\rm NL}}$ while accounting for its generic scale dependence. We then compute the angular power spectrum of CMB arising from such non-Gaussian corrections to explore possible imprints. We initially illustrate this method using the conventional templates of $f_{_{\rm NL}}$ such as local, equilateral and orthogonal types, with and without the running of the parameter. We further employ this method to an oscillatory form of $f_{_{\rm NL}}$ and lastly on a realistic model of inflation proposed by Starobinsky. Though this method does not improve much on the constraints on the first three templates of $f_{_{\rm NL}}$, it provides interesting insights on models that do not conform to these templates. We infer that the non-Gaussian correction to the spectrum can be sensitive to model parameters that are degenerate at the level of the original power spectrum. Hence, this method of computing indirect imprints of $f_{_{\rm NL}}$ on angular power spectrum of CMB provides a new avenue to explore primordial scalar non-Gaussianity and possibly constrain them effectively.

We investigated dusty and dust-free gas dynamics for a radiation-driven sub-pc scale outflow in an active galactic nucleus (AGN) associated with a supermassive black hole $10^7 M_\odot$ and bolometric luminosity $10^{44}$ erg s$^{-1}$ based on the two-dimensional radiation-hydrodynamic simulations. A radiation-driven ``lotus-like'' multi-shell outflow is launched from the inner part ($r \lesssim 0.04$ pc) of the geometrically thin disk, and it repeatedly and steadily produces shocks as mass accretion continues through the disk to the center. The shape of the dust sublimation radius is not spherical and depends on the angle ($\theta$) from the disk plane, reflecting the non-spherical radiation field and non-uniform dust-free gas. Moreover, we found that the sublimation radius of $\theta \sim 20$-$60$ deg varies on a timescale of several years. The ``inflow-induced outflow" contributes the obscuration of the nucleus in the sub-pc region. The column density of the dust-free gas is $N_{\rm H} \gtrsim 10^{22}$ cm$^{-2}$ for $r \lesssim 0.04$ pc. Gases near the disk plane ($\theta \lesssim 30$ degree) can be the origin of the Compton-thick component, which was suggested by the recent X-ray observations of AGNs. The dusty outflow from the sub-pc region can be also a source of material for the radiation-driven fountain for a larger scale.

The release of volatiles from comets is usually from direct sublimation of ices on the nucleus, but for very or hyper-active comets other sources have to be considered to account for the total production rates. In this work, we present new near-infrared imaging and spectroscopic observations of 46P/Wirtanen taken during its close approach to the Earth on 2018 December 19 with the MMIRS instrument at the MMT Observatory to search for signatures of icy or ice-rich grains in its inner coma that might explain its previously reported excess water production. The morphology of the images does not suggest any change in grain properties within the field of view, and the NIR spectra do not show the characteristic absorption features of water ice. Using a new MCMC-based implementation of the spectral modeling approach of Protopapa et al. (2018), we estimate the areal water ice fraction of the coma to be less than 0.6%. When combined with slit-corrected Afrho values for the J, H, and K bands and previously measured dust velocities for this comet, we estimate an icy grain production rate of less than 4.6 kg/s. This places a strict constraint on the water production rate from pure icy grains in the coma, and in turn we find that for the 2018-2019 apparition approximately 64% of 46P's surface was sublimating water near perihelion. WE then discuss 46P's modern properties within the context of other (formerly) hyper-active comets to understand how these complex objects evolve.

Carbon is an essential element for life but how much can be delivered to young planets is still an open question. The chemical characterization of planet-forming disks is a crucial step in our understanding of the diversity and habitability of exoplanets. Very low-mass stars ($<0.2~M_{\odot}$) are interesting targets because they host a rich population of terrestrial planets. Here we present the JWST detection of abundant hydrocarbons in the disk of a very low-mass star obtained as part of the MIRI mid-INfrared Disk Survey (MINDS). In addition to very strong and broad emission from C$_2$H$_2$ and its $^{13}$C$^{12}$CH$_2$ isotopologue, C$_4$H$_2$, benzene, and possibly CH$_4$ are identified, but water, PAH and silicate features are weak or absent. The lack of small silicate grains implies that we can look deep down into this disk. These detections testify to an active warm hydrocarbon chemistry with a high C/O ratio in the inner 0.1 au of this disk, perhaps due to destruction of carbonaceous grains. The exceptionally high C$_2$H$_2$/CO$_2$ and C$_2$H$_2$/H$_2$O column density ratios suggest that oxygen is locked up in icy pebbles and planetesimals outside the water iceline. This, in turn, will have significant consequences for the composition of forming exoplanets.

Recent high angular resolution ALMA observations have revealed numerous gaps in protoplanetary disks. A popular interpretation has been that planets open them. Most previous investigations of planet gap-opening have concentrated on viscous disks. Here, we carry out 2D (axisymmetric) global simulations of gap opening by a planet in a wind-launching non-ideal MHD disk with consistent thermochemistry. We find a strong concentration of poloidal magnetic flux in the planet-opened gap, where the gas dynamics are magnetically dominated. The magnetic field also drives a fast (nearly sonic) meridional gas circulation in the denser disk regions near the inner and outer edges of the gap, which may be observable through high-resolution molecular line observations. The gap is more ionized than its denser surrounding regions, with a better magnetic field-matter coupling. In particular, it has a much higher abundance of molecular ion HCO$^+$, consistent with ALMA observations of the well-studied AS 209 protoplanetary disk that has prominent gaps and fast meridional motions reaching the local sound speed. Finally, we provide fitting formulae for the ambipolar and Ohmic diffusivities as a function of the disk local density, which can be used for future 3D simulations of planet gap-opening in non-ideal MHD disks where thermochemistry is too computationally expensive to evolve self-consistently with the magneto-hydrodynamics.

We revisit the idea that the inflaton may have dissipated part of its energy into a thermal bath during inflation, considering monomial inflationary potentials and three different forms of dissipation rate. Using a numerical Fokker-Planck approach to describe the stochastic dynamics of inflationary fluctuations, we confront this scenario with current bounds on the spectrum of curvature fluctuations and primordial gravitational waves. We also obtain analytical approximations that outperform those frequently used in previous analyses. We show that only our numerical Fokker-Planck method is accurate, fast and precise enough to test these models against current data. We advocate its use in future studies of warm inflation. We also apply the stochastic inflation formalism to this scenario, finding that a commonly implemented large thermal correction to the primordial spectrum--that had been argued to become apparent with it--is actually not required. Improved bounds on the scalar spectral index will further constrain warm inflation in the near future.

The large-scale structure (LSS) of the Universe is comprised of galaxy filaments, tendrils, and voids. The majority of the Universe's volume is taken up by these voids, which exist as underdense, but not empty, regions. The galaxies found inside these voids are expected to be some of the most isolated objects in the Universe. This study, using the Galaxy and Mass Assembly (GAMA) and Galaxy Zoo surveys, aims to investigate basic physical properties and morphology of void galaxies versus field (filament and tendril) galaxies. We use void galaxies with stellar masses of $9.35 < log(M/M_\odot) < 11.25$, and this sample is split by identifying two redshift-limited regions, 0 < z < 0.075, and, $0.075 < z < 0.15$. To find comparable objects in the sample of field galaxies from GAMA and Galaxy Zoo, we identify "twins" of void galaxies as field galaxies within $\pm$0.05 dex and $\pm$0.15 dex of M and specific star formation rate. We determine the statistical significance of our results using the Kolmogorov-Smirnov (KS) test. We see that void galaxies, in contrast with field galaxies, seem to be disk-dominated and have predominantly round bulges (with > 50 percent of the Galaxy Zoo citizen scientists agreeing that bulges are present).

Recent work on metal-intermediate globular clusters (GCs) with [Fe/H]=$-1.5$ and $-0.75$ has illustrated the theoretical behavior of multiple populations in photometric diagrams obtained with the James Webb Space Telescope (JWST). These results are confirmed by observations of multiple populations among M-dwarfs of 47 Tucanae. Here, we explore the multiple populations in metal-poor GCs with [Fe/H]=$-$2.3. We take advantage of synthetic spectra and isochrones that account for the chemical composition of multiple populations to identify photometric diagrams that separate the distinct stellar populations of GCs. We derived high-precision photometry and proper motion for main-sequence stars in the metal-poor GC M 92 from JWST and Hubble Space Telescope (HST) images. We identified a first generation (1G) and two main groups of second-generation stars (2G$_{\rm A}$ and 2G$_{\rm B}$) and investigated their kinematics and chemical composition. We find isotropic motions with no differences among the distinct populations. The comparison between the observed colors of M 92 stars and the colors derived by synthetic spectra reveals that helium abundance of 2G$_{\rm A}$ and 2G$_{\rm B}$ stars are higher than that of the 1G by $\Delta Y \sim 0.01$ and $0.04$, respectively. The $m_{\rm F090W}$ vs. $m_{\rm F090W}-m_{\rm F277W}$ color-magnitude diagram shows that below the knee, MS stars exhibit a wide color broadening due to multiple populations. We constrain the amount of oxygen variation needed to reproduce the observed MS width, which is consistent with results on red-giant branch stars. We conclude that multiple populations with masses of $\sim$0.1-0.8$M_{\odot}$ share similar chemical compositions.

We demonstrate the potential of the state-of-the-art OpenAI GPT-4 large language model to engage in meaningful interactions with Astronomy papers using in-context prompting. To optimize for efficiency, we employ a distillation technique that effectively reduces the size of the original input paper by 50\%, while maintaining the paragraph structure and overall semantic integrity. We then explore the model's responses using a multi-document context (ten distilled documents). Our findings indicate that GPT-4 excels in the multi-document domain, providing detailed answers contextualized within the framework of related research findings. Our results showcase the potential of large language models for the astronomical community, offering a promising avenue for further exploration, particularly the possibility of utilizing the models for hypothesis generation.

The occurrence of singularities where spacetime curvature becomes infinite and geodesic evolution breaks down are inevitable events in classical general relativity (GR) unless one chooses an exotic matter violating weak energy condition. These singularities show up in various physical processes, such as the gravitational collapse, the birth of the universe in the standard cosmology as well as the classical solutions of the black hole spacetimes. In the last two decades, a rigorous understanding of the dynamics of quantum spacetime and the way it resolves singularities has been achieved in loop quantum cosmology (LQC) which applies the concepts and techniques of loop quantum gravity to the symmetry reduced cosmological spacetimes. Due to the fundamental discreteness of quantum geometry derived from the quantum theory, the big bang singularity has been robustly shown to be replaced by a big bounce. Strong curvature singularities intrinsic in the classical cosmology are generically resolved for a variety of cosmological spacetimes including anisotropic models and polarized Gowdy models. Using effective spacetime description the LQC universe also provides an ultra-violet complete description of the classical inflationary scenario as well as its alternatives such as the ekpyrotic and matter bounce scenarios. In this chapter we provide a summary of singularity resolution and its physical implications for various isotropic and anisotropic cosmological spacetimes in LQC and analyze robustness of results through variant models originating from different quantization prescriptions.

The confluence of major theoretical, experimental, and observational advances are providing a unique perspective on the equation of state of dense neutron-rich matter -- particularly its symmetry energy -- and its imprint on the mass-radius relation for neutron stars. In this contribution we organize these developments in an equation of state density ladder. Of particular relevance to this discussion is the impact of the various rungs on the equation of state and the identification of possible discrepancies among the various methods. A preliminary analysis identifies a possible tension between laboratory measurements and gravitational-wave detections that could indicate the emergence of a phase transition in the stellar core.

5D warped extra dimension models with multiple 3-branes can naturally realize multiple hierarchical mass scales which are ubiquitous in physics beyond the Standard Model. We discuss cosmological consequences of such multi-brane models with stabilized radions. It is confirmed that for temperatures below the scale of the IR brane at the end of the extra dimension, we recover the ordinary expansion of the Universe, with the Hubble expansion rate determined by sum of the physical energy densities on all 3-branes where they are localized. In addition, we explore the cosmology for temperatures above the scales of the intermediate and IR branes where the Universe is described by a spacetime with the 3-branes replaced by an event horizon. As the temperature of the Universe cools down, phase transitions are expected to take place, and the intermediate and IR branes come out from behind the event horizon. The Goldberger-Wise mechanism for radion stabilization has a well-known problem of having a supercooled phase transition, which typically does not get completed in time. This problem is even more severe when an intermediate brane is introduced, whose scale is well above TeV, as the corresponding Hubble rate is much larger. We circumvent the problem by employing an alternative mechanism for radion stabilization with dark Yang-Mills fields, which prevents a long supercooling epoch, but still allows the strong first order phase transitions. As a result, the phase transitions in our multi-brane Universe predict a stochastic gravitational wave background with a unique multi-peak signature, which is within the sensitivity reach of future space-based gravitational wave observers. We also show that there are $N-1$ radions for an $N$ 3-brane set-up, unlike a recent claim that there exists only one radion.

We revisit the Affleck-Dine leptogenesis via the $L H_u$ flat direction with a light slepton field. Although the light slepton field is favored in low-energy SUSY phenomenologies, such as the muon $g-2$ anomaly and bino-slepton coannihilation, it may cause a problem in the Affleck-Dine leptogenesis: it may create an unwanted charge-breaking vacuum in the Affleck-Dine field potential so that the Affleck-Dine field is trapped during the course of leptogenesis. We investigate the conditions under which such an unwanted vacuum exists and clarify that both thermal and quantum corrections are important for the (temporal) disappearance of the charge-breaking minimum. We also confirm that if the charge-breaking vacuum disappears due to the thermal or quantum correction, the correct baryon asymmetry can be produced while avoiding the cosmological gravitino problem.

We review a recently proposed definition of complexity of the structure of self--gravitating fluids \cite{ch1}, and the criterium to define the simplest mode of their evolution. We analyze the origin of these concepts and their possible applications in the study of gravitation collapse. We start by considering the static spherically symmetric case, extending next the study to static axially symmetric case. Afterward we consider the non--static spherically symmetric case. Two possible modes of evolution are proposed to be the simplest one. One is the homologous conditio,, however, as was shown later on, it may be useful to relax this last condition to enlarge the set of possible solutions, by adopting the so-called quasi-homologous condition. As another example of symmetry, we consider fluids endowed with hyperbolical symmetry. Exact solutions for static fluid distributions satisfying the condition of minimal complexity are presented.. An extension of the complexity factor to the vacuum solutions of the Einstein equations represented by the Bondi metric is discussed. A complexity hierarchy is established in this case, ranging from the Minkowski spacetime (the simplest one) to gravitationally radiating systems (the most complex). Finally we propose a list of questions which, we believe, deserve to be treated in the future

We present a review on the self-resonant dark matter scenarios where multiple components of dark matter give rise to a resonant condition in the $u$-channel diagrams for their comparable masses. In this case, there is no need of lighter mediators for enhancing the self-scattering and annihilation cross sections for dark matter. We discuss the velocity-dependent self-scattering for the small-scale problems, the relic density of self-resonant dark matter, and the observable signatures in indirect and detection experiments.