Papers for Tuesday, Aug 16 2022

With the rapid increase of fast radio burst (FRB) detections within the past few years, there is now a catalogue being developed for all-sky extragalactic dispersion measure (DM) observations in addition to the existing collection of all-sky extragalactic Faraday rotation measurements (RMs) of radio galaxies. We present a method of reconstructing all-sky information of the Galactic magnetic field component parallel to the line of sight, $B_{\parallel}$, using simulated observations of the RM and DM along lines of sight to radio galaxies and FRB populations, respectively. This technique is capable of distinguishing between different input Galactic magnetic field and thermal electron density models. Significant extragalactic contributions to the DM are the predominant impediment in accurately reconstructing the Galactic DM and $\left<B_{\parallel}\right>$ skies. We look at ways to improve the reconstruction by applying a filtering algorithm on the simulated DM lines of sight and we derive generalized corrections for DM observations at $|b|$ > 10 deg that help to disentangle Galactic and extragalactic DM contributions. Overall, we are able to reconstruct both large-scale Galactic structure and local features in the Milky Way's magnetic field from the assumed models. We discuss the application of this technique to future FRB observations and address possible differences between our simulated model and observed data, namely: adjusting the priors of the inference model, an unevenly distributed population of FRBs on the sky, and localized extragalactic DM structures.

JWST's Early Release Observations of the lensing cluster SMACS J0723.3-7327 have given an unprecedented spectroscopic look into the high-redshift universe. These observations reveal five galaxies at z > 5. All five have detectable [OIII]4363 line emission, indicating that these galaxies have high temperatures and low metallicities and that they are highly star forming. In recent work, the metallicities of these five galaxies have been studied using various techniques. Here we summarize and compare these previous results, as well as perform our own measurements of the metallicities using improved methodologies that optimize the extraction of the emission lines. In particular, we use simultaneous line fitting and a fixed Balmer decrement correction, as well as a novel footprint measurement of the emission lines in the 2D spectra, to produce higher fidelity line ratios that are less sensitive to calibration and systematic effects. We then compare our metallicities to those of z < 1 galaxies with high rest-frame equivalent widths of H-beta, finding that they may be good analogs. Finally, we estimate that the JWST galaxies out to z ~ 8 are young compared to the age of the universe.

We apply the radiative transfer (RT) code SKIRT on a sample of ~14000 low-redshift (z<= 0.1) galaxies extracted from the TNG50 simulation to enable an apples-to-apples comparison with observations. The RT procedure is calibrated via comparison of a subsample of TNG50 galaxies with the DustPedia observational sample: we compare several luminosity and colour scaling relations and spectral energy distributions in different specific SFR bins. We consistently derive galaxy luminosity functions for the TNG50 simulation in 14 broadband filters from UV to submillimetre wavelengths and investigate the effects of the aperture, orientation, radiative transfer recipe, and numerical resolution. We find that, while our TNG50+RT fiducial model agrees well with the observed luminosity functions at the knee (+/- 0.04 dex typical agreement), the TNG50+RT luminosity functions evaluated within 5R_1/2 are generally higher than observed at both the faint and bright ends, by 0.004 (total IR)-0.27 (UKIDSS H) dex and 0.12 (SPIRE250)-0.8 (GALEX FUV) dex, respectively. A change in the aperture does affect the bright end of the luminosity function, easily by up to 1 dex depending on the choice. However, we also find that the galaxy luminosity functions of a worse-resolution run of TNG50 (TNG50-2, with 8 times worse mass resolution than TNG50, similar to TNG100) are in better quantitative agreement with observational constraints. Finally, we publicly release the photometry for the TNG50 sample in 53 broadbands from FUV to submillimetre, in three orientations and four apertures, as well as galaxy spectral energy distributions.

Galaxy mergers are known to trigger both extended and central star formation. However, what remains to be understood is whether this triggered star formation is facilitated by enhanced star formation efficiencies, or an abundance of molecular gas fuel. This work presents spatially resolved measurements of CO emission collected with the Atacama Large Millimetre Array (ALMA) for 20 merging galaxies (either pairs or post-mergers) selected from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey. Eleven additional merging galaxies are selected from the ALMA MaNGA QUEnching and STar formation (ALMaQUEST) survey, resulting in a set of 31 mergers at various stages of interaction and covering a broad range of star formation rates (SFR). We investigate galaxy-to-galaxy variations in the resolved Kennicutt-Schmidt relation (rKS: $\Sigma_{H_2}$ vs. $\Sigma_{SFR}$), the resolved molecular gas main sequence (rMGMS: $\Sigma_{\star}$ vs. $\Sigma_{H_2}$), and the resolved star-forming main sequence (rSFMS: $\Sigma_{\star}$ vs. $\Sigma_{SFR}$). We quantify offsets from these resolved relations to determine if star formation rate, molecular gas fraction, and/or star formation efficiency (SFE) is enhanced in different regions of an individual galaxy. By comparing offsets in all three parameters we can discern whether gas fraction or SFE powers an enhanced $\Sigma_{SFR}$. We find that merger-induced star formation can be driven by a variety of mechanisms, both within a galaxy and between different mergers, regardless of interaction stage.

Observing habitable exoplanets that may resemble Earth is a key priority in astronomy that is dependent on not only detecting such worlds, but also ascertaining that apparent signatures of habitability are not due to other sources. Space telescopes designed to observe such worlds, such as that recommended by NASA's 2020 Astrophysics Decadal Survey, have a diffraction-limited resolution that effectively spreads light from a source in a region around the source point. In this letter, we show that the diffraction limit of a 6 meter space telescope results in a point spread function of an Earth-like planet that may contain additional unanticipated bodies for systems at distances relevant to proposed searches. These unexpected additional objects, such as other planets and moons, can influence obtained spectra for a putative habitable planet by producing spurious features and adding additional uncertainty in the spectra. A model of the Earth observed by a 6 meter space telescope as though it was an exoplanet shows that the light from the Earth would be blended with the Moon, Mercury, Venus and Mars in various combinations and at different times for numerous combinations of distance to the system and wavelength. Given the importance of extricating the true spectra of a potentially habitable planet in order to search for biosignatures, we highlight the need to account for this effect during the development of relevant telescopes and suggest some potential means of accounting for this photobombing effect.

We have studied the nuclear region of the previously detected dual AGN system in the galaxy pair IRAS 05589+2828 and 2MASX J06021107+2828382 through new optical spectroscopy observations, along with radio and X-ray archival data. Our multiwavelength data strongly suggest that the Sy1 \iras\, (z=0.0330$\pm$0.0002) conforms to a dual AGN system with the Sy2 \twomas\, (z=0.0334$\pm$0.0001) with a projected separation obtained from the radio data of 20.08\arcsec\, ($\sim$13.3\,kpc). Analysis of the optical spectra reveals a faint narrow extended emission from H$\alpha$ and [OIII] amidst the two AGN, supporting evidence for an ongoing merger. \iras\, is a double component narrow emission line AGN, with complex broad Balmer emission line profiles that clearly show a strong red-peaklet with a velocity shift of $\sim$3500\,km\,s$^{-1}$. The black hole mass estimates of \iras\, and \twomas\, are log\,M$\rm_{BH}$\,=\,8.59\,$\pm$\,0.14 (M$_\odot$) and log\,M$\rm_{BH}$\,=\,8.21$\pm$0.2 (M$_\odot$), respectively. In the X-ray bands, \iras\, is compatible with a Type 1 object, showing both spectral and flux variability. \chandra\, data of 2MASX\,J06021107+2828382 allowed us to measure a high hardness ratio in this source, providing evidence for a Type 2 AGN. The 22 GHz image obtained with the {\it Karl G. Jansky Very Large Array} has revealed that both AGN are compact radio objects with spectral indices -0.26$\pm$0.03 and -0.70$\pm$0.11, confirming for the first time its dual AGN nature in the radio bands.

We present an analysis of the kinematics of the Radcliffe Wave, a 2.7-kpc-long sinusoidal band of molecular clouds in the solar neighborhood recently detected via 3D dust mapping. With Gaia DR2 astrometry and spectroscopy, we analyze the 3D space velocities of $\sim 1500$ young stars along the Radcliffe Wave in action-angle space, using the motion of the wave's newly born stars as a proxy for its gas motion. We find that the vertical angle of young stars -- corresponding to their orbital phase perpendicular to the Galactic plane -- varies significantly as a function of position along the structure, in a pattern potentially consistent with a wave-like oscillation. This kind of oscillation is not seen in a control sample of older stars from Gaia occupying the same volume, disfavouring formation channels caused by long-lived physical processes. We use a ``wavy midplane'' model to try to account for the trend in vertical angles seen in young stars, and find that while the best-fit parameters for the wave's spatial period and amplitude are qualitatively consistent with the existing morphology defined by 3D dust, there is no evidence for additional velocity structure. These results support more recent and/or transitory processes in the formation of the Radcliffe Wave, which would primarily affect the motion of the wave's gaseous material. Comparisons of our results with new and upcoming simulations, in conjunction with new stellar radial velocity measurements in Gaia DR3, should allow us to further discriminate between various competing hypotheses.

We present HD-TESS, a catalog of 1,709 bright ($V\sim3-10$) red giants from the Henry Draper (HD) Catalog with asteroseismic measurements based on photometry from NASA's Transiting Exoplanet Survey Satellite (TESS). Using light curves spanning at least six months across a single TESS observing cycle, we provide measurements of global asteroseismic parameters ($\nu_{\mathrm{max}}$ and $\Delta\nu$) and evolutionary state for each star in the catalog. We adopt literature values of atmospheric stellar parameters to estimate the masses and radii of the giants in our catalog using asteroseismic scaling relations, and observe that HD-TESS giants on average have larger masses compared to Kepler red giants. Additionally, we present the discovery of oscillations in 99 red giants in astrometric binary systems, including those with subdwarf or white dwarf companions. Finally, we benchmark radii from asteroseismic scaling relations against those measured using long-baseline interferometry for 18 red giants and find that correction factors to the scaling relations improve the agreement between asteroseismic and interferometric radii to approximately 3%.

