Papers for Tuesday, Dec 13 2022

Ultra-diffuse galaxies (UDGs) are attractive candidates to probe cosmological models and test theories of galaxy formation at low masses; however, they are difficult to detect because of their low surface brightness. In the Local Group (LG) a handful of UDGs have been found to date, most of which are satellites of the Milky Way and M31, and only two are isolated galaxies. It is unclear whether so few UDGs are expected. We address this by studying the population of UDGs formed in hydrodynamic constrained simulations of the LG from the HESTIA suite. For a LG with mass $M_{\rm LG}\!\left(<2.5\, {\rm Mpc}\right)=8\times10^{12}{\rm M_\odot}$, we predict that there are $12\pm3$ UDGs (68 per cent confidence) with stellar masses $10^6 \leq M_\ast\, /\, {\rm M_\odot} < 10^9$, and effective radii $R_{\rm e} \geq 1.5\, {\rm kpc}$, in the field of the LG, of which $2^{+2}_{-1}$ (68 per cent confidence) are detectable in the footprint of the Sloan Digital Sky Survey (SDSS). Accounting for survey incompleteness, we find that up to 82, 90, and 100 per cent of all UDGs in the LG field would be observable in a future all-sky survey with a depth similar to the SDSS, the Dark Energy Survey, or the Legacy Survey of Space and Time, respectively. Our results suggest that there is a population of UDGs in the LG awaiting discovery.

We utilize $\sim17000$ bright Luminous Red Galaxies (LRGs) from the novel Dark Energy Spectroscopic Instrument Survey Validation spectroscopic sample, leveraging its deep ($\sim2.5$ hour/galaxy exposure time) spectra to characterize the contribution of recently quenched galaxies to the massive galaxy population at $0.4<z<1.3$. We use Prospector to infer non-parametric star formation histories and identify a significant population of post-starburst galaxies that have joined the quiescent population within the past $\sim1$ Gyr. The highest redshift subset (277 at $z>1$) of our sample of recently quenched galaxies represents the largest spectroscopic sample of post-starburst galaxies at that epoch. At $0.4<z<0.8$, we measure the number density of quiescent LRGs, finding that recently quenched galaxies constitute a growing fraction of the massive galaxy population with increasing lookback time. Finally, we quantify the importance of this population amongst massive ($\mathrm{log}(M_\star/M_\odot)>11.2$) LRGs by measuring the fraction of stellar mass each galaxy formed in the Gyr before observation, $f_{\mathrm{1 Gyr}}$. Although galaxies with $f_{\mathrm{1 Gyr}}>0.1$ are rare at $z\sim0.4$ ($\lesssim 0.5\%$ of the population), by $z\sim0.8$ they constitute $\sim3\%$ of massive galaxies. Relaxing this threshold, we find that galaxies with $f_\mathrm{1 Gyr}>5\%$ constitute $\sim10\%$ of the massive galaxy population at $z\sim0.8$. We also identify a small but significant sample of galaxies at $z=1.1-1.3$ that formed with $f_{\mathrm{1 Gyr}}>50\%$, implying that they may be analogues to high-redshift quiescent galaxies that formed on similar timescales. Future analysis of this unprecedented sample promises to illuminate the physical mechanisms that drive the quenching of massive galaxies after cosmic noon.

Early JWST spectroscopic campaigns have confirmed the presence of strong [O III] line-emitting galaxies in the redshift interval $7<z<9$. Although deduced earlier from Spitzer photometry as indicative of young stellar populations, some studies suggested the relevant photometric excesses attributed to [O III] emission could, in part, be due to Balmer breaks arising from older stars. We demonstrate that this is likely the case by exploiting medium-band near-infrared JWST photometry in the Hubble Ultra Deep Field. We locate a sample of 6 galaxies with redshifts 8.2$<z<$8.6 for which the relevant medium-band filters enables us to separate the contributions of [O III] emission and a Balmer break, thereby breaking earlier degeneracies of interpretation. The technique is particularly valuable since it provides photometric redshifts whose precision, $\Delta\,z\simeq\,\pm0.08$, approaches that of spectroscopic campaigns now underway with JWST. Although some sources are young, a third of our sample have prominent Balmer breaks consistent with stellar ages of $\geq$150 Myr. Our results indicate that even intense [O III] emitters experienced episodes of earlier star formation to $z\sim$10 and beyond, as is now being independently deduced from direct detection of the progenitors of similar systems.

Using early-release data from the JWST, Mowla et al. and Claeyssens et al. recently measured various properties for gravitationally lensed compact sources (`sparkles') around the `Sparkler' galaxy at a redshift of 1.378 (a look-back time of 9.1 Gyr). Here, we focus on the Mowla et al. as they were able to break the age-metallicity degeneracy and derive independent ages, metallicities and extinctions for each source. They identified 5 metal-rich, old GC candidates (with formation ages up to $\sim$13 Gyr). We examine the age--metallicity relation (AMR) for the GC candidates and other Sparkler compact sources. The Sparkler galaxy, which has a current estimated stellar mass of 10$^9$ M$_{\odot}$, is compared to the Large Magellanic Cloud (LMC), the disrupted dwarf galaxy Gaia--Enceladus and the Milky Way (MW). The Sparkler galaxy appears to have undergone very rapid chemical enrichment in the first few hundred Myr after formation, with its GC candidates similar to those of the MW's metal-rich subpopulation. We also compare the Sparkler to theoretical AMRs and formation ages from the E-MOSAICS simulation, finding the early formation age of its GCs to be in some tension with these predictions for MW-like galaxies. The metallicity of the Sparkler's star forming regions are more akin to a galaxy of stellar mass $\ge$ 10$^{10.5}$ M$_{\odot}$, i.e. at the top end of the expected mass growth over 9.1 Gyr of cosmic time. We conclude that the Sparkler galaxy may represent a progenitor of a MW-like galaxy, even including the ongoing accretion of a satellite galaxy.

TDEs have been observed in the optical and UV for more than a decade but the underlying emission mechanism still remains a puzzle. It has been suggested that viewing angle effects could potentially explain their large photometric and spectroscopic diversity. Polarization is indeed sensitive to the viewing angle and the first polarimetry studies of TDEs are now available, calling for a theoretical interpretation. In this study, we model the continuum polarization levels of TDEs using the radiative transfer code POSSIS and the collision-induced outflow (CIO) TDE emission scenario where unbound shocked gas originating from a debris stream intersection point offset from the black hole, reprocesses the hard emission from the accretion flow into UV and optical bands. We explore two different cases of peak mass fallback rates M'p (~3 and ~0.3 Msol/yr) while varying the following geometrical parameters: the distance R_int from the black hole (BH) to the intersection point, the radius of the photosphere around the BH R_ph, on the surface of which the photons are generated, and the opening angle Deltheta (anisotropic emission). For the high mass fallback rate case, we find for every viewing angle polarization levels below one (P<1%) and P<0.5% for 10/12 simulations. The absolute value of polarization reaches its maximum (P_max) for equatorial viewing angles. For the low mass fallback rate case, the maximum value predicted is P~8.8% and P_max is reached for intermediate viewing angles. We find that the polarization depends strongly on i) the optical depths at the central regions set by the different M'p values and ii) the viewing angle. Finally, by comparing our model predictions to polarization observations of a few TDEs, we attempt to constrain their observed viewing angles and we show that multi-epoch polarimetric observations can become a key factor in constraining the viewing angle of TDEs.

Entropy is an advantageous diagnostics to study the thermodynamic history of the intracluster plasma of galaxy clusters. We present the entropy profile of the Abell 2244 galaxy cluster derived both exclusively using X-ray data from the low-background Swift XRT telescope and also using Planck y data. The entropy profile derivation using X-rays only is robust at least to the virial radius because the cluster brightness is large compared to the X-ray background at low energies, temperature is strongly bounded by the lack of cluster X-ray photons at energies kT>3 keV, and the XRT background is low, stable and understood. In the observed solid angle, about one quadrant, the entropy radial profile deviates from a power-law at the virial radius, mainly because of a sharp drop of the cluster temperature. This bending of the entropy profile is confirmed when X-ray spectral information is replaced by the Compton map. Clumping and non-thermal pressure support are insufficient to restore a power law entropy profile because they are bound to be small by: a) the agreement between mass estimates from different tracers (gas and galaxies), b) the agreement between entropy profile determinations based on combinations of observables with different sensitivities and systematics, and c) the low value of clumping as estimated using the azimuthal scatter and the gas fraction. Based on numerical simulations, ion-electron equilibration is also insufficient to restore a linear entropy profile. Therefore, the bending of the entropy profiles seems to be robustly derived and witnesses the teoretically-predicted decrease in the inflow through the virial boundary.

