Papers for Wednesday, Mar 23 2022

We present a spectroscopic and photometric analysis of a sample of 416,288 galaxies from the Sloan Digital Sky Survey (SDSS) matched to mid-infrared (mid-IR) data from the Wide-Field Infrared Survey Explorer (WISE). By using a new spectroscopic fitting package, GELATO (Galaxy/AGN Emission Line Analysis TOol), we are able to retrieve emission line fluxes and uncertainties for SDSS spectra and robustly determine the presence of broad lines and outflowing components, enabling us to investigate WISE color space as a function of optical spectroscopic properties. In addition, we pursue SED template fitting to assess the relative AGN contribution and nuclear obscuration to compare to existing mid-IR selection criteria with WISE. We present a selection criterion in mid-IR color space to select Active Galactic Nuclei (AGNs) with a $\sim$80% accuracy and a completeness of $\sim$16%. This is the first mid-IR color selection defined by solely using the distribution of Type I and Type II optical spectroscopic AGNs in WISE mid-IR color space. Our selection is an improvement of $\sim$50% in the completeness of targeting spectroscopic AGNs with WISE down to an SDSS $r<17.77$ mag. In addition, our new criterion targets a less luminous population of AGNs, with on average lower [O III] luminosities by $\sim$30% ($>0.1$ dex) compared to typical WISE color-color selections. With upcoming large photometric surveys without corresponding spectroscopy, our method presents a way to select larger populations of AGNs at lower AGN luminosities and higher nuclear obscuration levels than traditional mid-IR color selections.

The Universe can inspire us to design communities that foster equity and inclusion.

The discrepancy between the value of the Hubble constant $H_0$ in the late, local universe and the one obtained from the Planck collaboration representing an all-sky value for the early universe reached the 5-$\sigma$ level. Approaches to alleviate the tension contain a wide range of ansatzes: increasing uncertainties in data acquisition, reducing biases in the astrophysical models that underly the probes, or taking into account observer-dependent variances in the parameters of the cosmological background model. Another way of interpreting the tension is to consider the density parameters in the Friedmann equations effective background fields which are probed at two different cosmic times, i.e. using the data of the early universe and the late, more complex universe for two different fits of a flat Friedmann-Lema\^itre-Robertson-Walker cosmology. The two sets of parameter values can differ and tensions arise if these background fits and perturbing biases are not consistently calibrated with respect to each other when they are linked to observables because these observables are all based on common standards. This consistent model-fitting calibration is lacking between the two $H_0$ values mentioned above, thus causing a tension. As shown here, this interpretation resolves the $H_0$ tension, if 15% of the matter-density parameter obtained from the fit to the cosmic microwave background, $\Omega_\mathrm{m}=0.315$, are assigned to decoupled perturbations yielding $\Omega_\mathrm{m}=0.267$ for the fit at redshifts of the supernova observations. Existing theoretical analyses and data evaluations which support this solution are given.

The Central Molecular Zone (CMZ) is a ring-like accumulation of molecular gas in the innermost few hundred parsecs of the Milky Way, generated by the inward transport of matter driven by the Galactic bar. The CMZ is the most extreme star-forming environment in the Galaxy. The unique combination of large-scale dynamics and extreme interstellar medium conditions, characterised by high densities, temperatures, pressures, turbulent motions, and strong magnetic fields, make the CMZ an ideal region for testing current star and planet formation theories. We review the recent observational and theoretical advances in the field, and combine these to draw a comprehensive, multi-scale and multi-physics picture of the cycle of matter and energy in the context of star formation in the closest galactic nucleus.

The engulfment of substellar bodies (SBs) such as brown dwarfs and planets has been invoked as a possible explanation for the presence of SBs orbiting subdwarfs and white dwarfs, rapidly rotating giants, and lithium-rich giants. We perform three-dimensional hydrodynamical simulations of the flow in the vicinity of an SB engulfed in a stellar envelope. We model the SB as a rigid body with a reflective boundary because it cannot accrete. This reflective boundary changes the flow morphology to resemble that of engulfed compact objects with outflows. We measure the drag coefficients for the ram pressure and gravitational drag forces acting on the SB, and use them to integrate its trajectory during engulfment. We find that SB engulfment can increase the stellar luminosity of a $1M_\odot$ star by up to a few orders of magnitude for timescales of up to a few thousand years when the star is $\approx10R_\odot$ and up to a few decades at the tip of the red giant branch. We find that no SBs can eject the envelope of a $1M_\odot$ star before it evolves to $\approx10R_\odot$. In contrast, SBs as small as $\approx10M_\text{Jup}$ can eject the envelope at the tip of the red giant branch, shrinking their orbits by several orders of magnitude in the process. The numerical framework we introduce here can be used to study the dynamics of planetary engulfment in a simplified setting that captures the physics of the flow at the scale of the SB.

The classification of the minor bodies of the Solar System based on observables has been continuously developed and iterated over the past 40 years. While prior iterations followed either the availability of large observational campaigns or new instrumental capabilities opening novel observational dimensions, we see the opportunity to improve primarily upon the established methodology. We develop an iteration of the asteroid taxonomy which allows for classification of partial and complete observations, i.e. visible, near-infrared, and visible-near-infrared spectrometry, and which re-introduces the visual albedo into the classification observables. The resulting class assignments are given probabilistically, enabling to quantify the uncertainty of a classification. We build the taxonomy based on 2983 observations of 2125 unique asteroids, representing an almost tenfold increase of sample size compared with the previous taxonomy. The asteroid classes are identified in a lower-dimensional representation of the observations using a Mixture of Common Factor Analyzers model. We identify 17 classes split into the three complexes C, M, and S, including the novel Z-class for extremely-red objects in the Main Belt. The visual albedo information resolves the spectral degeneracy of the X-complex and establishes the P-class as part of the C-complex. We present a classification tool which computes probabilistic class assignments within this taxonomic scheme from asteroid observations, intrinsically accounting for degeneracies between classes based on the observed wavelength region. The ability to classify partial observations and the re-introduction of the visual albedo into the classification provide a taxonomy which is well suited for the current and future datasets of asteroid observations, in particular provided by Gaia, the MITHNEOS-, and the SPHEREx surveys.

The detection of Ly$\alpha$ nebulae around $z\gtrsim 6$ quasars provides evidence for extended gas reservoirs around the first rapidly growing supermassive black holes. Observations of $z > 6$ quasars can be explained by cosmological models provided that the black holes by which they are powered evolve in rare, massive dark matter haloes. Whether these theoretical models also explain the observed extended Ly$\alpha$ emission remains an open question. We post-process a suite of cosmological, radiation-hydrodynamic simulations targeting a quasar host halo at $z>6$ with the Ly$\alpha$ radiative transfer code RASCAS. A combination of recombination radiation from photo-ionised hydrogen and emission from collisionally excited gas powers Ly$\alpha$ nebulae with a surface brightness profile in close agreement with observations. We also find that, even on its own, resonant scattering of the Ly$\alpha$ line associated to the quasar's broad line region can also generate Ly$\alpha$ emission on $\sim 100 \, \rm kpc$ scales, resulting in comparable agreement with observed surface brightness profiles. Even if powered by a broad quasar Ly$\alpha$ line, Ly$\alpha$ nebulae can have narrow line-widths $\lesssim 1000 \, \rm km \, s^{-1}$, consistent with observational constraints. Even if there is no quasar, we find that halo gas cooling produces a faint, extended Ly$\alpha$ glow. However, to light-up extended Ly$\alpha$ nebulae with properties in line with observations, our simulations unambiguously require quasar-powered outflows to clear out the galactic nucleus and allow the Ly$\alpha$ flux to escape and still remain resonant with halo gas. The close match between observations and simulations with quasar outflows suggests that AGN feedback already operates before $z \, = \, 6$ and confirms that high-$z$ quasars reside in massive haloes tracing overdensities.

Photochemistry is expected to change the chemical composition of the upper atmospheres of irradiated exoplanets through the dissociation of species, such as methane and ammonia, and the association of others, such as hydrogen cyanide. Although primarily the high altitude day side should be affected by photochemistry, it is still unclear how dynamical processes transport photochemical species throughout the atmosphere, and how these chemical disequilibrium effects scale with different parameters. In this work we investigate the influence of photochemistry in a two-dimensional context, by synthesizing a grid of photochemical models across a large range of temperatures. We find that photochemistry can strongly change the atmospheric composition, even up to depths of several bar in cool exoplanets. We further identify a sweet spot for the photochemical production of hydrogen cyanide and acetylene, two important haze precursors, between effective temperatures of 800 and 1400 K. The night sides of most cool planets (effective temperature < 1800 K) are shown to host photochemistry products, transported from the day side by horizontal advection. Synthetic transmission spectra are only marginally affected by photochemistry, but we suggest that observational studies probing higher altitudes, such as high-resolution spectroscopy, take photochemistry into account.

[Abridged] We present a novel method that takes advantage of the unique capabilities of the Gaia satellite to obtain large and reliable samples of dual or lensed AGN candidates with sub-arcsec separations, by selecting targets showing multiple peaks in the Gaia scans. With this method, we have identified 260 multiple system candidates, 26 of which have archival Hubble Space Telescope (HST) images. All these images show two or more point-sources at sub-arcsec separations, thus confirming that the proposed technique can be extremely efficient at small separations. Thirteen out of these 26 systems are known gravitational lenses, and two have been previously classified as dual AGNs. One of these two systems has an HST/STIS spatially resolved spectrum showing two distinct AGNs at z=2.95 with 3.6 kpc (0.46 arcsec) separation, confirming the previous identification as dual AGN. Dedicated adaptive-optics (AO) assisted high-resolution imaging at the Large Binocular Telescope (LBT) of 5 systems also detected multiple components in all the targets, further confirming the selection method. The measured separations are between 0.33 and 0.66 arcsec (2.7 - 5.6 kpc at the distance of the primary AGN). The nature of the unclassified sources is uncertain. Several tests demonstrate that contamination from foreground stars can only account for ~30% of these multiple systems, and the analysis of the observed colours suggests that most of the selected systems are expected to be dual/lensed AGNs with physical separations between 3 and 5 kpc. Our results show that this is a very efficient technique to select compact dual/lensed AGN systems with sub-arcsec separations, complementing other methods. This method samples separation down to 2 kpc at z>1 and thus allows us to probe the physical processes driving the inspiralling of the pairs of SMBH inside a single galaxy.