Linearly polarized emission from dust grains and molecular spectroscopy is an effective probe of the magnetic field topology in the interstellar medium and molecular clouds. The longstanding Davis-Chandrasekhar-Fermi (DCF) method and the recently developed Histogram of Relative Orientations (HRO) analysis and the polarization-intensity gradient (KTH) method are widely used to assess the dynamic role of magnetic fields in star formation based on the plane-of-sky component of field orientations inferred from the observations. We review the advances and limitations of these methods and summarize their applications to observations. Numerical tests of the DCF method, including its various variants, indicate that its largest uncertainty may come from the assumption of energy equipartition, which should be further calibrated with simulations and observations. We suggest that the ordered and turbulent magnetic fields of particular observations are local properties of the considered region. An analysis of the polarization observations using DCF estimations suggests that magnetically trans-to-super-critical and averagely trans-to-super-Alfv\'{e}nic clumps/cores form in sub-critical clouds. High-mass star-forming regions may be more gravity-dominant than their low-mass counterparts due to higher column density. The observational HRO studies clearly reveal that the preferential relative orientation between the magnetic field and density structures changes from parallel to perpendicular with increasing column densities, which, in conjunction with simulations, suggests that star formation is ongoing in trans-to-sub-Alfv\'{e}nic clouds. There is a possible transition back from perpendicular to random alignment at higher column densities. Results from observational studies using the KTH method broadly agree with those of the HRO and DCF studies.

Pulsar timing arrays (PTAs) are searching for gravitational waves from supermassive black hole binaries (SMBHBs). Here we show how future PTAs could use a detection of gravitational waves from individually resolved SMBHB sources to produce a purely gravitational wave-based measurement of the Hubble constant. This is achieved by measuring two separate distances to the same source from the gravitational wave signal in the timing residual: the luminosity distance $D_L$ through frequency evolution effects, and the parallax distance $D_\mathrm{par}$ through wavefront curvature (Fresnel) effects. We present a generalized timing residual model including these effects in an expanding universe. Of these two distances, $D_\mathrm{par}$ is challenging to measure due to the pulsar distance wrapping problem, a degeneracy in the Earth-pulsar distance and gravitational wave source parameters that requires highly precise, sub-parsec level, pulsar distance measurements to overcome. However, in this paper we demonstrate that combining the knowledge of two SMBHB sources in the timing residual largely removes the wrapping cycle degeneracy. Two sources simultaneously calibrate the PTA by identifying the distances to the pulsars, which is useful in its own right, and allow recovery of the source luminosity and parallax distances which results in a measurement of the Hubble constant. We find that, with optimistic PTAs in the era of the Square Kilometer Array, two SMBHB sources within a few hundred Mpc could be used to measure the Hubble constant with a relative uncertainty on the order of 10 per cent.

The deuteration of molecules forming in the ices such as methanol (CH$_3$OH) is sensitive to the physical conditions during their formation in dense cold clouds and can be probed through observations of deuterated methanol in hot cores. Observations with ALMA containing transitions of CH$_3$OH, CH$_2$DOH, CHD$_2$OH, $^{13}$CH$_3$OH, and CH$_3^{18}$OH are investigated. The column densities of CH$_2$DOH, CHD$_2$OH, and CH$_3$OH are determined for all sources, where the column density of CH$_3$OH is derived from optically thin $^{13}$C and $^{18}$O isotopologues. Consequently, the D/H ratio of methanol is derived taking statistical effects into account. Singly deuterated methanol (CH$_2$DOH) is detected toward 25 of the 99 sources in our sample of the high-mass protostars. Including upper limits, the $\rm (D/H)_{CH_3OH}$ ratio inferred from $N_\mathrm{CH_2DOH}/N_\mathrm{CH_3OH}$ was derived for 38 of the 99 sources and varies between $\sim10^{-3}-10^{-2}$. Including other high-mass hot cores from the literature, the mean methanol D/H ratio is $1.1\pm0.7\times10^{-3}$. This is more than one order of magnitude lower than what is seen for low-mass protostellar systems ($2.2\pm1.2\times10^{-2}$). Doubly deuterated methanol (CHD$_2$OH) is detected toward 11 of the 99 sources. Including upper limits for 15 sources, the $\rm (D/H)_{CH_2DOH}$ ratios derived from $N_\mathrm{CHD_2OH}/N_\mathrm{CH_2DOH}$ are more than two orders of magnitude higher than $\rm (D/H)_{CH_3OH}$ with an average of $2.0\pm0.8\times10^{-1}$ which is similar to what is found for low-mass sources. Comparison with literature GRAINOBLE models suggests that the high-mass prestellar phases are either warm ($>20$ K) or live shorter than the free-fall timescale. In contrast, for low-mass protostars, both a low temperature of $<15$ K and a prestellar phase timescale longer than the free-fall timescale are necessary.

We present a rest-frame optical morphological analysis of galaxies observed with the NIRCam imager on the James Webb Space Telescope (JWST) as part of the GLASS Early Release Science program. We select 217 sources at redshifts 0.8 < z < 5.4 and use the seven 0.9-5um NIRCam filters to generate rest-frame gri composite color images, and conduct visual morphological classification. Compared to HST-based work we find a higher incidence of disks and bulges than expected at z > 1.5, revealed by rest frame optical imaging. We detect 73 clear disks (43 at z > 1.5) of which 45 have bulges. No evolution of bulge fraction with redshift is evident: 60% at z < 2 (N = 26) versus 64% at z >= 2 (N = 19). A stellar mass dependence is evident, with bulges visible in 47% of all disk galaxies with mass 10^9.5 M_sun (N = 30) but only 9% at M < 10^9.5 M_sun (N = 15). We supplement visual morphologies with non-parametric measurements of Gini and Asymmetry coefficients in the rest-frame i-band. Our sources are more asymmetric than local galaxies, with slightly higher Gini values. When compared to high-z rest-frame ultraviolet measurements with Hubble Space Telescope, JWST shows more regular morphological types such as disks, bulges and spiral arms at z > 1.5, with smoother (i.e. lower Gini) and more symmetrical light distributions. The unexpected prevalence of detectable bulges and regular morphology at z > 1.5 will allow new tests of theoretical models of galaxy evolution.

We present a thermal evolution model coupled with a Henyey solver to study the circumstances under which a rocky planet could potentially host a dynamo in its liquid iron core and/or magma ocean. We calculate the evolution of planet thermal profiles by solving the energy balance equations for both the mantle and the core. We use a modified mixing length theory to model the convective heat flow in both the magma ocean and solid mantle. In addition, by including the Henyey solver, we self-consistently account for adjustments in the interior structure and heating (cooling) due to planet contraction (expansion). We evaluate whether a dynamo can operate using the critical magnetic Reynolds number. We run simulations to explore how planet mass ($M_{pl}$), core mass fraction (CMF) and equilibrium temperature ($T_{eq}$) affect the evolution and lifetime of possible dynamo sources. We find that the $T_{eq}$ determines the solidification regime of the magma ocean, and only layers with melt fraction greater than a critical value of 0.4 may contribute to the dynamo source region in the magma ocean. We find that the mantle mass, determined by $M_{pl}$ and CMF, controls the thermal isolating effect on the iron core. In addition, we show that the liquid core last longer with increasing planet mass. For a core thermal conductivity of 40$\ \mathrm{Wm^{-1}K^{-1}}$, the lifetime of the dynamo in the iron core is limited by the lifetime of the liquid core for 1$M_{\oplus}$ planets, and by the lack of thermal convection for 3$M_{\oplus}$ planets.

With a rapidly rising number of transients detected in astronomy, classification methods based on machine learning are increasingly being employed. Their goals are typically to obtain a definitive classification of transients, and for good performance they usually require the presence of a large set of observations. However, well-designed, targeted models can reach their classification goals with fewer computing resources. This paper presents SNGuess, a model designed to find young extragalactic nearby transients with high purity. SNGuess works with a set of features that can be efficiently calculated from astronomical alert data. Some of these features are static and associated with the alert metadata, while others must be calculated from the photometric observations contained in the alert. Most of the features are simple enough to be obtained or to be calculated already at the early stages in the lifetime of a transient after its detection. We calculate these features for a set of labeled public alert data obtained over a time span of 15 months from the Zwicky Transient Facility (ZTF). The core model of SNGuess consists of an ensemble of decision trees, which are trained via gradient boosting. Approximately 88% of the candidates suggested by SNGuess from a set of alerts from ZTF spanning from April 2020 to August 2021 were found to be true relevant supernovae (SNe). For alerts with bright detections, this number ranges between 92% and 98%. Since April 2020, transients identified by SNGuess as potential young SNe in the ZTF alert stream are being published to the Transient Name Server (TNS) under the AMPEL_ZTF_NEW group identifier. SNGuess scores for any transient observed by ZTF can be accessed via a web service. The source code of SNGuess is publicly available.

Ultra-light bosonic dark matter called fuzzy dark matter (FDM) has attracted much attention as an alternative to the cold dark matter. An intriguing feature of the FDM model is the presence of a soliton core, a stable dense core formed at the center of halos. In this paper, we analytically study the dependence of the soliton core properties on the halo characteristics by solving approximately the Schr\"odinger-Poisson equation. Focusing on the ground-state eigenfunction, we derive a key expression for the soliton core radius, from which we obtain the core-halo mass relations similar to those found in the early numerical work, but involving the factor dependent crucially on the halo concentration and cosmological parameters. Based on the new relations, we find that for a given cosmology, (i) there exist a theoretical bound on the radius and mass of soliton core for each halo mass (ii) incorporating the concentration-halo mass (c-M) relation into the predictions, the core-halo relations generally exhibit a non power-law behavior, and with the c-M relation suppressed at the low-mass scales, relevant to the FDM model, predictions tend to match the simulations well (iii) the scatter in the c-M relation produces a sizable dispersion in the core-halo relations, and can explain the results obtained from cosmological simulations. Finally, the validity of our analytical treatment are critically examined. A perturbative estimation suggests that the prediction of the core-halo relations is valid over a wide range of parameter space, and the impact of the approximation invoked in the analytical calculations is small, although it is not entirely negligible.