Supergiant Fast X-ray Transients (SFXT) are High Mass X-ray Binaries displaying X-ray outbursts reaching peak luminosities of 10$^{38}$ erg/s and spend most of their life in more quiescent states with luminosities as low as 10$^{32}$-10$^{33}$ erg/s. The main goal of our comprehensive and uniform analysis of the SFXT Swift triggers is to provide tools to predict whether a transient which has no known X-ray counterpart may be an SFXT candidate. These tools can be exploited for the development of future missions exploring the variable X-ray sky through large FoV instruments. We examined all available data on outbursts of SFXTs that triggered the Swift/BAT collected between 2005-08-30 and 2014-12-31, in particular those for which broad-band data, including the Swift/XRT ones, are also available. We processed all BAT and XRT data uniformly with the Swift Burst Analyser to produce spectral evolution dependent flux light curves for each outburst. The BAT data allowed us to infer useful diagnostics to set SFXT triggers apart from the general GRB population, showing that SFXTs give rise uniquely to image triggers and are simultaneously very long, faint, and `soft' hard-X-ray transients. The BAT data alone can discriminate very well the SFXTs from other fast transients such as anomalous X-ray pulsars and soft gamma repeaters. However, to distinguish SFXTs from, for instance, accreting millisecond X-ray pulsars and jetted tidal disruption events, the XRT data collected around the time of the BAT triggers are decisive. The XRT observations of 35/52 SFXT BAT triggers show that in the soft X-ray energy band, SFXTs display a decay in flux from the peak of the outburst of at least 3 orders of magnitude within a day and rarely undergo large re-brightening episodes, favouring in most cases a rapid decay down to the quiescent level within 3-5 days (at most). [Abridged]

Gamma-ray bursts (GRBs) are the most luminous transients in the universe and are utilized as probes of early stars, gravitational wave counterparts, and collisionless shock physics. In spite of studies on polarimetry of GRBs in individual wavelengths that characterized intriguing properties of prompt emission and afterglow, no coordinated multi-wavelength measurements have yet been performed. Here, we report the first coordinated simultaneous polarimetry in the optical and radio bands for the afterglow associated with the typical long GRB 191221B. Our observations successfully caught the radio emission, which is not affected by synchrotron self-absorption, and show that the emission is depolarized in the radio band compared to the optical one. Our simultaneous polarization angle measurement and temporal polarization monitoring indicate the existence of cool electrons that increase the estimate of jet kinetic energy by a factor of $>$ 4 for this GRB afterglow. Further coordinated multi-wavelength polarimetric campaigns would improve our understanding of the total jet energies and magnetic field configurations in the emission regions of various types of GRBs, which are required to comprehend the mass scales of their progenitor systems and the physics of collisionless shocks.

Line-intensity mapping (IM) experiments seek to perform statistical measurements of large-scale structure with spectral lines such as 21cm, CO, and Lyman-$\alpha$ (Ly$\alpha$). A challenge in these observations is to ensure that astrophysical foregrounds, such as galactic synchrotron emission in 21cm measurements, are properly removed. One method that has the potential to reduce foreground contamination is to cross correlate with a galaxy survey that overlaps with the IM volume. However, telescopes sensitive to high-redshift galaxies typically have small field of views (FOVs) compared to IM surveys. Thus, a galaxy survey for cross correlation would necessarily consist of pencil beams which sparsely fill the IM volume. In this paper, we develop the formalism to forecast the sensitivity of cross correlations between IM experiments and pencil-beam galaxy surveys. We find that a random distribution of pencil beams leads to very similar overall sensitivity as a lattice spaced across the IM survey and derive a simple formula for random configurations that agrees with the Fisher matrix formalism. We explore examples of combining high-redshift James Webb Space Telescope (JWST) observations with both a SPHEREx-like Ly$\alpha$ IM survey and a 21cm experiment based on the Hydrogen Epoch of Reionization Array (HERA). We find that the JWST-SPHEREx case is promising, leading to a total signal-to-noise of ${\sim}7$ after 100 total hours of JWST (at $z=7$). We find that HERA is not well-suited for this approach owing to its drift-scan strategy, but that a similar experiment that can integrate down on one field could be.

Understanding the nucleosynthetic origin of nitrogen and the evolution of the N/O ratio in the interstellar medium is crucial for a comprehensive picture of galaxy chemical evolution at high-redshift because most observational metallicity (O/H) estimates are implicitly dependent on the N/O ratio. The observed N/O at high-redshift shows an overall constancy with O/H, albeit with a large scatter. We show that these heretofore unexplained features can be explained by the pre-supernova wind yields from rotating massive stars (M$\gtrsim 10 \, \mathrm{M}_\odot$, $v/v_{\rm{crit}} \gtrsim 0.4$). Our models naturally produce the observed N/O plateau, as well as the scatter at low O/H. We find the scatter to arise from varying star formation efficiency. However, the models that have supernovae dominated yields produce a poor fit to the observed N/O at low O/H. This peculiar abundance pattern at low O/H suggests that dwarf galaxies are most likely to be devoid of SNe yields and are primarily enriched by pre-supernova wind abundances.

In this paper, we generalize the Weinberg's procedure to determine the comoving curvature perturbation $\cal R$ to non-attractor inflationary regimes. We show that both modes of $\cal R$ are related to a symmetry of the perturbative equations in the Newtonian gauge. As a byproduct, we clarify that adiabaticity does not generally imply constancy of $\cal R$, not even in the $k\rightarrow 0$ limit. Applying this knowledge to the separate Universe approach, we find that correlators of $\delta N$ {\it do not} generically correspond to comoving curvature perturbations correlators, even at the linear level, but rather to correlators of curvature perturbations at uniform density, at least at linear level. Thus, $\delta N$ formalism does not capture information about decaying (for slow-roll) or growing (beyond slow-roll) modes of $\cal R$. The latter being the only interesting mode for models of inflation related to primordial black holes formation.

We study the effect of deep heating (stellar irradiation deposited at the ~10^5 Pa level, as well as in the shallow region at the ~10^3 level) in hot-exoplanet atmospheres. Unlike with shallow heating (only), the atmosphere with deep heating exhibits a single equilibrium state, characterized by repeated generation of giant cyclonic vortices that move away westward from the point of emergence. The generation is accompanied by a burst of heightened turbulence activity, leading to the production of small-scale structures and large-scale mixing of temperature on a timescale of ~3 planetary rotations.

Approximately one-third of existing $\gamma$-ray sources identified by the $\textit{Fermi Gamma-Ray Space Telescope}$ are considered to be unassociated, with no known counterpart at other frequencies/wavelengths. These sources have been the subject of intense scrutiny and observational effort during the observatory's mission lifetime, and here we present a method of leveraging existing radio catalogs to examine these sources without the need for specific dedicated observations, which can be costly and complex. Via the inclusion of many sensitive low-frequency catalogs we specifically target steep spectrum sources such as pulsars. This work has found steep-spectrum radio sources contained inside 591 $\textit{Fermi}$ unassociated fields, with at least 21 of them being notable for having pulsar-like $\gamma$-ray properties as well. We also identify a number of other fields of interest based on various radio and $\gamma$-ray selections.

The proposed next generation Event Horizon Telescope (ngEHT) concept envisions the imaging of various astronomical sources on scales of microarcseconds in unprecedented detail with at least two orders of magnitude improvement in the image dynamic ranges by extending the Event Horizon Telescope (EHT). A key technical component of ngEHT is the utilization of large aperture telescopes to anchor the entire array, allowing the connection of less sensitive stations through highly sensitive fringe detections to form a dense network across the planet. Here, we introduce two projects for planned next generation large radio telescopes in the 2030s on the Chajnantor Plateau in the Atacama desert in northern Chile, the Large Submillimeter Telescope (LST) and the Atacama Large Aperture Submillimeter Telescope (AtLAST). Both are designed to have a 50-meter diameter and operate at the planned ngEHT frequency bands of 86, 230 and 345\,GHz. A large aperture of 50\,m that is co-located with two existing EHT stations, the Atacama Large Millimeter/Submillimeter Array (ALMA) and the Atacama Pathfinder Experiment (APEX) Telescope in the excellent observing site of the Chajnantor Plateau, will offer excellent capabilities for highly sensitive, multi-frequency, and time-agile millimeter very long baseline interferometry (VLBI) observations with accurate data calibration relevant to key science cases of ngEHT. In addition to ngEHT, its unique location in Chile will substantially improve angular resolutions of the planned Next Generation Very Large Array in North America or any future global millimeter VLBI arrays if combined. LST and AtLAST will be a key element enabling transformative science cases with next-generation millimeter/submillimeter VLBI arrays.

We develop a method to calibrate u-band photometry based on the observed color of blue galactic halo stars. The galactic halo stars belong to an old stellar population of the Milky Way and have relatively low metallicity. The "blue tip" of the halo population -- the main sequence turn-off (MSTO) stars -- is known to have a relatively uniform intrinsic edge u-g color with only slow spatial variation. In SDSS data, the observed variation is correlated with galactic latitude, which we attribute to contamination by higher-metallicity disk stars and fit with an empirical curve. This curve can then be used to calibrate u-band imaging if g-band imaging of matching depth is available. Our approach can be applied to single-field observations at $|b| > 30^\circ$, and removes the need for standard star observations or overlap with calibrated u-band imaging. We include in our method the calibration of g-band data with ATLAS-Refcat2. We test our approach on stars in KiDS DR 4, ATLAS DR 4, and DECam imaging from the NOIRLab Source Catalog (NSC DR2), and compare our calibration with SDSS. For this process, we use synthetic magnitudes to derive the color equations between these datasets, in order to improve zero-point accuracy. We find an improvement for all datasets, reaching a zero-point precision of 0.016 mag for KiDS (compared to the original 0.033 mag), 0.020 mag for ATLAS (originally 0.027 mag), and 0.016 mag for DECam (originally 0.041 mag). Thus, this method alone reaches the goal of 0.02 mag photometric precision in u-band for the Rubin Observatory's Legacy Survey of Space and Time (LSST).

A growing avenue for determining the prevalence of life beyond Earth is to search for "technosignatures" from extraterrestrial intelligences/agents. Technosignatures require significant energy to be visible across interstellar space and thus intentional signals might be concentrated in frequency, in time, or in space, to be found in mutually obvious places. Therefore, it could be advantageous to search for technosignatures in parts of parameter space that are mutually-derivable to an observer on Earth and a distant transmitter. In this work, we used the L-band (1.1-1.9 GHz) receiver on the Robert C. Byrd Green Bank Telescope (GBT) to perform the first technosignature search pre-synchronized with exoplanet transits, covering 12 Kepler systems. We used the Breakthrough Listen turboSETI pipeline to flag narrowband hits ($\sim$3 Hz) using a maximum drift rate of $\pm$614.4 Hz/s and a signal-to-noise threshold of 5 - the pipeline returned $\sim 3.4 \times 10^5$ apparently-localized features. Visual inspection by a team of citizen scientists ruled out 99.6% of them. Further analysis found 2 signals-of-interest that warrant follow-up, but no technosignatures. If the signals-of-interest are not re-detected in future work, it will imply that the 12 targets in the search are not producing transit-aligned signals from 1.1-1.9 GHz with transmitter powers $>$60 times that of the former Arecibo radar. This search debuts a range of innovative technosignature techniques: citizen science vetting of potential signals-of-interest, a sensitivity-aware search out to extremely high drift rates, a more flexible method of analyzing on-off cadences, and an extremely low signal-to-noise threshold.

We explore model performance for the Alfven Wave Solar atmosphere Model (AWSoM) with near-real-time (NRT) synoptic maps of the photospheric vector magnetic field. These maps, produced by assimilating data from the Helioseismic Magnetic Imager (HMI) onboard the Solar Dynamics Observatory (SDO), use a different method developed at the National Solar Observatory (NSO) to provide a near contemporaneous source of data to drive numerical models. Here, we apply these NSO-HMI-NRT maps to simulate three Carrington rotations (CRs): 2107-2108 (centered on 2011/03/07 20:12 CME event), 2123 (integer CR) and 2218--2219 (centered on 2019/07/2 solar eclipse), which together cover a wide range of activity level for solar cycle 24. We show simulation results, which reproduce both extreme ultraviolet emission (EUV) from the low corona while simultaneously matching in situ observations at 1 au as well as quantify the total unsigned open magnetic flux from these maps.

We study the non-equilibrium dynamics of Axion-like particles (ALP) coupled to Standard Model degrees of freedom in thermal equilibrium. The Quantum Master Equation (QME) for the (ALP) reduced density matrix is derived to leading order in the coupling of the (ALP) to the thermal bath, but to \emph{all} orders of the bath couplings to degrees of freedom within or beyond the Standard Model other than the (ALP). The (QME) describes the damped oscillation dynamics of an initial misaligned (ALP) condensate, thermalization with the bath, decoherence and entropy production within a unifying framework. The (ALP) energy density $\mathcal{E}(t)$ features two components: a ``cold'' component from the misaligned condensate and a ``hot'' component from thermalization with the bath, with $\mathcal{E}(t)= \mathcal{E}_{c}\,e^{-\gamma(T)\,t}+\mathcal{E}_{h}(1-e^{-\gamma(T)\,t})$ thus providing a ``mixed dark matter'' scenario. Relaxation of the (ALP) condensate, thermalization, decoherence and entropy production occur on similar time scales. An explicit example with (ALP)-photon coupling, valid post recombination yields a relaxation rate $\gamma(T)$ with a substantial enhancement from thermal emission and absorption. A misaligned condensate is decaying at least since recombination and on the same time scale thermalizing with the cosmic microwave background (CMB). Possible consequences for birefringence of the (CMB) and (ALP) contribution to the effective number of ultrarelativistic species and galaxy formation are discussed.

In this study, we report on a detailed single pulse analysis of the radio emission from the pulsar J1401$-$6357 (B1358$-$63) based on data observed with the ultrawideband low-frequency receiver on the Parkes radio telescope. In addition to a weak leading component, the integrated pulse profile features a single-humped structure with a slight asymmetry. The frequency evolution of the pulse profile is studied. Well-defined nulls, with an estimated nulling fraction greater than 2\%, are present across the whole frequency band. No emission is detected with significance above 3$\sigma$ in the average pulse profile integrated over all null pulses. Using fluctuation spectral analysis, we reveal the existence of temporal-dependent subpulse drifting in this pulsar for the first time. A clear double-peaked feature is present at exactly the alias border across the whole frequency band, which suggests that the apparent drift sense changes during the observation. Our observations provide further confirmation that the phenomena of pulse nulling and subpulse drifting are independent of observing frequency, which suggest that they invoke changes on the global magnetospheric scale.

We used the Immersion GRating Infrared Spectrometer (IGRINS) to determine fundamental parameters for 61 K- and M-type young stellar objects (YSOs) located in the Ophiuchus and Upper Scorpius star-forming regions. We employed synthetic spectra and a Markov chain Monte Carlo approach to fit specific K-band spectral regions and determine the photospheric temperature ($T_{\rm eff}$), surface gravity ($\log$ g), magnetic field strength (B), projected rotational velocity ($v \sin i$), and K-band veiling ($r_{\rm K}$). We determined B for $\sim$46% of our sample. Stellar parameters were compared to the results from Taurus-Auriga and the TW Hydrae Association (TWA) presented in Paper I of this series. We classified all the YSOs in the IGRINS survey with infrared spectral indices from 2MASS and WISE photometry between 2 and 24~$\mu$m. We found that Class II YSOs typically have lower $\log$ g and $v\sin i$, similar B, and higher K-band veiling than their Class III counterparts. Additionally, we determined the stellar parameters for a sample of K and M field stars also observed with IGRINS. We have identified intrinsic similarities and differences at different evolutionary stages with our homogeneous determination of stellar parameters in the IGRINS YSO Survey. Considering $\log$ g as a proxy for age, we found that the Ophiuchus and Taurus samples have a similar age. We also find that Upper Scorpius and TWA YSOs have similar ages, and are more evolved than Ophiuchus/Taurus YSOs.

Mutual events (MEs) are eclipses and occultations among planetary natural satellites. Most of the time, eclipses and occultations occur separately. However, the same satellite pair will exhibit an eclipse and an occultation quasi-simultaneously under particular orbital configurations. This kind of rare event is termed as a quasi-simultaneous mutual event (QSME). During the 2021 campaign of mutual events of jovian satellites, we observed a QSME between Europa and Ganymede. The present study aims to describe and study the event in detail. We observed the QSME with a CCD camera attached to a 300-mm telescope at the Hong Kong Space Museum Sai Kung iObservatory. We obtained the combined flux of Europa and Ganymede from aperture photometry. A geometric model was developed to explain the light curve observed. Our results are compared with theoretical predictions (O-C). We found that our simple geometric model can explain the QSME fairly accurately, and the QSME light curve is a superposition of the light curves of an eclipse and an occultation. Notably, the observed flux drops are within 2.6% of the theoretical predictions. The size of the event central time O-Cs ranges from -14.4 to 43.2 s. Both O-Cs of flux drop and timing are comparable to other studies adopting more complicated models. Given the event rarity, model simplicity and accuracy, we encourage more observations and analysis on QSMEs to improve Solar System ephemerides.

We study the statistical property of fast radio bursts (FRBs) based on a selected sample of 190 one-off FRBs in the first CHIME/FRB catalog. Three power law models are used in the analysis, and we find the cumulative distribution functions of energy can be well fitted by bent power law and thresholded power law models. And the distribution functions of fluctuations of energy well follow the Tsallis $q$-Gaussian distribution. The $q$ values in the Tsallis $q$-Gaussian distribution are constant with small fluctuations for different temporal scale intervals, indicating a scale-invariant structure of the bursts. The earthquakes and soft gamma repeaters show similar properties, which are consistent with the predictions of self-organized criticality systems.

About 900 fast radio burst (FRB) sources have been detected till now, among which 29 FRBs are found to burst out repeatedly. Although a firm connection between at least some FRBs and magnetars has been established, the trigger mechanism and radiation process in these enigma phenomena are still highly controversial. In this study, we build a sample of 16 repeating FRBs from which at least five bursts have been detected, including the most active four repeaters of FRBs 20121102A, 20180916B, 20190520B and 20201124A. Various key parameters of their bursts are collected from the literature, which include the trigger time, pulse width, dispersion measure (DM), Faraday rotation measure (RM), bandwidth, waiting time, peak flux and fluence. The distribution and time evolution of these parameters are investigated. Potential correlations between various parameter pairs are also extensively explored. The behaviors of different repeaters are then compared. It is found that the DM of FRB 20121102A seems to increase continuously on a long timescale. While the $DM$ of most repeaters varies in a narrow range of $\pm 3$ cm$^{-3}$ pc, FRB 20190520B is found to have a large variation range of $\pm 12$ cm$^{-3}$ pc. A linear correlation between $DM$ and the star formation rate is established. $RM$ evolves with time in a much more chaos behavior in different repeaters. Generally, the waiting time shows a bimodal distribution in each source. The implications of these features to the underlying physics are discussed.

We present the long term analysis of GS 1826-238, a neutron star X-ray binary known as the "Clocked Burster", using data from NuSTAR StrayCats. StrayCats, a catalogue of NuSTAR stray light data, contains data from bright, off-axis X-ray sources that have not been focused by the NuSTAR optics. We obtained stray light observations of the source from 2014-2021, reduced and analyzed the data using nustar-gen-utils Python tools, demonstrating the transition of source from the "island" atoll state to a "banana" branch. We also present the lightcurve analysis of Type I X-Ray bursts from the Clocked Burster and show that the bursts from the banana/soft state are systematically shorter in durations than those from the island/hard state and have a higher burst fluence. From our analysis, we note an increase in mass accretion rate of the source, and a decrease in burst frequency with the transition.

The standard technique for very low-frequency gravitational wave detection is mainly based on searching for a specific spatial correlation in the variation of the times of arrival of the radio pulses emitted by millisecond pulsars with respect to a timing model. This spatial correlation, which in the case of the gravitational wave background must have the form described by the Hellings and Downs function, has not yet been observed. Therefore, despite the numerous hints of a common red noise in the timing residuals of many millisecond pulsars compatible with that expected for the gravitational wave background, its detection has not yet been achieved. By now, the reason is not completely clear and, from some recent works, the urgency to adopt new detection techniques, possibly complementary to the standard one, is emerging clearly. Of course, this demand also applies to the detection of continuous gravitational waves emitted by supermassive black hole binaries populating the Universe. In the latter case, important information could, in principle, emerge from the millisecond pulsars considered individually in a single-pulsar search of continuous GWs. In this context, the surfing effect can then be exploited, helping to select the best pulsars to carry out such analysis. This paper aims to clarify when the surfing effect occurs and describe it exhaustively. A possible application to the case of the supermassive black hole binary candidate PKS 2131-021 and millisecond pulsar J2145-0750 is also analyzed.

X-ray spectroscopy of galactic black hole binaries serve as a powerful tool to gain an overall understanding of the system. Not only can the properties of the accretion disk be studied in detail, the fundamental properties of the black hole can also be inferred. The pursuit of these objectives also leads to an indirect validation of general relativity in strong field limit. In this work we carry out a comprehensive spectral analysis of the galactic X-ray binary MAXI J1631-479 using data from NICER and NuSTAR observatories. We trace the evolution of the accretion disk properties such as density, ionization, Fe abundance, etc as the source transitions from a disk dominated soft state to a power law dominated hard intermediate state. We provide strong constrains on the spin of the black hole and the inclination of the inner disk. We also use the soft state NICER observations to constrain the black hole mass using distance estimates from optical observations.

The splashback radius of a dark matter halo, which corresponds to the first apocenter radius reached by infalling matter and substructures, has been detected around galaxy clusters using a multitude of observational methods, including weak lensing measurements. In this manuscript, we present how the splashback feature in the halo density profile affects galaxy cluster masses derived through weak lensing measurements if it is not accounted for. We find that the splashback radius has an increasingly large effect on group-sized halos towards $M_{200m} \sim 10^{13.5} \mathrm{M_\odot}$. Depending on the model and the radial scale used, the cluster/group masses can be biased low by more than 0.1 dex. This bias, in turn, would result in a slightly lower $\Omega_m$ value if propagated into a cluster cosmology analysis. The splashback effect with group-sized dark matter halos may become important to consider, given the increasingly stringent cosmological constraints coming from optical wide-field surveys.

Water is crucial for the emergence and evolution of life on Earth. Recent studies of the water content in early forming planetary systems similar to our own show that water is an abundant and ubiquitous molecule, initially synthesized on the surfaces of tiny interstellar dust grains by the hydrogenation of frozen oxygen. Water then enters a cycle of sublimation/freezing throughout the successive phases of planetary system formation, namely, hot corinos and protoplanetary disks, eventually to be incorporated into planets, asteroids, and comets. The amount of heavy water measured on Earth and in early forming planetary systems suggests that a substantial fraction of terrestrial water was inherited from the very first phases of the Solar System formation and is 4.5 billion years old.

Observations of neutron stars may be used to study aspects of extremely dense matter, specifically a possibility of phase transitions to exotic states, such as de-confined quarks. We present a novel data analysis method for detecting signatures of dense-matter phase transitions in sets of mass-radius measurements, and study its sensitivity with respect to the size of observational errors and the number of observations. The method is based on machine learning anomaly detection coupled with normalizing flows technique: the algorithm trained on samples of astrophysical observations featuring no phase transition signatures interprets a phase transition sample as an ''anomaly''. For the sake of this study, we focus on dense-matter equations of state leading to detached branches of mass-radius sequences (strong phase transitions), use an astrophysically-informed neutron-star mass function, and various magnitudes of observational errors and sample sizes. The method is shown to reliably detect cases of mass-radius relations with phase transition signatures, while increasing its sensitivity with decreasing measurement errors and increasing number of observations. We discuss marginal cases, when the phase transition mass is located near the edges of the mass function range. Evaluated on the current state-of-art selection of real measurements of electromagnetic and gravitational-wave observations, the method gives inconclusive results, which we interpret as due to small available sample size, large observational errors and complex systematics.

Context. The globular cluster (GC) system of Circinus galaxy has not been probed previously partly because of the location of the galaxy at - 3.8$^\circ$ Galactic latitude which suffers severely from interstellar extinction, stellar crowding, and Galactic foreground contamination. However, the deep near-infrared (NIR) photometry by the VISTA Variables in the Via L\'actea Extended Survey (VVVX) in combination with the precise astrometry of Gaia EDR3 allow us to map GCs in this region. Aims. Our long-term goal is to study and characterise the distributions of GCs and Ultra-compact dwarfs of Circinus galaxy which is the nearest Seyfert II galaxy. Here we conduct the first pilot search for GCs in this galaxy. Methods. We use NIR VVVX photometry in combination with Gaia EDR3 astrometric features such as astrometric excess noise and BP/RP excess factor to build the first homogeneous catalogue of GCs in Circinus galaxy. A robust combination of selection criteria allows us to effectively clean interlopers from our sample. Results. We report the detection of$\sim$ 70 GC candidates in this galaxy at a 3 $\sigma$ confidence level. They show a bimodal colour distribution with the blue peak at (G-Ks)$_0$ = 0.985$\pm$0.127 mag with a dispersion of 0.211$\pm$0.091 mag and the red peak at (G-Ks)$_0$ = 1.625$\pm$0.177 mag with a dispersion of 0.482$\pm$0.114 mag. A GC specific frequency (S$_N$) of 1.3$\pm$0.2 was derived for the galaxy, and we estimated a total population of 120$\pm$40 GCs. Based on the projected radial distribution it appears that Circinus has a different distribution of GC candidates than MW and M31. Conclusions. We demonstrate that Circinus galaxy hosts a sizeable number of cluster candidates. This result is the first leap towards understanding the evolution of old stellar clusters in this galaxy.

Lobster eye telescopes are ideal monitors to detect X-ray transients, because they could observe celestial objects over a wide field of view in X-ray band. However, images obtained by lobster eye telescopes are modified by their unique point spread functions, making it hard to design a high efficiency target detection algorithm. In this paper, we integrate several machine learning algorithms to build a target detection framework for data obtained by lobster eye telescopes. Our framework would firstly generate two 2D images with different pixel scales according to positions of photons on the detector. Then an algorithm based on morphological operations and two neural networks would be used to detect candidates of celestial objects with different flux from these 2D images. At last, a random forest algorithm will be used to pick up final detection results from candidates obtained by previous steps. Tested with simulated data of the Wide-field X-ray Telescope onboard the Einstein Probe, our detection framework could achieve over 94% purity and over 90% completeness for targets with flux more than 3 mCrab (9.6 * 10-11 erg/cm2/s) and more than 94% purity and moderate completeness for targets with lower flux at acceptable time cost. The framework proposed in this paper could be used as references for data processing methods developed for other lobster eye X-ray telescopes.

Primordial black holes (PBHs) can be both candidates of dark matter and progenitors of binary black holes (BBHs) detected by the LIGO-Virgo-KAGRA collaboration. Since PBHs could form in the very early Universe through the gravitational collapse of primordial density perturbations, the population of BBHs detected by gravitational waves encodes much information on primordial curvature perturbation. In this work, we take a reliable and systematic approach to reconstruct the power spectrum of the primordial curvature perturbation from GWTC-3, under the hierarchical Bayesian inference framework, by accounting for the measurement uncertainties and selection effects. In addition to just considering the single PBH population model, we also report the results considering the multi-population model, i.e., the mixed PBH and astrophysical black hole binaries model. We find that the maximum amplitude of the reconstructed power spectrum of primordial curvature perturbation can be $\sim2.5\times10^{-2}$ at $\mathcal{O}(10^{5})~\rm Mpc^{-1}$ scales, which is consistent with the PBH formation scenario from inflation at small scales.

The chemical reservoir within protoplanetary disks has a direct impact on planetary compositions and the potential for life. A long-lived carbon-and nitrogen-rich chemistry at cold temperatures (<=50K) is observed within cold and evolved planet-forming disks. This is evidenced by bright emission from small organic radicals in 1-10 Myr aged systems that would otherwise have frozen out onto grains within 1 Myr. We explain how the chemistry of a planet-forming disk evolves from a cosmic-ray/X-ray-dominated regime to an ultraviolet-dominated chemical equilibrium. This, in turn, will bring about a temporal transition in the chemical reservoir from which planets will accrete. This photochemical dominated gas phase chemistry develops as dust evolves via growth, settling and drift, and the small grain population is depleted from the disk atmosphere. A higher gas-to-dust mass ratio allows for deeper penetration of ultraviolet photons is coupled with a carbon-rich gas (C/O > 1) to form carbon-bearing radicals and ions. This further results in gas phase formation of organic molecules, which then would be accreted by any actively forming planets present in the evolved disk.

SIPGI is a spectroscopic pipeline for the data reduction of optical/near-infrared data acquired by slit-based spectrographs. SIPGI is a complete spectroscopic data reduction environment retaining the high level of flexibility and accuracy typical of the standard "by-hand" reduction methods but with a significantly higher level of efficiency. This is obtained exploiting three main concepts: 1) a built-in data organiser to classify the data, together with a graphical interface; 2) the instrument model (analytic description of the main calibration relations); 3) the design and flexibility of the reduction recipes: the number of tasks required to perform a complete reduction is minimised, preserving the possibility to verify the accuracy of the main stages of data-reduction process. The current version of SIPGI manages data from the MODS and LUCI spectrographs mounted at the Large Binocular Telescope (LBT) with the idea to extend SIPGI to support other through-slit spectrographs.

Upcoming cosmological weak lensing surveys are expected to constrain cosmological parameters with unprecedented precision. In preparation for these surveys, large simulations with realistic galaxy populations are required to test and validate analysis pipelines. However, these simulations are computationally very costly -- and at the volumes and resolutions demanded by upcoming cosmological surveys, they are computationally infeasible. Here, we propose a Deep Generative Modeling approach to address the specific problem of emulating realistic 3D galaxy orientations in synthetic catalogs. For this purpose, we develop a novel Score-Based Diffusion Model specifically for the SO(3) manifold. The model accurately learns and reproduces correlated orientations of galaxies and dark matter halos that are statistically consistent with those of a reference high-resolution hydrodynamical simulation.

These lecture notes are based on the course "Gravitational waves from the early universe" given at the 27th W.E. Heraeus "Saalburg" Summer School 2021 by Valerie Domcke. Ongoing and future collaborations will probe different frequency ranges of the gravitational wave spectrum, allowing for probing different stages of the early universe and Beyond Standard Model physics. Due to the very high energies involved, accelerators cannot probe them. Therefore, current knowledge about new physics is limited and relies on bounds from CMB observations and theoretical assumptions about these energy scales. While some models are in tension with CMB data, others are unconstrained in shorter wavelength scales. Nonetheless, each one of these models has a gravitational wave density spectrum that can be compared to data. These lecture notes review the formalism of gravitational waves in General Relativity and introduce stochastic gravitational waves, primordial sources, and detection efforts.

The future astronomical imaging surveys are set to provide precise constraints on cosmological parameters, such as dark energy. However, production of synthetic data for these surveys, to test and validate analysis methods, suffers from a very high computational cost. In particular, generating mock galaxy catalogs at sufficiently large volume and high resolution will soon become computationally unreachable. In this paper, we address this problem with a Deep Generative Model to create robust mock galaxy catalogs that may be used to test and develop the analysis pipelines of future weak lensing surveys. We build our model on a custom built Graph Convolutional Networks, by placing each galaxy on a graph node and then connecting the graphs within each gravitationally bound system. We train our model on a cosmological simulation with realistic galaxy populations to capture the 2D and 3D orientations of galaxies. The samples from the model exhibit comparable statistical properties to those in the simulations. To the best of our knowledge, this is the first instance of a generative model on graphs in an astrophysical/cosmological context.

This study is based on observations of MWC560 during the last two observational seasons (2020/2021 and 2021/2022). Other than looking for flickering we were interested in following the variability of brightness in the same period. Looking for similarities in the spectra with other types of stars is also of great interest to us because it could help clarify the stellar configuration of such objects. Our observations during the last two observational seasons of MWC560 confirm the absence of flickering. From the similarities of the gathered spectra of XX Oph and MWC560 we assume that the components in XX Oph are a red giant and a white dwarf, which are also surrounded by a common shell.

KK242 is a LV dwarf of transition type residing in the void environment. Koda et al. present clear indications on its connection with Scd galaxy NGC6503. This implies the distance to KK242 of ~6.3 Mpc and its M_B = -10.5 mag. Its radial velocity, known from the Effelsberg radio telescope \HI\ observations, reveals, however, the difference with that of NGC6503, dV ~ 400 km/s. If real, this fact implies the substantial constraints on its origin. To clear-up the issue of KK242 radial velocity, we obtained with the SAO 6-m telescope spectra of its faint star-forming (SF) complex. H-alpha and H-beta emission is detected in two adjacent compact regions, the southern and northern, separated by ~2" (~60 pc). Their mean radial velocity is V_hel = -66 km/s, ~100 km/s lower than that of NGC6503. We use the HST Legacy Archive images and photometry of individual stars from the Extragalactic Distance Database, available for KK242, to identify in the SF complex the exciting hot stars, the probable BHeB and RHeB stars and a supernova remnant. We address, based on the possible range of its gas metallicity, the probable evolutionary paths of KK242. Using package Cloudy and parameters of the exciting B0V stars, we conclude that the observed flux ratio of [Sii] doublet to H-alpha is consistent with the value of 12+log(O/H) ~7.35+/-0.18 dex, expected for a stripped void dIrr galaxy.

Neutrino-cooled accretion flow around a spinning black hole, produced by a compact binary merger is a promising scenario for jet formation and launching magnetically-driven outflows. Based on GW170817 gravitational wave detection by LIGO and Virgo observatories followed by electromagnetic counterparts, this model can explain the central engine of the short duration gamma ray bursts (GRB) and kilonova radiations. Using the open-source GRMHD HARM-COOL code, we evolved several 2D magnetized accretion disk-black hole models with realistic equation of state in the fixed curved space-time background. We applied particle tracer technique to measure the properties of the outflows. The disk and black hole's initial parameters are chosen in a way to represent different possible post-merger scenarios of the merging compact objects. Our simulations show a strong correlation between black hole's spin and ejected mass. Generally, mergers producing massive disks and rapidly spinning black holes launch stronger outflows. We observed our models generate winds with moderate velocity ($v/c \sim 0.1-0.2$), and broad range of electron fraction. We use these results to estimate the luminosity and light curves of possible radioactively powered transients emitted by such systems. We found the luminosity peaks within the range of $10^{40}-10^{42}$ erg/s which agrees with previous studies for disk wind outflows.

Supersonic interacting flows occurring in phenomena such as protostellar jets give rise to strong shocks, and have been demonstrated in several laboratory experiments. To study such colliding flows, we use the AstroBEAR AMR code to conduct hydrodynamic simulations in three dimensions. We introduce variations in the flow parameters of density, velocity, and cross sectional radius of the colliding flows %radius in order to study the propagation and conical shape of the bow shock formed by collisions between two, not necessarily symmetric, hypersonic flows. We find that the motion of the interaction region is driven by imbalances in ram pressure between the two flows, while the conical structure of the bow shock is a result of shocked lateral outflows (SLOs) being deflected from the horizontal when the flows are of differing cross-section.

We present ALMA observations with a 800 au resolution and radiative-transfer modelling of the inner part ($r\approx6000$ au) of the ionized accretion flow around a compact star cluster in formation at the center of the luminous ultra-compact (UC) HII region G10.6-0.4. We modeled the flow with an ionized Keplerian disk with and without radial motions in its outer part, or with an external Ulrich envelope. The MCMC fits to the data give total stellar masses $M_\star$ from 120 to $200~M_\odot$, with much smaller ionized-gas masses $M_\mathrm{ion-gas} = 0.2$ to $0.25~M_\odot$. The stellar mass is distributed within the gravitational radius $R_g\approx 1000$ to 1500 au, where the ionized gas is bound. The viewing inclination angle from the face-on orientation is $i = 49$ to $56~\deg$. Radial motions at radii $r > R_g$ converge to $v_{r,0} \approx 8.7$ km/s, or about the speed of sound of ionized gas, indicating that this gas is marginally unbound at most. From additional constraints on the ionizing-photon rate and far-IR luminosity of the region, we conclude that the stellar cluster consists of a few massive stars with $M_\mathrm{star} = 32$ to $60~M_\odot$, or one star in this range of masses accompanied by a population of lower-mass stars. Any active accretion of ionized gas onto the massive (proto)stars is residual. The inferred cluster density is very large, comparable to that reported at similar scales in the Galactic Center. Stellar interactions are likely to occur within the next Myr.

Voids occupy about 3/4 of the volume of the Universe and contain about 15% of its mass. Due to various observational selection effects, these structure elements and galaxies populating voids, are highly under-explored. This especially relates to the lowest mass galaxies which comprise the main void population. Studying the nearby voids allows us to improve our understanding of the most elusive void objects. We present the brief overview of the current status and the prospects of the study of the nearest voids and their galaxies. First, we summarize the pioneer study of a hundred galaxies residing in the nearby Lynx-Cancer void which clearly evidence for the slower evolution of void galaxies and finds also the unusual very metal-poor and gas-rich dwarfs. Then we describe the recently defined sample of the nearby voids within the sphere with R = 25 Mpc and a sample of 1350 galaxies residing in these voids (~20% of all galaxies within this volume). We discuss the current results obtained for several directions of the study of this sample. They include: the search for Very Young Galaxies, the study of HI properties, the clustering of void galaxies and its relation to the void substructures, and the unbiased study of 260 void galaxies within the Local Volume (R < 11 Mpc). Altogether, this opens a perspective way to address the suggested peculiarities of the void galaxy formation and evolution. Finally, we briefly overview the expected advancements in the void galaxy studies related to the upcoming new facilities.

The multi-messenger anomalies, including spectral hardening or excess for nuclei, leptons, ratios of $\bar p/p$ and B/C, and anisotropic reversal, were observed in past years. AMS-02 experiment also revealed different spectral break for positron and electron at 284 GeV and beyond TeV respectively. It is natural to ask whether all those anomalies originate from one unified physical scenario. In this work, the spatially-dependent propagation (SDP) with a nearby SNR source is adopted to reproduce above mentioned anomalies. There possibly exists dense molecular cloud(DMC) around SNRs and the secondary particles can be produced by pp-collision or fragmentation between the accelerated primary cosmic rays and DMC. As a result, the spectral hardening for primary, secondary particles and ratios of $B/C$ and $\bar p/p$ can be well reproduced. Due to the energy loss at source age of 330 kyrs, the characteristic spectral break-off for primary electron is at about 1 TeV hinted from the measurements. The secondary positron and electron from charged pion take up $5\%$ energy from their mother particles, so the positron spectrum has a cut-off at $\sim$250 GeV. Therefore, the different spectral break for positron and electron together with other anomalies can be fulfilled in this unified physical scenario. More interesting is that we also obtain the featured structures as spectral break-off at 5 TV for secondary particles of Li, Be, B, which can be served to verify our model. We hope that those tagged structures can be observed by the new generation of space-borne experiment HERD in future.

We present a sample-variance-limited measurement of the temperature power spectrum ($TT$) of the cosmic microwave background (CMB) using observations of a $\sim\! 1500 \,\mathrm{deg}^2$ field made by SPT-3G in 2018. We report multifrequency power spectrum measurements at 95, 150, and 220GHz covering the angular multipole range $750 \leq \ell < 3000$. We combine this $TT$ measurement with the published polarization power spectrum measurements from the 2018 observing season and update their associated covariance matrix to complete the SPT-3G 2018 $TT/TE/EE$ data set. This is the first analysis to present cosmological constraints from SPT $TT$, $TE$, and $EE$ power spectrum measurements jointly. We blind the cosmological results and subject the data set to a series of consistency tests at the power spectrum and parameter level. We find excellent agreement between frequencies and spectrum types and our results are robust to the modeling of astrophysical foregrounds. We report results for $\Lambda$CDM and a series of extensions, drawing on the following parameters: the amplitude of the gravitational lensing effect on primary power spectra $A_\mathrm{L}$, the effective number of neutrino species $N_{\mathrm{eff}}$, the primordial helium abundance $Y_{\mathrm{P}}$, and the baryon clumping factor due to primordial magnetic fields $b$. We find that the SPT-3G 2018 $T/TE/EE$ data are well fit by $\Lambda$CDM with a probability-to-exceed of $15\%$. For $\Lambda$CDM, we constrain the expansion rate today to $H_0 = 68.3 \pm 1.5\,\mathrm{km\,s^{-1}\,Mpc^{-1}}$ and the combined structure growth parameter to $S_8 = 0.797 \pm 0.042$. The SPT-based results are effectively independent of Planck, and the cosmological parameter constraints from either data set are within $<1\,\sigma$ of each other. (abridged)

Sophisticated spectral energy distribution (SED) models describe dust attenuation and emission using geometry parameters. This treatment is natural since dust effects are driven by the underlying star-dust geometry in galaxies. An example is the Starduster SED model, which divides a galaxy into a stellar disk, a stellar bulge, and a dust disk. This work utilises the Starduster SED model to study the efficacy of inferring geometry parameters using spatially integrated SED fitting. Our method fits the SED model to mock photometry produced by combining a semi-analytic model with the same SED model. Our fitting results imply that the disk radius can be constrained, while the inclination angle, dust disk to stellar disk radius ratio, bulge radius and intrinsic bulge to total luminosity ratio are unconstrained, even though 21 filters from UV to FIR are used. We also study the impact of S/N, finding that the increase of S/N (up to 80) brings limited improvements to the results. We provide a detailed discussion to explain these findings, and point out the implications for models with more general geometry.

It has been suggested for a long time that dark matter would form a density spike around a black hole. However, no promising evidence has been observed so far to verify this theoretical suggestion. Here, we report the existence of a dark matter density spike around each of the two nearby stellar-mass black holes (A0620-00 and XTE J1118+480). The dynamical friction between dark matter and the companion stars can satisfactorily explain the abnormally fast orbital decays in the two binaries. The calculated spike index for A0620-00 and XTE J1118+480 are $\gamma=1.71^{+0.01}_{-0.02}$ and $\gamma=1.85^{+0.04}_{-0.04}$ respectively, which are close to the lower regime predicted by the stellar heating model. It may provide a possible indirect evidence for the existence of dark matter density spikes around stellar-mass black holes. We anticipate that analyzing observational data of nearby black hole X-ray binaries would be a new way to reveal the nature of dark matter.

We combined mapping-mode mid-infrared Spitzer spectra with complementary infrared imaging to perform a spatially resolved study of polycyclic aromatic hydrocarbons (PAHs) emission from the central regions of 66 nearby galaxies, roughly evenly divided into star-forming systems and low-luminosity active galactic nuclei (AGNs). In conjunction with similar measurements available for quasars, we aim to understand the physical properties of PAHs across a broad range of black hole accretion power, with the goal of identifying observational diagnostics that can be used to probe the effect of AGNs on the host galaxy. Whereas the PAH emission correlates tightly with far-ultraviolet luminosity in star-forming regions, the spatially resolved regions of AGNs tend to be PAH-deficient. Moreover, AGN regions exhibit on average smaller PAH 6.2 {\mu}m/7.7 {\mu}m and larger PAH 11.3 {\mu}m/7.7 {\mu}m band ratios. Although the current data are highly restrictive, they suggest that these anomalous PAH band ratios cannot be explained by the effects of the AGN radiation field alone. Instead, they hint that small grains may be destroyed by the combined effects of radiative processes and shocks, which are plausibly linked to jets and outflows preferentially associated with highly sub-Eddington, radiatively inefficient AGNs. While quasars also present a PAH deficit and unusual PAH band ratios, their characteristics differ in detail compared to those observed in more weakly accreting AGNs, a possible indicator of fundamental differences in their modes of energy feedback.

Emission from polycyclic aromatic hydrocarbons (PAHs), a commonly used indicator of star formation activity in galaxies, also has the potential to serve as an effective empirical tracer of molecular gas. We use a sample of 19 nearby galaxies with spatially resolved mid-infrared Spitzer spectroscopy, multi-wavelength optical and mid-infrared imaging, and millimeter interferometric CO(1-0) maps to investigate the feasibility of using PAH emission as an empirical proxy to estimate molecular gas mass. PAH emission correlates strongly with CO emission on sub-kpc scales over the diverse environments probed by our sample of star-forming galaxies and low-luminosity active galactic nuclei. The tight observed correlation, likely a consequence of photoelectronic heating of the diffuse interstellar gas by the PAHs, permits us to derive an empirical calibration to estimate molecular gas mass from the luminosity of PAH emission that has a total scatter of only {\sim} 0.2 - 0.25 dex. Mid-infrared bands sensitive to PAH emission (e.g., the Spitzer/IRAC4 and WISE/W3 filters) can also be used as a highly effective substitute for this purpose.

Emission from polycyclic aromatic hydrocarbons (PAHs) is a promising tool for estimating star formation rate (SFR) in galaxies, but the origin of its sources of excitation, which include not only young but possibly also old stars, remains uncertain. We analyze Spitzer mid-infrared mapping-mode spectroscopic observations of the nuclear and extra-nuclear regions of 33 nearby galaxies to study the contribution of evolved stars to PAH emission. In combination with photometric measurements derived from ultraviolet, H{\alpha}, and infrared images, the spatially resolved spectral decomposition enables us to characterize the PAH emission, SFR, and stellar mass of the sample galaxies on sub-kpc scales. We demonstrate that the traditional empirical correlation between PAH luminosity and SFR has a secondary dependence on specific SFR, or, equivalently, stellar mass. Ultraviolet-faint regions with lower specific SFRs and hence greater fraction of evolved stars emit stronger PAH emission at fixed SFR than ultraviolet-bright regions. We reformulate the PAH-based SFR estimator by explicitly introducing stellar mass as a second parameter to account for the contribution of evolved stars to PAH excitation. The influence of evolved stars can explain the sub-linear correlation between PAH emission and SFR, and it can partly account for the PAH deficit in dwarf galaxies and low-metallicity environments.

The origin of the inner Galactic emission, measured by COMPTEL with a flux of $\sim 10^{-2} ~{\rm MeV~ cm}^{-2}~ {\rm s}^{-1}~ {\rm sr}^{-1}$ in the 1-30 MeV energy range from the inner Galactic region, has remained unsettled since its discovery. In this paper, we elaborate on a model of individual MeV gamma-ray sources unresolved by COMPTEL. This is conducted for sources crossmatched between the Swift-BAT and Fermi-LAT catalogs by interpolating the energy spectra in the hard X-ray and GeV gamma-ray ranges, as well as unmatched sources between the two catalogs. We find that the source contribution to the COMPTEL emission would be at least ~20%. Combined with the Galactic diffuse emission, which is not well constrained, the COMPTEL emission can be roughly reproduced in some cases.

It was well known that most of gamma-ray bursts (GRBs) are dominated by positive spectral lags, while a small fraction of GRBs show negative lags. However, Wei et al. firstly identified a well-defined transition from positive lags to negative lags in GRB 160625B, and then got robust limits on possible violation of Lorentz Invariance (LIV) based on the observation. Recently, such a transition has been found in three different emission episodes in \thisgrb by Gunapati et al., which provides us a great opportunity to investigate whether the transition results from LIV-induced observed spectral lags. Our analysis shows that the LIV model can not be compatible with the current observations, whereas, only the spectral evolution induced spectral lags could responsible for the transition. So, spectral evolution can also explain the positive to negative lag in GRB 190530A.

Gas accretion of embedded stellar-mass black holes\,(sBHs) or stars in the accretion disk of active galactic nuclei\,(AGNs) will modify the mass distribution of these sBHs and stars, which will also affect the migration of the sBHs/stars. \textbf{With the introduction of the mass accretion effect, we simulate the evolution of the sBH/star distribution function in a consistent way by extending the Fokker-Planck equation of sBH/star distributions to the mass-varying scenario, and explore the mass distribution of sBHs in the nuclear region of the galaxy centre.} We find that the sBHs can grow up to several tens solar mass and form heavier sBH binaries, which will be helpful for us to understand the black-hole mass distribution as observed by the current and future ground-based gravitational wave detectors\,(e.g., LIGO/VIRGO, ET and Cosmic Explorer). We further estimate the event rate of extreme mass-ratio inspirals\,(EMRI) for sBH surrounding the massive black hole and calculate the stochastic gravitational wave\,(GW) background of the EMRIs. We find that the background can be detected in future space-borne GW detectors after considering the sBHs embedded in the AGN disk, while the mass accretion has a slight effect on the GW background.

Gamma-Ray Bursts (GRBs) can be used as a tool to probe the universe at high redshift. In this regard, the Amati relation which correlates the isotropic equivalent radiant energy $(E_{iso})$ and the spectral peak energy in the GRB rest frame $(E_p)$ allows us to use GRBs as distance indicators. However, the circularity issue that arises due to the lack of GRBs at low redshift has motivated several authors to come up with model-independent approaches to investigate this relation. For this same purpose, we use Hubble parameter measurements obtained from the differential age of the galaxies to circumvent the circularity problem. In this work, we apply a non-parametric approach namely Gaussian Process on the observational Hubble data (without assuming any cosmological model or parameters) to determine the luminosity distances needed to calculate $E_{iso}$. We find that the best fit values of the Amati relation parameters are in concordance with the earlier works.

This paper presents a detection of significant velocity offset between emission and absorption lines for a dual core system in SDSS~J155708.82+273518.74 (= SDSS~J1557). The photometric image of SDSS~J1557 exhibits clear two cores with a projected separation of $\sim$2.2 arcseconds (4.9 kpc) determined by GALFIT. Based on the applications of the commonly accepted pPXF code with 636 theoretical SSP templates, the host galaxy contribution can be well determined. Then, the emission line features of SDSS~J1557 can be well measured after subtraction of host starlight. It is found that the velocity offset of emission lines with respect to absorption lines reaches $458 \pm 13$ km/s. According to the Baldwin-Phillips-Terlevich (BPT) diagram, SDSS J1557 is a composite galaxy. In addition, SDSS J1557 can well fit the $M_{\rm BH}-\sigma_{\ast}$ relation of bulges and the galaxy merger would not change this relation. Two reasonable models (say, AGN-driven outflow vs. dual core system) have been discussed to explain this velocity offset. The model of AGN-driven outflow fails to interpret the systematic redshift of emission lines and similar velocity offsets for various emission lines in SDSS~J1557. A significant velocity offset between emission and absorption lines might be an effective indicator of dual core system.

The study of Active Galactic Nuclei (AGNs) is vital in order to understand their respective nuclei bounded activity primarily triggered by the accretion disk and the associated corona. The hard X-ray emission characterised by the spectral index $\Gamma$ emitted by these sources in 2-10 keV band ionizes the nearby physical features like disk, BLR, NLR and provides a radiative configuration. However, based on previous studies, there is degeneracy in the evolution of $\Gamma$ across the redshift, wherein few studies display a systematic trend in the evolution and others rule out the systematic variation. In the present work, we study the evolution of $\Gamma$ across the redshift for a very large sample and perform further studies on the variation of $\Gamma$ with optical parameters using SDSS data viz. H$_{\beta}$ \& Mg II and their luminosities. Based on analysis we find that there is no change in $\Gamma$ across the redshift and it does not show any correlation with the optical parameters. This strongly suggests that the $\Gamma$ is influenced by the soft excess and/or reflection component and is affecting the disk and BLR region at the same strength across the redshift.

On-going or soon to come cosmological large-scale structure surveys such as DESI, SPHEREx, Euclid, or the High-Latitude Spectroscopic Survey of the Nancy Grace Roman Space Telescope promise unprecedented measurement of the clustering of galaxies on large scales. When quantified with the Cartesian Fourier basis, the measurement of these large scales requires the introduction of so-called wide-angle corrections. By contrast, the measurement of the power spectrum in a spherical Fourier Bessel (SFB) basis does not require such corrections and naturally accounts for the spherical survey geometries. Here, we develop and implement a fast code to construct the SFB power spectrum and investigate how line of sight effects, physics such as non-Gaussianity, and differing survey geometries affect SFB power spectrum estimates. We then leverage our program to predict the tightness of cosmic growth constraints from realistic survey specifications using a Fisher matrix formalism.

A thermal relic WIMP remains a prime candidate for the nature of Dark Matter, particularly for the more poorly constrained case of a heavy ($\gtrsim$ 1 TeV) WIMP. The highest fluxes from WIMP annihilations are expected in the region of the Galactic Centre (GC) where current and near future gamma-ray observatories can be exploited to place tight limits on the WIMP paradigm. It is regularly noted that the annihilation flux of gammas will be accompanied by charged secondary particles which can produce 'delayed' inverse Compton (IC) gamma-ray emission, but this component is often neglected in indirect Dark Matter searches. In this work the inverse Compton emission is studied for the specific conditions of heavy WIMP annihilation in the GC. Using models for the magnetic and radiation fields of the region, and taking into consideration the transport of secondary particles, we find that for TeV WIMPs the IC component cannot be neglected in the GC, with the particles produced cooling within the region rather than propagating out in to the Galaxy. This effect changes the predicted spectral shape substantially and thus boosts the detection prospects for heavy WIMPs.

Evidence from the analysis of eclipsing binary systems revealed that late-type stars are larger and cooler than predicted by models, and that this is probably caused by stellar magnetic activity. In this work, we revisit this problem taking into account the advancements in the last decade. We provide and updated a list of 32 eclipsing binary or multiple systems, including at least one star with a mass $\lesssim 0.7$ M$_{\odot}$ and with mass and radius measured to an accuracy better than 3%. The~comparison with stellar structure and evolution theoretical models reveals an overall discrepancy of about 7% and -4% for the radius and effective temperature, respectively, and that it may be larger than previously found below the full convection boundary. Furthermore, the hypothesis of stellar activity is reinforced by the comparison of different systems with similar components. Further eclipsing binaries with accurately determined masses and radii, and with estimated activity levels, as well as the implementation of magnetic activity in theoretical models will help to improve our knowledge of low-mass stars, which are prime targets for exoplanet surveys.

Extreme ultraviolet imaging spectroscopic observations often show an increase in line width around the loop-top or above-loop-top (ALT) region of solar flares, suggestive of turbulence. In addition, recent spectroscopic observations found the oscillation in the Doppler velocity around the ALT region. We performed three-dimensional magnetohydrodynamic (MHD) simulations to investigate the dynamics in the ALT region, with a particular focus on the generation of turbulence and the excitation of the oscillatory motion. We found a rapid growth of MHD instabilities around the upper parts of the ALT region (arms of the magnetic tuning fork). The instabilities grow more rapidly than the magnetic Rayleigh-Taylor-type instabilities at the density interface beneath the reconnecting current sheet. Eventually, the ALT region is filled with turbulent flows. The arms of the magnetic tuning fork have bad-curvature and transonic flows. Therefore, we consider that the rapidly growing instabilities are combinations of pressure-driven and centrifugally driven Rayleigh-Taylor-type instabilities. Despite the presence of turbulent flows, the ALT region shows a coherent oscillation driven by the backflow of the reconnection jet. We examine the numerical results by re-analyzing the solar flare presented in Reeves et al. 2020. We find that the highest non-thermal velocity is always at the uppermost visible edge of the ALT region, where oscillations are present. This result is consistent with our models. We also argue that the turbulent magnetic field has a significant impact on the confinement of non-thermal electrons in the ALT region.

Bose-Einstein-condensed dark matter (BEC-DM), also called scalar field dark matter (SFDM), has become a popular alternative to the standard, collisionless cold dark matter (CDM) model, due to its long-held potential to resolve the small-scale crisis of CDM. Halos made of BEC-DM have been modelled using the Gross-Pitaevskii (GP) equation coupled to the Poisson equation; the so-called GPP equations of motion. These equations are based on fundamental microphysical conditions that need to be fulfilled in order for the equations to be valid in the first place, related to the diluteness of the DM gas and the nature of the scattering model. We use these conditions in order to derive the implications for the BEC-DM particle parameters, notably the self-interaction coupling strength $g$. We compare the bounds with the constraint that results from the assumption of virial equilibrium of the central cores of halos, deriving a relationship that connects the BEC-DM parameters $g$ and particle mass $m$. We find that the GPP conditions are greatly fulfilled, for BEC-DM particle masses of interest, if such models also obey the virial condition that turns out to be the strongest constraint. We also derive the implications for the elastic scattering cross section (per particle mass) in BEC-DM halos, based on the scattering model of GPP, and find a huge range of possible values, depending on the strength of self-interaction. We put our results into context to recent literature which predicts sub-kpc core size.

We study possible ways gravitational waves (GW) get sourced in a theory with minimal left-right symmetry breaking. First order phase transitions generically lead to gravitational waves sourced by bubble collisions, while second order phase transitions (SOPT) do not. Interesting variants on the standard classification of phase transitions occur due to the breaking of discrete parity combined with the limitation of light cone in the early Universe. If local effective potential signals SOPT or a cross over, breaking of discrete parity in conjunction with finiteness of the causal horizon leads to a \textsl{causal horizon limited} second order phase transition, which results in domain walls separating left-like and right-like domains. On the other hand for the case of usual first order phase transition (FOPT), we get the usual signal from spontaneously created bubbles, but also, as we argue from a lingering late time domain wall structure separating the two types of vacua. Thus both FOPT and putative SOPT give rise to distinct features in the spectrum of GW. The signatures are testable via experiments such as IPTA, and DECIGO and LISA. Finally we point out that a version of the left-right symmetric model which separates parity breaking from gauge symmetry breaking is also subject to domain wall formation and amenable to GW observations.

We present a study of optically selected dual AGN with projected separations of 3--97~kpc. Using multi-wavelength (MWL) information (optical, X-rays, mid-IR), we characterized the intrinsic nuclear properties of this sample and compared them with those of isolated systems. Among the 124 X-ray detected AGN candidates, 52 appear in pairs and 72 as single X-ray sources. Through MWL analysis, we confirmed the presence of the AGN in a fraction >80\% of the detected targets in pairs (42 over 52). X-ray spectral analysis confirms the trend of increasing AGN luminosity with decreasing separation, suggesting that mergers may have contributed in triggering more luminous AGN. Through X/mid-IR ratio $vs$ X-ray colors, we estimated a fraction of Compton-thin AGN (with 10$^{22}$ cm$^{-2}$ $<$ N$_{\rm H} <$10$^{24}$ cm$^{-2}$) of about 80\%, while about 16\% are Compton thick (CT, with N$_{\rm H}>$10$^{24}$ cm$^{-2}$) sources. These fractions of obscured sources are larger than those found in samples of isolated AGN, confirming that pairs of AGN show higher obscuration. This trend is further confirmed by comparing the de-reddened [O\ III] emission with the observed X-ray luminosity. However, the derived fraction of Compton-thick sources in this sample at early stage of merging is lower than reported for late-merging dual-AGN samples. Comparing N$_{\rm H}$ from X-rays with that derived from E(B-V) from Narrow Line Regions, we find that the absorbing material is likely associated with the torus or the Broad Line Regions. We also explored the X-ray detection efficiency of dual-AGN candidates, finding that, when observed properly (at on-axis positions and with long exposures), X-ray data represent a powerful way to confirm and investigate dual-AGN systems.

Photospheric C, N, and O abundances of 118 solar-analog stars were determined by applying the synthetic-fitting analysis to their spectra in the blue or near-UV region comprising lines of CH, NH, and OH molecules, with an aim of clarifying the behaviors of these abundances in comparison with [Fe/H]. It turned out that, in the range of -0.6<[Fe/H]<+0.3, [C/Fe] shows a marginally increasing tendency with decreasing [Fe/H] with a slight upturn around [Fe/H]~0, [N/Fe] tends to somewhat decrease towards lower [Fe/H], and [O/Fe] systematically increases (and thus [C/O] decreases) with a decrease in [Fe/H]. While these results are qualitatively consistent with previous determinations mostly based on atomic lines, the distribution centers of these [C/Fe], [N/Fe], and [O/Fe] at the near-solar metallicity are slightly negative by several hundredths dex, which is interpreted as due to unusual solar abundances possibly related to the planetary formation of our solar system. However, clear anomalies are not observed in the [C,N,O/Fe] ratios of planet-host stars. Three out of four very Be-deficient stars were found to show anomalous [C/Fe] or [N/Fe] which may be due to mass transfer from the evolved companion, though its relation to Be depletion mechanism is still unclear.

The main goal of the paper is to explore why observations of many astrospheres (or circumstellar bubbles) show quite stable and smooth structures of astropauses - the tangential discontinuities separating the stellar and interstellar winds, - while both theory and numerical simulations suggest that tangential discontinuities are unstable due to well known Kelvin-Helmholtz (K-H) instability. It was recognized before that magnetic fields may stabilize the astropauses. In this paper, we explore another mechanism to reduce the K-H instability of the astropauses. This mechanism is a periodic change of the stellar wind dynamic pressure. Fluctuations of the stellar wind parameters are quite expected. For example, the Sun has an 11-year cycle of global activity although there are also shorter periods of the solar wind fluctuations. We performed the parametric numerical study and demonstrate that the development of the K-H instability depends on the dimensionless parameter $\chi$ which is the ratio of the stellar wind terminal speed and interstellar flow speed. The larger the parameter $\chi$, the larger the fluctuations caused by the K-H instability. It has been shown that the K-H instability is convective which agrees with the previous linear analysis. The stabilization of the astropause by the periodic fluctuations in the stellar wind lead is demonstrated. It is shown that for the solar wind the most effective stabilization occurs when the period of stellar parameter change is about 1-4 years. For the eleven year solar cycle, the stabilization effect is weaker.

We present the 0.85-2.5-micron discovery spectrum and multi-epoch photometry of the new high-redshift quasar ULAS J081621.47+213442.6 obtained using the GNIRS spectrograph of the Gemini North and near-infrared wide-field camera of the 4-m UKIRT telescopes. The redshift of ULAS J081621.47+213442.6 measured from the MgII 2799 emission line is z=7.461. The absolute magnitude of the quasar is M1450=-25.33. The black hole mass estimated using the MgII 2799 line and Eddington accretion rate are ~5x10^8Msun and ~0.7. The spectrum of ULAS J081621.47+213442.6 exhibits strong NIII] 1750 emission line of a rest-frame equivalent width of ~12.5 A. The high abundance of nitrogen suggests that ULAS J081621.47+213442.6 may be at the peak of the nitrogen enrichment of the circumnuclear region by the asymptotic giant branch stars, which is expected ~0.25 Gyr after the bulk of star formation. The age of the starburst of ULAS J081621.47+213442.6 implied by the high nitrogen abundance, indicates that the active phase of the black hole growth of the quasar may have lasted only ~0.25 Gyr, favoring a massive initial black hole seed. We also observed the flux variations of the UV continuum of ULAS J081621.47+213442.6 caused by the variation in the line-of-sight absorbing column density on a rest-frame timescale of ~47 d. The estimated hydrogen column density of the gas cloud responsible for this variation is N~10^{23.5} cm^{-2}, consistent with the typical column density of mostly neutral, gravitationally bound clouds of the broad line region of quasars.

The Javalambre Photometric Local Universe Survey (J-PLUS) is a 12-band photometric survey using the 83-cm JAST telescope. Data Release 3 includes 47.4 million sources (29.8 million with $r \le 21$) on 3192 deg$^2$ (2881 deg$^2$ after masking). J-PLUS DR3 only provides star-galaxy classification so that quasars are not identified from the other sources. Given the size of the dataset, machine learning methods could provide a valid alternative classification and a solution to the classification of quasars. Our objective is to classify J-PLUS DR3 sources into galaxies, stars and quasars, outperforming the available classifiers in each class. We use an automated machine learning tool called {\tt TPOT} to find an optimized pipeline to perform the classification. The supervised machine learning algorithms are trained on the crossmatch with SDSS DR12, LAMOST DR7 and \textit{Gaia} DR3. We checked that the training set of about 570 thousand galaxies, one million stars and 220 thousand quasars is both representative and pure to a good degree. We considered 37 features: besides the twelve photometric bands with their errors, six colors, four morphological parameters, galactic extinction with its error and the PSF relative to the corresponding pointing. After exploring numerous pipeline possibilities through the TPOT genetic algorithm, we found that XGBoost provides the best performance: the AUC for galaxies, stars and quasars is above 0.99 and the average precision is above 0.99 for galaxies and stars and 0.94 for quasars. XGBoost outperforms the star-galaxy classifiers already provided in J-PLUS DR3 and also efficiently classifies quasars. We also found that photometry was very important in the classification of quasars, showing the relevance of narrow-band photometry.

Fast radio bursts (FRBs) are bright millisecond radio bursts at cosmological distances. Only a small fraction of FRBs apparently repeat. Polarization, a fundamental property of electromagnetic signals, often carries critical information about the radiation processes, the environment, and the intervening medium of FRBs. Here we report circular polarization detections of two active repeating FRBs, namely FRBs 20121102A and 20190520B, with the Five-hundred-meter Aperture Spherical radio Telescope. We detect circular polarization in both active repeating FRBs, which increases the number of repeating FRBs with circular polarization to three. In one of the bursts of FRB 20121102A, we detect 64% degree of circular polarization. The observed circular polarization is unlikely induced by multipath propagation. Our observations favor circular polarization induced by Faraday conversion or radiation mechanism intrinsic to the FRB source. The conditions to generate circular polarization have to be rare in either case.

We study the effect of baryons on the cosmic web -- halos, filaments, walls, and voids. To do so, we apply a modified version of NEXUS, a cosmic web morphological analysis algorithm, to the IllustrisTNG simulations. We find that halos lose more than $10\%$ of their mass due to baryons, mostly to filaments and a small portion to walls and voids. However, the mass transfer does not significantly shift the boundaries of structures, leaving the volume fractions of the cosmic structures largely unaffected. We quantify the effects of baryonic feedback on the power spectrum and the probability density function (PDF) of the density field for individual cosmic structures. For the power spectrum, most suppression due to feedback can be accounted for by including $M\ge10^{12}~M_\odot/h$ halos, without considering other cosmic structures. However, when examining the PDF of the density field, we find nearly $100\%$ suppression of the emptiest regions and $10\%$-level effects (boost or suppression) in the remaining regions of filaments, walls, and voids. Our results indicate the importance of modeling the effects of baryons in the whole cosmic web, not just halos, for cosmological analysis beyond two-point statistics or field-based inferences.

We present the results of a new reverberation mapping campaign for the broad-lined active galactic nucleus (AGN) in the edge-on spiral IC4329A. Monitoring of the optical continuum with $V-$band photometry and broad emission-line flux variability with moderate-resolution spectroscopy allowed emission-line light curves to be measured for H$\beta$, H$\gamma$, and HeII $\lambda 4686$. We find a time delay of $16.3^{+2.6}_{-2.3}$ days for H$\beta$, a similar time delay of $16.0^{+4.8}_{-2.6}$ days for H$\gamma$, and an unresolved time delay of $-0.6^{+3.9}_{-3.9}$ days for HeII. The time delay for H$\beta$ is consistent with the predicted value from the relationship between AGN luminosity and broad line region radius, after correction for the $\sim2.4$mag of intrinsic extinction at 5100A. Combining the measured time delay for H$\beta$ with the broad emission line width and an adopted value of $\langle f \rangle = 4.8$, we find a central supermassive black hole mass of $M_{\rm BH}=6.8^{+1.2}_{-1.1}\times10^7 M_{\rm \odot}$. Velocity-resolved time delays were measured across the broad H$\beta$ emission-line profile and may be consistent with an ''M''-like shape. Modeling of the full reverberation response of H$\beta$ was able to provide only modest constraints on some parameters, but does exhibit agreement with the black hole mass and average time delay. The models also suggest that the AGN structure is misaligned by a large amount from the edge-on galaxy disk. This is consistent with expectations from the unified model of AGNs, in which broad emission lines are expected to be visible only for AGNs that are viewed at relatively face-on inclinations.

We investigate deviations in the mean phase travel time of acoustic waves preceding the emergence of 46 large active regions observed by the Helioseismic and Magnetic Imager (HMI). In our investigation, we consider two different procedures for obtaining the mean phase travel time, by minimizing the difference between cross-correlations and a reference, as well as the Gabor wavelet fitting procedure. We cross-correlate the time series of mean phase travel time deviations with the surface magnetic field and determine the peak correlation time lag. We also compute the perturbation index--the area integrated mean phase travel time deviations exceeding quiet sun thresholds--and compare the time of peak perturbation index with the correlation time lag. We find that the lag times derived from the difference minimization procedure precede the flux emergence for 36 of the 46 active regions, and that this lag time has a noticeable correlation with the maximum flux rate. However, only 28 of the active regions have peak perturbation index times in the range of 24 to 48 hours prior to the flux emergence. Additionally, we examine the relationship between properties of the emerged active regions and the strength of helioseismic signals prior to their emergence.

The physics of baryons in halos, and their subsequent influence on the total matter phase space, has a rich phenomenology and must be well understood in order to pursue a vast set of questions in both cosmology and astrophysics. We use the CAMELS simulation suite to quantify the impact of four different galaxy formation parameters/processes (as well as two cosmological parameters) on the concentration-mass relation, $c_{\rm vir} - M_{\rm vir}$. We construct a simulation-informed nonlinear model for concentration as a function of halo mass, redshift, and 6 cosmological/astrophysical parameters. This is done for two galaxy formation models, IllustrisTNG and SIMBA, using 1000 simulations of each. We extract the imprints of galaxy formation across a wide range in mass $M_{\rm vir} \in [10^{11}, 10^{14.5}] M_{\rm \odot}/h$ and in redshift $z \in [0,6]$ finding many strong mass- and redshift-dependent features. Comparisons between the IllustrisTNG and SIMBA results show the astrophysical model choices cause significant differences in the mass and redshift dependence of these baryon imprints. Finally, we use existing observational measurements of $c_{\rm vir} - M_{\rm vir}$ to provide rough limits on the four astrophysical parameters. Our nonlinear model is made publicly available and can be used to include CAMELS-based baryon imprints in any halo model-based analysis.

Bouncing models of cosmology, as they arise e.g. in loop quantum cosmology, can generate close-to-scale-invariant fluctuation spectra as observed in the Cosmic Microwave Background (CMB). However, they are typically not Gaussian and also generate a bispectrum. It was proposed that these models can help to mitigate the large-scale anomalies of the CMB by considering large non-Gaussianities on very large scales, which decay exponentially on sub-horizon scales. It was therefore thought that this non-Gaussianity would not be visible in observations, which can only probe sub-horizon scales. In this letter we show that bouncing models with parameters such that they can mitigate the large-scale anomalies of the CMB are excluded by the Planck data with high significance of, depending on the specific model, $6.4$ or $14$ standard deviations.

The presence of neutral C60 fullerenes in circumstellar environments has been firmly established by astronomical observations as well as laboratory experiments and quantum-chemistry calculations. However, the large variations observed in the C60 17.4um/18.9um band ratios indicate that either additional emitters should contribute to the astronomical IR spectra or there exist unknown physical processes besides thermal and UV excitation. Fullerene-based molecules such as metallofullerenes and fullerene-adducts are natural candidate species as potential additional emitters, but no specific species has been identified to date. Here we report a model based on quantum-chemistry calculations and IR spectra simulation of neutral and charged endo(exo)hedral metallofullerenes, showing that they have a significant contribution to the four strongest IR bands commonly attributed to neutral C60. These simulations may explain the large range of 17.4um/18.9um band ratios observed in very different fullerene-rich circumstellar environments like those around planetary nebulae and chemically peculiar R Coronae Borealis stars. Our proposed model also reveals that the 17.4um/18.9um band ratio in the metallofullerenes simulated IR spectra mainly depends on the metal abundances, ionization level, and endo/exo concentration in the circumstellar envelopes. We conclude that metallofullerenes are potential emitters contributing to the observed IR spectra in fullerene-rich circumstellar envelopes. Our simulated IR spectra indicate also that the James Webb Space Telescope has the potential to confirm or refute the presence of metallofullerenes (or even other fullerene-based species) in circumstellar enviroments.

Magnetic fields permeate the entire Galaxy and are essential to, for example, the regulation of several stages of the star formation process and cosmic ray transportation. Unraveling its properties, such as intensity and topology, is an observational challenge that requires combining different and complementary techniques. The polarization of starlight due to the absorption by field-aligned non-spherical dust grains provides a unique source of information about the interstellar magnetic field in the optical band. This work introduces a first analysis of a new catalog of optical observations of linearly polarized starlight in the diffuse interstellar medium (ISM), the Interstellar Polarization Survey, General ISM (IPS-GI). We used data from the IPS-GI, focusing on 38 fields sampling lines of sight in the diffuse medium. The fields are about 0.3$^{\circ}$ by 0.3$^{\circ}$ in size and each of them contains $\sim1000$ stars on average. The IPS-GI catalog has polarimetric measurements of over 40000 stars, over 18000 of which have ${P}/\sigma_{P} > 5$. We added distances and other parameters from auxiliary catalogs to over 36000 of these stars. We analyzed parameter distributions and correlations between parameters of a high-quality subsample of 10516 stars (i.e. $\sim275$ stars per field). As expected, the degree of polarization tends to increase with the extinction, producing higher values of polarization at greater distances or at lower absolute Galactic latitudes. Furthermore, we find evidence for a large-scale ordered Galactic magnetic field.

We present images at 6 and 14 GHz of Source I in Orion-KL. At higher frequencies, from 43 to 340 GHz, images of this source are dominated by thermal emission from dust in a 100 AU diameter circumstellar disk, but at 6 and 14 GHz the emission is elongated along the minor axis of the disk, aligned with the SiO bipolar outflow from the central object. Gaussian fits to the 6, 14, 43, and 99 GHz images find a component along the disk minor axis whose flux and length vary with frequency consistent with free-free emission from an ionized outflow. The data favor a broad outflow from a disk wind, rather than a narrow ionized jet. Source I was undetected in higher resolution 5 GHz e-MERLIN observations obtained in 2021. The 5-6 GHz structure of SrcI may be resolved out by the high sidelobe structure of the e-MERLIN synthesized beam, or be time variable.

We have developed a pipeline called \mentari to generate the far-ultraviolet to far-infrared spectral energy distribution (SED) of galaxies from the \dustysage semi-analytic galaxy formation model (SAM). \dustysage incorporates dust-related processes directly on top of the basic ingredients of galaxy formation like gas infall, cooling, star formation, feedback, and mergers. We derive a physically motivated attenuation model from the computed dust properties in \dustysage, so each galaxy has a self-consistent set of attenuation parameters based on the complicated dust physics that occurred across the galaxy's assembly history. Then, we explore several dust emission templates to produce infrared spectra. Our results show that a physically-motivated attenuation model is better for obtaining a consistent multi-wavelength description of galaxy formation and evolution, compared to using a constant attenuation. We compare our predictions with a compilation of observations and find that the fiducial model is in reasonable agreement with: (i) the observed $z=0$ luminosity functions from the far-ultraviolet to far-infrared simultaneously, and hence (ii) the local cosmic SED in the same range, (iii) the rest-frame K-band luminosity function across $0 < z < 3$, and (iv) the rest-frame far-ultraviolet luminosity function across $0 < z < 1$. Our model underproduces the far-ultraviolet emission at $z=2$ and $z=3$, which can be improved by altering the AGN feedback and dust processes in \dustysage. However, this combination thus worses the agreement at $z=0$, which suggests that more detailed treatment of such processes is required.

We present the Galactic disk vertical velocity analysis using OB type stars (OB), Red Clump stars (RC), and Main-Sequence-Turn-Off stars (MSTO) with different average age populations crossed matched with LAMOST DR5 and Gaia DR3. We reveal the vertical velocities of the three populations varies clearly with the Galactocentric distance ($R$) and the younger stellar population has stronger increasing trend in general. The bending and breathing modes indicated by the vertical motions are dependent on the populations and they are varying with spatial locations. These vertical motions may be due to the Galactic warp, or minor mergers, or non-equilibrium of the disk. Assuming the warp is the dominant component, we find that the warp amplitude ($\gamma$, $Z_\omega$) for OB (younger population) is larger than that for RC (medium population) and the later one is also larger than that for MSTO (older population), which is in agreement with other independent analyses of stellar density distribution, and supports the warp is long-lived, non-steady structure and has time evolution. This conclusion is robust whether or not the line-of-nodes $\phi_w$ is fixed or as a free parameter (with $\phi_w$ is around 3$-$8.5$^{\circ}$ as best fit). Furthermore, we find that warp is lopsided with asymmetries along azimuthal angle ($\phi$).

Relativistic collisionless shocks are associated with efficient particle acceleration when propagating into weakly magnetized homogeneous media; as the magnetization increases, particle acceleration becomes suppressed. We demonstrate that this changes when the upstream carries kinetic-scale inhomogeneities, as is often the case in astrophysical environments. We use fully-kinetic simulations to study relativistic perpendicular shocks in magnetized pair plasmas interacting with upstream density perturbations. The upstream fluctuations are found to corrugate the shock front and generate large-scale turbulent shear motions in the downstream, which in turn are capable of accelerating particles. This can revive relativistic magnetized shocks as viable energization sites in astrophysical systems, such as jets and accretion disks. The generation of large-scale magnetic structures also has important implications for polarization signals from blazars.

Models invoking magnetic reconnection as the particle acceleration mechanism within relativistic jets often adopt a gradual energy dissipation profile within the jet. However, such a profile has yet to be reproduced in first-principles simulations. Here, we perform a suite of 3D general relativistic magnetohydrodynamic simulations of post-neutron star merger disks with an initially purely toroidal magnetic field. We explore the variations in both the microphysics (e.g., nuclear recombination, neutrino emission) and system parameters (e.g., disk mass). In all our simulations, we find the formation of magnetically striped jets. The stripes result from the reversals in the poloidal magnetic flux polarity generated in the accretion disk. The simulations display large variations in the distributions of stripe duration, $\tau$, and power, $\langle P_{\Phi} \rangle$. We find that more massive disks produce more powerful stripes, the most powerful of which reaches $\langle P_{\Phi} \rangle \sim 10^{49}$~erg~s$^{-1}$ at $\tau \sim 20$~ms. The power and variability that result from the magnetic reconnection of the stripes agree with those inferred in short duration gamma-ray bursts. We find that the dissipation profile of the cumulative energy is roughly a power-law in both radial distance, $z$, and $\tau$, with the slope in the range, $\sim 1.7-3$; more massive disks display larger slopes.

Understanding the statistical properties of synchrotron emission from our Galaxy is valuable from the perspective of observations targeting signals of cosmological origin, as well as for understanding physical processes in our Galaxy. In this work we extend the analysis of~\cite{Rahman:2021azv} to - (a) all-sky observed maps at different frequencies provided by WMAP, Planck and Stockert-Villa, (b) component separated synchrotron maps provided by WMAP, Planck and BeyondPlanck, and (c) component separated polarization maps provided by WMAP and Planck. Our main goals are to understand the variation of morphological properties with frequency, test the agreement between different component separation pipelines, and understand the nature of non-Gaussianity and statistical isotropy of synchrotron fluctuations towards smaller scales. The tools we use are Minkowski functionals and tensors. We find that the frequency maps exhibit large variation of the morphological features which persist after smoothing the maps at different scales. These differences are most likely due to the presence of different foreground components, besides synchrotron. For the component separated synchrotron temperature maps, we find that towards small scales BeyondPlanck and WMAP \texttt{MCMC-e} maps exhibit relatively good agreement with the Haslam map. We also find that all maps exhibit kurtosis-type non-Gaussianity in agreement with the Haslam map. However, the levels of deviations are much higher except for WMAP \texttt{MCMC-e} and BeyondPlanck which are comparable. We also find that the component separated maps tend towards isotropy on smaller scales, except for Planck and WMAP \texttt{MCMC-c} and \texttt{g} maps. The smaller scale non-Gaussian deviations of $E$ and $B$-mode polarization maps from WMAP and Planck are found to agree well with the kurtosis type non-Gaussianity.

Context. Cassiopeia A is one of the most extensively studied supernova remnants (SNRs) in our Galaxy. The analysis of its spectral features with the help of low frequency observations plays an important role for understanding the evolution of the radio source through the propagation of synchrotron emission to observers through the SNR environment and the interstellar medium. Aims. In this paper we present measurements of the integrated spectrum of Cas A to characterize the properties of free-free absorption towards this SNR. We also add new measurements to track its slowly evolving and decreasing integrated flux density. Methods. We use the Giant Ukrainian radio telescope (GURT) for measuring the continuum spectrum of Cassiopeia A within the frequency range of 16-72 MHz. The radio flux density of Cassiopeia A relative to the reference source of the radio galaxy Cygnus A has been measured on May-October, 2019 with two subarrays of the GURT, used as a two-element correlation interferometer. Results. We determine magnitudes of emission measure, electron temperature and an average number of charges of the ions for both internal and external absorbing ionized gas towards in Cassiopeia A. Generally, their values are close to the ones suggested by Arias et al. (2018), although for some there are slight differences. In the absence of clumping we find the unshocked ejecta of M = 2.61 solar mass at the electron density of 15.3 cm^-3 has a gas temperature of T=100 K. If the clumping factor is 0.67, then the unshocked ejecta of 0.96 solar mass the electron density of 18.7 cm^-3. Conclusions. The integrated flux density spectrum of Cassiopeia A obtained with the GURT interferometric observations is consistent with the theoretical model within measurement errors and also reasonably consistent with other recent results in the literature.

The sample of time-delay gravitational lenses appropriate for studying the geometry of the Universe continues to grow as dedicated campaigns, such as the Dark Energy Survey, the VST ATLAS survey, and the Large Synoptic Survey Telescope, complete their census of high-redshift sources. This catalog now includes hundreds of strong lensing systems, at least 31 of which have reasonably accurate time delay measurements. In this paper, we use them to compare the predictions of two competing Friedmann-Lemaitre-Robertson-Walker models: flat LCDM, characterized by two adjustable parameters (H_0 and Omega_m), and the R_h=ct universe (with H_0 as the single free variable). Over the past decade, the latter has accounted for the data better than the standard model, most recently the emergence of well-formed galaxies discovered by JWST at cosmic dawn. Here we show that the current sample of time-delay lenses favours R_h=ct with a likelihood of ~84% versus ~16% for the standard model. This level of accuracy will greatly improve as the ongoing surveys uncover many thousands additional lens systems over the next several years.

The symmetron field has an environment (density) dependent behavior which is a common feature of the models with the screening mechanism and results in a rich phenomenology. This model can produce domain walls between regions with different densities. We consider this aspect and study the physics of domain walls in between (underdensity) voids and (overdensity) halo structures. The (spherical) domain walls exert a repulsive force on a test mass outside of the wall while a test mass inside of the wall sees no force. This makes the structures outside the voids go further to a larger radius. Effectively, this means the voids are becoming larger in this scenario in comparison to the standard model of cosmology. Interestingly, this makes voids emptier which may shed light on Peebles' void phenomenon.

Motivated by the increasing number of gravitational wave detections of merging black holes (BHs) by LIGO-VIRGO-KAGRA, BH binary mergers in the discs of active galactic nuclei (AGN) is investigated as a possible merger channel. In this pathway, BHs in the large gas disc are expected to encounter one another, become mutually bound as a BH binary system through interaction with the gas in the disc and subsequently inspiral through gravitational torques induced by the local gas until merger. To more precisely determine the feasibility of this merger pathway, we present the first 3D global hydrodynamic simulations of both the formation and evolution of a stellar-mass BH binaries AGN discs with three different AGN disc masses and five different initial radial separations. These 15 different simulations show that binary capture can be successful in a range of local gas densities including cases well below that of a standard radiatively efficient alpha disc, and identify that the majority of these captured binaries are then subsequently hardened by the surrounding gas. Both prograde and retrograde binaries are formed in our simulations, which tend to be slowly hardened by the gas. The eccentricity evolution is found to depend strongly on the orbital rotation of the binary where prograde binaries are governed by gravitational torques form their circumbinary mini-disc, with eccentricities being damped over time, while for retrograde binaries the eccentricities are excited to >0.9 by accretion torques. In two cases, retrograde binaries ultimately undergo a close periapsis passage which results in a merger via gravitational waves after only a few thousand binary orbits. Thus, the merger timescale for retrograde binaries can be far shorter than the lifetime of an AGN disc. These simulations support the case that the AGN merger pathway could be highly efficient for merging BHs.

We analyze the effects of modified gravity on specific heats of electrons and ions, Debye temperature, crystallization process, and cooling mechanism in white dwarfs. We derive the Lane-Emden-Chandrasekhar equation and relate it to the cooling process equations for Palatini $f(R)$ gravity. Moreover, for the first time in the literature, we show that the gravity model plays a crucial role not only in the mass and size of the white dwarf, but also affects their internal properties. We further demonstrate that modified gravity can decrease the cooling age significantly.

Vector Dark Matter (VDM) that couples to lepton flavor ($L_e$, $L_{\mu}$, $L_{\tau}$) acts similarly to a chemical potential for the neutrino flavor eigenstates and modifies neutrino oscillations. VDM imparts unique signatures such as time and directional dependence with longer baselines giving better sensitivity. We use the non-observation of such a signal at Super-Kamiokande to rule out the existence of VDM in a region of parameter space several orders of magnitude beyond other constraints and show the projected reach of future experiments such as DUNE.

The dark matter relic density may be governed by the presence of new mediators that connect the dark matter field with the Standard Model particles. When the dark matter particle mass is larger than the mediator's, the pair production of mediators is kinematically open. This setup is known in the literature as secluded dark matter. Motivated by the appearance of secluded dark matter in several model building endeavours, we investigate the sensitivity of TeV gamma-ray instruments in the Southern Hemisphere namely, H.E.S.S., CTA, and SWGO to secluded dark matter annihilating in the Galactic Halo. We exploit the complementarity aspects of these detectors to find restrictive bounds on the annihilation cross-section for different annihilation channels. In particular, for a dark matter particle mass of $2$~TeV, H.E.S.S. is able to constraint $\langle \sigma v \rangle \geq 4 \times 10^{-26}\,\, {\rm cm}^3\, {\rm s}^{-1}$ at 95\% confidence level for the $4q$ and $4\tau$ channel, while CTA will be sensitive to $\langle \sigma v \rangle \geq 7 \times 10^{-27}\,\, {\rm cm}^3\, {\rm s}^{-1}$ and SWGO $\langle \sigma v \rangle \geq 6 \times 10^{-27}\,\, {\rm cm}^3\, {\rm s}^{-1}$ for the $4\tau$ channel, both well below the thermal relic cross-section. In fact, the combination of CTA and SWGO will be able to probe cross-sections below the thermal relic value for dark matter particles in the whole mass range between 100 GeV and 100 TeV in the $4q$ and $4\tau$ channels, and between 100 GeV and $\sim$40 TeV in the $4b$ channel.

We have explored a completely new and alternative way to restrict the parameter space of Horndeski theory of gravity. Using its Newtonian limit, it is possible to test the theory at a regime where, given its complexity and the small magnitude of the expected effects, it is poorly probed. At Newtonian level, it gives rise to a generalized Yukawa-like Newtonian potential which we have tested using S2 star orbit data. Our model adds five parameters to the General Relativity model, and the analysis constrains two of them with unprecedented precision to these energy scales, while only gives an exclusion region for the remaining parameters. We have shown the potential of weak-field tests to constrain Horndeski gravity opening, as a matter of fact, a new avenue that deserves to be further, and deeply, explored near in the future.

Over the last two decades, the possibility of using RPCs in outdoors systems has increased considerably. Our group has participated in this effort having installed several systems and continues to work on their optimization, while simultaneously studying and developing new approaches that can to use of RPCs in outdoor applications. In particular, some detectors were deployed in the field at the Pierre Auger Observatory in 2019 remained inactive, awaiting the commissioning of support systems. During the pandemic the detectors were left without gas flow for more than two years, but were recently reactivated with no major problems. The LouMu project combines particle physics and geophysics in order to map large geologic structures, using Muon Tomography. The development of the RPC system used and the data from the last two years will be presented. Finally, recent advances in a large area (1 m2) double gap-sealed RPC will be presented.

LISA is a space-based mHz gravitational-wave observatory, with a planned launch in 2034. It is expected to be the first detector of its kind, and will present unique challenges in instrumentation and data analysis. An accurate pre-flight simulation of LISA data is a vital part of the development of both the instrument and the analysis methods. The simulation must include a detailed model of the full measurement and analysis chain, capturing the main features that affect the instrument performance and processing algorithms. Here, we propose a new model that includes, for the first time, proper relativistic treatment of reference frames with realistic orbits; a model for onboard clocks and clock synchronization measurements; proper modeling of total laser frequencies, including laser locking, frequency planning and Doppler shifts; better treatment of onboard processing and updated noise models. We then introduce two implementations of this model, LISANode and LISA Instrument. We demonstrate that TDI processing successfully recovers gravitational-wave signals from the significantly more realistic and complex simulated data. LISANode and LISA Instrument are already widely used by the LISA community and, for example, currently provide the mock data for the LISA Data Challenges.

We calculate the expectation value of the energy density and pressure of a scalar field after its decoupling from a thermal bath in the spatially flat Friedman-Robertson-Walker space-time, within the framework of quantum statistical mechanics. By using the density operator determined by the condition of local thermodynamic equilibrium, we determine the mean value of the stress-energy tensor of a real scalar field by subtracting the vacuum expectation value at the time of the decoupling. The obtained expressions of energy density and pressure involve corrections with respect to the classical free-streaming solution of the relativistic Boltzmann equation, which may become relevant even at long times.

We consider numerical black hole solutions in the Weyl conformal geometry, and its associated conformally invariant Weyl quadratic gravity. In this model Einstein gravity (with a positive cosmological constant) is recovered in the spontaneously broken phase of Weyl gravity, after the Weyl gauge field ($\omega _{\mu}$) becomes massive through a Stueckelberg mechanism, and it decouples. As a first step in our investigations we write down the conformally invariant gravitational action, containing a scalar degree of freedom, and the Weyl vector. The field equations are derived from the variational principle in the absence of matter. By adopting a static spherically symmetric geometry, the vacuum field equations for the gravitational, scalar, and Weyl fields are obtained. After reformulating the field equations in a dimensionless form, and by introducing a suitable independent radial coordinate, we obtain their solutions numerically. We detect the formation of a black hole from the presence of a Killing horizon for the timelike Killing vector in the metric tensor components, indicating the existence of the singularity in the metric. Several models, corresponding to different functional forms of the Weyl vector, are considered. An exact black hole model, corresponding to a Weyl vector having only a radial spacelike component, is also obtained. The thermodynamic properties of the Weyl geometric type black holes (horizon temperature, specific heat, entropy and evaporation time due to Hawking luminosity) are also analyzed in detail.

Can modern cosmological observations be reconciled with a general-relativistic Universe without an anti-gravitating energy source? Usually, the answer to this question by cosmologists is in the negative, and it is commonly believed that the observed excess dimming of supernovae relative to that in the Milne model is evidence for dark energy. In this paper, we develop a theorem that clarifies the conditions for such an excess dimming, based on which we argue that the answer to the above question may counter-intuitively be `yes'.

