Papers for Tuesday, Mar 28 2023

The chemical fingerprints of the first stars are retained within the photospheres of ancient unevolved metal-poor stars. A significant fraction of these stellar fossils is represented by stars known as Carbon-Enhanced Metal-Poor (CEMP), $\rm [C/Fe]>+0.7$ and $\rm [Fe/H]<-2$, which are likely imprinted by low-energy primordial supernovae. These CEMP stars are largely observed in the Galactic halo and ultra-faint dwarf galaxies, with values reaching $\rm [C/Fe]=+4.5$. The Galactic bulge is predicted to host the oldest stars, but it shows a striking dearth of CEMP stars with $\rm [C/Fe]\gtrsim +2.0$. Here we explore the possible reasons for this anomaly by performing a statistical analysis of the observations of metal-poor stars in combination with the predictions of $\Lambda$CDM models. We suggest that the dearth of CEMP stars with high $\rm [C/Fe]$ is not due to the low statistics of observed metal-poor stars but is the result of the different formation process of the bulge. $N$-body simulations show that the first star-forming halos which end up in the bulge are characterized by the highest star-formation rates. These rates enable the formation of rare massive first stars exploding as pair-instability supernovae (PISNe), which wash out the signature of primordial faint supernovae. We demonstrate that the mean $\rm [C/Fe]$ of first stars polluted environments decreases with the increasing contribution of PISNe. We conclude that the dearth of CEMP stars in the Galactic bulge indirectly probes the existence of elusive PISNe, and propose a novel method which exploits this lack to constrain the mass distribution of the first stars.

Atmospheric metal enrichment (i.e., elements heavier than helium, also called "metallicity") is a key diagnostic of the formation of giant planets. The giant planets of the solar system exhibit an inverse relationship between mass and both their bulk metallicities and atmospheric metallicities. Extrasolar giant planets also display an inverse relationship between mass and bulk metallicity. However, there is significant scatter in the relationship and it is not known how atmospheric metallicity correlates with either planet mass or bulk metallicity. Here we show that the Saturn-mass exoplanet HD 149026b has an atmospheric metallicity 59 - 276 times solar (at 1 $\sigma$), which is greater than Saturn's atmospheric metallicity of ~7.5 times solar at >4 $\sigma$ confidence. This result is based on modeling CO$_2$ and H$_2$O absorption features in the thermal emission spectrum of the planet measured by JWST. HD 149026b is the most metal-rich giant planet known, with an estimated bulk heavy element abundance of 66 $\pm$ 2% by mass. We find that the atmospheric metallicities of both HD 149026b and the solar system giant planets are more correlated with bulk metallicity than planet mass.

We study the kinematics of stars both at their formation and today within 14 Milky Way (MW)-mass galaxies from the FIRE- 2 cosmological zoom-in simulations. We quantify the relative importance of cosmological disk settling and post-formation dynamical heating. We identify three eras: a Pre-Disk Era (typically >8 Gyr ago), when stars formed on dispersion-dominated orbits; an Early-Disk Era (~ 8 - 4 Gyr ago), when stars started to form on rotation-dominated orbits but with high velocity dispersion, sigma_form; and a Late-Disk Era (< 4 Gyr ago), when stars formed with low sigma_form. sigma_form increased with time during the Pre-Disk Era, peaking ~ 8 Gyr ago, then decreased throughout the Early-Disk Era as the disk settled and remained low throughout the Late-Disk Era. By contrast, the velocity dispersion measured today, sigma_now, increases monotonically with age because of stronger post-formation heating for Pre-Disk stars. Importantly, most of sigma_now was in place at formation, not added post-formation, for stars younger than ~ 10 Gyr. We compare the evolution of the three velocity components: at all times, sigma_R,form > sigma_phi,form > sigma_Z,form. Post-formation heating primarily increased sigma_R at ages < 4 Gyr but acted nearly isotropically for older stars. The lookback time that the disk began to settle correlates with its dynamical state today: earlier-settling galaxies currently form colder disks. Young stars in FIRE-2 are kinematically hotter than the MW but broadly agree with M31 and M33. Including stellar cosmic-ray feedback does not significantly change the amount of disk rotational support at fixed stellar mass.

Warm Dark Matter (WDM) particles with masses ($\sim$ kilo electronvolt) offer an attractive solution to the small-scale issues faced by the Cold Dark Matter (CDM) paradigm. The delay of structure formation in WDM models and the associated dearth of low-mass systems at high-redshifts makes this an ideal time to revisit WDM constraints in light of the unprecedented data-sets from the James Webb Space Telescope (JWST). Developing a phenomenological model based on the halo mass functions in CDM and WDM models, we calculate high-redshift ($z \gt 6$) the stellar mass functions (SMF) and the associated stellar mass density (SMD) and the maximum stellar mass allowed in a given volume. We find that: (i) WDM as light as 1.5 keV is already disfavoured by the low-mass end of the SMF (stellar mass $M_* \sim 10^7 \rm{M_\odot}$) although caution must be exerted given the impact of lensing uncertainties; (ii) 1.5 keV WDM models predict SMD values that show a steep decrease from $10^{8.8}$ to $10^{2} ~{\rm M_\odot ~cMpc^{-3}}$ from $z \sim 4$ to 17 for $M_* \gt 10^8 \rm{M_\odot}$; (iii) the 1.5 keV WDM model predicts a sharp and earlier cut-off in the maximum stellar masses for a given number density (or volume) as compared to CDM or heavier WDM models. For example, with a number density of $10^{-3} \rm {cMpc^{-3}}$, 1.5 (3) KeV WDM models do not predict bound objects at $z \gt 12$ (18). Forthcoming JWST observations of multiple blank fields can therefore be used as a strong probe of WDM at an epoch inaccessible by other means.

We present the discovery of a new highly compact quadruple star system, TIC 219006972, consisting of two eclipsing binary stars with orbital periods of 8.3 days and 13.7 days, and an outer orbital period of only 168 days. This period is a full factor of 2 shorter than the quadruple with the shortest outer period reported previously, VW LMi, where the two binary components orbit each other every 355 days. The target was observed by TESS in Full-Frame Images in sectors 14-16, 21-23, 41, 48 and 49, and produced two sets of primary and secondary eclipses. These show strongly non-linear eclipse timing variations (ETVs) with an amplitude of $\sim$0.1 days, where the ETVs of the primary and secondary eclipses, and of the two binaries are all largely positively correlated. This highlights the strong dynamical interactions between the two binaries and confirms the compact quadruple configuration of TIC 219006972. The two eclipsing binaries are nearly circular whereas the quadruple system has an outer eccentricity of about 0.25. The entire system is nearly edge-on, with a mutual orbital inclination between the two eclipsing binary star systems of about 1 degree.

A merger of binary neutron stars creates heavy unstable elements whose radioactive decay produces a thermal emission known as a kilonova. In this paper, we predict the photometric and polarimetric behaviour of this emission by performing 3-D Monte Carlo radiative transfer simulations. In particular, we choose three hydrodynamical models for merger ejecta, two including jets with different luminosities and one without a jet structure, to help decipher the impact of jets on the light curve and polarimetric behaviour. In terms of photometry, we find distinct color evolutions across the three models. Models without a jet show the highest variation in light curves for different viewing angles. In contrast, to previous studies, we find models with a jet to produce fainter kilonovae when viewed from orientations close to the jet axis, compared to a model without a jet. In terms of polarimetry, we predict relatively low levels (<~0.3-0.4%) at all orientations that, however, remain non-negligible until a few days after the merger and longer than previously found. Despite the low levels, we find that the presence of a jet enhances the degree of polarization at wavelengths ranging from 0.25 to 2.5\micron, an effect that is found to increase with the jet luminosity. Thus, future photometric and polarimetric campaigns should observe kilonovae in blue and red filters for a few days after the merger to help constrain the properties of the ejecta (e.g. composition) and jet.

Moons orbiting exoplanets ("exomoons") may hold clues about planet formation, migration, and habitability. We investigate the plausibility of exomoons orbiting the temperate ($T_\text{eq}=294$ K) giant ($R = 9.2$ $\text{R}_\oplus$) planet HIP 41378 f. Previous studies have suggested that HIP 41378 f has a low apparent bulk density of $0.09\,\text{g}\,\text{cm}^{-3}$ and a flat near-infrared transmission spectrum, suggesting that it may possess circumplanetary rings. Given the planet's long orbital period ($P\approx1.5$ yr), it may also host a large exomoon. Here, we consider a hypothetical exomoon orbiting HIP 41378 f with the same satellite-to-planet mass ratio as the Moon-Earth system and assess its orbital stability using a suite of N-body and tidal migration simulations. We find that satellites up to this size are largely stable against Hill sphere escape and collisions with the host planet, consistent with theoretical stability limits determined by previous studies. We then simulate the expected transit signal from the exomoon and show that current transit observations likely cannot constrain the presence of exomoons orbiting HIP 41378 f, though future observations may be capable of detecting exomoons in other systems. Finally, we simulate the combined transmission spectrum of HIP 41378 f and an exomoon with a low-metallicity atmosphere, and show that the total effective spectrum is contaminated at the $\sim$10 ppm level. Our work not only demonstrates the feasibility of exomoons orbiting HIP 41378 f, but also shows that large exomoons may be a source of uncertainty in future high-precision measurements of exoplanet systems.

Stellar atmospheres separate the hot and dense stellar interiors from the emptiness of space. Radiation escapes from the outermost layers of a star, carrying direct physical information. Underneath the atmosphere, the very high opacity keeps radiation thermalized and resembling a black body with the local temperature. In the atmosphere the opacity drops, and radiative energy leaks out, which is redistributed in wavelength according to the physical processes by which matter and radiation interact, in particular photoionization. In this article, I will evaluate the role of photoionization in shaping the stellar energy distribution of stars. To that end, I employ simple, state-of-the-art plane-parallel model atmospheres and a spectral synthesis code, dissecting the effects of photoionization from different chemical elements and species, for stars of different masses in the range of 0.3 to 2 M$_{\odot}$. I examine and interpret the changes in the observed spectral energy distributions of the stars as a function of the atmospheric parameters. The photoionization of atomic hydrogen and H$^-$ are the most relevant contributors to the continuum opacity in the optical and near-infrared regions, while heavier elements become important in the ultraviolet region. In the spectra of the coolest stars (spectral types M and later), the continuum shape from photoionization is no longer recognizable due to the accumulation of lines, mainly from molecules. These facts have been known for a long time, but the calculations presented provide an updated quantitative evaluation and insight into the role of photoionization on the structure of stellar atmospheres.