RX J0134.2-4258 is one of the most super-Eddington narrow-line Seyfert 1 (NLS1) galaxies, on which we conducted a monitoring campaign from radio to X-rays. In this paper, we present a detailed analysis of its optical/UV spectra and broadband spectral energy distribution (SED). Our study shows that the preferred black hole mass of RX J0134.2-4258 is $\sim 2~\times~10^{7}~M_{\odot}$, giving a mass accretion rate through the outer disc of $\dot{m}_{\rm out} \sim 20$ (assuming zero spin) compared to the observed luminosity ratio $L_{\rm bol}/L_{\rm Edd} \sim 6$. This reduction in radiative efficiency is expected for super-Eddington flows, as power can be lost via advection and/or winds. We find that the optical/UV lines of RX J0134.2-4258 resemble those from weak-like quasars (WLQs), as it has notably weak C IV and N V emission lines. It also has dramatic X-ray variability, again similar to that recently observed in some other WLQs. However, WLQs have systematically higher masses ($\gtrsim 10^8~M_{\odot}$), and lower Eddington ratios ($\dot{m}_{\rm out} \sim 1-3$) than RX J0134.2-4258. We compare instead to the most extreme NLS1s, with similarly large $\dot{m}_{\rm out}$ but smaller masses. These show similarly large reductions in radiative efficiency but their UV lines are not similarly wind-dominated. We suggest a new category of weak-line Seyfert (WLS) galaxies to describe sources like RX J0134.2-4258, and interpret its (so far unique) properties in a model where the lower disc temperature in the higher-mass black holes leads to UV line driving enhancing the super-Eddington radiation pressure driven wind.

The central part of the Galaxy host a multitude of stellar populations, including the spheroidal bulge stars, stars moved to the bulge through secular evolution of the bar, inner halo, inner thick disk, inner thin disk, as well as debris from past accretion events. We identified a sample of 58 candidate stars belonging to the stellar population of the spheroidal bulge, and analyse their abundances. The present calculations of Mg, Ca, and Si lines are in agreement with the APOGEE-ASPCAP abundances, whereas abundances of C, N, O, and Ce are re-examined. We find normal $\alpha$-element enhancements in oxygen, similar to magnesium, Si, and Ca abundances, which are typical of other bulge stars surveyed in the optical in Baade's Window. The enhancement of [O/Fe] in these stars suggests that they do not belong to accreted debris. No spread in N abundances is found, and none of the sample stars is N-rich, indicating that these stars are not second generation stars originated in globular clusters. Ce instead is enhanced in the sample stars, which points to an s-process origin such as due to enrichment from early generations of massive fast rotating stars, the so-called spinstars

The solutions of the \textit{diffuse reflection finite atmosphere problem} are very useful in the astrophysical context. Chandrasekhar was the first to solve this problem analytically, by considering atmospheric scattering. These results have wide applications in the modeling of planetary atmospheres. However, they cannot be used to model an atmosphere with emission. We solved this problem by including \textit{thermal emission effect} along with scattering.Here, our aim is to provide a complete picture of generalized finite atmosphere problem in presence of scattering and thermal emission, and to give a physical account of the same. For that, we take an analytical approach using the invariance principle method to solve the diffuse reflection finite atmosphere problem in the presence of atmospheric thermal emission. We established the general integral equations of modified scattering function $S(\tau; \mu, \phi; \mu_0, \phi_0)$, transmission function $T(\tau; \mu, \phi; \mu_0, \phi_0)$ and their derivatives with respect to $\tau$ for a thermally emitting atmosphere. We customize these equations for the case of isotropic scattering and introduce two new functions $V(\mu)$ and $W(\mu)$, analogous to Chandrasekhar $X(\mu)$, and $Y(\mu)$ functions respectively. We also derive a transformation relation between the modified S-T functions and give a physical account of $V(\mu)$ and $W(\mu)$ functions. Our final results are consistent with those of Chandrasekhar at low emission limit (i.e. only scattering). From the consistency of our results, we conclude that the consideration of thermal emission effect in diffuse reflection finite atmosphere problem gives more general and accurate results than considering only scattering.

Ultraviolet (UV) wavelength observations have made a significant contribution to our understanding of hot stellar populations of star clusters. Multi-wavelength spectral energy distributions (SEDs) of stars, including ultraviolet observations, have proven to be an excellent tool for discovering unresolved hot companions in exotic stars such as blue straggler stars (BSS), thereby providing helpful clues to constrain their formation mechanisms. Melotte 66 is a 3.4 Gyr old open cluster located at a distance of 4810 pc. We identify the cluster members by applying the ML-MOC algorithm on Gaia EDR3 data. Based on our membership identification, we find 1162 members, including 14 BSS candidates, 2 yellow straggler candidates (YSS), and one subdwarf B candidate (sdB). We generated SEDs for 11 BSS candidates and the sdB candidate using Swift/UVOT data combined with other archival data in the optical and IR wavelengths. We discover a hot companion of one BSS candidate, BSS3, with temperature of 38000$_{-6000}^{+7000}$ K, luminosity of 2.99$_{-1.86}^{+5.47}$ L$_\odot$, and radius of 0.04$_{-0.005}^{+0.008}$ R$_\odot$. This hot companion is a likely low-mass WD with an estimated mass of 0.24 $\text{-}$ 0.44 M$_\odot$. We report one BSS candidate, BSS6, as an Algol-type eclipsing binary with a period of 0.8006 days, based on the Gaia DR3 variability classification. We suggest that BSS3 is formed via either the Case A or Case B mass-transfer channel, whereas BSS6 is formed via the Case A mass transfer.

I have used high-precision photometry and astrometry from the third data release of Gaia (DR3) to identify candidate members of the 32 Ori association. Spectral types and radial velocities have been measured for subsets of the candidates using new and archival spectra. For the candidates that have radial velocity measurements, I have used UVW velocities to further constrain their membership, arriving at a final catalog of 169 candidates. I estimate that the completeness of the survey is ~90% for spectral types of <=M7 (>=0.06 Msun). The histogram of spectral types for the 32 Ori candidates exhibits a maximum at M5 (~0.15 Msun), resembling the distributions measured for other young clusters and associations in the solar neighborhood. The available UVW velocities indicate that the association is expanding, but they do not produce a well-defined kinematic age. Based on their sequences of low-mass stars in color-magnitude diagrams, the 32 Ori association and Upper Centaurus-Lupus/Lower Centaurus-Crux (UCL/LCC) are coeval to within +/-1.2 Myr, and they are younger than the Beta Pic moving group by ~3 Myr, which agrees with results from previous analysis based on the second data release of Gaia. Finally, I have used mid-infrared (IR) photometry from the Wide-field Infrared Survey Explorer to check for excess emission from circumstellar disks among the 32 Ori candidates. Disks are detected for 18 candidates, half of which are reported for the first time in this work. The fraction of candidates at <=M6 that have full, transitional, or evolved disks is 10/149=0.07+0.03/-0.02, which is consistent with the value for UCL/LCC.

We present a new pipeline developed to detect and characterize faint astronomical companions at small angular separation from the host star using sets of wide-field imaging observations not specifically designed for High Contrast Imaging analysis. The core of the pipeline relies on Karhunen-Lo\^eve truncated transformation of the reference PSF library to perform PSF subtraction and identify candidates. Tests of reliability of detections and characterization of companions are made through simulation of binaries and generation of Receiver Operating Characteristic curves for false positive/true positive analysis. The algorithm has been successfully tested on large HST/ACS and WFC3 datasets acquired for two HST Treasury Programs on the Orion Nebula Cluster. Based on these extensive numerical experiments we find that, despite being based on methods designed for observations of single star at a time, our pipeline performs very well on mosaic space based data. In fact, we are able to detect brown dwarf-mass companions almost down to the planetary mass limit. The pipeline is able to reliably detect signals at separations as close as $\gtrsim 0.1 "$ with a completeness of $\gtrsim 10\%$, or $\sim 0.2"$ with a completeness of $\sim 30\%$. This approach can potentially be applied to a wide variety of space based imaging surveys, starting with data in the existing HST archive, near-future JWST mosaics, and future wide-field Roman images.

Non-Maxwellian $\kappa$ electron energy distributions (EEDs) have been proposed in recent years to resolve the so-called ``electron temperature and abundance discrepancy problem'' in the study of planetary nebulae (PNe). Thus the need to develop diagnostic tools to determine from observations the EED of PNe is raised. Arising from high energy levels, the ultraviolet (UV) emission lines from PNe present intensities that depend sensitively on the high-energy tail of the EED. In this work, we investigate the feasibility of using the \ion{C}{2}]$\lambda$2326/\ion{C}{2}$\lambda$1335 intensity ratios as a diagnostic of the deviation of the EED from the Maxwellian distribution (as represented by the $\kappa$ index). We use a Maxwellian decomposition approach to derive the theoretical $\kappa$-EED-based collisionally excited coefficients of \ion{C}{2}, and then compute the \ion{C}{2} UV intensity ratio as a function of the $\kappa$ index. We analyze the archival spectra acquired by the {\it International Ultraviolet Explorer} and measure the intensities of \ion{C}{2} UV lines from 12 PNe. By comparing the observed line ratios and the theoretical predictions, we can infer their $\kappa$ values. With the Maxwellian-EED hypothesis, the observed \ion{C}{2}]$\lambda$2326/\ion{C}{2}$\lambda$1335 ratios are found to be generally lower than those predicted from the observed optical spectra. This discrepancy can be explained in terms of the $\kappa$ EED. Our results show that the $\kappa$ values inferred range from 15 to infinity, suggesting a mild or modest deviation from the Maxwellian distribution. However, the $\kappa$-distributed electrons are unlikely to exist throughout the whole nebulae. A toy model shows that if just about 1--5 percent of the free electrons in a PN had a $\kappa$-EED as small as $\kappa=3$, it would be sufficient to account for the observations.

