Papers for Tuesday, Apr 26 2022

The ionosphere introduces chromatic distortions on low frequency radio waves, and thus poses a hurdle for 21-cm cosmology. In this paper we introduce time-varying chromatic ionospheric effects on simulated antenna temperature data of a global 21-cm data analysis pipeline, and try to detect the injected global signal. We demonstrate that given turbulent ionospheric conditions, more than 5\% error in our knowledge of the ionospheric parameters could lead to comparatively low evidence and high root-mean-square error (RMSE), suggesting a false or null detection. When using a constant antenna beam for cases that include data at different times, the significance of the detection lowers as the number of time samples increases. It is also shown that for observations that include data at different times, readjusting beam configurations according to the time-varying ionospheric conditions should greatly improve the significance of a detection, yielding higher evidences and lower RMSE, and that it is a necessary procedure for a successful detection when the ionospheric conditions are not ideal.

A quantum-mechanical solution to the problem of radiative recombination of an electron in the Coulomb field has long been known. However, in astrophysics, the classical approach is sometimes used to describe similar processes in systems of magnetic monopoles or interacting dark matter particles. The importance of such problems is determined by the fact that recombination processes play a decisive role in the evolution of the large-scale structure of the Universe. It is shown that the applicability of the quantum and classical approaches to radiative recombination is closely related to the magnitude of the radiated angular momentum and its quantization. For situations where the classical approach is not suitable, a semiclassical approach based on consideration of the quantization of the angular momentum is proposed. It can be useful for an approximate description of radiative recombination in a variety of systems.

Polar faculae (PFe) are the footpoints of magnetic field lines near the Sun's poles that are seen as bright regions along the edges of granules. The time variation in the number of PFe has been shown to correlate with the strength of the polar magnetic field and to be a predictor of the subsequent solar cycle. Due to the small size and transient nature of these features, combined with different techniques and observational factors, previous counts of PFe differ in magnitude. Further, there were no scalable techniques to measure the statistical properties of faculae. Using data from the Helioseismic and Magnetic Imager (HMI) onboard the Solar Dynamics Observatory (SDO), we present two new methods for tracking faculae and measuring their properties. In the first, we calculate the standard deviation of the HMI images over one day, visualizing the faculae as streaks. The facular lifetime is found by dividing the angular streak length by the latitude-dependent rotation rate. We apply this method to 134 days of data over 11 years. The distribution of facular lifetimes has a mean of 6.0 hours, a FWHM of 5.4 hours, and a skew towards longer lifetimes, with some lasting up to 1 day. In the second method, we overlay images of the progressive standard deviation with the HMI magnetogram to show the close relationship between facular candidates and the magnetic field. The results allow us to distinguish between motion due to the Sun's rotation and "proper motion" due to faculae moving across the Sun's surface, confirming that PFe participate in convective motions at the poles. Counts of PFe using both methods agree with previous counts in their variation with the solar cycle and the polar magnetic field. These methods can be extended to automate the identification and measurement of other properties of of PFe, which would allow for daily measurements of all faculae since SDO began operation in 2010.

The three rung distance ladder, which calibrates Type Ia supernovae through stellar distances linked to geometric measurements, provides the highest precision direct measurement of the Hubble constant. In light of the Hubble tension, it is important to test the individual components of the distance ladder. For this purpose, we report a measurement of the Hubble constant from 35 extragalactic Cepheid hosts measured by the SH0ES team, using their distances and redshifts at cz < 3300 km /s , instead of any, more distant Type Ia supernovae, to measure the Hubble flow. The Cepheid distances are calibrated geometrically in the Milky Way, NGC 4258, and the Large Magellanic Cloud. Peculiar velocities are a significant source of systematic uncertainty at z $\sim$ 0.01, and we present a formalism for both mitigating and quantifying their effects, making use of external reconstructions of the density and velocity fields in the nearby universe. We identify a significant source of uncertainty originating from different assumptions about the selection criteria of this sample, whether distance or redshift limited, as it was assembled over three decades. Modeling these assumptions yields central values ranging from H0 = 71.8 to 77.0 km/s/Mpc. Combining the four best fitting selection models yields H0 = 73.1 (+2.6/-2.3) km/s/Mpc as a fiducial result, at $2.6\sigma$ tension with Planck. While Type Ia supernovae are essential for a precise measurement of H0, unknown systematics in these supernovae are unlikely to be the source of the Hubble tension

We study the intrinsic 3D shapes of quiescent galaxies over the last half of cosmic history based on their axial ratio distribution. To this end, we construct a sample of unprecedented size, exploiting multi-wavelength $u$-to-$K_s$ photometry from the deep wide area surveys KiDS+VIKING paired with high-quality $i$-band imaging from HSC-SSP. Dependencies of the shapes on mass, redshift, photometric bulge prominence and environment are considered. For comparison, the intrinsic shapes of quenched galaxies in the IllustrisTNG simulations are analyzed and contrasted to their formation history. We find that over the full $0<z<0.9$ range, and in both simulations and observations, spheroidal 3D shapes become more abundant at $M_* > 10^{11}\ M_{\odot}$, with the effect being most pronounced at lower redshifts. In TNG, the most massive galaxies feature the highest ex-situ stellar mass fractions, pointing to violent relaxation via mergers as the mechanism responsible for their 3D shape transformation. Larger differences between observed and simulated shapes are found at low to intermediate masses. At any mass, the most spheroidal quiescent galaxies in TNG feature the highest bulge mass fractions, and conversely observed quiescent galaxies with the highest bulge-to-total ratios are found to be intrinsically the roundest. Finally, we detect an environmental influence on galaxy shape, at least at the highest masses, such that at fixed mass and redshift quiescent galaxies tend to be rounder in denser environments.

Recent work investigating the impact of winds and outflows from active galactic nuclei (AGN) on the habitability of exoplanets suggests that such activity could be deleterious for the long-term survival of planetary atmospheres and the habitability of planets subject to such winds. Here, we discuss the relative importance of the effect of AGN winds compared to stellar winds and the effect of the planet's magnetosphere and stellar irradiation and conclude that AGN winds are not likely to play a significant role in the evolution of atmospheric conditions in planets under conditions otherwise favorable for habitability.

We present an improved inverse ray-shooting code based on GPUs for generating microlensing magnification maps. In addition to introducing GPUs for acceleration, we put the efforts in two aspects: (i) A standard circular lens plane is replaced by a rectangular one to reduce the number of unnecessary lenses as a result of an extremely prolate rectangular image plane. (ii) Interpolation method is applied in our implementation which has achieved an significant acceleration when dealing with large number of lenses and light rays required by high resolution maps. With these applications, we have greatly reduced the running time while maintaining high accuracy: the speed has been increased by about 100 times compared with ordinary GPU based IRS code and GPU-D code when handling large number of lenses. If encountered the high resolution situation up to $10000^2$ pixels, resulting in almost $10^{11}$ light rays, the running time can also be reduced by two orders of magnitude.

The new Gaia data release (EDR3) with improved astrometry has opened a new era in studying our Milky Way in fine detail. We use Gaia EDR3 astrometry together with 2MASS and WISE photometry to study two of the most massive molecular clouds in the solar vicinity: Orion A and California. Despite having remarkable similarities in the plane of the sky in terms of shape, size, and extinction, California has an order of magnitude lower star formation efficiency. We use our state-of-the-art dust mapping technique to derive the detailed three-dimensional (3D) structure of the two clouds, taking into account both distance and extinction uncertainties, and a full 3D spatial correlation between neighbouring points. We discover that, despite the apparent filamentary structure in the plane of the sky, California is a flat 120-pc-long sheet extending from 410 to 530 $pc$. We show that not only Orion A and California differ substantially in their 3D shapes, but also Orion A has considerably higher density substructures in 3D than California. This result presents a compelling reason why the two clouds have different star formation activities. We also demonstrate how the viewing angle of California can substantially change the cloud's position in the Kennicutt-Schmidt relation. This underlines the importance of 3D information in interpreting star formation relations and challenges studies that rely solely on the column density thresholds to determine star formation activities in molecular clouds. Finally, we provide accurate distance estimates to multiple lines of sight towards various parts of the two clouds.

The coupling between radio and X-ray luminosity is an important diagnostic tool to study the connection between the accretion inflow and jet outflow for low-mass X-ray binaries (LMXBs). When comparing NS- and BH-LMXBs, we find that the radio/X-ray correlation for individual NS-LMXBs is scattered, whereas for individual BH-LMXBs a more consistent correlation is generally found. Furthermore, we observe jet quenching for both types of LMXBs, but it is unclear what exactly causes this, and if jets in NS-LMXBs quench as strongly as those in BH-LMXBs. While additional soft X-ray spectral components can be present for NS-LMXBs due to the presence of the neutron star's surface, disentangling the individual X-ray spectral components has thus far not been considered when studying the radio/X-ray coupling. Here we present eleven epochs of Swift/XRT observations matched with quasi-simultaneous archival radio observations of the 2009 November outburst of Aql X-1. We decompose the thermal and Comptonised spectral components in the Swift/XRT spectra, with the aim of studying whether the presence of additional thermal emission affects the coupling of the radio/X-ray luminosity. We find that there is no evidence of a significant thermal contribution in Swift/XRT spectra that could cause scatter in the radio/X-ray coupling. However, a noteworthy finding is that the X-ray observation with the strongest statistically significant thermal component occurs around the same epoch as a bright radio detection. Follow-up research using more sensitive X-ray observations combined with densely sampled near-simultaneous radio observations is required to explore the radio/X-ray coupling for NS-LMXBs in more detail.