The radio jets of an active galactic nucleus (AGN) can heat up the gas around a host galaxy and quench star formation activity. The presence of a radio jet could be related to the evolutionary path of the host galaxy and may be imprinted in the morphology and kinematics of the galaxy. In this work, we use data from the Sloan Digital Sky Survey's Mapping Nearby Galaxies at Apache Point Observatory survey and the Low Frequency Array (LOFAR) Two-Metre Sky Survey as well as the National Radio Astronomy Observatory (NRAO) the Karl G. Jansky Very Large Array (VLA) Sky Survey and the Faint Images of the Radio Sky at Twenty Centimeter survey. We combine these integral field spectroscopic data and radio data to study the link between stellar kinematics and radio AGNs. We find that the luminosity-weighted stellar angular momentum $\lambda_{Re}$ is tightly related to the range of radio luminosity and the fraction of radio AGNs F radio present in galaxies, as high-luminosity radio AGNs are only in galaxies with a small $\lambda_{Re}$, and the $F_{radio}$ at a fixed stellar mass decreases with $\lambda_{Re}$. These results indicate that galaxies with stronger random stellar motions with respect to the ordered motions might be better breeding grounds for powerful radio AGNs. This would also imply that the merger events of galaxies are important in the triggering of powerful radio jets in our sample.

The emission of gamma-ray burst (GRB) 221009A at 18 TeV has been detected by the large high-altitude air shower observatory (LHAASO). We suggest jitter radiation as a possible explanation for the TeV emission for this energetic GRB. In our scenario, the radiation field is linked to the perturbation field, and the perturbation field is dominated by kinetic turbulence. Kinetic turbulence takes a vital role in both magnetic field generation and particle acceleration. The jitter radiation can reach the TeV energy band when we consider either electron cooling or Landau damping. We further suggest that the jitter radiation in the very high-energy band is coherent emission. Our modeling results can be constrained by the observational results of GRB 221009A in the TeV energy band. This radiation mechanism is expected to have wide applications in the high-energy astrophysical research field.

The determination of the absolute and relative position of a spacecraft is critical for its operation, observations, data analysis, scientific studies, as well as deep space exploration in general. A spacecraft that can determine its own absolute position autonomously may perform more than that must rely on transmission solutions. In this work, we report an absolute navigation accuracy of $\sim$ 20 km using 16-day Crab pulsar data observed with $Fermi$ Gamma ray Burst Monitor (GBM). In addition, we propose a new method with the inverse process of the triangulation for joint navigation using repeated bursts like that from the magnetar SGR J1935+2154 observed by the Gravitational wave high-energy Electromagnetic Counterpart All-sky Monitor (GECAM) and GBM.

We present here the East Asia to Italy Nearly Global VLBI (EATING VLBI) project. How this project started and the evolution of the international collaboration between Korean, Japanese, and Italian researchers to study compact sources with VLBI observations is reported. Problems related to the synchronization of the very different arrays and technical details of the telescopes involved are presented and discussed. The relatively high observation frequency (22 and 43 GHz) and the long baselines between Italy and East Asia produced high-resolution images. We present example images to demonstrate the typical performance of the EATING VLBI array. The results attracted international researchers and the collaboration is growing, now including Chinese and Russian stations. New in progress projects are discussed and future possibilities with a larger number of telescopes and a better frequency coverage are briefly discussed herein.

We have performed the near-infrared photometric monitoring observations of two TeV gamma-ray binaries with O-stars (LS 5039 and 1FGL J1018.6-5856), using IRSF/SIRIUS at SAAO, in order to study the stellar parameters and their perturbations caused by the binary interactions. The whole orbital phase was observed multiple times and no significant variabilities including orbital modulations are detected for both targets. Assuming that the two systems are colliding wind binaries, we estimate the amplitude of flux variation caused by the difference in the optical depth of O-star wind at inferior conjunction, where the star is seen through the cavity created by pulsar wind, and other orbital phases without pulsar-wind intervention. The derived amplitude is <0.001 mag, which is about two orders of magnitude smaller than the observed upper limit. Also using the upper limits of the near-infrared variability, we for the first time obtain the upper limit of the dust formation rate resulting from wind-wind collision in O-star gamma-ray binaries.

Tidally locked worlds provide a unique opportunity for constraining the probable climates of certain exoplanets. They are unique in that few exoplanet spin and obliquity states are known or will be determined in the near future: both of which are critical in modeling climate. A recent study shows the dynamical conditions present in the TRAPPIST-1 system make rotation and large librations of the substellar point possible for these planets, which are usually assumed to be tidally locked. We independently confirm the tendency for planets in TRAPPIST-1-like systems to sporadically transition from tidally locked libration to slow rotation using N-body simulations. We examine the nature and frequency of these spin states to best inform energy balance models which predict the temperature profile of the planet's surface. Our findings show that tidally locked planets with sporadic rotation are able to be in both long-term persistent states and chaotic states: where rapid transitions between behaviors are present. Quasi-stable spin regimes, where the planet exhibits one spin behavior for up to hundreds of millennia, are likely able to form stable climate systems while the spin behavior is constant. 1D energy balance models show that tidally locked planets with sporadic rotation around M-dwarfs will experience a relatively small change in substellar temperature due to the lower albedo of ice in an infrared dominant stellar spectrum. The exact effects of large changes in temperature profiles on these planets as they rotate require more robust climate models, like 3D global circulation models, to better examine.

MASTER OT J004527.52+503213.8 (hereafter MASTER J004527) is a dwarf nova discovered by the MASTER project in 2013. At 18:20 UTC on 24 October 2020, brightening of this object was reported to vsnet-alert (24843 by Denisenko). This was the second report of a superoutburst after its discovery. Photometric observations were made using the 23.5-cm Schmidt-Cassegrain telescope at Okayama University of Science observatory soon after the alert through 4 November 2020. In this work, we present the photometric data from our observation, and the analysis of the light curves of MASTER J004527 during the 2020 outburst. We propose a method to determine the period of superhumps by polynomial fitting, which can be applied to a light curve with many missing data. In addition to our own data, we incorporate other all sky survey data of the outburst to better understand the properties of the superhumps. Based on our observations, we conclude that MASTER J004527 is an SU UMa-type dwarf nova, since no early superhumps occurred.

Context: Given the current big data era in Astronomy, machine learning based methods have being applied over the last years to identify or classify objects like quasars, galaxies and stars from full sky photometric surveys. Aims: Here we systematically evaluate the performance of Random Forests (RF) in classifying quasars using either magnitudes or colours, both from broad and narrow-band filters, as features. Methods: The working data consists of photometry from the ALHAMBRA Gold Catalogue that we cross-matched with the Sloan Digital Sky Survey (SDSS) and with the Million Quasars Catalogue (Milliquas) for objects labelled as quasars, galaxies or stars. A RF classifier is trained and tested to evaluate the effect on final accuracy and precision of varying the free parameters and the effect of using narrow or broad-band magnitudes or colours. Results: Best performances of the classifier yielded global accuracy and quasar precision around 0.9. Varying model free parameters (within reasonable ranges of values) has no significant effects on the final classification. Using colours instead of magnitudes as features results in better performances of the classifier, especially using colours from the ALHAMBRA Survey. Colours that contribute the most to the classification are those containing the near-infrared $JHK$ bands.

With dedicated exoplanet surveys underway for multiple extreme precision radial velocity (EPRV) instruments, the near-future prospects of RV exoplanet science are promising. These surveys' generous time allocations are expected to facilitate the discovery of Earth analogs around bright, nearby Sun-like stars. But survey success will depend critically on the choice of observing strategy, which will determine the survey's ability to mitigate known sources of noise and extract low-amplitude exoplanet signals. Here, we present an analysis of the Fisher information content of simulated EPRV surveys, accounting for the most recent advances in our understanding of stellar variability on both short and long timescales (i.e., oscillations and granulation within individual nights, and activity-induced variations across multiple nights). In this analysis, we capture the correlated nature of stellar variability by parameterizing these signals with Gaussian Process kernels. We describe the underlying simulation framework as well as the physical interpretation of the Fisher information content, and we evaluate the efficacy of EPRV survey strategies that have been presented in the literature. We explore and compare strategies for scheduling observations over various timescales and we make recommendations to optimize survey performance for the detection of Earth-like exoplanets.

Planet formation imprints signatures on the physical structures of disks. In this paper, we present high-resolution ($\sim$50 mas, 8 au) Atacama Large Millimeter/submillimeter Array (ALMA) observations of 1.3 mm dust continuum and CO line emission toward the disk around the M3.5 star 2MASS J04124068+2438157. The dust disk consists only of two narrow rings at radial distances of 0.47 and 0.78 arcsec ($\sim$70 and 116 au), with Gaussian $\sigma$ widths of 5.6 and 8.5 au, respectively. The width of the outer ring is smaller than the estimated pressure scale height by $\sim25\%$, suggesting dust trapping in a radial pressure bump. The dust disk size, set by the location of the outermost ring, is significantly larger (by $3\sigma$) than other disks with similar millimeter luminosity, which can be explained by an early formation of local pressure bump to stop radial drift of millimeter dust grains. After considering the disk's physical structure and accretion properties, we prefer planet--disk interaction over dead zone or photoevaporation models to explain the observed dust disk morphology. We carry out high-contrast imaging at $L'$ band using Keck/NIRC2 to search for potential young planets, but do not identify any source above $5\sigma$. Within the dust gap between the two rings, we reach a contrast level of $\sim$7 mag, constraining the possible planet below $\sim$2--4 $M_{\rm Jup}$. Analyses of the gap/ring properties suggest a $\sim$Saturn mass planet at $\sim$90 au is likely responsible for the formation of the outer ring, which can be potentially revealed with JWST.

The possible time variation of the fundamental constants of nature has been an active subject of research in modern physics. In this paper, we propose a new method to investigate such possible time variation of the speed of light $c$ using the updated Hubble diagram of high-redshift standard candles including Type Ia Supernovae (SNe Ia) and high-redshift quasars (based on UV-X relation). Our findings show that the SNe Ia Pantheon sample, combined with currently available sample of cosmic chronometers, would produce robust constraints on the speed of light at the level of $c/c_0=1.03\pm0.03$. For the Hubble diagram of UV+X ray quasars acting as a new type of standard candles, we obtain $c/c_0=1.19\pm0.07$. Therefore, our results confirm that there is no strong evidence for the deviation from the constant speed of light up to $z\sim 2$. Moreover, we discuss how our technique might be improved at much higher redshifts ($z\sim5$), focusing on future measurements of the acceleration parameter $X(z)$ with gravitational waves (GWs) from binary neutron star mergers. In particular, in the framework of the second-generation space-based GW detector, DECi-hertz Interferometer Gravitational-wave Observatory (DECIGO), the speed of light is expected to be constrained with the precision of $\Delta{c}/c=10^{-3}$.

A large number of fast radio bursts (FRBs) detected with the CHIME telescope enable us to investigate their energy distributions in different redshift intervals, incorporating with the consideration of the selection effects of CHIME. As a result, we obtain a non-evolving energy function (EF), which is in a power law form with a low-energy exponential cutoff, for the high-energy FRBs (HEFRBs) of energies $E\gtrsim2\times0^{38}$ erg. On the contrary, the energy distribution of the low-energy FRBs (LEFRBs) obviously cannot be described by the same EF. Including the lowest dispersion measures (DMs) samples, the LEFRBs are concentrated towards the Galactic plane and their latitude distribution is similar to that of Galactic rotational radio transients (RRATs). These indications hint that LEFRBs might compose a special type of RRATs with relatively higher DMs and energies (i.e., $\sim10^{28-31}$ erg for a reference distance of $\sim10$ kpc if they belong to the Milky Way). Finally, we revisit the redshift-dependent event rate of HEFRBs and confirm that they could be produced by the remnants of cosmological compact binary mergers.