Of all the solar fundamental parameters (mass, diameter, gravity at the surface,...), the gravitational moments have been quite often ignored in the past, mainly due to the great difficulty to get a reliable estimate. Even though the order of magnitude of the solar quadrupole moment $J_2$ is now known to be $10^{-7}$, its accurate value is still discussed. Indeed, the expansion in multipoles $J_{(l,~ l = 2, ...)}$ of the gravitational potential of a rotating body affects the orbital motion of planets at a relativistic level. We will recall here the recent progresses made in testing General Relativity through the contribution of the first solar quadrupole moment. Using the Eddington-Robertson parameters, we recall the constraints both on a theoretical and experimental point of view. Together with $\gamma$, which encodes the amount of curvature of space-time per unit rest-mass, the Post--Newtonian Parameter $\beta$ contributes to the relativistic precession of planets. The latter parameter encodes the amount of non-linearity in the superposition law of gravitation. Even though in principle, it would be possible to extract $J_2$ from planetary ephemerides, we observe that it is significantly correlated with other solution parameters (semi-major axis of planets, mass of asteroids...). Focusing on the $J_2$ correlations, we show that in general, when ~$\beta$ and ~$\gamma$ are freed, the correlations ~[$\beta, J_2$] and ~[$\gamma, J_2$] are $\approx$ 45\% and $\approx$ 55\% respectively. Moreover, all the planetary dynamics-based values are biased by the Lense--Thiring effect, which has never been modeled and solved for so far, but can be estimated to $\approx$ 7\%. It is thus possible to get a good estimate of the solar quadrupole moment:$1.66\times10^{-7}$$\leq$$J_2$$\leq$$2.32\times10^{-7}$.

In the framework of Randall -- Sundrum theory with extra dimension Reissner -- Nordstrom black hole solutions with a tidal charge have been found. The shadow around the supermassive black hole in M87 was reconstructed in 2019 based on observations with the Event Horizon Telescope (EHT) in April 2017. In May 2022 the EHT Collaboration presented results of a shadow reconstruction for our Galactic Center. Earlier, for Reissner -- Nordstr\"om metric we derived analytical expressions for shadow size as a function of charge and later generalized these results for a tidal charge case. We discuss opportunities to evaluate parameters of alternative theories of gravity with shadow size estimates done by the EHT Collaboration, in particular, a tidal charge could be estimated from these observations.

Dark matter is believed to constitute the majority of the matter content of the universe, but virtually nothing is known about its nature. Physical properties of a candidate particle can be probed via indirect detection by observing the decay and/or annihilation products. While this has previously been done primarily through gamma-ray studies, the increased sensitivity of new radio interferometers means that searches via the radio bandwidth are the new frontrunners. MeerKAT's high sensitivity, ranging from 3 $\mu$Jy beam$^{-1} $ for an 8 arcsecond beam to 10 $\mu$Jy beam$^{-1} $ for an 15 arcsecond beam, make it a prime candidate for radio dark matter searches. Using MeerKAT Galaxy Cluster Legacy Survey (MGCLS) data to obtain diffuse synchrotron emission within galaxy clusters, we are able to probe the properties of a dark matter model. In this work we consider both generic WIMP annihilation channels as well as the 2HDM+S model. The latter was developed to explain various anomalies observed in Large Hadron Collider (LHC) data from runs 1 and 2. The use of public MeerKAT data allows us to present the first WIMP dark matter constraints produced using this instrument.

Cosmological shock waves are essential to understanding the formation of cosmological structures. To study them, scientists run computationally expensive high-resolution 3D hydrodynamic simulations. Interpreting the simulation results is challenging because the resulting data sets are enormous, and the shock wave surfaces are hard to separate and classify due to their complex morphologies and multiple shock fronts intersecting. We introduce a novel pipeline, Virgo, combining physical motivation, scalability, and probabilistic robustness to tackle this unsolved unsupervised classification problem. To this end, we employ kernel principal component analysis with low-rank matrix approximations to denoise data sets of shocked particles and create labeled subsets. We perform supervised classification to recover full data resolution with stochastic variational deep kernel learning. We evaluate on three state-of-the-art data sets with varying complexity and achieve good results. The proposed pipeline runs automatically, has only a few hyperparameters, and performs well on all tested data sets. Our results are promising for large-scale applications, and we highlight now enabled future scientific work.

Early results of JWST observations have delivered bright $z\gtrsim 10$ galaxy candidates in greater numbers than expected, enabling construction of the rest-frame UV luminosity functions (LFs). The LFs contain key information on the galaxy assembly history, star formation activity, and stellar population in the distant universe. Given an upper bound of the total baryonic mass inflow rate to galaxies from their parent halos estimated from abundance matching, we derive a lower bound on the product of the star formation and UV-photon production efficiency in galaxies at each redshift. This stringent constraint requires a high efficiency ($\gtrsim 10-30\%$) converting gas into stars, assuming a normal stellar population with a Salpeter-like mass distribution. The efficiency is substantially higher than those of typical nearby galaxies, but is consistent with those seen in starburst galaxies and super star clusters observed in the nearby universe. Alternatively, the star formation efficiency may be as low as a few percent, which is the average value for the entire galaxy population at $z\simeq 6$, if the stellar population is metal-free and drawn from a top-heavy mass distribution that produces more intense UV radiation. We discuss several other possible scenarios to achieve the constraint, for instance, energetic radiation produced from compact stellar-remnants and quasars, and propose ways to distinguish the scenarios by forthcoming observations.

In this work, we investigate how the central stellar metallicity ([Z/H]) of 1363 galaxies from the SAMI galaxy survey is related to their stellar mass and a proxy for the gravitational potential, $\Phi$ = log10(M/M*) - log10($r_e$/kpc). In agreement with previous studies, we find that passive and star-forming galaxies occupy different areas of the [Z/H]-M* plane, with passive galaxies having higher [Z/H] than star-forming galaxies at fixed mass (a difference of 0.23 dex at log10(M/M*)=10.3). We show for the first time that all galaxies lie on the same relation between [Z/H] and $\Phi$, and show that the offset in [Z/H] between passive and star-forming galaxies at fixed $\Phi$ is smaller than or equal to the offset in [Z/H] at fixed mass (an average $\Delta$[Z/H] of 0.11 dex at fixed $\Phi$ compared to 0.21 dex at fixed mass). We then build a simple model of galaxy evolution to explain and understand our results. By assuming that [Z/H] traces $\Phi$ over cosmic time and that the probability that a galaxy quenches depends on both its mass and size, we are able to reproduce these offsets in stellar metallicity with a model containing instantaneous quenching. We therefore conclude that an offset in metallicity at fixed mass cannot by itself be used as evidence of slow quenching processes, in contrast to previous studies. Instead, our model implies that metal-rich galaxies have always been the smallest objects for their mass in a population. Our findings reiterate the need to consider galaxy size when studying stellar populations.

The partial spatial separation of cold dark matter (DM) and gas is a ubiquitous feature in the formation of cosmic large-scale structure. This separation, termed dissociation, is prominent in galaxy clusters that formed through collisions of massive progenitors, such as the famous `Bullet' cluster. A direct comparison of the incidence of such dissociated structures with theoretical predictions is challenged by the rarity of strongly dissociated systems and the difficulty to quantify dissociation. This paper introduces a well-defined dimensionless dissociation index $S\in[-1,1]$ that encodes the quadrupole difference between DM and gas in a custom region. Using a simulation of cosmic large-scale structure with cold DM and ideal non-radiating gas, in $\Lambda$CDM cosmology, we find that 90 per cent of the haloes are positively dissociated ($S>0$), meaning their DM is more elongated than their gas. The spatial density of highly dissociated massive structures appears consistent with observations. Through idealised $N$-body+SPH simulations of colliding gaseous DM haloes, we further explore the details of how ram-pressure causes dissociation in binary collisions. A suite of 300 such simulations reveals a scale-free relation between the orbital parameters of binary collisions and the resulting dissociation. Building on this relation, we conclude that the frequency of dissociated structures in non-radiative cosmological simulations is nearly fully accounted for by the major (mass ratio $>1:10$) binary collisions predicted by such simulations. In principle, our results allow us to constrain the orbital parameters that produced specific observed dissociated clusters.

Intrinsic alignments between galaxies and the large-scale structure contaminate galaxy clustering analyses and impact constraints on galaxy bias and the growth rate of structure in the Universe. This is the result of alignments inducing a selection effect on spectroscopic samples which is correlated with the large-scale structure. In this work, we quantify the biases on galaxy bias and the growth rate when alignments are neglected. We also examine different options for the mitigation of alignments by considering external priors on the effect and different probe combinations. We find that conservative analyses that restrict to $k_{\rm max}=0.1$ Mpc$^{-1}$ are not significantly affected. However, analyses that aim to go to higher wave numbers could evidence a significant contamination from alignments. In those cases, including a prior on alignment amplitude, or combining clustering with the position-intrinsic shape correlation of galaxies, can recover the same expected constraining power, or even inform bias and growth rate measurements.

Models of `modified-inertia' formulation of MOND are described and applied to nonrelativistic many-body systems. They involve time-nonlocal equations of motion. Momentum, angular momentum, and energy are (nonlocally) defined, whose total values are conserved for isolated systems. The models make all the salient MOND predictions. Yet, they differ from existing `modified-gravity' formulations in some second-tier predictions. The models describe correctly the motion of a composite body in a low-acceleration field even when the internal accelerations of its constituents are high. They exhibit a MOND external field effect (EFE) that shows some important differences from what we have come to expect from modified-gravity versions: In one, simple example of the models, what determines the EFE, in the case of a dominant external field, is $\mu(\theta\langle a_{ex}\rangle/a_0)$, where $\mu(x)$ is the MOND `interpolating function' that describes rotation curves, compared with $\mu(a_{ex}/a_0)$ for presently-known `modified-gravity' formulations. The two main differences are that while $a_{ex}$ is the momentary value of the external acceleration, $\langle a_{ex}\rangle$ is a certain time average of it, and that $\theta>1$ is an extra factor that depends on the frequency ratio of the external- and internal-field variations. Only ratios of frequencies enter, and $a_0$ remains the only new dimensioned constant. For a system on a circular orbit in a galaxy (such as the vertical dynamics in a disc galaxy), the first difference disappears, since $\langle a_{ex}\rangle=a_{ex}$. But the $\theta$ factor can appreciably enhance the EFE in quenching MOND effects, over what is deduced in modified gravity. Some exact solutions are also described, such as for rotation curves, for an harmonic force, and the general, two-body problem, which in the deep-MOND regime reduces to a single-body problem.