A key item of interest for planetary scientists and astronomers is the habitable zone, or the distance from a host star where a terrestrial planet can maintain necessary temperatures in order to retain liquid water on its surface. However, when observing a system's habitable zone, it is possible that one may instead observe a Venus-like planet. We define "Venus-like" as greenhouse-gas-dominated atmosphere occurring when incoming solar radiation exceeds infrared radiation emitted from the planet at the top of the atmosphere, resulting in a runaway greenhouse. Our definition of Venus-like includes both incipient and post-runaway greenhouse states. Both the possibility of observing a Venus-like world and the possibility that Venus could represent an end-state of evolution for habitable worlds, requires an improved understanding of the Venus-like planet; specifically, the distances where these planets can exist. Understanding this helps us define a "Venus zone", or the region in which Venus-like planets could exist, and assess the overlap with the aforementioned "Habitable Zone". In this study, we use a 1D radiative-convective climate model to determine the outer edge of the Venus zone for F0V, G2V, K5V, and M3V and M5V stellar spectral types. Our results show that the outer edge of the Venus zone resides at 3.01, 1.36, 0.68, 0.23, and 0.1 AU, respectively. These correspond to incident stellar fluxes of 0.8, 0.55, 0.38, 0.32, and 0.3 S, respectively, where stellar flux is relative to Earth (1.0). These results indicate that there may be considerable overlap between the habitable zone and the Venus zone.

It is widely believed that at least some fast radio bursts (FRBs) are produced by magnetars. Even though mounting observational evidence points towards a magnetospheric origin of FRB emission, the question of the location for FRB generation continues to be debated. One argument recently suggested against the magnetospheric origin of bright FRBs is that the radio waves associated with an FRB may lose most of their energy before escaping the magnetosphere because the cross-section for $e^\pm$ to scatter large-amplitude EM waves in the presence of a strong magnetic field is much larger than the Thompson cross-section. We have investigated this suggestion in this work, and find that FRB radiation traveling through the open field line region of a magnetar's magnetosphere does not suffer much loss due to two previously ignored factors. First, the plasma in the outer magnetosphere ($r \gta 10^9 \ $cm), where the losses are potentially most severe, is likely to be flowing outward at a high Lorentz factor $\gamma_p \geq 10^3$. Second, the angle between the wave vector and the magnetic field vector, $\theta_B$, in the outer magnetosphere is likely of the order of 0.1 radian or smaller due in part to the intense FRB pulse that tilts open magnetic field lines so that they get aligned with the pulse propagation direction. Both these effects reduce the interaction between the FRB pulse and the plasma substantially. We find that a bright FRB with an isotropic luminosity $L_{\rm frb} \gta 10^{42} \ {\rm erg \ s^{-1}}$ can escape the magnetosphere unscathed for a large section of the $\gamma_p-\theta_B$ parameter space, and therefore conclude that the generation of FRBs in magnetar magnetosphere passes this test.

We investigate the impact of strong initial magnetic fields in core-collapse supernovae of non-rotating progenitors by simulating the collapse and explosion of a 16.9 Msun star for a strong- and weak-field case assuming a twisted-torus field with initial central field strengths of ~1012 G and ~106 G. The strong-field model has been set up with a view to the fossil-field scenario for magnetar formation and emulates a pre-collapse field configuration that may occur in massive stars formed by a merger. This model undergoes shock revival already 100 s after bounce and reaches an explosion energy of 9.3 x 10^50 erg at 310 ms, in contrast to a more delayed and less energetic explosion in the weak-field model. The strong magnetic fields help trigger a neutrino-driven explosion early on, which results in a rapid rise and saturation of the explosion energy. Dynamically, the strong initial field leads to a fast build-up of magnetic fields in the gain region to 40% of kinetic equipartition and also creates sizable pre-shock ram pressure perturbations that are known to be conducive to asymmetric shock expansion. For the strong-field model, we find an extrapolated neutron star kick of ~350 km/s, a spin period of ~70 ms, and no spin-kick alignment. The dipole field strength of the proto-neutron star is 2 x 10^14 G by the end of the simulation with a declining trend. Surprisingly, the surface dipole field in the weak-field model is stronger, which argues against a straightforward connection between pre-collapse fields and the birth magnetic fields of neutron stars.

We investigate the properties of stars participating in double compact star merger events considering interacting model of stable strange quark matter. We model the matter making it compatible with the recent astrophysical observations of compact star mass-radius and gravitational wave events. In this context we consider modified MIT bag model and vector bag model with and without self interaction. We find new upper bound on tidal deformability of $1.4~M_\odot$ strange star corresponding to the upper bound of effective tidal deformability inferred from gravitational wave event. Range of compactness of $1.4~M_\odot$ strange star is obtained as ${0.175}\leq{C_{1.4}}\leq{0.199}$. Radius range of $1.5M_\odot$ primary star is deduced to be ${10.57}\leq{R_{1.5}}\leq{12.04}$ km, following stringent GW170817 constraints. GW190425 constraints provide with upper limit on radius of $1.7$ solar mass strange star that it should be less than $13.41$ $\text{km}$.

Changes in the supercycle lengths of some SU UMa-type dwarf novae have been detected by other studies, and indicate that the mass transfer rates noticeably decrease over time. We investigated the supercycle lengths of three SU UMa-type dwarf novae: AR Pic, QW Ser and V521 Peg, to determine if they have detectable changes in their supercycles. We present the results of optical spectroscopic and photometric observations of these sources. Our observations were conducted in 2016 and 2017 at the Boyden Observatory and the Sutherland station of the South African Astronomical Observatory. The quiescent results indicated that all three sources are typical SU UMa-type dwarf novae. We also present results of AR Pic and QW Ser in outburst and of V521 Peg during a precursor outburst and superoutburst. Light curves were supplemented by the Catalina Real-Time Transient Survey, the ASAS-3 and ASAS-SN archives, and the AAVSO International database in order to investigate the long-term behavior of these sources. Our results combined with catalogued properties for all short-period dwarf novae, shows a possible relationship between the supercycle time in SU UMa systems and their orbital periods, which is interpreted as the decline in the mass transfer rate as systems evolve toward and away from the 'period minimum'. At the shortest orbital periods, SU UMa systems are almost indistinguishable from WZ Sge systems. However, we propose that the scaleheight between the secondary's photosphere and L1 may be a factor that distinguish the SU UMa subclasses.

We present the results of fast spectrophotometry of flares on the RS CVn-type star AR Lac with a time resolution of 34 s and a spectroscopic resolution R $\sim $ 1300. The observations were performed on July 21-22, 2021 with the 2.0 m Karl Zeiss telescope at the Terskol Observatory. During the flares, an additional emission appeared in the spectrum of AR Lac at wavelengths nearby Ca II H\&K ($\lambda $ = 3933, 3968 \AA), nearby Si IV lines ($\lambda $ = 4089, 4116 \AA). Variations in these lines range from 2 to 5\%. We have estimated the UBV magnitudes from spectrograms through a mathematical convolution of the spectra with the filter transmission curves. The flare amplitudes in the UBV band up to 0.8 magnitudes were discovered. A detailed colorimetric analysis has allowed important parameters of the flares on AR Lac to be estimated: the temperatures at maximum light and their sizes. The color-color $(U-B)$ - $(B-V)$ diagrams confirm that all the flares at maximum light radiate as a blackbody. The temperatures at maximum light were up to 12000 $\pm $ 300 K. Based on our colorimetric analysis, we have estimated the linear sizes of the flares at maximum light. The linear size of the flares at the maximum luminosity is approximately 2\% of the radius of the star.

Whiting 1 is a faint and young globular cluster in the halo of the Milky Way, and was suggested to have originated in the Sagittarius spherical dwarf galaxy (Sgr dSph). In this paper, we use the deep DESI Legacy Imaging Surveys to explore tentative spatial connection between Whiting 1 and the Sgr dSph. We redetermine the fundamental parameters of Whiting 1 and use the best-fitting isochrone (age $\tau$=6.5 Gyr, metalicity Z=0.005 and $\rm d_{\odot}$=26.9 kpc) to construct a theoretical matched filter for the extra-tidal features searching. Without any smooth technique to the matched filter density map, we detect a round-shape feature with possible leading and trailing tails on either side of the cluster. This raw image is not totally new compared to old discoveries, but confirms that no more large-scale features can be detected under a depth of r<=22.5 mag. In our results, the whole feature stretches 0.1-0.2 degree along the orbit of Whiting 1, which gives a much larger area than the cluster core. The tails on both sides of the cluster align along the orbital direction of the Sgr dSph as well as the cluster itself, which implies that these debris are probably stripped remnants of Whiting 1 by the Milky Way.

We model long-term variations in the scintillation of binary pulsar PSR J1603$-$7202, observed by the 64 m Parkes radio telescope (Murriyang) between 2004 and 2016. We find that the time variation in the scintillation arc curvature is well-modelled by scattering from an anisotropic thin screen of plasma between the Earth and the pulsar. Using our scintillation model, we measure the inclination angle and longitude of ascending node of the orbit, yielding a significant improvement over the constraints from pulsar timing. From our measurement of the inclination angle, we place a lower bound on the mass of J1603$-$7202's companion of $\gtrsim 0.5\,\text{M}_\odot$ assuming a pulsar mass of $\gtrsim1.2\,\text{M}_\odot$. We find that the scintillation arcs are most pronounced when the electron column density along the line of sight is increased, and that arcs are present during a known extreme scattering event. We measure the distance to the interstellar plasma and its velocity, and we discuss some structures seen in individual scintillation arcs within the context of our model.

The presence of forming planets embedded in their protoplanetary disks has been inferred from the detection of multiring structures in such disks. Most of these suspected planets are undetectable by direct imaging observations at current measurement sensitivities. Inward migration and accretion might make these putative planets accessible to the Doppler method, but the actual extent of growth and orbital evolution remains unconstrained. Under the premise that the gaps in the disk around HD 163296 originate from new-born planets, we investigate if and under which circumstances the gap-opening planets could represent progenitors of the exoplanet population detected around A-type stars. In particular, we study the dependence of final planetary masses and orbital parameters on the viscosity of the disk. The evolution of the embedded planets was simulated throughout the disk lifetime and up to 100 Myr after the dispersal of the disk, taking the evolving disk structure and a likely range of disk lifetimes into account. We find that the final configuration of the planets is largely determined by the $\alpha$ viscosity parameter of the disk and less dependent on the choice for the disk lifetime and the initial planetary parameters. If we assume that planets such as those in HD 163296 evolve to form the observed exoplanet population of A-type stars, a $\alpha$ parameter on the order of $3.16 \times 10^{-4} \lesssim \alpha \lesssim 10^{-3}$ is required for the disks to induce sufficiently high migration rates. Depending on whether or not future direct imaging surveys will uncover a larger number of planets with $m_\mathrm{pl} \lesssim 3 M_\mathrm{Jup}$ and $a_\mathrm{pl} \gtrsim 10 \mathrm{AU}$ we expect the $\alpha$ parameter to be at the lower or upper end of this range, always under the assumption that such disks indeed harbor wide orbit planets.

