Papers for Wednesday, Nov 16 2022

Stars collect most of their mass during the protostellar stage, yet the accretion luminosity and stellar parameters, which are needed to compute the mass accretion rate, are poorly constrained for the youngest sources. The aim of this work is to fill this gap, computing the stellar properties and the accretion rates for a large sample of Class I protostars located in nearby (< 500 pc) star-forming regions and analysing their interplay. We used a self-consistent method to provide accretion and stellar parameters using SED modeling and veiling information from near-IR observations, when possible. We calculated accretion and stellar properties for the first time for 50 young stars. We focused our analysis on the 39 confirmed protostars, finding that their mass accretion rate varies between about 10^(-8) and about 10^(-4) Msun/yr in a stellar mass range between about 0.1 and 3 Msun. We find systematically larger mass accretion rates for our Class I sample than in Class II objects. Although the mass accretion rate we found is high, it still suggests that either stars collect most of its mass before Class I stage, or eruptive accretion is needed during the overall protostellar phase. Indeed, our results suggest that for a large number of protostars the disk can be unstable, which can result in accretion bursts and disk fragmentation in the past or in the future.

Enhanced ionizing radiation in close proximity to redshift $z\gtrsim 6$ quasars creates short windows of intergalactic Ly$\alpha$ transmission blueward of the quasar Ly$\alpha$ emission lines. The majority of these Ly$\alpha$ near-zones are consistent with quasars that have optically/UV bright lifetimes of $t_{\rm Q}\sim 10^{5}-10^{7}\rm\,yr$. However, lifetimes as short as $t_{\rm Q}\lesssim 10^{4}\rm\,yr$ appear to be required by the smallest Ly$\alpha$ near-zones. These short lifetimes present an apparent challenge for the growth of $\sim 10^{9}\rm\,M_{\odot}$ black holes at $z\gtrsim 6$. Accretion over longer timescales is only possible if black holes grow primarily in an obscured phase, or if the quasars are variable on timescales comparable to the equilibriation time for ionized hydrogen. Distinguishing between very young quasars and older quasars that have experienced episodic accretion with Ly$\alpha$ absorption alone is challenging, however. We therefore predict the signature of proximate 21-cm absorption around $z\gtrsim 6$ radio-loud quasars. For modest pre-heating of intergalactic hydrogen by the X-ray background, where the spin temperature $T_{\rm S} \lesssim 10^{2}\rm\,K$ prior to any quasar heating, we find proximate 21-cm absorption should be observable in the spectra of radio-loud quasars. The extent of the proximate 21-cm absorption is sensitive to the integrated lifetime of the quasar. Evidence for proximate 21-cm absorption from the diffuse intergalactic medium within $2-3\rm\,pMpc$ of a (radio-loud) quasar would be consistent with a short quasar lifetime, $t_{\rm Q}\lesssim 10^{5}\rm\,yr$, and would provide a complementary constraint on models for high redshift black hole growth.

The local distance ladder estimate of the Hubble constant ($H_0$) is important in cosmology, given the recent tension with the early universe inference. We estimate $H_0$ from the Type Ia supernova (SN Ia) distance ladder, inferring SN Ia distances with the hierarchical Bayesian SED model, BayeSN. This method has a notable advantage of being able to continuously model the optical and near-infrared (NIR) SN Ia light curves simultaneously. We use two independent distance indicators, Cepheids or the tip of the red giant branch (TRGB), to calibrate a Hubble-flow sample of 67 SNe Ia with optical and NIR data. We estimate $H_0 = 74.82 \pm 0.97$ (stat) $\pm\, 0.84$ (sys) km s$^{-1}$ Mpc$^{-1}$ when using the calibration with Cepheid distances to 37 host galaxies of 41 SNe Ia, and $70.92 \pm 1.14$ (stat) $\pm\,1.49$ (sys) km s$^{-1}$ Mpc$^{-1}$ when using the calibration with TRGB distances to 15 host galaxies of 18 SNe Ia. For both methods, we find a low intrinsic scatter $\sigma_{\rm int} \lesssim 0.1$ mag. We test various selection criteria and do not find significant shifts in the estimate of $H_0$. Simultaneous modelling of the optical and NIR yields up to $\sim$15% reduction in $H_0$ uncertainty compared to the equivalent optical-only cases. With improvements expected in other rungs of the distance ladder, leveraging joint optical-NIR SN Ia data can be critical to reducing the $H_0$ error budget.

We explore three sets of cosmological hydrodynamical simulations, IllustrisTNG, EAGLE, and SIMBA, to investigate the physical processes impacting the distribution of baryons in and around haloes across an unprecedented mass range of $10^8<M_{\rm 200c}/{\rm M_{\odot}}<10^{15}$, from the halo centre out to scales as large as $30\,R_{\rm 200c}$. We demonstrate that baryonic feedback mechanisms significantly redistribute gas, lowering the baryon fractions inside haloes while simultaneously accumulating this material outside the virial radius. To understand this large-scale baryonic redistribution and identify the dominant physical processes responsible, we examine several variants of TNG that selectively exclude stellar and AGN feedback, cooling, and radiation. We find that heating from the UV background in low-mass haloes, stellar feedback in intermediate-mass haloes, and AGN feedback in groups ($10^{12} \leq M_{\rm 200c}/{\rm M_{\odot}}<10^{14}$) are the dominant processes. Galaxy clusters are the least influenced by these processes on large scales. We introduce a new halo mass-dependent characteristic scale, the closure radius $R_{\rm c}$, within which all baryons associated with haloes are found. For groups and clusters, we introduce a universal relation between this scale and the halo baryon fraction: $R_{\rm c}/R_{\rm 200c,500c}-1=\beta(z)(1-f_{\rm b}(<R_{\rm 200c,500c})/f_{\rm b,cosmic})$, where $\beta(z)=\alpha\,(1+z)^\gamma$, and $\alpha$ and $\gamma$ are free parameters fit using the simulations. Accordingly, we predict that all baryons associated with observed X-ray haloes can be found within $R_{\rm c}\sim 1.5-2.5 R_{\rm 200c}$. Our results can be used to constrain theoretical models, particularly the physics of supernova and AGN feedback, as well as their interplay with environmental processes, through comparison with current and future X-ray and SZ observations.

We present the first large-scale 3D map of interstellar dust temperature. We build upon existing 3D reddening maps derived from starlight absorption (Bayestar19), covering 3/4 of the sky. Starting with the column density for each of 500 million 3D voxels, we propose a temperature and emissivity power-law slope ($\beta$) for each of them, and integrate along the line of sight to synthesize an emission map in five frequency bands observed by \emph{Planck} and \emph{IRAS}. The reconstructed emission map is constrained to match observations on a $10'$ scale, and does so with good fidelity. We produce 3D temperature maps at resolutions of $110', 55', $and $27'$. We assess performance on Cepheus, a dust cloud with two distinct components along the line of sight, and find distinct temperatures for the two components. We thus show that this methodology has enough precision to constrain clouds with different temperature along the line of sight up to $1-\sigma$ error. This would be an important result for dust frequency decorrelation foreground analysis for cosmic microwave background experiments, which would be impacted by a line-of-sight with varying temperature and magnetic field components. In addition to $T$ and $\beta$, we constrain the conversion factor between emission optical depth and reddening. This conversion factor is assumed to be constant in commonly used emission-based reddening maps. However, this work shows a factor of two variation that may prove significant for some applications.

We show that the 216.8$\pm$0.1 ms periodicity reported for the fast radio burst (FRB) 20191221A is very constraining for burst models. The high accuracy of burst periodicity (better than one part in 10$^3$), and the 2\% duty cycle (ratio of burst-duration and inter-burst interval), suggest a pulsar-like rotating beam model for the observed activity; the radio waves are produced along open field lines within $\sim 10^7$ cm of the neutron star surface, and the beam periodically sweeps across the observer as the star spins. According to this picture, FRB 20191221A is a factor $\sim 10^{12}$ scaled up version of galactic pulsars with one major difference: whereas pulsars convert rotational kinetic energy to EM waves, the outbursts of 20191221A require conversion of magnetic energy to radiation.

The epoch in which galactic discs settle is a major benchmark to test models of galaxy formation and evolution but is as yet largely unknown. Once discs settle and become self-gravitating enough, stellar bars are able to form; therefore, determining the ages of bars can shed light on the epoch of disc settling, and on the onset of secular evolution. Nevertheless, until now, timing when the bar formed has proven challenging. In this work, we present a new methodology for obtaining the bar age, using the star formation history of nuclear discs. Nuclear discs are rotation-supported structures, built by gas pushed to the centre via bar-induced torques, and their formation is thus coincident with bar formation. In particular, we use integral field spectroscopic (IFS) data from the TIMER survey to disentangle the star formation history of the nuclear disc from that of the underlying main disc, which enables us to more accurately determine when the nuclear disc forms. We demonstrate the methodology on the galaxy NGC 1433 -- which we find to host an old bar that is $8.0^{+1.6}_{-1.1}\rm{(sys)}^{+0.2}_{-0.5}\rm{(stat)}$ Gyr old -- and describe a number of tests carried out on both the observational data and numerical simulations. In addition, we present evidence that the nuclear disc of NGC 1433 grows in accordance with an inside-out formation scenario. This methodology is applicable to high-resolution IFS data of barred galaxies with nuclear discs, making it ideally suited for the TIMER survey sample. In the future we will thus be able to determine the bar age for a large sample of galaxies, shedding light on the epoch of disc settling and bar formation.