It is now clear that the stars in the Solar neighbourhood display large-scale coherent vertical breathing motions. At the same time, Milky Way-like galaxies experience tidal interactions with satellites/companions during their evolution. While these tidal interactions can excite vertical oscillations, it is still not clear whether vertical breathing motions are excited \textit{directly} by the tidal encounters or are driven by the tidally-induced spirals. We test whether excitation of breathing motions are directly linked to tidal interactions by constructing a set of $N$-body models (with mass ratio 5:1) of unbound, single fly-by interactions with varying orbital configurations. We first reproduce the well-known result that such fly-by interactions can excite strong transient spirals (lasting for $\sim 2.9-4.2$ Gyr) in the outer disc of the host galaxy. The generation and strength of the spirals are shown to vary with the orbital parameters (the angle of interaction, and the orbital spin vector). Furthermore, we demonstrate that our fly-by models exhibit coherent breathing motions whose amplitude increases with height. The amplitudes of breathing motions show characteristic modulation along the azimuthal direction, with compressing breathing motions coinciding with the peaks of the spirals and expanding breathing motions falling in the inter-arm regions -- a signature of a spiral-driven breathing motion. These breathing motions in our models end when the strong tidally-induced spiral arms fade away. Thus, it is the tidally-induced spirals which drive the large-scale breathing motions in our fly-by models, and the dynamical role of the tidal interaction in this context is indirect.

A ubiquitous presence of weak energy releases is one of the most promising hypotheses to explain coronal heating, referred to as the nanoflare hypothesis. The accelerated electrons associated with such weak heating events are also expected to give rise to coherent impulsive emission via plasma instabilities in the meterwave radio band, making this a promising spectral window to look for their presence. Recently \citet{Mondal2020b} reported the presence of weak impulsive emissions from quiet Sun regions which seem to meet the requirements of being radio counterparts of the hypothesized nanoflares. Detection of such low-contrast weak emission from the quiet Sun is challenging and, given their implications, it is important to confirm their presence. In this work, using data from the Murchison Widefield Array, we explore the use of an independent robust approach for their detection by separating the dominant slowly varying component of emission from the weak impulsive one in the visibility domain. We detect milli-SFU level bursts taking place all over the Sun and characterize their brightness temperatures, distributions, morphologies, durations and association with features seen in EUV images. We also attempt to constraint the energies of the nonthermal particles using inputs from the FORWARD coronal model along with some reasonable assumptions and find them to lie in the sub-pico flare ($\sim 10^{19}-10^{21}$ ergs) range. In the process, we also discover perhaps the weakest type III radio burst and another one that shows clear signatures of weakest quasi-periodic pulsations.

This chapter describes the history, principles, and recent developments of large field of view X-ray optics based on lobster eye designs. Most of grazing incidence (reflective) X-ray imaging systems used in astronomy and other applications, are based on the Wolter 1 (or modified) arrangement. But there are also other designs and configurations proposed for future applications for both laboratory and space environments. Kirkpatrick-Baez (K-B) based lenses as well as various types of lobster eye optics serve as an example. Analogously to Wolter lenses, all these systems use the principle that the X-rays are reflected twice to create focal images. Various future projects in X-ray astronomy and astrophysics will require large optics with wide fields of view. Both large Kirkpatrick-Baez modules and lobster eye X-ray telescopes may serve as solutions as these can offer innovations such as wide fields of view, low mass and reduced costs. The basic workings of lobster eye optics using Micro Pore Optics (MPOs) and their various uses are discussed. The issues and limiting factors of these optics are evaluated and current missions using lobster eye optics to fulfil their science objectives are reviewed. The Multi Foil Optics (MFO) approach represents a promising alternative. These arrangements can also be widely applied in laboratory devices. The chapter also examines the details of alternative applications for non-Wolter systems in other areas of science, where some of these systems have already demonstrated their advantages such as the K-B systems which have already found wide applications in laboratories and synchrotrons.

We present the first observational study of pulsars performed with the second-generation precursor stations to the low-frequency component of the Square Kilometre Array (SKA-Low): the Aperture Array Verification System 2 (AAVS2) and the Engineering Development Array 2 (EDA2). Using the SKA-Low stations, we have observed 100 southern-sky pulsars between 70-350 MHz, including follow-up observations at multiple frequencies for a selected sample of bright pulsars. These observations have yielded detections of 22 pulsars, including the lowest-frequency detections ever published for 6 pulsars, despite the modest sensitivity of initial system where the recording bandwidth is limited to ~1 MHz. By comparing simultaneous flux density measurements obtained with the SKA-Low stations and performing rigorous electromagnetic simulations, we verify the accuracy of the SKA-Low sensitivity simulation code presented in Sokolowski et al. (2022). Furthermore, we perform model fits to the radio spectra of the detected pulsars using the method developed by Jankowski et al. (2018), including 9 pulsars which were not fitted in the original work. We robustly classify the spectra into 5 morphological classes and find that all but one pulsar exhibit deviations from simple power-law behaviour. These findings suggest that pulsars with well-determined spectra are more likely to show spectral flattening or turn-over than average. Our work demonstrates how SKA-Low stations can be meaningfully used for scientifically useful measurements and analysis of pulsar radio spectra, which are important inputs for informing pulsar surveys and science planned with the SKA-Low.

The Etherington distance duality relation is well-established for metric theories of gravity, and confirms the duality between the luminosity distance and the angular diameter distance through the conservation of surface brightness. A violation of the Etherington distance duality due to lensing in a non-metric spacetime would lead to a fluctuations in surface brightness of galaxies. Likewise, fluctuations of the surface brightness can arise in classical astrophysics as a consequence of intrinsic tidal interaction of galaxies with their environment. Therefore, we study these in two cases in detail: Firstly, for intrinsic size fluctuations and the resulting changes in surface brightness, and secondly, for an area-metric spacetime as an example of a non-metric spacetime where the distance duality relation itself acquires modifications. We compare the auto- and cross-correlations of the angular spectra in these two cases and show that the fluctuations in intrinsic brightness can potentially be measured with a cumulative signal-to-noise ratio $\Sigma(\ell) \geq 3$ in an Euclid-like survey. The measurement in area-metric spacetimes, however, depends on the specific parameter choices, which also determine the shape and amplitude of the spectra. While lensing surveys do have sensitivity to lensing-induced surface brightness fluctuations in area-metric spacetimes, the measurement does not seem to be possible for natural values of the Etherington-breaking parameters.

Extragalactic Background Light (EBL) studies have revealed a significant discrepancy between direct measurements -- via instruments measuring bare sky from which Zodiacal and Galactic light models are subtracted -- and measurements of the Integrated Galaxy Light (IGL). This discrepancy could lie in either method, whether it be an incomplete Zodiacal model or missed faint galaxies in the IGL calculations. It has been proposed that the discrepancy is due to deep galaxy surveys, such as those with the Hubble Space Telescope, missing up to half of the faint galaxies with $24 \le m_{AB} \le 29$ mag. We address this possibility by creating three replications of the Hubble UltraDeep Field (HUDF) through successive rotations and additions to the original. SourceExtractor is used to analyze these replications to compare the recovered counts and photometry to the original HUDF, allowing us to assess how many galaxies were missed due to confusion, i.e., due to blending with neighboring faint galaxies. This exercise reveals that, while up to 50% of faint galaxies with $28 \le m_{AB} \le 29$ mag were missed or blended with neighboring objects in certain filters, not enough were missed to account for the EBL discrepancy alone in any of the replications.

Digital radio detection of cosmic-ray air showers has emerged as an alternative technique in high-energy astroparticle physics. Estimation of the detection efficiency of cosmic-ray radio arrays is one of the few remaining challenges regarding this technique. To address this problem, we developed a model based on the explicit probabilistic treatment of key elements of the radio technique for air showers: the footprint of the radio signal on ground, the detection of the signal in an individual antenna, and the detection criterion on the level of the entire array. The model allows for estimation of sky regions of full efficiency and can be used to compute the aperture of the array, which is essential to measure the absolute flux of cosmic rays. We also present a semi-analytical method that we apply to the generic model, to calculate the efficiency and aperture with high accuracy and reasonable calculation time. The model in this paper is applied to the Tunka-Rex array as example instrument and validated against Monte Carlo simulations. The validation shows that the model performs well, in particular, in the prediction of regions with full efficiency. It can thus be applied to other antenna arrays to facilitate the measurement of absolute cosmic-ray fluxes and to minimize a selection bias in cosmic-ray studies.

Using the convolution of seeing and diffraction, the relation between seeing and aperture in the visibility of sunspots is explored. It is shown that even telescopes with apertures smaller than 5 centimetres are significantly affected by seeing. Although larger aperture instruments suffer more from seeing than smaller ones, their level of detail always remains better under the same atmospheric conditions.