A millisecond magnetar engine has been widely suggested to exist in gamma-ray burst (GRB) phenomena, in view of its substantial influences on the GRB afterglow emission. In this paper, we investigate the effects of the magnetar engine on the supernova (SN) emission which is associated with long GRBs and, specifically, confront the model with the observational data of SN 2006aj/GRB 060218. SN 2006aj is featured by its remarkable double-peaked ultraviolet-optical (UV-opt) light curves. By fitting these light curves, we demonstrate that the first peak can be well accounted for by the breakout emission of the shock driven by the magnetar wind, while the primary supernova emission is also partly powered by the energy injection from the magnetar. The magnetic field strength of the magnetar is constrained to be $\sim 10^{15}$ G, which is in good agreement with the common results inferred from the afterglow emission of long GRBs. In more detail, it is further suggested that the UV excess in the late emission of the supernova could also be due to the leakage of the non-thermal emission of the pulsar wind nebula (PWN), if some special conditions can be satisfied. The consistency between the model and the SN 2006aj observation indicates that the magnetar engine is likely to be ubiquitous in the GRB phenomena and even further intensify their connection with the phenomena of superluminous supernovae.

The Fix and Pair techniques were designed to generate simulations with reduced variance in the 2-point statistics by modifying the Initial Conditions (ICs). In this paper we show that this technique is also valid when the initial conditions have local Primordial non-Gaussianities (PNG), parametrised by $f_{\rm NL}$, without biasing the 2-point statistics but reducing significantly their variance. We show how to quantitatively use these techniques to test the accuracy of galaxy/halo clustering models down to a much reduced uncertainty and we apply them to test the standard model for halo clustering in the presence of PNG. Additionally, we show that by Matching the stochastic part of the ICs for two different cosmologies (Gaussian and non-Gaussian) we obtain a large correlation between the (2-point) statistics that can explicitly be used to further reduce the uncertainty of the model testing. For our reference analysis ($f_{\rm NL}=100$, $V=1 [h^{-1}{\rm Gpc}]^3$, $n= 2.5\times 10^{-4}[h^{-1}{\rm Mpc}]^{-3}$, $b=2.32$), we obtain an uncertainty of $\sigma(f_{\rm NL})=60$ with a standard simulation, whereas using Fixed [Fixed-Paired] initial conditions it reduces to $\sigma(f_{\rm NL})=12$ [$\sigma(f_{\rm NL})=12$]. When also Matching the ICs we obtain $\sigma (f_{\rm NL})=18$ for the standard case, and $\sigma (f_{\rm NL})=8$ [$\sigma (f_{\rm NL})=7$] for Fixed [Fixed-Paired]. The combination of the Fix, Pair and Match techniques can be used in the context of PNG to create simulations with an effective volume incremented by a factor $\sim 70$ at given computational resources.

Properties of the extragalactic magnetic field (EGMF) outside of clusters and filaments of the large-scale structure are largely unknown. The EGMF could be probed with $\gamma$-ray observations of distant (redshift $z > 0.1$) blazars. TeV $\gamma$ rays from these sources strongly absorb on extragalactic background light photons; secondary electrons and positrons produce cascade $\gamma$ rays with the observable flux dependent on EGMF parameters. We put constraints on the EGMF strength using 145 months of Fermi-LAT observations of the blazars 1ES 1218+304, 1ES 1101-232, and 1ES 0347-121, and imaging atmospheric Cherenkov telescope observations of the same sources. We perform a series of full direct Monte Carlo simulations of intergalactic electromagnetic cascades with the ELMAG 3.01 code and calculate the observable spectra inside the point spread functions of the observing instruments for a range of EGMF strengths. We compare the combined observed spectra with the models for various values of the EGMF strength $B$ and calculate the exclusion statistical significance for every value of $B$. We find that the values of the EGMF strength $B \le 10^{-17}$ G are excluded at a high level of statistical significance $Z > 4 \sigma$ for all the four options of the intrinsic spectral shape (power-law, power-law with exponential cutoff, log-parabola, log-parabola with exponential cutoff). On the other hand, $B = 10^{-16}$ G is still a viable option of the EGMF strength. These results were obtained for the case of stable sources.

After studying the design geometry of the Antikythera Mechanism Saros spiral, new critical geometrical/mechanical characteristics of the Back plate design were detected. The geometrical characteristics related to the symmetry of the Antikythera Mechanism design, are independent to the present irregular deformation of the Mechanism parts and were used as calibration points for the Saros spiral cells positional measurements. The Saros cells numbering was recalculated using the calibration points position. A correction of minus one to the currently accepted numbering of the Saros cells was applied. Following the new numbering, a new proper position for the (displaced) Saros pointer axis-g, in graphic design environment was calculated. The measurements were tested on a bronze reconstruction of the Back plate, by the authors. This research leads to a new important result that the Saros does not start in a random or arbitrary date but only when a solar eclipse occurs within a month. Additional results were also calculated regarding the symmetry of the eclipse events/sequence. The new Saros cell numbering strongly affects the calculations for the initial starting date of the Saros spiral and the eclipse events scheme of the Antikythera Mechanism.

COIN-Gaia 13 is a newly discovered open cluster revealed by Gaia DR2 data. It is a nearby open cluster with a distance of about 513 pc. Combined with the five-dimensional astrometric data of Gaia EDR3 with higher accuracy, we use the membership assignment algorithm (pyUPMASK) to determine the membership of COIN-Gaia 13 in a large extended spatial region. The cluster has found 478 candidate members. After obtaining reliable cluster members, we further study its basic properties and spatial distribution. Our results show that there is an obvious extended structure of the cluster in the X-Y plane. This elongated structure is distributed along the spiral arm, and the whole length is about 270 pc. The cluster age is 250 Myr, the total mass is about 439 M$_\odot$, and the tidal radius of the cluster is about 11 pc. Since more than half of the member stars (352 stars) are located outside twice the tidal radius, it is suspected that this cluster is undergoing the dynamic dissolution process. Furthermore, the spatial distribution and kinematic analysis indicate that the extended structure in COIN-Gaia 13 is more likely to be caused by the differential rotation of the Galaxy.

Using photometric data from the Dark Energy Survey and the Wide-field Infrared Survey Explorer, we estimate photometric redshifts for 105 million galaxies using the nearest-neighbour algorithm. From such a large database, 151,244 clusters of galaxies are identified in the redshift range of 0.1<z<1.5 based on the overdensity of the total stellar mass of galaxies within a given photometric redshift slice, among which 76,826 clusters are newly identified and 30,477 clusters have a redshift z>1. We cross-match these clusters with those in the catalogues identified from the X-ray surveys and the Sunyaev-Zel'dovich (SZ) effect by the Planck, South Pole Telescope and Atacama Cosmology Telescope surveys, and get the redshifts for 45 X-ray clusters and 56 SZ clusters. More than 95% SZ clusters in the sky region have counterparts in our catalogue. We find multiple optical clusters in the line of sight towards about 15% of SZ clusters.

We present a numerical evidence supporting the scenario that the peculiar alignments of the galaxy stellar spins with the major principal axes of the local tidal tensors are produced during the quiescent evolution period when the galaxies experience no recent merger events. Analyzing the merger tree from the TNG300-1 simulation of the IllustrisTNG project, we find the latest merger epochs, $a(z_{m})$, of the galaxies, and create four $a(z_{m})$-selected samples that are controlled to share the identical mass and density distributions. For each sample, we determine the spin and shape vectors of the galaxy stellar, cold and hot gas, and dark matter components separately, and compute the average strengths of their alignments with the principal directions of the local tidal fields as well as their mutual alignment tendencies. It is found that the stellar (cold gas) spin axes of the galaxies whose latest merger events occur at earlier epochs are more strongly aligned (weakly anti-aligned) with the major principal axes of the tidal fields. It is also shown that although the mass-dependent transition of the galaxy DM spins have little connection with the merger events, the morphologies, spin-shape and shape-shear alignment strengths of the galaxy four components sensitively depend on $a(z_{m})$. Noting that the stellar components of the galaxies which undergo long quiescent evolution have distinctively oblate shapes and very strong spin-shape alignments, we suggest that the local tidal field might be traced by using the stellar shapes of galaxies without signatures of mergers as a proxy of their stellar spins.

Context. The chemical building blocks of life contain a large proportion of nitrogen, an essential element. Titan, the largest moon of Saturn, with its dense atmosphere of molecular nitrogen and methane, offers an exceptional opportunity to explore how this element is incorporated into carbon chains through atmospheric chemistry in our Solar System. A brownish dense haze is consistently produced in the atmosphere and accumulates on the surface on the moon. This solid material is nitrogen-rich and may contain prebiotic molecules carrying nitrogen. Aims. To date, our knowledge of the processes leading to the incorporation of nitrogen into organic chains has been rather limited. In the present work, we investigate the formation of nitrogen-bearing ions in an experiment simulating Titan s upper atmosphere, with strong implications for the incorporation of nitrogen into organic matter on Titan. Methods. By combining experiments and theoretical calculations, we show that the abundant N2+ ion, produced at high altitude by extreme-ultraviolet solar radiation, is able to form nitrogen-rich organic species. Results. An unexpected and important formation of CH3N2+ and CH2N2+ diazo-ions is experimentally observed when exposing a gas mixture composed of molecular nitrogen and methane to extreme-ultraviolet radiation. Our theoretical calculations show that these diazo-ions are mainly produced by the reaction of N2+ with CH3 radicals. These small nitrogen-rich diazo-ions, with a N/C ratio of two, appear to be a missing link that could explain the high nitrogen content in Titan s organic matter. More generally, this work highlights the importance of reactions between ions and radicals, which have rarely been studied thus far, opening up new perspectives in astrochemistry.