We present a $L$-band (2.98--3.96$\mu$m) spectroscopic study of 8 young L dwarfs with spectral types ranging from L2 to L7. Our spectra (${\lambda}/{\Delta \lambda}\approx$ 250 to 600) were collected using the Gemini Near-InfraRed Spectrograph. We first examine the young $L$-band spectral sequence, most notably analyzing the evolution of the $Q$-branch of methane absorption feature at 3.3 $\mu$m. We find the $Q$-branch feature first appears between L3 and L6, as previously seen in older field dwarfs. Secondly, we analyze how well various atmospheric models reproduce the $L$-band and published near-IR (0.7--2.5 $\mu$m) spectra of our objects by fitting five different grids of model spectra to the data. Best-fit parameters for the combined near-IR and $L$-band data are compared to best-fit parameters for just the near-IR data, isolating the impact that the addition of the $L$-band has on the results. This addition notably causes a $\sim$100 K drop in the best-fit effective temperature. Also, when clouds and a vertical mixing rate ($K_{\mathrm{zz}}$) are included in the models, thick clouds and higher $K_{\mathrm{zz}}$ values are preferred. Five of our objects also have previously published effective temperatures and surface gravities derived using evolutionary models, age estimates, and bolometric luminosities. Comparing model spectra matching these parameters to our spectra, we find disequilibrium chemistry and clouds are needed to match these published effective temperatures. Three of these objects are members of AB Dor, allowing us to show the temperature dependence of the $Q$-branch of methane.

Hyper-luminous Quasi-Stellar Objects (QSOs) represent the ideal laboratory to investigate Active Galactic Nuclei (AGN) feedback mechanism since their formidable energy release causes powerful winds at all scales and thus the maximum feedback is expected. We aim at deriving the mean Spectral Energy Distribution (SED) of a sample of 85 WISE-SDSS Selected Hyper-luminous (WISSH) quasars. Since the SED provides a direct way to investigate the AGN structure, our goal is to understand if quasars at the bright end of the luminosity function have peculiar properties compared to the bulk of the population. We built a mean intrinsic SED after correcting for the dust extinction, absorption and emission lines and intergalactic medium absorption. We also derived bolometric, IR band and monochromatic luminosities together with bolometric corrections at lambda = 5100 A and 3 micron. We define a new relation for the 3 micron bolometric correction. We find that the mean SED of hyper-luminous WISSH QSOs is different from that of less luminous sources, i.e. a relatively lower X-ray emission and a near and mid IR excess which can be explained assuming a larger dust contribution. WISSH QSOs have stronger emission from both warm and very hot dust, the latter being responsible for shifting the typical dip of the AGN SED from 1.3 to 1.1 micron. We also derived the mean SEDs of two sub-samples created according to the presence of Broad Absorption Lines and equivalent width of CIV line. We confirm that BALs are X-ray weak and that they have a reddened UV-optical continuum. We also find that BALs tend to have stronger emission from the hot dust component. This analysis suggests that hyper-luminous QSOs have a peculiar SED compared to less luminous objects. It is therefore critical to use SED templates constructed exclusively from very bright quasars samples when dealing with particularly luminous sources.

In this work, we use ~500 low-redshift (z ~ 0.1) X-ray AGNs observed by XMM-Newton and SDSS to investigate the prevalence and nature of AGNs that apparently lack optical emission lines (``optically dull AGNs''). Although 1/4 of spectra appear absorption-line dominated in visual assessment, line extraction with robust continuum subtraction from the MPA/JHU catalog reveals usable [OIII] measurements in 98% of the sample, allowing us to study [OIII]-underluminous AGNs together with more typical AGNs in the context of the L$_{\mathrm{[OIII]}}$--L$_{X}$ relation. We find that ``optically dull AGNs'' do not constitute a distinct population of AGNs. Instead, they are the [OIII]-underluminous tail of a single, unimodal L$_{\mathrm{[OIII]}}$--L$_{X}$ relation that has substantial scatter (0.6 dex). We find the degree to which an AGN is underluminous in [OIII] correlates with the specific SFR or D$_{4000}$ index of the host, which are both linked to the molecular gas fraction. Thus the emerging physical picture for the large scatter seems to involve the gas content of the narrow-line region. We find no significant role for previously proposed scenarios for the presence of optically dull AGNs, such as host dilution or dust obscuration. Despite occasionally weak lines in SDSS spectra, >80% of X-ray AGNs are identified as such with the BPT diagram. >90% are classified as AGNs based only on [NII]/H$\alpha$, providing more complete AGN samples when [OIII] or H$\beta$ are weak. X-ray AGNs with LINER spectra obey essentially the same \lxo\ relation as Seyfert 2s, suggesting their line emission is produced by AGN activity.

We present a Chandra study of surface brightness fluctuations in the diffuse intracluster medium of the Perseus Cluster. Our study utilizes deep, archival imaging of the cluster core as well as a new mosaic of 29 short 5 ks observations extending in 8 different directions out to radii of r_500 ~ 2.2r_2500. Under the assumption that the distribution of densities at a given radius is log-normally distributed, two important quantities can be derived from the width of the log-normal density distribution on a given spatial scale: the density bias, which is equal to the square root of the clumping factor C; and the one-component turbulent velocity, v_(k, 1D). We forward-model all contributions to the measured surface brightness, including astrophysical and particle background components, and account for the Poisson nature of the measured signal. Measuring the distribution of surface brightness fluctuations in 1 arcmin^2 regions, spanning the radial range 0.3-2.2 r_2500 (7.8-57.3 arcmin), we find a small to moderate average density bias of around 3% at radii below 1.6r_2500. We also infer an average turbulent velocity at these radii of v_1D <400 km s^-1. Direct confirmation of our results on turbulent velocities inferred from surface brightness fluctuations should be possible using the X-ray calorimeter spectrometers to be flown aboard the XRISM and Athena. observatories.

The characteristic variability of cool giants and supergiants is attributed to a combination of stellar pulsation and large-scale convective flows. Full 3D radiation-hydrodynamical modelling is an essential tool for understanding the nature of these dynamical processes. The parameter space in our 3D model grid of red giants has expanded in recent years. These models can provide many insights on the nature and properties of the pulsations, including the interplay between convection and pulsations. We treat 3D dynamical models of asymptotic giant branch (AGB) stars and red supergiants similar to observational data. We aim to explore the relation between stellar parameters and the properties of the self-excited pulsations. Output from global 'star-in-a-box' models computed with the CO5BOLD radiation-hydrodynamics code were analysed, particularly in regards to the pulsation properties, to find possible correlations with input and emergent stellar parameters. The fast Fourier transform was applied to spherically averaged mass flows to identify possible radial pulsation periods beneath the photosphere of the modelled stars. Stellar parameters were investigated for correlations with the extracted pulsation periods. We find that the pulsation periods varied with the stellar parameters in good agreement with the current expectations. The pulsation periods follow Ritter's period-mean density relation well and our AGB models agree with period-luminosity relations derived from observations. A mass estimate formula was derived from the 3D models, relating the stellar mass to the fundamental mode pulsation period and the stellar radius. While the non-linearity of the interplay between the self-excited pulsations and the self-consistent convection complicates analyses, the resulting correlations are in good agreement with respect to current theoretical and observational understandings.

Ultraluminous X-ray Sources (ULXs) in globular clusters are low mass X-ray binaries that achieve high X-ray luminosities through a currently uncertain accretion mechanism. Using archival Chandra and Hubble Space Telescope observations, we perform a volume-limited search ($\lesssim$ 70 Mpc) of 21 of the most massive ($>10^{11.5} M_\odot$) early-type galaxies to identify ULXs hosted by globular cluster (GC) candidates. We find a total of 34 ULX candidates above the expected background within 5 times the effective radius of each galaxy, with 10 of these ($\sim29.4\%$) potentially hosted by a GC. A comparison of the spatial and luminosity distributions of these new candidate GC ULXs with previously identified GC ULXs shows that they are similar: both samples peak at $L_X \sim$ a few $\times 10^{39}$ erg/s and are typically located within a few effective radii of their host galaxies.

Young exoplanets are attractive targets for atmospheric characterization to explore the early phase of planetary evolution and the surrounding environment. Recent observations of the 10 Myr young Neptune-sized exoplanet K2-33b revealed that the planet's transit depth drastically decreases from the optical to near-infrared wavelengths. Thao et al. (2022) suggested that a thick planetary haze and/or stellar spots may be the cause; however, even the best-fit model only barely explains the data. Here, we propose that the peculiar transmission spectrum may indicate that K2-33b possess a circumplanetary dust ring, an analog of Jupiter's dust ring. We demonstrate that the ring could produce a steep slope in the transmission spectrum even if its optical depth is as low as $\sim {10}^{-2}$. We then apply a novel joint atmosphere-ring retrieval to K2-33b and find that the ring scenario could well explain the observed spectrum for various possible ring compositions. Importantly, the dust ring also exhibits prominent absorption features of ring particles around $\sim$10 $\rm{\mu}m$, whose shape and strength depend on the composition of the ring. Thus, future observations by JWST-MIRI would be able to test not only the ring hypothesis but also, if it indeed exists, to constrain the composition of the ring -- providing a unique opportunity to explore the origins of the dust ring around its parent planet, soon after the planetary system's formation.

Although all-sky surveys have led to the discovery of dozens of young planets, little is known about their atmospheres. Here, we present multi-wavelength transit data for the super Neptune-sized exoplanet, K2-33b -- the youngest (~10 Myr) transiting exoplanet to-date. We combined photometric observations of K2-33 covering a total of 33 transits spanning >2 years, taken from K2, MEarth, Hubble, and Spitzer. The transit photometry spanned from the optical to the near-infrared (0.6-4.5$\mu$m), enabling us to construct a transmission spectrum of the planet. We find that the optical transit depths are nearly a factor of two deeper than those from the near-infrared. This difference holds across multiple datasets taken over years, ruling out issues of data analysis and unconstrained systematics. Surface inhomogeneities on the young star can reproduce some of the difference, but required spot coverage fractions (>60%) are ruled out by the observed stellar spectrum(<20%). We find a better fit to the transmission spectrum using photochemical hazes, which were predicted to be strong in young, moderate-temperature, and large-radius planets like K2-33b. A tholin haze with CO as the dominant gaseous carbon carrier in the atmosphere can reasonably reproduce the data with small or no stellar surface inhomogeneities, consistent with the stellar spectrum. The HST data quality is insufficient for the detection of any molecular features. More observations would be required to fully characterize the hazes and spot properties and confirm the presence of CO suggested by current data.