Near-Earth asteroid population models predict the existence of bodies located inside the orbit of Venus. Despite searches up to the end of 2019, none had been found. We report discovery and follow-up observations of (594913) 'Ayl\'o'chaxnim, an asteroid with an orbit entirely interior to Venus. (594913) 'Ayl\'o'chaxnim has an aphelion distance of ~0.65 au, is ~2 km in diameter and is red in colour. The detection of such a large asteroid inside the orbit of Venus is surprising given their rarity according to near-Earth asteroid population models. As the first officially numbered and named asteroid located entirely within the orbit of Venus, we propose that the class of interior to Venus asteroids be referred to as 'Ayl\'o'chaxnim asteroids.

Surface brightness -- colour relations (SBCRs) are very useful tools for predicting the angular diameters of stars. They offer the possibility to calculate very precise spectrophotometric distances by the eclipsing binary method or the Baade-Wesselink method. Double-lined Detached Eclipsing Binary stars (SB2 DEBs) with precisely known trigonometric parallaxes allow for a calibration of SBCRs with unprecedented precision. In order to improve such calibrations, it is important to enlarge the calibration sample of suitable eclipsing binaries with very precisely determined physical parameters. We carefully chose a sample of ten SB2 DEBs in the solar neighbourhood which contain inactive main-sequence components. The components have spectral types from early A to early K. All systems have high-precision parallaxes from the Gaia mission. We analysed high precision ground- and space-based photometry simultaneously with the radial velocity curves derived from HARPS spectra. We used spectral disentangling to obtain the individual spectra of the components and used these to derive precise atmospheric parameters and chemical abundances. For almost all components, we derived precise surface temperatures and metallicities. We derived absolute dimensions for 20 stars with an average precision of 0.2% and 0.5% for masses and radii, respectively. Three systems show slow apsidal motion. One system, HD 32129, is most likely a triple system with a much fainter K6V companion. Also three systems contain metallic-line components and show strong enhancements of barium and ittrium. The components of all systems compare well to the SBCR derived before from the detached eclipsing binary stars. With a possible exception of HD 32129, they can be used to calibrate SBCRs with a precision better than 1% with available Gaia DR3 parallaxes.

Context. The so-called "light echoes" of supernovae - the apparent motion of outburst-illuminated interstellar dust - can be detected in astronomical difference images; however, light echoes are extremely rare which makes manual detection an arduous task. Surveys for centuries-old supernova light echoes can involve hundreds of pointings of wide-field imagers wherein the subimages from each CCD amplifier require examination. Aims. We introduce ALED, a Python package that implements (i) a capsule network trained to automatically identify images with a high probability of containing at least one supernova light echo, and (ii) routing path visualization to localize light echoes and/or light echo-like features in the identified images. Methods. We compare the performance of the capsule network implemented in ALED (ALED-m) to several capsule and convolutional neural networks of different architectures. We also apply ALED to a large catalogue of astronomical difference images and manually inspect candidate light echo images for human verification. Results. ALED-m, was found to achieve 90% classification accuracy on the test set, and to precisely localize the identified light echoes via routing path visualization. From a set of 13,000+ astronomical images, ALED identified a set of light echoes that had been overlooked in manual classification. ALED is available via github.com/LightEchoDetection/ALED.

The observed presence of water molecules in the dayside lunar regolith was an unexpected discovery and remains poorly understood. Standard thermophysical models predict temperatures that are too high for adsorbed water to be stable. We propose that this problem can be caused by the assumption of locally flat surfaces that is common in such models. Here we apply a model that explicitly considers surface roughness, and accounts for solar illumination, shadows cast by topography, self-heating, thermal reradiation, and heat conduction. We couple the thermophysical model to a model of first-order desorption of lunar surface water and demonstrate that surface roughness substantially increases the capability of the Moon to retain water on its sunlit hemisphere at any latitude, and within $45^{\circ}$ of the poles, at any time of the lunar day. Hence, we show that lunar surface roughness has a strong influence on lunar water adsorption and desorption. Therefore, it is of critical importance to take account of surface roughness to get an accurate picture of the amount of water on the Moon's surface and in its exosphere.

The Hubble law (HL) governs the low-redshift (low-z) evolution of the distance of an object. However, there is a lack of an investigation of its validity and effective radius for a long time, since the low-z background data with a high precision is scarce. The latest Type Ia supernovae sample Pantheon+ having a significant increase of low-z data provides an excellent opportunity to test the HL. We propose a generalized HL and implement the first modern test of the HL with Pantheon+. We obtain the constraint on the deviation parameter $\alpha=1.00118\pm0.00044$, confirm the validity of linear HL with a $0.04\%$ precision and give the transition redshift $z_t=0.03$ and luminosity distance $D_{L,t}=123.13\pm1.75$ Mpc, which means that HL holds when $z<0.03$ and breaks down at a distance of $D_L>123.13$ Mpc. Comparing the ability of Type Ia supernovae and HII galaxies in testing the HL, we stress the uniqueness and strong power of Type Ia supernovae in probing the low-z physics.

Context. Despite being a prominent subset of the exoplanet population discovered in the past three decades, the nature and provenance of sub-Neptune-sized planets are still one of the open questions in exoplanet science. Aims. For planets orbiting bright stars, precisely measuring the orbital and planet parameters of the system is the best approach to distinguish between competing theories regarding their formation and evolution. Methods. We obtained 69 new radial velocity observations of the mid-M dwarf G 9-40 with the CARMENES instrument to measure for the first time the mass of its transiting sub-Neptune planet, G 9-40 b, discovered in data from the K2 mission. Results. Combined with new observations from the TESS mission during Sectors 44, 45, and 46, we are able to measure the radius of the planet to an uncertainty of 3.4% (Rb = 1.900 +- 0.065 Re) and determine its mass with a precision of 16% (Mb = 4.00 +- 0.63 Me). The resulting bulk density of the planet is inconsistent with a terrestrial composition and suggests the presence of either a water-rich core or a significant hydrogen-rich envelope. Conclusions. G 9-40 b is referred to as a keystone planet due to its location in period-radius space within the radius valley. Several theories offer explanations for the origin and properties of this population and this planet is a valuable target for testing the dependence of those models on stellar host mass. By virtue of its brightness and small size of the host, it joins L 98-59 d as one of the two best warm (Teq ~ 400 K) sub-Neptunes for atmospheric characterization with JWST, which will probe cloud formation in sub-Neptune-sized planets and break the degeneracies of internal composition models.

We present a status update for MagAO-X, a 2000 actuator, 3.6 kHz adaptive optics and coronagraph system for the Magellan Clay 6.5 m telescope. MagAO-X is optimized for high contrast imaging at visible wavelengths. Our primary science goals are detection and characterization of Solar System-like exoplanets, ranging from very young, still-accreting planets detected at H-alpha, to older temperate planets which will be characterized using reflected starlight. First light was in Dec, 2019, but subsequent commissioning runs were canceled due to COVID-19. In the interim, MagAO-X has served as a lab testbed. Highlights include implementation of several focal plane and low-order wavefront sensing algorithms, development of a new predictive control algorithm, and the addition of an IFU module. MagAO-X also serves as the AO system for the Giant Magellan Telescope High Contrast Adaptive Optics Testbed. We will provide an overview of these projects, and report the results of our commissioning and science run in April, 2022. Finally, we will present the status of a comprehensive upgrade to MagAO-X to enable extreme-contrast characterization of exoplanets in reflected light. These upgrades include a new post-AO 1000-actuator deformable mirror inside the coronagraph, latest generation sCMOS detectors for wavefront sensing, optimized PIAACMC coronagraphs, and computing system upgrades. When these Phase II upgrades are complete we plan to conduct a survey of nearby exoplanets in reflected light.

We present the conceptual design of GMagAO-X, an extreme adaptive optics system for the 25 m Giant Magellan Telescope (GMT). We are developing GMagAO-X to be available at or shortly after first-light of the GMT, to enable early high contrast exoplanet science in response to the Astro2020 recommendations. A key science goal is the characterization of nearby potentially habitable terrestrial worlds. GMagAO-Xis a woofer-tweeter system, with integrated segment phasing control. The tweeter is a 21,000 actuator segmented deformable mirror, composed of seven 3000 actuator segments. A multi-stage wavefront sensing system provides for bootstrapping, phasing, and high order sensing. The entire instrument is mounted in a rotator to provide gravity invariance. After the main AO system, visible (g to y) and near-IR (Y to H) science channels contain integrated coronagraphic wavefront control systems. The fully corrected and, optionally, coronagraphically filtered beams will then be fed to a suite of focal plane instrumentation including imagers and spectrographs. This will include existing facility instruments at GMT via fiber feeds. To assess the design we have developed an end-to-end frequency-domain modeling framework for assessing the performance of GMagAO-X. The dynamics of the many closed-loop feedback control systems are then modeled. Finally, we employ a frequency-domain model of post-processing algorithms to analyze the final post-processed sensitivity. The CoDR for GMagAO-X was held in September, 2021. Here we present an overview of the science cases, instrument design, expected performance, and concept of operations for GMagAO-X.

GMagAO-X is the near first light ExAO coronagraphic instrument for the 25.4m GMT. It is designed for a slot on the folded port of the GMT. To meet the strict ExAO fitting and servo error requirement (<90nm rms WFE), GMagAO-X must have 21,000 actuator DM capable of >2KHz correction speeds. To minimize wavefront/segment piston error GMagAO-X has an interferometric beam combiner on a vibration isolated table, as part of this "21,000 actuator parallel DM". Piston errors are sensed by a Holographic Dispersed Fringe Sensor (HDFS). In addition to a coronagraph, it has a post-coronagraphic Lyot Low Order WFS (LLOWFS) to sense non-common path (NCP) errors. The LLOWFS drives a non-common path DM (NCP DM) to correct those NCP errors. GMagAO-X obtains high-contrast science and wavefront sensing in the visible and/or the NIR. Here we present our successful externally reviewed (Sept. 2021) CoDR optical-mechanical design that satisfies GMagAO-X's top-level science requirements and is compliant with the GMT instrument requirements and only requires COTS parts.