Determinants of the second-rank tensors stand useful in forming generally invariant terms as in the case of the volume element of the gravitational actions. Here, we extend the action of the matter fields by an arbitrary function $f(D)$ of the determinants of their energy-momentum, and the metric, $D=|\textbf{det}.T|/|\textbf{det}.g|$. We derive the gravitational field equations and examine the nonlinear terms induced by the determinant, specifically, the inverse of the energy-momentum tensor. We also show that these extensions require a nonzero stress-energy tensor for the vacuum. We propose a scale-free model, $f(D)=\lambda D^{1/4}$, and show how it induces the familiar invariant terms formed by the trace of the energy-momentum tensor by expanding the action around the stress-energy of the vacuum. We study the hydrostatic equilibrium equations for a neutron star by providing relevant values of the dimensionless constant $\lambda$. We show that the differences from the predictions of general relativity, in the mass-radius relations, which are sensitive to the equations of state, are conspicuous for $\lambda \sim -10^{-2}$. We also show that the model does not affect the predictions on the primordial nucleosynthesis when it is applied to the early radiation era. This novel and unfamiliar type of gravity-matter coupling can lead to a rich phenomenology in gravitational physics.

We combine experimental constraints from direct detection searches and from neutrino telescopes looking for WIMP annihilations in the Sun to derive halo-independent bounds on each of the 28 WIMP-proton and WIMP-neutron couplings of the effective non-relativistic Hamiltonian that drives the scattering process off nuclei of a WIMP of spin 1/2. The method assumes that the velocity distribution is normalized to one and homogeneous at the the solar system scale, as well as equilibrium between WIMP capture and annihilation in the Sun, and requires to fix the WIMP annihilation channels (we assume $b\bar{b}$). We consider a single non-vanishing coupling at a time, and find that for most of the couplings the degree of relaxation of the halo-independent bounds compared to those obtained by assuming the Standard Halo Model is with few exceptions relatively moderate in the low and high WIMP mass regimes, where it can be as small as a factor of $\simeq 2$, while in the intermediate mass range between 10 GeV and 200 GeV it can be as large as $\sim 10^3$. An exception to this general pattern, with more moderate values of the bound relaxation, is observed in the case of spin-dependent WIMP-proton couplings with no or a comparatively small momentum suppression, for which WIMP capture is strongly enhanced because it is driven by scattering events off $^1H$, which is the most abundant target in the Sun. Within this class of operators the relaxation is particularly small for interactions that are driven by only the velocity-dependent term, for which the solar capture signal is enhanced compared to the direct detection one, thanks to the highest speed of scattering WIMPs within the Sun due to the larger gravitational acceleration.

