Papers for Tuesday, Mar 12 2024

Polycyclic aromatic hydrocarbon (PAH) molecules have long been suggested to be present in the interstellar medium (ISM). Nevertheless, despite their expected ubiquity and sustained searching efforts, identifying specific interstellar PAH molecules from their infrared (IR) spectroscopy has so far been unsuccessful. However, due to its unprecedented sensitivity, the advent of the James Webb Space Telescope (JWST) may change this. Meanwhile, recent years have witnessed breakthroughs in detecting specific PAH molecules (e.g., indene, cyanoindene, and cyanonaphthalene) through their rotational lines in the radio frequencies. As JWST holds great promise for identifying specific PAH molecules in the ISM based on their vibrational spectra in the IR, in this work we model the vibrational excitation of indene, a molecule composed of a six-membered benzene ring fused with a five-membered cyclopentene ring, and calculate its IR emission spectra for a number of representative astrophysical regions. This will facilitate JWST to search for and identify indene in space through its vibrational bands and to quantitatively determine or place an upper limit on its abundance.

Super-puffs are planets with exceptionally low densities ($\rho \lesssim 0.1$~g~cm$^{-3}$) and core masses ($M_c \lesssim 5 M_\oplus$). Many lower-mass ($M_p\lesssim10M_\oplus$) super-puffs are expected to be unstable to catastrophic mass loss via photoevaporation and/or boil-off, whereas the larger gravitational potentials of higher-mass ($M_p\gtrsim10M_\oplus$) super-puffs should make them more stable to these processes. We test this expectation by studying atmospheric loss in the warm, higher-mass super-puff TOI-1420b ($M = 25.1M_\oplus$, $R = 11.9R_\oplus$, $\rho = 0.08$~g~cm$^{-3}$, $T_\mathrm{eq} = 960$~K). We observed one full transit and one partial transit of this planet using the metastable helium filter on Palomar/WIRC and found that the helium transits were $0.671\pm0.079\%$ (8.5$\sigma$) deeper than the TESS transits, indicating an outflowing atmosphere. We modeled the excess helium absorption using a self-consistent 1D hydrodynamics code to constrain the thermal structure of the outflow given different assumptions for the stellar XUV spectrum. These calculations then informed a 3D simulation which provided a good match to the observations with a modest planetary mass-loss rate of $10^{10.82}$~g~s$^{-1}$ ($M_p/\dot{M}\approx70$~Gyr). Super-puffs with $M_p\gtrsim10M_\oplus$, like TOI-1420b and WASP-107b, appear perfectly capable of retaining atmospheres over long timescales; therefore, these planets may have formed with the unusually large envelope mass fractions they appear to possess today. Alternatively, tidal circularization could have plausibly heated and inflated these planets, which would bring their envelope mass fractions into better agreement with expectations from core-nucleated accretion.

When orbiting hotter stars, hot Jupiters are often highly inclined relative to their host star equator planes. By contrast, hot Jupiters orbiting cooler stars are more aligned. Prior attempts to explain this correlation between stellar obliquity and effective temperature have proven problematic. We show how resonance locking -- the coupling of the planet's orbit to a stellar gravity mode (g mode) -- can solve this mystery. Cooler stars with their radiative cores are more likely to be found with g-mode frequencies increased substantially by core hydrogen burning. Strong frequency evolution in resonance lock drives strong tidal evolution; locking to an axisymmetric g mode damps semi-major axes, eccentricities, and as we show for the first time, obliquities. Around cooler stars, hot Jupiters evolve into spin-orbit alignment and avoid engulfment. Hotter stars lack radiative cores, and therefore preserve congenital spin-orbit misalignments. We focus on resonance locks with axisymmetric modes, supplementing our technical results with simple physical interpretations, and show that non-axisymmetric modes also damp obliquity.

The pervasive presence of warm gas in galaxy halos suggests that the circumgalactic medium (CGM) is multiphase in its ionization structure and complex in its kinematics. Some recent state-of-the-art cosmological galaxy simulations predict an azimuthal dependence of CGM metallicities. We investigate the presence of such a trend by analyzing the distribution of gas properties in the CGM around 47 $z <$ 0.7 galaxies from the Multiphase Galaxy Halos Survey determined using a cloud-by-cloud, multiphase, ionization modelling approach. We identify three distinct populations of absorbers: cool clouds ($T \sim$ 10$^{4.1}$ K) in photoionization equilibrium, warm-hot collisionally ionized clouds ($T \sim$ 10$^{4.5-5}$ K) affected by time-dependent photoionization, and hotter clouds ($T \sim$ 10$^{5.4-6}$ K) with broad OVI and Lya absorption consistent with collisional ionization. We find that fragmentation can play a role in the origin of cool clouds, that warm-hot clouds are out of equilibrium due to rapid cooling, and that hotter clouds are representative of virialized halo gas in all but the lowest mass galaxies. The metallicities of clouds do not depend on the azimuthal angle or other galaxy properties for any of these populations. At face value, this disagrees with the simplistic model of the CGM with bipolar outflows and cold-mode planar accretion. However, the number of clouds per sightline is significantly larger close to the minor and major axes. This implies that the processes of outflows and accretion are contributing to these CGM cloud populations, but that they are well mixed in these low redshift galaxies.

Cosmological observables are particularly sensitive to key ratios of energy densities and rates, both today and at earlier epochs of the Universe. Well-known examples include the photon-to-baryon and the matter-to-radiation ratios. Equally important, though less publicized, are the ratios of pressure-supported to pressureless matter and the Thomson scattering rate to the Hubble rate around recombination, both of which observations tightly constrain. Preserving these key ratios in theories beyond the $\Lambda$ Cold-Dark-Matter ($\Lambda$CDM) model ensures broad concordance with a large swath of datasets when addressing cosmological tensions. We demonstrate that a mirror dark sector, reflecting a partial $\mathbb{Z}_2$ symmetry with the Standard Model, in conjunction with percent level changes to the visible fine-structure constant and electron mass which represent a \textit{phenomenological} change to the Thomson scattering rate, maintains essential cosmological ratios. Incorporating this ratio preserving approach into a cosmological framework significantly improves agreement to observational data ($\Delta\chi^2=-35.72$) and completely eliminates the Hubble tension with a cosmologically inferred $H_0 = 73.80 \pm 1.02$ km/s/Mpc when including the S$H_0$ES calibration in our analysis. While our approach is certainly nonminimal, it emphasizes the importance of keeping key ratios constant when exploring models beyond $\Lambda$CDM.

