Papers for Monday, Jan 30 2023

Searching for as yet undetected gamma-ray sources is a major target of the Fermi LAT Collaboration. We present an algorithm capable of identifying such type of sources by non-parametrically clustering the directions of arrival of the high-energy photons detected by the telescope onboard the Fermi spacecraft. n particular, the sources will be identified using a von Mises-Fisher kernel estimate of the photon count density on the unit sphere via an adjustment of the mean-shift algorithm to account for the directional nature of data. This choice entails a number of desirable benefits. It allows us to by-pass the difficulties inherent on the borders of any projection of the photon directions onto a 2-dimensional plane, while guaranteeing high flexibility. The smoothing parameter will be chosen adaptively, by combining scientific input with optimal selection guidelines, as known from the literature. Using statistical tools from hypothesis testing and classification, we furthermore present an automatic way to skim off sound candidate sources from the gamma-ray emitting diffuse background and to quantify their significance. The algorithm was calibrated on simulated data provided by the Fermi LAT Collaboration and will be illustrated on a real Fermi LAT case-study.

Charge exchange is an atomic process that primarily occurs at interfaces between the neutral and ionized gas. The study of the process has been carried out on three levels: the theoretical calculation of the cross sections, the laboratory measurements of reaction rates and line strengths, and the observational constraints using celestial objects. For a long time in the past, the status of astrophysical observations in the X-ray band lagged behind the other two aspects until the discovery of X-ray from a comet was made in 1996, which changed the research landscape. Recent observational evidence suggests that charge exchange has been seen or can be expected from a surprisingly broad range of locations, from the Earth's exosphere to the large-scale structures of the Universe. The rapid development of high-resolution X-ray spectroscopy, in particular the non-dispersive micro-calorimeters, is paving the way to revolutionary new science possibilities both in the laboratory and astrophysics. This chapter summarizes the current knowledge of charge exchange and its relevance on astrophysics, especially X-ray spectroscopy.

We present the first homogeneous catalog of $Kepler$, $K2$, and $TESS$ host stars and the corresponding catalog of exoplanet properties, which contain 7993 stars and 9324 planets, respectively. We used isochrone fitting and $Gaia$ DR3 photometry, parallaxes, and spectrophotometric metallicities to compute precise, homogeneous $T_{\mathrm{eff}}$, $\log g$, masses, radii, mean stellar densities, luminosities, ages, distances, and V-band extinctions for 3248, 565, and 4180 $Kepler$, $K2$, and $TESS$ stars, respectively. We compared our stellar properties to studies using fundamental and precise constraints, such as interferometry and asteroseismology, and find residual scatters of 2.8%, 5.6%, 5.0%, and 31%, with offsets of 0.2%, 1.0%, 1.2%, and 0.7% between our $T_{\mathrm{eff}}$, radii, masses, and ages and those in the literature, respectively. In addition, we compute planet radii, semimajor axes, and incident fluxes for 4281, 676, and 4367 $Kepler$, $K2$, and $TESS$ planets, respectively, and find that the exoplanet radius gap is less prominent in the $K2$, $TESS$, and combined samples than it is in the $Kepler$ sample alone. We suspect this difference is largely due to heterogeneous planet-to-star radius ratios, shorter time baselines of $K2$ and $TESS$, and smaller sample sizes. Finally, we identify a clear radius inflation trend in our large sample of hot Jupiters and find 150 hot sub-Neptunian desert planets, in addition to a population of over 400 young host stars as potential opportunities for testing theories of planet formation and evolution.

Apsidal precession in stellar binaries is the main non-Keplerian dynamical effect impacting the radial-velocities of a binary star system. Its presence can notably hide the presence of orbiting circumbinary planets because many fitting algorithms assume perfectly Keplerian motion. To first order, apsidal precession ($\dot{\omega}$) can be accounted for by adding a linear term to the usual Keplerian model. We include apsidal precession in the kima package, an orbital fitter designed to detect and characterise planets from radial velocity data. In this paper, we detail this and other additions to kima that improve fitting for stellar binaries and circumbinary planets including corrections from general relativity. We then demonstrate that fitting for $\dot{\omega}$ can improve the detection sensitivity to circumbinary exoplanets by up to an order of magnitude in some circumstances, particularly in the case of multi-planetary systems. In addition, we apply the algorithm to several real systems, producing a new measurement of aspidal precession in KOI-126 (a tight triple system), and a detection of $\dot{\omega}$ in the Kepler-16 circumbinary system. Although apsidal precession is detected for Kepler-16, it does not have a large effect on the detection limit or the planetary parameters. We also derive an expression for the precession an outer planet would induce on the inner binary and compare the value this predicts with the one we detect.

Electron impact ionization is critical in producing the ionospheres on many planetary bodies and, as discussed here, is critical for interpreting spacecraft and telescopic observations of the tenuous atmospheres of the icy Galilean satellites of Jupiter (Europa, Ganymede, and Callisto), which form an interesting planetary system. Fortunately, laboratory measurements, extrapolated by theoretical models, were developed and published over a number of years by K. Becker and colleagues (see Deutsch et al. 2009) to provide accurate electron impact ionization cross sections for atoms and molecules, which are crucial to correctly interpret these measurements. Because of their relevance for the Jovian icy satellites we provide useful fits to the complex, semi-empirical Deutsch-Mark formula for energy-dependent electron impact ionization cross-sections of gas-phase water products (i.e., H2O, H2, O2, H, O). These are then used with measurements of the thermal plasma in the Jovian magnetosphere to produce ionization rates for comparison with solar photo-ionization rates at the icy Galilean satellites.