Recent studies using New Horizons LORRI images have returned the most precise measurement of the cosmic optical background to date, yielding a flux that exceeds that expected from deep galaxy counts by roughly a factor of two. We investigate whether this excess, detected at $\sim 4\sigma$ significance, is due to dark matter that decays to a monoenergetic photon with a rest-frame energy in the range $0.5-10$ eV. We derive the spectral energy distribution from such decays and the contribution to the flux measured by LORRI. The parameter space that explains the measured excess with decays to photons with energies $E\gtrsim 4$ eV is unconstrained to date. If the excess arises from dark-matter decay to a photon line, there will be a significant signal in forthcoming line-intensity mapping measurements. Moreover, the ultraviolet instrument aboard New Horizons (which will have better sensitivity and probe a different range of the spectrum) and future studies of very high-energy $\gamma$-ray attenuation will also test this hypothesis and expand the search for dark matter to a wider range of frequencies.

The tidal interactions between a planet and moon can provide insight into the properties of the host planet. The recent exomoon candidates Kepler-1708 b-i and Kepler-1625 b-i are Neptune-sized satellites orbiting Jupiter-like planets and provide an opportunity to apply such methods. We show that if the tidal migration time is roughly equal to the age of these systems, then the tidal dissipation factor Q for the planets Kepler-1708 b and Kepler-1625 b have values of ~$3\times10^5-3\times10^6$ and ~$1.5\times10^5-4\times10^5$, respectively. In each case, these are consistent with estimates for gas giant planets. Even though some work suggests an especially large semimajor axis for Kepler-1625 b-i, we find that this would imply a surprisingly low Q~2000 for a gas giant unless the moon formed at essentially its current position. More detailed predictions for the moons' initial semimajor axis could provide even better constraints on Q, and we discuss the formation scenarios for a moon in this context. Similar arguments can be used as more exomoons are discovered in the future to constrain exoplanet interior properties. This could be especially useful for exoplanets near the sub-Neptune/super-Earth radius gap where the planet structure is uncertain.