We introduce the Hawai`i Supernova Flows project and present summary statistics of the first 1218 astronomical transients observed, 669 of which are spectroscopically classified Type Ia Supernovae (SNe Ia). Our project is designed to obtain systematics-limited distances to SNe Ia while consuming minimal dedicated observational resources. This growing sample will provide increasing resolution into peculiar velocities as a function of position on the sky and redshift, allowing us to more accurately map the structure of dark matter. This can be used to derive cosmological parameters such as $\sigma_8$ and can be compared with large scale flow maps from other methods such as luminosity-line width or luminosity-velocity dispersion correlations in galaxies. Additionally, our photometry will provide a valuable test bed for analyses of SNe Ia incorporating near-infrared data. In this survey paper, we describe the methodology used to select targets, collect and reduce data, and calculate distances.

We present a method that calibrates a semi-analytic model to the Renaissance Simulations, a suite of cosmological hydrodynamical simulations with high-redshift galaxy formation. This approach combines the strengths of semi-analytic techniques and hydrodynamical simulations, enabling the extension to larger volumes and lower redshifts inaccessible to simulations due to computational expense. Using a sample of Renaissance star formation histories (SFHs) from an average density region of the Universe, we construct a four parameter prescription for metal-enriched star formation characterized by an initial bursty stage followed by a steady stage where stars are formed at constant efficiencies. Our model also includes a treatment of Population III star formation where a minimum halo mass and log-normal distribution of stellar mass are adopted to match the numerical simulations. Star formation is generally well reproduced for halos with masses $\lesssim$$10^{9} M_{\mathrm{\odot}}$. Between $11<z<25$ our model produces metal-enriched star formation rate densities (SFRDs) that typically agree with Renaissance within a factor of $\sim$2 for the average density region. Additionally, the total metal-enriched stellar mass only differs from Renaissance by about $10\%$ at $z \sim 11$. For regions that are either more overdense or rarefied not included in the calibration, we produce metal-enriched SFRDs that agree with Renaissance within a factor of $\sim$2 at high-$z$, but eventually differ by higher factors for later times. This is likely due to environmental dependencies not included in the model. Our star formation prescriptions can easily be adopted in other analytic or semi-analytic works to match our calibration to Renaissance.

We quantify for the first time the gravitational wave (GW) phase shift appearing in the waveform of eccentric binary black hole (BBH) mergers formed dynamically in three-body systems. For this, we have developed a novel numerical method where we construct a reference binary, by evolving the post-Newtonian (PN) evolution equations backwards from a point near merger without the inclusion of the third object, that can be compared to the real binary that evolves under the influence from the third BH. From this we quantify how the interplay between dynamical tides, PN-effects, and the time-dependent Doppler shift of the eccentric GW source results in unique observable GW phase shifts that can be mapped to the gravitational dynamics taking place at formation. We further find a new analytical expression for the GW phase shift, which surprisingly has a universal functional form that only depends on the time-evolving BBH eccentricity. The normalization scales with the BH masses and initial separation, which can be linked to the underlying astrophysical environment. GW phase shifts from a chaotic 3-body BH scattering taking place in a cluster, and from a BBH inspiraling in a disk migration trap near a super-massive BH, are also shown for illustration. When current and future GW detectors start to observe eccentric GW sources with high enough signal-to-noise-ratio, we propose this to be among the only ways of directly probing the dynamical origin of individual BBH mergers using GWs alone.

We present new simulations investigating the impact of mass transfer on the asteroseismic signals of slowly pulsating B stars. We use MESA to simulate the evolution of a binary star system and GYRE to compute the asteroseismic properties of the accretor star. We show that, compared to a single star of the same final mass, a star that has undergone accretion (of non-enriched material) has a significantly different internal structure, evident in both the hydrogen abundance profile and Brunt-V\"ais\"al\"a frequency profile. These differences result in significant changes in the observed period spacing patterns, implying that one may use this as a diagnostic to test whether a star's core has been rejuvenated as a result of accretion. We show that it is essential to consider the full multimodal posterior distributions when fitting stellar properties of mass-gainers to avoid drawing misleading conclusions. Even with these considerations, stellar ages will be significantly underestimated when assuming single star evolution for a mass-gainer. We find that future detectors with improved uncertainties would rule out single star models with the correct mass and central hydrogen fraction. Our proof of principle analysis demonstrates the need to further investigate the impact of binary interactions on stellar asteroseismic signals for a wide range of parameters, such as initial mass, amount of mass transferred and the age of the accretor star at the onset of mass transfer.

We analyze accretion-rate time series for equal-mass binaries in co-planar gaseous disks spanning a continuous range of orbital eccentricities up to 0.8, for both prograde and retrograde systems. The dominant variability timescales match that of previous investigations; the binary orbital period is dominant for prograde binaries with $e \gtrsim 0.1$, with a 5 times longer "lump" period taking over for $e\lesssim 0.1$. This lump period fades and drops from 5 times to 4.5 times the binary period as $e$ approaches 0.1, where it vanishes. For retrograde orbits, the binary orbital period dominates at $e \lesssim 0.55$ and is accompanied by a 2 times longer-timescale periodicity at higher eccentricities. The shape of the accretion-rate time series varies with binary eccentricity. For prograde systems, the orientation of an eccentric disk causes periodic trading of accretion between the binary components in a ratio that we report as a function of binary eccentricity. We present a publicly available tool, binlite, that can rapidly ($\lesssim 0.01$~sec) generate templates for the accretion-rate time series, onto either binary component, for choice of binary eccentricity below 0.8. As an example use-case, we build lightcurve models where the accretion rate through the circumbinary disk and onto each binary component sets contributions to the emitted specific flux. We combine these rest-frame, accretion-variability lightcurves with observer-dependent Doppler boosting and binary self-lensing. This allows a flexible approach to generating lightcurves over a wide range of binary and observer parameter space. We envision binlite as the access point to a living database that will be updated with state-of-the-art hydrodynamical calculations as they advance.

One major goal in fast radio burst science is to detect fast radio bursts (FRBs) over a wide field of view without sacrificing the angular resolution required to pinpoint them to their host galaxies. Wide-field detection and localization capabilities have already been demonstrated using connected-element interferometry; the CHIME/FRB Outriggers project will push this further using widefield cylindrical telescopes as widefield outriggers for very long baseline interferometry (VLBI). This paper describes an offline VLBI software correlator written in Python for the CHIME/FRB Outriggers project. It includes features well-suited to modern widefield instruments like multibeaming/multiple phase center correlation, pulse gating including coherent dedispersion, and a novel correlation algorithm based on the quadratic estimator formalism. This algorithm mitigates sensitivity loss which arises in instruments where the windowing and channelization is done outside the VLBI correlator at each station, which accounts for a 30 percent sensitivity drop away from the phase center. Our correlation algorithm recovers this sensitivity on both simulated and real data. As an end to end check of our software, we have written a preliminary pipeline for VLBI calibration and single-pulse localization, which we use in Lanman et al. (2024) to verify the astrometric accuracy of the CHIME/FRB Outriggers array.

We report on the discovery of Swift J010902.6-723710, a rare eclipsing Be/X-ray Binary system by the Swift SMC Survey (S-CUBED). Swift J010902.6-723710 was discovered via weekly S-CUBED monitoring observations when it was observed to enter a state of X-ray outburst on 10 October 2023. X-ray emission was found to be modulated by a 182s period. Optical spectroscopy is used to confirm the presence of a highly-inclined circumstellar disk surrounding a B0-0.5Ve optical companion. Historical UV and IR photometry are then used to identify strong eclipse-like features re-occurring in both light curves with a 60.623 day period, which is adopted as the orbital period of the system. Eclipsing behavior is found to be the result of a large accretion disk surrounding the neutron star. Eclipses are produced when the disk passes in front of the OBe companion, blocking light from both the stellar surface and circumstellar disk. This is only the third Be/X-ray Binary to have confirmed eclipses. We note that this rare behavior provides an important opportunity to constrain the physical parameters of a Be/X-ray Binary with greater accuracy than is possible in non-eclipsing systems.

Giant planets are expected to form within circumplanetary disks, which shape their formation history and the local environment. Here, we consider the formation and structure of circumplanetary disks that arise during the late stages of giant planet formation. During this phase, when most of the final mass is accumulated, incoming material enters the Hill sphere and falls toward the planet. In the absence of torques, the falling parcels of gas conserve their specific angular momentum and collect into a circumplanetary disk. Generalizing previous work, we consider a range of possible geometries for the flow entering the sphere of influence of the planet. Specifically, we consider five geometric patterns for the inward flow, ranging from concentration toward the rotational poles of the system to isotropic flow to concentration along the equatorial plane. For each case, we derive analytic descriptions for the density field of the infall region, the disk surface density in the absence of viscosity, and steady-state solutions for viscous disks. These results, in turn, specify the luminosity contributions of the planet, the circumplanetary disk, and the envelope. These power sources, in conjunction with the surrounding material, collectively determine the observational appearance of the forming planet. We conclude with an approximate determination of these radiative signatures.

This paper presents the multi-scale temperature structures in the Milky Way (MW) hot gas, as part of the XMM-Newton Line Emission Analysis Program (X-LEAP), surveying the O VII, O VIII, and Fe-L band emission features in the XMM-Newton archive. In particular, we define two temperature tracers, $I_{\rm OVIII}/I_{\rm OVII}$ (O87) and $I_{\rm FeL}/(I_{\rm OVII}+I_{\rm OVIII})$ (FeO). These two ratios cannot be explained simultaneously using single-temperature collisional ionization models, which indicates the need for multi-temperature structures in hot gas. In addition, we show three large-scale features in the hot gas: the eROSITA bubbles around the Galactic center (GC); the disk; and the halo. In the eROSITA bubbles, the observed line ratios can be explained by a log-normal temperature distribution with a median of $\log T/{\rm K} \approx 6.4$ and a scatter of $\sigma_T \approx 0.2$ dex. Beyond the bubbles, the line ratio dependence on the Galactic latitude suggests higher temperatures around the midplane of the MW disk. The scale height of the temperature variation is estimated to be $\approx$2 kpc assuming an average distance of $5$ kpc for the hot gas. The halo component is characterized by the dependence on the distance to the GC, showing a temperature decline from $\log\,T/{\rm K}\,\approx\, 6.3$ to $5.8$. Furthermore, we extract the auto-correlation and cross-correlation functions to investigate the small-scale structures. O87 and FeO ratios show a consistent auto-correlation scale of $\approx$$ 5^\circ$ (i.e., $\approx$$ 400$ pc at 5 kpc), which is consistent with expected physical sizes of X-ray bubbles associated with star-forming regions or supernova remnants. Finally, we examine the cross-correlation between the hot and UV-detected warm gas, and show an intriguing anti-correlation.

The true nature of primordial magnetic fields (PMFs) and their role in the formation of galaxies still remains elusive. To shed light on these unknowns, we investigate their impact by varying two sets of properties: (i) accounting for the effect of PMFs on the initial matter power spectrum, and (ii) accounting for their magneto-hydrodynamical effects on the formation of galaxies. By comparing both we can determine the dominant agent in shaping galaxy evolution. We use the magneto-hydrodynamics code RAMSES, to generate multiple zoom-in simulations for eight different host halos of dwarf galaxies across a wide luminosity range of $10^3-10^6\,L_{\odot}$. We explore a variety of primordial magnetic field (comoving) strengths ranging from $0.05$ to $0.50\,\mathrm{nG}$. We find magnetic fields in the interstellar medium not only modify star formation in dwarf spheroidal galaxies but also completely prevent the formation of stars in less compact ultra-faints with halo mass and stellar mass below $\sim 2.5\cdot10^9$ and $3\cdot10^6\,M_{\odot}$, respectively. At high redshifts, the impact of PMFs on host halos of dwarf galaxies through the modification of the matter power spectrum is more dominant than the influence of magneto-hydrodynamics in shaping their gaseous structure. Through the amplification of small perturbations ranging in mass from $10^7$ to $10^9\,M_{\odot}$ in the $\Lambda$CDM$+$PMFs matter power spectrum, primordial fields expedite the formation of the first dark matter halos, leading to an earlier onset and a higher star formation rate at redshifts $z>12$.

A new fine grid of nonlinear convective pulsation models for the so-called "bump Cepheids" is presented to investigate the Hertzprung progression (HP) phenomenon shown by their light and radial pulsation velocity curves. The period corresponding to the center of the HP is investigated as a function of various model assumptions, such as the efficiency of super-adiabatic convection, the mass-luminosity relation, and the metal and helium abundances. The assumed mass-luminosity relation is found to significantly affect the phenomenon but variations in the chemical composition as well as in the stellar mass (at fixed mass-luminosity relation) also play a key role in determining the value of the HP center period. Finally, the predictive capability of the presented theoretical scenario is tested against observed light curves of bump Cepheids in the ESA Gaia database, also considering the variation of the pulsation amplitudes and of the Fourier parameters $R_{21}$ and $\Phi_{21}$ with the pulsation period. A qualitative agreement between theory and observations is found for what concerns the evolution of the light curve morphology as the period moves across the HP center, as well for the pattern in period-amplitude, period-$R21$ and period-$\Phi_{21}$ planes. A larger sample of observed Cepheids with accurate light curves and metallicities is required in order to derive more quantitative conclusions.

Starburst galaxies are undergoing intense episodes of star formation. In these galaxies, gas is ejected into the surrounding environment through winds created by the effect of hot stars and supernova explosions. When interacting with the intergalactic medium, these winds can produce strong shocks capable of accelerating cosmic rays. The radiation from these cosmic rays mainly occurs in radio and gamma rays. The radio halo can be characterized using the scale height. We searched for the presence of radio halos in a sample of edge-on starburst galaxies gathered from the MeerKAT 1.28 GHz Atlas of Southern Sources in the IRAS Revised Bright Galaxy Sample. We selected a sample of 25 edge-on galaxies from the original sample and modeled their disk and halo contributions. We have detected and characterized 11 radio halos, seven of which are reported here for the first time. We found that the halo scale heights increase linearly with the radio diameters and this relation does not depend on the star formation rate. All galaxies in our sample follow the radio-infrared relation with a q parameter value of $2.5\pm0.1$. The dependence of the halo luminosity on the star formation rate and the infrared luminosity supports the hypothesis that the radio halos are the result of synchrotron radiation produced by relativistic electrons and points toward the fact that the star formation activity plays a crucial role in halo creation. The average scale height of 1 kpc implies a dynamical range of 4 Myr, several orders of magnitude greater than the synchrotron losses for electrons of 10 TeV. This suggests that some process must exist to reaccelerate cosmic rays in the halo if gamma-ray emission of a leptonic origin is detected from the halo. According to the relation between the radio and gamma-ray luminosities, we found that NGC 4666 is a potential gamma-ray source for future observations.

We conducted an all-sky imaging transient search with the Owens Valley Radio Observatory Long Wavelength Array (OVRO-LWA) data collected during the Perseid meteor shower in 2018. The data collection during the meteor shower was motivated to conduct a search for intrinsic radio emission from meteors below 60 MHz known as the meteor radio afterglows (MRAs). The data collected were calibrated and imaged using the core array to obtain lower angular resolution images of the sky. These images were input to a pre-existing LWA transient search pipeline to search for MRAs as well as cosmic radio transients. This search detected 5 MRAs and did not find any cosmic transients. We further conducted peeling of bright sources, near-field correction, visibility differencing and higher angular resolution imaging using the full array for these 5 MRAs. These higher angular resolution images were used to study their plasma emission structures and monitor their evolution as a function of frequency and time. With higher angular resolution imaging, we resolved the radio emission size scales to less than 1 km physical size at 100 km heights. The spectral index mapping of one of the long duration event showed signs of diffusion of plasma within the meteor trails. The unpolarized emission from the resolved radio components suggest resonant transition radiation as the possible radiation mechanism of MRAs.

Spatial template is important to study the nearby supernova remnant (SNR). For SNR G332.5-5.6, we report a gaussian disk with radius of about 1.06 degrees to be a potential good spatial model in the gamma-ray band. Employing this new gaussian disk, its GeV lightcurve shows a significant variability of about 7 sigma. The $\gamma$-ray observations of this SNR could be explained well either by a leptonic model or a hadronic model, in which a flat spectrum for the ejected electrons/protons.

The interaction between a Coronal Mass Ejection (CME) and a comet has been observed several times by in-situ observations from the Rosetta Plasma Consortium (RPC), which is designed to investigate the cometary magnetosphere of comet 67P/Churyumov-Gerasimenko (CG). Goetz et al. (2019) reported a magnetic field of up to 300 nT measured in the inner coma, which is among the largest interplanetary magnetic fields observed in the solar system. They suggested the large magnetic field observations in the inner coma come from magnetic field pile-up regions, which are generated by the interaction between a CME and/or corotating interaction region and the cometary magnetosphere. However, the detailed interaction between a CME and the cometary magnetosphere of comet CG in the inner coma has not been investigated by numerical simulations yet. In this manuscript, we will use a numerical model to simulate the interaction between comet CG and a Halloween class CME and investigate its magnetospheric response to the CME. We find that the plasma structures change significantly during the CME event, and the maximum value of the magnetic field strength is more than 500nT close to the nucleus. Virtual satellites at similar distances as Rosetta show that the magnetic field strength can be as large as 250nT, which is slightly less than what Goetz et al. (2019) reported.

Modeling the brightness of satellites in large Low-Earth Orbit (LEO) constellations can not only assist the astronomical community in assessing the impact of reflected light from satellites, optimizing observing schedules and guiding data processing, but also motivate satellite operators to improve their satellite designs, thus facilitating cooperation and consensus among different stakeholders. This work presents a photometric model of the Starlink satellites based on the Bidirectional Reflectance Distribution Function (BRDF) using millions of photometric observations. To enhance model accuracy and computational efficiency, data filtering and reduction are employed, and chassis blocking on the solar array and the earthshine effect are taken into account. The assumptions of the model are also validated by showing that the satellite attitude is as expected, the solar array is nearly perpendicular to the chassis, and both the solar array pseudo-specular reflection and the chassis earthshine should be included in the model. Reflectance characteristics of the satellites and the apparent magnitude distributions over station are finally discussed based on the photometric predictions from the model. In addition to assessing the light pollution and guiding the development of response measures, accurate photometric models of satellites can also play an important role in areas such as space situational awareness.

Context. In Gaia era, atmospheric turbulence, which causes stochastic wander of a star image, is a fundamental limitation to the astrometric accuracy of ground-based optical imaging. However, the positional bias caused by turbulence (called turbulence error here) can be effectively reduced by measuring a target relative to another reference (a star or a fast-moving target) which locates in the range of only several tens of arcsec, since they suffer from similar turbulence errors. This phenomenon is called the precision premium and has been effectively applied to the astrometry of solar system. Further investigation for the precision premium shows that, the precision premium works at less than about 100 arcsec for two specific objects and the relative positional precision as a function of their angular seperation can be well fitted by a sigmoidal function, called the precision premium curve (PPC). Aims. We want to reduce the turbulence error of a target if it is imaged in an area of high stellar density of a ground-based observation by taking advantage of more Gaia reference stars. Methods. Based on the PPC, we proposed a high-precision astrometric solution called precision premium transformation (PPT) in this paper, which takes advantage of high similarity of turbulence errors in a small region and the dense Gaia reference stars in the region to reduce the turbulence errors on the observation, through a weighted solution. Results. Through systematic analysis, the PPT method exhibits significant advantages in terms of not only precision but also applicability when a target is imaged in an area of high stellar density. The PPT method is also applied to the determination of the proper motion of an open cluster, and the results demonstrate and quantify benefits that the PPT method bestows on ground-based astrometry.

In eclipsing X-ray binary systems, the direct X-ray emission is blocked by the companion star during the eclipse. We observe only reprocessed emission that contains clues about the environment of the compact object and its chemical composition, ionization levels, etc. We have found flares in some X-ray binaries during their eclipses. The study of eclipse flares provides additional clues regarding the size of the reprocessing region and helps distinguish between different components of the X-ray spectrum observed during the eclipse. In the archival data, we searched for flares during eclipses of high-mass X-ray binaries and found flares in three sources: Vela X-1, LMC X-4, and 4U 1700-37. Comparing spectral properties of the eclipse flare and non-flare data, we found changes in the power-law photon index in all three sources and multiple emission lines in Vela X-1 and 4U 1700-37. The fluxes of prominent emission lines showed a similar increase as the overall X-ray flux during the eclipse flare, suggesting the lines originate in the binary environment and not in the interstellar medium. We also observed a soft excess in 4U 1700-37 that remains unchanged during both eclipse flare and non-flare states. Our analysis suggests that this emission originates from the extremely thin shell of the stellar wind surrounding the photosphere of its companion star. The detection of short (100-200 seconds) count-rate doubling timescale in 4U 1700-37 and LMC X-4 indicates that the eclipse reprocessing occurs in a region larger than, but comparable to the size of the companion star.

Because the conventional method of creating type-III radiation by coalescence of counter-propagating Langmuir waves has not been verified with in-situ data, Parker Solar Probe data was examined in search of such in-situ evidence, which was not found. Instead, a new mechanism for creating type-III radiation has been found as a result of observing slow Langmuir waves (~2200 km/sec) with electric fields as large as 300 mV/m during a developing type-III burst on March 21, 2023. Because of their slow phase velocities, these Langmuir waves had short wavelengths, several times the Debye length of 2.65 meters, and, as a result, k{\lambda}d~0.93. Such waves may be strongly damped to be replaced by new growing bursts of waves that create the characteristic Langmuir waveform that is composed of peaks and valleys of a few milliseconds duration. The average electron current that produced these Langmuir waves is estimated from the Generalized Ohms law to be ~15 microamps/m2 and, from the strahl, 8 microamps/m2. Peak currents were at least twice these averages. These Langmuir waves, acting as antennas, produced electrostatic harmonics having slow phase velocities (~2000 km/sec) at frequencies of n{\omega}p, where {\omega}p = 6.28*125 kHz is the Langmuir wave frequency and n = 2, 3, 4, 5, 6, and 7. Such waves are not the type-III emission. As at least the first harmonic wave evolved through the huge density irregularities, its wave number decreased and it became the electromagnetic wave that was the type-III radiation.

The ROME/REA (Robotic Observations of Microlensing Events/Reactive Event Assessment) Survey was a Key Project at Las Cumbres Observatory (hereafter LCO) which continuously monitored 20 selected fields (3.76 sq.deg.) in the Galactic Bulge throughout their seasonal visibility window over a three-year period, between March 2017 and March 2020. Observations were made in three optical passbands (SDSS-g', -r', -i'), and LCO's multi-site telescope network enabled the survey to achieve a typical cadence of $\sim$10\,hrs in i' and ~15 hrs in g' and r'. In addition, intervals of higher cadence (<1 hr) data were obtained during monitoring of key microlensing events within the fields. This paper describes the Difference Image Analysis data reduction pipeline developed to process these data, and the process for combining the photometry from LCO's three observing sites in the Southern Hemisphere. The full timeseries photometry for all 8 million stars, down to a limiting magnitude of i~18 mag is provided in the data release accompanying this paper, and samples of the data are presented for exemplar microlensing events, illustrating how the tri-band data are used to derive constraints on the microlensing source star parameters, a necessary step in determining the physical properties of the lensing object. The timeseries data also enables a wealth of additional science, for example in characterizing long-timescale stellar variability, and a few examples of the data for known variables are presented.

In this study we present the analysis of six solar flare events that occurred in 2022, using new data from the third-generation Miniature X-Ray Solar Spectrometer (MinXSS), also known as the Dual-zone Aperture X-ray Solar Spectrometer (DAXSS). The primary focus of this study is on the flare's "onset phase", which is characterized by elevated soft X-ray emissions even before the flare's impulsive phase. We analyze the temporal evolution of plasma temperature, emission measure, and elemental abundance factors during the flare onset phase, by fitting the DAXSS spectra with the Astrophysical Plasma Emission Code (APEC) model. The model fitting results indicate that the flaring-plasma is already at a high temperature (10-15 MK) during the onset period. The temperature rises during the onset phase, followed by a decrease and subsequent increase during the impulsive phase. Elemental abundance factors show a trend of falling below pre-flare values during the onset phase, with some recovery before the impulsive phase. During the impulsive phase, the abundance factors decrease from elevated coronal values to about photospheric values. We also analyze images from the 193 Angstrom channel of the Atmospheric Imaging Assembly (AIA), highlighting the formation or brightening of coronal loop structures during the onset phase. Two distinct onset loop configurations are observed which are referred to as 1-loop and 2-loop onsets. Both DAXSS and AIA observations indicate that the flare onset phase exhibits similar hot coronal plasma properties as the impulsive phase, suggesting that the onset phase may act as a preconditioning effect for some flares.

We report a detection of GeV $\gamma$-ray emission potentially originating from the pulsar wind nebula in CTA 1 by analyzing about 15 yr of Fermi Large Area Telescope data. By selecting an energy range from 50 GeV to 1 TeV to remove contamination from the $\gamma$-ray pulsar PSR J0007+7303, we have discovered an extended $\gamma$-ray source with a TS value of $\sim$ 44.94 in the region of CTA 1. The obtained flux is measured to be 6.71 $\pm$ 2.60 $\times$ $10^{-12}$ erg $\mathrm{cm}^{-2}$ $\mathrm{s}^{-1}$ with a spectral index of 1.61 $\pm$ 0.36, which allows for a smooth connection with the flux in the TeV band. CTA 1 is also considered to be associated with 1LHAASO J0007+7303u, which is an Ultra-High-Energy source listed in the recently published catalog of the Large High Altitude Air Shower Observatory. We assume that the radiation originates from the pulsar wind nebula and that its multi-wavelength spectral energy distribution can be explained well with a time-dependent one-zone model.

We consider the broadband spectral energy distribution of the high energy peaked (HBL) blazar Mkn 501 using $\textit{Swift}$-XRT/UVOT, NuSTAR and $\textit{Fermi}$-LAT observations taken between 2013 and 2022. The spectra were fitted with a one-zone leptonic model using synchrotron and synchrotron self-Compton emission from different particle energy distributions such as a broken power-law, log-parabola, as well as distributions expected when the diffusion or the acceleration time scale are energy dependent. The jet power estimated for a broken power-law distribution was $ \sim 10^{47} (10^{44})$ erg s$^{-1}$ for a minimum electron energy $\gamma_\text{min} \sim 10 (10^3)$. However, for electron energy distributions with intrinsic curvature (such as the log-parabola form), the jet power is significantly lower at a few times $ 10^{42}$ erg s$^{-1}$ which is a few percent of the Eddington luminosity of a $10^7$ M$_\odot$ black hole, suggesting that the jet may be powered by accretion processes. We discuss the implications of these results.

We use a Kernel Density Estimator method to evaluate the stellar velocity dispersion in the open cluster NGC 2571. We derive the 3-D velocity dispersion using both proper motions as extracted from Gaia DR3 and single epoch radial velocities as obtained with the instrument FLAMES at ESO VLT. The mean-square velocity along the line-of-sight is found to be larger than the one in the tangential direction by a factor in the interval [6,8]. We argue that the most likely explanation for such an occurrence is the presence of a significant quantity of unresolved binary and multiple stars in the radial velocity sample. Special attention should be paid to single line spectroscopic binaries (SB1) since in this case we observe the spectral lines of the primary component only, and therefore the derived radial velocity is not the velocity of the binary system center of mass. To investigate this scenario, we performed numerical experiments at varying the fractional abundance of SB1 in the observed sample. These experiments show that the increase of the mean-square radial velocity depends actually on the fractional abundance of SB1 to a power in the range of [0.39,0.45]. We used the 3-D velocity dispersion obtained by the dispersions in the tangential directions and the assumption that the radial velocity dispersion is the same as a tangential one to estimate the virial cluster mass and the cluster mass taking into account the gravitational field of the Galaxy and the non-stationarity of the cluster. These estimates are $650\pm30 \; M_\odot$ and $310\pm80 \; M_\odot$, respectively, and they are in substantial agreement with the photometric cluster mass.

The solar dynamo is essentially a cyclic process in which the toroidal component of the magnetic field is converted into the poloidal one and vice versa. This cyclic loop is disturbed by some nonlinear and stochastic processes mainly operating in the toroidal to poloidal part. Hence, the memory of the polar field decreases in every cycle. On the other hand, the dynamo efficiency and, thus, the supercriticality of the dynamo decreases with the Sun's age. Previous studies have shown that the memory of the polar magnetic field decreases with the increase of supercriticality of the dynamo. In this study, we employ popular techniques of time series analysis, namely, compute Higuchi's fractal dimension, Hurst exponent, and Multi-Fractal Detrended Fluctuation Analysis, to the amplitude of the solar magnetic cycle obtained from dynamo models operating at near-critical and supercritical regimes. We show that the magnetic field in the near-critical regime is governed by strong memory, less stochasticity, intermittency, and breakdown of self-similarity. On the contrary, the magnetic field in the supercritical region has less memory, strong stochasticity, and shows a good amount of self-similarity. Finally, applying the same time series analysis techniques in the reconstructed sunspot data of 85 cycles and comparing their results with that from models, we conclude that the solar dynamo is possibly operating near the critical regime and not too much supercritical regime. Thus Sun may not be too far from the critical dynamo transition.

The estimation of the bulge and disk massses, the main baryonic components of a galaxy, can be performed using various approaches, but their implementation tend to be challenging as they often rely on strong assumptions about either the baryon dynamics or the dark matter model. In this work, we present an alternative method for predicting the masses of galactic components, including the disk, bulge, stellar and total mass, using a set of machine learning algorithms: KNN-neighbours (KNN), Linear Regression (LR), Random Forest (RF) and Neural Network (NN). The rest-frame absolute magnitudes in the ugriz-photometric system were selected as input features, and the training was performed using a sample of spiral galaxies hosting a bulge from Guo's mock catalogue \citep{Guo-Catalog} derived from the Millennium simulation. In general, all the algorithms provide good predictions for the galaxy's mass components ranging from $10^9\,M_\odot$ to $10^{11}\,M_\odot$, corresponding to the central region of the training mass domain; however, the NN give rise to the most precise predictions in comparison to other methods. Additionally, to test the performance of the NN architecture, we used a sample of observed galaxies from the SDSS survey whose mass components are known. We found that the NN can predict the luminous masses of disk-dominant galaxies within the same range of magnitudes that for the synthetic sample up to a $99\%$ level of confidence, while mass components of galaxies hosting larger bulges are well predicted up to $95\%$ level of confidence. The NN algorithm can also bring up scaling relations between masses of different components and magnitudes.

We study the pulsar energy-dependent $\gamma$-ray light curves and spectra from curvature radiation in the dissipative magnetospheres. The dissipative magnetospheres with the combined force-free (FFE) and Aristotelian (AE) are computed by a pseudo-spectral method with the high-resolution simulation in the rotating coordinate system, which produces a near force-free field structure with the dissipative region only near the equatorial current sheet outside the light cylinder (LC). We use the test particle trajectory method to compute the energy-dependent $\gamma$-ray light curves, phase-average and phase-resolved spectra by including both the accelerating electric field and radiation reaction. The predicted energy-dependent $\gamma$-ray light curves and spectra are then compared with those of the Vela pulsar observed by Fermi. Our results can generally reproduce the observed trends of the energy-dependent $\gamma$-ray light curves and spectra for the Vela pulsar.

We review the recent advances in the pulsar high-energy $\gamma$-ray observation and the electrodynamics of the pulsar magnetospheres from the early vacuum model to the recent plasma-filled models by the numerical simulations. The numerical simulations have made the significant progresses toward the self-consistent modeling of the plasma-filled magnetosphere by including the particle acceleration and radiation. The current numerical simulations confirm a near force-free magnetosphere with the particle acceleration in the separatrix near the light cylinder and the current sheet outside the light cylinder, which can provide a good match to the recent high-energy $\gamma$-ray observations. The modeling of the combined multi-wavelength light curves, spectra, and polarization are expected to provide a stronger constrain on the geometry of the magnetic field lines, the location of the particle acceleration and the emission region, and the emission mechanism in the pulsar magnetospheres.

We present a new method - called HINORA (HIgh-NOise RANdom SAmple Consensus) - for the identification of regular structures in 3D point distributions. Motivated by the possible existence of the so called Council of Giants, i.e. a ring of twelve massive galaxies surrounding the Local Group in the Local Sheet with a radius of 3.75 Mpc, we apply HINORA to the Local Volume Galaxy catalogue confirming its existence. When varying the lower limit of K-band luminosity of the galaxy entering the catalogue, we further report on the existence of another ring-like structure in the Local Volume that now contains the Milky Way and M31. However, this newly found structure is dominated by low-mass (satellite) galaxies. While we here simply present the novel method as well as its first application to observational data, follow-up work using numerical simulations of cosmic structure formation shall shed light into the origin of such regular patterns in the galaxy distribution. Further, the method is equally suited to identify similar (or even different) structures in various kinds of astrophysical data (e.g. locating the actual 'baryonic-acoustic oscillation spheres' in galaxy redshift surveys).

We apply the Convolutional Neural Networks (CNNs) to the mock 21cm maps from the post-reionization epoch to show that the $\Lambda$CDM and warm dark matter (WDM) model can be distinguished for WDM particle masses $m_{FD}<3$ keV, under the assumption of thermal production of WDM following the Fermi-Dirac (FD) distribution. We demonstrate that the CNN is a potent tool in distinguishing the dark matter masses, highlighting its sensitivity to the subtle differences in the 21cm maps produced by varying dark matter masses. Furthermore, we extend our analysis to encompass different WDM production mechanisms, recognizing that the dark matter production mechanism in the early Universe is among the sources of the most significant uncertainty for the dark matter model building. In this work, given the mass of the dark matter, we discuss the feasibility of discriminating four different WDM models: Fermi-Dirac (FD) distribution model, neutrino Minimal Standard Model ($\nu$MSM), Dodelson-Widrow (DW), and Shi-Fuller (SF) model. For instance, when the WDM mass is 2 keV, we show that one can differentiate between CDM, FD, $\nu$MSM, and DW models while discerning between the DW and SF models turns out to be challenging. Our results reinforce the viability of the CNN as a robust analysis for 21cm maps and shed light on its potential to unravel the features associated with different dark matter production mechanisms.

Among Neptunian mass planets exoplanets ($20-50$ M$_\oplus$), puffy hot Neptunes are extremely rare, and their unique combination of low mass and extended radii implies very low density ($\rho < 0.3$ g cm$^{-3}$). Over the last decade, only a few puffy planets have been detected and precisely characterized with both transit and radial velocity observations, most notably including WASP-107 $b$, TOI-1420 $b$, and WASP-193 $b$. In this paper, we report the discovery of TOI-1173 A $b$, a low-density ($\rho = 0.269_{-0.024}^{+0.028}$ g cm$^{-3}$) super-Neptune with $P = 7.06$ days in a nearly circular orbit around the primary G-dwarf star in the wide binary system TOI-1173 A/B. Using radial velocity observations with the MAROON-X spectrograph and transit photometry from TESS, we determined a planet mass of $M_{\rm{p}} = 26.1\pm1.9\ M_{\oplus}$ and radius of $R_{\rm{p}} = 8.10\pm0.17\ R_{\oplus}$. TOI-1173 A $b$ is the first puffy Super-Neptune planet detected in a wide binary system (separation $\sim 11,400$ AU). We explored several mechanisms to understand the puffy nature of TOI-1173 A $b$, and showed that tidal heating is the most promising explanation. Furthermore, we demonstrate that TOI-1173 A $b$ likely has maintained its orbital stability over time and may have undergone von-Zeipel-Lidov-Kozai migration followed by tidal circularization given its present-day architecture, with important implications for planet migration theory and induced engulfment into the host star. Further investigation of the atmosphere of TOI-1173 A $b$ will shed light on the origin of close-in low-density Neptunian planets in field and binary systems, while spin-orbit analyses may elucidate the dynamical evolution of the system.

Galaxy clusters as gravitational lenses play a unique role in astrophysics and cosmology: they permit mapping the dark matter distribution on a range of scales; they reveal the properties of high and intermediate redshift background galaxies that would otherwise be unreachable with telescopes; they constrain the particle nature of dark matter and are a powerful probe of global cosmological parameters, like the Hubble constant. In this review we summarize the current status of cluster lensing observations and the insights they provide, and offer a glimpse into the capabilities that ongoing, and the upcoming next generation of telescopes and surveys will deliver. While many open questions remain, cluster lensing promises to remain at the forefront of discoveries in astrophysics and cosmology.

We explore the multiwavelength radiation properties of the light curves and energy spectra in the dissipative magnetospheres of pulsars. The dissipative magnetospheres are simulated by the pseudo-spectral method with the combined force-free and Aristotelian electrodynamics, which can produce self-consistent accelerating electric fields mainly distributed in the equatorial current sheet outside the light cylinder. The multiwavelength light curves and spectra are computed by using the multiple emission mechanisms of both the primary particles accelerated by the accelerating electric fields in the equatorial current sheet and the secondary pairs with an assumed distribution spectrum. We then compare the predicted multiwavelength light curves and spectra with the observed data from the Crab, Vela, and Geminga pulsars. Our modeling results can systematically well reproduce the observed trends of the multiwavelength light curves and the spectra for these three pulsars.

This study investigates accelerated cosmic expansion using the Viscous Modified Chaplygin Gas (VMMG) and Generalized Cosmic Chaplygin Gas (GCCM) within Horava-Lifshitz gravity. Our primary objective is to constrain essential cosmological parameters, such as the Hubble Parameter ($H_{0}$) and Sound Horizon ($r_{d}$). We incorporate recent datasets comprising 17 Baryon Acoustic Oscillation observations, 33 Cosmic Chronometer measurements, 40 Type Ia Supernovae data points, 24 quasar Hubble diagram data points, and 162 Gamma Ray Bursts data points. Additionally, we integrate the most recent determination of the Hubble constant (R22). Treating $r_{d}$ as a free parameter helps alleviate bias, enhance precision, and improve dataset compatibility. By simulating random correlations in the covariance matrix, errors are reduced. The obtained values of ($H_{0}$) and ($r_{d}$) are consistent with those from Planck and SDSS. Cosmography tests offer insights into the dynamics of both models, aiding our understanding of cosmological evolution. Statefinder diagnostics deepen our understanding of cosmic expansion dynamics and aid in distinguishing between different cosmological models. The $o_{m}$ Diagnostic test reveals that at late times, VMMG falls in the phantom region and GCCM falls in the phantom quintessence region. The Akaike Information Criterion (AIC) and Bayesian Information Criterion (BIC) provide support for all models under consideration, indicating that each model offers a plausible explanation. Notably, the $\Lambda$CDM model emerges with the lowest AIC score, suggesting its relatively superior fit compared to others. Additionally, validation through the reduced $\chi_{\text{red}}^{2}$ statistic confirms satisfactory fits across all models, further reinforcing their credibility in explaining the observed data.

The solar cycle is a complex phenomenon, a comprehensive understanding of which requires the study of various tracers. Here, we consider the solar cycle as manifested in the harmonics of the solar large-scale surface magnetic field, including zonal, sectorial and tesseral harmonics, divided into odd and even relative to the solar equator. In addition to considering the amplitudes of the harmonics, we analyze their contribution to the magnetic energy. It turns out that the relative contribution of different types of harmonics to the magnetic energy is virtually independent of the cycle height. We identify different phases of the activity cycle using harmonics of different symmetries. A possible way to incorporate the obtained result into the solar dynamo theory is proposed.

We investigate the effect of a dark matter caustic passing through the Solar System. We find, confirming a previous result, that the Sun tracks the caustic surface for some time. We integrate numerically the equations of motion of the Sun and a comet for a large number of initial conditions and of caustic passage properties. We calculate the probability for the comet to escape the Solar System and the probability for it to fall within 50 A.U. of the Sun, given the initial semi-major axis and eccentricity of its orbit. We find that the average probability for a comet to fall within 50 A.U. of the Sun is of order $3 \cdot 10^{-4}$ and that comets which are initially at a distance larger than about $10^5$ A.U. have a probability of order one to be ejected from the Solar System.

Gravitational Waves (GWs) provide a unique way to explore our Universe. The ongoing ground-based detectors, e.g., LIGO, Virgo, and KAGRA, and the upcoming next-generation detectors, e.g., Cosmic Explorer and Einstein Telescope, as well as the future space-borne GW antennas, e.g., LISA, TianQin, and TaiJi, cover a wide range of GW frequencies {from $\sim 10^{-4}\;\rm Hz$ to $\sim 10^3\;\rm Hz$} and almost all types of compact objects in close orbits serve as the potential target sources for these GW detectors. The synergistic multi-band GW and EM observations would allow us to study fundamental physics from stars to cosmology. {The formation of stellar GW sources has been extensively explored in recent years, and progress on physical processes in binary interaction has been made as well. Furthermore, some studies have shown that the progress in binary evolution may significantly affect the properties of the stellar GW sources.}

We explore the morphological features and star formation activities of [OII] emitters in the COSMOS UltraDeep field at $z \sim 1.5$ using JWST NIRCam data from the COSMOS-Web survey and Subaru Hyper Suprime-Cam. We also report the discovery of large filamentary structures traced by [OII] emitters, surrounding an extremely overdense core with a galaxy number density $\sim11\times$ higher than the field average. These structures span over 50 cMpc, underscoring their large scale in the cosmic web at this epoch. After matching the stellar mass distributions, the core galaxies show a higher frequency of disturbances (50$\%$ $ \pm$ 9$\%$) than those in outskirts (41$\%$ $ \pm$ 9$\%$) and the field (21$\%$ $ \pm$ 5$\%$), indicative of more frequent mergers and interactions in the innermost $\lesssim1.5 $ arcmin region. Additionally, we observe that specific star formation rates are elevated in denser environments. A Kolmogorov-Smirnov (KS) test comparing the distribution of specific star formation rates of core and field galaxies yields a $\textit{p}$-value of 0.02, suggesting an enhancement of star-formation activity driven by the dense environment. Our findings underscore the environmental impact on galaxy evolution during a pivotal cosmic epoch and set the stage for further investigation with the increasing larger data from upcoming surveys.

We investigated the gas obscuration and host galaxy properties of active galactic nuclei (AGN) during the peak of cosmic accretion growth of supermassive black holes (SMBHs) at redshift 0.8-1.8 using X-ray detected AGN with mid-infrared and far-infrared detection. The sample was classified as type-1 and type-2 AGN using optical spectral and morphological classification while the host galaxy properties were estimated with multiwavelength SED fitting. For type-1 AGN, the black hole mass was determined from MgII emission lines while the black hole mass of type-2 AGN was inferred from the host galaxy's stellar mass. Based on the derived parameters, the distribution of the sample in the absorption hydrogen column density ($N_{\rm H}$) vs. Eddington ratio diagram is examined. Among the type-2 AGN, $28\pm5$\% are in the forbidden zone, where the obscuration by dust torus cannot be maintained due to radiation pressure on dusty material. The fraction is higher than that observed in the local universe from the BAT AGN Spectroscopic Survey (BASS) data release 2 ($11\pm3$\%). The higher fraction implies that the obscuration of the majority of AGN is consistent with the radiation pressure regulated unified model but with an increased incidence of interstellar matter (ISM) obscured AGN. We discuss the possibility of dust-free absorption in type-1 AGN and heavy ISM absorption in type-2 AGN. We also find no statistical difference in the star-formation activity between type-1 and type-2 AGN which may suggest that obscuration triggered by a gas-rich merging is not common among X-ray detected AGN in this epoch.

The origin of supermassive black holes (SMBHs) residing in the centers of most galaxies remains a mystery. Various growth models, such as accretion and hierarchical mergers, have been proposed to explain the existence and cosmological evolution of these SMBHs, but no direct observational evidence is available to test these models. The Event Horizon Telescope (EHT) offered direct imaging of nearby SMBHs, in particular, the one at the center of the Milky Way Galaxy named Sgr~A*. Measurements suggest that the Sgr~A* BH spins rapidly with significant spin axis misalignment relative to the angular momentum of the Galactic plane. Through investigating various SMBH growth models, here we show that the spin properties of Sgr~A* provides strong evidence of a past SMBH merger. Inspired by the merger between the Milky Way and Gaia-Enceladus, which has a 4:1 mass ratio as inferred from Gaia data, we have discovered that a 4:1 major merger of SMBH with a binary angular momentum inclination angle of 15-45 degrees with respect to the line of sight (LOS), can successfully replicate the measured spin properties of Sgr A*. This merger event in our galaxy provides observational support for the theory of hierarchical BH mergers in the formation and growth of SMBHs. The inferred merger rate, consistent with theoretical predictions, suggests a promising detection rate of SMBH mergers for space-borne gravitational wave detectors expected to operate in 2030s.

Exploring the concept of a massive photon has been an important area in astronomy and physics. If photons have mass, their propagation in non-vacuum space would be affected by both the non-zero mass $m_{\gamma}$ and the presence of a plasma medium. This would lead to a delay time proportional to $m_{\gamma}^2\nu^{-4}$, which deviates from the classical dispersion relation (proportional to $\nu^{-2}$). For the first time, we have derived the dispersion relation of a photon with a non-zero mass propagating in plasma. To reduce the impact of variations in the dispersion measure (DM), we employed the high-precision timing data to constrain the upper bound of the photon mass. Specifically, the DM/time of arrival (TOA) uncertainties derived from ultra-wide bandwidth (UWB) observations conducted by the Parkes Pulsar Timing Array (PPTA) are used. The de-dispersed pulses from fast radio bursts (FRBs) with minimal scattering effects are also used to constrain the upper bound of photon mass. The stringent limit on the photon mass is determined by uncertainties of the TOA of pulsars, with an optimum value of $9.52\times 10^{-46} \, \rm kg \,\,(5.34 \times 10^{-10}\, \rm eV/c^2$). In the future, it is essential to investigate the photon mass, as pulsar timing data are collected by PTA and UWB receivers, or FRBs with wide-band spectra are detected by UWB receivers.

We study the incidence and spatial distribution of galaxies that are currently undergoing gravitational merging (M) or that have signs of a post merger (PM) in six galaxy clusters (A754, A2399, A2670, A3558, A3562, and A3716) within the redshift range, 0.05$\lesssim$$z$$\lesssim$0.08. To this aim, we obtained Dark Energy Camera (DECam) mosaics in $u^{\prime}$, $g^{\prime}$, and $r^{\prime}$-bands covering up to $3\times R_{200}$ of the clusters, reaching 28 mag/arcsec$^2$ surface brightness limits. We visually inspect $u^{\prime}$$g^{\prime}$$r^{\prime}$ color-composite images of volume-limited ($M_r < -20$) cluster-member galaxies to identify whether galaxies are of M or PM types. We find 4% M-type and 7% PM-type galaxies in the galaxy clusters studied. By adding spectroscopic data and studying the projected phase space diagram (PPSD) of the projected clustocentric radius and the line-of-sight velocity, we find that PM-type galaxies are more virialized than M-type galaxies, having 1--5% point higher fraction within the escape-velocity region, while the fraction of M-type was $\sim$10% point higher than PM-type in the intermediate environment. Similarly, on a substructure analysis, M types were found in the outskirt groups, while PM types populated groups in ubiquitous regions of the PPSD. Adopting literature-derived dynamical state indicator values, we observed a higher abundance of M types in dynamically relaxed clusters. This finding suggests that galaxies displaying post-merging features within clusters likely merged in low-velocity environments, including cluster outskirts and dynamically relaxed clusters.

We aim at searching for secular evolution of the orbital period in the short-period binary system WR 127 (WN3b+O9.5V, $P = 9.555^d$). We performed new low-resolution spectroscopic observations of WR 127 on 2.5-m CMO SAI telescope to construct the radial velocity curves of the components suggesting the component masses $M_\mathrm{WR}\sin^3(i) = 11.8\pm1.4$ $M_{\odot}$, $M_\mathrm{O}\sin^3(i)=17.2\pm1.4$ $M_{\odot}$. The comparison with archival radial velocity curves enabled us to calculate the $(O-C)$ plot with accuracy sufficient to search for the orbital period change in WR 127. We report on the reliable detection of a secular increase in the orbital period of WR 127 at a rate of $\dot{P} = 0.83\pm0.14~\mbox{s~yr}^{-1}$ corresponding to the dynamical mass-loss rate from the WR star $\dot{M}_\mathrm{WR} = (2.6\pm0.5)\times 10^{-5}$ $M_{\odot}$ yr$^{-1}$. The mass-loss rate from WR stars in three WR+OB binaries (WR 127, CX Cep and V444 Cyg) as inferred from spectroscopic and photometric measurements suggests a preliminary empirical correlation between the WR star mass and the dynamical mass-loss rate $\dot M_\mathrm{WR}\sim M_\mathrm{WR}^{1.8}$. This relation is important for the understanding of the evolution of massive close binaries with WR stars -- precursors of gravitational-wave binary merging events with neutron stars and black holes.

Velocity fields of molecular clouds (MCs) can provide crucial information on the merger and split between clouds, as well as their internal kinematics and maintenance, energy injection and redistribution, even star formation within clouds. Using the CO spectral lines data from the Milky Way Imaging Scroll Painting (MWISP) survey, we measure the relative velocities along the line of sight ($\Delta$V$_{\rm LOS}$) between $^{13}$CO structures within $^{12}$CO MCs. Emphasizing MCs with double and triple $^{13}$CO structures, we find that approximately 70$\%$ of $\Delta$V$_{\rm LOS}$ values are less than $\sim$ 1 km s$^{-1}$, and roughly 10$\%$ of values exceed 2 km s$^{-1}$, with a maximum of $\sim$ 5 km s$^{-1}$. Additionally, we compare $\Delta$V$_{\rm LOS}$ with the internal velocity dispersion of $^{13}$CO structures ($\sigma_{\rm ^{13}CO,in}$) and find that about 40$\%$ of samples in either double or triple regime display distinct velocity discontinuities, i.e. the relative velocities between $^{13}$CO structures are larger than the internal linewidths of $^{13}$CO structures. Among these 40$\%$ samples in the triple regime, 33$\%$ exhibit signatures of combinations through the two-body motion, whereas the remaining 7$\%$ show features of configurations through the multiple-body motion. The $\Delta$V$_{\rm LOS}$ distributions for MCs with double and triple $^{13}$CO structures are similar, as well as their $\Delta$V$_{\rm LOS}$/$\sigma_{\rm ^{13}CO,in}$ distributions. This suggests that relative motions of $^{13}$CO structures within MCs are random and independent of cloud complexities and scales.

We present first-order models for tilt-to-length (TTL) coupling in LISA, both for the individual interferometers as well as in the time-delay interferometry (TDI) Michelson observables. These models include the noise contributions from angular and lateral jitter coupling of the six test masses, six movable optical subassemblies (MOSAs), and three spacecraft. We briefly discuss which terms are considered to be dominant and reduce the TTL model for the second-generation TDI Michelson X observable to these primary noise contributions to estimate the resulting noise level. We show that the expected TTL noise will initially violate the entire mission displacement noise budget, resulting in the known necessity to fit and subtract TTL noise in data post-processing. By comparing the noise levels for different assumptions prior to subtraction, we show why noise mitigation by realignment prior to subtraction is favorable. We then discuss that the TTL coupling in the individual interferometers will have noise contributions that will not be present in the TDI observables. Models for TTL coupling noise in TDI and in the individual interferometers are therefore different, and commonly made assumptions are valid as such only for TDI but not for the individual interferometers. Finally, we analyze what implications can be drawn from the presented models for the subsequent fit-and-subtraction in post-processing. We show that noise contributions from the test mass and inter-satellite interferometers are indistinguishable, such that only the combined coefficients can be fit and used for subtraction. However, a distinction is considered not necessary. Additionally, we show a correlation between coefficients for transmitter and receiver jitter couplings in each individual TDI Michelson observable. This full correlation can be resolved by using all three Michelson observables for fitting the TTL coefficients.

We present a new determination of the star-forming main sequence (MS), obtained through stacking 100k K-band-selected galaxies in the far-IR Herschel and James Clerk Maxwell Telescope (JCMT) imaging. By fitting the dust emission curve to the stacked far-IR photometry, we derive the IR luminosities (LIR) and, hence, star formation rates (SFR) out to z~7. The functional form of the MS is found, with the linear SFR-M* relation that flattens at high stellar masses and the normalization that increases exponentially with redshift. We derive the corresponding redshift evolution of the specific star formation rate (sSFR) and compare our findings with the recent literature. We find our MS to be exhibiting slightly lower normalization at z<=2 and to flatten at larger stellar masses at high redshifts. By deriving the relationship between the peak dust temperature (Td) and redshift, where Td increases linearly from ~20K at z=0.5 to ~50 K at z=6, we conclude that the apparent inconsistencies in the shapes of the MS are most likely caused by the different dust temperatures assumed when deriving SFRs in the absence of far-IR data. Finally, we investigate the derived shape of the star-forming MS by simulating the time evolution of the observed galaxy stellar mass function (GSMF). While the simulated GSMF is in good agreement with the observed one, some inconsistencies persist. In particular, we find the simulated GSMF to be somewhat overpredicting the number density of low-mass galaxies at z>2.

The geometry of a coronal hole (CH) affects the density profile of the reflected part of an incoming global coronal wave (CW). In this study, we perform for the first time magnetohydrodynamic (MHD) simulations of fast-mode MHD waves interacting with CHs of different geometries, such as circular, elliptic, convex, and concave shapes. We analyse the influence these geometries have on the density profiles of the reflected waves and we generate the corresponding simulation-based time-distance plots. Within these time-distance plots we determine regions that exhibit specific density features, such as large reflected density amplitudes. In a further step, these interaction features can be compared to actual observed CW-CH interaction events which makes it possible to explain interaction parameters from the observed interaction events, such as the density structure of the reflected wave, which are usually difficult to comprehensively understand by only analysing the measurements. Moreover, we show that the interaction between a concave shaped CH and CWs, whose density profile include an enhanced as well as a depleted wave, can lead to reflected density amplitudes that are more than two times larger than the incoming ones. Another effect of the interplay between the constructive and destructive interference of the reflected wave parts is a strongly depleted region in the middle of the CW-CH interaction process. In addition, we show how important the choice of the path is that is used to generate the time-distance plots and how this choice affects the interpretation of the CW-CH interaction results.

We discuss the possibility to observe neutrinos emitted from a failed core collapse Supernova in the various experiments at Laboratori Nazionali del Gran Sasso. We show that the veto regions of dark matter and neutrinoless double beta decay experiments can be used as a network of small detectors to measure Supernova neutrinos. In addition we show that this network can measure very precisely the moment of black hole formation, which can be then used in the nearby VIRGO detector and future Einstein Telescope to look for the gravitational wave counterpart to the neutrino signal.

We propagate relativistic test particles in the field of a steady 3D MHD simulations of the solar wind. We use the MPI-AMRVAC code for the wind simulations and integrate the relativistic guiding center equations using a new third-order accurate time integration scheme to solve the particle trajectories. Diffusion in velocity space, given a particle-turbulence mean free path $\lambda_\parallel$ along the magnetic field, is also included. Preliminary results for $81\:{\rm keV}$ electrons injected at 0.139 AU heliocentric distance and mean free path $\lambda_\parallel =0.5\:{\rm AU}$ are in a good qualitative agreement with measurements at 1 AU.

The composition of the solar corona differs from that of the photosphere, with the plasma thought to fractionate in the solar chromosphere according to the First Ionisation Potential (FIP) of the different elements. This produces a FIP bias, wherein elements with a low FIP are preferentially enhanced in the corona compared to their photospheric abundance, but direct observations of this process remain elusive. Here we use a series of spectroscopic observations of Active Region AR 12759 as it transited the solar disc over a period of 6 days from 2-7 April 2020 taken using the Hinode Extreme ultraviolet Imaging Spectrometer (EIS) and Interface Region Imaging Spectrograph (IRIS) instruments to look for signatures of plasma fractionation in the solar chromosphere. Using the Si X/S X and Ca XIV/Ar XIV diagnostics, we find distinct differences between the FIP bias of the leading and following polarities of the active region. The widths of the IRIS Si IV lines exhibited clear differences between the leading and following polarity regions, indicating increased unresolved wave activity in the following polarity region compared to the leading polarity region, with the chromospheric velocities derived using the Mg II lines exhibiting comparable, albeit much weaker, behaviour. These results are consistent with plasma fractionation via resonant/non-resonant waves at different locations in the solar chromosphere following the ponderomotive force model, and indicate that IRIS could be used to further study this fundamental physical process.

Precise measurements of energy spectra of different cosmic ray species were obtained in recent years, by particularly the AMS-02 experiment on the International Space Station. It has been shown that apparent differences exist in different groups of the primary cosmic rays. However, it is not straightforward to conclude that the source spectra of different particle groups are different since they will experience different propagation processes (e.g., energy losses and fragmentations) either. In this work, we study the injection spectra of different nuclear species using the measurements from Voyager-1 outside the solar system, and ACR-CRIS and AMS-02 on top of atmosphere, in a physical framework of cosmic ray transportation. Two types of injection spectra are assumed, the broken power-law and the non-parametric spline interpolation form. The non-parametric form fits the data better than the broken power-law form, implying that potential structures beyond the constrained spectral shape of broken power-law may exist. For different nuclei the injection spectra are overall similar in shape but do show some differences among each other. For the non-parametric spectral form, the helium injection spectrum is the softest at low energies and the hardest at high energies. For both spectral shapes, the low-energy injection spectrum of neon is the hardest among all these species, and the carbon and oxygen spectra have more prominent bumps in 1-10 GV in the R2dN/dR presentation. Such differences suggest the existence of differences in the sources or acceleration processes of various nuclei of cosmic rays.

We super-resolve the seeing-limited Subaru Hyper Suprime-Cam (HSC) images for 32,187 galaxies at z~2-5 in three techniques, namely, the classical Richardson-Lucy (RL) point spread function (PSF) deconvolution, sparse modeling, and generative adversarial networks to investigate the environmental dependence of galaxy mergers. These three techniques generate overall similar high spatial resolution images but with some slight differences in galaxy structures, for example, more residual noises are seen in the classical RL PSF deconvolution. To alleviate disadvantages of each technique, we create combined images by averaging over the three types of super-resolution images, which result in galaxy sub-structures resembling those seen in the Hubble Space Telescope images. Using the combined super-resolution images, we measure the relative galaxy major merger fraction corrected for the chance projection effect, f_merg, for galaxies in the ~300 deg^2-area data of the HSC Strategic Survey Program and the CFHT Large Area U-band Survey. Our f_merg measurements at z~3 validate previous findings showing that f_merg is higher in regions with a higher galaxy overdensity delta at z~2-3. Thanks to the large galaxy sample, we identify a nearly linear increase in f_merg with increasing delta at z~4-5, providing the highest-z observational evidence that galaxy mergers are related to delta. In addition to our f_merg measurements, we find that the galaxy merger fractions in the literature also broadly align with the linear f_merg-delta relation across a wide redshift range of z~2-5. This alignment suggests that the linear f_merg-delta relation can serve as a valuable tool for quantitatively estimating the contributions of galaxy mergers to various environmental dependences. This super-resolution analysis can be readily applied to datasets from wide field-of-view space telescopes such as Euclid and Roman.

We present a Wideband Tapered Gridded Estimator (TGE), which incorporates baseline migration and variation of the primary beam pattern for neutral hydrogen (${\rm H\hspace{0.5mm}}{\scriptsize {\rm I}}$) 21-cm intensity mapping (IM) with large frequency bandwidth radio-interferometric observations. Here we have analysed $394-494 \, {\rm MHz}$ $(z = 1.9 - 2.6)$ uGMRT data to estimate the Multi-frequency Angular Power Spectrum (MAPS) $C_\ell(\Delta\nu)$ from which we have removed the foregrounds using the polynomial fitting (PF) and Gaussian Process Regression (GPR) methods developed in our earlier work. Using the residual $C_\ell(\Delta\nu)$ to estimate the mean squared 21-cm brightness temperature fluctuation $\Delta^2(k)$, we find that this is consistent with $0 \pm 2 \sigma$ in several $k$ bins. The resulting $2\sigma$ upper limit $\Delta^2(k) < (4.68)^2 \, \rm{mK^2}$ at $k=0.219\,\rm{Mpc^{-1}}$ is nearly $15$ times tighter than earlier limits obtained from a smaller bandwidth ($24.4 \, {\rm MHz}$) of the same data. The $2\sigma$ upper limit $[\Omega_{{\rm H\hspace{0.5mm}}{\scriptsize {\rm I}}} b_{{\rm H\hspace{0.5mm}}{\scriptsize {\rm I}}}] < 1.01 \times 10^{-2}$ is within an order of magnitude of the value expected from independent estimates of the ${\rm H\hspace{0.5mm}}{\scriptsize {\rm I}}$ mass density $\Omega_{{\rm H\hspace{0.5mm}}{\scriptsize {\rm I}}}$ and the ${\rm H\hspace{0.5mm}}{\scriptsize {\rm I}}$ bias $b_{{\rm H\hspace{0.5mm}}{\scriptsize {\rm I}}}$. The techniques used here can be applied to other telescopes and frequencies, including $\sim 150 \, {\rm MHz}$ Epoch of Reionization observations.

We present point-source photometry from the Spitzer Space Telescope's final survey of the Small Magellanic Cloud (SMC). We mapped 30 square degrees in two epochs in 2017, with the second extending to early 2018 at 3.6 and 4.5 microns using the Infrared Array Camera. This survey duplicates the footprint from the SAGE-SMC program in 2008. Together, these surveys cover a nearly 10 yr temporal baseline in the SMC. We performed aperture photometry on the mosaicked maps produced from the new data. We did not use any prior catalogs as inputs for the extractor in order to be sensitive to any moving objects (e.g., foreground brown dwarfs) and other transient phenomena (e.g., cataclysmic variables or FU Ori-type eruptions). We produced a point-source catalog with high-confidence sources for each epoch as well as combined-epoch catalog. For each epoch and the combined-epoch data, we also produced a more complete archive with lower-confidence sources. All of these data products will be available to the community at the Infrared Science Archive.

Among multi-messenger observations of the next galactic core-collapse supernova, Super-Kamiokande (SK) plays a critical role in detecting the emitted supernova neutrinos, determining the direction to the supernova (SN), and notifying the astronomical community of these observations in advance of the optical signal. On 2022, SK has increased the gadolinium dissolved in its water target (SK-Gd) and has achieved a Gd concentration of 0.033%, resulting in enhanced neutron detection capability, which in turn enables more accurate determination of the supernova direction. Accordingly, SK-Gd's real-time supernova monitoring system (Abe te al. 2016b) has been upgraded. SK_SN Notice, a warning system that works together with this monitoring system, was released on December 13, 2021, and is available through GCN Notices (Barthelmy et al. 2000). When the monitoring system detects an SN-like burst of events, SK_SN Notice will automatically distribute an alarm with the reconstructed direction to the supernova candidate within a few minutes. In this paper, we present a systematic study of SK-Gd's response to a simulated galactic SN. Assuming a supernova situated at 10 kpc, neutrino fluxes from six supernova models are used to characterize SK-Gd's pointing accuracy using the same tools as the online monitoring system. The pointing accuracy is found to vary from 3-7$^\circ$ depending on the models. However, if the supernova is closer than 10 kpc, SK_SN Notice can issue an alarm with three-degree accuracy, which will benefit follow-up observations by optical telescopes with large fields of view.

Computing the tidal deformations of Mars, we explored various Mars internal structures by examining profiles that include or exclude a basal molten layer within the mantle and a solid inner core. By assessing their compatibility with a diverse set of geophysical observations we show that despite the very short periods of excitation, tidal deformation is very efficient to constraint the Mars interior. We calculated densities and thicknesses for Martian lithosphere, mantle, {core-mantle boundary} layer and core and {found them} coherent with preexisting results from other methods. We also estimated new viscosities for these layers. We demonstrated that the geodetic record associated with thermal constraint is very sensitive to the {presence} of a basal molten layer in deep Martian mantle, and less sensitive to the solid core. Our results also indicate that the existence of the basal molten layer necessarily comes together with an inversion of viscosity between the lithosphere and the mantle. In this case, we could attribute this reverse viscosity contrast to the poor hydration state of Martian mantle and we underlined that this result prevents a strict Earth-like plate tectonics on Mars. The existence of the basal molten layer is also associated with a non-inversion of viscosity between the core-mantle boundary layer and the liquid core. Finally, in our results, a basal molten layer is incompatible with the existence of a solid inner core. Efforts to detect basal molten layer are then of prime importance to decipher the Martian interior. {Inversely, viscosity profiles appear to be very good tools for probing the existence of such molten layer at the base of the Mars mantle.

Tidal disruptions of stars on the equatorial plane orbiting Kerr black holes have been widely studied. However thus far, there have been fewer studies of stars in inclined precessing orbits around a Kerr black hole. In this paper, by using tensor virial equations, we show the presence of possible resonances in these systems for typical physical parameters of black hole-neutron star binaries in close orbits or of a white dwarf/an ordinary star orbiting a supermassive black hole. This suggests the presence of a new instability before the tidal disruption limit is encountered in such systems.

The variations in metallicity and spatial patterns within star-forming regions of galaxies result from diverse physical processes unfolding throughout their evolutionary history, with a particular emphasis in recent events. Analysing MaNGA and \textsc{eagle} galaxies, we discovered an additional dependence of the mass-metallicity relation (MZR) on metallicity gradients ($\nabla_{{\rm (O/H)}}$). Two regimes emerged for low and high stellar mass galaxies, distinctly separated at approximately ${\rm M_{\star}} >10^{9.75}$. Low-mass galaxies with strong positive $\nabla_{{\rm (O/H)}}$ appear less enriched than the MZR median, while those with strong negative gradients are consistently more enriched in both simulated and observed samples. Interestingly, low-mass galaxies with strong negative $\nabla_{{\rm (O/H)}}$ exhibit high star-forming activity, regardless of stellar surface density or $\nabla_{{\rm (O/H)}}$. In contrast, a discrepancy arises for massive galaxies between MaNGA and \textsc{eagle} datasets. The latter exhibit a notable anticorrelation between specific star formation rate and stellar surface density, independent of $\nabla_{{\rm (O/H)}}$, while MaNGA galaxies show this trend mainly for strong positive $\nabla_{{\rm (O/H)}}$. Further investigation indicates that galaxies with strong negative gradients tend to host smaller central black holes in observed datasets, a trend not replicated in simulations. These findings suggest disparities in metallicity recycling and mixing history between observations and simulations, particularly in massive galaxies with varying metallicity gradients. These distinctions could contribute to a more comprehensive understanding of the underlying physics.

The computational fluid dynamics on a sphere is relevant to global simulations of geophysical fluid dynamics. Using the conventional spherical-polar (or lat-lon) grid results in a singularity at the poles, with orders of magnitude smaller cell sizes at the poles in comparison to the equator. To address this problem, we developed a general circulation model (dynamic core) with a gnomonic equiangular cubed-sphere configuration. This model is developed based on the Simulating Nonhydrostatic Atmospheres on Planets (SNAP) model, using a finite volume numerical scheme with a Riemann-solver-based dynamic core and the vertical implicit correction (VIC) scheme. This change of the horizontal configuration gives a 20-time acceleration of global simulations compared to the lat-lon grid with a similar number of cells at medium resolution. We presented standard tests ranging from 2D shallow-water models to 3D general circulation tests, including earth-like planets and shallow hot Jupiters, to validate the accuracy of the model. The method described in this article is generic to transform any existing finite-volume hydrodynamic model in the Cartesian geometry to the spherical geometry.

Utilizing the well-established Radial Tully-Fisher (RTF) relation observed in a `large' (843) sample of local galaxies, we report the maximum allowed variance in the Hubble parameter, $H_0$. We estimate the total intrinsic scatter in the magnitude of the RTF relation(s) implementing a cosmological model-independent cosmographic expansion. We find that the maximum allowed local variation in our baseline analysis, using 4 RTF relations in the galaxy sample is $\Delta H_0/H_0 \lesssim 3 \%$ at a $95\%$ C.L. significance. Which is implied form a constraint of $\Delta H_0/H_0 = 0.54^{+1.32}_{-1.37} \%$ estimated at $D_{\rm{L}}\sim 10\, [\rm{Mpc}]$. Using only one `best-constrained' radial bin we report a conservative $95\%$ C.L. limit of $\Delta H_0/H_0 \lesssim 4 \%$. Through our estimate of maximum variation, we propose a novel method to validate several late-time/local modifications put forth to alleviate the $H_0$ tension. We find that within the range of the current galaxy sample redshift distribution $10 \, [\rm{Mpc}] \le D_{\rm{L}} \le 140\, [\rm{Mpc}]$, it is highly unlikely to obtain a variation of $\Delta H_0/H_0 \sim 9\%$, necessary to alleviate the $H_0$-tension. However, we also elaborate on the possible alternative inferences when the innermost radial bin is included in the analysis. Alongside the primary analysis of fitting the individual RTF relations independently, we propose and perform a joint analysis of the RTF relations useful to create a pseudo-standardizable sample of galaxies. We also test for the spatial variation of $H_0$, finding that the current samples' galaxies distributed only in the southern hemisphere support the null hypothesis of isotropy, within the allowed noise levels.