The element abundances in the solar corona and solar wind are often different from those of the solar photosphere, typically with a relative enrichment of elements with low first ionization potential (FIP effect). Here we study the spatial distribution and temporal evolution of the coronal chemical composition in an active region (AR) over about 10 days, using Hinode/EIS spectra, and we also analyze coordinated IRIS observations of the chromospheric and transition region emission to investigate any evidence of the footprints of the FIP effect in the lower atmosphere. To derive the coronal abundances we use a spectral inversion method recently developed for the MUSE investigation (Cheung et al. 2019, De Pontieu et al. 2020). We find that in the studied active region (AR 12738) the coronal FIP bias presents significant spatial variations, with its highest values (~ 2.5-3.5) in the outflow regions at the boundary of the AR, but typically modest temporal variability. Some moss regions and some regions around the AR sunspot show enhanced FIP (~ 2-2.5) with respect to the AR core, which has only a small FIP bias of ~1.5. The FIP bias appears most variable in some of these moss regions. The IRIS observations reveal that the chromospheric turbulence, as derived from IRIS2 inversions of the Mg II spectra, is enhanced in the outflow regions characterized by the high FIP bias, providing significant new constraints to both models aimed at explaining the formation of AR outflows and models of chemical fractionation.

We present a new diagnostic diagram for local ultraluminous infrared galaxies (ULIRGs) and quasars, analysing particularly the Spitzer Space Telescope's Infrared Spectrograph (IRS) spectra of 102 local ULIRGs and 37 Palomar Green quasars. Our diagram is based on a special non-linear mapping of these data, employing the Kernel Principal Component Analysis method. The novelty of this map lies in the fact that it distributes the galaxies under study on the surface of a well-defined ellipsoid, which, in turn, links basic concepts from geometry to physical properties of the galaxies. Particularly, we have found that the equatorial direction of the ellipsoid corresponds to the evolution of the power source of ULIRGs, starting from the pre-merger phase, moving through the starburst-dominated coalescing stage towards the active galactic nucleus (AGN)-dominated phase, and finally terminating with the post-merger quasar phase. On the other hand, the meridian directions distinguish deeply obscured power sources of the galaxies from unobscured ones. These observations have also been verified by comparison with simulated ULIRGs and quasars using radiative transfer models. The diagram correctly identifies unique galaxies with extreme features that lie distinctly away from the main distribution of the galaxies. Furthermore, special two-dimensional projections of the ellipsoid recover almost monotonic variations of the two main physical properties of the galaxies, the silicate and PAH features. This suggests that our diagram naturally extends the well-known Spoon diagram and it can serve as a diagnostic tool for existing and future infrared spectroscopic data, such as those provided by the James Webb Space Telescope.

While giant planet occurrence rates increase with stellar mass, occurrence rates of close-in super-Earths decrease. This is in contradiction to the expectation that the total mass of the planets in a system scale with the protoplanetary disc mass and hence the stellar mass. Since the snow line plays an important role in the planet formation process we examine differences in the temperature structure of protoplanetary gas discs around stars of different mass. Protoplanetary discs likely contain a dead zone at the midplane that is sufficiently cold and dense for the magneto-rotational instability to be suppressed. As material builds up, the outer parts of the dead zone may be heated by self-gravity. The temperature in the disc can be below the snow line temperature far from the star and in the inner parts of a dead zone. The inner icy region has a larger radial extent around smaller mass stars. The increased mass of solid icy material may allow for the in situ formation of larger and more numerous planets close to a low-mass star. Super-Earths that form in the inner icy region may have a composition that includes a significant fraction of volatiles.

This paper identifies spectral lines from X-ray to infrared wavelengths which are optimally suited to measuring vector magnetic fields as high as possible in the solar atmosphere. Instrumental and Earth's atmospheric properties, as well as solar abundances, atmospheric properties and elementary atomic physics are considered without bias towards particular wavelengths or diagnostic techniques. While narrowly-focused investigations of individual lines have been reported in detail, no assessment of the comparative merits of all lines has ever been published. Although in the UV, on balance the Mg+ h and k lines near 2800 Angstroms are optimally suited to polarimetry of plasma near the base of the solar corona. This result was unanticipated, given that longer-wavelength lines offer greater sensitivity to the Zeeman effect. While these lines sample optical depths photosphere to the coronal base, we argue that cores of multiple spectral lines provide a far more discriminating probe of magnetic structure as a function of optical depth than the core and inner wings of a strong line. Thus, together with many chromospheric lines of Fe+ between 2585 and the h line at 2803 Angstrom, this UV region promises new discoveries concerning how the magnetic fields emerge, heat, and accelerate plasma as they battle to dominate the force and energy balance within the poorly-understood chromosphere.

We are of the opinion that several anomalies in the atmosphere of Venus provide evidence of yet-unknown processes and systems that are out of equilibrium. The investigation of these anomalies on Venus should be open to the wide range of explanations, including unknown biological activity. We provide an overview of two anomalies, the tentative detection of ammonia and phosphine in Venus's atmosphere. These anomalies fly in the face of the tacit assumption that the atmosphere of Venus must be in chemical redox equilibrium, an assumption connected to the belief that Venus is lifeless. We then discuss several major past discoveries in astronomy, biology and geology, which lead to the abandonment of certain assumptions held by many scientists as though they were well-established principles. The anomalies of ammonia and phosphine in the atmosphere of Venus are placed in the context of these historical discoveries. This context supports our opinion that persistence by the community in the exploration of these anomalies with a skeptical eye towards tacit assumptions will increase the chances of making profound discoveries about the atmosphere of Venus and the diverse and often strange nature of planetary environments.

The Small Magellanic Cloud (SMC) is the closest low-metallicity galaxy to the Milky Way where the dynamical state of molecular clouds can be analyzed. We present a CO(2-1) survey at 9 pc resolution obtained with the APEX telescope in an extensive region of the SMC and characterize the properties of the molecular clouds. We study the dynamical state and stability of these clouds uniformly. We identify 177 molecular clouds within the SMC by using CPROPS, of which 124 clouds are fully resolved with signal-to-noise ratio $>$ 5. The scaling relationships show that the SMC clouds are (on average) less turbulent and less luminous than their inner Milky Way counterparts of similar size by a factor of 2 and 3, respectively, while for a fixed linewidth, the SMC clouds are over-luminous by a factor of 3.5. Using the virial masses, we derive a CO-to-H2 conversion factor for the SMC CO clouds of 10.5 M$_{\odot}(K km s^{-1} pc^{2})^{-1}$, measured at 9 pc resolution. We also determine a dust-based conversion factor of 28 M$_{\odot}(K km s^{-1} pc^{2})^{-1}$, obtained at 12 pc resolution. We find that the SMC clouds appear to be in approximate gravitational virial equilibrium. We find that the cumulative mass functions based on both the luminous mass and the virial mass are steeper than $\frac{dN}{dM} \propto M^{-2}$, suggesting that most of the molecular mass of the SMC is contained in low-mass clouds.

Trajectory optimization of low-thrust perturbed orbit rendezvous is a crucial technology for space missions in low Earth orbits, which is difficult to solve due to its initial value sensitivity, especially when the transfer trajectory has many revolutions. This paper investigated the time-fixed perturbed orbit rendezvous between low-eccentricity orbits and proposed a priori quasi-optimal thrust strategy to simplify the problem into a parametric optimization problem, which significantly reduces the complexity. The optimal trajectory is divided into three stages including transfer to a certain intermediate orbit, thrust-off drifting and transfer from intermediate orbit to the target orbit. In the two transfer stages, the spacecraft is assumed to use a parametric law of thrust. Then, the optimization model can be then obtained using very few unknowns. Finally, a differential evolution algorithm is adopted to solve the simplified optimization model and an analytical correction process is proposed to eliminate the numerical errors. Simulation results and comparisons with previous methods proved this new method's efficiency and high precision for low-eccentricity orbits. The method can be well applied to premilitary analysis and high-precision trajectory optimization of missions such as in-orbit service and active debris removal in low Earth orbits.

ALMA observations of the long wavelength dust continuum are used to estimate the gas masses in a sample of 708 star-forming (SF) galaxies at z = 0.3 to 4.5. We determine the dependence of gas masses and star formation efficiencies (SFE=SFR per unit gass mass). We find that 70 percent of the increase in SFRs of the MS is due to the increased gas masses at earlier epochs while 30 percent is due to increased efficiency of SF. For galaxies above the MS this is reversed with 70 percent of the increased SFR relative to the MS being due to elevated SFEs. Thus, the major evolution of star formation activity at early epochs is driven by increased gas masses, while the starburst activity taking galaxies above the MS is due to enhanced triggering of star formation (likely due to galactic merging). The interstellar gas peaks at z = 2 and dominates the stellar mass down to z = 1.2. Accretion rates needed to maintain continuity of the MS evolution exceed 100 Msun per yr at z > 2. The galactic gas contents are likely the driving determinant for both the rise in SF and AGN activity from z = 5 to their peak at z = 2 and subsequent fall to lower z. We suggest that for self-gravitating clouds with supersonic turbulence, cloud collisions and the filamentary structure of the clouds regulate the star formation activity.