The previous decade saw the discovery of the first four known interstellar objects due to advances in astronomical viewing equipment. Future sky surveys with greater sensitivity will allow for more frequent detections of such objects, including increasingly small objects. We consider the capabilities of the Legacy Survey of Space and Time (LSST) of the Vera C. Rubin Observatory to detect interstellar objects of small sizes during its period of operation over the next decade. We use LSST's detection capabilities and simulate populations of interstellar objects in the range of 1-50m in diameter to calculate the expected number of small interstellar objects that will be detected. We use previous detections of interstellar objects to calibrate our object density estimates. We also consider the impact of the population's albedo on detection rates by considering populations with two separate albedo distributions: a constant albedo of 0.06 and an albedo distribution that resembles near earth asteroids. We find that the number of detections increases with the diameter over the range of diameters we consider. We estimate a detection rate of up to a small ISO every two years of LSST's operation with an increase by a factor of ten for future surveys that extend a magnitude deeper.

We present a routinized and reliable method to obtain source catalogs from the $\textit{Nuclear Spectroscopic Telescope Array}$ ($\textit{NuSTAR}$) extragalactic surveys of the Extended $\textit{Chandra}$ Deep Field-South (E-CDF-S) and $\textit{Chandra}$ Deep Field-North (CDF-N). The $\textit{NuSTAR}$ E-CDF-S survey covers a sky area of $\approx30'\times30'$ to a maximum depth of $\sim$ 230 ks corrected for vignetting in the 3--24 keV band, with a total of 58 sources detected in our E-CDF-S catalog; the $\textit{NuSTAR}$ CDF-N survey covers a sky area of $\approx7'\times10'$ to a maximum depth of $\sim$ 440 ks corrected for vignetting in the 3--24 keV band, with a total of 42 sources detected in our CDF-N catalog that is produced for the first time. We verify the reliability of our two catalogs by crossmatching them with the relevant catalogs from the $\textit{Chandra}$ X-ray observatory, and find that the fluxes of our $\textit{NuSTAR}$ sources are generally consistent with that of their $\textit{Chandra}$ counterparts. Our two catalogs are produced following the exactly same method and made publicly available, thereby providing a uniform platform that facilitates further studies involving these two fields. Our source-detection method provides a systematic approach for source cataloging in other $\textit{NuSTAR}$ extragalactic surveys.

The present study collects 456 new times of maximum light of the classical Cepheid RT Aur, covering the period from 1897 to 2022. The O-C diagram resulting from these observations shows that the period given by the GCVS has to be corrected. It results that no strong period variation is found. However, the observed O-C residuals show a long term periodic trend. In the hypothesis of RT Aur being in a binary system, an orbit cannot be deduced from the available astrophysical data.

Strong gravitationally lensed arcs and arclets produced by the mass distribution in galaxy clusters have been observationally detected for several decades now. These strong lensing constraints provided high-fidelity mass models for cluster lenses that include a detailed census of the substructure down to $10^{9-10}\,\mathrm{M}_\odot$. Optimizing lens models, where the cluster mass distribution is modeled by a smooth component and subhalos associated with the locations of individual cluster galaxies, has enabled deriving the subhalo mass function, providing important constraints on the nature and granularity of dark matter. In this work, we explore and present a novel method to detect and measure individual perturbers (subhalos, line-of-sight halos, and wandering supermassive black holes) by exploiting their proximity to highly distorted lensed arcs in galaxy clusters, and by modeling the local lensing distortions with curved arc bases. This method offers the possibility of detecting individual low-mass perturber subhalos in clusters and halos along the line-of-sight down to a mass resolution of $10^8\, \mathrm{M}_\odot$. We quantify our sensitivity to low-mass perturbers with masses $M \sim 10^{7-9}\,\mathrm{M}_\odot$ in clusters with masses $M \sim 10^{14-15}\mathrm{M}_\odot$, by creating realistic mock data. Using three lensed images of a background galaxy in the cluster SMACS J0723, as seen by the $\textit{James Webb Space Telescope}$, we study the retrieval of the properties of potential perturbers with masses $M = 10^{7-9}\,\mathrm{M}_\odot$. From the derived posterior probability distributions for the perturber, we constrain its concentration, redshift, and ellipticity. By allowing us to probe lower-mass substructures, the use of the curved arc bases can lead to powerful constraints on the nature of dark matter as discrimination between dark matter models appears on smaller scales.

We present a method for testing the predictions of hadronic interaction models and improving their consistency with observed two-dimensional distributions of the depth of shower maximum, $X_\text{max}$, and signal at the ground level as a function of zenith angle. The method relies on the assumption that the mass composition is the same at all zenith angles, while the atmospheric shower development and attenuation depend on composition in a correlated way. In the present work, for each of the three leading LHC-tuned hadronic interaction models, we allow a global shift $\Delta X_\text{max}$ of the predicted shower maximum, which is the same for every mass and energy, and a rescaling $R_\text{Had}$ of the hadronic component at the ground level which is constant with the zenith angle. We apply the analysis to 2297 events reconstructed with both the fluorescence and surface detectors of the Pierre Auger Observatory with energies $10^{18.5-19.0}$ eV and zenith angles below 60$^\circ$. Given the modeling assumptions made in this analysis, the best fit reaches its optimum value when shifting the $X_\text{max}$ predictions of hadronic interaction models to deeper values and increasing the hadronic signal. This change in the predicted $X_\text{max}$ scale alleviates the previously identified model deficit in the hadronic signal (commonly called the muon puzzle) but does not fully remove it. Because of the size of the adjustments $\Delta X_\text{max}$ and $R_\text{Had}$ and the large number of events in the sample, the statistical significance of need for these adjustments is large, greater than 5$\sigma_\text{stat}$, even for the combination of the systematic experimental shifts within 1$\sigma_\text{sys}$ that are the most favorable for the models.

Although significant progress has been achieved in recent surveys of the magnetism in massive stars, the origin of the detected magnetic fields remains to be the least understood topic in their studies. We present an analysis of 61 high-resolution spectropolarimetric observations of 36 systems with O-type primaries, among them ten known particle-accelerating colliding-wind binaries exhibiting synchrotron radio emission. Our sample consists of multiple systems with components at different evolutionary stages with wide and tight orbits and different types of interactions. For the treatment of the complex composite spectra of the multiple systems, we used a special procedure involving different line masks populated for each element separately. Out of the 36 systems, 22 exhibit in their LSD Stokes V profiles definitely detected Zeeman features, among them seven systems with colliding winds. For fourteen systems the detected Zeeman features are most likely associated with O-type components whereas for three systems we suggest an association with an early B-type component. For the remaining five systems the source of the field is unclear. Marginal evidence for the detection of a Zeeman feature is reported for eleven systems and non-detection for three systems. The large number of systems with definitely detected Zeeman features presents a mystery, but probably indicates that multiplicity plays a definite role in the generation of magnetic fields in massive stars. The newly found magnetic systems are supreme candidates for spectropolarimetric monitoring over their orbital and rotation periods to obtain trustworthy statistics on the magnetic field geometry and the distribution of field strength.

This paper presents a study of the use of numerical simulation and Bayesian optimisation techniques to investigate the dynamics of celestial systems. Initially, the study focuses on Lagrange points in restricted three-body systems where a 2D three-body system simulator is employed to locate the five Lagrange points. An appropriate loss function is developed to capture the gravitational stability of the system, and the stability properties of the different Lagrange points are explored. Additionally, the study investigates how varying the number of variables for the satellite impacts the search for the Lagrange points. Finally, the scope of the study is expanded to explore stable configurations in multi-star systems represented by regular convex n-gons. In this case, Bayesian optimisation is used to find suitable settings for the n-gon's radius and the stars' velocity vectors, such that the overall system is stable.

The TRAPPIST-1 system is remarkable for its seven planets that are similar in size, mass, density, and stellar heating to the rocky planets Venus, Earth, and Mars in our own Solar System (Gillon et al. 2017). All TRAPPIST-1 planets have been observed with the transmission spectroscopy technique using the Hubble or Spitzer Space Telescopes, but no atmospheric features have been detected or strongly constrained (Ducrot et al. 2018; de Wit et al. 2018; Zhang et al. 2018; Garcia et al. 2022). TRAPPIST-1 b is the closest planet to the system's M dwarf star, and it receives 4 times as much irradiation as Earth receives from the Sun. This relatively large amount of stellar heating suggests that its thermal emission may be measurable. Here we present photometric secondary eclipse observations of the Earth-sized TRAPPIST-1 b exoplanet using the F1500W filter of the MIRI instrument on JWST. We detect the secondary eclipse in each of five separate observations with 8.7-sigma confidence when all data are combined. These measurements are most consistent with the re-radiation of the TRAPPIST-1 star's incident flux from only the dayside hemisphere of the planet. The most straightforward interpretation is that there is little or no planetary atmosphere redistributing radiation from the host star and also no detectable atmospheric absorption from carbon dioxide (CO$_2$) or other species.

Gravitational polarization is examined for equilibrium self-gravitating polytropic sheets perturbed by gravitational field due to test mass sheet. We find equilibrium solutions to the corresponding perturbed Lane-Emden equations for non-negative polytropic indexes. It is shown that gravitational polarization may be observed even in a finite extent of self-gravitating systems in addition to previously discussed infinite systems. In the polytropic sheets, the maximum gravitational amplification gets greater with a higher polytropic index while the height at which the maximum amplification occurs gets lower. The ratio of height change to the original height increases with polytropic index. The last result constrains the linear approximation method used for the present perturbation method.

We present an extensive catalog of 5405 early-type dwarf (dE) galaxies located in the various environments, i.e., clusters, groups and fields, of the local universe ($z$ $<$ 0.01). The dEs are selected through visual inspection of the Legacy survey's $g$-$r$-$z$ combined tri-color images. The inspected area, covering a total sky area of 7643 deg$^{2}$, encompasses two local clusters, Virgo and Fornax, 265 groups, and the regions around 586 field galaxies of $M_{K}$ $<$ $-$21 mag. \trev{The catalog aims to be one of the most extensive and publicly accessible collections of data on dE, despite its complex completeness limits that may not accurately represent its statistical completeness.} The strength of the catalog lies in the morphological characteristics, including nucleated, tidal, and ultradiffuse dE. The two clusters contribute nearly half (2437 out of 5405) dEs, and the 265 groups contribute 2103 dEs. There are 864 dEs in 586 fields, i.e., $\sim$\,1.47 dEs per field. Using a standard definition commonly used in literature, we identify 100 ultra-diffuse galaxies (UDGs), which take $\sim$\,2\,\% of the dE population. We find that 40\,\% of our sample dEs harbor a central nucleus, and among the UDG population, a majority, 79\%, are nonnucleated. About 1.3\,\% of dEs suffer from ongoing tidal disturbance by nearby massive galaxies, and only 0.03\,\% show the sign of recent dwarf-dwarf mergers. The association between dEs and their nearest bright neighbor galaxies suggests that dEs are more likely created where their neighbors are non-star-forming ones.