We establish a dynamical mechanism to explain the origin of the asymmetry between the arms observed in some barred disk galaxies, where one of the two arms emanating from the bar ends is very well defined, while the second one displays a ragged structure, extending between its ridge and the bar. To this purpose, we study the invariant manifolds associated to the Lyapunov periodic orbits around the unstable equilibrium points at the ends of the bar. Matter from the galaxy center is transported along these manifolds to the periphery, forming this way the spiral arms that emanate from the bar ends. If the mass distribution in the galaxy center is not homogeneous, because of an asymmetric bar with one side stronger than the other, or because of a non-centered bulge, the dynamics about the two unstable Lagrange points at the ends of the bar will not be symmetric as well. One of their invariant manifolds becomes more extended than the other, enclosing a smaller section and the escaping orbits on it are fewer and dispersed in a wider region. The result is a weaker arm, and more ragged than the one at the other end of the bar.

We report new measurements of branching fractions for 20 UV and blue lines in the spectrum of neutral silicon (Si I) originating in the 3$s^{2}$3$p$4$s$ $^{3}$P$^{\rm o}_{1,2}$, $^{1}$P$^{\rm o}_{1}$ and 3$s$3$p^{3}$ $^{1}$D$^{\rm o}_{1,2}$ upper levels. Transitions studied include both strong, nearly pure LS multiplets as well as very weak spin-forbidden transitions connected to these upper levels. We also report a new branching fraction measurement of the $^{4}$P$_{1/2}$ - $^{2}$P$^{\rm o}_{1/2,3/2}$ intercombination lines in the spectrum of singly-ionized silicon (Si II). The weak spin-forbidden lines of Si I and Si II provide a stringent test on recent theoretical calculations, to which we make comparison. The branching fractions from this study are combined with previously reported radiative lifetimes to yield transition probabilities and log($gf$)s for these lines. We apply these new measurements to abundance determinations in five metal-poor stars.