This is a report on some highlights of some research on the rare light elements, lithium (Li), beryllium (Be), and boron (B), that I presented in my Henry Norris Russell Lecture in January, 2020. It is not a comprehensive review of work on these light elements, but contains sections on Big Bang nucleosynthesis of Li and the rarity of these light elements. It includes information on how they are observed, both historically and currently, and the difficulties entailed in determining their abundances. %It traces the evolution of these elements through the history of the Galaxy. The production of Li, Be, and B is ongoing so the youngest stars contain the most Li, and they have had less time to destroy it. All three elements are readily destroyed in stellar interiors, but have differing degrees of susceptibility to the particular nuclear fusion reactions which deplete their surface content. This feature makes them remarkably good probes into the otherwise unobservable interiors of stars. It also enhances the use of two or more of the three in sorting out the various processes at work in the insides of stars.

There are a handful of spectral windows in the near-IR through which we can see down to Venus' surface on the night side of the planet. The surface of our sister planet has thus been imaged by sensors on Venus-orbiting platforms (Venus Express, Akatsuki) and during fly-by with missions to other planets (Galileo, Cassini). The most tantalizing finding, so far, is the hint of possible active volcanism. However, the thermal radiation emitted by the searing (c. 475 degrees C) surface of Venus has to get through the opaque clouds between 50 and 70 km altitude, as well as the sub-cloud atmosphere. In the clouds, the light is not absorbed but scattered, many times. This results in blurring the surface imagery to the point where the smallest discernible feature is roughly 100 km in size, full-width half-max (FWHM), and this has been explained using numerical models. We describe a new analytical modeling framework for predicting the width of the atmospheric point-spread function (APSF), which is what determines the effective resolution of surface imaging from space. Our best estimates of the APSF width for the 1-to-1.2 micron spectral range are clustered around 130 km FWHM. Interestingly, this is somewhat larger than the accepted value of about 100 km, which is based on visual image inspection and numerical simulations.

We first discuss and determine the isospin mixing of the two $2^-$ states (12.53 MeV and 12.97 MeV) of $^{16}$O nucleus using the inelastic electron scattering data. We then evaluate the cross section of 4.4-MeV $\gamma$ rays produced in the neutrino neutral-current (NC) reaction $^{16}$O($\nu, \nu^{\prime}$)$^{16}$O$(12.97~{\rm MeV}, 2^-$) in a water Cherenkov detector at the low energy below 100 MeV. The detection of $\gamma$ rays for $E_{\gamma}>5$~MeV from the NC reaction $^{16}$O($\nu, \nu^{\prime}$)$^{16}$O$(E_x>16\ {\rm MeV}, T=1$) with a water Cherenkov detector in the supernova neutrino bursts has been proposed and discussed by several authors previously. In this article, we discuss a new NC reaction channel from $^{16}$O(12.97 ${\rm MeV}, 2^-$) producing a 4.4-MeV $\gamma$ ray, the cross section of which is more robust and even larger at the low energy ($E_{\nu}<25$ MeV) than the NC cross section from $^{16}$O$(E_x>16\ {\rm MeV}, T=1$). We also evaluate the number of such events induced by neutrinos from supernova explosion which can be observed by the Super-Kamiokande, a 32 kton water Cherenkov detector in the Earth.

We present the survey design, implementation, and outlook for COSMOS-Web, a 255 hour treasury program conducted by the James Webb Space Telescope in its first cycle of observations. COSMOS-Web is a contiguous 0.54 deg$^2$ NIRCam imaging survey in four filters (F115W, F150W, F277W, and F444W) that will reach 5$\sigma$ point source depths ranging $\sim$27.5-28.2 magnitudes. In parallel, we will obtain 0.19 deg$^2$ of MIRI imaging in one filter (F770W). COSMOS-Web will build on the rich heritage of multiwavelength observations and data products available in the COSMOS field. The design of COSMOS-Web is motivated by three primary science goals: (1) to discover thousands of galaxies in the Epoch of Reionization ($6<z<11$) and map reionization's spatial distribution, environments, and drivers on scales sufficiently large to mitigate cosmic variance, (2) to identify hundreds of rare quiescent galaxies at $z>4$ and place constraints on the formation of the Universe's most massive galaxies ($M_\star>10^{10}$\,M$_\odot$), and (3) directly measure the evolution of the stellar mass to halo mass relation using weak gravitational lensing out to $z\sim2.5$ and measure its variance with galaxies' star formation histories and morphologies. In addition, we anticipate COSMOS-Web's legacy value to reach far beyond these scientific goals, touching many other areas of astrophysics, such as the identification of the first direct collapse black hole candidates, ultracool sub-dwarf stars in the Galactic halo, and possibly the identification of $z>10$ pair-instability supernovae. In this paper we provide an overview of the survey's key measurements, specifications, goals, and prospects for new discovery.

It has been suggested in the literature that the sound velocity of the nuclear matter $v_s$ violates the so-called sound velocity bound $v_s \le c/\sqrt{3}$ at high density, where $c$ is the speed of light. In this paper, we revisit this issue and confront the current measurements of mass, radius, and tidal deformability of neutron stars with $10^5$ different equations of state which are parametrized at low density and saturates the sound velocity bound beyond twice the saturation density where the equation of state has not been constrained yet, by which we can conservatively obtain the maximum mass of the neutron stars compatible both with the observed properties of neutron stars and the sound velocity bound. We find that majority of the models are eliminated by the incompatibility with the observations and, especially, the recently detected massive pulsar ($2.35\pm 0.17 M_\odot$) is hardly realized by our simulations. Our study strongly supports the violation of the sound velocity bound.

Based upon a metallicity selection, we identify a significant sub-population of the inner halo globular clusters in the Andromeda Galaxy which we name the Dulais Structure. It is distinguished as a co-rotating group of 10-20 globular clusters which appear to be kinematically distinct from, and on average more metal-poor than, the majority of the inner halo population. Intriguingly, the orbital axis of this Dulais Structure is closely aligned with that of the younger accretion event recently identified using a sub-population of globular clusters in the outer halo of Andromeda, and this is strongly suggestive of a causal relationship between the two. If this connection is confirmed, a natural explanation for the kinematics of the globular clusters in the Dulais Structure is that they trace the accretion of a substantial progenitor (~10^11 Msun) into the halo of Andromeda during the last few billion years, that may have occurred as part of a larger group infall.

We perform a novel multi-messenger analysis for the identification and parameter estimation of the Standing Accretion Shock Instability (SASI) in a core collapse supernova with neutrino and gravitational wave (GW) signals. In the neutrino channel, this method performs a likelihood ratio test for the presence of SASI in the frequency domain. For gravitational wave signals we process an event with a modified constrained likelihood method. Using simulated supernova signals, the properties of the Hyper-Kamiokande neutrino detector, and O3 LIGO Interferometric data, we produce the two-dimensional probability density distribution (PDF) of the SASI activity indicator and calculate the probability of detection $P_\mathrm{D}$ as well as the false identification probability $P_\mathrm{FI}$. We discuss the probability to establish the presence of the SASI as a function of the source distance in each observational channel, as well as jointly. Compared to a single-messenger approach, the joint analysis results in $P_\mathrm{D}$ (at $P_\mathrm{FI}=0.1$) of SASI activities that is larger by up to $\approx~40\%$ for a distance to the supernova of 5 kpc. We also discuss how accurately the frequency and duration of the SASI activity can be estimated in each channel separately. Our methodology is suitable for implementation in a realistic data analysis and a multi-messenger setting.

The climate of a planet can be strongly affected by its eccentricity due to variations in the stellar flux. There are two limits for the dependence of the inner habitable zone boundary (IHZ) on eccentricity: (1) the mean-stellar flux approximation ($S_{\mbox{IHZ}} \propto \sqrt{1-e^2}$), in which the temperature is approximately constant throughout the orbit, and (2) the maximum-stellar flux approximation ($S_{\mbox{IHZ}} \propto (1-e)^2$), in which the temperature adjusts instantaneously to the stellar flux. Which limit is appropriate is determined by the dimensionless parameter $\Pi = \frac{C}{BP}$, where $C$ is the heat capacity of the planet, $P$ is the orbital period, and $B=\frac{\partial \Omega}{\partial T_s}$, where $\Omega$ is the outgoing longwave radiation and $T_s$ is the surface temperature. We use the Buckingham $\Pi$ theorem to derive an analytical function for the IHZ in terms of eccentricity and $\Pi$. We then build a time-dependent energy balance model to resolve the surface temperature evolution and constrain our analytical result. We find that $\Pi$ must be greater than about $\sim 1$ for the mean-stellar flux approximation to be nearly exact and less than about $\sim 0.01$ for the maximum-stellar flux approximation to be nearly exact. In addition to assuming a constant heat capacity, we also consider the effective heat capacity including latent heat (evaporation and precipitation). We find that for planets with an Earth-like ocean, the IHZ should follow the mean-stellar flux limit for all eccentricities. This work will aid in the prioritization of potentially habitable exoplanets with non-zero eccentricity for follow-up characterization.

We report the discovery of 2M1222-57 as a low-mass, pre-main-sequence (PMS) eclipsing binary (EB) in the Lower Centaurus Crux (LCC) association for which, using Gaia parallaxes and proper motions with a neural-net age estimator, we determine an age of 16.2$\pm$2.2 Myr. The broadband spectral energy distribution (SED) shows clear excess at ~10 um indicative of a circumbinary disk, and new speckle-imaging observations reveal a faint, tertiary companion separated by ~100 AU. H-alpha emission is modulated on the orbital period, consistent with theoretical models of orbitally pulsed accretion streams reaching from the inner disk edge to the central stars. From a joint analysis of spectroscopically determined radial velocities and TESS light curves, together with additional tight constraints provided by the SED and the Gaia parallax, we measure masses for the eclipsing stars of 0.74 Msun and 0.67 Msun; radii of 0.98 Rsun and 0.94 Rsun; and effective temperatures of 3750 K and 3645 K. The masses and radii of both stars are measured to an accuracy of ~1%. The measured radii are inflated, and the temperatures suppressed, relative to predictions of standard PMS evolutionary models at the age of LCC; also, the Li abundances are ~2 dex less depleted than predicted by those models. However, models that account for the global and internal effects of surface magnetic fields are able to simultaneously reproduce the measured radii, temperatures, and Li abundances at an age of 17.0$\pm$0.5 Myr. Altogether, the 2M1222-57 system presents very strong evidence that magnetic activity in young stars alters both their global properties and the physics of their interiors.