We present the discovery of two exoplanets transiting TOI-836 (TIC 440887364) using data from TESS Sector 11 and Sector 38. TOI-836 is a bright ($T = 8.5$ mag), high proper motion ($\sim\,200$ mas yr$^{-1}$), low metallicity ([Fe/H]$\approx\,-0.28$) K-dwarf with a mass of $0.68\pm0.05$ M$_{\odot}$ and a radius of $0.67\pm0.01$ R$_{\odot}$. We obtain photometric follow-up observations with a variety of facilities, and we use these data-sets to determine that the inner planet, TOI-836 b, is a $1.70\pm0.07$ R$_{\oplus}$ super-Earth in a 3.82 day orbit, placing it directly within the so-called 'radius valley'. The outer planet, TOI-836 c, is a $2.59\pm0.09$ R$_{\oplus}$ mini-Neptune in an 8.60 day orbit. Radial velocity measurements reveal that TOI-836 b has a mass of $4.5\pm0.9$ M$_{\oplus}$ , while TOI-836 c has a mass of $9.6\pm2.6$ M$_{\oplus}$. Photometric observations show Transit Timing Variations (TTVs) on the order of 20 minutes for TOI-836 c, although there are no detectable TTVs for TOI-836 b. The TTVs of planet TOI-836 c may be caused by an undetected exterior planet.

The Extremely Large Telescopes will require hundreds of actuators across the pupil for high Strehl in the visible. We envision a triple-stage AO (TSAO) system for GMT/GMagAO-X to achieve this. The first stage is a 4K DM controlled by an IR pyramid wavefront sensor that provides the first order correction. The second stage contains the high-order parallel DM of GMagAO-X that has 21000 actuators and contains an interferometric delay line for phasing of each mirror segment. This stage uses a Zernike wavefront sensor for high-order modes and a Holographic Dispersed Fringe Sensor for segment piston control. Finally, the third stage uses a dedicated 3K dm for non-common path aberration control and the coronagraphic wavefront control by using focal plane wavefront sensing and control. The triple stage architecture has been chosen to create simpler decoupled control loops. This work describes the performance of the proposed triple-stage AO architecture for ExAO with GMagAO-X.

Future Extremely Large Telescopes (ELTs) will require advances in Adaptive Optics (AO) systems to fully realize their potential. In addition to separate, dedicated wavefront sensors, it is recognized that wavefront sensing within the science focal plane itself will also be needed for many new instruments. One approach is to use On-Detector Guide Windows (ODGWs), whereby a small sub-window of a science detector is read-out continuously (~10s-100s of Hz) in parallel with slower reads of the full chip (>10 s). Guide star centroids from these windows can be used to correct for vibrations and flexure. Another potential use for these windows is to perform localized resets at high cadence to prevent saturation and to minimize persistence from bright sources. We have prototyped an ODGW system using a 5-um cutoff Teledyne HAWAII-2RG infrared detector, and the new Astronomical Research Cameras Gen-4 controller. We describe our implementation of an ODGW mode, and science image artifacts that were observed.

The search for exoplanets is pushing adaptive optics systems on ground-based telescopes to their limits. Currently, we are limited by two sources of noise: the temporal control error and non-common path aberrations. First, the temporal control error of the AO system leads to a strong residual halo. This halo can be reduced by applying predictive control. We will show and described the performance of predictive control with the 2K BMC DM in MagAO-X. After reducing the temporal control error, we can target non-common path wavefront aberrations. During the past year, we have developed a new model-free focal-plane wavefront control technique that can reach deep contrast (<1e-7 at 5 $\lambda$/D) on MagAO-X. We will describe the performance and discuss the on-sky implementation details and how this will push MagAO-X towards imaging planets in reflected light. The new data-driven predictive controller and the focal plane wavefront controller will be tested on-sky in April 2022.

The 25.4m Giant Magellan Telescope (GMT) will be amongst the first in a new series of segmented extremely large telescopes (ELTs). The 25.4 m pupil is segmented into seven 8.4 m circular segments in a flower petal pattern. At the University of Arizona we have developed a novel pupil slicer that will be used for ELT extreme adaptive optics (ExAO) on the up and coming ExAO instrument, GMagAO-X. This comes in the form of a six-sided reflective pyramid with a hole through the center known as a "hexpyramid". By passing the GMT pupil onto this reflective optic, the six outer petals will be sent outward in six different directions while the central segment passes through the center. Each segment will travel to its own polarization independent flat fold mirror mounted on a piezoelectric piston/tip/tilt controller then onto its own commercial 3,000 actuator deformable mirror (DM) that will be employed for extreme wavefront control. This scheme of seven DMs working in parallel to produce a 21,000 actuator DM is a new ExAO architecture that we named a "parallel DM," in which the hexpyramid is a key optical component. This significantly surpasses any current or near future actuator count for any monolithic DM architecture. The optical system is designed for high-quality wavefront (lambda/10 surface PV) with no polarization errors and no vignetting. The design and fabrication of the invar mechanical mounting structure for this complex optical system is described in this paper.

To take advantage of high-resolution optics sensitive to a broad energy range, future X-ray imaging instruments will require thick detectors with small pixels. This pixel aspect ratio affects spectral response in the soft X-ray band, vital for many science goals, as charge produced by the photon interaction near the entrance window diffuses across multiple pixels by the time it is collected, and is potentially lost below the imposed noise threshold. In an effort to understand these subtle but significant effects and inform the design and requirements of future detectors, we present simulations of charge diffusion using a variety of detector characteristics and operational settings, assessing spectral response at a range of X-ray energies. We validate the simulations by comparing the performance to that of real CCD detectors tested in the lab and deployed in space, spanning a range of thickness, pixel size, and other characteristics. The simulations show that while larger pixels, higher bias voltage, and optimal backside passivation improve performance, reducing the readout noise has a dominant effect in all cases. We finally show how high-pixel-aspect-ratio devices present challenges for measuring the backside passivation performance due to the magnitude of other processes that degrade spectral response, and present a method for utilizing the simulations to qualitatively assess this performance. Since compelling science requirements often compete technically with each other (high spatial resolution, soft X-ray response, hard X-ray response), these results can be used to find the proper balance for a future high-spatial-resolution X-ray instrument.

MagAO-X is an extreme adaptive optics (ExAO) instrument for the Magellan Clay 6.5-meter telescope at Las Campanas Observatory in Chile. Its high spatial and temporal resolution can produce data rates of 1 TB/hr or more, including all AO system telemetry and science images. We describe the tools and architecture we use for commanding, telemetry, and science data transmission and storage. The high data volumes require a distributed approach to data processing, and we have developed a pipeline that can scale from a single laptop to dozens of HPC nodes. The same codebase can then be used for both quick-look functionality at the telescope and for post-processing. We present the software and infrastructure we have developed for ExAO data post-processing, and illustrate their use with recently acquired direct-imaging data.

We promote the teaching of mass functions as an integral part of an interdisciplinary science education. Mass functions characterize the frequency distributions of objects with different masses on all cosmic scales. We intend to enhance experiential learning of this concept with a creative LEGO brick experiment for a diverse student audience. To our surprise, the LEGO mass function is not only qualitatively but also quantitatively comparable to mass functions found across the Universe. We also discuss the relation between gravitation and mass distributions as a possible explanation for the continuity of the universal mass function.

We study the impact of the muon pair production and double pair production processes induced by ultra-high energy photons on the cosmic microwave background. Although the muon pair production cross section is smaller than the electron pair production one, the associated energy loss length is comparable or shorter than the latter (followed by inverse Compton in the deep Klein-Nishina regime) at high-redshift, where the effect of the astrophysical radio background is expected to be negligible. By performing a simulation taking into account the details of $e/\gamma$ interactions at high energies, we show that a significant fraction of the electromagnetic energy injected at $E\gtrsim 10^{19}\,$eV at redshift $z\gtrsim 5$ is channeled into neutrinos. The double pair production plays a crucial role in enhancing the multiplicity of muon production in these electromagnetic cascades. The ultra-high energy neutrino spectrum, yet to be detected, can in principle harbour information on ultra-high energy sources in the young universe, either conventional or exotic ones, with weaker constraints from the diffuse gamma ray flux compared to their low redshift counterparts.

NASA's Gateway to Astronaut Photography of Earth Contains over 30000 photos of 2500 cataloged urban lightscapes (urban night lights) taken from the International Space Station. Over 100 of these multispectral DSLR photos are of sufficient spatial resolution, sharpness and exposure to be used for broadband spectral characterization of urban lightscapes. Analysis of simulated atmospheric transmissivity from the MODTRAN radiative transfer model shows that spectral slopes of transmissivity spectra are relatively insensitive to choice of model atmopshere, with variations in atmospheric path length and aerosol optical depth primarily affecting the bias of the spectrum rather than the slope. This suggests that color tempterature-calibrated RGB channels can be corrected for relative differences in atmospheric scattering and absorption to allow for quantitative intercomparison. A mosaic of 18 intercalibrated RGB photos renders a spectral feature space with four clearly defined spectral endmembers corresponding to white, yellow and red light sources, with brightness modulated by a dark background endmember. These four spectral endmembers form the basis of a linear spectral mixture model which can be inverted to provide estimates of the areal fraction of each endmember present within every pixel instantaneous field of view. The resulting spectral feature space shows two distinct mixing trends extending from the dark endmember to near flat spectrum (white-yellow) and warm spectrum (orange) sources. The distribution of illuminated pixels is strongly skewed toward a lower luminance background of warm spectrum street lighting, with brighter lights generally corresponding to point sources and major thoroughfares. Intercomparison of 18 individual urban lightscape spectral feature spaces show consistent topology, despite variations in exposure and interior mixing trends.