Horizon thermodynamics is expected to be related to the effective energy based on the energy density calculated from the Friedmann equation for a Friedmann--Robertson--Walker (FRW) universe. In the present study, the effective energy and thermostatistical quantities on a cosmological horizon are examined to clarify the holographic-like connection between them, with a focus on a de Sitter universe. To this end, the Helmholtz free energy on the horizon is derived from horizon thermodynamics. The free energy is found to be equivalent to the effective energy calculated from the Friedmann equation. This consistency is interpreted as a kind of holographic-like connection. To examine this connection, Padmanabhan's holographic equipartition law, which is related to the origin of spacetime dynamics, is applied to a de Sitter universe. It is found that the law should lead to a holographic-like connection. The holographic-like connection is considered to be a bridge between thermostatistical quantities on the horizon and in the bulk. For example, cosmological equations for a flat FRW universe can be derived from horizon thermodynamics by accepting the connection as a viable scenario. In addition, a thermal entropy equivalent to the Bekenstein--Hawking entropy is obtained from the Friedmann equation using the concept of a canonical ensemble in statistical physics. The present study provides new insight into the discussion of horizon thermodynamics and cosmological equations.

We investigate the use of spinors to describe the secular evolution of quasi-Keplerian systems. Evaluating their Poisson brackets, we show that the components of a properly-chosen spinor are canonical variables. We illustrate this formalism with a satellite's motion around an oblate body.