We present here a validation of sub-Saturn exoplanet TOI-181b orbiting a K spectral type star TOI-181 (Mass: 0.781 $\pm$ 0.042 M$_{\odot}$, Radius: 0.74 $\pm$ 0.13 R$_{\odot}$, Temperature: 4994 $\pm$ 50 K) as a part of Validation of Transiting Exoplanets using Statistical Tools (VaTEST) project. TOI-181b is a planet with radius 0.631521 $\pm$ 0.0086 R$_J$, mass 0.127891 $\pm$ 0.01929 M$_J$, orbiting in a slightly eccentric orbit with eccentricity 0.193 $\pm$ 0.066 and semi-major axis of 0.05578 $\pm$ 0.00473 AU, with an orbital period of 4.532058 $\pm$ 0.000002 days. The transit photometry data was collected using Transiting Exoplanet Survey Satellite (TESS) and spectroscopic data for radial velocity analysis was collected using The European Southern Observatory's (ESO) High Accuracy Radial Velocity Planet Searcher (HARPS) telescope. Based on the radial velocity best-fit model we measured RV semi-amplitude to be 18.8615 $\pm$ 2.1312 ms$^{-1}$. We carried out LATTE tests to make sure the transit signal was caused by a planet and not by any instrument fault or false positive scenarios. Additionally, we used VESPA and TRICERATOPS to compute the False Positive Probability (FPP), and the findings were FPP values of $1.68\times10^{-14}$ and $3.812\times10^{-04}$, respectively, which are significantly lower than the 1\% threshold. The finding of TOI-181b is significant in the perspective of future work on the formation and migration history of analogous planetary systems since warm sub-Saturns are uncommon in the known sample of exoplanets.

The nonlinear Lagrangian displacement field and initial linear density field are highly correlated. Therefore, reconstructing the nonlinear displacement field could better extract the primordial cosmological information from the late time density field. Continuing from Ref. [1], we investigate to what extent the iterative displacement reconstruction in Ref. [2] can recover the true displacement field with a particular emphasis on improving the numerical discreteness effect and improving the perturbation theory model for the postreconstructed field. We propose two new displacement field estimators to overcome the numerical discreteness effect in measuring the displacement field; the normalized momentum estimator (NME) and the rescaled resumed estimator (RRE). Furthermore, by adequately accounting for the UV-sensitive term in a reconstruction procedure with an effective field theory approach, we manage to decrease the discrepancy between the theoretical model and the simulation result almost five times, to the level of a few \% at $k\sim 0.2h/{\rm Mpc}$ for the redshift $z=0.6$. We conclude that the previous displacement field reconstruction cannot fully reconstruct the shape of the nonlinear displacement field at $k\sim 0.2h/{\rm Mpc}$ at $z=0.6$, while it is still an efficient BAO reconstruction method.

The rotational states of the members in the dwarf planet - satellite systems in the transneptunian region are determined by the formation conditions and the tidal interaction between the components, and these rotational characteristics are the prime tracers of their evolution. Previously a number of authors claimed highly diverse values for the rotation period for the dwarf planet Eris, ranging from a few hours to a rotation (nearly) synchronous with the orbital period (15.8 d) of its satellite, Dysnomia. In this letter we present new light curve data of Eris, taken with $\sim$1-2m-class ground based telescopes, and with the TESS and Gaia space telescopes. TESS data could not provide a well-defined light curve period, but could constrain light curve variations to a maximum possible light curve amplitude of $\Delta m$ $\leq$ 0.03 mag (1-$\sigma$) for P $\leq$ 24 h periods. Both the combined ground-based data and the Gaia measurements unambiguously point to a light curve period equal to the orbital period of Dysnomia, P = 15.8 d, with a light curve amplitude of $\Delta m$ $\approx$ 0.03 mag, i.e. the rotation of Eris is tidally locked. Assuming that Dysnomia has a collisional origin, calculations with a simple tidal evolution model show that Dysnomia has to be relatively massive (mass ratio of q = 0.01--0.03) and large (radius of $R_s$ $\geq$ 300 km) to slow down Eris to synchronized rotation. These simulations also indicate that -- assuming tidal parameters usually considered for transneptunian objects -- the density of Dysnomia should be 1.8-2.4 $g cm^{-3}$, an exceptionally high value among similarly sized transneptunian objects, putting important constraints on the formation conditions.

The HH80-81 system is one of the most powerful jets driven by a massive protostar. We present new near-infrared (NIR) line imaging observations of the HH80-81 jet in the H$_2$ (2.122 $\mu$m) and [Fe II] (1.644 $\mu$m) lines. These lines trace not only the jet close to the exciting source but also the knots located farther away. We have detected nine groups of knot-like structures in the jet including HH80 and HH81 spaced $0.2-0.9$ pc apart. The knots in the northern arm of the jet show only [Fe II] emission closer to the exciting source, a combination of [Fe II] and H$_2$ at intermediate distances, and solely H$_2$ emission farther outwards. Towards the southern arm, all the knots exhibit both H$_2$ and [Fe II] emission. The nature of the shocks is inferred by assimilating the NIR observations with radio and X-ray observations from literature. In the northern arm, we infer the presence of strong dissociative shocks, in the knots located close to the exciting source. The knots in the southern arm that include HH80 and HH81 are explicable as a combination of strong and weak shocks. The mass-loss rates of the knots determined from [Fe II] luminosities are in the range $\sim 3.0\times 10^{-7}-5.2\times 10^{-5}$ M$_{\odot}$ yr$^{-1}$, consistent with those from massive protostars. Towards the central region, close to the driving source of the jet, we have observed various arcs in H$_2$ emission which resemble bow shocks, and strings of H$_2$ knots which reveal traces of multiple outflows.

The number density of extragalactic 21-cm radio sources as a function of their spectral line-widths -- the HI width function (HIWF) -- is in principle a sensitive tracer of the dark matter halo mass function (HMF). The $\Lambda$ cold dark matter model predicts that the HMF should be identical everywhere provided it is sampled in sufficiently large volumes, implying that the same should be true of the HIWF. The ALFALFA 21-cm survey measured the HIWF in two separate, northern ('spring') and southern ('fall') Galactic fields and found a systematically higher number density of sources in the spring field. Taken at face value, this is in tension with theoretical predictions. We use the Sibelius-DARK N-body simulation and the semi-analytical galaxy formation model GALFORM to create a mock ALFALFA survey. We find that the offset in number density likely has two origins: the sensitivity of the survey is different in the two survey fields, which has not been correctly accounted for in previous measurements; and the $1/V_{\mathrm{eff}}$ algorithm used for completeness corrections does not fully account for biases arising from spatial clustering in the galaxy distribution. The latter bias is primarily driven by a foreground overdensity in the spring field within a distance of $30\,\mathrm{Mpc}$, but more distant structure also plays a role. We provide an updated measurement of the ALFALFA HIWF (and HIMF) correcting for the variations in survey sensitivity. Only when systematic effects such as these are understood and corrected for can the HIWF fulfil its potential as a test of cosmological models.

We show that the Friedmann-Lema\^{i}tre-Robertson-Walker (FLRW) framework has an instability towards the growth of fluid flow anisotropies, even if the Universe is accelerating. This flow (tilt) instability in the matter sector is invisible to Cosmic No-Hair Theorem-like arguments, which typically only flag shear anisotropies in the metric. We illustrate our claims in the setting of ``dipole cosmology'', the maximally Copernican generalization of FLRW that can accommodate a flow. Simple models are sufficient to show that the cosmic flow need not track the shear, even in the presence of a positive cosmological constant. We also emphasize that the growth of the tilt hair is fairly generic if the total equation of state $w(t) \rightarrow -1$ at late times (as it does in standard cosmology), irrespective of the precise model of dark energy.

The rotation of the fluid cell generates additional pressure and exhibits properties similar to temperature in isotrotropic expansion, affecting the convection criteria in the form of size, with small fluid cells in a natural convection state and large fluid cells in a forced convection state. This is verified in observational data on solar granulation, which also infer that the critical size of the granule is negatively correlated with the value of the local average vortex.

It has been theorized that Dynamical Dark Energy (DDE) could be a possible solution to the Hubble tension. To avoid the degeneracy between the Hubble parameter $H_0$ and the sound horizon scale $r_d$, in this article we use their multiplication as one parameter $c/\left(H_0 r_d\right)$ and we use it to infer cosmological parameters for 6 different models - $\Lambda$CDM and 5 DDE parametrizations -- the Chevallier-Polarski-Linder (CPL), the Barboza-Alcaniz (BA), the Low correlation (LC), the Jassal-Bagla-Padmanabhan (JBP) and the Feng-Shen-Li-Li model. We choose a dataset that treats this combination as one parameter, that includes the Baryon Acoustic Oscillation (BAO) data $0.11 \le z \le 2.40$ and additional points from the Cosmic Microwave Background (CMB) Peaks ($z \approx 1090$). To them, we add the marginalized Pantehon dataset and GRB dataset. We see that the tension is moved from $H_0$ and $r_d$ to $c/\left(H_0 r_d\right)$ and $\Omega_m$. There is only one model that satisfies the Planck 2018 constraints on both parameters and this is LC with a huge error. The rest cannot fit into both constraints. $\Lambda$CDM is preferred with respect to the statistical measures.