In order to assess the multiplicity statistics of stars across spectral types and populations in a volume-limited sample, we censused nearby stars for companions with Robo-AO. We report on observations of 1157 stars of all spectral types within 25 pc with decl. $>-13^{\circ}$ searching for tight companions. We detected 154 companion candidates with separations ranging from $\sim$0.15$''$ to 4.0$''$ and magnitude differences up to $\Delta$m$_{\textit{i'}}\le$7 using the robotic adaptive optics instrument Robo-AO. We confirmed physical association from Gaia EDR3 astrometry for 53 of the companion candidates, 99 remain to be confirmed, and 2 were ruled out as background objects. We complemented the high-resolution imaging companion search with a search for co-moving objects with separations out to 10,000 AU in Gaia EDR3, which resulted in an additional 147 companions registered. Of the 301 total companions reported in this study, 49 of them are new discoveries. Out of the 191 stars with significant acceleration measurements in the Hipparcos-Gaia catalog of accelerations, we detect companions around 115 of them, with the significance of the acceleration increasing as the companion separation decreases. From this survey, we report the following multiplicity fractions (compared to literature values): 40.9%$\pm$3.0% (44%) for FGK stars and 28.2%$\pm$2.3% (27%) for M stars, as well as higher-order fractions of 5.5%$\pm$1.1% (11%) and 3.9%$\pm$0.9% (5%) for FGK stars and M type stars, respectively.

We report the first direct measurement of the helium isotope ratio, 3He/4He, outside of the Local Interstellar Cloud, as part of science verification observations with the upgraded CRyogenic InfraRed Echelle Spectrograph (CRIRES). Our determination of 3He/4He is based on metastable HeI* absorption along the line-of-sight towards Tet02 Ori A in the Orion Nebula. We measure a value 3He/4He=(1.77+/-0.13)x10^{-4}, which is just ~40 per cent above the primordial relative abundance of these isotopes, assuming the Standard Model of particle physics and cosmology, (3He/4He)_p = (1.257+/-0.017)x10^-4. We calculate a suite of galactic chemical evolution simulations to study the Galactic build up of these isotopes, using the yields from Limongi & Chieffi (2018) for stars in the mass range M=8-100 M_sun and Lagarde (2011,2012) for M=0.8-8 M_sun. We find that these simulations simultaneously reproduce the Orion and protosolar 3He/4He values if the calculations are initialized with a primordial ratio (3He/4He)_p=(1.043+/-0.089)x10^-4. Even though the quoted error does not include the model uncertainty, this determination agrees with the Standard Model value to within ~2sigma. We also use the present-day Galactic abundance of deuterium (D/H), helium (He/H), and 3He/4He to infer an empirical limit on the primordial 3He abundance, (3He/H)_p < (1.09+/-0.18)x10^-5, which also agrees with the Standard Model value. We point out that it is becoming increasingly difficult to explain the discrepant primordial 7Li/H abundance with non-standard physics, without breaking the remarkable simultaneous agreement of three primordial element ratios (D/H, 4He/H, and 3He/4He) with the Standard Model values.

Variable accretion in young stellar objects reveals itself photometrically and spectroscopically over a continuum of timescales and amplitudes. Most dramatic are the large outbursts (e.g., FU Ori, V1647 Ori, and EX Lup type events), but more frequent are the less coherent, smaller burst-like variations in accretion rate. Improving our understanding of time-variable accretion directly addresses the fundamental question of how stars gain their masses. We review variability phenomena, as characterized from observations across the wavelength spectrum, and how those observations probe underlying physical conditions. The diversity of observed lightcurves and spectra at optical and infrared wavelengths defies a simple classification of outbursts and bursts into well-defined categories. Mid-infrared and sub-millimeter wavelengths are sensitive to lower-temperature phenomena, and it is currently unclear if observed flux variations probe similar or distinct physics relative to the shorter wavelengths. We highlight unresolved issues and emphasize the value of spectroscopy, multiwavelength studies, and ultimately patience in using variable accretion to understand stellar mass assembly.

The maximum mass of a star that can produce a white dwarf (WD) is an important astrophysical quantity. One of the best approaches to establishing this limit is to search for WDs in young star clusters in which only massive stars have had time to evolve and where the mass of the progenitor can be established from the cooling time of the WD together with the age of the cluster. Searches in young Milky Way clusters have not thus far yielded WD members more massive than about 1.1Msun, well below the Chandrasekhar mass of 1.38Msun, nor progenitors with masses in excess of about 6Msun. However, the hunt for potentially massive WDs that escaped their cluster environs is yielding interesting candidates. To expand the cluster sample further, we used HST to survey four young and massive star clusters in the Magellanic Clouds for bright WDs that could have evolved from stars as massive as 10Msun. We located five potential WD candidates in the oldest of the four clusters examined, the first extragalactic single WDs thus far discovered. As these hot WDs are very faint at optical wavelengths, final confirmation will likely have to await spectroscopy with 30-metre class telescopes.

Globular clusters display an anti-correlation between the fraction of the first generation of stars ($N({\rm G1})/N({\rm tot})$) and the slope of the present-day mass function of the clusters ($\alpha_{pd}$), which is particularly significant for massive clusters. In the framework of the binary-mediated collision scenario for the formation of the second generation stars in globular clusters, we test the effect of a varying initial stellar mass function (IMF) of the G1 stars on the $(N({\rm G1})/N({\rm tot}))-\alpha_{pd}$ anti-correlation. We use a simple collision model which has only two inputs, the IMF of G1 stars and the fraction of G1 stars that coalesce to form G2 stars. We show that a variable efficiency of the collision process is necessary in order to explain the $(N({\rm G1})/N({\rm tot}))-\alpha_{pd}$ anti-correlation, however, the scatter in the anti-correlation can only be explained by variations in the IMF, and in particular by variations of the slope in the mass interval $\approx$ (0.1-0.5) M$_{\odot}$. Our results indicate that in order to explain the scatter in the $(N({\rm G1})/N({\rm tot}))-\alpha_{pd}$ relation, it is necessary to invoke variations of the slope in this mass range between $\approx -0.9$ and $\approx -1.9$. Interpreted in terms of a Kroupa-like broken power law, this translates into variations of the mean mass between $\approx 0.2$ and $0.55$ M$_{\odot}$. This level of variations is consistent with what is observed for young stellar clusters in the Milky Way and may reflect variations in the physical conditions of the globular clusters progenitor clouds at the time the G1 population has formed or to the occurrence of collisions between protostellar embryos before stars settle on the Main Sequence.

The minimum - maximum method, belonging to the precursor class of the solar activity forecasting methods, is based on a linear relationship between relative sunspot number in the minimum and maximum epochs of solar cycles. In the present analysis we apply a modified version of this method using data not only from the minimum year, but also from a couple of years before and after the minimum. The revised 13-month smoothed monthly total sunspot number data set from SILSO/SIDC is used. Using data for solar cycle nos. 1-24 the largest correlation coefficient (CC) is obtained when correlating activity level 3 years before solar cycle minimum with the subsequent maximum (CC = 0.82), independent of inclusion or exclusion of the solar cycle no. 19. For the next solar maximum of the cycle no. 25 we predict: Rmax = 121 +- 33. Our results indicate that the next solar maximum (of the cycle no. 25) will be of the similar amplitude as the previous one, or even something lower. This is in accordance with the general middle-term lowering of the solar activity after the secular maximum in the 20th century and consistent with the Gleissberg period of the solar activity. The reliability of the 3 years before the minimum predictor is experimentally justified by the largest correlation coefficient and verified with the Student t-test. It is satisfactorily explained with the two empirical well-known findings: the extended solar cycle and the Waldmeier effect. Finally, we successfully reproduced the maxima of the last four solar cycles, nos. 21-25, using the 3 years before the minimum method.

Modelling ultra-low-noise far-infrared grating spectrometers has become crucial for the next generation of far-infrared space observatories. Conventional techniques are awkward to apply because of the partially coherent form of the incident spectral field, and the few-mode response of the optics and detectors. We present a modal technique for modelling the behaviour of spectrometers, which allows for the propagation and detection of partially coherent fields, and the inclusion of straylight radiated by warm internal surfaces. We illustrate the technique by modelling the behaviour of the Long Wavelength Band of the proposed SAFARI instrument on the well-studied SPICA mission.

The CMB lensing signal from cosmic voids and superclusters probes the growth of structure in the low-redshift cosmic web. In this analysis, we cross-correlated the Planck CMB lensing map with voids detected in the Dark Energy Survey Year 3 (Y3) data set ($\sim$5,000 deg$^{2}$), extending previous measurements using Y1 catalogues ($\sim$1,300 deg$^{2}$). Given the increased statistical power compared to Y1 data, we report a $6.6\sigma$ detection of negative CMB convergence ($\kappa$) imprints using approximately 3,600 voids detected from a redMaGiC luminous red galaxy sample. However, the measured signal is lower than expected from the MICE N-body simulation that is based on the $\Lambda$CDM model (parameters $\Omega_{\rm m} = 0.25$, $\sigma_8 = 0.8$). The discrepancy is associated mostly with the void centre region. Considering the full void lensing profile, we fit an amplitude $A_{\kappa}=\kappa_{\rm DES}/\kappa_{\rm MICE}$ to a simulation-based template with fixed shape and found a moderate $2\sigma$ deviation in the signal with $A_{\kappa}\approx0.79\pm0.12$. We also examined the WebSky simulation that is based on a Planck 2018 $\Lambda$CDM cosmology, but the results were even less consistent given the slightly higher matter density fluctuations than in MICE. We then identified superclusters in the DES and the MICE catalogues, and detected their imprints at the $8.4\sigma$ level; again with a lower-than-expected $A_{\kappa}=0.84\pm0.10$ amplitude. The combination of voids and superclusters yields a $10.3\sigma$ detection with an $A_{\kappa}=0.82\pm0.08$ constraint on the CMB lensing amplitude, thus the overall signal is $2.3\sigma$ weaker than expected from MICE.

We present an updated analysis of systematics in the Atacama Large Millimeter/submillimeter Array (ALMA) proposal ranks from Carpenter (2020) to include the last two ALMA cycles, when significant changes were introduced in the proposal review process. In Cycle 7, the investigator list on the proposal cover sheet was randomized such that the reviewers were aware of the overall proposal team but did not know the identity of the principal investigator (PI). In Cycle 8, ALMA adopted distributed peer review for most proposals and implemented dual-anonymous review for all proposals, in which the identity of the proposal team was not revealed to the reviewers. The most significant change in the systematics in Cycles 7 and 8 compared to previous cycles is related to the experience of PIs in submitting ALMA proposals. PIs that submit a proposal every cycle tend to have ranks that are consistent with average in Cycles 7 and 8 whereas previously they had the best overall ranks. Also, PIs who submitted a proposal for the second time show improved ranks over previous cycles. These results suggest some biases related to the relative prominence of the PI have been present in the ALMA review process. Systematics related to regional affiliation remain largely unchanged in that PIs from Chile, East Asia, and non-ALMA regions tend to have poorer overall ranks than PIs from Europe and North America. The systematics of how one region ranks proposals from another region are also investigated. No significant differences in the overall ranks based on gender of the PI are observed.

China Space Station Telescope (CSST) is a forthcoming powerful Stage IV space-based optical survey equipment. It is expected to explore a number of important cosmological problems in extremely high precision. In this work, we focus on investigating the constraints on neutrino mass and other cosmological parameters under the model of cold dark matter with a constant equation of state of dark energy ($w$CDM), using the mock data from the CSST photometric galaxy clustering and cosmic shear surveys (i.e. 3$\times$2pt). The systematics from galaxy bias, photometric redshift uncertainties, intrinsic alignment, shear calibration, baryonic feedback, non-linear, and instrumental effects are also included in the analysis. We generate the mock data based on the COSMOS catalog considering the instrumental and observational effects of the CSST, and make use of the Markov Chain Monte Carlo (MCMC) method to perform the constraints. Comparing to the results from current similar measurements, we find that CSST 3$\times$2pt surveys can improve the constraints on the cosmological parameters by one order of magnitude at least. We can obtain an upper limit for the sum of neutrino mass $\Sigma m_{\nu} \lesssim 0.36$ (0.56) eV at 68\% (95\%) confidence level, and $\Sigma m_{\nu} \lesssim 0.23$ (0.29) eV at 68\% (95\%) confidence level if ignore the baryonic effect, which is comparable to the {\it Planck} results and much better than the current photometric surveys. This indicates that the CSST photometric surveys can provide stringent constraints on the neutrino mass and other cosmological parameters, and the results also can be further improved by including data from other kinds of CSST cosmological surveys.

Polarized Galactic synchrotron emission is an undesirable foreground for cosmic microwave background (CMB) experiments observing at frequencies $< 150$ GHz. We perform a combined analysis of observational data at 1.4, 2.3, 23, 30 and 33 GHz to quantify the spatial variation of the polarized synchrotron spectral index, $\beta^{pol}$, on $\sim3.5^\circ$ scales. We compare results from different data combinations to address limitations and inconsistencies present in these public data, and form a composite map of $\beta^{pol}$. Data quality masking leaves 44% sky coverage (75% for $|b|> 45^\circ$). Generally $-3.2 < \beta^{pol} \lesssim -3$ in the inner Galactic plane and spurs, but the Fan Region in the outer Galaxy has a flatter index. We find a clear spectral index steepening with increasing latitude south of the Galactic plane with $\Delta \beta^{pol}=0.4$, and a smaller steepening of $0.25$ in the north. Near the south Galactic pole the polarized synchrotron spectral index is $\beta^{pol} \approx -3.4$. Longitudinal spectral index variations of $\Delta \beta^{pol} \sim 0.1$ about the latitudinal mean are also detected. Within the BICEP2/Keck survey footprint, we find consistency with a constant value, $\beta^{pol} = -3.25 \pm 0.04$ (statistical) $\pm 0.02$ (systematic). We compute a map of the frequency at which synchrotron and thermal dust emission contribute equally to the total polarized foreground. The limitations and inconsistencies among datasets encountered in this work make clear the value of additional independent surveys at multiple frequencies, especially between $10-20$ GHz, provided these surveys have sufficient sensitivity and control of instrumental systematic errors.

We searched for correlations between the number of satellites and fundamental galactic properties for the Milky Way-like host galaxies in order to better understand their diverse satellite populations. We specifically aim to understand why galaxies that are very similar in stellar mass content, star formation rate, and local environment have very different numbers of satellites. Satellites of Galactic Analogs (SAGA) spectroscopic survey has completed spectroscopic observations of 36 Milky Way-like galaxies within their virial radii down to the luminosity of Leo I dwarf galaxy. All the available galactic properties of SAGA galaxies from the literature along with measured stellar masses were correlated with the number of satellites and no significant correlation was found. However, when we considered the "expected" number of satellites based on the correlation between the baryonic bulge-to-total ratio and the number of satellites confirmed for several nearby galaxies then strong correlations emerge between this number and (1) the mass of the bulge, and (2) the total specific angular momentum. The first correlation is positive, implying that galaxies with more massive bulges have more satellites, as already confirmed. The second correlation with the angular momentum is negative, meaning that, the smaller the angular momentum, the greater the number of expected satellites. This would imply that either satellites cannot form if galaxy angular momentum is too high, or that satellites form inside-out, so that angular momentum is being transferred to the outer parts of the galaxies. However, deeper spectroscopic observations are needed to confirm these findings, because they rely on the expected rather than detected number of satellites. There was a luminosity limit to the SAGA survey equivalent to the luminosity of Leo I dwarf satellite of the Milky Way galaxy.

There exists an inevitable scatter in intrinsic luminosity of Gamma Ray Bursts(GRBs). If there is relativistic beaming in the source, viewing angle variation necessarily introduces variation in the intrinsic luminosity function (ILF). Scatter in the ILF can cause a selection bias where distant sources that are detected have a larger median luminosity than those detected close by. Median luminosity divides any given population into equal halves. When the functional form of a distribution is unknown, it can be a more robust diagnostic than any that use trial functional forms. In this work we employ a statistical test based on median luminosity and apply it to test a class of models for GRBs. We assume that the GRB jet has a finite opening angle and that the orientation of the GRB jet is random relative to the observer.We calculate $L_{median}$ as a function of redshift by simulating GRBs empirically, theoretically and use the luminosity vs redshift {\it Swift} data in order to compare the theoretical results with the observed ones. The method accounts for the fact that at some redshifts there may be some GRBs that go undetected. We find that $L_{median}$ is extremely insensitive to the on-axis (i.e. maximal) luminosity of the jet.

Characterising spectral variability of radio sources is a technique that offers the ability to determine the astrophysics of the intervening media, source structure, emission and absorption processes. We present broadband (0.072--10 GHz) spectral variability of 15 peaked-spectrum (PS) sources with the Australia Telescope Compact Array (ATCA) and the Murchison Widefield Array (MWA). These 15 PS sources were observed quasi-contemporaneously with ATCA and the MWA four to six times during 2020 with approximately a monthly cadence. Variability was not detected at 1--10GHz frequencies but 13 of the 15 targets show significant variability with the MWA at megahertz frequencies. We conclude the majority of variability seen at megahertz frequencies is due to refractive interstellar scintillation of a compact component ~25 mas across. We also identify four PS sources that show a change in their spectral shape at megahertz frequencies. Three of these sources are consistent with a variable optical depth from an inhomogeneous free-free absorbing cloud around the source. One PS source with a variable spectral shape at megahertz frequencies is consistent with an ejection travelling along the jet. We present spectral variability as a method for determining the physical origins of observed variability and for providing further evidence to support absorption models for PS sources where spectral modelling alone is insufficient.

Prediction of the solar cycle is challenging but essential because it drives space weather. Several predictions with varying amplitudes of the ongoing Cycle~25 have been made. We show that an aspect of the Waldmeier effect, i.e., a strong positive correlation between the rise rate and the amplitude of the cycle, has a physical link with the build-up of the previous cycle's polar field after its reversal. We find that the rise rate of the polar field is highly correlated with the rise rate and the amplitude of the next solar cycle. Thus, the prediction of the amplitude of the solar cycle can be made just a few years after the reversal of the previous cycle's polar field, thereby extending the scope of the solar cycle prediction to much earlier than the usual time. Our prediction of Cycle 25 based on the rise rate of the previous polar field is $137\pm 23$, which is quite close to the prediction $138\pm 26$ based on the WE2 computed from the available 2 years sunspot data of the ongoing cycle.

The Crab Nebula emits bright non-thermal radiation from radio to the most energetic photons. The underlying physical model of a relativistic wind from the pulsar terminating in a hydrodynamic standing shock remains unchanged since the early 1970s. In this model, an increase of the toroidal magnetic field downstream from the shock is expected. We introduce a detailed radiative model to calculate non-thermal synchrotron and inverse Compton as well as thermal dust emission self-consistently to compare quantitatively with observational data. Special care is given to the radial dependence of electron and seed field density. The radiative model is used to estimate the parameters of electrons and dust in the nebula. A combined fit based upon a $\chi^2$ minimisation reproduces successfully the complete data set used. For the best-fitting model, the energy density of the magnetic field dominates over the particle energy density up to a distance of $\approx 1.3~r_s$ ($r_s$: distance of the termination shock from the pulsar). The very high energy (VHE: $E>100$ GeV) gamma-ray spectra set the strongest constraints on the radial dependence of the magnetic field: $B(r)=(264\pm9)~\mu\mathrm{G} (r/r_s)^{-0.51\pm0.03}$. The reconstructed magnetic field and its radial dependence indicates a ratio of Poynting to kinetic energy flux $\sigma\approx 0.1$ at the termination shock, $\approx 30$ times larger than estimated up to now. Consequently, the confinement of the nebula would require additional mechanisms to slow down the flow through e.g. excitation of small-scale turbulence with possible dissipation of magnetic field.

We present a probabilistic model built for the optical cross-match between the SRG/eROSITA X-ray sources and photometric data from the DESI Legacy Imaging Surveys. The model relies both on positional and photometric information on optical objects nearby X-ray sources and allows performing selection with precision and recall $\approx94$% (for $F_{\rm X,0.5-2}>10^{-14}$ erg/s/cm$^2$). With this model, we calibrated positional error of the SRG/eROSITA sources detected in the Lockman Hole: $\sigma_{\rm corr} = 0.87\sqrt{ \sigma_{\rm det}^{2.53} + 1.12^2}$. The model will become a part of the SRGz system for data analysis of the X-ray data obtained from the all-sky SRG/eROSITA survey.

Using light curves with $\sim$3 min cadence and a duration of 2 hrs made using the OmegaWhite survey, we present the results of a search for short-period variable stars in the field of 20 open clusters. We identified 92 variable stars in these fields. Using a range of cluster member catalogues and Gaia EDR3 data, we have determined that 10 are cluster members and 2 more are probable members. Based on their position on the Gaia HRD and their photometric periods, we find that most of these are $\delta$ Sct stars. We obtained low-resolution optical spectroscopy of some of these cluster members and field stars. We discuss the cluster variable stars in the context of $\delta$ Sct stars in other open clusters.

Galaxy formation and evolution simulations are essential tools to probe poorly known astrophysics processes, but particular care is needed to compare simulations with galaxy observations, as observed data need to be modelled as well. We present a method to generate mock galaxies from the hydro-dynamical IllustrisTNG simulations which are suited to compare with integral field spectroscopic observation of galaxies from the SDSS-IV/MaNGA survey. Firstly, we include the same instrumental effects and procedures as adopted in the acquisition and analysis of real data. Furthermore, we generate the galaxy spectra from the simulations using new stellar population models based on the MaNGA stellar library (MaStar). In this way, our mock data cubes have the same spatial sampling, cover the same wavelength range (3600-10300 \r{A}), and share the same spectral resolution ($R\approx1800$) and flux calibration of real MaNGA galaxy spectra. In this first paper, we demonstrate the method over an early-type and a late-type simulated galaxy from TNG50. We analyse the correspondent mock MaNGA-like data cubes with the same full spectral fitting code, \textsc{FIREFLY}, which was used for the observed spectra. We find that the intrinsic and recovered age and metallicity gradients are consistent within 1${\sigma}$, with residuals over all tassels consistent with $0$ at the 68$\%$ confidence level. We also perform the challenging test at comparing intrinsic and recovered star formation histories, finding a close resemblance between input and output. In follow-up papers, we will present a full simulated MaNGA-like catalogue ($\approx10,000$ galaxies) with a comprehensive comparison of TNG50 simulations to MaNGA observational results.

Young stellar objects in their pre-main sequence phase are characterized by irregular changes in brightness, generally attributed to an increase of the mass accretion rate due to various kind of instabilities occurring in the circumstellar disk. In the era of large surveys aimed to monitor the sky, we present a pipeline to detect irregular bursts, in particular EXors-like ( EX Lupi type eruptive variables), in the light curves. The procedure follows a heuristic approach and is tested against the light curves already collected for a few objects presently recognized as bona fide or candidate EXors.

Nanoflares in quiet-Sun regions during solar cycle 24 are studied with the best available plasma diagnostics to derive their energy distribution and contribution to coronal heating during different levels of solar activity. Extreme ultraviolet (EUV) filters of the Atmospheric Imaging Assembly (AIA) onboard the Solar Dynamics Observatory (SDO) are used. We analyze 30 AIA/SDO image series between 2011 and 2018, each covering a $400 \times 400$ arcsec quiet-Sun field-of-view over two hours with a 12-second cadence. Differential emission measure (DEM) analysis is used to derive the emission measure (EM) and temperature evolution for each pixel. We detect nanoflares as EM-enhancements using a threshold-based algorithm and derive their thermal energy from the DEM observations. Nanoflare energy distributions follow power-laws that show slight variations in steepness ($\alpha =$ 2.02 to 2.47) but no correlation to the solar activity level. The combined nanoflare distribution of all data sets covers five orders of magnitude in event energies ($10^{24} \mathrm{~to~} 10^{29} \mathrm{~erg}$) with a power-law index $\alpha=2.28 \pm0.03$. The derived mean energy flux of $(3.7\pm 1.6)\times 10^4\mathrm{~erg~cm^{-2}~s^{-1}}$ is one order of magnitude smaller than the coronal heating requirement. We find no correlation between the derived energy flux and solar activity. Analysis of the spatial distribution reveals clusters of high energy flux (up to $3\times 10^5 \mathrm{~erg~cm^{-2}~s^{-1}}$) surrounded by extended regions with lower activity. Comparisons with magnetograms from the Helioseismic and Magnetic Imager (HMI) demonstrate that high-activity clusters are located preferentially in the magnetic network and above regions of enhanced magnetic flux density.

We present a new reconstruction of the distribution of atomic hydrogen in the inner Galaxy that is based on explicit radiation-transport modelling of line and continuum emission and a gas-flow model in the barred Galaxy that provides distance resolution for lines of sight toward the Galactic Center. The main benefits of the new gas model are, a), the ability to reproduce the negative line signals seen with the H$I$4PI survey and, b), the accounting for gas that primarily manifests itself through absorption. We apply the new model of Galactic atomic hydrogen to an analysis of the diffuse gamma-ray emission from the inner Galaxy, for which an excess at a few GeV was reported that may be related to dark matter. We find with high significance an improved fit to the diffuse gamma-ray emission observed with the Fermi-LAT, if our new H$I$ model is used to estimate the cosmic-ray induced diffuse gamma-ray emission. The fit still requires a nuclear bulge at high significance. Once this is included there is no evidence for a dark-matter signal, be it cuspy or cored. But an additional so-called boxy bulge is still favoured by the data. This finding is robust under the variation of various parameters, for example the excitation temperature of atomic hydrogen, and a number of tests for systematic issues.