Observations of the Milky Way's low-$\alpha$ disk show that at fixed metallicity, [Fe/H], several element abundance, [X/Fe], correlate with age, with unique slopes and small scatters around the age-[X/Fe] relations. In this study, we turn to simulations to explore the age-[X/Fe] relations for the elements C, N, O, Mg, Si, S, and Ca that are traced in a FIRE-2 cosmological zoom-in simulation of a Milky Way-like galaxy, m12i, and understand what physical conditions give rise to the observed age-[X/Fe] trends. We first explore the distributions of mono-age populations in their birth and current locations, [Fe/H], and [X/Fe], and find evidence for inside-out radial growth for stars with ages < 7 Gyr. We then examine the age-[X/Fe] relations across m12i's disk and find that the direction of the trends agree with observations, apart from C, O, and Ca, with remarkably small intrinsic scatters, $\sigma_{int}$ (0.01-0.04 dex). This $\sigma_{int}$ measured in the simulations is also metallicity-dependent, with $\sigma_{int}$ $\approx$ 0.025 dex at [Fe/H]=-0.25 dex versus $\sigma_{int}$ $\approx$ 0.015 dex at [Fe/H]=0 dex, and a similar metallicity dependence is seen in the GALAH survey for the elements in common. Additionally, we find that $\sigma_{int}$ is higher in the inner galaxy, where stars are older and formed in less chemically-homogeneous environments. The age-[X/Fe] relations and the small scatter around them indicate that simulations capture similar chemical enrichment variance as observed in the Milky Way, arising from stars sharing similar element abundances at a given birth place and time.

The light curve of the microlensing event KMT-2021-BLG-1898 exhibits a short-term central anomaly with double-bump features that cannot be explained by the usual binary-lens or binary-source interpretations. With the aim of interpreting the anomaly, we analyze the lensing light curve under various sophisticated models. We find that the anomaly is explained by a model, in which both the lens and source are binaries (2L2S model). For this interpretation, the lens is a planetary system with a planet/host mass ratio of $q\sim 1.5\times 10^{-3}$, and the source is a binary composed of a turn off or a subgiant star and a mid K dwarf. The double-bump feature of the anomaly can also be depicted by a triple-lens model (3L1S model), in which the lens is a planetary system containing two planets. Among the two interpretations, the 2L2S model is favored over the 3L1S model not only because it yields a better fit to the data, by $\Delta\chi^2=[14.3$--18.5], but also the Einstein radii derived independently from the two stars of the binary source result in consistent values. According to the 2L2S interpretation, KMT-2021-BLG-1898 is the third planetary lensing event occurring on a binary stellar system, following MOA-2010-BLG-117 and KMT-2018-BLG-1743. Under the 2L2S interpretation, we identify two solutions resulting from the close-wide degeneracy in determining the planet-host separation. From a Bayesian analysis, we estimate that the planet has a mass of $\sim 0.7$--0.8~$M_{\rm J}$, and it orbits an early M dwarf host with a mass of $\sim 0.5~M_\odot$. The projected planet-host separation is $\sim 1.9$~AU and $\sim 3.0$~AU according to the close and wide solutions, respectively.

The light curve of the microlensing event KMT-2021-BLG-0240 exhibits a short-lasting anomaly with complex features near the peak at the 0.1~mag level from a single-lens single-source model. We conducted modeling of the lensing light curve under various interpretations to reveal the nature of the anomaly. It is found that the anomaly cannot be explained with the usual model based on a binary-lens (2L1S) or a binary-source (1L2S) interpretation. However, a 2L1S model with a planet companion can describe part of the anomaly, suggesting that the anomaly may be deformed by a tertiary lens component or a close companion to the source. From the additional modeling, we find that all the features of the anomaly can be explained with either a triple-lens (3L1S) model or a binary-lens binary-source (2L2S) model obtained under the 3L1S interpretation. However, it is difficult to validate the 2L2S model because the light curve does not exhibit signatures induced by the source orbital motion and the ellipsoidal variations expected by the close separation between the source stars according to the model. We, therefore, conclude that the two interpretations cannot be distinguished with the available data, and either can be correct. According to the 3L1S solution, the lens is a planetary system with two sub-Jovian-mass planets in which the planets have masses of 0.32--0.47~$M_{\rm J}$ and 0.44--0.93~$M_{\rm J}$, and they orbit an M dwarf host. According to the 2L2S solution, on the other hand, the lens is a single planet system with a mass of $\sim 0.21~M_{\rm J}$ orbiting a late K-dwarf host, and the source is a binary composed of a primary of a subgiant or a turnoff star and a secondary of a late G dwarf. The distance to the planetary system varies depending on the solution: $\sim 7.0$~kpc according to the 3L1S solution and $\sim 6.6$~kpc according to the 2L2S solution.

We present the multi-wavelength counterparts of 850-$\mu$m selected submillimetre sources over a 2-deg$^2$ field centred on the North Ecliptic Pole. In order to overcome the large beam size (15 arcsec) of the 850-$\mu$m images, deep optical to near-infrared (NIR) photometric data and arcsecond-resolution 20-cm images are used to identify counterparts of submillimetre sources. Among 647 sources, we identify 514 reliable counterparts for 449 sources (69 per cent in number), based either on probabilities of chance associations calculated from positional offsets or offsets combined with the optical-to-NIR colours. In the radio imaging, the fraction of 850-$\mu$m sources having multiple counterparts is 7 per cent. The photometric redshift, infrared luminosity, stellar mass, star-formation rate (SFR), and the AGN contribution to the total infrared luminosity of the identified counterparts are investigated through spectral energy distribution fitting. The SMGs are infrared-luminous galaxies at an average $\langle z\rangle=2.5$ with $\mathrm{log}_{10} (L_\mathrm{IR}/\mathrm{L}_\odot)=11.5-13.5$, with a mean stellar mass of $\mathrm{log}_{10} (M_\mathrm{star}/\mathrm{M}_\odot)=10.90$ and SFR of $\mathrm{log}_{10} (\mathrm{SFR/M_\odot\,yr^{-1}})=2.34$. The SMGs show twice as large SFR as galaxies on the star-forming main sequence, and about 40 per cent of the SMGs are classified as objects with bursty star formation. At $z\ge4$, the contribution of AGN luminosity to total luminosity for most SMGs is larger than 30 per cent. The FIR-to-radio correlation coefficient of SMGs is consistent with that of main-sequence galaxies at $z\simeq2$.

We present a multi-zone galactic chemical evolution model for the Milky Way that takes into account the updated yields of major nucleosynthesis channels. It incorporates physical processes like radial flow of gas in the disk and infall of fresh gas, along with stellar scattering processes like radial migration. For the first time using a physically motivated model, we qualitatively reproduce the observed $([{\rm Fe/H}], [\alpha/{\rm Fe}])$ distribution of stars at different Galactic locations, along with the age-$[\alpha/{\rm Fe}]$ and age-${\rm [Fe/H]}$ relations. The model successfully reproduces the bifurcated high- and low-$[\alpha/{\rm Fe}]$ sequences and sheds light on the origin of the thick disc. We analyse the effect of physical processes individually on the observed properties of the Galaxy. We show that radial flow of gas plays an important role in establishing the radial gradient for both ${\rm [Fe/H]}$ and $[\alpha/{\rm Fe}]$. The dichotomy in $[\alpha/{\rm Fe}]$ is primarily due to the sharp fall of $[\alpha/{\rm Fe}]$ with time and influenced by star formation history. We identify at least three independent factors that govern the rate of \alphafe{} declining with time: the fraction of white dwarfs that form SNe Ia, the timescale associated with the delay time distribution of Type Ia supernovae (SNe Ia), and star formation efficiency. ${\rm [Fe/H]}$ and $[\alpha/{\rm Fe}]$ evolve rapidly at early times but change gradually for at least the last 6 Gyr as they approach equilibrium values. The final equilibrium value of $[\alpha/{\rm Fe}]$ is governed by the ratio of SNe Ia to core-collapse supernovae (CCSN).

(Abridged) In this paper, we investigate the physical conditions of [CI]-traced gas in high-mass star-forming regions by analyzing APEX [CI] 492 GHz single-pointing observations of the ATLASGAL Top100 sources along with other multi-wavelength data. Our 98 sources are clearly detected in [CI] 492 GHz emission, and the observed integrated intensities and line widths tend to increase toward evolved stages of star formation. In addition to these "main" components that are associated with the Top100 sample, 41 emission and two absorption features are identified by their velocities toward 28 and two lines of sight respectively as "secondary" components. The secondary components have systematically smaller integrated intensities and line widths than the main components. We found that [CI] 492 GHz and 13CO(2-1) are well correlated with the 13CO(2-1)-to-[CI] 492 GHz integrated intensity ratio varying from 0.2 to 5.3. In addition, we derived the H2-to-[CI] conversion factor, X(CI), by dividing 870 micron-based H2 column densities by the observed [CI] 492 GHz integrated intensities and found that X(CI) ranges from 2.3e20 to 1.3e22 with a median of 1.7e21. In contrast to the strong correlation with 13CO(2-1), [CI] 492 GHz has a scattered relation with the 870 micron-traced molecular gas. Finally, we performed LTE and non-LTE analyses of the [CI] 492 GHz and 809 GHz data for a subset of the Top100 sample and inferred that [CI] emission likely originates from warm (kinetic temperature > 60 K), optically thin (opacity < 0.5), and highly pressurized (thermal pressure ~ e5 to e8 K/cm3) regions.