Ge et al. reported a peculiar large glitch observed from SGR 1935 recently~\citep{2022arXiv221103246G}. % Interestingly, this glitch occurred about $3.1\pm 2.5~\rm{days}$ before FRB 200428, accompanied by a delayed spin up (DSU) process with a timescale of $8\pm 1~\rm{days}$. % This DSU is the first one detected from magnetars unambiguously, previously, DSU has only been detected from large glitches of the young Crab pulsar. % Strikingly, this DSU even resulted in torque reversal, i.e., the spin down state of SGR 1935 was turned into a spin up one within $8\pm 1~\rm{days}$, the first case observed from isolated neutron stars as far as we know. % What's more, the average positive torque is consistent in value and order of magnitude with those deduced from DSUs observed from the Crab pulsar. % The coincidences of DSU appearance and consistent positive torque in the Crab pulsar and SGR 1935 suggest that, these DSUs have the same physical origin and should be analyzed simultaneously. % Joint analysis on DSUs from the Crab pulsar and SGR 1935 suggests a physical origin independent of (at least not directly dependent on) rotational parameters. % If this new DSU and torque reversal phenomena be further confirmed and accretion be excluded convincingly, nearly all models regarding DSU phenomenon under the framework of vortex motion should be reconsidered, new physics will be required. % We stress that, joint analysis on DSUs from different pulsars with different rotational parameters may provide pivotal clue to reveal their physical origin, even neutron star equation of state and the triggering mechanism(s) of (galactic) fast radio bursts ultimately.

Context. The gamma-ray emitting source WISE J141046.00+740511.2 has been associated with a Fermi-LAT detection by crossmatching with Swift/XRT data. It has shown all the canonical observational characteristics of a BL Lac source, including a power-law, featureless optical spectrum. However, it was only recently detected at radio frequencies and its radio flux is significantly low. Aims. Given that a radio detection is fundamental to associate lower-energy counterparts to Fermi-LAT sources, we aim to unambiguously classify this source by performing a multiwavelength analysis based on contemporaneous data. Methods. By using multifrequency observations at the Jansky Very Large Array, Giant Metrewave Radio Telescope, Gran Telescopio Canarias, Gemini, William Herschel Telescope and Liverpool observatories, together with Fermi-LAT and Swift data, we carried out two kinds of analyses. On one hand, we studied several known parameters that account for the radio loudness or weakness characterization and their application to blazars (in general) and to our source (in particular). And, on the other hand, we built and analyzed the observed spectral energy distribution (SED) of this source to try to explain its peculiar characteristics. Results. The multiwavelength analysis indicates that WISE J141046.00+740511.2 is a blazar of the high-frequency peaked (HBL) type that emits highly polarized light and that is likely located at a low redshift. In addition, the one-zone model parameters that best fit its SED are those of an extreme HBL (EHBL); this blazar type has been extensively predicted in theory to be lacking in the radio emission that is otherwise typical of canonical gamma-ray blazars. Conclusions. We confirm that WISE J141046.00+740511.2 is indeed a highly polarized BL Lac of the HBL type. Further studies will be conducted to explain the atypical low radio flux detected for this source.

We present a first sample of 117 [OIII]$\lambda\lambda$4960,5008-selected star-forming galaxies at $5.33 < z < 6.93$ detected in JWST/NIRCam 3.5$\mu$m slitless spectroscopy of a $6.5 \times 3.4$ arcmin$^2$ field centered on the hyperluminous quasar SDSS J0100+2802, obtained as part of the EIGER (Emission-line galaxies and Intergalactic Gas in the Epoch of Reionization) survey. Three prominent galaxy overdensities are observed, one of them at the redshift of the quasar. Galaxies are found within 200 pkpc and 105 km s$^{-1}$ of four known metal absorption-line systems in this redshift range. We focus on the role of the galaxies in ionizing the surrounding intergalactic medium (IGM) during the later stages of cosmic reionization and construct the mean Ly$\alpha$ and Ly$\beta$ transmission as a function of distance from the galaxies. At the lowest redshifts in our study, $5.3 < z < 5.7$, the IGM transmission rises monotonically with distance from the galaxies. This is as expected when galaxies reside at peaks in the overdensity field of an IGM that is ionized by more or less uniform ionizing background, and has been seen at lower redshifts. In contrast, at $5.7 < z < 6.14$, the transmission of both Ly$\alpha$ and Ly$\beta$ first increases with distance, but then peaks at a distance of 5 cMpc before declining. This peak in transmission is qualitatively similar to that seen (albeit at smaller distances and higher redshifts) in the THESAN simulations. Finally, in the region $6.15 < z < 6.26$ where the additional ionizing radiation from the quasar dominates, the monotonic increase in transmission with distance is re-established. This result is interpreted to represent evidence that the transmission of the IGM at $z \sim 5.9$ towards J0100+2802 results from the ``local'' ionizing radiation of galaxies that dominates over the much reduced cosmic background.

We present emission-line measurements and physical interpretations for a sample of 117 [OIII] emitting galaxies at $z=5.33-6.93$, using the first deep JWST/NIRCam wide field slitless spectroscopic observations. Our 9.7-hour integration is centered upon the $z=6.3$ quasar J0100+2802 -- the first of six fields targeted by the EIGER survey -- and covers $\lambda=3-4$ microns. We detect 133 [OIII] doublets, but merge pairs within $\approx$10 kpc and 600 km s$^{-1}$, motivated by their small scale clustering excess. We detect H$\beta$ in 68 and H$\gamma$ in two galaxies. The galaxies are characterised by a UV luminosity M$_{\rm UV}\sim-19.7$ ($-17.7$ to $-22.3$), mass ~$10^8$ $(10^{6.8-9.8})$ M$_{\odot}$, H$\beta$ and [OIII] EWs $\approx$ 1000 Angstrom (up to 3000 Angstrom), young ages (~5 Myr), a highly excited interstellar medium ([OIII]/H$\beta\approx6$) and low dust attenuations. These high EWs are very rare in the local Universe, but we show they are ubiquitous at $z\sim6$ based on the measured number densities. The stacked spectrum reveals H$\gamma$ and [OIII]$_{4364}$ which shows that the galaxies are typically dust and metal poor (E(B-V)=0.1, 12+log(O/H)=7.4) with a high electron temperature ($2\times10^4$ K) and a high production efficiency of ionising photons ($\xi_{\rm ion}=10^{25.6}$ Hz erg$^{-1}$). We further show the existence of a strong mass-metallicity relation as expected by simulations. The young highly ionising stellar populations, moderately low metallicities, low dust attenuations and high ionisation state in z~6 galaxies conspire to maximise the [OIII] output from galaxies, yielding an [OIII] luminosity density at z~6 that is significantly higher than at z~2, despite the order of magnitude decline in cosmic star formation. Thus, [OIII] emission-line surveys with JWST prove a highly efficient method to trace the galaxy density in the epoch of reionization.

Aiming to detect planetary companions to young stars with debris disks via the radial velocity method, we observed HD114082 during April 2018 - August 2022 as one of the targets of our RVSPY program (Radial Velocity Survey for Planets around Young stars). We used the FEROS spectrograph, mounted to the MPG/ESO 2.2 m telescope in Chile, to obtain high signal-to-noise spectra and time series of precise radial velocities (RVs). Additionally, we analyzed archival HARPS spectra and TESS photometric data. We used the CERES, CERES++ and SERVAL pipelines to derive RVs and activity indicators and ExoStriker for the independent and combined analysis of the RVs and TESS photometry. We report the discovery of a warm super-Jovian companion around HD114082 based on a 109.8$\pm$0.4 day signal in the combined RV data from FEROS and HARPS, and on one transit event in the TESS photometry. The best-fit model indicates a 8.0$\pm$1.0 Mjup companion with a radius of 1.00$\pm$0.03 Rjup in an orbit with a semi-major axis of 0.51$\pm$0.01 au and an eccentricity of 0.4$\pm$0.04. The companions orbit is in agreement with the known near edge-on debris disk located at about 28 au. HD114082b is possibly the youngest (15$\pm$6 Myr), and one of only three younger than 100 Myr giant planetary companions for which both their mass and radius have been determined observationally. It is probably the first properly model-constraining giant planet that allows distinguishing between hot and cold-start models. It is significantly more compatible with the cold-start model.

Rapid cooling of planetesimal cores has been inferred for several iron meteorite parent bodies based on metallographic cooling rates, and linked to the loss of their insulating mantles during impacts. However, the timing of these disruptive events is poorly constrained. Here, we used the short-lived 107Pd / 107Ag decay system to date rapid core cooling by determining Pd-Ag ages for iron meteorites. We show closure times for the iron meteorites equate to cooling in the timeframe ~7.8 to 11.7 Myr after CAI, and indicate that an energetic inner Solar System persisted at this time. This likely results from the dissipation of gas in the protoplanetary disk, after which the damping effect of gas drag ceases. An early giant planet instability between 5 and 14 Myr after CAI could have reinforced this effect. This correlates well with the timing of impacts recorded by the Pd Ag system for iron meteorites.

We compare the size distributions of {\sl self-organized criticality (SOC)} systems in the solar photosphere and the transition region, using magnetogram data from {\sl Helioseismic and Magnetic Imager (HMI)} and {\sl Interface Region Imaging Spectrograph (IRIS)} data. For each dataset we fit a combination of a Gaussian and a power law size distribution function, which yields information on four different physical processes: (i) photosopheric granulation convection dynamics (explaining the Gaussian random noise distribution in IRIS data); (ii) spicular plage events in the transition region (explaining the power law size distribution in IRIS data); (iii) salt-and-pepper small-scale magnetic structures (explaining the random noise distributions in HMI magnetograms); and (iv) magnetic reconnection processes in flares and nanoflares (explaining the power law size distribution in HMI data). We find a high correlation (CCC=0.97) between IRIS and HMI data. Datasets with magnetic flux balance are found to match the SOC-predicted power law slope $\alpha_F=1.80$ for mean fluxes, which confirms the universality of SOC-inferred flux size distributions, and agrees with the results of Parnell et al.~(2009), $\alpha_F=1.85\pm0.14$.