We study the propagation of cosmological gravitational wave (GW) backgrounds from the early radiation era until the present day in modified theories of gravity. Comparing to general relativity (GR), we study the effects that Horndeski parameters, such as the run rate of the effective Planck mass $\alpha_{\rm M}$ and the tensor speed excess $\alpha_{\rm T}$, have on the present-day GW spectrum. We use both the WKB estimate, which provides an analytical description but fails at superhorizon scales, and numerical simulations that allow us to go beyond the WKB approximation. We show that $\alpha_{\rm T}$ makes relatively insignificant changes to the GR solution, especially taking into account the constraints on its value from GW observations by the LIGO-Virgo collaboration, while $\alpha_{\rm M}$ can introduce modifications to the spectral slopes of the GW energy spectrum in the low-frequency regime depending on the considered time evolution of $\alpha_{\rm M}$. The latter effect is additional to the damping or growth occurring equally at all scales that can be predicted by the WKB approximation. In light of the recent observations by pulsar timing array collaborations and future detectors such as SKA, LISA, DECIGO, BBO, or ET, we show that, in most of the cases, constraints can not be placed on the effects of $\alpha_{\rm M}$ and the initial GW energy density $\mathcal{E}_{\rm GW}^*$ separately, but only on the combined effects of the two.

We review various aspects of de Sitter spacetime in string theory: its status as an effective field theory spacetime solution, its relation to the vacuum energy problem in string theory, its (global) holographic definition in terms of two entangled and non-canonical conformal field theories, as well as a realization of a realistic de Sitter universe endowed with the observed visible matter and the necessary dark sector in order to reproduce the realistic cosmological structure. In particular, based on the new insight regarding the cosmological constant problem in string theory, we argue that in a doubled, T-duality-symmetric, phase-space-like and non-commutative generalized-geometric formulation, string theory can naturally lead to a small and positive cosmological constant that is radiatively stable and technically natural. Such a formulation is fundamentally based on a quantum spacetime, but in an effective spacetime description of this general formulation of string theory, the curvature of the dual spacetime is the cosmological constant of the observed spacetime, while the size of the dual spacetime is the gravitational constant of the same observed spacetime. Also, the three scales associated with intrinsic non-commutativity of string theory, the cosmological constant scale and the Planck scale, as well as the Higgs scale, can be arranged to satisfy various seesaw-like formulae. Along the way, we show that these new features of string theory can be implemented in a particular deformation of cosmic-string-like models.