We present the multi-epoch analysis of 13 variable, nearby (z<0.1), Compton-thin (22<logN_H<24) active galactic nuclei (AGN) selected from the 105-month BAT catalog. Analyzing all available archival soft and hard X-ray observations, we investigate the line-of-sight hydrogen column density (N_H) variability on timescales ranging from a few days to approximately 20 years. Each source is analyzed by simultaneously modeling the data with three physical torus models, providing tight constraints on torus properties, including the covering factor, the cloud dispersion, and the torus average hydrogen column density (N_H,av). For each epoch, we measure the N_H and categorize the source as `N_H Variable', `Non-variable in N_H', or `Undetermined' based on the degree of variability. Our final sample includes 27 variable, Compton-thin AGN after implementing another 14 AGN analyzed in our previous work. We find that all sources require either flux or N_H variability. We classify 37% of them as `N_H Variable', 44% as `Non-variable in N_H', and 19% as `Undetermined'. Noticeably, there is no discernible difference between geometrical and intrinsic properties among the three variability classes, suggesting no intrinsic differences between the N_H-variable and non-variable sources. We measure the median variation in N_H between any observation pair of the same source to be 25% with respect to the lowest N_H measure in the pair. Furthermore, 48% of the analyzed sources require the inclusion of a Compton-thick reflector in the spectral fitting. Among these, the 30% exhibits recorded 22 GHz water megamaser emission, suggesting a potential shared nature between the two structures.