We present recent 2-port vector network analyzer (VNA) measurements of the complete set of scattering parameters for the antenna used within the Long Wavelength Array (LWA) and the associated front end electronics (FEEs). Full scattering parameter measurements of the antenna yield not only the reflection coefficient for each polarization, S11 and S22, but also the coupling between polarizations, S12 and S21. These had been previously modeled using simulations, but direct measurements had not been obtained until now. The measurements are used to derive a frequency dependent impedance mismatch factor (IMF) which represents the fraction of power that is passed through the antenna-FEE interface and not reflected due to a mismatch between the impedance of the antenna and the impedance of the FEE. We also present results from a two antenna experiment where each antenna is hooked up to a separate port on the VNA. This allows for cross-antenna coupling to be measured for all four possible polarization combinations. Finally, we apply the newly measured IMF and FEE forward gain corrections to LWA data to investigate how well they remove instrumental effects.

Recently, we witnessed how the synergy of small satellite technology and solar sailing propulsion enables new missions. Together, small satellites with lightweight instruments and solar sails offer affordable access to deep regions of the solar system, also making it possible to realize hard-to-reach trajectories that are not constrained to the ecliptic plane. Combining these two technologies can drastically reduce travel times within the solar system, while delivering robust science. With solar sailing propulsion capable of reaching the velocities of ~5-10 AU/yr, missions using a rideshare launch may reach the Jovian system in two years, Saturn in three. The same technologies could allow reaching solar polar orbits in less than two years. Fast, cost-effective, and maneuverable sailcraft that may travel outside the ecliptic plane open new opportunities for affordable solar system exploration, with great promise for heliophysics, planetary science, and astrophysics. Such missions could be modularized to reach different destinations with different sets of instruments. Benefiting from this progress, we present the "Sundiver" concept, offering novel possibilities for the science community. We discuss some of the key technologies, the current design of the Sundiver sailcraft vehicle and innovative instruments, along with unique science opportunities that these technologies enable, especially as this exploration paradigm evolves. We formulate policy recommendations to allow national space agencies, industry, and other stakeholders to establish a strong scientific, programmatic, and commercial focus, enrich and deepen the space enterprise and broaden its advocacy base by including the Sundiver paradigm as a part of broader space exploration efforts.

HESS J1809-193 is an extended TeV $\gamma$-ray source and the origin of its $\gamma$-ray emission remains ambiguous. Pulsar wind nebula (PWN) of PSR J1809-1917 laying inside the extended $\gamma$-ray emission is a possible candidate. Powered by the central pulsar, ultrarelativistic electrons in PWN can produce radio to X-ray emission through synchrotron and $\gamma$-ray emission by inverse Compton (IC) scattering. To check whether this PWN is the counterpart of HESS J1809-193, we analyzed Chandra X-ray radial intensity profile and the spectral index profile of this PWN. We then adopt a one-zone isotropic diffusion model to fit the keV and the TeV data. We find diffuse nonthermal X-ray emission extending beyond PWN, which is likely an X-ray halo radiated by escaping electron/positron pairs from the PWN. A relatively strong magnetic field of $21\,\mu$G is required to explain the spatial evolution of the X-ray spectrum (i.e., the significant softening of the spectrum with increasing distance from the pulsar), which, however, would suppress the IC radiation of pairs. Our result implies that a hadronic component may be needed to explain HESS J1809-193.

Asteroid 4 Vesta has a set of parallel troughs aligned with its equator. Although previous evaluations suggest that it is of shock fracturing tectonic origin, we propose that the equatorial troughs can be created by secondary cratering from the largest impact basin, Rheasilvia. We calculated the trajectories of ejecta particles from Rheasilvia by considering Vesta's rapid rotation. As a result, we found that secondary craters should be parallel to the latitude. In particular, if we assume that ejecta particles are launched at an initial launch velocity of approximately 350-380 m/s and a launch angle of 25 degree, the parallel equatorial troughs, the Divalia Fossae, can be suitably explained by secondary cratering. This model works well on objects, such as Haumea, Salacia, and Chariklo, but not on Mercury, the Moon, and regular satellites.

We present a new 3D MHD heliospheric model for space-weather forecasting driven by boundary conditions defined from white-light observations of the solar corona. The model is based on the MHD code PLUTO, constrained by an empirical derivation of the solar wind background properties at 0.1au. This empirical method uses white-light observations to estimate the position of the heliospheric current sheet. The boundary conditions necessary to run HelioCast are then defined from pre-defined relations between the necessary MHD properties (speed, density and temperature) and the distance to the current sheet. We assess the accuracy of the model over six Carrington rotations during the first semester of 2018. Using point-by-point metrics and event based analysis, we evaluate the performances of our model varying the angular width of the slow solar wind layer surrounding the heliospheric current sheet. We also compare our empirical technique with two well tested models of the corona: Multi-VP and WindPredict-AW. We find that our method is well suited to reproduce high speed streams, and does -- for well chosen parameters -- better than full MHD models. The model shows, nonetheless, limitations that could worsen for rising and maximum solar activity.

I describe the selection and initial characterisation of 20 eclipsing binary stars that are suitable for calibration and testing of stellar models and data analysis algorithms used by the PLATO mission and spectroscopic surveys. The binary stars selected are F-/G-type dwarf stars with M-type dwarf companions that contribute less than 2% of the flux at optical wavelengths. The light curves typically show well-defined total eclipses with very little variability between the eclipses. I have used near-infrared spectra obtained by the APOGEE survey to measure the spectroscopic orbit for both stars in HD22064. Combined with an analysis of the TESS light curve, I derive the following masses and radii: $M_1 = 1.35 \pm 0.03 M_{\odot}$, $M_2 = 0.58 \pm 0.01 M_{\odot}$, $R_1 = 1.554 \pm 0.014 R_{\odot}$, $R_2 = 0.595 \pm 0.008 R_{\odot}$. Using $R_1$ and the parallax from Gaia EDR3, I find that the primary star's angular diameter is $\theta = 0.1035 \pm 0.0009 $ mas. The apparent bolometric flux of the primary star is ${\mathcal F}_{\oplus,0} = (7.51\pm 0.09)\times10^{-9}$ erg cm$^{-2}$ s$^{-1}$. Hence, this F2V star has an effective temperature $T_{\rm eff,1} = 6763{\rm\,K} \pm 39{\rm \,K}$. HD22064 is an ideal benchmark star that can be used for ``end-to-end'' tests of the stellar parameters measured by large-scale spectroscopic surveys, or stellar parameters derived from asteroseismology with PLATO. The techniques described here for HD22064 can be applied to the other eclipsing binaries in the sample in order to create an all-sky network of such benchmark stars.

Astronomers have been measuring the separations and position angles between the two components of binary stars since William Herschel began his observations in 1781. In 1970, Anton Labeyrie pioneered a method, speckle interferometry, that overcomes the usual resolution limits induced by atmospheric turbulence by taking hundreds or thousands of short exposures and reducing them in Fourier space. Our 2022 automation of speckle interferometry allowed us to use a fully robotic 1.0-meter PlaneWave Instruments telescope, located at the El Sauce Observatory in the Atacama Desert of Chile, to obtain observations of many known binaries with established orbits. The long-term objective of these observations is to establish the precision, accuracy, and limitations of this telescope's automated speckle interferometry measurements. This paper provides an early overview of the Known Binaries Project and provide example results on a small-separation (0.27") binary, WDS 12274-2843 B 228.

Determining accurate effective temperatures of stars buried in the dust-obscured Galactic regions is extremely difficult from photometry. Fortunately, high-resolution infrared spectroscopy is a powerful tool for determining the temperatures of stars with no dependence on interstellar extinction. It has long been known that the depth ratios of temperature-sensitive and relatively insensitive spectral lines are excellent temperature indices. In this work, we provide the first extensive line depth ratio (LDR) method application in the infrared region that encompasses both \H\ and \K\ bands (1.48 $\mu$m - 2.48 $\mu$m). We applied the LDR method to high-resolution (R $\simeq$ 45,000) \H\ and \K-band spectra of 110 stars obtained with the Immersion Grating Infrared Spectrograph (IGRINS). Our sample contained stars with 3200 $<$ \teff\ (K) $<$ 5500, 0.20 $\leq$ log g $<$ 4.6, and $-$1.5 $<$ [M/H] $<$ 0.5. Application of this method in the \K-band yielded 21 new LDR$-$\teff\ relations. We also report five new LDR$-$\teff\ relations found in the \H-band region, augmenting the relations already published by other groups. The temperatures found from our calibrations provide reliable temperatures within $\sim$70 K accuracy compared to spectral \teff\ values from the literature.

MOdified Gravity (MoG)) is widely constrained in different astrophysical and astronomical systems. Since these different systems are based on different scales it is not trivial to get a combined constraint that is based on different phenomenology. Here, for the first time (to the best of our knowledge), we combine constraints for MoG from late time Cosmology and the orbital motion of the stars around the galactic center. MoG give different potentials that are tested directly in the galactic center. The cosmological data set includes the type Ia supernova and baryon acoustic oscillations. For the galactic star center data set we use the published orbital measurements of the S2 star. The constraints on the universal parameter $\beta$ from the combined system give: $\beta_{HS}=0.154 \pm 0.109$ for the Hu-Sawicki model, while $\beta_{St}= 0.309 \pm 0.19 $ for the Starobinsky model. These results improve on the cosmological results we obtain. The results show that {{\it combined constraint}} from different systems yields a stronger constraint for different theories under consideration. Future measurements from the galactic center and from cosmology will give better constraints on MoG.

Measurements of the characteristic length scale $r_s$ of the baryon acoustic oscillations (BAO) provide a robust determination of the distance-redshift relation. Currently, the best (sub-per cent) estimate of $r_s$ at the drag epoch is provided by Cosmic Microwave Background (CMB) observations assuming the validity of the standard $\Lambda$CDM model at $z \sim 1000$. Therefore, inferring $r_s$ from low-$z$ observations in a model-independent way and comparing its value with CMB estimates provides a consistency test of the standard cosmology and its assumptions at high-$z$. In this paper, we address this question and estimate the absolute BAO scale combining angular BAO measurements and type Ia Supernovae data. Our analysis uses two different methods to connect these data sets and finds a good agreement between the low-$z$ estimates of $r_{s}$ with the CMB sound horizon at drag epoch, regardless of the value of the Hubble constant $H_0$ considered. These results highlight the robustness of the standard cosmology at the same time that they also reinforce the need for more precise cosmological observations at low-$z$.