Optimal filtering is the crucial technique for the data analysis of transition-edge-sensor (TES) calorimeters to achieve their state-of-the-art energy resolutions. Filtering out the `bad' data from the dataset is important because it otherwise leads to the degradation of energy resolutions, while it is not a trivial task. We propose a neural network-based technique for the automatic goodness tagging of TES pulses, which is fast and automatic and does not require bad data for training.

In the context of Open Science, provenance has become a decisive piece of information to provide along with astronomical data. Provenance is explicitly cited in the FAIR principles, that aims to make research data Findable, Accessible, Interoperable and Reusable. The IVOA Provenance Data Model, published in 2020, puts in place the foundations for structuring and managing detailed provenance information, from the acquisition of raw data, to the dissemination of final products. The ambition is to provide for each astronomical dataset a sufficiently fine grained and detailed provenance information so that end-users understand the quality, reliability and trustworthiness of the data. This would ensure that the Reusable principle is respected.

The distribution of galaxies provides an ideal laboratory to test for deviations from General Relativity. In particular, redshift-space distortions are commonly used to constrain modifications to the Poisson equation, which governs the strength of dark matter clustering. Here, we show that these constraints rely on the validity of the weak equivalence principle, which has never been tested for the dark matter component. Relaxing this restrictive assumption leads to modifications in the growth of structure that are fully degenerate with modifications induced by the Poisson equation. This in turns strongly degrades the constraining power of redshift-space distortions. Such degeneracies can however be broken and tight constraints on modified gravity can be recovered using relativistic distortions in the galaxy distribution, which will be observable by the coming generation of large-scale structure surveys.

One of the ways to understand the genesis and evolution of the universe is to know how galaxies have formed and evolved. In this regard, the study of star formation history (SFH) plays an important role in the accurate understanding of galaxies. In this paper, we used long-period variable stars (LPVs) to estimate the SFH in the Andromeda galaxy (M31). These cool stars reach their peak luminosity in the final stage of their evolution; their birth mass is directly related to their luminosity. Therefore, we construct the mass function and the star formation history using stellar evolution models.

The Galactic Center black hole and the nuclear star cluster are surrounded by a clumpy ring of gas and dust (the circumnuclear disk, CND) that rotates about them at a standoff distance of ~1.5 pc. The mass and density of individual clumps in the CND are disputed. We seek to use H$_2$ to characterize the clump size distribution and to investigate the morphology and dynamics of the interface between the ionized interior layer of the CND and the molecular reservoir lying further out (corresponding to the inner rim of the CND, illuminated in ultraviolet light by the central star cluster). We have observed two fields of approximately 20"x20" in the CND at near-infrared wavelengths with the OSIRIS spectro-imager at the Keck Observatory. These two fields, located at the approaching and receding nodes of the CND, best display this interface. Our data cover two H$_2$ lines as well as the Br$\gamma$ line (tracing H ii). We have developed the tool CubeFit, an original method to extract maps of continuous physical parameters (such as velocity field and velocity dispersion) from integral-field spectroscopy data, using regularization to largely preserve spatial resolution in regions of low signal-to-noise ratio. This original method enables us to isolate compact, bright features in the interstellar medium of the CND. Several clumps in the southwestern field assume the appearance of filaments, many of which are parallel to each other. We conclude that these clumps cannot be self-gravitating.

The composition of forming planets is strongly affected by the protoplanetary disc's thermal structure. This thermal structure is predominantly set by dust radiative transfer and viscous (accretional) heating and can be impacted by gaps - regions of low dust and gas density that can occur when planets form. The effect of variations in dust surface density on disc temperature has been poorly understood until now. In this work, we use the radiative transfer code MCMax to model the 2D dust thermal structure with individual gaps corresponding to planets with masses of 0.1 M$_J$ - 5 M$_J$ and orbital radii of 3, 5, and 10 AU. Low dust opacity in the gap allows radiation to penetrate deeper into the disc and warm the midplane by up to 16 K, but only for gaps located in the region of the disc where stellar irradiation is the dominant source of heat (here, a$\gtrsim$4 AU). In viscously-heated regions (a$\lesssim$4 AU), the midplane of the gap is relatively cooler by up to 100 K. Outside of the gap, broad radial oscillations in heating and cooling are present due to changes in disc flaring. These thermal features affect local segregation of volatile elements (H$_2$O, CH$_4$, CO2, CO) between the dust and gas. We find that icelines experience dramatic shifts relative to gapless models: up to 6.5 AU towards the star and 4.3 AU towards the midplane. While quantitative predictions of iceline deviations will require more sophisticated models which include transport and sublimation/condensation kinetics, our results provide evidence that planet-induced iceline variations represent a potential feedback from the planet onto the composition of material it is accreting.

The transition of high-mass galaxies from being blue and star forming to being red and dead is a crucial step in galaxy evolution, yet not fully understood. In this work, we use the NIHAO suite of galaxy simulations to investigate the relation between the transition time through the green valley and other galaxy properties. The typical green valley crossing time of our galaxies is approximately 400 Myr, somewhat shorter than observational estimates. The crossing of the green valley is triggered by the onset of AGN feedback and the subsequent shut down of star formation. Interestingly the time spent in the green valley is not related to any other galaxy properties, such as stellar age or metallicity, or the time at which the star formation quenching takes place. The crossing time is set by two main contributions: the ageing of the current stellar population and the residual star formation in the green valley. These effects are of comparable magnitude, while major and minor mergers have a negligible contribution. Most interestingly, we find the time that a galaxy spends to travel through the green valley is twice the $e$-folding time of the star formation quenching. This result is stable against galaxy properties and the exact numerical implementation of AGN feedback in the simulation. Assuming a typical crossing time of about one Gyr inferred from observations, our results imply that any mechanism or process aiming to quench star formation, must do it on a typical timescale of 500 Myr.

The standard model of cosmology assumes that our Universe began 14 Gyrs (billion years) ago from a singular Big Bang creation. This can explain a vast range of different astrophysical data from a handful of free cosmological parameters. However, we have no direct evidence or fundamental understanding of some key assumptions: Inflation, Dark Matter and Dark Energy. Here we review the idea that cosmic expansion originates instead from gravitational collapse and bounce. The collapse generates a Black Hole (BH) of mass $ M \simeq 5 \times 10^{22} M_{\odot}$ that formed 25~Gyrs ago. As there is no pressure support, the cold collapse can continue inside in free fall until it reaches atomic nuclear saturation (GeV), when is halted by Quantum Mechanics, as two particles cannot occupy the same quantum state. The collapse then bounces like a core-collapse supernovae, producing the Big Bang expansion. Cosmic acceleration results from the BH event horizon. During collapse, perturbations exit the horizon to re-enter during expansion, giving rise to the observed universe without the need for Inflation or Dark Energy. Using Ockham's razor, this makes the BH Universe (BHU) model more compelling than the standard singular Big Bang creation.

We investigate the habitability of hypothethical moons orbiting known exoplanets. This study focuses on big, rocky exomoons that are capable of maintaining a significant atmosphere. To determine their habitability, we calculate the incident stellar radiation and the tidal heating flux arising in the moons as the two main contributors to the energy budget. We use the runaway greenhouse and the maximum greenhouse flux limits as a definition of habitability. For each exoplanet we run our calculations for plausible ranges of physical and orbital parameters for the moons and the planet using a Monte Carlo approach. We calculate the moon habitability probability for each planet which is the fraction of the investigated cases that lead to habitable conditions. Based on our results, we provide a target list for observations of known exoplanets of which the top 10 planets have more than 50\% chance for hosting habitable moons on stable orbits. Two especially promising candidates are Kepler-62~f and Kepler-16~b, both of them with known masses and radii. Our target list can help to detect the first habitable exomoon.

We propose a simple analytical jet model of magnetic jets, utilizing conservation laws and some results of published GRMHD jet simulations, and consider radially-averaged profiles of the physical quantities. We take into account conversion of the magnetic energy flux to bulk acceleration in jets formed around rotating black holes assuming the mass continuity equation and constant jet power, leading to the Bernoulli equation. For assumed profiles of the bulk Lorentz factor and the radius, this gives the profile of toroidal magnetic field component along the jet. We then consider the case where the poloidal field component is connected to a rotating black hole surrounded by an accretion disc. The formalism recovers the standard formula for the power extracted from a rotating black hole. We find that the poloidal field strength dominates over the toroidal one in the comoving frame up to large distances, which means that jets should be more stable to current-driven kink modes. The resulting magnetic field profiles can then be used to calculate the jet synchrotron emission.

We investigate the formation of multiple images as the radio signals from fast radio bursts (FRBs) pass through the plane of a plasma clump. The exponential model for the plasma clump is adopted to analyze the properties of the multiple images. By comparing with the classical dispersion relations, we find that one image has exhibited specific inverse properties to others, such as their delay times at high frequency is higher than that at low frequency, owing to the lensing effects of the plasma clump. We demonstrate that these inverse effects should be observable in some repeating FRBs. Our results predict deviation in the estimated DM across multiple images, consistent with the observations of FRB 121102 and FRB 180916.J0158+65. If other plasma lenses have effects similar to an exponential lens, we find that they should also give rise to the similar dispersion relation in the multiple images. For some repeating FRBs, analysis of the differences in time delay and in DM between multiple images at different frequencies can serve as a method to reveal the plasma distribution.

Recent observations of coherent radiation from the Crab pulsar (Bij et al 2021) suggest the emission is driven by an ultra - relativistic ($\gamma \sim 10^4$), cold plasma flow. A relativistically expanding plasma shell can compress the ambient magnetic field, like a moving mirror, and thus produce coherent radiation whose wavelength is shorter than that of the ambient medium by $\gamma^2$. This mechanism has been studied in the past by Colgate and Noerdelinger (1971), in the context of radio loud supernova explosions. In this work we propose that a similar mechanism drives the coherent emission in fast radio bursts. The high Lorenz factors dramatically lower the implied energy and magnetic field requirements, allowing the spin down energy of regular (or even recycled), fast spinning pulsars, rather than slow spinning magnetars, to explain FRBs. We show that this model can explain the frequency and the time evolution of observed FRBs, as well as their duration, energetics and absence of panchromatic counterparts. We also predict that the peak frequency of sub pulses decline with observation time as $\omega_{\rm obs} \propto t_{\rm obs}^{-1/2}$. Unfortunately, with current capabilities it is not possible to constrain the shape of the curve $\omega_{\rm obs} \left(t_{\rm obs} \right)$. Finally, we find that a variation of this model can explain weaker radio transients, such as the one observed from a galactic magnetar. In this variant, the shock wave produces low frequency photons which are then Compton scattered to the GHz range.

During solar flares plasma is typically heated to very high temperatures, and the resulting redistribution of energy via thermal conduction is a primary mechanism transporting energy throughout the flaring solar atmosphere. The thermal flux is usually modeled using Spitzer's theory, which is based on local Coulomb collisions between the electrons carrying the thermal flux and those in the background. However, often during flares, temperature gradients become sufficiently steep that the collisional mean free path exceeds the temperature gradient scale size, so that thermal conduction becomes inherently non-local. Further, turbulent angular scattering, which is detectable in nonthermal widths of atomic emission lines, can also act to increase the collision frequency and so suppress the heat flux. Recent work by Emslie & Bian (2018) extended Spitzer's theory of thermal conduction to account for both non-locality and turbulent suppression. We have implemented their theoretical expression for the heat flux (which is a convolution of the Spitzer flux with a kernel function) into the RADYN flare-modeling code and performed a parameter study to understand how the resulting changes in thermal conduction affect flare dynamics and hence the radiation produced. We find that models with reduced heat fluxes predict slower bulk flows, less intense line emission, and longer cooling times. By comparing features of atomic emission lines predicted by the models with Doppler velocities and nonthermal line widths deduced from a particular flare observation, we find that models with suppression factors between 0.3 to 0.5 relative to the Spitzer value best reproduce observed Doppler velocities across emission lines forming over a wide range of temperatures. Interestingly, the model that best matches observed nonthermal line widths has a kappa-type velocity distribution function.

The CARMENES instrument is searching for periodic radial-velocity (RV) variations of M dwarfs, which are induced by orbiting planets. However, there are other potential sources of such variations, including rotational modulation caused by stellar activity. We aim to investigate four M dwarfs (Ross 318, YZ CMi, TYC 3529-1437-1, and EV Lac) with different activity levels and spectral sub-types. Our goal is to compare the periodicities seen in 22 activity indicators and the stellar RVs, and to examine their stability over time. For each star, we calculated GLS periodograms of pseudo-equivalent widths of chromospheric lines, indices of photospheric bands, the differential line width as a measure of the width of the average photospheric absorption line, the RV, the chromatic index that describes the wavelength dependence of the RV, and CCF parameters. We also calculated periodograms for subsets of the data and compared our results to TESS photometry. We find the rotation periods of all four stars to manifest themselves in the RV and photospheric indicators, particularly the TiO~7050 index, whereas the chromospheric lines show clear signals only at lower activity levels. For EV Lac and TYC 3529-1437-1, we find episodes during which indicators vary with the rotation period, and episodes during which they vary with half the rotation period, similarly to photometric light curves. The changing periodicities reflect the evolution of stellar activity features on the stellar surface. We therefore conclude that our results not only emphasise the importance of carefully analysing indicators complementary to the RV in RV surveys, but they also suggest that it is also useful to search for signals in activity indicators in subsets of the dataset, because an activity signal that is present in the RV may not be visible in the activity indicators all the time, in particular for the most active stars.

(abridged) Star formation in the outer Galaxy, i.e., outside of the Solar circle, has been lightly studied in part due to low CO brightness of molecular clouds linked with the negative metallicity gradient. Recent infrared surveys provide an overview of dust emission in large sections of the Galaxy, but suffer from cloud confusion and poor spatial resolution at far-infrared wavelengths. We aim to develop a methodology to identify and classify young stellar objects (YSOs) in star-forming regions in the outer Galaxy, and use it to solve a long-standing confusion with the distance and evolutionary status of IRAS 22147+5948. We use Support Vector Machine learning algorithm to complement standard color-color and color-magnitude diagrams in search for YSOs in the IRAS 22147 region using publicly available data from the `Spitzer Mapping of the Outer Galaxy' survey. The agglomerative hierarchical clustering algorithm is used to identify clusters, along with the Robitaille et al. (2017) code to calculate physical properties of individual YSOs. We identify 13 Class I and 13 Class II YSO candidates using the color-color diagrams, and additional 2 and 21 sources, respectively, using the machine learning techniques. Spectral energy distributions of 23 sources are modelled with a star and a passive disk, corresponding to Class II objects. Models of 3 sources include envelopes typical for Class I objects. The objects are grouped in 2 clusters located at the distance of ~2.2 kpc, and 5 clusters at ~5.6 kpc. The spatial extent of CO, radio continuum, and dust emission confirms the origin of YSOs in two distinct star-forming regions along a similar line-of-sight. The outer Galaxy might serve as a unique laboratory of star formation across environments on condition that complementary methods and ancillary data are used to properly account for cloud confusion and distance uncertainties.