Using a 3D General Circulation Model, the Unified Model, we present results from simulations of a tidally-locked TRAPPIST-1e with varying carbon dioxide CO2 and methane CH4 gas concentrations, and their corresponding prescribed spherical haze profiles. Our results show that the presence of CO2 leads to a warmer atmosphere globally due to its greenhouse effect, with the increase of surface temperature on the dayside surface reaching up to ~14.1 K, and on the nightside up to ~21.2 K. Increasing presence of CH4 first elevates the surface temperature on the dayside, followed by a decrease due to the balance of tropospheric warming and stratospheric cooling. A thin layer of haze, formed when the partial pressures of CH4 to CO2 (pCH4/pCO2) = 0.1, leads to a dayside warming of ~4.9K due to a change in the water vapour H2O distribution. The presence of a haze layer that formed beyond the ratio of 0.1 leads to dayside cooling. The haze reaches an optical threshold thickness when pCH4/pCO2 ~0.4 beyond which the dayside mean surface temperature does not vary much. The planet is more favourable to maintaining liquid water on the surface (mean surface temperature above 273.15 K) when pCO2 is high, pCH4 is low and the haze layer is thin. The effect of CO2, CH4 and haze on the dayside is similar to that for a rapidly-rotating planet. On the contrary, their effect on the nightside depends on the wind structure and the wind speed in the simulation.