We investigate the spin state of a protoplanet during the pebble accretion influenced by the gas flow in the gravitational potential of the protoplanet and how it depends on the planetary mass, the headwind speed, the distance from the host star, and the pebble size. We perform nonisothermal three-dimensional hydrodynamical simulations in a local frame to obtain the gas flow around the planet. We then numerically integrate three-dimensional orbits of pebbles under the obtained gas flow. Finally, assuming uniform spatial distribution of incoming pebbles, we calculate net spin by summing up specific angular momentum that individual pebbles transfer to the protoplanet at impacts. We find that a protoplanet with the envelope acquires prograde net spin rotation regardless of the planetary mass, the pebble size, and the headwind speed of the gas. This is because accreting pebbles are dragged by the envelope that commonly has prograde rotation. As the planetary mass or orbital radius increases, the envelope is thicker and the prograde rotation is faster, resulting in faster net prograde spin. When the dimensionless thermal mass of the planet, $m = R_{\mathrm{Bondi}} / H$, where $R_{\mathrm{Bondi}}$ and $H$ are the Bondi radius and the disk gas scale height, is larger than a certain critical mass ($m \gtrsim 0.3$ at $0.1 \, \mathrm{au}$ or $m \gtrsim 0.1$ at $1 \, \mathrm{au}$), the spin rotation exceeds the breakup one. The predicted spin frequency reaches the breakup one at the planetary mass $m_{\mathrm{iso,rot}} \sim 0.1 \, (a / 1 \, \mathrm{au})^{-1/2}$ (where $a$ is the orbital radius), suggesting that the protoplanet cannot grow beyond $m_{\mathrm{iso,rot}}$. It is consistent with the Earth's current mass and could help the formation of the Moon by a giant impact on fast-spinning proto-Earth.

Turbulent processes at work in the intracluster medium perturb this environment, displacing gas, and creating local density fluctuations that can be quantified via X-ray surface brightness fluctuation analyses. Improved knowledge of these phenomena would allow for a better determination of the mass of galaxy clusters, as well as a better understanding of their dynamic assembly. In this work, we aim to set constraints on the structure of turbulence using X-ray surface brightness fluctuations. We seek to consider the stochastic nature of this observable and to constrain the structure of the underlying power spectrum. We propose a new Bayesian approach, relying on simulation-based inference to account for the whole error budget. We used the X-COP cluster sample to individually constrain the power spectrum in four regions and within $R_{500}$. We spread the analysis on the 12 systems to alleviate the sample variance. We then interpreted the density fluctuations as the result of either gas clumping or turbulence. For each cluster considered individually, the normalisation of density fluctuations correlates positively with the Zernike moment and centroid shift, but negatively with the concentration and the Gini coefficient. The spectral index within $R_{500}$ and evaluated over all clusters is consistent with a Kolmogorov cascade. The normalisation of density fluctuations, when interpreted in terms of clumping, is consistent within $0.5 R_{500}$ with the literature results and numerical simulations; however, it is higher between 0.5 and $1 R_{500}$. Conversely, when interpreted on the basis of turbulence, we deduce a non-thermal pressure profile that is lower than the predictions of the simulations within 0.5 $R_{500}$, but still in agreement in the outer regions. We explain these results by the presence of central structural residues that are remnants of the dynamic assembly of the clusters.

We present a new technique to generate the light curves of RRab stars in different photometric bands ($I$ and $V$ bands) using Artificial Neural Networks (ANN). A pre-computed grid of models was used to train the ANN, and the architecture was tuned using the $I$ band light curves. The best-performing network was adopted to make the final interpolators in the $I$ and $V$ bands. The trained interpolators were used to predict the light curve of RRab stars in the Magellanic Clouds, and the distances to the LMC and SMC were determined based on the reddening independent Wesenheit index. The estimated distances are in good agreement with the literature. The comparison of the predicted and observed amplitudes, and Fourier amplitude ratios showed good agreement, but the Fourier phase parameters displayed a few discrepancies. To showcase the utility of the interpolators, the light curve of the RRab star EZ Cnc was generated and compared with the observed light curve from the Kepler mission. The reported distance to EZ Cnc was found to be in excellent agreement with the updated parallax measurement from Gaia EDR3. Our ANN interpolator provides a fast and efficient technique to generate a smooth grid of model light curves for a wide range of physical parameters, which is computationally expensive and time-consuming using stellar pulsation codes.

The radial density of planets increases with depth due to compressibility, leading to impacts on their convective dynamics. To account for these effects, including the presence of a quasi-adiabatic temperature profile and entropy sources due to dissipation, the compressibility is expressed through a dissipation number, $\mathcal{D}$, proportional to the planet's radius and gravity. In Earth's mantle, compressibility effects are moderate, but in large rocky or liquid exoplanets (Super-Earths), the dissipation number can become very large. This paper explores the properties of compressible convection when the dissipation number is significant. We start by selecting a simple Murnaghan equation of state that embodies the fundamental properties of condensed matter at planetary conditions. Next, we analyze the characteristics of adiabatic profiles and demonstrate that the ratio between the bottom and top adiabatic temperatures is relatively small and probably less than 2. We examine the marginal stability of compressible mantles and reveal that they can undergo convection with either positive or negative superadiabatic Rayleigh numbers. Lastly, we delve into simulations of convection performed using the exact equations of mechanics, neglecting inertia (infinite Prandtl number case), and examine their consequences for Super-Earths dynamics.

[abridged] We analyzed mid-infrared high-contrast coronagraphic images of the beta Pictoris system, taking advantage of the NEAR experiment using the VLT/VISIR instrument. The goal of our analysis is to investigate both the detection of the planet beta Pictoris b and of the disk features at mid-IR wavelengths. In addition, by combining several epochs of observation, we expect to constrain the position of the known clumps and improve our knowledge on the dynamics of the disk. To evaluate the planet b flux contribution, we extracted the photometry and compared it to the flux published in the literature. In addition, we used previous data from T-ReCS and VISIR, to study the evolution of the position of the southwest clump that was initially observed in the planetary disk back in 2003. While we did not detect the planet b, we were able to put constraints on the presence of circumplanetary material, ruling out the equivalent of a Saturn-like planetary ring around the planet. The disk presents several noticeable structures, including the known southwest clump. Using a 16-year baseline, sampled with five epochs of observations, we were able to examine the evolution of the clump: the clump orbits in a Keplerian motion with an sma of 56.1+-0.4 au. In addition to the known clump, the images clearly show the presence of a second clump on the northeast side of the disk and fainter and closer structures that are yet to be confirmed. We found correlations between the CO clumps detected with ALMA and the mid-IR images. If the circumplanetary material were located at the Roche radius, the maximum amount of dust determined from the flux upper limit around beta Pictoris b would correspond to the mass of an asteroid of 5 km in diameter. Finally, the Keplerian motion of the southwestern clump is possibly indicative of a yet-to-be-detected planet or signals the presence of a vortex.

We present a sample of 19,583 ultracool dwarf candidates brighter than z $\leq 23$ selected from the Dark Energy Survey DR2 coadd data matched to VHS DR6, VIKING DR5 and AllWISE covering $\sim$ 4,800 $deg^2$. The ultracool candidates were first pre-selected based on their (i-z), (z-Y), and (Y-J) colours. They were further classified using a method that compares their optical, near-infrared and mid-infrared colours against templates of M, L and T dwarfs. 14,099 objects are presented as new L and T candidates and the remaining objects are from the literature, including 5,342 candidates from our previous work. Using this new and deeper sample of ultracool dwarf candidates we also present: 20 new candidate members to nearby young moving groups (YMG) and associations, variable candidate sources and four new wide binary systems composed of two ultracool dwarfs. Finally, we also show the spectra of twelve new ultracool dwarfs discovered by our group and presented here for the first time. These spectroscopically confirmed objects are a sanity check of our selection of ultracool dwarfs and photometric classification method.

We demonstrate that planet formation via pebble accretion is sensitive to external photoevaporation of the outer disc. In pebble accretion, planets grow by accreting from a flux of solids (pebbles) that radially drift inwards from the pebble production front. If external photoevaporation truncates the outer disc fast enough, it can shorten the time before the pebble production front reaches the disc outer edge, cutting off the supply of pebble flux for accretion, hence limiting the pebble mass reservoir for planet growth. Conversely, cloud shielding can protect the disc from strong external photoevaporation and preserve the pebble reservoir. Because grain growth and drift can occur quickly, shielding even on a short time-scale (<1 Myr) can have a non-linear impact on the properties of planets growing by pebble accretion. For example a $10^{-3} M_\oplus$ planetary seed at 25 au stays at 25 au with a lunar mass if the disc is immediately irradiated by a $10^3$ G$_0$ field, but grows and migrates to be approximately Earth-like in both mass and orbital radius if the disc is shielded for just 1 Myr. In NGC 2024, external photoevaporation is thought to happen to discs that are <0.5 Myr old, which coupled with the results here suggests that the exact planetary parameters can be very sensitive to the star forming environment. Universal shielding for time-scales of at least $\sim1.5$ Myr would be required to completely nullify the environmental impact on planetary architectures.

Catastrophically evaporating planets are observed through their dusty tails, formed through rocky material evaporated from their highly irradiated molten surfaces. The composition of these tails offers an avenue for studying the composition of rocky exoplanets, but only if the evolution of the underlying interior is understood. This is because it is the interior evolution that sets the composition at the base of the mass outflow. In this work, we present a model of the evolution of the interiors of catastrophically evaporating planets. Its basis is a one-dimensional code that takes into account energy flow through conduction and convection as well as melting. We find that these planets are likely to be entirely solid when significant mass loss occurs, other than a thin magma pool on the day side. Consequently, the outflows from the planets, and thus the dust tails, sample material only from the surface of the planet. We also use our model to investigate the occurrence rate of planets that can catastrophically evaporate, and find that, on average, a star in the Kepler sample has approximately one such planet. Our value is above, but within an order of magnitude of, the occurrence rate inferred from the Kepler statistics for Super-Earths, implying that exotic mechanisms to produce the catastrophically evaporating planet population may not be required. We also find that the range of substellar temperatures of the observed systems are well explained by recent theoretical models which only produce dust in a limited temperature region.

We review the status of $f(R,T)$ theories, where $T$ is the trace of the energy momentum tensor $T^{\mu\nu}$, concerning the evolution of the cosmological flat Friedmann-Lema\^itre-Robertson-Walker (FLRW) background expansion. We start focusing on the modified Friedmann equations for the case of a minimally coupled gravitational Lagrangian of the type $f(R,T)=R +\alpha e^{\beta T} + \gamma_{n} T^{n}$. With this choice one is allowed to cover all existing proposals in the literature via four free parameters and all relevant $f(R,T)$ models as well as the $\Lambda$CDM model can be achieved in the appropriate limit. We show that in such minimally coupled case there exists a useful constraining relation between the effective fractionary total matter density with arbitrary equation of state parameter and the modified gravity parameters. Then, with this association the modified gravity sector can be independently constrained using estimations of the gas mass fraction in galaxy clusters. Using cosmological background data and demanding the universe is old enough to accommodate the existence of Galactic globular clusters with estimated age of at least $\sim 13$ Gyrs we find a narrow range of the modified gravity free parameter space in which this class of theories remains cosmologically viable. As expected, this preferred parameter space region accommodates the $\Lambda$CDM limit of $f(R,T)$ models. We also work out the non-minimally coupled case in the metric-affine formalism and find that there are no viable cosmologies in the latter situation.