We have analyzed the retrieval of the relation between mass and absolute Gaia G magnitude as obtained from theoretical models by \cite{2012MNRAS.427..127B} for the stars of different luminosity classes. For subgiant and main-sequence stars, we provide approximate analytical direct and inverse relations based on the most probable value of G magnitude in the course of evolution within the given stage. A comparison with (albeit few) observational data confirms that our results can be used for the estimation of the stellar mass from Gaia photometry in the range of 1 to 10 solar mass. We argue that similar relations for other luminosity classes are not informative.

We present SHINeS, a space simulator which can be used to replicate the thermal environment in the immediate neighborhood of the Sun down to a heliocentric distance r~0.06 au. The system consists of three main parts: the solar simulator which was designed and constructed in-house, a vacuum chamber, and the probing and recording equipment needed to monitor the experimental procedures. Our motivation for building this experimental system was to study the effect of intense solar radiation on the surfaces of asteroids when their perihelion distances become smaller than the semi-major axis of the orbit of Mercury. Comparisons between observational data and recent orbit and size-frequency models of the population of near-Earth asteroids suggest that asteroids are super-catastrophically destroyed when they approach the Sun. Whereas the current models are agnostic about the disruption mechanism, SHINeS was developed to study the mechanism or mechanisms responsible. The system can, however, be used for other applications that need to study the effects of high solar radiation on other natural or artificial objects.

Knowing the Galactic 3D dust distribution is relevant for understanding many processes in the interstellar medium and for correcting many astronomical observations for dust absorption and emission. Here, we aim for a 3D reconstruction of the Galactic dust distribution with an increase in the number of meaningful resolution elements by orders of magnitude with respect to previous reconstructions, while taking advantage of the dust's spatial correlations to inform the dust map. We use iterative grid refinement to define a log-normal process in spherical coordinates. This log-normal process assumes a fixed correlation structure, which was inferred in an earlier reconstruction of Galactic dust. Our map is informed through 111 Million data points, combining data of PANSTARRS, 2MASS, Gaia DR2 and ALLWISE. The log-normal process is discretized to 122 Billion degrees of freedom, a factor of 400 more than our previous map. We derive the most probable posterior map and an uncertainty estimate using natural gradient descent and the Fisher-Laplace approximation. The dust reconstruction covers a quarter of the volume of our Galaxy, with a maximum coordinate distance of $16\,\text{kpc}$, and meaningful information can be found up to at distances of $4\,$kpc, still improving upon our earlier map by a factor of 5 in maximal distance, of $900$ in volume, and of about eighteen in angular grid resolution. Unfortunately, the maximum posterior approach chosen to make the reconstruction computational affordable introduces artifacts and reduces the accuracy of our uncertainty estimate. Despite of the apparent limitations of the presented 3D dust map, a good part of the reconstructed structures are confirmed by independent maser observations. Thus, the map is a step towards reliable 3D Galactic cartography and already can serve for a number of tasks, if used with care.

The propagation direction and true velocity of a solar coronal mass ejection, which are among the most decisive factors for its geo-effectiveness, are difficult to determine through single-perspective imaging observations. Here we show that Sun-as-a-star spectroscopic observations, together with imaging observations, could allow us to solve this problem. Using observations of the Extreme-ultraviolet Variability Experiment onboard the Solar Dynamics Observatory, we found clear blue-shifted secondary emission components in extreme ultraviolet spectral lines during a solar eruption on October 28, 2021. From simultaneous imaging observations, we found that the secondary components are caused by a mass ejection from the flare site. We estimated the line-of-sight (LOS) velocity of the ejecta from both the double Gaussian fitting method and the red-blue asymmetry analysis. The results of both methods agree well with each other, giving an average LOS velocity of the plasma of $\sim 423~\rm{km~s^{-1}}$. From the $304$ \AA~image series taken by the Extreme Ultraviolet Imager onboard the Solar Terrestrial Relation Observatory-A (STEREO-A) spacecraft, we estimated the plane-of-sky (POS) velocity from the STEREO-A viewpoint {to be around $587~\rm{km~s^{-1}}$}. The full velocity of the bulk motion of the ejecta was then computed by combining the imaging and spectroscopic observations, which turns out to be around $596~\rm{km~s^{-1}}$ with an angle of $42.4^\circ$ to the west of the Sun-Earth line and $16.0^\circ$ south to the ecliptic plane.

Population studies of exoplanets are key to unlocking their statistical properties. So far the inferred properties have been mostly limited to planetary, orbital and stellar parameters extracted from, e.g., Kepler, radial velocity, and GAIA data. More recently an increasing number of exoplanet atmospheres have been observed in detail from space and the ground. Generally, however, these atmospheric studies have focused on individual planets, with the exception of a couple of works which have detected the presence of water vapor and clouds in populations of gaseous planets via transmission spectroscopy. Here, using a suite of retrieval tools, we analyse spectroscopic and photometric data of 25 hot Jupiters, obtained with the Hubble and Spitzer Space Telescopes via the eclipse technique. By applying the tools uniformly across the entire set of 25 planets, we extract robust trends in the thermal structure and chemical properties of hot Jupiters not obtained in past studies. With the recent launch of JWST and the upcoming missions Twinkle, and Ariel, population based studies of exoplanet atmospheres, such as the one presented here, will be a key approach to understanding planet characteristics, formation, and evolution in our galaxy.

The LIGO-Virgo-KAGRA (LVK) Collaboration has reported ~100 BH-BH mergers to date. The LVK provides estimates of rates, masses, effective spins, and redshifts for these mergers. Yet, the formation channel(s) of the mergers remains uncertain. One way to search for a formation site is to contrast properties of detected BH-BH mergers with different models of BH-BH merger formation. Here, we assess the usefulness of the average total mass of BH-BH mergers and its change with redshift to help identify the formation site of BH-BH mergers. We find that the average intrinsic BH-BH total merger mass shows exceptionally different behavior for the models that we adopt for our analysis. In the local universe (z=0) the average merger mass changes from Mtot=25Msun for CE binary evolution and open clusters channels, to Mtot=30Msun for the RLOF binary channel, to Mtot=45Msun for the globular cluster channel. These differences are even more pronounced at larger redshifts. However, these differences are diminished when considering LVK O3 detector sensitivity. Comparison with LVK O3 data shows that none of our adopted models can match the data despite large errors on BH-BH masses and redshifts. We emphasize that our conclusions are derived from a small set of 5 specific models that are subject to numerous known uncertainties. We also note that BH-BH mergers may originate from a mix of several channels and that other BH-BH formation channels may exist. Our study was designed to investigate the usefulness of the total BH-BH merger mass in establishing the origin of gravitational-wave sources.