We present a catalog that we named MODEST containing depth-dependent information on the atmospheric conditions inside sunspot groups of all types. The catalog is currently composed of 942 observations of 117 individual active regions with sunspots that cover all types of features observed in the solar photosphere. We use the SPINOR-2D code to perform spatially coupled inversions of the Stokes profiles observed by Hinode/SOT-SP at high spatial resolution. SPINOR-2D accounts for the unavoidable degradation of the spatial information due to the point spread function of the telescope. The sunspot sample focuses on complex sunspot groups, but simple sunspots are also part of the catalog for completeness. Sunspots were observed from 2006 to 2019, covering parts of solar cycles 23 and 24. The catalog is a living resource, as with time, more sunspot groups will be included.

Scalar-induced Gravitational Waves (SIGWs) represent a particular class of primordial signals which are sourced at second-order in perturbation theory whenever a scalar fluctuation of the metric is present. They form a guaranteed Stochastic Gravitational Wave Background (SGWB) that, depending on the amplification of primordial scalar fluctuations, can be detected by GW detectors. The amplitude and the frequency shape of the scalar-induced SGWB can be influenced by the statistical properties of the scalar density perturbations. In this work we study the intuitive physics behind SIGWs and we analyze the imprints of local non-Gaussianity of the primordial curvature perturbation on the GW spectrum. We consider all the relevant non-Gaussian contributions up to fifth-order in the scalar seeds without any hierarchy, and we derive the related GW energy density $\Omega_{\rm GW}(f)$. We perform a Fisher matrix analysis to understand to which accuracy non-Gaussianity can be constrained with the LISA detector, which will be sensitive in the milli-Hertz frequency band. We find that LISA, neglecting the impact of astrophysical foregrounds, will be able to measure the amplitude, the width and the peak of the spectrum with an accuracy up to $\mathcal{O}(10^{-4})$, while non-Gaussianity can be measured up to $\mathcal{O}(10^{-3})$. Finally, we discuss the implications of our non-Gaussianity expansion on the fraction of Primordial Black Holes.