We have studied the spectral time variations of candidate luminous blue variable stars (cLBV) in two low-metallicity blue compact dwarf galaxies, DDO 68 and PHL 293B. The LBV in DDO 68, located in HII region #3, shows an outburst, with an increase of more than 1000 times in Halpha luminosity during the period 2008-2010. The broad emission of the HI and HeI lines display a P Cygni profile, with a relatively constant terminal velocity of ~800 km/s, reaching a maximum luminosity L(Halpha) of ~2x10^38 erg/s, with a FWHM of ~1000-1200 km/s. On the other hand, since the discovery of a cLBV in 2001 in PHL 293B, the fluxes of the broad components and the broad-to-narrow flux ratios of the HI and HeI emission lines in this galaxy have remained nearly constant over 16 years, with small variations. The luminosity of the broad Halpha component varies between ~2x10^38 erg/s and ~10^39 erg/s, with the FWHM varying in the range ~500-1500 km/s. Unusually persistent P Cygni features are clearly visible until the end of 2020 despite a decrease of the broad-to-narrow flux ratio in the most recent years. A terminal velocity of ~800 km/s is measured from the P Cygni profile, similar to the one in DDO 68, although the latter is 3.7 more metal-deficient than PHL 293B. The relative constancy of the broad Halpha luminosity in PHL 293B suggests that it is due to a long-lived stellar transient of type LBV/SN IIn.