Context: High contrast imaging is a powerful technique to detect and characterize planetary companions at large orbital separations from their parent stars. Aims: We aim at studying the limiting magnitude of the VLT/SPHERE Adaptive Optics system and the corresponding instrument performance for faint targets (G $\ge$ 11.0 mag). Methods: We computed coronagraphic H-band raw contrast at 300 [mas] and FWHM of the non-coronagraphic PSF, for a total of 111 different stars observed between 2016 and 2022 with IRDIS. For this, we processed a large number of individual frames that were obtained under different atmospheric conditions. We then compared the resulting raw contrast and the PSF shape as a function of the visible wave front sensor instant flux which scales with the G-band stellar magnitude. We repeated this analysis for the top 10\% and 30\% best turbulence conditions in Cerro Paranal. Results: We found a strong decrease in the coronagraphic achievable contrast for star fainter than G $\sim$ 12.5 mag, even under the best atmospheric conditions. In this regime, the AO correction is dominated by the read-out noise of the WFS detector. In particular we found roughly a factor ten decrease in the raw contrast ratio between stars with G $\sim$ 12.5 and G $\sim$ 14.0 mag. Similarly, we observe a sharp increase in the FWHM of the non-coronagraphic PSF beyond G $\sim$ 12.5 mag, and a corresponding decrease in the strehl ratio from $\sim$ 50\% to $\sim$ 20\% for the faintest stars. Although these trend are observed for the two turbulence categories, the decrease in the contrast ratio and PSF sharpness is more pronounced for poorer conditions.

We observed the nearly edge-on debris disk system HD 111520 at $J$, $H$, & $K1$ near infrared (NIR) bands using both the spectral and polarization modes of the Gemini Planet Imager (GPI). With these new observations, we have performed an empirical analysis in order to better understand the disk morphology and its highly asymmetrical nature. We find that the disk features a large brightness and radial asymmetry, most prominent at shorter wavelengths. We also find that the radial location of the peak polarized intensity differs on either side of the star by 11 AU, suggesting that the disk may be eccentric, although, such an eccentricity does not fully explain the large brightness and radial asymmetry observed. Observations of the disk halo with HST also show the disk to be warped at larger separations, with a bifurcation feature in the northwest, further suggesting that there may be a planet in this system creating an asymmetrical disk structure. Measuring the disk color shows that the brighter extension is bluer compared to the dimmer extension, suggesting that the two sides have different dust grain properties. This finding, along with the large brightness asymmetry, are consistent with the hypothesis that a giant impact occurred between two large bodies in the northern extension of the disk, although confirming this based on NIR observations alone is not feasible. Follow-up imaging with ALMA to resolve the asymmetry in the dust mass distribution is essential in order to confirm this scenario.

In this study, we model a sequence of a confined and a full eruption, employing the relaxed end state of the confined eruption of a kink-unstable flux rope as the initial condition for the ejective one. The full eruption, a model of a coronal mass ejection, develops as a result of converging motions imposed at the photospheric boundary, which drive flux cancellation. In this process, parts of the positive and negative external flux converge toward the polarity inversion line, reconnect, and cancel each other. Flux of the same amount as the canceled flux transfers to a flux rope, increasing the free magnetic energy of the coronal field. With sustained flux cancellation and the associated progressive weakening of the magnetic tension of the overlying flux, we find that a flux reduction of $\approx\!11\%$ initiates the torus instability of the flux rope, which leads to a full eruption. These results demonstrate that a homologous full eruption, following a confined one, can be driven by flux cancellation.

Both leptonic and hadronic emission processes may contribute to blazar jet emission; which dominates in blazars's high energy emission component remains an open question. Some intermediate synchrotron peaked blazars transition from their low to high energy emission components in the X-ray band making them excellent laboratories to probe both components simultaneously, and good targets for the newly launched Imaging X-ray Polarimetry Explorer. We characterize the spectral energy distributions for three such blazars: CGRaBS~J0211+1051, TXS~0506+056, and S5~0716+714, predicting their X-ray polarization behavior by fitting a multizone polarized leptonic jet model. We find that a significant detection of electron synchrotron dominated polarization is possible with a 300~ks observation for S5~0716+714 and CGRaBS~J0211+1051 in their flaring states, while even 500~ks observations are unlikely to measure synchrotron self-Compton polarization. Importantly, non-leptonic emission processes like proton synchrotron are marginally detectable for our brightest ISP, S5~0716+714, during a flaring state. Improved {\it IXPE} data reduction methods or next generation telescopes like {\it eXTP} are needed to confidently measure SSC polarization.

In the near-future, dedicated telescopes observe Earth-like exoplanets in reflected light, allowing their characterization. Because of the huge distances, every exoplanet will be a single pixel, but temporal variations in its spectral flux hold information about the planet's surface and atmosphere. We test convolutional neural networks for retrieving a planet's rotation axis, surface and cloud map from simulated single-pixel flux and polarization observations. We investigate the assumption that the planets reflect Lambertian in the retrieval while their actual reflection is bidirectional, and of including polarization in retrievals. We simulate observations along a planet's orbit using a radiative transfer algorithm that includes polarization and bidirectional reflection by vegetation, desert, oceans, water clouds, and Rayleigh scattering in 6 spectral bands from 400 to 800 nm, at various photon noise levels. The surface-types and cloud patterns of the facets covering a model planet are based on probability distributions. Our networks are trained with simulated observations of millions of planets before retrieving maps of test planets. The neural networks can constrain rotation axes with a mean squared error (MSE) as small as 0.0097, depending on the orbital inclination. On a bidirectionally reflecting planet, 92% of ocean and 85% of vegetation, desert, and cloud facets are correctly retrieved, in the absence of noise. With realistic noise, it should still be possible to retrieve the main map features with a dedicated telescope. Except for face-on orbits, a network trained with Lambertian reflecting planets, yields significant retrieval errors when given observations of bidirectionally reflecting planets, in particular, brightness artefacts around a planet's pole. Including polarization improves retrieving the rotation axis and the accuracy of the retrieval of ocean and cloud facets.

Radar and spacecraft observations show the permanently shadowed regions around Mercury's North Pole to contain water ice and complex organic material. One possible source of this material are impacts by interplanetary dust particles (IDPs), asteroids, and comets. We have performed numerical simulations of the dynamical evolution of asteroids and comets over the few Myr and checked for their impacts with Mercury. We use the N-body integrator RMVS/Swifter to propagate the Sun and the eight planets from their current positions. We add comets and asteroids to the simulations as massless test particles, based on their current orbital distributions. Asteroid impactors are assigned a probability of being water-rich (C-class) based on the measured distribution of taxonomic types. For comets, we assume a constant water fraction. For IDPs, we use a dynamical meteoroid model to compute the dust flux on Mercury. Relative to previous work on asteroid and comet impacts (Moses et al. 1999), we leverage 20 years of progress in minor body surveys. Immediate post-impact ejection of impactor material into outer space is taken into account as is the migration efficiency of water across Mercury's surface to the polar cold traps. We find that asteroids deliver $\sim 1 \times 10^{3}$ kg/yr of water to Mercury, comets deliver $\sim 1 \times 10^{3}$ kg/yr and IDPs deliver $\sim 16 \times 10^{3}$ kg/yr within a factor of several. Over a timescale of $\sim 1$ Gyr, this is enough to deliver the minimum amount of water required by the radar and MESSENGER observations. While other sources of water on Mercury are not ruled out by our analysis, we show that they are not required to explain the currently available observational lower limits.

We demonstrate how to use persistent homology for cosmological parameter inference in a tomographic cosmic shear survey. We obtain the first cosmological parameter constraints from persistent homology by applying our method to the first-year data of the Dark Energy Survey. To obtain these constraints, we analyse the topological structure of the matter distribution by extracting persistence diagrams from signal-to-noise maps of aperture masses. This presents a natural extension to the widely used peak count statistics. Extracting the persistence diagrams from the cosmo-SLICS, a suite of $N$-body simulations with variable cosmological parameters, we interpolate the signal using Gaussian Processes and marginalise over the most relevant systematic effects, including intrinsic alignments and baryonic effects. We find for the structure growth parameter $S_8=0.747^{+0.025}_{-0.031}$, which is in full agreement with other late-time probes. We also constrain the intrinsic alignment parameter to $A=1.54\pm 0.52$, ruling out the case of no intrinsic alignments at a $3\sigma$-level.

Massive vector fields feature in several areas of particle physics, e.g., as carriers of weak interactions, dark matter candidates, or as an effective description of photons in a plasma. Here we investigate vector fields with self-interactions by replacing the mass term in the Proca equation with a general potential. We show that this seemingly benign modification inevitably introduces ghost instabilities of the same kind as those recently identified for vector-tensor theories of modified gravity (but in this simpler, minimally coupled theory). It has been suggested that nonperturbative dynamics may drive systems away from such instabilities. We demonstrate that this is not the case by evolving a self-interacting Proca field on a black hole background, where it grows due to the superradiant instability. The system initially evolves as in the massive case, but instabilities are triggered in a finite time once the self-interaction becomes significant. These instabilities have implications for the formation of condensates of massive, self-interacting vector bosons, the possibility of spin-one bosenovae, vector dark matter models, and effective models for interacting photons in a plasma.

We investigate the linear stability of a sinusoidal shear flow with an initially uniform streamwise magnetic field in the framework of incompressible magnetohydrodynamics (MHD) with finite resistivity and viscosity. This flow is known to be unstable to the Kelvin-Helmholtz instability in the hydrodynamic case. The same is true in ideal MHD, where dissipation is neglected, provided the magnetic field strength does not exceed a critical threshold beyond which magnetic tension stabilizes the flow. Here, we demonstrate that including viscosity and resistivity introduces two new modes of instability. One of these modes, which we call a resistively-unstable Alfv\'en wave due to its connection to shear Alfv\'en waves, exists for any nonzero magnetic field strength as long as the magnetic Prandtl number $Pm < 1$. We present a reduced model for this instability that reveals its excitation mechanism to be the negative eddy viscosity of periodic shear flows described by Dubrulle & Frisch (1991). Finally, we demonstrate numerically that this mode saturates in a quasi-stationary state dominated by counter-propagating solitons.