Studies of massive stars using space photometry have revealed that they commonly exhibit stochastic low-frequency (SLF) variability. This has been interpreted as being caused by internal gravity waves (IGWs) excited at the interface of convective and radiative regions within stellar interiors, such as the convective core or sub-surface convection zones. We aim to compare the properties of SLF variability in massive main-sequence stars observed by the TESS mission determined by different statistical methods, and confirm the correlation between the morphology of SLF variability and a star's location in the HR diagram. From a sample of 30 massive stars observed by TESS, we compare the resultant parameters of SLF variability, in particular the characteristic frequency, obtained from fitting the amplitude spectrum of the light curve with those inferred from fitting the covariance structure of the light curve using the celerite2 Gaussian process (GP) regression software and a damped SHO kernel. We find a difference in the characteristic frequency obtained from the amplitude spectrum fitting and from fitting the co-variance structure of the light curve using a GP regression with celerite2 for only a minority of the considered sample. However, the trends among mass, age and the properties of SLF variability previously reported remain unaffected. GP regression is a useful and novel methodology to efficiently characterise SLF variability in massive stars compared to previous techniques used in the literature. We conclude that the correlation between a star's SLF variability, in particular the characteristic frequency, and its location in the HR diagram is robust for main-sequence massive stars. There also exists a distribution in the stochasticity of SLF variability in massive stars, which indicates that the coherency of SLF variability is also a function of mass and age in massive stars.

We present a non-perturbative method for calculating the abundance of primordial black holes given an arbitrary one-point probability distribution function for the primordial curvature perturbation, $P(\zeta)$. A non-perturbative method is essential when considering non-Gaussianities that cannot be treated using a conventional perturbative expansion. To determine the full statistics of the density field, we relate $\zeta$ to a Gaussian field by equating the cumulative distribution functions. We consider two examples: a specific local-type non-Gaussian distribution arising from ultra slow roll models, and a general piecewise model for $P(\zeta)$ with an exponential tail. We demonstrate that the enhancement of primordial black hole formation is due to the intermediate regime, rather than the far tail. We also show that non-Gaussianity can have a significant impact on the shape of the primordial black hole mass distribution.

Monte Carlo (MC) method played an essential role in many areas of human activity and has found application in many branches of physical science. This paper proposes a computational technique based on MC algorithms to solve Lane-Emden (LE) type equations. We analyze four LE equations arising in astrophysics: the positive and negative indices of the polytropic gas spheres, the isothermal gas sphere, and the white dwarf equation. We calculated eleven models (i.e., eleven LE equations) of the positive index polytropes, nine for the negative index polytrope, the isothermal gas sphere, and the white dwarf equation. Comparing the MC and numerical/analytical models gives good agreement for the four LE equations under study.

Magnetic reconnection in the quiet Sun is a phenomenon that is consistently observed, and it has recently become feasible to simulate via 3D numerical models of realistically stratified and convection-driven reconnection. We aim to illustrate ways by which quiet Sun fields may contribute to solar atmospheric heating via magnetic reconnection that is driven by convective motion. We also aim to compare our stratified model to earlier idealized coronal models in terms of reconnection drivers and topological conditions. We analyzed a simulation of the quiet Sun in which a complex coronal magnetic field is self-consistently driven by the underlying convection. We employed a selection of Lagrangian markers to trace the spatiotemporal behavior of specific magnetic features that are relevant to magnetic reconnection and atmospheric heating. A large-scale reconnection-driven heating event occurs in the simulated corona, in a flattened X-shaped feature characterized by a weak field and high current. Relevant features include a smooth overlying horizontal field, an arcade, and a horizontal flux rope which eventually reconnect with the overlying field, raising coronal plasma temperatures up to 1.47 MK. We find that our results are in good agreement with idealized coronal flare models, which demonstrates that the same physical concepts are valid. We also find that the reconnecting flux rope and arcade are neither formed by any obvious coherent flux emergence, nor by any ordered photospheric motion or flux cancellation. Instead, they seem to develop merely from the self-consistent convective driving of pre-existing tangled field lines. This gradual ordering suggests an inverse cascade of magnetic helicity via smaller reconnection events, located at or above photospheric flux concentrations. We suggest that this case is representative of heating events that may be ubiquitous in the real quiet Sun.

We use a set of about 300 simulated clusters from The Three Hundred Project to calculate their hydrostatic masses and evaluate the associated bias by comparing them with the true cluster mass. Over a redshift range from 0.07 to 1.3, we study the dependence of the hydrostatic bias on redshift, concentration, mass growth, dynamical state, mass, and halo shapes. We find almost no correlation between the bias and any of these parameters. However, there is a clear evidence that the scatter of the mass-bias distribution is larger for low-concentrated objects, high mass growth, and more generically for disturbed systems. Moreover, we carefully study the evolution of the bias of twelve clusters throughout a major-merger event. We find that the hydrostatic-mass bias follows a particular evolution track along the merger process: to an initial significant increase of the bias recorded at the begin of merger, a constant plateaus follows until the end of merge, when there is a dramatic decrease in the bias before the cluster finally become relaxed again. This large variation of the bias is in agreement with the large scatter of the hydrostatic bias for dynamical disturbed clusters. These objects should be avoided in cosmological studies because their exact relaxation phase is difficult to predict, hence their mass bias cannot be trivially accounted for.

We present MargNet, a deep learning-based classifier for identifying stars, quasars and compact galaxies using photometric parameters and images from the Sloan Digital Sky Survey (SDSS) Data Release 16 (DR16) catalogue. MargNet consists of a combination of Convolutional Neural Network (CNN) and Artificial Neural Network (ANN) architectures. Using a carefully curated dataset consisting of 240,000 compact objects and an additional 150,000 faint objects, the machine learns classification directly from the data, minimising the need for human intervention. MargNet is the first classifier focusing exclusively on compact galaxies and performs better than other methods to classify compact galaxies from stars and quasars, even at fainter magnitudes. This model and feature engineering in such deep learning architectures will provide greater success in identifying objects in the ongoing and upcoming surveys, such as Dark Energy Survey (DES) and images from the Vera C. Rubin Observatory.

Recent observations indicate that galactic outflows are ubiquitous in high redshift galaxies, including normal star forming galaxies, quasar hosts, and dusty star forming galaxies (DSFGs). However, the impact of outflows on the evolution of their hosts is still an open question. Here, we analyse the star formation histories (SFH) and galactic outflow properties of galaxies in massive haloes ($10^{12}M_{\odot}<M_{\rm vir} <5\times 10^{12}M_{\odot}$) at $z\gtrsim5.5$ in three zoom-in cosmological simulations from the MassiveFIRE suite, as part of the Feedback In Realistic Environments (FIRE) project. The simulations were run with the FIRE-2 model, which does not include feedback from active galactic nuclei (AGN). The simulated galaxies resemble $z>4$ DSFGs, with SFRs of $\sim 1000\ M_{\odot}\rm yr^{-1}$ and molecular gas masses of $M_{\rm mol}\sim 10^{10}\ M_{\odot}$. However, the simulated galaxies are characterised by higher circular velocities than those observed in high-z DSFGs. The mass loading factors from stellar feedback are of the order of $\sim 0.1$, implying that stellar feedback is inefficient in driving galactic outflows and gas is consumed by star formation on much shorter time-scales than it is expelled from the interstellar medium (ISM). We also find that stellar feedback is highly inefficient in self-regulating star formation in this regime, with an average integrated star formation efficiency (SFE) per dynamical time of $30\%$. Finally, compared to FIRE-2 galaxies hosted in similarly massive haloes at lower redshift, we find lower mass loading factors and higher SFEs in the high redshift sample. We argue that both effects originate from the higher total and gas surface densities that characterise high$-z$ massive systems.

We study the possibility that the axion Peccei-Quinn symmetry is spontaneously broken after the beginning of inflation. This scenario interpolates between pre-inflationary and post-inflationary axion cosmology with significant phenomenological differences from both. Since the axion is not present at the early stages of inflation large inflationary fluctuations are produced only at scales not constrained by CMB, avoiding the strongest isocurvature constraints. The energy density in isocurvature perturbations at short scales however can be comparable with the adiabatic contribution from misalignment. These large overdensities can lead to the formation of axion mini-clusters and also to constraints from Lyman-$\alpha$ forest and future CMB spectral distortion measurements. If Peccei-Quinn symmetry is broken during the first O(25) e-foldings of inflation no axions are produced from the string network but contributions from the annihilation of domain walls can further boost the abundance. This scenario is minimally realized if the Hubble scale during inflation drops below the Peccei-Quinn symmetry breaking scale but other realizations are possible.

The positive hyperon production asymmetries that have been measured at LHC are real demonstrations of string junction role in the baryon charge transfer at baryon production in HE proton-proton interactions. In order to invent the neutral heavy particle with zero baryon charge as a candidate for Dark Matter, it is necessary to turn back to the times of dual topological representation of hadron physics and pomeron theory. The topological presentation of pomeron exchange at the HE proton-proton collision is cylinder that is covered with quark-gluon net. Taking into account that the junction of three gluons (SJ) has the positive baryon number, as well as the antijunction {antiSJ) is of negative baryon charge, the neutral self-connected baryonium configuration with zero baryon charge is torus that is covered by a discrete number of hexagons with 3 string junction and 3 antijunction vertices each. The possibilities to observe such baryonium tori are discussed in this paper for three cases: collisions at p-p colliders, the cosmic particles interactions with atmosphere and in the relativistic jets at the Active Galaxy Nuclei as well.