This paper continues our study of radio pulsar emission-beam configurations with the primary intent of extending study to the lowest possible frequencies. Here we focus on a group of 133 more recently discovered pulsars, most of which were included in the (100-200 MHz) LOFAR High Band Survey, observed with Arecibo at 1.4 GHz and 327 MHz, and some observed at decameter wavelengths. Our analysis framework is the core/double-cone beam model, and we took opportunity to apply it as widely as possible, both conceptually and quantitatively, while highlighting situations where modeling is difficult, or where its premises may be violated. In the great majority of pulsars, beam forms consistent with the core/double-cone model were identified. Moreover, we found that each pulsar's beam structure remained largely constant over the frequency range available; where profile variations were observed, they were attributable to different component spectra and in some instances to varying conal beam sizes. As an Arecibo population, many or most of the objects tend to fall in the Galactic anticenter region and/or at high Galactic latitudes, so overall it includes a number of nearer, older pulsars. We found a number of interesting or unusual characteristics in some of the pulsars that would benefit from additional study. The scattering levels encountered for this group are low to moderate, apart from a few pulsars lying in directions more toward the inner Galaxy.

For more than two years, we monitored with the HARPS-N spectrograph the 400 Myr-old star HD\,63433, which hosts two close-in (orbital periods $P_b\sim7.1$ and $P_c\sim20.5$ days) sub-Neptunes detected by the TESS space telescope, and it was announced in 2020. Using radial velocities and additional TESS photometry, we aim to provide the first measurement of their masses, improve the measure of their size and orbital parameters, and study the evolution of the atmospheric mass-loss rate due to photoevaporation. We tested state-of-the-art analysis techniques and different models to mitigate the dominant signals due to stellar activity that are detected in the radial velocity time series. We used a hydro-based analytical description of the atmospheric mass-loss rate, coupled with a core-envelope model and stellar evolutionary tracks, to study the past and future evolution of the planetary masses and radii. We derived new measurements of the planetary orbital periods and radii ($P_b=7.10794\pm0.000009$ d, $r_b=2.02^{+0.06}_{-0.05}$ $R_{\oplus}$; $P_c=20.54379\pm0.00002$ d, $r_c=2.44\pm0.07$ $R_{\oplus}$), and determined mass upper limits ($m_b\lesssim$11 $M_{\oplus}$; $m_c\lesssim$31 $M_{\oplus}$; 95$\%$ confidence level), with evidence at a 2.1--2.7$\sigma$ significance level that HD\,63433\,c might be a dense mini-Neptune with a Neptune-like mass. For a grid of test masses below our derived dynamical upper limits, we found that HD\,63433\,b has very likely lost any gaseous H-He envelope, supporting HST-based observations that are indicative of there being no ongoing atmospheric evaporation. HD\,63433\,c will keep evaporating over the next $\sim$5 Gyr if its current mass is $m_c\lesssim$15 $M_{\oplus}$, while it should be hydrodynamically stable for higher masses.

In the context of atmospheric shower arrays designed for $\gamma$-ray astronomy and in the context of the ALTO project, we present: a study of the impact of heavier nuclei in the cosmic-ray background on the estimated $\gamma$-ray detection performance on the basis of dedicated Monte Carlo simulations, a method to calculate the sensitivity to a point-like source, and finally the required observation times to reach a firm detection on a list of known point-like sources.

Several observational studies have shown that many Galactic globular clusters (GCs) are characterised by internal rotation. Theoretical studies of the dynamical evolution of rotating clusters have predicted that, during their long-term evolution, these stellar systems should develop a dependence of the rotational velocity around the cluster's centre on the mass of stars, with the internal rotation increasing for more massive stars. In this paper we present the first observational evidence of the predicted rotation-mass trend. In our investigation, we exploited the $\mathit{Gaia}$ Data Release 3 catalogue of three GCs: NGC 104 (47 Tuc), NGC 5139 ($\omega$ Cen) and NGC 5904 (M 5). We found clear evidence of a cluster rotation-mass relation in 47 Tuc and M 5, while in $\omega$ Cen, the dynamically youngest system among the three clusters studied here, no such trend was detected.

Much of the research in supernova cosmology is based on an assumption that the peak luminosity of type Ia supernovae (SNe Ia), after a standardization process, is independent of the galactic environment. A series of recent studies suggested that there is a significant correlation between the standardized luminosity and the progenitor age of SNe Ia. The correlation found in the most recent work by Lee et al. is strong enough to explain the extra dimming of distant SNe Ia, and therefore casts doubts on the direct evidence of cosmic acceleration. The present work improves the previous work by incorporating the uncertainties of progenitor ages, which were ignored in Lee et al., into a fully Bayesian inference framework. We find a weaker dependence of supernova standardized luminosity on the progenitor age, but the detection of correlation remains significant (3.3$\sigma$). Assuming that such correlation can be extended to high redshift and applying it to the Pantheon SN Ia data set, we find that the correlation cannot be the primary cause of the observed extra dimming of distant SNe Ia. Further more, we use the PAge formalism, which is a good approximation to many dark energy and modified gravity models, to do a model comparison. The concordance Lambda cold dark matter model remains a good fit when the progenitor-age dependence of SN Ia luminosity is corrected. The best-fit parameters, however, are in $\sim 2\sigma$ tension with the standard values inferred from cosmic microwave background observations.

This survey is devoted to recent developments in the statistical analysis of spherical data, with a view to applications in Cosmology. We will start from a brief discussion of Cosmological questions and motivations, arguing that most Cosmological observables are spherical random fields. Then, we will introduce some mathematical background on spherical random fields, including spectral representations and the construction of needlet and wavelet frames. We will then focus on some specific issues, including tools and algorithms for map reconstruction (\textit{i.e.}, separating the different physical components which contribute to the observed field), geometric tools for testing the assumptions of Gaussianity and isotropy, and multiple testing methods to detect contamination in the field due to point sources. Although these tools are introduced in the Cosmological context, they can be applied to other situations dealing with spherical data. Finally, we will discuss more recent and challenging issues such as the analysis of polarization data, which can be viewed as realizations of random fields taking values in spin fiber bundles.

Several energetic disturbances have been identified as triggers of the large-amplitude oscillations (LAOs) in prominences. However, the mechanisms for LAOs excitation are not well understood. We aim to study these mechanisms, performing time-dependent numerical simulations in 2.5D and 2D setups using magnetohydrodynamic (MHD) code MANCHA3D. Two types of disturbances are applied to excite prominence oscillations, such as a perturbation associated with an eruption and the waves caused by an artificial energy release. In the simulation with the eruption, we obtain that it does not produce LAOs in the prominence located in its vicinity. While the erupting flux rope rises, an elongated current sheet forms behind it, which becomes unstable and breaks into plasmoids. The downward-moving plasmoids cause perturbations in the velocity field by merging with the post-reconnection loops. This velocity perturbation propagates in the surroundings and perturbs the nearby prominence. The analysis of the oscillatory motions of the prominence plasma reveals the excitation of small-amplitude oscillations (SAOs), which are a mixture of longitudinal and vertical oscillations. In the simulation with a distant artificial perturbation, a fast-mode shock wave is produced, and it gradually reaches two flux rope prominences at different distances. This shock wave excites vertical LAOs and longitudinal SAOs with similar amplitudes, periods, and damping times in both prominences. Finally, in the experiment with the external triggering of LAOs in a dipped arcade prominence model, we find that, although the vector normal to the front of a fast-mode shock wave is parallel to the spine of the dipped arcade well before the contact, this wave does not excite longitudinal LAOs. When the wave front approaches the prominence, it pushes the dense plasma down, establishing vertical LAOs.

We review the current literature on the formation of Coherent Structures (CoSs) in strongly turbulent 3D magnetized plasmas. CoSs (Current Sheets (CS), magnetic filaments, large amplitude magnetic disturbances, vortices, and shocklets) appear intermittently inside a turbulent plasma and are collectively the locus of magnetic energy transfer (dissipation) into particle kinetic energy, leading to heating and/or acceleration of the latter. CoSs and especially CSs are also evolving and fragmenting, becoming locally the source of new clusters of CoSs. Strong turbulence can be generated by the nonlinear coupling of large amplitude unstable plasma modes, by the explosive reorganization of large scale magnetic fields, or by the fragmentation of CoSs. A small fraction of CSs inside a strongly turbulent plasma will end up reconnecting. Magnetic Reconnection (MR) is one of the potential forms of energy dissipation of a turbulent plasma. Analysing the evolution of CSs and MR in isolation from the surrounding CoSs and plasma flows may be convenient for 2D numerical studies, but it is far from a realistic modeling of 3D astrophysical, space and laboratory environments, where strong turbulence can be exited, as e.g. in the solar wind, the solar atmosphere, solar flares and Coronal Mass Ejections (CMEs), large scale space and astrophysical shocks, the magnetosheath, the magnetotail, astrophysical jets, Edge Localized Modes (ELMs) in confined laboratory plasmas (TOKAMAKS), etc.

Polarization of the cosmic microwave background (CMB) is sensitive to new physics violating parity symmetry, such as the presence of a pseudoscalar "axionlike" field. Such a field may be responsible for early dark energy (EDE), which is active prior to recombination and provides a solution to the so-called Hubble tension. The EDE field coupled to photons in a parity-violating manner would rotate the plane of linear polarization of the CMB and produce a cross-correlation power spectrum of $E$- and $B$-mode polarization fields with opposite parities. In this paper, we fit the $EB$ power spectrum predicted by the photon-axion coupling of the EDE model with a potential $V(\phi)\propto [1-\cos(\phi/f)]^3$ to polarization data from Planck. We find that the unique shape of the predicted $EB$ power spectrum is not favored by the data and obtain a first constraint on the photon-axion coupling constant, $g=(0.04\pm 0.16)M_{\text{Pl}}^{-1}$ (68% CL), for the EDE model that best fits the CMB and galaxy clustering data. This constraint is independent of the miscalibration of polarization angles of the instrument or the polarized Galactic foreground emission. Our limit on $g$ may have important implications for embedding EDE in fundamental physics, such as string theory.