X-ray spectroscopy of heavily obscured Active Galactic Nuclei (AGN) offers a unique opportunity to study the circum-nuclear environment of accreting supermassive black holes (SMBHs). However, individual models describing the obscurer have unique parameter spaces that give distinct parameter posterior distributions when fit to the same data. To assess the impact of model-specific parameter dependencies, we present a case study of the nearby heavily obscured low-luminosity AGN NGC 3982, which has a variety of column density estimations reported in the literature. We fit the same broadband XMM-Newton + NuSTAR spectra of the source with five unique obscuration models and generate posterior parameter distributions for each. By using global parameter exploration, we traverse the full prior-defined parameter space to accurately reproduce complex posterior shapes and inter-parameter degeneracies. The unique model posteriors for the line-of-sight column density are broadly consistent, predicting Compton-thick $N_{\rm H}$ $>1.5\times10^{24}\rm cm^{-2}$ at the 3$\sigma$ confidence level. The posterior median intrinsic X-ray luminosity in the 2-10 keV band however was found to differ substantially, with values in the range log $L_{ 2-10\,{\rm keV}}$ergs$^{-1}$ = 40.9-42.1 for the individual models. We additionally show that the posterior distributions for each model occupy unique regions of their respective multi-dimensional parameters spaces, and how such differences can propagate into the inferred properties of the central engine. We conclude by showcasing the improvement in parameter inference attainable with the High Energy X-ray Probe (HEX-P) with a uniquely broad simultaneous and high-sensitivity bandpass of 0.2-80 keV.