GJ 229B was the first T-dwarf to be discovered in 1995, and its spectrum has been more thoroughly observed than most other brown dwarfs. Yet a full spectroscopic analysis of its atmosphere has not been done with modern techniques. This spectrum has several peculiar features, and recent dynamical estimates of GJ 229B's mass and orbit have disagreed widely, both of which warrant closer investigation. With a separation of tens of AU from its host star, GJ 229B falls near the border of the planet and stellar population formation regimes, so its atmosphere could provide clues to formation processes for intermediate objects of this type. In an effort to resolve these questions, we performed retrievals on published spectra of GJ 229B over a wide range of wavelengths (0.5-5.1 ${\rm \mu m}$) using the open-source APOLLO code. Based on these retrievals, we present a more precise mass estimate of $41.6\pm3.3\, M_J$ and an effective temperature estimate of $869_{-7}^{+5}$ K, which are more consistent with evolutionary models for brown dwarfs and suggest an older age for the system of $>$1.0 Gyr. We also present retrieved molecular abundances for the atmosphere, including replicating the previously-observed high CO abundance, and discuss their implications for the formation and evolution of this object. This retrieval effort will give us insight into how to study other brown dwarfs and directly-imaged planets, including those observed with JWST and other next-generation telescopes.

The problem of orbit flips caused by eccentric von Zeipel-Lidov-Kozai effects is systematically investigated by means of three approaches, including Poincar\'e sections, dynamical system theory (periodic orbits and invariant manifolds) and perturbation treatments. Poincar\'e sections show that orbit flips are due to the existence of islands of libration centered at inclination of $90^{\circ}$, dynamical system theory shows that orbit flips are due to the existence of polar periodic orbits and invariant manifolds, and perturbative treatments indicate that orbit flips are due to the libration of a certain critical argument. Using these approaches, the boundaries of flipping regions in the entire parameter space are produced and they are in excellent agreement with each other. Through analysis, the essence of flipping orbits is reached: (a) flipping orbits are a kind of quasi-periodic trajectories around polar periodic orbits and invariant manifolds at the same level of Hamiltonian provide boundaries of flipping regions, and (b) flipping orbits are a kind of resonant trajectories and resonant width measures the size of flipping regions.

A stellar occultation occurs when a Solar System object passes in front of a star for an observer. This technique allows the determination of sizes and shapes of the occulting body with kilometer precision. Also, this technique constrains the occulting body's positions, albedos, densities, etc. In the context of the Galilean moons, these events can provide their best ground-based astrometry, with uncertainties in the order of 1 mas ($\sim$ 3 km at Jupiter's distance during opposition). We organized campaigns and successfully observed a stellar occultation by Io (JI) in 2021, one by Ganymede (JIII) in 2020, and one by Europa (JII) in 2019, with stations in North and South America. Also, we re-analyzed two previously published events, one by Europa in 2016 and another by Ganymede in 2017. Then, we fit the known 3D shape of the occulting satellite and determine its center of figure. That resulted in astrometric positions with uncertainties in the milliarcsecond level. The positions obtained from these stellar occultations can be used together with dynamical models to ensure highly accurate orbits of the Galilean moons. These orbits can help plan future space probes aiming at the Jovian system, such as JUICE by ESA and Europa Clipper by NASA, and allow more efficient planning of flyby maneuvers.

From estimates of the near-surface heat capacity of planets it is shown that the thermal time scale is larger than the orbital period in the presence of a global ocean that is well-mixed to a depth of 100 m, or of an atmosphere with a pressure of several tens of bars. As a consequence, the temperature fluctuations of such planets on eccentric orbits are damped. The average temperature should be calculated by taking the temporal mean of the irradiation over an orbit, which increases with $1/\sqrt{1-e^2}$. This conclusion is independent of the orbital distance and valid for Sun-like stars; the damping is even stronger for low-mass main sequence hosts.

Ru{\dj}er Bo\v{s}kovi\'c developed methods for determination of solar rotation elements: the solar equator inclination i, the longitude of the node {\Omega} and the period of solar rotation. In his last work Opera pertinentia ad opticam et astronomiam, published in 1785, in the chapter Opuscule II he described his methods, the formulae with figure descriptions and an example for calculation of the solar rotation elements with detailed numerical explanation using his own observations performed in September 1777. The original numerical procedure was performed using logarithmic formulae. In present work we give a description of the original results of R. Bo\v{s}kovi\'c and compare them with our recalculated values.

Images of the Sun at millimeter wavelengths obtained by ALMA show a significant correspondence with the magnetograms. In this paper, we investigate this correspondence by comparing ALMA full-disk solar image taken at 1.2 mm with a SDO/HMI magnetogram and analyze their correlation. It is found that chromospheric network and active regions show a positive correlation where brightness temperature is increasing with the line-of-sight magnetic field strength, while sunspots have a negative correlation. Quiet Sun regions do not show any dependence of the brightness temperature with the magnetic field. Thermal bremsstrahlung is given as the best explanation for the observed correlations.

We develop a method for implementing a proposal on utilizing knowledge of neutron star (NS) equation of state (EoS) for inferring the Hubble constant from a population of binary neutron star (BNS) mergers. This method is useful in exploiting BNSs as standard sirens when their redshifts are unavailable. Gravitational wave (GW) signals from compact object binaries provide a direct measurement of their luminosity distances, but not the redshifts. Unlike in the past, we employ a realistic EoS parameterization in a Bayesian framework to simultaneously measure the Hubble constant and refine the constraints on the EoS parameters. The uncertainty in the redshift depends on the uncertainty in EoS and mass parameters estimated from GW data. Combining the inferred BNS redshifts with the corresponding luminosity distances, one constructs a redshift-distance relationship and deduces the Hubble constant from it. Here, we show that in the Cosmic Explorer era, one can measure the Hubble constant to a precision of $\sim 3\%$ (with a $90\%$ credible interval) with a realistic distribution of thousand BNSs, while allowing for uncertainties in their EoS parameters. The methodology implemented in this work demonstrates a comprehensive algorithm to infer NS EoS and the Hubble constant by simultaneously combining GW observations from merging NSs, choosing a simple population model of NS masses and keeping the merger rate of NSs constant. This method can be immediately extended to incorporate merger rate, population properties, and different cosmological parameters.

The sparse layouts of radio interferometers result in an incomplete sampling of the sky in Fourier space which leads to artifacts in the reconstructed images. Cleaning these systematic effects is essential for the scientific use of radiointerferometric images. Established reconstruction methods are often time-consuming, require expert-knowledge, and suffer from a lack of reproducibility. We have developed a prototype Deep Learning-based method that generates reproducible images in an expedient fashion. To this end, we take advantage of the efficiency of Convolutional Neural Networks to reconstruct image data from incomplete information in Fourier space. The Neural Network architecture is inspired by super-resolution models that utilize residual blocks. Using simulated data of radio galaxies that are composed of Gaussian components we train Deep Learning models whose reconstruction capability is quantified using various measures. The reconstruction performance is evaluated on clean and noisy input data by comparing the resulting predictions with the true source images. We find that source angles and sizes are well reproduced, while the recovered fluxes show substantial scatter, albeit not worse than existing methods without fine-tuning. Finally, we propose more advanced approaches using Deep Learning that include uncertainty estimates and a concept to analyze larger images.

We search for a first-order phase transition (PT) gravitational wave (GW) signal from Advanced LIGO and Advanced Virgo's first three observing runs. Due to the large theoretical uncertainties, four GW energy spectral shapes from bubble and sound wave collisions widely adopted in literature are investigated, separately. Our results indicate that there is no evidence for the existence of such GW signals, and therefore we give the upper limits on the amplitude of GW energy spectrum $h^2\Omega_\text{pt}(f_*)$ in the peak frequency range of $f_*\in [0.5,500]$ Hz for these four theoretical models, separately. We find that $h^2\Omega_\text{pt}(f_*\simeq 40\ \text{Hz})<(1.3\sim 1.4)\times10^{-8}$ at $95\%$ credible level, and roughly $H_*/\beta\lesssim 0.1$ and $\alpha\lesssim 1$ at $68\%$ credible level in the peak frequency range of $f_*\in [10,100]$ Hz corresponding to the most sensitive frequency band of Advanced LIGO and Advanced Virgo's first three observing runs, where $H_*$ is the Hubble parameter when PT happens, $\beta$ is the bubble nucleation rate and $\alpha$ is the ratio of vacuum and relativistic energy density.

The Fermi Gamma-ray Burst Monitor (GBM) triggers on-board in response to $\sim$ 40 short gamma-ray bursts (SGRBs) per year; however, their large localization regions have made the search for optical counterparts a challenging endeavour. We have developed and executed an extensive program with the wide field of view of the Zwicky Transient Facility (ZTF) camera, mounted on the Palomar 48 inch Oschin telescope (P48), to perform target-of-opportunity (ToO) observations on 10 Fermi-GBM SGRBs during 2018 and 2020-2021. Bridging the large sky areas with small field of view optical telescopes in order to track the evolution of potential candidates, we look for the elusive SGRB afterglows and kilonovae (KNe) associated with these high-energy events. No counterpart has yet been found, even though more than 10 ground based telescopes, part of the Global Relay of Observatories Watching Transients Happen (GROWTH) network, have taken part in these efforts. The candidate selection procedure and the follow-up strategy have shown that ZTF is an efficient instrument for searching for poorly localized SGRBs, retrieving a reasonable number of candidates to follow-up and showing promising capabilities as the community approaches the multi-messenger era. Based on the median limiting magnitude of ZTF, our searches would have been able to retrieve a GW170817-like event up to $\sim$ 200 Mpc and SGRB afterglows to z = 0.16 or 0.4, depending on the assumed underlying energy model. Future ToOs will expand the horizon to z = 0.2 and 0.7 respectively.