We consider cosmic chronometer (CC) data for the Hubble parameter, quasar (QSO) luminosities data of X-rays and ultraviolet rays emission, and the latest measurements of the present value of the Hubble parameter from 2018 Planck mission (PL18), and SH0ES observations (SHOES) to constrain the present value of cosmic curvature density parameter. We consider three kinds of dark energy models: the $\Lambda$CDM model, the wCDM model, and the CPL parametrization. In all these three models, we find higher values of the matter-energy density parameter, $\Omega_{\rm m0}$ compared to the one obtained from the Planck 2018 mission of CMB observation. Also, we find evidence for a nonflat and closed Universe at 0.5$\sigma$ to 3$\sigma$ confidence levels. The flat Universe is almost 2 to 3$\sigma$, 1 to 1.5$\sigma$, and 0.5 to 1$\sigma$ away from the corresponding mean values, obtained in $\Lambda$CDM model, wCDM model, and CPL parametrization respectively obtained from different combinations of datasets. The evidence for nonzero cosmic curvature is lesser in dynamical dark energy models compared to the $\Lambda$CDM model. That means the evidence of nonzero cosmic curvature depends on the behavior of the equation of state of the dark energy. Since the values of the cosmic curvature are degenerate to the equation of state of the dark energy, we also consider a model independent analysis to constrain the cosmic curvature using the combination of Gaussian process regression analysis and artificial neural networks analysis. In the model independent analysis, we also find evidence for a closed Universe, and the flat Universe is almost 1$\sigma$ away. So, both the model dependent and independent analyses favor a closed Universe from the combinations of CC, QSO, and $H_0$ observations.

Although the number of exoplanets reported in the literature exceeds 5000 so far, only a few dozen of them are young planets ($\le$900 Myr). However, a complete characterization of these young planets is key to understanding the current properties of the entire population. Hence, it is necessary to constrain the planetary formation processes and the timescales of dynamical evolution by measuring the masses of exoplanets transiting young stars. We characterize and measure the masses of two transiting planets orbiting the 400 Myr old solar-type star HD\,63433, which is a member of the Ursa Major moving group. We analysed precise photometric light curves of five sectors of the TESS mission with a baseline of $\sim$750 days and obtained $\sim$150 precise radial velocity measurements with the visible and infrared arms of the CARMENES instrument at the Calar Alto 3.5 m telescope in two different campaigns of $\sim$500 days. We performed a combined photometric and spectroscopic analysis to retrieve the planetary properties of two young planets. The strong stellar activity signal was modelled by Gaussian regression processes. We have updated the transit parameters of HD\,63433\,b and c and obtained planet radii of R$_p^b$\,=\,2.140\,$\pm$\,0.087 R$_\oplus$ and R$_p^c$\,=\,2.692\,$\pm$\,0.108 R$_\oplus$. Our analysis allowed us to determine the dynamical mass of the outer planet with a 4$\sigma$ significance ($M_p^c$\,=\,15.54\,$\pm$\,3.86 M$_\oplus$) and set an upper limit on the mass of the inner planet at 3$\sigma$ ($M_p^b$\,$<$\,21.76 M$_\oplus$). According to theoretical models, both planets are expected to be sub-Neptunes, whose interiors mostly consist of silicates and water with no dominant composition of iron, and whose gas envelopes are lower than 2\% in the case of HD\,63433\,c. The envelope is unconstrained in HD\,63433\,b.

Transmission spectroscopy is still the preferred characterization technique for exoplanet atmospheres, although it presents unique challenges which translate into characterization bottlenecks when robust mitigation strategies are missing. Stellar contamination is one of such challenges that can overpower the planetary signal by up to an order of magnitude, and thus not accounting for stellar contamination can lead to significant biases in the derived atmospheric properties. Yet, accounting for stellar contamination may not be straightforward, as important discrepancies exist between state-of-the-art stellar models and measured spectra and between models themselves. Here we explore the extent to which stellar models can be used to reliably correct for stellar contamination and yield a planet's uncontaminated transmission spectrum. We find that (1) discrepancies between stellar models can dominate the noise budget of JWST transmission spectra of planets around stars with heterogeneous photospheres; (2) the true number of unique photospheric spectral components and their properties can only be accurately retrieved when the stellar models have a sufficient fidelity; and (3) under such optimistic circumstances the contribution of stellar contamination to the noise budget of a transmission spectrum is considerably below that of the photon noise for the standard transit observation setup. Therefore, we suggest (1) increased efforts towards development of model spectra of stars and their active regions in a data-driven manner; and (2) the development of empirical approaches for deriving spectra of photospheric components using the observatories with which the atmospheric explorations are carried out.

During the first 500 million years of cosmic history, the first stars and galaxies formed and seeded the cosmos with heavy elements. These early galaxies illuminated the transition from the cosmic "dark ages" to the reionization of the intergalactic medium. This transitional period has been largely inaccessible to direct observation until the recent commissioning of JWST, which has extended our observational reach into that epoch. Excitingly, the first JWST science observations uncovered a surprisingly high abundance of early star-forming galaxies. However, the distances (redshifts) of these galaxies were, by necessity, estimated from multi-band photometry. Photometric redshifts, while generally robust, can suffer from uncertainties and/or degeneracies. Spectroscopic measurements of the precise redshifts are required to validate these sources and to reliably quantify their space densities, stellar masses, and star formation rates, which provide powerful constraints on galaxy formation models and cosmology. Here we present the results of JWST follow-up spectroscopy of a small sample of galaxies suspected to be amongst the most distant yet observed. We confirm redshifts z > 10 for two galaxies, including one of the first bright JWST-discovered candidates with z = 11.4, and show that another galaxy with suggested z ~ 16 instead has z = 4.9, with strong emission lines that mimic the expected colors of more distant objects. These results reinforce the evidence for the rapid production of luminous galaxies in the very young Universe, while also highlighting the necessity of spectroscopic verification for remarkable candidates.

Strong gravitational lensing of supernovae is exceedingly rare. To date, only a handful of lensed supernovae are known. Despite their rarity, lensed supernovae have emerged as one of the most promising methods for measuring the current expansion rate of the Universe and breaking the Hubble tension. We present an extensive search for gravitationally lensed supernovae within the Zwicky Transient Facility (ZTF) public survey, covering 12,524 transients with good light curves discovered during four years of observations. We crossmatch a catalogue of known and candidate lens galaxies with our transient sample and find only one coincident source, which was due to chance alignment. To search for supernovae magnified by unknown lens galaxies, we test multiple methods that have been suggested in the literature, for the first time on real data. This includes selecting objects with extremely red colours and those that appear inconsistent with the host galaxy redshift. In both cases, we find a few hundred candidates, most of which are due to contamination from activate galactic nuclei, bogus detections, or unlensed supernovae. The false positive rate from these methods presents significant challenges for future surveys. In total, 65 unique transients were identified across all of our selection methods that required detailed manual rejection, which would be infeasible for larger samples. Overall, we do not find any compelling candidates for lensed supernovae, which is broadly consistent with previous estimates for the rate of lensed supernovae in the ZTF public survey and the number expected to pass the selection cuts we apply.

Current pulsar timing array (PTA) techniques for characterizing the spectrum of a nanohertz-frequency stochastic gravitational-wave background (SGWB) begin at the stage of timing data. This can be slow and memory intensive with computational scaling that will worsen PTA analysis times as more pulsars and observations are added. Given recent evidence for a common-spectrum process in PTA data sets and the need to understand present and future PTA capabilities to characterize the SGWB through large-scale simulations, we have developed efficient and rapid approaches that operate on intermediate SGWB analysis products. These methods refit SGWB spectral models to previously-computed Bayesian posterior estimations of the timing power spectra. We test our new methods on simulated PTA data sets and the NANOGrav $12.5$-year data set, where in the latter our refit posterior achieves a Hellinger distance from the current full production-level pipeline that is $\lesssim 0.1$. Our methods are $\sim10^2$--$10^4$ times faster than the production-level likelihood and scale sub-linearly as a PTA is expanded with new pulsars or observations. Our methods also demonstrate that SGWB spectral characterization in PTA data sets is driven by the longest-timed pulsars with the best-measured power spectral densities which is not necessarily the case for SGWB detection that is predicated on correlating many pulsars. Indeed, the common-process spectral properties found in the NANOGrav $12.5$-year data set are given by analyzing only the $\sim10$ longest-timed pulsars out of the full $45$ pulsar array, and we find that the ``shallowing'' of the common-process power-law model occurs when gravitational-wave frequencies higher than $\sim 50$~nanohertz are included. The implementation of our methods is openly available as a software suite to allow fast and flexible PTA SGWB spectral characterization and model selection.

Within the framework of modified gravity, the quasi-static and sub-horizon approximations are widely used in analyses aiming to identify departures from the concordance model at late-times. Under these approximations, it is generally assumed that time derivatives are subdominant with respect to spatial derivatives given that the relevant physical modes are those well inside the Hubble radius. Here, in the context of the effective fluid approach applied to $f(R)$ theories, we put forward a new parameterization which allows us to obtain analytical expressions for the gravitational potentials, whence for the effective dark energy perturbations. In order to track the validity of the two aforementioned approximations, we compare our results and the standard results found in the literature against full numerical solutions for two well-known toy-models; namely, the designer ($f$DES) model and the Hu-Sawicki (HS) model. We find that: $i)$ the sub-horizon approximation can be safely applied for scales $k \gtrsim 200 H_0$, $ii)$ in this ``safety region'', the quasi-static approximation is a very accurate description of the late-time dynamics even when dark energy significantly contribute to the cosmic budget, $iii)$ some relevant terms were neglected in the standard procedure yielding to inaccurate results in some of the dark energy effective fluid quantities; e.g. the sound speed. Instead, our expressions show overall agreement with respect to the full solutions. Therefore, our results indirectly indicate that the effective fluid expressions derived for more general modified gravity theories, such as Horndeski, should be revisited.

In many areas of the observational and experimental sciences data is scarce. Data observation in high-energy astrophysics is disrupted by celestial occlusions and limited telescope time while data derived from laboratory experiments in synthetic chemistry and materials science is time and cost-intensive to collect. On the other hand, knowledge about the data-generation mechanism is often available in the sciences, such as the measurement error of a piece of laboratory apparatus. Both characteristics, small data and knowledge of the underlying physics, make Gaussian processes (GPs) ideal candidates for fitting such datasets. GPs can make predictions with consideration of uncertainty, for example in the virtual screening of molecules and materials, and can also make inferences about incomplete data such as the latent emission signature from a black hole accretion disc. Furthermore, GPs are currently the workhorse model for Bayesian optimisation, a methodology foreseen to be a guide for laboratory experiments in scientific discovery campaigns. The first contribution of this thesis is to use GP modelling to reason about the latent emission signature from the Seyfert galaxy Markarian 335, and by extension, to reason about the applicability of various theoretical models of black hole accretion discs. The second contribution is to extend the GP framework to molecular and chemical reaction representations and to provide an open-source software library to enable the framework to be used by scientists. The third contribution is to leverage GPs to discover novel and performant photoswitch molecules. The fourth contribution is to introduce a Bayesian optimisation scheme capable of modelling aleatoric uncertainty to facilitate the identification of material compositions that possess intrinsic robustness to large scale fabrication processes.