Globular clusters exhibit abundance variations, defining `multiple populations', which have prompted a protracted search for their origin. Properties requiring explanation include: the high fraction of polluted stars ($\sim 40{-}90$~percent, correlated with cluster mass), the absence of pollution in young clusters and the lower pollution rate with binarity and distance from the cluster centre. We present a novel mechanism for late delivery of pollutants into stars via accretion of sub-stellar companions. In this scenario, stars move through a medium polluted with AGB and massive star ejecta, accreting material to produce companions with typical mass ratio $q\sim 0.1$. These companions undergo eccentricity excitation due to dynamical perturbations by passing stars, culminating in a merger with their host star. The accretion of the companion alters surface abundances via injected pollutant. Alongside other self-enrichment models, the companion accretion model can explain the dilution of pollutant and correlation with intra-cluster location. The model also explains the ubiquity and discreteness of the populations and correlations of enrichment rates with cluster mass, cluster age and stellar binarity. Abundance variations in some clusters can be broadly reproduced using AGB and massive binary ejecta abundances from the literature. In other clusters, some high companion mass ratios ($q\gtrsim 1$) are required. In these cases, the available mass budget necessitates a variable degree of mixing of the polluted material with the primary star, deviations from model ejecta abundances or mixing of internal burning products. We highlight the avenues of further investigation which are required to explore some of the key processes invoked in this model.

We report the spectroscopic confirmation of a massive quiescent galaxy, GS-9209 at a new redshift record of $z=4.658$, just 1.25 Gyr after the Big Bang, using new deep continuum observations from JWST NIRSpec. From our full-spectral-fitting analysis, we find that this galaxy formed its stellar population over a $\simeq200$ Myr period, approximately $600-800$ Myr after the Big Bang ($z_\rm{form}=7.3\pm0.2$), before quenching at $z_\rm{quench}=6.7\pm0.3$. GS-9209 demonstrates unambiguously that massive galaxy formation was already well underway within the first billion years of cosmic history, with this object having reached a stellar mass of log$_{10}(M_*/\rm{M}_\odot) > 10.3$ by $z=7$. This galaxy also clearly demonstrates that the earliest onset of galaxy quenching was no later than $\simeq800$ Myr after the Big Bang. We estimate the iron abundance and $\alpha$-enhancement, finding [Fe/H] $= -0.97^{+0.06}_{-0.07}$ and [$\alpha$/Fe] $= 0.67^{+0.25}_{-0.15}$, suggesting the stellar mass vs iron abundance relation at $z\simeq7$, when this object formed most of its stars, was $\simeq0.4$ dex lower than at $z\simeq3.5$. Whilst its spectrum is dominated by stellar emission, GS-9209 also exhibits broad H$\alpha$ emission, indicating that it hosts an active galactic nucleus (AGN), for which we measure a black-hole mass of log$_{10}(M_\bullet/\rm{M}_\odot) = 8.7\pm0.1$. Although large-scale star formation in GS-9209 has been quenched for almost half a billion years, the significant integrated quantity of accretion implied by this large black-hole mass suggests AGN feedback plausibly played a significant role in quenching star formation in this galaxy. GS-9209 is also extremely compact, with an effective radius of just $215\pm20$ parsecs. This intriguing object offers perhaps our deepest insight yet into massive galaxy formation and quenching during the first billion years of cosmic history.

We present 18 years of OGLE photometry together with spectra obtained over 12 years, revealing that the early Oe star AzV 493 shows strong photometric (Delta I < 1.2 mag) and spectroscopic variability with a dominant, 14.6-year pattern and ~40-day oscillations. We estimate stellar parameters T_eff = 42000 K, log L/L_sun = 5.83 +/- 0.15, M/M_sun = 50 +/- 9, and vsini = 370 +/- 40 km/s. Direct spectroscopic evidence shows episodes of both gas ejection and infall. There is no X-ray detection, and it is likely a runaway star. AzV 493 may have an unseen companion on a highly eccentric (e > 0.93) orbit. We propose that close interaction at periastron excites ejection of the decretion disk, whose variable emission-line spectrum suggests separate inner and outer components, with an optically thick outer component obscuring both the stellar photosphere and the emission-line spectrum of the inner disk at early phases in the photometric cycle. It is plausible that AzV 493's mass and rotation have been enhanced by binary interaction followed by the core-collapse supernova explosion of the companion, which now could be either a black hole or neutron star. This system in the Small Magellanic Cloud can potentially shed light on OBe decretion disk formation and evolution, massive binary evolution, and compact binary progenitors.

The indirect detection of cosmic rays via the radio signal of extensive air showers is gaining a lot of ground. Many new arrays of radio antennas are under construction or in the phase of development. Calibrating these arrays is important for the reconstruction of observed events and for the comparability between observatories. Using reference antennas in calibration campaigns is not ideal because of uncertainties on their signal output strength that are large or difficult to assess. In a different approach the arrays can be calibrated against the Galactic radio emission as the dominant source of background. This so-called Galactic Calibration relies on predictions of the diffuse Galactic radio emission, for which models are publicly available. We present a comparison of these models in the frequency range from 10 to 408 MHz in order to estimate the systematic uncertainties on the strength of the Galactic background. We do this comparison on a global level as well as adapted for selected radio arrays and discuss implications for applying the Galactic calibration method. Furthermore we study the influence of the quiet Sun as an additional source of radio emission in the sky.

We use the RIOTS4 sample of SMC field OB stars to determine the origin of massive runaways in this low-metallicity galaxy using Gaia proper motions, together with stellar masses obtained from RIOTS4 data. These data allow us to estimate the relative contributions of stars accelerated by the dynamical ejection vs binary supernova mechanisms, since dynamical ejection favors faster, more massive runaways, while SN ejection favors the opposite trend. In addition, we use the frequencies of classical OBe stars, high-mass X-ray binaries, and non-compact binaries to discriminate between the mechanisms. Our results show that the dynamical mechanism dominates by a factor of 2 - 3. This also implies a significant contribution from two-step acceleration that occurs when dynamically ejected binaries are followed by SN kicks. We update our published quantitative results from Gaia DR2 proper motions with new data from DR3.

We observed HD 19467 B with JWST's NIRCam in six filters spanning 2.5-4.6 $\mu m$ with the Long Wavelength Bar coronagraph. The brown dwarf HD 19467 B was initially identified through a long-period trend in the radial velocity of G3V star HD 19467. HD 19467 B was subsequently detected via coronagraphic imaging and spectroscopy, and characterized as a late-T type brown dwarf with approximate temperature $\sim1000$K. We observed HD 19467 B as a part of the NIRCam GTO science program, demonstrating the first use of the NIRCam Long Wavelength Bar coronagraphic mask. The object was detected in all 6 filters (contrast levels of $2\times10^{-4}$ to $2\times10^{-5}$) at a separation of 1.6 arcsec using Angular Differential Imaging (ADI) and Synthetic Reference Differential Imaging (SynRDI). Due to a guidestar failure during acquisition of a pre-selected reference star, no reference star data was available for post-processing. However, RDI was successfully applied using synthetic Point Spread Functions (PSFs) developed from contemporaneous maps of the telescope's optical configuration. Additional radial velocity data (from Keck/HIRES) are used to constrain the orbit of HD 19467 B. Photometric data from TESS are used to constrain the properties of the host star, particularly its age. NIRCam photometry, spectra and photometry from literature, and improved stellar parameters are used in conjunction with recent spectral and evolutionary substellar models to derive physical properties for HD 19467 B. Using an age of 9.4$\pm$0.9 Gyr inferred from spectroscopy, Gaia astrometry, and TESS asteroseismology, we obtain a model-derived mass of 62$\pm 1M_{J}$, which is consistent within 2-$\sigma$ with the dynamically derived mass of 81$^{+14}_{-12}M_{J}$.

We present predictions for cosmic evolution of populations of supermassive black holes (SMBHs) forming from Population III.1 seeds, i.e., early, metal-free dark matter minihalos forming far from other sources, parameterized by isolation distance, $d_{\rm{iso}}$. Extending previous work that explored this scenario to $z=10$, we follow evolution of a $(60\:{\rm{Mpc}})^3$ volume to $z=0$. We focus on evolution of SMBH comoving number densities, halo occupation fractions, angular clustering and 3D clustering, exploring a range of $d_{\rm{iso}}$ constrained by observed local number densities of SMBHs. We also compute synthetic projected observational fields, in particular a case comparable to the Hubble Ultra Deep Field. We compare Pop III.1 seeding to a simple halo mass threshold model, commonly adopted in cosmological simulations of galaxy formation. Major predictions of the Pop III.1 model include that all SMBHs form by $z\sim25$, after which their comoving number densities are near-constant, with low merger rates. Occupation fractions evolve to concentrate SMBHs in the most massive halos by $z=0$, but with rare cases in halos down to $\sim10^8\:M_\odot$. The $d_{\rm{iso}}$ scale at epoch of formation, e.g., $100\:$kpc-proper at $z\sim30$, i.e., $\sim3\:$Mpc-comoving, is imprinted in the SMBH two-point angular correlation function, remaining discernible as a low-amplitude feature to $z\sim1$. The SMBH 3D two-point correlation function at $z=0$ also shows lower amplitude compared to equivalently massive halos. We discuss prospects for testing these predictions with observational surveys of SMBH populations.

Galaxy formation and evolution are regulated by the feedback from galactic winds. Absorption lines provide the most widely available probe of winds. However, since most data only provide information integrated along the line-of-sight, they do not directly constrain the radial structure of the outflows. In this paper, we present a method to directly measure the gas electron density in outflows (ne), which in turn yields estimates of outflow cloud properties (e.g., density, volume filling-factor, and sizes/masses). We also estimate the distance (r) from the starburst at which the observed densities are found. We focus on 22 local star-forming galaxies primarily from the COS Legacy Archive Spectroscopic SurveY (CLASSY). In half of them, we detect absorption lines from fine structure excited transitions of Si II (i.e., Si II*). We determine ne from relative column densities of Si II and Si II*, given Si II* originates from collisional excitation by free electrons. We find that the derived ne correlates well with the galaxy's star-formation rate per unit area. From photoionization models or assuming the outflow is in pressure equilibrium with the wind fluid, we get r ~ 1 to 2 * rstar or ~ 5 * rstar, respectively, where rstar is the starburst radius. Based on comparisons to theoretical models of multi-phase outflows, nearly all of the outflows have cloud sizes large enough for the clouds to survive their interaction with the hot wind fluid. Most of these measurements are the first-ever for galactic winds detected in absorption lines and, thus, will provide important constraints for future models of galactic winds.

The Seoul Radio Astronomy Observatory (SRAO) operates a 6.1-meter radio telescope on the Gwanak campus of Seoul National University. We present the efforts to reform SRAO to a Very Long Baseline Interferometry (VLBI) station, motivated by recent achievements by millimeter interferometer networks such as Event Horizon Telescope, East Asia VLBI Network, and Korean VLBI Network (KVN). For this goal, we installed a receiver that had been used in the Combined Array for Research in Millimeter-wave Astronomy and a digital backend, including an H-maser clock. The existing hardware and software were also revised, which had been dedicated only to single-dish operations. After several years of preparations and test observations in 1 and 3-millimeter bands, a fringe was successfully detected toward 3C 84 in 86 GHz in June 2022 for a baseline between SRAO and KVN Ulsan station separated by 300 km. Thanks to the dual frequency operation of the receiver, the VLBI observations will soon be extended to the 1 mm band and verify the frequency phase referencing technique between 1 and 3-millimeter bands.

Surfaces of planetary bodies can have strong electric fields, subjecting conductive grains to repulsive electrostatic forces. This has been proposed as mechanism to eject grains from the ground. To quantify this process, we study mm-sized basalt aggregates consisting of micrometer constituents exposed to an electric field in drop tower experiments. The dust aggregates acquire high charges on sub-second timescales while sticking to the electrodes according to the field polarity. Charging at the electrodes results in a repulsive (lifting) force and continues until repulsion overcomes adhesion and particles are lifted, moving towards the opposite electrode. Some aggregates remain attached, which is consistent with a maximum charge limit being reached, providing an electrostatic force too small to counteract adhesion. All observations are in agreement with a model of moderately conductive grains with a small but varying number of adhesive contacts to the electrodes. This supports the idea that on planetary surfaces with atmospheres, electrostatic repulsion can significantly contribute to airborne dust and sand, i.e. decrease the threshold wind speed that is required for saltation and increase the particle flux as suggested before.

Methods. We use the code MCFOST to solve the non-LTE problem of line formation in non-axisymmetric accreting magnetospheres. We compute the Br{\gamma} line profile originating from accretion columns for models with different magnetic obliquities. We also derive monochromatic synthetic images of the Br{\gamma} line emitting region across the line profile. This spectral line is a prime diagnostics of magnetospheric accretion in young stars and is accessible with the long baseline near-infrared interferometer GRAVITY installed at the ESO Very Large Telescope Interferometer. Results. We derive Br{\gamma} line profiles as a function of rotational phase and compute interferometric observables, visibilities and phases, from synthetic images. The line profile shape is modulated along the rotational cycle, exhibiting inverse P Cygni profiles at the time the accretion shock faces the observer. The size of the line's emission region decreases as the magnetic obliquity increases, which is reflected in a lower line flux. We apply interferometric models to the synthetic visibilities in order to derive the size of the line-emitting region. We find the derived interferometric size to be more compact than the actual size of the magnetosphere, ranging from 50 to 90\% of the truncation radius. Additionally, we show that the rotation of the non-axisymmetric magnetosphere is recovered from the rotational modulation of the Br{\gamma}-to-continuum photo-centre shifts, as measured by the differential phase of interferometric visibilities.

Very long baseline interferometry (VLBI) is a radio-astronomical technique in which the correlated signal from various baselines is combined into an image of highest angular resolution. Due to sparsity of the measurements, this imaging procedure constitutes an ill-posed inverse problem. For decades the CLEAN algorithm was the standard choice in VLBI studies, although having some serious disadvantages and pathologies that are challenged by the requirements of modern frontline VLBI applications. We develop a novel multi-scale CLEAN deconvolution method (DoB-CLEAN) based on continuous wavelet transforms that address several pathologies in CLEAN imaging. We benchmark this novel algorithm against CLEAN reconstructions on synthetic data and reanalyze BL Lac observations of RadioAstron with DoB-CLEAN. DoB-CLEAN approaches the image by multi-scalar and multi-directional wavelet dictionaries. Two different dictionaries are used. Firstly, a difference of elliptical spherical Bessel functions dictionary fitted to the uv-coverage of the observation that is used to sparsely represent the features in the dirty image. Secondly, a difference of elliptical Gaussian wavelet dictionary that is well suited to represent relevant image features cleanly. The deconvolution is performed by switching between the dictionaries. DoB-CLEAN achieves super-resolution compared to CLEAN and remedies the spurious regularization properties of CLEAN. In contrast to CLEAN, the representation by basis functions has a physical meaning. Hence, the computed deconvolved image still fits the observed visibilities, opposed to CLEAN. State-of-the-art multi-scalar imaging approaches seem to outperform single-scalar standard approaches in VLBI and are well suited to maximize the extraction of information in ongoing frontline VLBI applications.

It has been shown that during the outburst of accretion activity observed in UX Ori type star RZ Psc in 2013, the accretion rate increased approximately by an order of magnitude. This means that the accretion process at the late stages of the Pre-Main Sequence evolution is very unstable. Using the spectra obtained during this episode we have studied the magnetospheric emission in the H$\alpha$ line. Models of magnetospheric accretion are calculated to obtain the parameters of the magnetosphere from this observation. In present work we have taken into account the influence of the recombination delay effect during gas motion in the stellar magnetosphere. The accounting for this effect and the presence of the magnetospheric absorption in the IR CaII triplet lines and its absence in D Na I resonance lines allowed us to place a lower limit on the temperature in the magnetosphere at $\approx$ 10000 K, which significantly improved precision of our estimate of accretion rate. According to the best fit model the logarithm of accretion rate is $\log\dot{M} = -10.1\pm0.3$ ($\dot{M} \approx 7\times10^{-11}$ M$_\odot$yr$^{-1}$) and the inclination angle of RZ Psc is $43\pm 3^\circ$. It is less than the inclination, typical for the UX Ori stars (about 70$^\circ$), that explains the weak photometric variability of this star. Using the obtained accretion rate and magnetosphere radius we estimate the strength of the dipole component of the magnetic field of RZ Psc $\approx$ 0.1 kGs.

We investigate the flow of material from highly misaligned and polar circumbinary discs that feed the formation of circumstellar discs around each binary component. With three-dimensional hydrodynamic simulations we consider equal mass binaries with low eccentricity. We also simulate inclined test particles and highly-misaligned circumstellar discs around one binary component for comparison. During Kozai-Lidov (KL) cycles, the circumstellar disc structure is altered through exchanges of disc eccentricity with disc tilt. Highly inclined circumstellar discs and test particles around individual binary components can experience very strong KL oscillations. The continuous accretion of highly misaligned material from the circumbinary disc allows the KL oscillations of circumstellar discs to be long-lived. In this process, the circumbinary material is continuously delivered with a high inclination to the lower inclination circumstellar discs. We find that the simulation resolution is important for modeling the longevity of the KL oscillations. An initially polar circumbinary disc forms nearly polar, circumstellar discs that undergo KL cycles. The gas steams accreting onto the polar circumstellar discs vary in tilt during each binary orbital period, which determines how much material is accreted onto the discs. The long-lived KL cycles in polar circumstellar discs may lead to the formation of polar S-type planets in binary star systems.

Extending previous results [JHEP 11 (2021) 091], we explore aspects of the reheating mechanism for non-minimal Higgs inflation in the strong coupling regime. We constrain the radiative corrections for the inflaton's potential by considering the Coleman-Weinberg approximation and use the Renormalization Group Equations for the Higgs field to derive an upper limit on the quark top mass, $m_t$. Using the current Cosmic Microwave Background, Baryon Acoustic Oscillation, and Supernova data, we obtain $m_t \leq 170.44$ GeV, confirming the observational compatibility of the model with recent $m_t$ estimates reported by the CMS collaboration. We also analyze the breakdown of the well-known correlation involving the Hubble constant $H_0$ and the clustering parameter $\sigma_8$, which makes the model interesting in light of the cosmological tensions discussed over the last decade.

Particle cascades induced by ultra-high-energy (UHE) cosmic rays and neutrinos impacting on the lunar regolith usually radiate Cherenkov radio emissions due to the presence of excess negative charge, which is known as Askaryan effect. Several experiments have been carried out to detect the Cherenkov radio emissions in the lunar regolith. To prepare for future lunar Ultra-Long Wavelength (ULW, frequencies below 30 MHz) radio astronomy missions, we study the detection of the Cherenkov radio emissions with the ULW radio telescope that are operating at the lunar orbit. We have carried out instrument modelling and analytic calculations for the analysis of aperture, flux and event rate, and the analyses show the detectability of the Cherenkov radiation. Based on the properties of the Cherenkov radiation, we have demonstrated that the cosmic ray and neutrino events could be reconstructed with the three ULW vector antennas onboard the lunar satellites via measurements of the Askaryan radio pulse intensity, polarizations, etc. The results obtained by this study would be useful for future lunar radio explorer mission, where the detections of UHE cosmic rays and neutrinos could be successfully attempted.

Convection and mass loss by stellar winds are two dynamical processes that shape asymptotic giant branch (AGB) stars and their evolution. Observations and earlier 3D models indicate that giant convection cells cause high-contrast surface intensity patterns, and contribute to the origin of clumpy dust clouds. We study the formation and resulting properties of dust-driven winds from AGB stars, using new global 3D simulations. The dynamical stellar interiors, atmospheres, and wind acceleration zones of two M-type AGB stars were modeled with the CO5BOLD code. These first global 3D simulations are based on frequency-dependent gas opacities, and they feature time-dependent condensation and evaporation of silicate grains. Convection and pulsations emerge self-consistently, allowing us to derive wind properties (e.g., mass-loss rates and outflow velocities), without relying on parameterized descriptions of these processes. In contrast to 1D models with purely radial pulsations, the shocks induced by convection and pulsation in the 3D models cover large parts, but not the entirety, of the sphere, leading to a patchy, nonspherical structure of the atmosphere. Since dust condensation critically depends on gas density, new dust clouds form mostly in the dense wakes of atmospheric shocks, where the grains can grow efficiently. The resulting clumpy distribution of newly formed dust leads to a complex 3D morphology of the extended atmosphere and wind-acceleration zone, with simultaneous infall and outflow regions close to the star. Highly nonspherical isotherms and short-lived cool pockets of gas in the stellar vicinity are prominent features. Efficient dust formation sets in closer to the star than spherical averages of the temperature indicate, in dense regions where grain growth rates are higher than average.

We present the merging of the Particle-Mesh (PM) relativistic Gevolution code with the TreePM Gadget-4 code, with the aim of studying general relativity effects in cosmology. Our code, called GrGadget, is able to track the evolution of metric perturbations in the weak field limit by using Gevolution's implementation of a relativistic PM in the Poisson gauge. To achieve this, starting from Gevolution we have written a C++ library called libgevolution, that allows a code to access and use the same abstractions and resources that Gevolution uses for its PM-only N-body simulations. The code works under the assumption that particle interactions at short distances can be approximated as Newtonian, so that we can combine the forces computed with a Newtonian Tree with those computed with a relativistic PM. The result is a TreePM simulation code that represents metric perturbations at the scales where they are relevant, while resolving non-linear structures. We validate our code by closely matching Gadget-4 forces, computed with the Tree switched off, with those computed with libgevolution in the Newtonian limit. With GrGadget we obtain a matter power spectrum that is compatible with Newtonian Gadget at small scales and contains GR features at large scales that are consistent with results obtained with Gevolution. We demonstrate that, due to the better resolution of the highly non-linear regime, the representation of the relativistic fields sampled on the mesh improves with respect to the PM-only simulations.

We present results from new self-consistent 3D MHD simulations of the magnetospheres from massive stars with a dipole magnetic axis that has a non-zero obliquity angle ($\beta$) to the star's rotation axis. As an initial direct application, we compare the global structure of co-rotating disks for nearly aligned ($\beta=5^o$) versus half-oblique ($\beta=45^o$) models, both with moderately rapid rotation ($\sim$ 0.5 critical). We find that accumulation surfaces broadly resemble the forms predicted by the analytic Rigidly Rotating Magnetosphere (RRM) model, but the mass buildup to near the critical level for centrifugal breakout against magnetic confinement distorts the field from the imposed initial dipole. This leads to an associated warping of the accumulation surface toward the rotational equator, with the highest density concentrated in {\em wings} centered on the intersection between the magnetic and rotational equators. These MHD models can be used to synthesize rotational modulation of photometric absorption and H$\alpha$ emission for a direct comparison with observations.

While the sample of optical Type Ia Supernova (SN Ia) light curves (LCs) usable for cosmological parameter measurements surpasses 2000, the sample of published, cosmologically viable near-infrared (NIR) SN Ia LCs, which have been shown to be good "standard candles," is still $\lesssim$ 200. Here, we present high-quality NIR LCs for 83 SNe Ia ranging from $0.002 < z < 0.09$ as a part of the Dark Energy, H$_0$, and peculiar Velocities using Infrared Light from Supernovae (DEHVILS) survey. Observations are taken using UKIRT's WFCAM, where the median depth of the images is 20.7, 20.1, and 19.3 mag (Vega) for $Y$, $J$, and $H$-bands, respectively. The median number of epochs per SN Ia is 18 for all three bands ($YJH$) combined and 6 for each band individually. We fit 47 SN Ia LCs that pass strict quality cuts using three LC models, SALT3, SNooPy, and BayeSN and find scatter on the Hubble diagram to be comparable to or better than scatter from optical-only fits in the literature. Fitting NIR-only LCs, we obtain standard deviations ranging from 0.128-0.135 mag. Additionally, we present a refined calibration method for transforming 2MASS magnitudes to WFCAM magnitudes using HST CALSPEC stars that results in a 0.03 mag shift in the WFCAM $Y$-band magnitudes.

Sagittarius A*, the supermassive black hole at the center of our galaxy, exhibits episodic near-infrared flares. The recent monitoring of three such events by the GRAVITY instrument has shown that some flares are associated with orbital motions in the close environment of the black hole with super Keplerian velocity. We develop a semi-analytic model of Sagittarius~A* flares based on an ejected large plasmoid, inspired by recent particle-in-cell global simulations of black hole magnetospheres. We model the infrared astrometric and photometric signatures associated to this model. We consider a spherical large plasmoid ejected along a conical orbit around the black hole. This plasmoid is assumed to be formed by successive mergers of smaller plasmoids produced through magnetic reconnection. Non-thermal electrons are injected in the plasmoid. We compute the evolution of the electron-distribution under the influence of synchrotron cooling. We solve the radiative transfer problem and transport the radiation along null geodesics of the Schwarzschild spacetime. We also take into account the quiescent radiation of the accretion flow, on top of which the flare evolves. For the first time, we successfully account for the astrometric and flux variations of the GRAVITY data with a flare model that incorporates an explicit modeling of the emission mechanism. We find good agreement between the prediction of our model and the recent data. In particular, the azimuthal velocity is set by the magnetic field line it belongs to, which is anchored in the inner parts of the accretion flow, hence the super-Keplerian motion. The astrometric track is also shifted with respect to the center of mass due to the quiescent radiation, in agreement with the difference measured with the GRAVITY data. These results support the picture of magnetic reconnection as a viable model for Sagittarius~A* infrared flares.

The analysis of photometric large-scale structure data is often complicated by the need to account for many observational and astrophysical systematics. The elaborate models needed to describe them often introduce many ``nuisance parameters'', which can be a major inhibitor of an efficient parameter inference. In this paper we introduce an approximate method to analytically marginalise over a large number of nuisance parameters based on the Laplace approximation. We discuss the mathematics of the method, its relation to concepts such as volume effects and profile likelihood, and show that it can be further simplified for calibratable systematics by linearising the dependence of the theory on the associated parameters. We quantify the accuracy of this approach by comparing it with traditional sampling methods in the context of existing data from the Dark Energy Survey, as well as futuristic Stage-IV photometric data. The linearised version of the method is able to obtain parameter constraints that are virtually equivalent to those found by exploring the full parameter space for a large number of calibratable nuisance parameters, while reducing the computation time by a factor 3-10. Furthermore, the non-linearised approach is able to analytically marginalise over a large number of parameters, returning constraints that are virtually indistinguishable from the brute-force method in most cases, accurately reproducing both the marginalised uncertainty on cosmological parameters, and the impact of volume effects associated with this marginalisation. We provide simple recipes to diagnose when the approximations made by the method fail and one should thus resort to traditional methods. The gains in sampling efficiency associated with this method enable the joint analysis of multiple surveys, typically hindered by the large number of nuisance parameters needed to describe them.

Studying solar wind conditions is central to forecasting impact of space weather on Earth. Under the assumption that the structure of this wind is constant in time and corotates with the Sun, solar wind and thereby space weather forecasts have been made quite effectively. Such corotation forecasts are well studied with decades of observations from STEREO and near-Earth spacecrafts. Forecast accuracy depends upon the latitudinal separation (or offset $\Delta \theta$) between source and spacecraft, forecast lead time ($\Delta t$) and the solar cycle via the sunspot number (SSN). The precise dependencies factoring in uncertain- ties however, are a mixture of influences from each of these factors. And for high precision forecasts, it is important to understand what drives the forecast accuracy and its uncertainty. Here we present a causal inference approach based on information theoretic measures to do this. Our framework can compute not only the direct (linear and non-linear) dependencies of the forecast mean absolute error (MAE) on SSN, $\Delta t$ and $\Delta t$, but also how these individual variables combine to enhance or diminish the MAE. We provide an initial assessment of this with potential of aiding data assimilation in the future.

Strong magnetic fields are observed in a substantial fraction of upper main sequence stars and white dwarfs. Many such stars are observed to exhibit photometric modulations as the magnetic poles rotate in and out of view, which could be a consequence of magnetic perturbations to the star's thermal structure. The magnetic pressure is typically larger than the gas pressure at the star's photosphere, but much smaller than the gas pressure in the star's interior, so the expected surface flux perturbations are not clear. We compute magnetically perturbed stellar structures of young $3 \, M_\odot$ stars that are in both hydrostatic and thermal equilibrium, and which contain both poloidal and toroidal components of a dipolar magnetic field as expected for stable fossil fields. This provides semi-analytical models of such fields in baroclinic stably stratified regions. The star's internal pressure, temperature, and flux perturbations can have a range of magnitudes, though we argue the most likely configurations exhibit flux perturbations much smaller than the ratio of surface magnetic pressure to surface gas pressure, but much larger than the ratio of surface magnetic pressure to central gas pressure. The magnetic pole is hotter than the equator in our models, but a cooler magnetic pole is possible depending on the magnetic field configuration. The expected flux variations for observed field strengths are $\delta L/L \! \lesssim \! 10^{-6}$, much smaller than those observed in magnetic stars, suggesting that observed perturbations stem from changes to the emergent spectrum rather than changes to the bolometric flux.

When the scale factor of expansion of the universe is written as $ a(t)\equiv Ae^{\alpha(t)}$, with $A$ as some real constant and $\alpha(t)$ a real function, the gravitational action $I_G$ appears in the same form as the matter action $I_M$ for a homogeneous and isotropic scalar field with a specific scale factor-dependent potential. We observe that by making analytic continuation of the Lagrangian function in this $I_G$ to the complex plane of $\alpha$, one can obtain terms corresponding to both parts of a total action that describes a viable cosmological model. The solution gives a bouncing and coasting universe, which first contracts and then expands linearly, with a smooth bounce in between them. Such a bounce, called a Tolman wormhole, is due to a Casimir-like negative energy that appears naturally in the solution. In a novel approach to quantum cosmology, we perform canonical quantization of the above model using the operator method and show that it can circumvent the operator ordering ambiguity in the conventional formalism. It leads to a quantum wave equation for the universe, solving which we get the interesting result that the universe is in a ground state with nonzero energy. This solution is different from the Wheeler-DeWitt wave function and possesses exact quantum-classical correspondence during all epochs.

We revisit graviton production via Bremsstrahlung from the decay of the inflaton during inflationary reheating. Using two complementary computational techniques, we first show that such 3-body differential decay rates differ from previously reported results in the literature. We then compute the stochastic gravitational wave (GW) background that forms during the period of reheating, when the inflaton perturbatively decays with the radiative emission of gravitons. By computing the number of relativistic degrees of freedom in terms of $\Delta N_\text{eff}$, we constrain the resulting GW energy density from BBN and CMB. Finally, we project current and future GW detector sensitivities in probing such a stochastic GW background, which typically peaks in the GHz to THz ballpark, opening up the opportunity to be detected with microwave cavities and space-based GW detectors.

Primordial black holes form in the early Universe and constitute one of the most viable candidates for dark matter. The study of their formation process requires the determination of a critical energy density perturbation threshold $\delta_\mathrm{c}$, which in general depends on the underlying gravity theory. Up to now, the majority of analytic and numerical techniques calculate $\delta_\mathrm{c}$ within the framework of general relativity. In this work, using simple physical arguments we estimate semi-analytically the PBH formation threshold within the framework of quantum gravity, working for concreteness within loop quantum gravity (LQG), which constitutes a non-perturbative and background-independent quantization of general relativity. In particular, for low mass PBHs formed close to the quantum bounce, we find a reduction in the value of $\delta_\mathrm{c}$ up to $50\%$ compared to the general relativistic regime quantifying for the first time to the best of our knowledge how quantum effects can influence PBH formation within a quantum gravity framework. Finally, by varying the Barbero-Immirzi parameter $\gamma$ of LQG we show its effect on the value of $\delta_\mathrm{c}$ while using the observational/phenomenological signatures associated to ultra-light PBHs, namely the ones affected by LQG effects, we propose the PBH portal as a novel probe to constrain the potential quantum nature of gravity.

Wave-like dark matter made of spin-1 particles (dark photons) is expected to form ground state clumps called "vector solitons", which can have different polarizations. In this work, we consider the interaction of dark photons with photons, expressed as dimension-6 operators, and study the electromagnetic radiation that arises from an isolated vector soliton due to parametric resonant amplification of the ambient electromagnetic field. We characterize the directional dependence and polarization of the outgoing radiation, which depends on the operator as well as the polarization state of the underlying vector soliton. We discuss the implications of this radiation for the stability of solitons and as a possible channel for detecting mergers of vector solitons through astrophysical observations.

We perform an exhaustive follow-up analysis of a subsolar-mass black hole candidate from the second observing run of Advanced LIGO, reported by Phukon et al. in 2021. The origin of this trigger is unclear, because the reported signal-to-noise ratio (SNR) of $8.6$ and inverse false alarm rate of about $0.5$ yr are too low to claim a gravitational-wave origin, but large enough to be intriguing. When using more precise waveforms, extending the frequency range down to 20 Hz, removing a prominent blip glitch and marginalizing over all the model parameters, we find that the network signal-to-noise ratio posterior distribution lies mostly below the search value, with the 90\% confidence interval being $7.94^{+0.70}_{-1.05}$. If one assumes that the signal comes from a real gravitational-wave merger event, we find a light component $m_2 = 0.76^{+0.50}_{-0.14} M_\odot$, suggesting a compact object of mass below one solar mass at $83.8\%$ confidence level. Such a low mass for a compact object would suggest an unexpectedly light neutron star or a black hole of primordial or exotic origin. The primary mass would be $m_1 = 4.71^{+1.57}_{-2.18} M_\odot$, likely in the lower mass gap, for a mass ratio of $q=0.16^{+0.34}_{-0.06}$, at a distance of $D_{\rm L}=124^{+82}_{-48}$ Mpc. The improved sensitivity of the next observing runs would make it possible to observe similar signals with a higher SNR and to distinguish a sub-solar mass component.

The problem of the speed of the objects inside the Schwarzschild black hole is considered. The general result is that the value of the relative speed of the objects following their non-zero angular momentum trajectories, both of geodesic and non-geodesic character, when approaching the ultimate singularity, tends to the value of speed of light. There is only one exception when both objects move in the same plane and have parallel angular momenta. This outcome appears to have a deeper sense: it reflects the anisotropic character of the dynamics of interior of this particular black hole. The result in question means that near the singularity, collisions of two particles lead to an indefinitely large energy in the center of mass frame. Aforementioned properties have their counterpart in the phenomenon of an indefinitely large blueshift near the singularity.

We compute the power spectrum of super-horizon curvature perturbations generated during a late period of thermal inflation, taking into account fluctuation-dissipation effects resulting from the scalar flaton field's interactions with the ambient radiation bath. We find that, at the onset of thermal inflation, the flaton field may reach an equilibrium with the radiation bath even for relatively small coupling constants, maintaining a spectrum of thermal fluctuations until the critical temperature $T_c$, below which thermal effects stop holding the field at the false potential minimum. This enhances the field variance compared to purely quantum fluctuations, therefore increasing the average energy density during thermal inflation and damping the induced curvature perturbations. In particular, we find that this inhibits the later formation of primordial black holes, at least on scales that leave the horizon for $T>T_c$. The larger thermal field variance also reduces the duration of a period of fast-roll inflation below $T_c$, as the field rolls to the true potential minimum, which should also affect the generation of (large) curvature perturbations on even smaller scales.

The gravitational memory effect and its electromagnetic (EM) analog are potential probes in the strong gravity regime. In the literature, this effect is derived for static observers at asymptotic infinity. While this is a physically consistent approach, it restricts the space-time geometries for which one can obtain the EM memory effect. To circumvent this, we evaluate the EM memory effect for comoving observers (defined by the 4-velocity $u_{\mu}$) in arbitrary curved space-times. Using the covariant approach, we split Maxwell's equations into two parts -- projected parallel to the 4-velocity $u_{\mu}$ and into the 3-space orthogonal to $u_{\mu}$. Further splitting the equations into $1+1+2$-form, we obtain \emph{master equation} for the EM memory in an arbitrary curved space-time. We provide a geometrical understanding of the contributions to the memory effect. We then obtain EM memory for specific space-time geometries and discuss the salient features.

Neutrino masses and mixings produce vacuum oscillations, an established quantum mechanical phenomenon. In matter, the Mikheev-Smirnov-Wolfenstein effect, due to neutrino interactions with the background particles, triggers resonant flavor modification. In dense environments, sizable neutrino-neutrino interactions, shock waves and turbulence impact the neutrino flavor content under a variety of phenomena. Theoretical approaches of neutrino propagation range from the mean-field approximation to the full quantum kinetic equations. Intriguing connections have been uncovered between weakly interacting dense neutrino gases and other many-body systems and domains, from condensed matter and nuclear physics to quantum computing. Besides the intrinsic theoretical interest, establishing how neutrinos change flavor contributes to answer the longstanding open questions of how massive stars explode and of the r-process sites. It is also important for future observations of core-collapse supernova neutrinos and of the diffuse supernova neutrino background that should be discovered in the foreseeable future.