We report the discovery of Pegasus IV, an ultra-faint dwarf galaxy found in archival data from the Dark Energy Camera processed by the DECam Local Volume Exploration Survey. Pegasus IV is a compact, ultra-faint stellar system ($r_{1/2} = 41^{+8}_{-6}$ pc; $M_V = -4.25 \pm 0.2$ mag) located at a heliocentric distance of $90^{+4}_{-6}$ kpc. Based on spectra of seven non-variable member stars observed with Magellan/IMACS, we confidently resolve Pegasus IV's velocity dispersion, measuring $\sigma_{v} = 3.3^{+1.7}_{-1.1} \text{ km s}^{-1}$ (after excluding three velocity outliers); this implies a mass-to-light ratio of $M_{1/2}/L_{V,1/2} = 167^{+224}_{-99} M_{\odot}/L_{\odot}$ for the system. From the five stars with the highest signal-to-noise spectra, we also measure a systemic metallicity of $\rm [Fe/H] = -2.67^{+0.25}_{-0.29}$ dex, making Pegasus IV one of the most metal-poor ultra-faint dwarfs. We tentatively resolve a non-zero metallicity dispersion for the system. These measurements provide strong evidence that Pegasus IV is a dark-matter-dominated dwarf galaxy, rather than a star cluster. We measure Pegasus IV's proper motion using data from Gaia Early Data Release 3, finding ($\mu_{\alpha*}, \mu_{\delta}) = (0.33\pm 0.07, -0.21 \pm 0.08) \text{ mas yr}^{-1}$. When combined with our measured systemic velocity, this proper motion suggests that Pegasus IV is on an elliptical, retrograde orbit, and is currently near its orbital apocenter. Lastly, we identify three potential RR Lyrae variable stars within Pegasus IV, including one candidate member located more than ten half-light radii away from the system's centroid. The discovery of yet another ultra-faint dwarf galaxy strongly suggests that the census of Milky Way satellites is still incomplete, even within 100 kpc.

Blue large-amplitude pulsators (BLAPs) form a small group of hot objects pulsating in a fundamental radial mode with periods of the order of 30 minutes. Proposed evolutionary scenarios explain them as evolved low-mass stars: either ~0.3 M$_\odot$ shell-hydrogen-burning objects with a degenerated helium core, or more massive (0.5 - 0.8) M$_\odot$ core-helium-burning stars, or ~0.7 M$_\odot$ surviving companions of type Ia supernovae. Therefore, their origin remains to be established. Using data from Transiting Exoplanet Survey Satellite, we discovered that HD 133729 is a binary consisting of a late B-type main-sequence star and a BLAP. The BLAP pulsates with a period of 32.37 min decreasing at a rate of $(-7.11 \pm 0.33)\times 10^{-11}$. Due to light dilution by a brighter companion, the observed amplitude of pulsation is much smaller than in other BLAPs. From available photometry, we derived times of maximum light, which revealed the binary nature of the star via O-C diagram. The diagram shows variations with a period of 23.08433 d that we attribute to the light-travel-time effect in the system. The analysis of these variations allowed to derive the spectroscopic parameters of the BLAP's orbit around the center of the mass of the binary. The presence of a hot companion in the system was confirmed by the analysis of its spectral energy distribution, which was also used to place the components in the H-R diagram. The obtained position of the BLAP fully agrees with the location of the other members of the class. With the estimated V~11 mag and the Gaia distance of less than 0.5 kpc, the BLAP is the brightest and the nearest of all known BLAPs. It may become a clue object in the verification of the evolutionary scenarios for this class of variable. We argue that low-mass progenitors of the BLAP are excluded if the components are coeval and no mass transfer between the components took place.

Context: The Helmi streams are kinematic substructures whose progenitor is likely a dwarf galaxy. Although 20 years have past since their discovery, it is still unclear whether their members are chemically distinguishable from other halo stars in the Milky Way. Aim: We aim to precisely characterize the chemical properties of the Helmi streams. Methods: We analysed high-resolution, high-signal-to-noise ratio spectra for eleven Helmi streams' stars through a line-by-line abundance analysis. We compared the derived abundances to homogenized literature abundances of the other halo stars, including those belonging to other kinematic substructures, such as Gaia-Enceladus and Sequoia. Results: Compared to typical halo stars, the Helmi streams' members clearly show low values of [X/Fe] in elements produced by massive stars, such as Na and $\alpha$-elements. This tendency is seen down to metallicities of at least $[\mathrm{Fe/H}]\sim -2.2$, suggesting type~Ia supernovae already started to contribute to the chemical evolution at this metallicity. We found that the [$\alpha$/Fe] ratio does not evolve much with metallicity, making the Helmi streams' stars less distinguishable from Gaia-Enceladus stars at $[\mathrm{Fe/H}]\gtrsim -1.5$. The almost constant but low value of [$\alpha$/Fe] might be indicative of quiescent star formation with low efficiency at the beginning and bursty star formation at later times. We also found extremely low values of [Y/Fe] at low metallicity, providing further support that light neutron capture-elements are deficient in Helmi streams. While Zn is deficient at low metallicity, it shows a large spread at high metallicity. The origin of the extremely low Y abundances and Zn variations remain unclear. Conclusion: The Helmi streams' stars are distinguishable from the majority of the halo stars if homogeneously derived abundances are compared.

Massive stars exhibit a variety of instabilities, many of which are poorly understood. We explore instabilities induced by centrifugal forces and angular momentum transport in massive rotating stars. First, we derive and numerically solve linearized oscillation equations for adiabatic radial modes in polytropic stellar models. In the presence of differential rotation, we show that centrifugal and Coriolis forces combined with viscous angular momentum transport can excite stellar pulsation modes, under both low- or high-viscosity conditions. In the low-viscosity limit, which is common in real stars, we demonstrate how to compute mode growth/damping rates via a work integral. Finally, we build realistic rotating $30\,M_\odot$ star models and show that overstable (growing) radial modes are predicted to exist for most of the star's life, in the absence of non-adiabatic effects. Peak growth rates are predicted to occur while the star is crossing the Hertzsprung-Russell gap, though non-adiabatic damping may dominate over viscous driving, depending on the effective viscosity produced by convective and/or magnetic torques. Viscous instability could be a new mechanism to drive massive star pulsations and is possibly related to instabilities of luminous blue variable stars.

The Cassini spacecraft discovered many close-in small satellites in Saturnian system, and several of them exhibit exotic orbital states due to interactions with Mimas and the oblateness of the planet. This work is devoted to Methone, which is currently involved in a 15:14 Mean-Motion Resonance with Mimas. We give an in deep study the current orbit of Methone by analyzing and identifying the short, resonant and long-term gravitational perturbations on its orbit. In addition, we perform numerical integrations of full equations of motion of ensembles of close-in small bodies orbiting the non-central field of Saturn. Spectral analyses of the orbits and interpretation of them in dynamical maps allow us to describe the orbit and the dynamics of Methone in view of resonant and long-term dynamics. We show that the current geometric orbit of Methone is aligned with Mimas' due to a forced resonant component in eccentricity, leading to simultaneous oscillations of several critical angles of the expanded disturbing function. Thus, we explain the simultaneous oscillations of four critical arguments associated to the resonance. The mapping of the Mimas-Methone resonance shows that the domains of the 15:14 Mimas-Methone resonance are dominated by regular motions associated to the Corotation resonance located at eccentricities lower than $\sim 0.015$ and osculating semi-major axis in the interval 194,660-194,730 km. Methone is currently located deeply within this site.

We explore the evolution of the stellar mass-size relation of galaxies of different morphological types and specifically bulge and disk components. We use a sample of $\sim35,000$ galaxies within a redshift range $0 < z < 1$, and stellar mass $\log_{10}(\mathrm{M}_*/\mathrm{M}_\odot) \geq 9.5$ volume-limited sample drawn from the combined DEVILS and HST-COSMOS region for which we presented a morphological classification into sub-classes of double-component (BD), pure-disk (pD), elliptical (E), and compact (C) in Paper-I and a structural decomposition into disk (D), diffuse bulge (dB), and compact bulge (cB) in Paper-II. We find that compared to disks, ellipticals and bulges follow steeper $M_*-R_e$ relations, likely indicating distinct evolutionary mechanisms. Ellipticals and disk structures follow consistently unchanged slopes of $\sim0.5$ and $\sim0.3$, respectively, at all redshifts. We quantify that disks follow a redshift independent $M_*-R_e$ slope regardless of the presence or absence of a bulge component (i.e., BD or pD) suggesting a similar origin and evolutionary pathway for all disks. Since $z = 1$ compact-bulges present a steepening relation which do not follow that of Es whilst diffuse-bulges experience a modest flattening. Overall, we find a close-to-no variation in the $M_*-R_e$ relations over the last $\sim8$ Gyr suggesting that despite ongoing although declining star-formation, mass evolution, morphological transitions and mergers, evolution moves galaxies along their $M_*-R_e$ trails. This seems to be consistent with an inside-out growth and evolution picture in which galaxies grow in size as they do in stellar mass. Besides, minor mergers are likely to be responsible for the growth of Es, at least in $z < 1$.