We use physics-informed neural networks (PINNs) to compute the first quasi-normal modes of the Kerr geometry via the Teukolsky equation. This technique allows us to extract the complex frequencies and separation constants of the equation without the need for sophisticated numerical techniques, and with an almost immediate implementation under the \texttt{PyTorch} framework. We are able to compute the oscillation frequencies and damping times for arbitrary black hole spins and masses, with accuracy typically below the percentual level as compared to the accepted values in the literature. We find that PINN-computed quasi-normal modes are indistinguishable from those obtained through existing methods at signal-to-noise ratios (SNRs) larger than 100, making the former reliable for gravitational-wave data analysis in the mid term, before the arrival of third-generation detectors like LISA or the Einstein Telescope, where SNRs of ${\cal O}(1000)$ might be achieved.

A scenario of the inflaton evolution from cosmic inflation and matter creation to dark energy/dark matter today is presented. To start with, a model of the inflationary phase of the inflaton is discussed. The inflaton rolls down a simple quadratic hilltop potential along with matter creation, following an exact tracking solution of its dynamics. Being dragged down by the presence of matter, it rolls down slowly and naturally ending inflation as the universe stops accelerating due to the presence of matter. The model predictions for the standard metrics such as scalar/tensor spectral indexes and tensor to scalar ratio are fully consistent with the current CMB limits. The quadratic potential could be extended to complete a potential hill subsequent to inflation. The evolution of inflaton discussed recently is consistent with this picture of its journey from cosmic inflation to dark matter today.