A recently published study of dark matter in distant galaxies has found a direct interaction between unknown dark matter particles and ordinary galactic baryonic matter. It is proposed here that the dark matter is Rydberg Matter composed of water nanoclusters ejected from abundant amorphous water-ice-coated cosmic dust, leading to non-gravitational interactions with ordinary matter. This may also clarify recently observed CMB birefringence, suggesting possible new physics such as quintessence.

We discuss the fundamentals of classical black hole (BH) thermodynamics in a new framework determined by light surfaces and their frequencies. This new approach allows us to study BH transitions inside the Kerr geometry. In the case of BHs, we introduce a new parametrization of the metric in terms of the maximum extractable rotational energy or, correspondingly, the irreducible mass, which is an alternative to the spin parametrization. It turns out that BH spacetimes with spins $a/M= \sqrt {8/9}$ and $a/M=1/\sqrt{2}$ show anomalies in the rotational energy extraction and surface gravity whereas the case $a/M=\sqrt{3}/2$ is of particular relevance to study the variations of the horizon area. We find the general conditions under which BH transitions can occur and express them in terms of the masses of the initial and final states. This shows that BH transitions in the Kerr geometry are not arbitrary but depend on the relationship between the mass and spin of the initial and final states. From an observational point of view, we argue that near the BH poles it is possible to detect photon orbits with frequencies that characterize the light surfaces analyzed in this work.

Interacting massive spin-1 fields have been widely used in cosmology and particle physics. We obtain a new condition on the validity of the classical limit of these theories related to the non-trivial constraints that exist for vector field components. A violation of this consistency condition causes a singularity in the time derivative of the auxiliary component and could impact, for example, the field's cosmic history and superradiance around black holes. Such a condition is expected to exist generically in many other non-trivially constrained systems.

Scalar-tensor theories of gravity that embrace conformal coupling to the scalar curvature are the focal point of cosmology on discussions of inflation and late-time accelerating universe. Although there exists a stringent nucleo-synthesis constraint on conformal gravity, one can formulate how to evade this difficulty by modifying the standard particle theory action consistently with the principles of gauge invariant quantum field theory. It is shown that stronger gravity at early epochs of cosmological evolution than previously thought of is inevitable in a class of conformal gravity models. This enhances discovery potentials of primordial gravitational wave emission and primordial black hole formation. The strong gravity effect may be enormous if massive clumps are energetically dominated by cold dark matter made of inflaton field, and created black holes may become a major candidate of cold dark matter.

To undergo diffusive shock acceleration, electrons need to be pre-accelerated to increase their energies by several orders of magnitude, else their gyro-radii are smaller than the finite width of the shock. In oblique shocks, where the upstream magnetic field orientation is neither parallel or perpendicular to the shock normal, electrons can escape to the shock upstream, modifying the shock foot to a region called the electron foreshock. To determine the pre-acceleration in this region, we undertake PIC simulations of oblique shocks while varying the obliquity and in-plane angles. We show that while the proportion of reflected electrons is negligible for $\theta_{\rm Bn} = 74.3^\circ$, it increases to $R \sim 5\%$ for $\theta_{\rm Bn} = 30^\circ$, and that, via the electron acoustic instability, these electrons power electrostatic waves upstream with energy density proportional to $R^{0.6}$ and a wavelength $\approx 2 \lambda_{\rm{se}}$, where $\lambda_{\rm{se}}$ is the electron skin length. While the initial reflection mechanism is typically a combination of shock surfing acceleration and magnetic mirroring, we show that once the electrostatic waves have been generated upstream they themselves can increase the momenta of upstream electrons parallel to the magnetic field. In $\lesssim 1\%$ of cases, upstream electrons are prematurely turned away from the shock and never injected downstream. In contrast, a similar fraction are re-scattered back towards the shock after reflection, re-interact with the shock with energies much greater than thermal, and cross into the downstream.

We consider a cosmological model of non-minimal derivative coupling (NMDC) to gravity with holographic effect from Bekenstein-Hawking entropy using Hubble horizon IR cutoff. Holographic parameter $c$ is constant in a range, $0 \leq c < 1$. NMDC effect allows gravitational constant to be time-varying. Definition of holographic density include time-varying part of the gravitational constant. NMDC part reduces strength of gravitational constant for $\k > 0$ and opposite for $\k < 0$. The holographic part enhances gravitational strength. We use spectral index and tensor-to-scalar ratio to test the model against CMB constraint. Number of e-folding is chosen to be $N \geq 60$. Potentials, $V = V_0 \phi^n $ with $n = 2, 4$, and $V = V_0 \exp{(-\beta \phi)}$ are considered. Combined parametric plots of $\k$ and $\phi$ show that the allowed regions of the power spectrum index and of the tensor-to-scalar ratio are not overlapping. NMDC inflation is ruled out and the holographic NMDC inflation is also ruled out for $0 < c < 1$. NMDC significantly changes major anatomy of the dynamics, i.e. it gives new late-time attractor trajectories in acceleration regions. The holographic part clearly affects pattern of trajectories. However, for the holographic part to affect shape of the acceleration region, the NMDC field must be in presence. To constrain the model at late time, variation of gravitational constant is considered. Gravitational-wave standard sirens and supernovae data give a constraint, $\dot{G}/G|_{t_0} \lesssim 3\times10^{-12} \, \text{year}^{-1}$ \cite{Zhao:2018gwk} which, for this model, results in $ 10^{-12} \, \text{year}^{-1} \, \gtrsim \, {- \kappa} \dot{\phi}\ddot{\phi}/{M^2_{\p}}\,. $ Positive $\k$ is favored and greater $c^2$ results in lifting up lower bound of $\k$.

Interacting dark energy-dark matter models have been widely analyzed in the literature in an attempt to find traces of new physics beyond the usual cosmological ($\Lambda$CDM) models. Such a coupling between both dark components is usually introduced in a phenomenological way through a flux in the continuity equation. However, models with a Lagrangian formulation are also possible. A class of the latter assumes a conformal/disformal coupling that leads to a fifth force on the dark matter component, which consequently does not follow the same geodesics as the other (baryonic, radiation, and dark energy) matter sources. Here we analyze how the usual cosmological singularities of the standard matter frame are seen from the dark matter one, concluding that by choosing an appropriate coupling, dark matter observers will see no singularities but a non-beginning, non-ending universe. By considering two simple phenomenological models we show that such a type of coupling can fit observational data as well as the usual $\Lambda$CDM model.

Cosmological fine-tuning has traditionally been associated with the narrowness of the intervals in which the parameters of the physical models must be located to make life possible. A more thorough approach focuses on the probability of the interval, not on its size. Most attempts to measure the probability of the life-permitting interval for a given parameter rely on a Bayesian statistical approach for which the prior distribution of the parameter is uniform. However, the parameters in these models often take values in spaces of infinite size, so that a uniformity assumption is not possible. This is known as the normalization problem. This paper explains a framework to measure tuning that, among others, deals with normalization, assuming that the prior distribution belongs to a class of maximum entropy (maxent) distributions. By analyzing an upper bound of the tuning probability for this class of distributions the method solves the so-called weak anthropic principle, and offer a solution, at least in this context, to the well-known lack of invariance of maxent distributions. The implication of this approach is that, since all mathematical models need parameters, tuning is not only a question of natural science, but also a problem of mathematical modeling. Cosmological tuning is thus a particular instantiation of a more general scenario. Therefore, whenever a mathematical model is used to describe nature, not only in physics but in all of science, tuning is present. And the question of whether the tuning is fine or coarse for a given parameter -- if the interval in which the parameter is located has low or high probability, respectively -- depends crucially not only on the interval but also on the assumed class of prior distributions. Novel upper bounds for tuning probabilities are presented.

The faithful inclusion of the effects of bulk viscosity induced by the presence of chemical reactions is an important issue for simulations of core-collapse supernovae, binary neutron star mergers and neutron star oscillations, where particle abundances are locally pushed out of chemical equilibrium by rarefaction and compression of the fluid elements. In this work, we discuss three different approaches that can be used to implement bulk viscosity in general relativistic hydrodynamic simulations of neutron stars: the exact multi-component reacting fluid, and two M\"uller-Israel-Stewart theories, namely the second order Hiscock-Lindblom model and its linear limit, the Maxwell-Cattaneo model. After discussing the theory behind the three approaches, we specialize their dynamics equations to spherical symmetry in the radial gauge-polar slicing (i.e., Schwarzschild) coordinates. We also discuss a particular choice for the equation of state of the fluid and the associated neutrino emission rates, which are used in a companion paper for the numerical comparison of the three frameworks, and we obtain the effective sound speed for the Hiscock-Lindblom theory in the non-linear regime.

Out-of-equilibrium reactions between different particle species are the main process contributing to bulk viscosity in neutron stars. In this work, we numerically compare three different approaches to the modeling of bulk viscosity: the multi-component fluid with reacting particle species and two bulk stress formalism based on the M\"uller-Israel-Stewart theory, namely the Hiscock-Lindblom and the Maxwell-Cattaneo models, whose flux-conservative formulation in radial gauge-polar slicing coordinates and spherical symmetry is derived in a companion paper. In our knowledge, this is the first time that a neutron star is simulated with the complete Hiscock-Lindblom model of bulk viscosity. We find that the Hiscock-Lindblom and Maxwell-Cattaneo models are good approximations of the multi-component fluid for small perturbations and when the non-equilibrium equation of state of the fluid depends on only one independent particle fraction. For more than one independent particle fraction and for large perturbations, the bulk stress approximation is still valid but less accurate. In addition, we include the energy loss due to the luminosity of the reactions in the bulk stress formulation. We find that the energy loss due to bulk viscosity has a larger effect on the dynamics than the bulk stress or the variation in particle composition per se. The new one-dimensional, general-relativistic hydrodynamic code developed for this work, hydro-bulk-1D, is publicly available.