Taking advantage of the unprecedented statistical power of upcoming cosmic shear surveys will require exquisite knowledge of the matter power spectrum over a wide range of scales. Analytical methods can achieve such precision only up to quasi-linear scales. For smaller non-linear scales we must resort to $N$-body and hydrodynamical simulations, which despite recent technological advances and improved algorithms remain computationally expensive. Over the past decade machine learning and the advent of emulators have propelled our ability to hit the target accuracy with impressive computing efficiency. Yet, realistically these techniques will be able to produce high-precision non-linear matter power spectra only for a restricted sub-set of extensions to the "vanilla" ${\Lambda}$CDM cosmology. I will present a promising alternative to alleviate these shortcomings that draws strength from the combination of halo model, perturbation theory and emulators--the reaction framework. I will show how a power spectrum evolved in the standard cosmology can be readily adjusted to account for physics beyond ${\Lambda}$CDM, and discuss the accuracy of the reaction for well-known modifications to gravity, dark energy parametrizations and massive neutrino cosmologies.

In this research we test the ability of a three Washington filter combination, (C-T1)-(T1-T2), compared with that of the traditional C-T1 color to find multiple populations on two globular clusters: NGC 7099 and NGC 1851, types I and II Globular clusters, respectively. Our improved photometry and membership selection, now using Gaia proper motions, finds that the second population stars are more centrally concentrated than first population stars, as expected and contrary to our previous findings for NGC 7099. We find that multiple populations are more easily detected in both clusters using the new (C-T1)-(T1-T2) color, although C-T1 conserves the best width/error ratio. We also search for differences of both colors while splitting the red-RGB and the blue-RGB in NGC 1851, but no significant improvement was found.

We derive the first detailed chemical abundances of three star clusters in the Large Magellanic Cloud (LMC), NGC1831 (436+/-22 Myr), NGC1856 (350+/-18 Myr) and [SL63]268 (1230+/-62 Myr) using integrated-light spectroscopic observations obtained with the Magellan Echelle spectrograph on Magellan Baade telescope. We derive [Fe/H], [Mg/Fe], [Ti/Fe], [Ca/Fe], [Ni/Fe], [Mn/Fe], [Cr/Fe] and [Na/Fe] for the three clusters. Overall, our results match the LMC abundances obtained in the literature as well as those predicted by detailed chemical evolution models. For clusters NGC1831 and NGC1856, the [Mg/Fe] ratios appear to be slightly depleted compared to [Ca/Fe] and [Ti/Fe]. This could be hinting at the well-known Mg-Al abundance anti-correlation observed in several Milky Way globular clusters. We note, however, that higher signal-to-noise observations are needed to confirm such a scenario, particularly for NGC 1831. We also find a slightly enhanced integrated-light [Na/Fe] ratio for cluster [SL63]268 compared to those from the LMC field stars, possibly supporting a scenario of intracluster abundance variations. We stress that detailed abundance analysis of individual stars in these LMC clusters is required to confirm the presence or absence of MSPs.

We analyse the aperiodic flaring features, also known as shots, observed in Cyg X-1 in the 0.1-80 keV energy band using a 6.39 ks simultaneous observation with AstroSat and NICER. We detect 49 simultaneous shots in the soft and hard X-ray bands with NICER and AstroSat-LAXPC, respectively. We observe the shot profile for the first time in soft X-rays (0.1-3 keV), which shows a spectral peak at $\sim$2 keV. Using time-averaged spectroscopy, we measured the truncation of the inner accretion disk at $6.7\pm0.2$ gravitational radii. The shot-phase resolved spectroscopy allowed us to identify the origin of some of the brightest aperiodic peaks in the soft X-rays. We find that the accretion rate is consistent with a constant during the shots while the inner edge of the accretion disk moves inwards/outwards as these shots rise/decay. We discuss the possible mechanisms causing the swing in the inner radius.

We describe how analytic solutions for linear hydromagnetic waves can be used for testing cosmological magnetohydrodynamic (MHD) codes. We start from the comoving MHD equations and derive analytic solutions for the amplitude evolution of linear hydromagnetic waves in a matter-dominated, flat Einstein-de-Sitter (EdS) universe. The waves considered are comoving, linearly polarized Alfv\'en waves and comoving, magnetosonic (fast) waves modified by selfgravity. The solution for compressible waves is found for a general adiabatic index and we consider the limits of hydrodynamics without selfgravity in addition to the full solution. In addition to these analytic solutions, the linearized equations are solved numerically for a $\Lambda$CDM cosmology. We use the analytic and numeric solutions to compare with results obtained using the cosmological MHD code AREPO and find good agreement in all cases. We envision that our examples could be useful when developing a new cosmological MHD code or for regression testing of existing codes.

We derive an estimator for the lensing potential from galaxy number counts which contains a linear and a quadratic term. We show that this estimator has a much larger signal-to-noise ratio than the corresponding estimator from intensity mapping. This is due to the additional lensing term in the number count angular power spectrum which is present already at linear order. We estimate the signal-to-noise ratio for future photometric surveys. Particularly at high redshifts, $z\gtrsim 1.5$, the signal to noise ratio can become of order 30. Therefore, the number counts in photometric surveys would be an excellent means to measure tomographic lensing spectra.

Besides their causal connection with long and short-term magnetic variability, solar bipolar magnetic regions are our chief source of insight into the location, size, and properties of large-scale toroidal magnetic structures in the solar interior. The great majority of these regions (~95%) follow a systematic east-west polarity orientation (Hale's law) that reverses in opposite hemispheres and across even and odd cycles. These regions also present a systematic north-south polarity orientation (Joy's law) that helps build the poloidal field that seeds the new cycle. Exceptions to Hale's law are rare and difficult to study due to their low numbers. Here, we present a statistical analysis of the inclination (tilt) with respect to the equator of Hale versus anti-Hale regions spanning four solar cycles, considering two complementary tilt definitions adopted in previous studies. Our results show that anti-Hale regions belong to a separate population than Hale regions, suggesting a different originating mechanism. However, we find that anti-Hale region tilts present similar systematic tilt properties and similar latitudinal distributions to Hale regions, implying a strong connection between the two. We see this as evidence that they belong to a common toroidal flux system. We speculate that anti-Hale regions originate from poloidal field sheared and strengthened on the spot after the emergence of Hale regions with very strong poloidal contribution. Thus, they are not in contradiction with the idea of largely coherent toroidal flux systems inside the solar interior.

Formation of dark matter halos is sensitive to the expansion rate of the Universe and to the growth of structures under gravitational collapse. Virialization of halos heats the gaseous intra-cluster medium to high temperatures, leading to copious emission of photons at X-ray wavelengths. We summarize the progress of X-ray surveys in determining cosmology using galaxy clusters. We review recent cosmological results based on cluster volume abundance, clustering, standard candles, extreme object statistics, and present relevant theoretical considerations. We discuss clusters as gravitation theory probes and present an outlook on future developments.

Bipolar magnetic regions (BMRs) are the cornerstone of solar variability. They are tracers of the large-scale magnetic processes that give rise to the solar cycle, shapers of the solar corona, building blocks of the large-scale solar magnetic field, and significant contributors to the free-energetic budget that gives rise to flares and coronal mass ejections. Surprisingly, no homogeneous catalog of BMRs exists today, in spite of the existence of systematic measurements of the magnetic field since the early 1970's. The purpose of this work is to address this deficiency by creating a homogenous catalog of BMRs from the 1970's until the present. For this purpose, in this paper we discuss the strengths and weaknesses of the automatic and manual detection of BMRs and how both methods can be combined to form the basis of our Bipolar Active Region Detection (BARD) code and its supporting human supervision module. At present, the BARD catalog contains more than 10,000 unique BMRs tracked and characterized during every day of their observation. Here we also discuss our future plans for the creation of an extended multi-scale magnetic catalog combining the SWAMIS and BARD catalogs.

Photon counting detector systems on sounding rocket payloads often require interfacing asynchronous outputs with a synchronously clocked telemetry (TM) stream. Though this can be handled with an on-board computer, there are several low cost alternatives including custom hardware, microcontrollers and field-programmable gate arrays (FPGAs). This paper outlines how a TM interface (TMIF) for detectors on a sounding rocket with asynchronous parallel digital output can be implemented using low cost FPGAs and minimal custom hardware. Low power consumption and high speed FPGAs are available as commercial off-the-shelf (COTS) products and can be used to develop the main component of the TMIF. Then, only a small amount of additional hardware is required for signal buffering and level translating. This paper also discusses how this system can be tested with a simulated TM chain in the small laboratory setting using FPGAs and COTS specialized data acquisition products.

The solar cycle periodically reshapes the magnetic structure and radiative output of the Sun and determines its impact on the heliosphere roughly every 11 years. Besides this main periodicity, it shows century-long variations (including periods of abnormally low solar activity called grand minima). The Maunder Minimum (1645-1715) has generated significant interest as the archetype of a grand minimum in magnetic activity for the Sun and other stars, suggesting a potential link between the Sun and changes in terrestrial climate. Recent reanalyses of sunspot observations have yielded a conflicted view on the evolution of solar activity during the past 400 years (a steady increase versus a constant level). This has ignited a concerted community-wide effort to understand the depth of the Maunder Minimum and the subsequent secular evolution of solar activity. The goal of this Perspective is to review recent work that uses historical data to estimate long-term solar variability, and to provide context to users of these estimates that may not be aware of their limitations. We propose a clear visual guide than can be used to easily assess observational coverage for different periods, as well as the level of disagreement between currently proposed sunspot group number series.