In recent work we showed that the relic dark matter abundance of QCD axions can be altered when the Peccei-Quinn (PQ) field is coupled to very light scalar/s, rendering the effective axion mass dynamical in the early universe. In this work we develop this framework further, by introducing a new extended symmetry group to protect the new particles' mass. We find that with a new global $SO(N)$ symmetry, with large $N$, we can not only account for the lightness of the new scalars, but we can reduce the axion relic abundance in a technically natural way. This opens up the possibility of large PQ scales, including approaching the GUT scale, and still naturally producing the correct relic abundance of axions. Also, in these models the effective PQ scale is relatively small in the very early universe, and so this can help towards alleviating the isocurvature problem from inflation. Furthermore, instead of possible over-closure from cosmic strings, the extended symmetry implies the formation of non-topological textures which provide a relatively small abundance.

Cold interstellar gas systems have been used to constrain dark matter (DM) models by the condition that the heating rate from DM must be lower than the astrophysical cooling rate of the gas. Following the methodology of Wadekar and Farrar (2021), we use the interstellar medium of a gas-rich dwarf galaxy, Leo T, and a Milky Way-environment gas cloud, G33.4-8.0 to constrain DM. Leo T is a particularly strong system as its gas can have the lowest cooling rate among all the objects in the late Universe (owing to the low volume density and metallicity of the gas). Milky Way clouds, in some cases, provide complementary limits as the DM-gas relative velocity in them is much larger than that in Leo T. We derive constraints on the following scenarios in which DM can heat the gas: $(i)$ interaction of axions with hydrogen atoms or free electrons in the gas, $(ii)$ deceleration of relic magnetically charged DM in gas plasma, $(iii)$ dynamical friction from compact DM, $(iv)$ hard sphere scattering of composite DM with gas. Our limits are complementary to DM direct detection searches. Detection of more gas-rich low-mass dwarfs like Leo T from upcoming 21cm and optical surveys can improve our bounds.

Extreme conditions present in the interiors of the core-collapse supernovae make neutrino-neutrino interactions not only feasible but dominant in specific regions, leading to the non-linear evolution of the neutrino flavor. Results obtained when such collective neutrino oscillations are treated in the mean-field approximation deviate from the results using the many-body picture because of the ignored quantum correlations. We present the first three flavor many-body calculations of the collective neutrino oscillations. The entanglement is quantified in terms of the entanglement entropy and the components of the polarization vector. We propose a qualitative measure of entanglement in terms of flavor-lepton number conserved quantities. We find that in the cases considered in the present work, the entanglement can be underestimated in two flavor approximation. The dependence of the entanglement on mass ordering is also investigated. We also explore the mixing of mass eigenstates in different mass orderings.

We investigate Natural Inflation with non-minimal coupling to gravity, characterized either by a quadratic or a periodic term, within the Warm Inflation paradigm during the slow roll stage, in both strong and weak dissipation limits, and show that it can accommodate the spectral index $n_s$ and tensor-to-scalar ratio $r$ observables given by Planck 2018 constraints, albeit with a too small value of the e-folding number to solve the horizon problem, providing thus only a partial solution to Natural Inflation issues.

We compute bounds and forecasts on screened modified gravity theories, specialising to the chameleon model in Casimir force experiments. In particular, we investigate the classical interaction between a plate and sphere subject to a screened interaction of the chameleon type. We compare numerical simulations of the field profile and the classical pressure exerted on the sphere to analytical approximations for these non-linear field theories. In particular, we focus on the proximity force approximation (PFA) and show that, within the range of sphere sizes $R$ and plate-sphere distance $D$ simulated numerically, the PFA does not reproduce the numerical results. This differs from the case of linear field theories such as Newtonian gravity and a Yukawa model where the PFA coincides with the exact results. We show that for chameleon theories, the screening factor approximation (SFA) whereby the sphere is modelled as a screened sphere embedded in the external field due to the plates, fares better and can be used in the regime $D\gtrsim R$ to extract constraints and forecasts from existing and forthcoming data. In particular, we forecast that future Casimir experiments would corroborate the closing of the parameter space for the simplest of chameleon models at the dark energy scale.

Coronal mass ejections (CMEs) can trigger geomagnetic storms and induce geoelectric currents that degrade the performance of terrestrial power grid operations; in particular, CMEs are known for causing large-scale outages in electrical grids. Submarine internet cables are powered through copper conductors spanning thousands of kilometers and are vulnerable to damage from CMEs, raising the possibility of a large-scale and long-lived internet outage. To better understand the magnitude of these risks, we monitor voltage changes in the cable power supply of four different transoceanic cables during time periods of high solar activity. We find a strong correlation between the strength of the high-frequency geomagnetic field at the landing sites of the systems and the line voltage change. We also uncover that these two quantities exhibit a near-linear power law scaling behavior that allows us to estimate the effects of once-in-a-century CME events. Our findings reveal that long-haul submarine cables, regardless of their length and orientation, will not be damaged during a solar superstorm, even one as large as the 1859 Carrington event.

The problem of detecting a sinusoidal signal with randomly varying frequency has a long history. It is one of the core problems in signal processing, arising in many applications including, for example, underwater acoustic frequency line tracking, demodulation of FM radio communications, laser phase drift in optical communications and, recently, continuous gravitational wave astronomy. In this paper we describe a Markov Chain Monte Carlo based procedure to compute a specific detection posterior density. We demonstrate via simulation that our approach results in an up to $25$ percent higher detection rate than Hidden Markov Model based solutions, which are generally considered to be the leading techniques for these problems.

By keeping the cosmological principle i.e., an isotropic and homogeneous universe, we consider the cosmology of a vector-tensor theory of gravitation known as the \textit{bumblebee} model. In this model a single Lorentz-violating timelike vector field with a nonzero vacuum expectation value (VEV) couples to the Ricci tensor and scalar, as well. Taking the ansatz $B(t)\sim t^\beta$ for the time evolution of the vector field we derive the relevant dynamic equations of the Universe, where $\beta$ is a free parameter. In particular, by employing observational data coming from the Big Bang Nucleosynthesis (BBN) and the matter-antimatter asymmetry in the Baryogenesis era, we impose some constraints on the VEV of the bumblebee timelike vector field i.e., $\xi b^2$, and the exponent parameter $\beta$. The former and the latter limit the size of Lorentz violation, and the rate of the time evolution of the background Lorentz-violating bumblebee field, respectively.

One of the most pressing questions for modern physics is the nature of dark matter (DM). Several efforts have been made to model this elusive kind of matter. The largest fraction of DM cannot be made of any of the known particles of the Standard Model (SM). We focus on brane world theory as a prospective framework for DM candidates beyond the SM of particle physics. The new degrees of freedom that appear in flexible brane world models, corresponding to brane fluctuations, are called branons. They behave as weakly interacting massive particles (WIMPs), which are one of the most favored candidates for DM. We present a multi-target DM search in dwarf spheroidal galaxies for branon DM annihilation signatures with the ground-based gamma-ray telescope MAGIC leading to the most constraining branon DM limits in the TeV mass range.

In this paper, we would like to figure out whether a k-inflation model admits the Bianchi type I metric as its inflationary solution under a constant-roll condition in the presence of the supergravity motivated coupling between scalar and vector fields, $f^2(\phi)F_{\mu\nu}F^{\mu\nu}$. As a result, some novel anisotropic inflationary solutions are shown to appear along with a power-law one in this scenario. Furthermore, these solutions are numerically confirmed to be attractive, in contrast to the prediction of the Hawking's cosmic no-hair conjecture.

We demonstrate analytically and numerically that "optimized partial dressing" (OPD) thermal mass resummation, which uses gap equation solutions inserted into the tadpole, efficiently tames finite temperature perturbation theory calculations of the effective thermal potential, without necessitating use of the high-temperature approximation. An analytical estimate of the scale dependence for OPD resummation, standard Parwani Daisy-resummation and dimensional reduction shows that OPD has similar scale dependence to dimensional reduction, greatly improved over Parwani resummation. We also elucidate how to construct and solve the gap equation for realistic numerical calculations, and demonstrate OPD's improved accuracy and precision for a toy scalar model. An example of the physical significance of OPD's improved accuracy is the maximal gravitational wave amplitude that a model is capable of generating, which Parwani resummation underestimates by two orders of magnitude. This highlights the need to bring theoretical uncertainties under control even when analysing broad features of a model. Given the simplicity of the OPD compared to two loop dimensional reduction, as well as the ease with which this scheme handles departures from the high temperature expansion, we argue this scheme has great potential in analyzing the parameter space of realistic BSM models.

The investigation of relativistic effects in the most general way requires a unified treatment of light propagation in cosmology. This goal can be achieved with the new interpretation of the geodesic deviation equation in terms of the bilocal geodesic operators (BGO). The BGO formalism extends the standard formulation, providing a unified framework to describe all possible optical phenomena due to the interaction between light and spacetime curvature. In my dissertation, I present {\tt BiGONLight}, a {\tt Mathematica} package that applies the BGO formalism to study light propagation in numerical relativity. The package encodes the 3+1 bilocal geodesic operators framework as a collection of {\tt Mathematica} functions. The inputs are the spacetime metric plus the kinematics of the observer and the source in the form of the 3+1 quantities, which may come directly from a numerical simulation or can be provided by the user as analytical components. These data are then used for ray tracing and computing the BGOs in a completely general way, i.e. without relying on symmetries or specific coordinate choices. The primary purpose of the package is the computation of optical observables in arbitrary spacetimes. The uniform theoretical framework of the BGO formalism allows for the extraction of multiple observables within a single computation, while the {\tt Wolfram} language provides a flexible computational framework that makes the package highly adaptable to perform both numerical and analytical studies of light propagation.

We consider generic freeze-in processes for generation of Dark Matter, together with the consequent re-thermalization of the Standard Model fluid. We find that Dark Matter inherits the Standard Model adiabatic inhomogeneities on the cosmological scales probed by current observations, that were super-horizon during freeze-in. Thereby, freeze-in satisfies the bounds on iso-curvature perturbations.