Tidal torques can alter the spins of tidally interacting stars and planets, usually over shorter timescales than the tidal damping of orbital separations or eccentricities. Simple tidal models predict that in eccentric binary or planetary systems, rotation periods will evolve toward a "pseudosynchronous" ratio with the orbital period. However, this prediction does not account for "inertial" waves that are present in stars or gaseous planets with (i) convective envelopes, and (ii) even very slow rotation. We demonstrate that tidal driving of inertial oscillations in eccentric systems generically produces a network of stable "synchronization traps" at ratios of orbital to rotation period that are simple to predict, but can deviate significantly from pseudosynchronization. The mechanism underlying spin synchronization trapping is similar to tidal resonance locking, involving a balance between torques that is maintained automatically by the scaling of inertial mode frequencies with the rotation rate. In contrast with many resonance locking scenarios, however, the torque balance required for synchronization trapping need not drive mode amplitudes to nonlinearity. Synchronization traps may provide an explanation for low-mass stars and hot Jupiters with observed rotation rates that deviate from pseudosynchronous or synchronous expectations.

We present a James Webb Space Telescope (JWST) imaging survey of I Zw 18, the archetypal extremely metal-poor, star-forming, blue compact dwarf galaxy. With an oxygen abundance of only $\sim$3% $Z_{\odot}$, it is among the lowest-metallicity systems known in the local universe, and is, therefore, an excellent accessible analog for the galactic building blocks which existed at early epochs of ionization and star formation. These JWST data provide a comprehensive infrared (IR) view of I Zw 18 with eight filters utilizing both NIRCam (F115W, F200W, F356W, and F444W) and MIRI (F770W, F1000W, F1500W, and F1800W) photometry, which we have used to identify key stellar populations that are bright in the near- and mid-IR. These data allow for a better understanding of the origins of dust and dust-production mechanisms in metal-poor environments by characterizing the population of massive, evolved stars in the red supergiant (RSG) and asymptotic giant branch (AGB) phases. In addition, it enables the identification of the brightest dust-enshrouded young stellar objects (YSOs), which provide insight into the formation of massive stars at extremely low metallicities typical of the very early universe. This paper provides an overview of the observational strategy and data processing, and presents first science results, including identifications of dusty AGB star, RSG, and bright YSO candidates. These first results assess the scientific quality of JWST data and provide a guide for obtaining and interpreting future observations of the dusty and evolved stars inhabiting compact dwarf star-forming galaxies in the local universe.