We investigate the differences in the small-scale structure of vector dark matter (VDM) and scalar dark matter (SDM) using 3+1 dimensional simulations of single/multicomponent Schr\"{o}dinger-Poisson system. We find that the amount of wave interference, core to halo mass ratio (and its scatter), spin of the core, as well as the shape of the central regions of dark matter halos can distinguish VDM and SDM. Starting with a collection of idealized halos (self-gravitating solitons) as an initial condition, we show that the system dynamically evolves to an approximately spherically symmetric configuration that has a core surrounded by a halo of interference patterns in the mass density. In the vector case, the central soliton in less dense and has a smoother transition to an $r^{-3}$ tail compared to the scalar case. Wave interference is $\sim 1/\sqrt{3}$ times smaller in VDM compared to SDM, resulting in fewer low and high density regions in VDM compared to SDM, with more diffuse granules in the halo. The ratio of VDM core mass to the total halo mass is lower than that in SDM, with a steeper dependence on the total energy of the system and a slightly larger scatter. Finally, we also initiate a study of the evolution of intrinsic spin angular momentum in the VDM case. We see a positive correlation between the total intrinsic spin in the simulation and the spin of the final central core, with significant scatter. We see large intrinsic spin in the core being possible even with vanishing amounts total angular momentum in the initial conditions. Our results point towards the possibility of distinguishing VDM from SDM using astrophysical and terrestrial observations.

The detection of an astrophysical flux of neutrinos in the TeV-PeV energy range by the IceCube observatory has opened new possibilities for the study of extreme cosmic accelerators. The apparent isotropy of the neutrino arrival directions favors an extragalactic origin for this flux, potentially created by a large population of distant sources. Recent evidence for the detection of neutrino emission from extragalactic sources include the active galaxies TXS 0506+056 and NGC 1068. We here review the current status of the search for the sources of the high-energy neutrino flux, concentrating on its extragalactic contribution. We discuss the implications of these observations for multimessenger studies of cosmic sources and present an outlook for how additional observations by current and future instruments will help address fundamental questions in the emerging field of high-energy neutrino astronomy.

Quasi-periodic eruptions (QPEs) are recurrent X-ray bursts found so far in the nuclei of low-mass galaxies. Their trigger mechanism is still unknown, but recent models involving one or two stellar-mass companions around the central massive ($\approx10^5-10^6$ solar masses) black hole have gathered significant attention. While these have been compared only qualitatively with observations, the phenomenology of QPEs is developing at a fast pace, with the potential to reveal new insights. Here we report two new observational results found in eRO-QPE1, the brightest QPE source discovered so far: i) the eruptions in eRO-QPE1 occur sometimes as single isolated bursts, and at others as chaotic mixtures of multiple overlapping bursts with very different amplitudes; ii) we confirm that QPEs peak at later times and are broader at lower energies, with respect to higher energies while, for the first time, we find that QPEs also start earlier at lower energies. Furthermore, eruptions appear to undergo an anti-clockwise hysteresis cycle in a plane of hardness ratio versus total count rate. Behavior i) was not found before in any other QPE source and implies that if a common trigger mechanism is in place for all QPEs, it must be able to produce both types of timing properties, regular and complex. Result ii) implies that the X-ray emitting component does not have an achromatic evolution even during the start of QPEs, and that the rise is harder than the decay at a given total count rate. This specific energy dependence could be qualitatively compatible with inward radial propagation during the rise within a compact accretion flow, the presence of which is suggested by the stable quiescence spectrum observed in general for QPE sources.

We estimate the stochastic gravitational wave spectrum emitted from a network of cosmic strings in which the latter are effectively stable against breaking by monopole pair creation. The monopoles are produced at a higher scale from an earlier symmetry breaking and experience significant inflation before reentering the horizon. This gives rise to monopole-antimonopole pairs connected by string segments and the string loop formation essentially ceases. As a consequence, the lower frequency portion of the gravitational wave spectrum is suppressed relative to the no-inflation case with stable strings, which evades the stringent PPTA bound on the dimensionless string tension $G\mu$. We display the modified spectrum, accessible in the ongoing and future experiments, for $G\mu$ values in the range $10^{-10} - 10^{-15}$. We show how this `quasi-stable' string network is realized in realistic grand unified theories.

We argue that standard tools of holography can be used to describe fully non-perturbative microscopic models of cosmology in which a period of accelerated expansion may result from the positive potential energy of time-dependent scalar fields evolving towards a region with negative potential. In these models, the fundamental cosmological constant is negative, and the universe eventually recollapses in a time-reversal symmetric way. The microscopic description naturally selects a special state for the cosmology. In this framework, physics in the cosmological spacetime is dual to the vacuum physics in a static planar asymptotically AdS Lorentzian wormhole spacetime, in the sense that the background spacetimes and observables are related by analytic continuation. The dual spacetime is weakly curved everywhere, so any cosmological observables can be computed in the dual picture via effective field theory without detailed knowledge of the UV completion or the physics near the big bang. In particular, while inflation may explain the origin of perturbations in the cosmology picture, the perturbations can be deduced from the dual picture without any knowledge of the inflationary potential.

Repulsive gravity is a well known characteristic of naked singularities. In this work, we explore light surfaces and find new effects of repulsive gravity. We compare Kerr naked singularities with the corresponding black hole counterparts and find certain structures that are identified as horizon remnants. We argue that these features might be significant for the comprehension of processes that lead to the formation or eventually destruction of black hole Killing horizons. These features can be detected by observing photon orbits, particularly close to the rotation axis, which can be used to distinguish naked singularities from black hole.

In the context of Effective Field Theory, the Hilbert space of states increases in an expanding universe. Hence, the time evolution cannot be unitary. The formation of structure is usually studied using effective field theory techniques. We study the constraints on effective field theory analyses of early universe models which come from demanding that the factor of the space of states corresponding to length scales where the primordial fluctuations are manifest does not suffer from the unitarity problem. For bouncing and emergent cosmologies, no constraints arise provided that the energy scale of the bounce or emergent phases is smaller than the ultraviolet (UV) cutoff scale. On the other hand, in the case of the inflationary scenario, non-trivial upper bounds on the energy scale of inflation arise.

With the improved accuracy of neutron star observational data, it is necessary to derive new equation of state where the crust and the core are consistently calculated within a unified approach. For this purpose we describe non-uniform matter in the crust of neutron stars employing a compressible liquid-drop model, where the bulk and the neutron fluid terms are given from the same model as the one describing uniform matter present in the core. We then generate a set of fifteen unified equations of state for cold catalyzed neutron stars built on realistic modelings of the nuclear interaction, which belongs to two main groups: the first one derives from the phenomenological Skyrme interaction and the second one from $\chi_{EFT}$ Hamiltonians. The confrontation of these model predictions allows us to investigate the model dependence for the crust properties, and in particular the effect of neutron matter at low density. The new set of unified equations of state is available at the CompOSE repository.

In this proceedings, we are interested in dark gravitational wave standard sirens and their use for cosmology and for constraining modified gravity theories. Due to the extra friction term introduced in their propagation equation those theories predict different luminosity distances for electromagnetic and gravitational waves (GWs). This effect can be parametrized by the two variables $\Xi_0$ and $n$, that can be measured from gravitational wave observations, and specifically from the binary black hole (BBH) mergers detected by LIGO and Virgo. By fitting jointly BBH population models in mass and redshift, the cosmological parameters, and the modified GW luminosity distance to $\sim$ 60 signals observed during the first three LIGO/Virgo observation runs, we conclude that general relativity is consistently the preferred model. The future observation runs O4 and O5 are also considered. Using the same approach, we forecast a measurement uncertainty on the modified gravity parameter $\Xi_0$ of $51\%$ with O4, and $20\%$ with O4 and O5, respectively if GR is the correct theory of gravity. However, we underline that there are strong correlations between astrophysical, cosmological and modified gravity parameters, possibly leading to bias if wrong priors are assumed.

We consider gravitational wave signals produced by a first-order phase transition in a theory with a generic renormalizable thermal effective potential of power law form. We find the frequency and amplitude of the gravitational wave signal can be related in a straightforward manner to the parameters of the thermal effective potential. This leads to a general conclusion; if the mass of the dark Higgs is less than 1% of the dark Higgs vacuum expectation value, then the gravitational wave signal will be unobservable at all upcoming and planned gravitational wave observatories.

Non-Standard Interactions (NSI) between neutrinos and matter affect the neutrino flavor oscillations. Due to the high matter density in the core of the Sun, solar neutrinos are suited to probe these interactions. Using the $277$ kton-yr exposure of Super-Kamiokande to $^{8}$B solar neutrinos, we search for the presence of NSI. Our data favors the presence of NSI with down quarks at 1.8$\sigma$, and with up quarks at 1.6$\sigma$, with the best fit NSI parameters being ($\epsilon_{11}^{d},\epsilon_{12}^{d}$) = (-3.3, -3.1) for $d$-quarks and ($\epsilon_{11}^{u},\epsilon_{12}^{u}$) = (-2.5, -3.1) for $u$-quarks. After combining with data from the Sudbury Neutrino Observatory and Borexino, the significance increases by 0.1$\sigma$.

In this paper, we introduce a self-consistent mean field approximation to study the QCD phase transition and the structure of hybrid stars within the framework of NJL model. In our practice, a phenomenological parameter $\alpha$ is introduced, which reflects the weights of "direct" channel and "exchange" channel under a finite chemical potential. The mass-radius relation is obtained by solving the Tolman-Oppenheimer-Volkoff equation using a crossover equation of state (EOS). We calculate the density distribution in a two solar-mass hybrid star to show the effects of different parameters. We also calculate the tidal Love number $k_2$ and the deformability $(\Lambda)$. It is found that the stiffness of the EOS increases with $\alpha $, which allows us to obtain a hybrid star with a maximum mass of 2.40 solar-mass through our model. The observation of over 2.06 solar-mass neutron stars may indicates that the chiral transition may be a crossover on the whole $T-\mu$ plane.