We consider a model of a random media with fixed and finite memory time with abrupt losses of memory (renovation model). Within the memory intervals we can observe either amplification or oscillation of the vector field in a given particle. The cumulative effect of amplifications in many subsequent intervals leads to amplification of the mean field and mean energy. Similarly, the cumulative effect of intermittent amplifications or oscillations also leads to amplification of the mean field and mean energy, however, at a lower rate. Finally, the random oscillations alone can resonate and yield the growth of the mean field and energy. These are the three mechanisms that we investigate and compute analytically and numerically the growth rates based on the Jacobi equation with the random curvature parameter.

We analyse the dynamics of a light scalar field responsible for the $\mu$ term of the Higgs potential and coupled to matter via the Higgs-portal mechanism. We find that this dilaton model is stable under radiative corrections induced by the standard model particle masses. When the background value of the scalar field is stabilised at the minimum of the scalar potential, the scalar field fluctuations only couple quadratically to the massive fields of the standard model preventing the scalar direct decay into standard model particles. Cosmologically and prior to the electroweak symmetry breaking, the scalar field rolls down along its effective potential before eventually oscillating and settling down at the electroweak minimum. These oscillations can be at the origin of dark matter due to the initial misalignment of the scalar field compared to the electroweak minimum, and we find that, when the mass of the scalar field is less than the eV scale and acts as a condensate behaving like dark matter on large scales, the scalar particles cannot thermalise with the standard model thermal bath. As matter couples in a composition-dependent manner to the oscillating scalar, this could lead to a violation of the equivalence principle aboard satellites such as the MICROSCOPE experiment and the next generation of tests of the equivalence principle. Local gravitational tests are evaded thanks to the weakness of the quadratic coupling in the dark matter halo, and we find that, around other sources, these dilaton models could be subject to a screening akin to the symmetron mechanism.

Black holes with dyonic charges in Einstein-Maxwell-dilaton-axion supergravity theory are revisited in the context of black hole shadows. We consider static as well as rotating (namely the dyonic Kerr-Sen) black holes. The matter stress-energy tensor components, sourced by the Maxwell, axion and dilaton fields satisfy the standard energy conditions. The analytical expressions for the horizon and the shadow radius of the static spacetimes demonstrate their dependence on $P^2+Q^2$ ($P$, $Q$ the magnetic and electric charges, respectively) and the mass parameter $M$. The shadow radius lies in the range $2M <R_{shadow}<3\sqrt{3} M$ and there is no stable photon orbit outside the horizon. Further, shadows cast by the rotating dyonic Kerr-Sen black holes are also studied and compared graphically with their Kerr-Newman and Kerr-Sen counterparts. Deviation of the shadow boundary is prominent with the variation of the magnetic charge, for the relatively slowly rotating dyonic Kerr-Sen spacetimes. We test any possible presence of a magnetic monopole charge in the backdrop of recent EHT observations for the supermassive black holes M87$^*$ and Sgr A$^*$. Deviation from circularity of the shadow boundary ($\Delta C$) and deviation of the average shadow radius from the Schwarzschild shadow radius (quantified as the fractional deviation parameter $\delta$) are the two observables used here. Observational bound on $\Delta C$ (available only for M87$^*$) is satisfied for all theoretically allowed regions of parameter space and thus cannot constrain the parameters. The observational bound on $\delta$ available for Sgr A$^*$ translates into an upper limit on any possible magnetic monopole charge linked to Sgr A$^*$ and is given as $P\lesssim 0.873\, M$. Such a constraint on $P$ is however expected to be far more stringent for other astrophysical tests.

In this article we present an alternative formalism for the inflationary phenomenology of rescaled Einstein-Gauss-Bonnet models which are in agreement with the GW170817 event. By constraining the propagation velocity of primordial tensor perturbations, an approximate form for the time derivative of the scalar field coupled to the Gauss-Bonnet density is extracted. In turn, the overall degrees of freedom decrease and similar to the case of the canonical scalar field, only one scalar function needs to be designated, while the other is extracted from the continuity equation of the scalar field. We showcase explicitly that the slow-roll indices can be written in a closed form as functions of three dimensionless parameters, namely $x=\frac{1}{2\alpha}\bigg(\frac{\kappa\xi'}{\xi''}\bigg)^2$, $\beta=8H^2\xi''$ and $\gamma=\frac{\xi'\xi'''}{\xi''^2}$ and in turn, we prove that the Einstein-Gauss-Bonnet model can in fact produce a blue-tilted tensor spectral index if the condition $\beta\geq1$ is satisfied, which is possible only for Einstein-Gauss-Bonnet models with $\xi''(\phi_k)>0$. Afterwards, a brief comment on the running of the spectral indices is made where it is shown that $a_{\mathcal{S}}(k_*)$ and $a_{\mathcal{T}}(k_*)$ in the constrained case are approximately of the order $\mathcal{O}(10^{-3})$, if not smaller. Last but not least, we examine the conditions under which the Swampland criteria are satisfied. We connect the tracking condition related to scalar field theories with the present models, and we highlight the important feature of the models we propose that the tracking condition can be satisfied only if the Swampland criteria are simultaneously satisfied, however the cases with $\xi\sim1/V$ and $\xi\sim V$ are excluded, as they cannot describe the inflationary era properly.

We study a radiative $p^{15}$N capture on the ground state of $^{16}$O at stellar energies within the framework of a modified potential cluster model (MPCM) with forbidden states, including low lying resonances. The investigation of the $^{15}$N($p,\gamma $)$^{16}$O reaction includes the consideration of $^{3}S_{1}$ resonances due to $E1$ transitions and contribution of $^{3}P_{1}$ scattering wave in $p$ + $^{15}$N channel due to $^{3}P_{1}\longrightarrow $ $^{3}P_{0}$ $M1$ transition. We calculate the astrophysical low-energy $S-$factor and extrapolated $S(0)$ turned out to be within $34.7-40.4$ keV$\cdot $b. It is elucidated the important role of the asymptotic constant (AC) for the $^{15}$N($p,\gamma $)$^{16}$O process with interfering $^{3}S_{1}$(312) and $^{3}S_{1}$(962) resonances. A comparison of our calculation for $S-$factor with existing experimental and theoretical data is addressed and the reasonable agreement is found. The reaction rate is calculated and compared with the existing rates. It has negligible dependence on the variation of AC, but shows strong impact of the interference of $^{3}S_{1}$(312) and $^{3}S_{1}$(962) resonances, especially at $T_{9}$ referring to the CNO Gamow windows. We present a stellar temperature dependence on the Gamow energy and a comparison of rates for radiative proton capture reactions for CNO cycle on nitrogen isotopes obtained in the framework of the MPCM and give temperature windows, prevalence, and significance of each process.

In this paper, we worked in the framework of an inflationary $f(R,T)$ theory, in the presence of a canonical scalar field. More specifically, the $f(R,T)=\gamma R+2\kappa\alpha T$ gravity. The values of the dimensionless parameters $\alpha$ and $\gamma$ are taken to be $\alpha \geq 0$ and $0 < \gamma \leq 1$. The motivation for that study was the striking similarities between the slow-roll parameters of the inflationary model used in this work and the ones obtained by the rescaled Einstein-Hilbert gravity inflation $f(R)=\alpha R$. We examined a variety of potentials to determine if they agree with the current Planck Constraints. In addition, we checked whether these models satisfy the Swampland Criteria and we specified the exact region of the parameter space that produces viable results for each model. As we mention in Section IV the inflationary $f(R,T)$ theory used in this work can not produce a positive $n_T$ which implies that the stochastic gravitational wave background will not be detectable.

We report on the first search for nuclear recoils from dark matter in the form of weakly interacting massive particles (WIMPs) with the XENONnT experiment which is based on a two-phase time projection chamber with a sensitive liquid xenon mass of $5.9$~t. During the approximately 1.1 tonne-year exposure used for this search, the intrinsic $^{85}$Kr and $^{222}$Rn concentrations in the liquid target were reduced to unprecedentedly low levels, giving an electronic recoil background rate of $(15.8\pm1.3)~\mathrm{events}/(\mathrm{t\cdot y \cdot keV})$ in the region of interest. A blind analysis of nuclear recoil events with energies between $3.3$~keV and $60.5$~keV finds no significant excess. This leads to a minimum upper limit on the spin-independent WIMP-nucleon cross section of $2.58\times 10^{-47}~\mathrm{cm}^2$ for a WIMP mass of $28~\mathrm{GeV}/c^2$ at $90\%$ confidence level. Limits for spin-dependent interactions are also provided. Both the limit and the sensitivity for the full range of WIMP masses analyzed here improve on previous results obtained with the XENON1T experiment for the same exposure.

In this paper, we implement a generalized pseudo-Newtonian potential and prescribe a numerical fitting formalism, to study the off-equatorial orbits inclined at a certain angle with the equatorial plane around both Schwarzschild and Kerr-like compact object primaries surrounded by a dipolar halo of matter. The chaotic dynamics of the orbits are detailed for both non-relativistic and special-relativistic test particles. The dependence of the degree of chaos on the rotation parameter $a$ and the inclination angle $i$ is established individually using widely used indicators, such as the Poincar\'e Map and the Lyapunov Characteristic Number. We find that although the chaoticity of the orbits has a positive correlation with $i$, the growth in the chaotic behaviour is not systematic. There exists a threshold value of the inclination angle $i_{\text{c}}$, after which the degree of chaos shows a sharp increase. On the other hand, the chaoticity of the inclined orbits anti-correlates with $a$ at the lower inclination angles. At higher values of $i$, the degree of chaos is maximum for the maximally counter-rotating compact objects, though it has a weak negative, sometimes positive, correlation with $a$ at its higher values. The studies performed with several initial conditions and orbital parameters reveal the intricate nature of the system.

We study the linear polarization from the accretion disk around weakly and strongly naked Janis-Newman-Winicour singularities. We consider an analytical toy model of thin magnetized fluid ring orbiting in the equatorial plane and emitting synchrotron radiation. The observable polarized images are calculated and compared to the Schwarzschild black hole for physical parameters compatible with the radio source M87. For small inclination angles the direct images of the weakly naked singularities closely mimic the Schwarzschild black hole. The deviation in the polarization properties increases if we consider larger inclination angles or higher order images as for indirect images the polarization intensity grows several times in magnitude compared to black holes. Strongly naked singularities produce significant observational signatures already in the direct images. They create a second image of the fluid ring with times larger polarization intensity and characteristic twist of the polarization direction. Due to this additional structure they can be distinguished in polarimetric experiments.

Oscillating scalar field condensates induce small amplitude oscillations of the Hubble parameter which can induce a decay of the condensate due to a parametric resonance instability [1]. We show that this instability can lead to the decay of the coherence of the condensate of axion-like particle (ALP) fields during the radiation phase of standard cosmology for rather generic ALP parameter values, with possible implications for certain experiments aiming to search for ALP candidates. As an example, we study the application of this instability to the QCD axion. We also study the magnitude of the induced entropy fluctuations.