We challenge the traditional wisdom that cosmological (big bang relic) neutrinos can only be hot Dark Matter. We provide a critical review of the concepts, derivations and arguments in foundational books and recent publications that led respected researchers to proclaim that "[Dark Matter] cannot be neutrinos". We then provide the physics resulting in relic neutrino's significant power loss from the interaction of its anomalous magnetic moment with a high-intensity primordial magnetic fields, resulting in subsequent condensation into Condensed Neutrino Objects (CNOs). Finally, the experimental degenerate mass bounds that would rule out condensed cosmological neutrinos as the Dark Matter (unless there is new physics that would require a modification to the CNO Equation of State) are provided. We conclude with a discussion on new directions for research.

For decades, agronomists have used remote sensing to monitor key crop parameters like biomass, fractional cover, and plant health. Vegetation indices (VIs) are popular for this purpose, primarily leveraging the spectral red edge in multispectral imagery. In contrast, spectral mixture models use the full reflectance spectrum to simultaneously estimate area fractions of multiple endmember materials present within a mixed pixel. Here, we characterize the relationships between hyperspectral endmember fractions and 6 common multispectral VIs in crops & soils of California agriculture. Fractional area of green vegetation (Fv) was estimated directly from 64,000,000 5 nm, 3 to 5 m reflectance spectra compiled from a mosaic of 15 AVIRIS-ng flightlines. Simulated Planet SuperDove reflectance spectra were then derived from the AVIRIS-ng, and used to compute 6 popular VIs (NDVI, NIRv, EVI, EVI2, SR, DVI). Multispectral VIs were compared to hyperspectral Fv using parametric (Pearson correlation, r) and nonparametric (Mutual Information, MI) similarity metrics. 4 VIs (NIRv, DVI, EVI, EVI2) showed strong linear relationships to Fv (r > 0.94; MI > 1.2). NIRv & DVI showed strong interrelation (r > 0.99, MI > 2.4), but deviated significantly from 1:1 relative to Fv. EVI & EVI2 were also strongly interrelated (r > 0.99, MI > 2.3) and more closely followed a 1:1 relation with Fv. In contrast, NDVI & SR showed weaker, nonlinear, heteroskedastic relation to Fv (r < 0.84, MI = 0.69). NDVI showed especially severe sensitivity to substrate background reflectance (-0.05 < NDVI < +0.6 for unvegetated spectra) and saturation (0.2 < Fv < 0.8 for NDVI = 0.7). These direct observational constraints on multispectral VI and hyperspectral mixture model comparability can serve as a quantitative benchmark for agronomic applications in the coming era of increasing spatial & spectral resolution Earth observation.

We propose using trapped electrons as high-$Q$ resonators for detecting meV dark photon dark matter. When the rest energy of the dark photon matches the energy splitting of the two lowest cyclotron levels, the first excited state of the electron cyclotron will be resonantly excited. A proof-of-principle measurement, carried out with one electron, demonstrates that the method is background-free over a 7.4 day search. It sets a limit on dark photon dark matter at 148 GHz (0.6 meV) that is around 75 times better than previous constraints. Dark photon dark matter in the 0.1-1 meV mass range (20-200 GHz) could likely be detected at a similar sensitivity in an apparatus designed for dark photon detection.

The phenomena of neutrino spin flavour precession in the presence of an extraneous magnetic field is a repercussion of neutrino magnetic moment which is consociated with the physics beyond the standard model of electroweak interactions. Ultra high energy neutrinos are spawned from a number of sources in the universe including the highly energetic astrophysical objects such as active galactic nuclei, blazar or supermassive black holes. When such high energy neutrinos pass through any compact stellar objects like neutron stars or white dwarfs, their flux can significantly reduce due to the exorbitant magnetic field provided by these compact objects. For Dirac neutrinos, such phenomena occur due to the conversion of neutrinos to their sterile counterparts. In this work, we consider a neutron star possessing a spatially varying magnetic field which may or may not decay with time. We find that, for the non-decaying magnetic field, the flux of high energy Dirac neutrinos becomes nearly half after passing through the neutron star. The flux is further enfeebled by $\sim 10\%$ in the presence of muons inside the neutron star. For decaying magnetic field, the flux reduction is abated by $\sim 5\%$ as compared to the temporally static magnetic field. In the case of a white dwarf, the depletion of flux is lesser as compared to the neutron stars.

In this paper the influence of large-scale decreasing and increasing gradients of the density of magnetized plasma on the relaxation process of a continuously injected relativistic electron beam with an energy of 611 keV ($v_b=0.9c$) and a pitch-angle distribution is studied using particle-in-cell numerical simulations. It is found that for the selected parameters in the case of a smoothly decreasing gradient and in a homogeneous plasma the formation of spatially limited plasma oscillations of large amplitude occurs. In such cases, modulation instability develops and a long-wave longitudinal modulation of the ion density is formed. In addition, the large amplitude of plasma waves accelerates plasma electrons to energies on the order of the beam energy. In the case of increasing and sharply decreasing gradients, a significant decrease in the amplitude of plasma oscillations and the formation of a turbulent ion density spectrum are observed. The possibility of acceleration of beam electrons to energies more than 2 times higher than the initial energy of the beam particles is also demonstrated. This process takes place not only during beam propagation in growing plasma density, but also in homogeneous plasma due to interaction of beam particles with plasma oscillations of large amplitude.

We revisit the constraint from the recently reported cosmic birefringence on axion-like particles with a general decay constant. A particular attention is paid to the naturalness of the model parameter space, which has been overlooked in the literature. We show that the observed cosmic birefringence is naturally explained by the electroweak axion with a string-theory inspired decay constant $F_A\simeq 10^{16}$ GeV.

If the neutrino has a large magnetic moment, then a phase effect may appear in its spin-flavor precession inside the supernova. It differs from the ordinary flavor oscillation phase effect in two aspects: It can develop even if there is only one partially adiabatic resonance and it effects a large part of the neutrino energy spectrum. We examine the spin-flavor precession phase effect both analytically and numerically for the Majorana neutrinos in a core collapse supernova. Our analytical approach is based on the assumption that spin-flavor precession and Mikheev-Smirnov-Wolfenstein resonances are completely decoupled. Where this decoupling assumption fails, we present only numerical results. We show that the sensitive phase dependence of the survival probabilities can be treated as an uncertainty which smears the final neutrino energy spectra to be observed at Earth.

The concept of dark energy can be a candidate for preventing the gravitational collapse of compact objects to singularities. According to the usefulness of gravity's rainbow in UV completion of general relativity (by providing a new description of spacetime), it can be an excellent option to study the behavior of compact objects near phase transition regions. In this work, we obtain a modified Tolman-Openheimer-Volkof (TOV) equation for anisotropic dark energy as a fluid by solving the field equations in gravity's rainbow. Next, to compare the results with general relativity, we use a generalized Tolman-Matese-Whitman mass function to determine the physical quantities such as energy density, radial pressure, transverse pressure, gravity profile, and anisotropy factor of the dark energy star. We evaluate the junction condition and investigate the dynamical stability of dark energy star thin shell in gravity's rainbow. We also study the energy conditions for the interior region of this star. We show that the coefficients of gravity's rainbow can significantly affect this non-singular compact object and modify the model near the phase transition region.

The cosmological horizon problem is the question why spatial domains that were never in causal contact with each other now appear in precise symphony. We propose a solution to the horizon problem in which a globally synchronized early state is reached as the $\omega$-limit point of a transient, inhomogeneous Mixmaster universe. We show that the $\alpha$-limit set of the latter is a Kasner circle which represents a synchronized initial state of minimal entropy. Accordingly, unless the evolution is disrupted by quantum gravitational effects so that the initial state is not attained, Planck size domains emerge as causally disconnected, albeit in complete synchrony, as the universe enters an `isotropic' state to remain so in future unison.

In the first-order phase transitions (PTs) colliding bubble is an important gravitational wave (GW) source. Following bubble collision, domain walls can be formed when degenerate vacua occur as a result of the breaking of a discrete symmetry relevant to new physics at electroweak or higher scales. Using lattice simulations, we study the dynamical evolution of domain walls and find that the networks of the domain wall are formed around the completion of PTs and the lifetime of the wall networks largely depends on whether or not the degeneracy of true vacua is broken. Our numerical results indicate that domain wall networks continue to produce GWs in the aftermath of PTs, leading to dramatically changing the spectral shape and enhancing the magnitude by about one order. The resulting GW power spectra are peaked at $kR_* \simeq \pi$, above the peak wavenumber it has a decaying power law close to $k^{-1.2}$ followed by a slowly decreasing plateau with the UV cutoff at $kR_* \sim \mathcal{O}(10^2)$

Pulsar timing arrays (PTAs) detect gravitational waves (GWs) via the correlations that the waves induce in the arrival times of pulses from different pulsars. The mean correlation $\mu_{\rm u}(\gamma)$ as a function of the angle $\gamma$ between the directions to two pulsars was calculated by Hellings and Downs in 1983. The variance $\sigma^2_{\rm tot}(\gamma)$ in this correlation was recently calculated for a single pulsar pair at angle $\gamma$. Averaging over many such pairs, uniformly distributed on the sky, reduces this to an intrinsic cosmic variance $\sigma^2_{\rm cos}(\gamma)$. We extend that analysis to an arbitrary finite set of pulsars, distributed at specific sky locations, for which the pulsar pairs are grouped into finite-width bins in $\gamma$. Given (measurements or calculations of) the correlations for any set of pulsars, we find the best way to estimate the mean in each bin. The optimal estimator of the correlation takes into account correlations among all of the pulsars that contribute to that angular bin. We also compute the variance in the binned estimate. For narrow bins, as the number of pulsar pairs grows, the variance drops to the cosmic variance. For wider bins, by sacrificing angular resolution in $\gamma$, the variance can even be reduced below the cosmic variance. Our calculations assume that the GW signals are described by a Gaussian ensemble, which provides a good description of the confusion noise produced by expected PTA sources. We illustrate our methods with plots of the GW variance for the sets of pulsars currently monitored by several PTA collaborations. The methods can also be applied to future PTAs, where the improved telescopes will provide larger pulsar populations and higher-precision timing.