Anti-quark nuggets (AQNs) have been suggested to solve the dark matter (DM) and the missing antimatter problem in the universe and have been proposed as an explanation of various observations. Their size is in the {\mu}m range and their density is about equal to the nuclear density with an expected flux of about $0.4 / km^2 / year$. For the typical velocity of DM constituents ($\sim$250 km/s), the solar system bodies act as highly performing gravitational lenses. Here we assume that DM streams or clusters are impinging, e.g., on the Earth, as it was worked out for DM axions and Weakly Interacting Massive Particles (WIMPs). Interestingly, in the LHC beam, unforeseen beam losses are triggered by so-called Unidentified Falling Objects (UFOs), which are believed to be constituted of dust particles with a size in the {\mu}m range and a density of several orders of magnitude lower than AQNs. Prezeau suggested that streaming DM constituents incident on the Earth should result in jet-like structures ("hairs") exiting the Earth, or a kind of caustics. Such ideas open novel directions in the search for DM. This work suggests a new analysis of the UFO results at the Large Hadron Collider (LHC), assuming that they are eventually, at least partly, due to AQNs. Firstly, a reanalysis of the existing data from the 4000 beam monitors since the beginning of the LHC is proposed, arguing that dust and AQNs should behave differently. The feasibility of this idea has been discussed with CERN accelerator people and potential collaborators.

Relativistic bubble walls from cosmological phase transitions (PT) necessarily accumulate expanding shells of particles. We systematically characterize shell properties, and identify and calculate the processes that prevent them from free streaming: phase-space saturation effects, out-of-equilibrium $2\to2$ and $3\to2$ shell-shell and shell-bath interactions, and shell interactions with bubble walls. We find that shells do not free stream in scenarios widely studied in the literature, where standard predictions will need to be reevaluated, including those of bubble wall velocities, gravitational waves (GW) and particle production. Our results support the use of bulk-flow GW predictions in all regions where shells free stream, irrespectively of whether or not the latent heat is mostly converted in the scalar field gradient.

We explore the possibility of %dynamically producing the observed matter-antimatter asymmetry of the Universe uniquely from the evaporation of primordial black holes (PBH) that are formed in an inflaton-dominated background. Considering the inflaton $(\phi)$ to oscillate in a monomial potential $V(\phi)\propto\phi^n$, we show, it is possible to obtain the desired baryon asymmetry via vanilla leptogenesis from evaporating PBHs of initial mass $\lesssim 10$ g. We find that the allowed parameter space is heavily dependent on the shape of the inflaton potential during reheating (determined by the exponent of the potential $n$), the energy density of PBHs (determined by $\beta$), and the nature of the coupling between the inflaton and the Standard Model (SM). To complete the minimal gravitational framework, we also include in our analysis the gravitational leptogenesis set-up through inflaton scattering via exchange of graviton, which opens up an even larger window for PBH mass, depending on the background equation of state. We finally illustrate that such gravitational leptogenesis scenarios can be tested with upcoming gravitational wave (GW) detectors, courtesy of the blue-tilted primordial GW with inflationary origin, thus paving a way to probe a PBH-induced reheating together with leptogenesis.

Binary neutron star mergers produce massive, hot, rapidly differentially rotating neutron star remnants; electromagnetic and gravitational wave signals associated with the subsequent evolution depend on the stability of these remnants. Stability of relativistic stars has previously been studied for uniform rotation and for a class of differential rotation with monotonic angular velocity profiles. Stability of those equilibria to axisymmetric perturbations was found to respect a turning point criterion: along a constant angular momentum sequence, the onset of unstable stars is found at maximum density less than but close to the density of maximum mass. In this paper, we test this turning point criterion for non-monotonic angular velocity profiles and non-isentropic entropy profiles, both chosen to more realistically model post-merger equilibria. Stability is assessed by evolving perturbed equilibria in 2D using the Spectral Einstein Code. We present tests of the code's new capability for axisymmetric metric evolution. We confirm the turning point theorem and determine the region of our rotation law parameter space that provides highest maximum mass for a given angular momentum.

The dark matter halo has non-negligible effects on the gravitational lensing of supermassive black hole in the galaxy center. Our work gives a study on the time-delay of light in gravitational lensings of black holes enclosed by dark matter halos. To provide a precise description on the distribution of dark matter in galaxies, we choose several famous phenomenological dark matter halo models in astrophysics, including the NFW, Beta, Burkert and Moore models, to carry out the present study. Through numerically calculating the time-delay of light in gravitational lensing, a comparative analysis of the dark matter effects within different halo models has been performed. Assuming typical length scales associated with the galactic gravitational lensing, numerical results indicate that the NFW, Beta, Moore dark matter halos can significantly enhance the time delay of light. Conversely, the dark matter in Burkert model exhibits a negative contribution on the time delay, resulting in smaller time delay compared with those for pure black hole cases. Keywords: Black Hole; Gravitational Lensing; Time Delay; Dark Matter Halo

The violations of parity and Lorentz symmetries in gravity can change the propagating properties of gravitational waves (GWs) in the cosmological background, which can arise from a large number of parity- and Lorentz-violating theories. In this paper, through a systematic parametrization for characterizing possible derivations from the standard GW propagation in general relativity, we study both the parity- and Lorentz-violating effects on the power spectra and the polarization of the primordial gravitational waves (PGWs) during the slow-roll inflation. To this end, we calculate explicitly the power spectrum and the corresponding circular polarization of the PGWs analytically by using the uniform asymptotic approximation. It is shown that the new contributions to power spectra contain two parts, one from the parity-violating terms and the other from the Lorentz-violating terms. While the Lorentz-violating terms can only affect the overall amplitudes of PGWs, the parity-violating terms induce nonzero circular polarization of PGWs, i.e., the left-hand and right-hand polarization modes of GWs have different amplitudes.

The Barbero-Immirzi parameter $\gamma$ appears as a coupling constant in the spinfoam dynamics of loop quantum gravity and can be understood as a measure of gravitational parity violation via a duality rotation. We investigate an effective field theory for gravity and a scalar field, with dynamics given by a $\gamma$-dual action obtained via a duality rotation of a parity-non-violating one. The resulting relation between the coupling constants of parity-even and parity-odd higher-curvature terms is determined by $\gamma$, opening the possibility of its measurement in the semiclassical regime. For a choice of $\gamma$-dual effective action, we study cosmic inflation and show that the observation of a primordial tensor polarization, together with the tensor tilt and the tensor-to-scalar ratio, provides a measurement of the Barbero-Immirzi parameter and, therefore, of the scale of discreteness of the quantum geometry of space.

We propose new slow-roll approximations for inflationary models with the Gauss-Bonnet term. We find more accurate expressions of the standard slow-roll parameters as functions of the scalar field. To check the accuracy of approximations considered we construct inflationary models with quadratic and quartic monomial potentials and the Gauss-Bonnet term. Numerical analysis of these models indicates that the proposed inflationary scenarios do not contradict to the observation data. New slow-roll approximations show that the constructed inflationary models are in agreement with the observation data, whereas one does not get allowed observational parameters at the same values of parameters of the constructed models in the standard slow-roll approximation.

We propose a setup for the origin of dark matter based on spacetime with a warped extra dimension and three branes: the Planck brane, the TeV brane, at a (few) TeV scale $\rho_T$, and a dark brane, at a (sub)-GeV scale $\rho_1\lesssim 100$ GeV $\ll\rho_T$. The Standard Model is localized in the TeV brane, thus solving the Higgs hierarchy problem, while the dark matter $\chi$, a Dirac fermion with mass $m_\chi<\rho_1$, is localized in the dark brane. The radion, with mass $m_r<m_\chi$, interacts strongly ($\sim m_\chi/\rho_1\sim\mathcal O(1)$) with dark matter and very weakly ($\sim m_{f}\rho_1/\rho_T^2\ll 1$) with the Standard Model matter $f$. The generic conflict between the bounds on its detection signatures and its proper relic abundance is avoided as dark matter annihilation is $p$-wave suppressed. The former is determined by its very weak interactions with the SM and the latter by its much stronger annihilation into radions. Therefore, there is a vast range in the Dark Matter's parameter space where the correct relic abundance is achieved consistently with the existing bounds. Moreover, for the dark brane with $\rho_1\lesssim 10$ GeV, a confinement/deconfinement first order phase transition, where the radion condensates, produces a stochastic gravitational waves background at the nanoHz frequencies, which can be identified with the signal detected by the Pulsar Timing Array (PTA) experiments. In the PTA window, for $0.15 \textrm{ GeV}\lesssim m_\chi\lesssim 0.5$ GeV the relic abundance is reproduced and all constraints are satisfied.

In this paper, we compute departures in the black hole thermodynamics induced by either geometric or topological corrections to general relativity. Specifically, we analyze the spherically symmetric spacetime solutions of two modified gravity scenarios with Lagrangians $\mathcal{L}\sim R^{1+\epsilon}$ and $\mathcal{L}\sim R+\epsilon\, \mathcal{G}^2$, where $\mathcal{G}$ is the Euler density in four dimensions, while $ 0<\epsilon\ll 1$ measures the perturbation around the Hilbert-Einstein action. Accordingly, we find the expressions of the Bekenstein-Hawking entropy by the Penrose formula, and the black hole temperature and horizon of the obtained solutions. We then investigate the heat capacities in terms of the free parameters of the theories under study. In doing so, we show that healing the problem of negative heat capacities can be possible under particular choices of the free constants, albeit with limitations on the masses allowed for the black hole solutions.

Incorporating first-order QED effects, we explore the shadows of Kerr-Newman black holes with a magnetic charge through the numerical backward ray-tracing method. Our investigation accounts for both the direct influence of the electromagnetic field on light rays and the distortion of the background spacetime metric due to QED corrections. We notice that the area of the shadow increases with the QED effect, mainly due to the fact that the photons move more slowly in the effective medium and become easier to be trapped by the black hole.