Papers for Wednesday, Feb 15 2023

Europa's subsurface ocean is a potential candidate for life in the outer solar system. It is thought that plumes may exist which eject ocean material out into space, which may be detected by a spacecraft flyby. Previous work on the feasibility of these detections has assumed a collisionless model of the plume particles. New models of the plumes including particle collisions have shown that a shock can develop in the plume interior as rising particles collide with particles falling back to the moon's surface, limiting the plume's altitude. Results show that the region over Europa's surface within which plumes would be separable from the H$_2$O atmosphere by JUICE (the region of separability) is reduced by up to a half with the collisional model compared to the collisionless model. Putative plume sources which are on the border of the region of separability for the collisionless model cannot be separated from the atmosphere when the shock is considered for a mass flux case of 100kg/s. Increasing the flyby altitude by 100km such that the spacecraft passes above the shock canopy results in a reduction in region of separability by a third, whilst decreasing the flyby altitude by 100km increases the region of separability by the same amount. We recommend flybys pass through or as close to the shock as possible to sample the most high-density region. If the spacecraft flies close to the shock, the structure of the plume could be resolvable using the neutral mass spectrometer on JUICE, allowing us to test models of the plume physics and understand the underlying physics of Europa's plumes. As the altitude of the shock is uncertain and dependent on unpredictable plume parameters, we recommend flybys be lowered where possible to reduce the risk of passing above the shock and losing detection coverage, density and duration.

Context. The existence of magnetic fields in the circumgalactic medium (CGM) is largely unconstrained. Their detection is important as magnetic fields can have a significant impact on the evolution of the CGM and, in turn, the fields can serve as tracers for dynamical processes in the CGM. Aims. With Faraday rotation of polarised background sources, we aim to detect a possible excess of the rotation measure in the surrounding area of nearby galaxies. Methods. We use 2,461 residual rotation measures (RRMs) observed with the LOw Frequency ARray (LOFAR), where the foreground contribution from the Milky Way is subtracted. The RRMs are then studied around a subset of 183 nearby galaxies that was selected by apparent $B$-band magnitude. Results. We find that, in general, the RRMs show no significant excess for small impact parameters (i.e. the perpendicular distance to the line of sight). However, if we only consider galaxies at higher inclination angles and sight lines that pass close to the minor axis of the galaxies, we find significant excess at impact parameters of less than 100 kpc. The excess in |RRM| is 3.7 $\rm rad\,m^{-2}$ with an uncertainty between $\pm 0.9~\rm rad\,m^{-2}$ and $\pm 1.3~\rm rad\,m^{-2}$ depending on the statistical properties of the background (2.8$\sigma$-4.1$\sigma$). With electron densities of ~$10^{-4}~\rm cm^{-3}$ this suggests magnetic field strengths of a few tenths of a micro Gauss. Conclusions. Our results suggest a slow decrease of the magnetic field strength with distance from the galactic disc such as expected if the CGM is magnetised by galactic winds and outflows.

The mergers of neutron stars expel a heavy-element enriched fireball which can be observed as a kilonova. The kilonova's geometry is a key diagnostic of the merger and is dictated by the properties of ultra-dense matter and the energetics of the collapse to a black hole. Current hydrodynamical merger models typically show aspherical ejecta. Previously, Sr$^+$ was identified in the spectrum of the the only well-studied kilonova AT2017gfo, associated with the gravitational wave event GW170817. Here we combine the strong Sr$^+$ P Cygni absorption-emission spectral feature and the blackbody nature of kilonova spectrum, to determine that the kilonova is highly spherical at early epochs. Line shape analysis combined with the known inclination angle of the source also shows the same sphericity independently. We conclude that energy injection by radioactive decay is insufficient to make the ejecta spherical. A magnetar wind or jet from the black hole disk could inject enough energy to induce a more spherical distribution in the overall ejecta, however an additional process seems necessary to make the element distribution uniform

The 21-cm power spectrum of reionization is a promising probe for cosmology and fundamental physics. Exploiting this new observable, however, requires fast predictors capable of efficiently scanning the very large parameter space of cosmological and astrophysical uncertainties. In this paper, we introduce the halo model of reionization (HMreio), a new analytical tool that combines the halo model of the cosmic dawn with the excursion-set bubble model for reionization, assuming an empirical correction factor to deal with overlapping ionization bubbles. First, HMreio is validated against results from the well-known semi-numerical code 21cmFAST, showing a good overall agreement for wave-modes of $k\lesssim 1$ h/Mpc. Based on this result, we perform a Monte-Carlo Markov-Chain (MCMC) forecast analysis assuming mock data from 1000-hour observations with the low-frequency part of the Square Kilometre Array (SKA) observatory. We simultaneously vary the six standard cosmological parameters together with seven astrophysical nuisance parameters quantifying the abundance and spectral properties of sources. Depending on the assumed theory error, we find very competitive constraints on cosmological parameters. In particular, it will be possible to conclusively test current cosmological tensions related to the Hubble parameter ($H_0$-tension) and the matter clustering amplitude ($S_8$-tension). Furthermore, the sum of the neutrino masses can be strongly constrained, making it possible to determine the neutrino mass hierarchy at the $\sim 90$ percent confidence level. However, these goals can only be achieved if the current modelling uncertainties are substantially reduced to below $\sim 3$ percent.

We present a high sensitivity spectral line survey of CO(1-0), CO(2-1), CO(3-2) and [CI](1-0) in 40 local (ultra) luminous infrared galaxies ((U)LIRGs), all with previous Herschel OH119 $\mu$m observations. We use single-dish observations (PI and archival) conducted with APEX, complemented with ALMA and ACA data. We study the total emission and pay special attention to the extended low-surface brightness components. We find a tight correlation between low-J CO and [CI] line luminosities suggesting their emission arise from similar regions, at least when averaged over galactic scales. We estimate a median CO-to-H$_2$ conversion factor of $1.7\pm 0.5$ M$_{\odot}$ (K km s$^{-1}$ pc$^2)^{-1}$ for ULIRGs, using [CI] as an independent tracer. We derive median galaxy-integrated CO line ratios ($r_{21}$, $r_{31}$ and $r_{32}$), as well as $r_{CICO}$, significantly higher than normal star forming galaxies, confirming the exceptional molecular gas properties of ULIRGs. We find that $r_{21}$ and $r_{32}$ are poor tracers of CO excitation in ULIRGs, while $r_{31}$ shows a positive trend with $L_{IR}$ and SFR, and a negative trend with the H$_2$ gas depletion timescales ($\tau_{dep}$). When studying CO line ratios as a function of gas kinematics, we find a positive relation between $r_{21}$ and $\sigma_v$, which can be explained by CO opacity effects. We find that the linewidths of [CI] lines are ~10% narrower than CO lines, which may suggest that the low optical depth of [CI] can challenge its detection in diffuse, low-surface brightness outflows, and so its use as a tracer of CO-dark H$_2$ gas in these components. Finally, we find that higher $L_{AGN}$ are associated to longer $\tau_{dep}$, consistent with the hypothesis that AGN feedback may reduce the efficiency of star formation.

How elliptical galaxies form is a key question in observational cosmology. While the formation of massive ellipticals is strongly linked to mergers, the low mass (Mstar < 10^9.5 MSun) regime remains less well explored. In particular, studying elliptical populations when they are blue, and therefore rapidly building stellar mass, offers strong constraints on their formation. Here, we study 108 blue, low-mass ellipticals (which have a median stellar mass of 10^8.7 MSun) at z < 0.3 in the COSMOS field. Visual inspection of extremely deep optical HSC images indicates that less than 3 per cent of these systems have visible tidal features, a factor of 2 less than the incidence of tidal features in a control sample of galaxies with the same distribution of stellar mass and redshift. This suggests that the star formation activity in these objects is not driven by mergers or interactions but by secular gas accretion. We combine accurate physical parameters from the COSMOS2020 catalog, with measurements of local density and the locations of galaxies in the cosmic web, to show that our blue ellipticals reside in low-density environments, further away from nodes and large-scale filaments than other galaxies. At similar stellar masses and environments, blue ellipticals outnumber their normal (red) counterparts by a factor of 2. Thus, these systems are likely progenitors of not only normal ellipticals at similar stellar mass but, given their high star formation rates, also of ellipticals at higher stellar masses. Secular gas accretion, therefore, likely plays a significant (and possibly dominant) role in the stellar assembly of elliptical galaxies in the low mass regime.

Placing the architecture of the Solar System within the broader context of planetary architectures is one of the primary topics of interest within planetary science. Exoplanet discoveries have revealed a large range of system architectures, many of which differ substantially from the Solar System model. One particular feature of exoplanet demographics is the relative prevalence of super-Earth planets, for which the Solar System lacks a suitable analog, presenting a challenge to modeling their interiors and atmospheres. Here we present the results of a large suite of dynamical simulations that insert a hypothetical planet in the mass range 1-10 $M_\oplus$ within the semi-major axis range 2-4 AU, between the orbits of Mars and Jupiter. We show that, although the system dynamics remain largely unaffected when the additional planet is placed near 3 AU, Mercury experiences substantial instability when the additional planet lies in the range 3.1-4.0 AU, and perturbations to the Martian orbit primarily result when the additional planet lies in the range 2.0-2.7 AU. We further show that, although Jupiter and Saturn experience relatively small orbital perturbations, the angular momentum transferred to the ice giants can result in their ejection from the system at key resonance locations of the additional planet. We discuss the implications of these results for the architecture of the inner and outer solar system planets, and for exoplanetary systems.

The assembly history of the Milky Way (MW) is a rapidly evolving subject, with numerous small accretion events and at least one major merger proposed in the MW's history. Accreted alongside these dwarf galaxies are globular clusters (GCs), which act as spatially coherent remnants of these past events. Using high precision differential abundance measurements from our recently published study, we investigate the likelihood that the MW clusters NGC 362 and NGC 288 are galactic siblings, accreted as part of the Gaia-Sausage-Enceladus (GSE) merger. To do this, we compare the two GCs at the 0.01 dex level for 20+ elements for the first time. Strong similarities are found, with the two showing chemical similarity on the same order as those seen between the three LMC GCs, NGC 1786, NGC 2210 and NGC 2257. However, when comparing GC abundances directly to GSE stars, marked differences are observed. NGC 362 shows good agreement with GSE stars in the ratio of Eu to Mg and Si, as well as a clear dominance in the r- compared to the s-process, while NGC 288 exhibits only a slight r-process dominance. When fitting the two GC abundances with a GSE-like galactic chemical evolution model, NGC 362 shows agreement with both the model predictions and GSE abundance ratios (considering Si, Ni, Ba and Eu) at the same metallicity. This is not the case for NGC 288. We propose that the two are either not galactic siblings, or GSE was chemically inhomogeneous enough to birth two similar, but not identical clusters with distinct chemistry relative to constituent stars.

We present JWST/NIRSpec Integral Field Spectrograph rest-frame optical data of the compact $z=5.55$ galaxy GS_3073. Its prominent broad components in several hydrogen and helium lines (while absent in the forbidden lines), and the detection of a large equivalent width of He II $\lambda4686$, EW(He II) $\sim20$ Angstrom, unambiguously identify it as an active galactic nucleus (AGN). We measure a gas-phase metallicity of $Z_{\rm gas}/Z_\odot\sim0.21^{+0.08}_{-0.04}$, lower than what has been inferred for both more luminous AGN at similar redshift and lower redshift AGN. We empirically show that classical emission line ratio diagnostic diagrams cannot be used to distinguish between the primary ionisation source (AGN or star formation) for such low-metallicity systems, whereas different diagnostic diagrams involving He II$\lambda4686$ prove very useful, independent of metallicity. We measure the central black hole mass to be $\log(M_{\rm BH}/M_\odot)\sim8.20^{+0.11}_{-0.16}$. While this places GS_3073 at the lower end of known high-redshift black hole masses, it still appears to be over-massive compared to its host galaxy properties. We detect an outflow with projected velocity $\gtrsim700$~km/s and an ionised gas mass outflow rate of about $100\ M_\odot/$yr, suggesting that GS_3073 is able to enrich the intergalactic medium with metals one billion years after the Big Bang.

A ray-theoretic phase space description of linear waves in a two-fluid (charges and neutrals) magnetized plasma is used to calculate analytic decay rates and mode transmission and conversion coefficients between fast and slow waves in two dimensions due to finite ion-neutral collision frequencies at arbitrary ionization fraction. This is relevant to partially ionized astrophysical plasmas, in particular solar and stellar atmospheres. The most important parameter governing collisional effects is the ratio of the wave frequency to the neutral-charges collision frequency, $\epsilon=\omega/\nu_{nc}$, with secondary dependence on ionization fraction and wave attack angle. Comparison is made to the one-fluid magnetohydrodynamic (MHD) case, and it is found that acoustic-to-acoustic and magnetic-to-magnetic transmission through the Alfv\'en-acoustic equipartition layer is decreased by a term of $\mathcal{O}(\epsilon^2)$ relative to one-fluid (infinite collision frequency), and correspondingly acoustic-to-magnetic and magnetic-to-acoustic conversion is increased. The neutral acoustic mode is shown to dissipate rapidly as $\nu_{nc}\to\infty$. Away from the mode conversion region, dissipative decay along the remaining magneto-acoustic rays scales as $\mathcal{O}(\epsilon)$ and is found to be much more effective on magnetically dominated rays compared to acoustically dominated rays. This produces a steep jump in dissipation in mode conversion regions, where the rays change character, and can produce localized heating there and beyond. Applications to the solar chromosphere are discussed.

We perform population synthesis of massive binaries to study the mergers of neutron stars (NSs) and black holes (BHs) with the cores of their giant secondaries during common envelope evolution (CEE). We use different values of the efficiency parameter $\alpha_{\rm CE}$ in the framework of the energy formalism for traditional CEE ($\alpha_{\rm CE} \leq 1$) and including additional energy sources to unbind the envelope ($\alpha_{\rm CE} > 1$). We constrain the possible values of $\alpha_{\rm CE}$ by comparing the results of our simulations with local rate densities of binary compact object mergers as inferred from gravitational waves observations. We find two primary evolutionary pathways of binary systems that result in NS-core mergers, while only one of them can also lead to the merger of a BH with the core of the giant star. We explore the zero age main sequence (ZAMS) statistical properties of systems that result in NS/BH-core mergers and find that the two evolutionary channels correspond to a bimodal distribution of orbital separations. We estimate the percentage of the mergers' event rates relative to core collapse supernovae (CCSNe). We include the effect of mass accreted by the NS/BH during CEE in a separate set of simulations and find it does not affect the mergers' event rates.

X-ray polarimetry is a unique way to probe geometrical configuration of highly-magnetized accreting neutron stars (X-ray pulsars). GRO J1008$-$57 is the first transient X-ray pulsar observed at two different flux levels by the Imaging X-ray Polarimetry Explorer (IXPE) during its outburst in November 2022. The polarization properties were found to be independent of the source luminosity, with the polarization degree varying between non-detection to about 15% over the pulse phase. Fitting the phase-resolved spectro-polarimetric data with the rotating vector model allowed us to estimate the pulsar inclination (130 deg, which is in good agreement with the orbital inclination), the position angle (75 deg) of the pulsar spin axis, and the magnetic obliquity (74 deg). This makes GRO J1008$-$57 the first confidently identified X-ray pulsar as a nearly orthogonal rotator. The results are discussed in the context of the neutron star atmosphere models and theories of pulsars' axis alignment.

Superluminous supernovae (SLSNe) are a rare class of stellar explosions with luminosities ~10-100 times greater than ordinary core-collapse supernovae. One popular model to explain the enhanced optical output of hydrogen-poor (Type I) SLSNe invokes energy injection from a rapidly spinning magnetar. A prediction in this case is that high-energy gamma rays, generated in the wind nebula of the magnetar, could escape through the expanding supernova ejecta at late times (months or more after optical peak). This paper presents a search for gamma-ray emission in the broad energy band from 100 MeV to 30 TeV from two Type I SLSNe, SN2015bn, and SN2017egm, using observations from Fermi-LAT and VERITAS. Although no gamma-ray emission was detected from either source, the derived upper limits approach the putative magnetar's spin-down luminosity. Prospects are explored for detecting very-high-energy (VHE; 100 GeV - 100 TeV) emission from SLSNe-I with existing and planned facilities such as VERITAS and CTA.

The observed properties of galaxies are strongly dependent on both their total stellar mass and their morphology. Furthermore, the environment is known to play a strong role in shaping them. The galaxy population in the local universe that is located in virialized clusters is found to be red, poorly star-forming, and mostly composed of early morphological types. Towards a holistic understanding of the mechanisms that drive galaxy evolution, we exploit the spectrophotometric data from the WINGS and OmegaWINGS local galaxy cluster surveys, and study the role of both the local and the large-scale environments. We attempt to disentangle their effects from the intrinsic characteristics of the galaxies, in shaping the star formation activity at fixed morphological type and stellar mass. Using a sample of field galaxies from the same surveys for comparison, we analyse the effects of the environment, embodied by the local density, clustercentric distance, and close neighbours, respectively, on the star formation histories of cluster galaxies. We find that local effects have a more relevant impact on galaxy stellar properties than the large-scale environment, and that morphology needs to be taken into account to pinpoint the mechanisms that are driving the influence of clusters in galaxy evolution.

Nearly one million light curves from the TESS Year 1 southern hemisphere extracted from Full Frame Images with the DIAmante pipeline are processed through the AutoRegressive Planet Search statistical procedure. ARIMA models remove trends and lingering autocorrelated noise, the Transit Comb Filter identifies the strongest periodic signal in the light curve, and a Random Forest machine learning classifier is trained and applied to identify the best potential candidates. Classifier training sets include injections of both planetary transit signals and contaminating eclipsing binaries. The optimized classifier has a True Positive Rate of 92.8% and a False Positive Rate of 0.37% from the labeled training set. The result of this DIAmante TESS autoregressive planet search (DTARPS) analysis is a list of 7,377 potential exoplanet candidates. The classifier has a False Positive Rate of 0.3%, a 64% recall rate for previously confirmed exoplanets, and a 78% negative recall rate for known False Positives. The completeness map of the injected planetary signals shows high recall rates for planets with 8 - 30 R(Earth) radii and periods 0.6-13 days and poor completeness for planets with radii < 2 R(Earth) or periods < 1 day. The list has many False Alarms and False Positives that need to be culled with multifaceted vetting operations (Paper II).

We explore local and global dynamical differences between the kinematics of ionised gas and stars in a sample of galaxies from Data Release 3 of the SAMI Galaxy Survey. We find better agreement between local (i.e., comparing on a spaxel-to-spaxel basis) velocities and dispersion of gas and stars in younger systems as with previous work on the asymmetric drift in galaxies, suggesting that the dynamics of stars and ionised gas are initially coupled. The intrinsic scatter around the velocity and dispersion relations increases with increasing stellar age and mass, suggesting that subsequent mechanisms such as internal processes, divergent star formation and assembly histories also play a role in setting and altering the dynamics of galaxies. The global (flux-weighted) dynamical support of older galaxies is hotter than in younger systems. We find that the ionised gas in galaxies is almost always dynamically colder than the stars with a steeper velocity gradient. In absolute terms, the local difference in velocity dispersion is more pronounced than the local difference in velocity, possibly reflecting inherent differences in the impact of turbulence, inflow and/or feedback on gas compared to stars. We suggest how these findings may be taken into account when comparing high and low redshift galaxy samples to infer dynamical evolution.

The DIAmante TESS AutoRegressive Planet Search (DTARPS) project seeks to identify photometric transiting planets from 976,814 southern hemisphere stars observed in Year 1 of the TESS mission. This paper follows the methodology developed by Melton et al. (Paper I) using light curves extracted and pre-processed by the DIAmante project (Montalto et al. 2020). Paper I emerged with a list of 7,377 light curves with statistical properties characteristic of transiting planets but dominated by False Alarms and False Positives. Here a multistage vetting procedure is applied including: centroid motion and crowding metrics, False Alarm and False Positive reduction, photometric binary elimination, and ephemeris match removal. The vetting produces a catalog of 462 DTARPS Candidates across the southern ecliptic hemisphere and 310 objects in a spatially incomplete Galactic Plane list. Fifty-eight percent were not previously identified as transiting systems. Candidates are flagged for possible blending from nearby stars based on Zwicky Transient Facility data and for possible radial velocity variations based on Gaia satellite data. Orbital periods and planetary radii are refined using astrophysical modeling; the resulting parameters closely match published values for Confirmed Planets. Their properties are discussed in Paper III.

In this work we study the potentialities of machine learning models in reconstructing the solar wind speed observations gathered in the first Lagrangian point by the ACE satellite in 2016--2017 using as input data galactic cosmic-ray flux variations measured with particle detectors hosted onboard the LISA Pathfinder mission also orbiting around L1 during the same years. We show that ensemble models composed of heterogeneous weak regressors are able to outperform weak regressors in terms of predictive accuracy. Machine learning and other powerful predictive algorithms open a window on the possibility of substituting dedicated instrumentation with software models acting as surrogates for diagnostics of space missions such as LISA and space weather science.

An algorithm is presented for solving the cold-start problem using observations of X-ray pulsars. Using a norm-minimization-based approach, the algorithm extends Lohan's banded-error intersection model to 3-dimensional space while reducing compute time by an order of magnitude. Higher-fidelity X-ray pulsar signal models, including the parallax effect, Shapiro delay, time dilation, and higher-order pulsar timing models, are considered. The feasibility of solving the cold-start problem using X-ray pulsar navigation is revisited with the improved models and prior knowledge requirements are discussed. Monte Carlo simulations are used to establish upper bounds on uncertainty and determine the accuracy of the algorithm. Results indicate that it is necessary to account for the parallax effect, time dilation, and higher-order pulsar timing models in order to successfully determine the position of the spacecraft in a cold-start scenario. The algorithm can uniquely identify a candidate spacecraft position within a 10 AU $\times$ 10 AU $\times$ 0.01 AU spheroid domain by observing eight to nine pulsars. The median position error of the algorithm is on the order of 15 km. Prior knowledge of spacecraft position is technically required, but only to an accuracy of 100 AU, making it practically unnecessary for navigation within the Solar System. Results further indicate that choosing lower-frequency pulsars increases the maximum domain size but also increases position error.

The DIAmante TESS AutoRegressive Planet Search (DTARPS) project, using novel statistical methods, has identified several hundred candidates for transiting planetary systems obtained from 0.9 million Full Frame Image light curves obtained in the TESS Year 1 southern hemisphere survey (Melton et al. 2022a and 2022b). Several lines of evidence, including limited reconnaissance spectroscopy, indicate that at least half are true planets rather than False Positives. Here various population properties of these objects are examined. Half of the DTARPS candidates are hot Neptunes, populating the 'Neptune desert' found in Kepler planet samples. The DTARPS samples also identify dozens of Ultra Short Period planets with orbital periods down to 5 hours, high priority systems for atmospheric transimssion spectroscopy, and planets orbiting low-mass M stars. DTARPS methodology is sufficiently well-characterized at each step that preliminary planet occurrence rates can be estimated. Except for the increase in hot Neptunes, DTARPS planet occurrence rates are consistent with Kepler rates. Overall, DTARPS provides one of the largest and most reliable catalog of TESS exoplanet candidates that can be tapped to improve our understanding of various exoplanetary populations and astrophysical processes.

We investigate the cosmological implications of an effective gravitational action, inspired by Sakharov's idea of induced gravity, containing non-local contributions from the operator $\left(\Box +\beta \right)^{-1} R$. The $\beta$ term is a novel feature in the panorama of non-local models of gravity, and arises naturally within Sakharov theory from the potential of a non-minimally coupled scalar field after a spontaneous symmetry breaking takes place. In this class of models the non-local contribution can acquire an oscillatory behaviour, thereby avoiding the divergence of the Hubble parameter at late times, which is another characteristic feature of the non-local models treated in the literature. Furthermore, the effective gravitational coupling $G_{\rm eff}$ inherits the oscillatory behaviour, resulting in alternating epochs of stronger and weaker gravity. Interestingly, this framework can describe a higher value for the Hubble factor today whilst potentially suppressing the growth of structures during the matter dominated epoch, providing a tool to catch \textit{two birds with one stone}, i.e. addressing at the same time both the $H_0$ and the $\sigma_8$ tensions.

We investigate the thermal condensation caused by a massive object that passes through the interstellar medium with high velocity, and propose a mechanism for creating a filamentary gaseous object, or interstellar contrail. Our main result shows that a long interstellar contrail can form with a certain parameter; a compact object more massive than $10^4\ {\rm M_\odot}$ can make a filament whose length is larger than $100\ {\rm pc}$. Observation of interstellar contrails may provide information on the number, masses, and velocities of fast-moving massive objects, and can be a new method for probing invisible gravitating sources such as intermediate-mass black holes.

Understanding the acceleration of Ultra High Energy Cosmic Rays is one of the great challenges of contemporary astrophysics. In this short review, we summarize the general observational constraints on their composition, spectrum and isotropy which indicate that nuclei heavier than single protons dominate their spectra above $\sim 5\,{\rm EeV}$, that they are strongly suppressed above energies $\sim50\,{\rm EeV}$, and that the only significant departure from isotropy is a dipole. Constraints based upon photopion and photodisintegration losses allow their ranges and luminosity density to be estimated. Three general classes of source model are discussed - magnetospheric models (including neutron stars and black holes), jet models (including Gamma Ray Bursts, Active Galactic Nuclei and Tidal Disruption Events) and Diffusive Shock Acceleration models (involving large accretion shocks around rich clusters of galaxies). The value of constructing larger and more capable arrays to measure individual masses at the highest energies and probably identifying their sources is emphasized.

Fast radio bursts (FRBs) are brief, luminous pulses with unknown physical origin. The repetition pattern of FRBs contains essential information about their physical nature and emission mechanisms. Using the two largest samples of FRB 20121102A and FRB 20201124, we report that the sources of the two FRBs reveal memory over a large range of timescales, from a few minutes to about an hour. This suggests that bright bursts are very likely to be followed by bursts of similar or larger magnitude, allowing for the prediction of intense bursts. The memory is detected from the coherent growths in burst-rate structures and the Hurst exponent. The waiting time distribution displays a power-law tail, which is consistent with a Poisson model with a time-varying rate. From cellular automaton simulations, we find that these characteristics can be well understood within the physical framework of a self-organized criticality system driven in a correlation way, such as random walk functions. These properties indicate that the triggers of bursts are correlated, preferring the crustal failure mechanism of neutron stars.

We report a new test of modified gravity theories using the large-scale structure of the Universe. This paper is the first attempt to (1) apply a joint analysis of the anisotropic components of galaxy two- and three-point correlation functions (2 and 3PCFs) to actual galaxy data and (2) constrain the nonlinear effects of degenerate higher-order scalar-tensor (DHOST) theories on cosmological scales. Applying this analysis to the Baryon Oscillation Spectroscopic Survey (BOSS) data release 12, we obtain the lower bounds of $-1.655 < \xi_{\rm t}$ and $-0.504 < \xi_{\rm s}$ at the $95\%$ confidence level on the parameters characterising the time evolution of the tidal and shift terms of the second-order velocity field. These constraints are consistent with GR predictions of $\xi_{\rm t}=15/1144$ and $\xi_{\rm s}=0$. Moreover, they represent a $35$-fold and $20$-fold improvement, respectively, over the joint analysis with only the isotropic 3PCF. We ensure the validity of our results by investigating various quantities, including theoretical models of the 3PCF, window function corrections, cumulative ${\rm S/N}$, Fisher matrices, and statistical scattering effects of mock simulation data. We also find statistically significant discrepancies between the BOSS data and the Patchy mocks for the 3PCF measurement. Finally, we package all of our 3PCF analysis codes under the name \textsc{HITOMI} and make them publicly available so that readers can reproduce all the results of this paper and easily apply them to ongoing future galaxy surveys.

We propose a non-exotic electromagnetic solution to the cosmological $^7$Li problem based upon a 2 MeV photo-emission line from color superconducting quark nuggets (CSCQN) that destroys Big-Bang nucleosynthesis (BBN) $^7$Be without affecting other abundances or Cosmic Microwave Background (CMB) physics. Conversion of CSCQN gluonic (rest-mass) energy to 2 MeV photons in the radiation-dominated post-BBN epoch reduces the primordial abundance of $^7$Be (and thus cosmological $^7$Li) by 2/3, provided the combined mass of the nuggets is greater than the total baryonic mass in the Universe. CSCQNs in our model are colorless, charge neutral, optically thin and decouple from the strong interaction (i.e. have minimal interaction with baryons and stars) making them a viable cold dark matter (CDM) candidate. The drainage (i.e. quantum tunnelling) of CSCQN Quantum-ChromoDynamics (QCD) vacuum to the external space-time (trivial QCD) vacuum offers a natural explanation of dark energy in our model and allows for a cosmology which evolves into a Lambda-CDM universe at low redshift with a possible resolution of the Hubble tension. The connection between CDM and DE in our model supports the notion that in-hadron confinement contains QCD condensates and disagrees with the conventional view of (space-filling) QCD condensates.

We present a multi-instrument study of the two precursor brightenings prior to the M6.5 flare (SOL2015-06-22T18:23) in the NOAA Active Region 12371, with a focus on the temperature (T), electron number density (n), and emission measure (EM). The data used in this study were obtained from four instruments with a variety of wavelengths, i.e., the Solar Dynamics Observatory's Atmospheric Imaging Assembly (AIA), in six extreme ultraviolet (EUV) passbands; the Expanded Owens Valley Solar Array (EOVSA) in microwave (MW); the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) in hard X-rays (HXR); and the Geostationary Operational Environmental Satellite (GOES) in soft X-rays (SXR). We compare the temporal variations of T, n, and EM derived from the different data sets. Here are the key results. (1) GOES SXR and AIA EUV have almost identical EM variations (1.5-3x10^48 per cm^3) and very similar T variations, from 8 to 15 million Kelvin (MK). (2) Listed from highest to lowest, EOVSA MW provides the highest temperature variations (15-60 MK), followed by RHESSI HXR (10-24 MK), then GOES SXR and AIA EUV (8-15 MK). (3) The EM variation from the RHESSI HXR measurements is always less than the values from AIA EUV and GOES SXR by at most 20 times. The number density variation from EOVSA MW is greater than the value from AIA EUV by at most 100 times. The results quantitatively describe the differences in the thermal parameters at the precursor phase, as measured by different instruments operating at different wavelength regimes and for different emission mechanisms.

At low ambient pressure, temperature gradients in porous soil lead to a gas flow, called thermal creep. With this regard, Mars is a unique as the conditions for thermal creep to occur in natural soil only exist on this planet in the solar system. Known as Knudsen compressor, thermal creep induces pressure variations. In the case of Mars, there might be a pressure maximum below the very top dust particle layers of the soil, which would support particle lift and might decrease threshold wind velocities necessary to trigger saltation or reduce angles of repose on certain slopes. In laboratory experiments, we applied diffusing wave spectroscopy (DWS) to trace minute motions of grains on the nm-scale in an illuminated simulated soil. This way, DWS visualizes pressure variations. We observe a minimum of motion which we attribute to the pressure maximum ~ 2 mm below the surface. The motion above but especially below that depth characteristically depends on the ambient pressure with a peak at an ambient pressure of about 3 mbar for our sample. This is consistent with earlier work on ejection of particle layers and is in agreement to a thermal creep origin. It underlines the supporting nature of thermal creep for particle lift which might be especially important on Mars.

A perturbation on the background inflaton potential can lead inflation into the ultraslow-roll stage and can thus remarkably enhance the power spectrum ${\cal P}_{\cal R}(k)$ of the primordial curvature perturbation on small scales. Such an enhanced ${\cal P}_{\cal R}(k)$ will result in primordial black holes (PBHs), contributing a significant fraction of dark matter, and will simultaneously generate sizable scalar-induced gravitational waves (SIGWs) as a secondorder effect. In this work, we calculate the PBH abundances $f_{\rm PBH}(M)$ and SIGW spectra $\Omega_{\rm GW}(f)$ in peak theory. We obtain the PBHs with desirable abundances in one or two typical mass windows at $10^{-17}\, M_\odot$, $10^{-13}\, M_\odot$, and $30\, M_\odot$, respectively. At the same time, the relevant SIGWs are expected to be observed by the next-generation gravitational wave detectors, without spoiling the current constraint. Especially, the SIGW associated with the PBH of $30\, M_\odot$ can also interpret the potential isotropic stochastic gravitational wave background from the NANOGrav 12.5-year dataset.

Extraction of characteristics of the complex light curve of pulsating stars is essential for stellar physics. We investigate the complex network (natural and horizontal visibility graphs) properties of the \dsct\ stars light curves observed by \tess. We find that the average shortest path length of \dsct\ light curves is a linear function of the logarithm of the network sizes, indicating the small-world and non-random properties. The small-world property signifies the connection of significant peaks of the light curve with small nearing peaks and other important peaks. The lognormal behavior of nodes' degree is evidence for non-random processes for stellar pulsations. This may be rooted in the different mechanisms of stellar dynamics, such as rotation, excitation of modes, and magnetic activity. The PageRank and nodes' degree distributions of \dsct\ stars collect in two HADS and non-HADS groups. The lower clustering for HADS than non-HADS indicates a more straightforward light curve (containing one or two independent modes) than a more complex light curve (including various oscillation modes) that might be a signature of surface gravity as an indication of stellar evolution. We show that reducing the window size of a light curve to about 5\% of the original one based on the network ranking preserves most of the star modes information. In this case, we also observe that the amplitude of most natural modes amplifies compared to the noise background in the power spectrum. These results indicate the capability of the network approach for interpreting pulsating stars' light curves.

In 2017 the Konkoly Observatory in Budapest published its first call for application inviting university students to carry out financially supported supervised research work and observing duties. The initiative quickly became popular and, so far, the program has supported 37 students. Five years later, is now time to summarize the experience gathered from both the institute and the participants. Notable results include numerous OTDK (student project) prizes awarded, first papers published, and acceptances into MSc and PhD programs both domestically and abroad, thus laying the foundations for the careers of several students. Among the feedback we have received from the students is the need for a more complex mentoring program, over and above of the funded research opportunities. A survey we conducted among the students indicates that communal and educational events are in the greatest demand, probably also induced by the lockdown restrictions experienced in the last few years. Through such events the students would not only build their community and start professional collaborations, but also learn more about various aspects of academia. In light of these results, we review possible avenues to improve the program.

Infrared measurements of asteroids are crucial for the determination of physical and thermal properties of individual objects, and for the understanding of the small-body populations in the solar system as a whole. But standard radiometric methods can only be applied if the orbit of an object is known, hence its position at the time of the observation. We present MIRI observations of the outer-belt asteroid 10920 and an unknown object, detected in all 9 MIRI bands in close proximity to 10920. We developed a new method "STM-ORBIT" to interpret the multi-band measurements without knowing the object's true location. The method leads to a confirmation of radiometric size-albedo solution for 10920 and puts constraints on the asteroid's location and orbit in agreement with its true orbit. Groundbased lightcurve observations of 10920, combined with Gaia data, indicate a very elongated object (a/b >= 1.5), with a spin-pole at (l, b) = (178{\deg}, 81{\deg}), and a rotation period of 4.861191 h. A thermophysical study leads to a size of 14.5 - 16.5 km, a geometric albedo between 0.05 and 0.10, and a thermal inertia in the range 9 to 35 Jm-2s-0.5K-1. For the newly discovered MIRI object, the STM-ORBIT method revealed a size of 100-230 m. The new asteroid must be on a very low-inclination orbit and it was located in the inner main-belt region during JWST observations. A beaming parameter {\eta} larger than 1.0 would push the size even below 100 meter, a main-belt regime which escaped IR detections so far. These kind of MIRI observations can therefore contribute to formation and evolution studies via classical size-frequency studies which are currently limited to objects larger than about one kilometer in size. We estimate that MIRI frames with pointings close to the ecliptic and only short integration times of a few seconds will always include a few asteroids, most of them will be unknown objects.

The Solar Mean Magnetic Field (SMMF) is the mean value of the line of sight (LOS) component of the solar vector magnetic field averaged over the visible hemisphere of the Sun. So far, the studies on SMMF have mostly been confined to the magnetic field measurements at the photosphere. In this study, we calculate and analyse the SMMF using magnetic field measurements at the chromosphere, in conjunction with that of photospheric measurements. For this purpose, we have used full disk LOS magnetograms derived from spectropolarimetric observations carried out in Fe I 630.15 nm and Ca II 854.2 nm by the Synoptic Optical Long term Investigations of the Sun (SOLIS)/Vector Spectromagnetograph (VSM) instrument during 2010 to 2017. It is found from this study that the SMMF at the chromosphere is weaker by a factor of 0.60 compared to the SMMF at the upper photosphere. The correlation analysis between them gives a Pearson correlation coefficient of 0.80. The similarity and reduced intensity of the chromospheric SMMF with respect to the photospheric SMMF corroborate the idea that it is the source of the Interplanetary Magnetic Field (IMF).

The abundance, temperature, and clustering of metals in the intergalactic medium are important parameters for understanding their cosmic evolution and quantifying their impact on cosmological analysis with the Ly $\alpha$ forest. The properties of these systems are typically measured from individual quasar spectra redward of the quasar's Ly $\alpha$ emission line, yet that approach may provide biased results due to selection effects. We present an alternative approach to measure these properties in an unbiased manner with the two-point statistics commonly employed to quantify large-scale structure. Our model treats the observed flux of a large sample of quasar spectra as a continuous field and describes the one-dimensional, two-point statistics of this field with three parameters per ion: the abundance (column density distribution), temperature (Doppler parameter) and clustering (cloud-cloud correlation function). We demonstrate this approach on multiple ions (e.g., C IV, Si IV, Mg II) with early data from the Dark Energy Spectroscopic Instrument (DESI) and high-resolution spectra from the literature. Our initial results show some evidence that the C IV abundance is higher than previous measurements and evidence for abundance evolution over time. The first full year of DESI observations will have over an order of magnitude more quasar spectra than this study. In a future paper we will use those data to measure the growth of clustering and its impact on the Ly $\alpha$ forest, as well as test other DESI analysis infrastructure such as the pipeline noise estimates and the resolution matrix.

We present a unified approach to (bi-)orthogonal basis sets for gravitating systems. Central to our discussion is the notion of mutual gravitational energy, which gives rise to the self-energy inner product on mass densities. We consider a first-order differential operator that is self-adjoint with respect to this inner product, and prove a general theorem that gives the conditions under which a (bi-)orthogonal basis set arises by repeated application of this differential operator. We then show that these conditions are fulfilled by all the families of analytical basis sets with infinite extent that have been discovered to date. The new theoretical framework turns out to be closely connected to Fourier-Mellin transforms, and it is a powerful tool for constructing general basis sets. We demonstrate this by deriving a basis set for the isochrone model and demonstrating its numerical reliability by reproducing a known result concerning unstable radial modes.

We discuss the gravitational wave spectrum produced by first-order phase transitions seeded by domain wall networks. This setup is important for many two-step phase transitions as seen for example in the singlet extension of the standard model. Whenever the correlation length of the domain wall network is larger than the typical bubble size, this setup leads to a gravitational wave signal that is shifted to lower frequencies and with an enhanced amplitude compared to homogeneous phase transitions without domain walls. We discuss our results in light of the recent PTA hints for gravitational waves.

(Abridged) In this article, we discuss the state of ``AGN feedback'' in radio-quiet (RQ) AGN. This study involves heterogeneous samples of nearby Seyfert and LINER galaxies as well as QSOs that have been observed at low radio frequencies (few ~100 MHz) with the GMRT and ~GHz frequencies with the VLA and VLBA. These multi-frequency, multi-resolution observations detect a range of arcsecond-scale radio spectral indices that are consistent with the presence of multiple contributors including starburst winds and AGN jets or winds; steep spectrum ``relic'' emission is observed as well. Polarization-sensitive data from the VLA and GMRT suggest that the radio outflows are stratified (e.g., in IIIZw2, Mrk231); distinct polarization signatures suggest that there could either be a ``spine + sheath'' structure in the radio outflow, or there could be a ``jet + wind'' structure. Similar nested biconical outflows can also explain the VLBA and SDSS emission-line data in the KISSR sample of double-peaked emission-line Seyfert and LINER galaxies. Furthermore, the modeling of the emission-lines with plasma modeling codes such as MAPPINGS indicates that parsec-scale jets and winds in these sources can disturb or move the narrow-line region gas clouds via the ``shock + precursor'' mechanism. Apart from the presence of ``relic'' emission, several Seyfert and LINER galaxies show clear morphological signatures of episodic jet activity. In one such source, NGC2639, at least four distinct episodes of jets are observed, the largest one of which was only detectable at 735 MHz with the GMRT. Additionally, a ~6 kpc hole in the CO molecular gas along with a dearth of young stars in the center of its host galaxy is observed. This suggests a link between episodic jet activity in RQ AGN and ``AGN feedback'' influencing the evolution of their host galaxies.

We briefly discuss a possible cosmological implication of the observed binary black hole mergings detected by LIGO-Virgo-Kagra collaboration (GWTC-3 catalogue) for the primordial black hole (PBH) formation in the early Universe. We show that the bumpy chirp mass distribution of the LVK BH+BH binaries can be fit with two distinct and almost equal populations: (1) astrophysical mergings from BH+BH formed in the modern Universe from evolution of massive binaries and (2) mergings of binary PBHs with initial log-normal mass distribution. We find that the PBH central mass ($M_c\simeq 30 M_\odot$) and distribution width derived from the observed LVK chirp masses are almost insensitive to the assumed double PBH formation model. To comply with the observed LVK BH+BH merging rate, the CDM PBH mass fraction should be $f_{pbh}\sim 10^{-3}$ but can be higher if PBH clustering is taken into account.

Precise measurements of stellar parameters are required in order to develop our theoretical understanding of stellar structure. These measurements enable errors and uncertainties to be quantified in theoretical models and constrain the physical interpretation of observed phenomena, such as the inflated radii of low-mass stars. We use newly-available TESS light curves combined with published radial velocity measurements to improve the characterization of 12 low mass eclipsing binaries composed of an M~dwarf accompanied by a brighter F/G star. We present and analyse ground-based simultaneous four-colour photometry for two targets. Our results include the first measurements of the fundamental properties of two of the systems. Light curve and radial velocity information were converted into the physical parameters of each component of the systems using an isochrone fitting method. We also derive the effective temperatures of the M~dwarfs, almost tripling the number of such measurements. The results are discussed in the context of radius inflation. We find that exquisite precision in the age estimation of young objects is required to determine their inflation status. However, all but three of the objects are securely located among the main sequence, demonstrating radius inflation and the necessity to develop our understanding of the complex physical processes governing the evolution of low-mass stars. We investigated the hypothesis that luminosity is unaffected by the inflation problem but the findings were not conclusive.

In November 2020, a new, bright object, eRASSt J234402.9$-$352640, was discovered in the second all-sky survey of SRG/eROSITA. The object brightened by a factor of at least 150 in 0.2--2.0 keV flux compared to an upper limit found six months previous, reaching an observed peak of $1.76_{-0.24}^{+0.03} \times 10^{-11}$ erg cm$^{-2}$ s$^{-1}$. The X-ray ignition is associated with a galaxy at $z=0.10$, making the peak luminosity log$_{10}(L_{\rm 0.2-2keV}/[\textrm{erg s}^{-1}])$=$44.7\pm0.1$. Around the time of the rise in X-ray flux, the nucleus of the galaxy brightened by approximately 3 mag. in optical photometry, after correcting for the host. We present data from Swift, XMM-Newton, and NICER, which reveal a very soft spectrum as well as strong 0.2--2.0 keV flux variability on multiple timescales. Optical spectra taken in the weeks after the ignition event show a blue continuum with broad, asymmetric Balmer emission lines, and high-ionisation ([OIII]$\lambda\lambda$4959,5007) and low-ionisation ([NII]$\lambda$6585, [SII]$\lambda\lambda$6716,6731) narrow emission lines. Following the peak in the optical light curve, the X-ray, UV, and optical photometry all show a rapid decline. The X-ray light curve shows a decrease in luminosity of $\sim$0.45 over 33 days and the UV shows a drop of $\sim$0.35. eRASSt J234402.9$-$352640 also shows a brightening in the mid-infrared, likely powered by a dust echo of the luminous ignition. We find no evidence in Fermi-LAT $\gamma$-ray data for jet-like emission. The event displays characteristics of a tidal disruption event (TDE) as well as of an active galactic nucleus (AGN), complicating its classification. Based on the softness of the X-ray spectrum, the presence of high-ionisation optical emission lines, and the likely infrared echo, we find that a TDE within a turned-off AGN best matches our observations.

Intermediate-Luminosity Red Transients (ILRTs) are a class of observed transient posited to arise from the production of an electron-capture supernova from a super-asymptotic giant branch star within a dusty cocoon. In this paper, we present a systematic analysis of narrow Na I D absorption as a means of probing the circumstellar environment of these events. We find a wide diversity of evolution in ILRTs in terms of line strength, time-scale, and shape. We present a simple toy model designed to predict this evolution as arising from ejecta from a central supernova passing through a circumstellar environment wherein Na II is recombining to Na I over time. We find that while our toy model can qualitatively explain the evolution of a number of ILRTs, the majority of our sample undergoes evolution more complex than predicted. The success of using the Na I D doublet as a diagnostic tool for studying circumstellar material will rely on the availability of regular high-resolution spectral observations of multiple ILRTs, and more detailed spectral modelling will be required to produce models capable of explaining the diverse range of behaviours exhibited by ILRTs. In addition, the strength of the Na I D absorption feature has been used as a means of estimating the extinction of sources, and we suggest that the variability visible in ILRTs would prevent such methods from being used for this class of transient, and any others showing evidence of variability

The Gaia Data Release 3 (DR3), published in June 2022, delivers a diverse set of astrometric, photometric, and spectroscopic measurements for more than a billion stars. The wealth and complexity of the data makes traditional approaches for estimating stellar parameters for the full Gaia dataset almost prohibitive. We have explored different supervised learning methods for extracting basic stellar parameters as well as distances and line-of-sight extinctions, given spectro-photo-astrometric data (including also the new Gaia XP spectra). For training we use an enhanced high-quality dataset compiled from Gaia DR3 and ground-based spectroscopic survey data covering the whole sky and all Galactic components. We show that even with a simple neural-network architecture or tree-based algorithm (and in the absence of Gaia XP spectra), we succeed in predicting competitive results (compared to Bayesian isochrone fitting) down to faint magnitudes. We will present a new Gaia DR3 stellar-parameter catalogue obtained using the currently best-performing machine-learning algorithm for tabular data, XGBoost, in the near future.

It is expected that cosmic rays (CRs) escape from high-redshift galaxies at redshift $z\sim 10 \, - \, 20$ because CRs are accelerated by supernova remnants of the first stars. Although ultraviolet and X-ray photons are widely considered the main source of heating of the intergalactic medium, CRs can also contribute to it. When the CRs propagate in the intergalactic medium, in addition to the heating process due to CR ionization, resistive heating occurs due to the electron return current induced by the streaming CRs. We evaluate the heating rate around a galaxy as a function of the distance from the galaxy. We find that the resistive heating induced by CRs dominates over the other heating processes in the vicinity of the galaxy $r \lesssim 10^2 \, \mathrm{kpc}$ until the temperature reaches $T\sim 10^4 \, \mathrm{K}$. We also recalculate the strength of the magnetic field generated by streaming CRs under the presence of X-ray heating and show that achieved strength can be about $1$ order of magnitude smaller when the X-ray heating is included. The presence of the "first" CRs could be confirmed from the characteristic signature of CR heating imprinted on the $21$-$\mathrm{cm}$ line map in future radio observations.

Blind cleaning methods are currently the preferred strategy for handling foreground contamination in single-dish HI intensity mapping surveys. Despite the increasing sophistication of blind techniques, some signal loss will be inevitable across all scales. Constructing a corrective transfer function using mock signal injection into the contaminated data has been a practice relied on for HI intensity mapping experiments. However, assessing whether this approach is viable for future intensity mapping surveys where precision cosmology is the aim, remains unexplored. In this work, using simulations, we validate for the first time the use of a foreground transfer function to reconstruct power spectra of foreground-cleaned low-redshift intensity maps and look to expose any limitations. We reveal that even when aggressive foreground cleaning is required, which causes ${>}\,50\%$ negative bias on the largest scales, the power spectrum can be reconstructed using a transfer function to within sub-percent accuracy. We specifically outline the recipe for constructing an unbiased transfer function, highlighting the pitfalls if one deviates from this recipe, and also correctly identify how a transfer function should be applied in an auto-correlation power spectrum. We validate a method that utilises the transfer function variance for error estimation in foreground-cleaned power spectra. Finally, we demonstrate how incorrect fiducial parameter assumptions (up to ${\pm}100\%$ bias) in the generation of mocks, used in the construction of the transfer function, do not significantly bias signal reconstruction or parameter inference (inducing ${<}\,5\%$ bias in recovered values).

Extreme mass ratio inspirals (EMRIs) are among the primary targets for the Laser Interferometer Space Antenna (LISA). The extreme mass ratios of these systems result in relatively weak GW signals, that can be individually resolved only for cosmologically nearby sources (up to $z\approx2$). The incoherent piling up of the signal emitted by unresolved EMRIs generate a confusion noise, that can be formally treated as a stochastic GW background (GWB). In this paper, we estimate the level of this background considering a collection of astrophysically motivated EMRI models, spanning the range of uncertainties affecting EMRI formation. To this end, we employed the innovative Augmented Analytic Kludge waveforms and used the full LISA response function. For each model, we compute the GWB SNR and the number of resolvable sources. Compared to simplified computations of the EMRI signals from the literature, we find that for a given model the GWB SNR is lower by a factor of $\approx 2$ whereas the number of resolvable sources drops by a factor 3-to-5. Nonetheless, the vast majority of the models result in potentially detectable GWB which can also significantly contribute to the overall LISA noise budget in the 1-10 mHz frequency range.

The halo occupation distribution (HOD) framework is an empirical method to describe the connection between dark matter halos and galaxies, which is constrained by small scale clustering data. Efficient fitting procedures are required to scan the HOD parameter space. This paper describes such a method based on Gaussian Processes to iteratively build a surrogate model of the posterior of the likelihood surface from a reasonable amount of likelihood computations, typically two orders of magnitude less than standard Monte Carlo Markov chain algorithms. Errors in the likelihood computation due to stochastic HOD modelling are also accounted for in the method we propose. We report results of reproducibility, accuracy and stability tests of the method derived from simulation, taking as a test case star-forming emission line galaxies, which constitute the main tracer of the Dark Energy Spectroscopic Instrument and have so far a poorly constrained galaxy-halo connection from observational data.

Debris disks are the signposts of collisionally eroding planetesimal circumstellar belts, whose study can put important constraints on the structure of extrasolar planetary systems. The best constraints on the morphology of disks are often obtained from spatially resolved observations in scattered light. Here, we investigate the young (~16 Myr) bright gas-rich debris disk around HD121617. We use new scattered-light observations with VLT/SPHERE to characterize the morphology and the dust properties of this disk. From these properties we can then derive constraints on the physical and dynamical environment of this system, for which significant amounts of gas have been detected. The disk morphology is constrained by linear-polarimetric observations in the J band. Based on our modeling results and archival photometry, we also model the SED to put constraints on the total dust mass and the dust size distribution. We explore different scenarios that could explain these new constraints. We present the first resolved image in scattered light of the debris disk HD121617. We fit the morphology of the disk, finding a semi-major axis of 78.3$\pm$0.2 au, an inclination of 43.1$\pm$0.2{\deg} and a position angle of the major axis with respect to north, of 239.8$\pm$0.3{\deg}, compatible with the previous continuum and CO detection with ALMA. Our analysis shows that the disk has a very sharp inner edge, possibly sculpted by a yet-undetected planet or gas drag. While less sharp, its outer edge is steeper than expected for unperturbed disks, which could also be due to a planet or gas drag, but future observations probing the system farther from the main belt would help explore this further. The SED analysis leads to a dust mass of 0.21$\pm$0.02 M$_{\oplus}$ and a minimum grain size of 0.87$\pm$0.12 $\mu$m, smaller than the blowout size by radiation pressure, which is not unexpected for very bright col...

High-energy cosmic rays interact in the Earth's atmosphere and produce extensive air showers (EASs) which can be measured with large detector arrays at the ground. The interpretation of these measurements relies on models of the EAS development which represents a challenge as well as an opportunity to test quantum chromodynamics (QCD) under extreme conditions. The EAS development is driven by hadron-ion collisions under low momentum transfer in the non-perturbative regime of QCD. Under these conditions, hadron production cannot be described using first principles and these interactions cannot be probed with existing collider experiments. Thus, accurate measurements of the EAS development provide a unique probe of multi-particle production in hadronic interactions.

In 2018, images were released of a planet being formed around the star PDS 70, offering a tantalizing glimpse into how planets come into being. However, many questions remain about how dust evolves into planets, and astrophysical observations are unable to provide all the answers. It is therefore necessary to perform experiments to reveal key details and, to avoid unwanted effects from the Earth's gravitational pull, it is often necessary to perform such experiments in microgravity platforms. This Review sketches current models of planet formation and describes the experiments needed to test the models.

We report initial results from an ongoing spectroscopic survey of central stars of faint planetary nebulae (PNe), obtained with the Low-Resolution Spectrograph on the Hobby-Eberly Telescope. The six PN nuclei (PNNi) discussed here all have strong emission at the O VI 3811-3834 A doublet, indicative of very high temperatures. Five of them--the nuclei of Ou 2, Kn 61, Kn 15, Abell 72, and Kn 130--belong to the hydrogen-deficient PG 1159 class, showing a strong absorption feature of He II and C IV at 4650-4690 A. Based on exploratory comparisons with synthetic model-atmosphere spectra, and the presence of Ne VIII emission lines, we estimate them to have effective temperatures of order 170,000 K. The central star of Kn 15 has a Wolf-Rayet-like spectrum, with strong and broad emission lines of He II, C IV, N V, and O V-VI. We classify it [WO2], but we note that the N V 4604-4620 A emission doublet is extremely strong, indicating a relatively high nitrogen abundance. Several of the emission lines in Kn 15 vary in equivalent width by factors as large as 1.5 among our four observations from 2019 to 2022, implying significant variations in the stellar mass-loss rate. We encourage spectroscopic monitoring. Follow-up high-time-resolution photometry of these stars would be of interest, given the large fraction of pulsating variables seen among PG 1159 and [WO] PNNi.

We present JWST NIRCam 9-band near-infrared imaging of the luminous $z=10.6$ galaxy GN-z11 from the JWST Advanced Deep Extragalactic Survey (JADES) of the GOODS-N field. We find a spectral energy distribution (SED) entirely consistent with the expected form of a high-redshift galaxy: a clear blue continuum from 1.5 to 4 microns with a complete dropout in F115W. The core of GN-z11 is extremely compact in JWST imaging. We analyze the image with a two-component model, using a point source and a S\'{e}rsic profile that fits to a half-light radius of 200 pc and an index $n=0.9$. We find a low-surface brightness haze about $0.4''$ to the northeast of the galaxy, which is most likely a foreground object but might be a more extended component of GN-z11. At a spectroscopic redshift of 10.60 (Bunker et al. 2023), the comparison of the NIRCam F410M and F444W images spans the Balmer jump. From population synthesis modeling, here assuming no light from an active galactic nucleus, we reproduce the SED of GN-z11, finding a stellar mass of $\sim$$10^{9}~M_{\odot}$, a star-formation rate of $\sim$$20~M_{\odot}~\mathrm{yr}^{-1}$ and a young stellar age of $\sim$$20~\mathrm{Myr}$. As massive galaxies at high redshift are likely to be highly clustered, we search for faint neighbors of GN-z11, finding 9 galaxies out to $\sim$5 comoving Mpc transverse with photometric redshifts consistent with $z=10.6$, and a 10$^{\rm th}$ more tentative dropout only $3''$ away.

Stars are fossils that retain the history of their host galaxies. Elements heavier than helium are created inside stars and are ejected when they die. From the spatial distribution of elements in galaxies, it is therefore possible to constrain the physical processes during galaxy formation and evolution. This approach, Galac- tic archaeology, has been popularly used for our Milky Way Galaxy with a vast amount of data from Gaia satellite and multi-object spectrographs to understand the origins of sub-structures of the Milky Way. Thanks to integral field units, this ap- proach can also be applied to external galaxies from nearby to distant universe with the James Webb Space Telescope. In order to interpret these observational data, it is necessary to compare with theoretical predictions, namely chemodynamical simula- tions of galaxies, which include detailed chemical enrichment into hydrodynamical simulations from cosmological initial conditions. These simulations can predict the evolution of internal structures (e.g., metallicity radial gradients) as well as that of scaling relations (e.g., the mass-metallicity relations). After explaining the formula and assumptions, we will show some example results, and discuss future prospects.

We present JADES JWST/NIRSpec spectroscopy of GN-z11, the most luminous candidate $z>10$ Lyman break galaxy in the GOODS-North field with $M_{UV}=-21.5$. We derive a redshift of $z=10.603$ (lower than previous determinations) based on multiple emission lines in our low and medium resolution spectra over $0.8-5.3\,\mu$m. We significantly detect the continuum and measure a blue rest-UV spectral slope of $\beta=-2.4$. Remarkably, we see spatially-extended Lyman-$\alpha$ in emission (despite the highly-neutral IGM expected at this early epoch), offset 555 km/s redward of the systemic redshift. From our measurements of collisionally-excited lines of both low- and high-ionization (including [O II] $\lambda3727$, [Ne III] $\lambda 3869$ and C III] $\lambda1909$) we infer a high ionization parameter ($\log U\sim -2$). We detect the rarely-seen N IV] $\lambda1486$ and N III]$\lambda1748$ lines in both our low and medium resolution spectra, with other high ionization lines seen in low resolution spectrum such as He II (blended with O III]) and C IV (with a possible P-Cygni profile). Based on the observed rest-UV line ratios, we cannot conclusively rule out photoionization from AGN. The high C III]/He II ratios, however, suggest a likely star-formation explanation. If the observed emission lines are powered by star formation, then the strong N III] $\lambda1748$ observed may imply an unusually high $N/O$ abundance. Balmer emission lines (H$\gamma$, H$\delta$) are also detected, and if powered by star formation rather than an AGN we infer a star formation rate of $\sim 20-30 M_{\odot}\,\rm yr^{-1}$ (depending on the IMF) and low dust attenuation. Our NIRSpec spectroscopy confirms that GN-z11 is a remarkable galaxy with extreme properties seen 430 Myr after the Big Bang.

This work presents the results from extending the long-term monitoring program of stellar motions within the Galactic Center to include stars with separations of 2 - 7 arcseconds from the compact radio source, Sgr A*. In comparison to the well studied inner 2 arcsec, a longer baseline in time is required to study these stars. With 17 years of data, a sufficient number of positions along the orbits of these outer stars can now be measured. This was achieved by designing a source finder to track the positions of ~ 2000 stars in NACO/VLT adaptive-optics-assisted images of the Galactic Center from 2002 to 2019. Of the studied stars, 54 exhibit significant accelerations toward Sgr A*, most of which have separations of between 2 and 3 arcseconds from the black hole. A further 20 of these stars have measurable radial velocities from SINFONI/VLT stellar spectra, which allows for the calculation of the orbital elements for these stars, thus increasing the number of known orbits in the Galactic Center by ~ 40 %. With orbits, we can consider which structural features within the Galactic Center nuclear star cluster these stars belong to. Most of the stars have orbital solutions that are consistent with the known clockwise rotating disk feature. Further, by employing Monte Carlo sampling for stars without radial velocity measurements, we show that many stars have a subset of possible orbits that are consistent with one of the known disk features within the Galactic Center.

Colour-magnitude diagrams reveal a population of blue (hot) sub-luminous objects with respect to the main sequence. These hot sub-luminous stars are the result of evolutionary processes that require stars to expel their obscuring, hydrogen-rich envelopes to reveal the hot helium core. As such, these objects offer a direct window into the hearts of stars that are otherwise inaccessible to direct observation. We showcase MeerLICHT's capabilities of detecting faint hot subdwarfs and identifying the dominant frequency in the photometric variability of these compact hot stars, in comparison to their $Gaia$ DR3 data. We hunt for oscillations, which will be an essential ingredient for accurately probing stellar interiors in future asteroseismology. Comparative MeerLICHT and $Gaia$ colour-magnitude diagrams are presented as a way to select hot subdwarfs from our sample. A dedicated frequency determination technique is developed and applied to the selected candidates to determine their dominant variability using time-series data from MeerLICHT and $Gaia$ DR3. We explore the power of both datasets in determining the dominant frequency. Using the $g-i$ colour, MeerLICHT offers a colour-magnitude diagram that is comparable in quality to that of $Gaia$ DR3. The MeerLICHT colour-colour diagrams allow for the study of different stellar populations. The frequency analysis of MeerLICHT and $Gaia$ DR3 data demonstrates the superiority of our MeerLICHT multi-colour photometry in estimating the dominant frequency compared to the sparse $Gaia$ DR3 data. MeerLICHT's multi-band photometry leads to the discovery of high-frequency faint subdwarfs. Our MeerLICHT results are a proof-of-concept of the capacity of the BlackGEM instrument currently in the commissioning stage at ESO's La Silla Observatory in Chile.

According to the inflationary theory of cosmology, most elementary particles in the current universe were created during a period of reheating after inflation. In this work we self-consistently couple the Einstein-inflaton equations to a strongly coupled quantum field theory (QFT) as described by holography. We show that this leads to an inflating universe, a reheating phase and finally a universe dominated by the QFT in thermal equilibrium.

Scalar boson stars have attracted attention as simple models for exploring the nonlinear dynamics of a large class of ultra compact and black hole mimicking objects. Here, we study the impact of interactions in the scalar matter making up these stars. In particular, we show the pivotal role the scalar phase and vortex structure play during the late inspiral, merger, and post-merger oscillations of a binary boson star, as well as their impact on the properties of the merger remnant. To that end, we construct constraint satisfying binary boson star initial data and numerically evolve the nonlinear set of Einstein-Klein-Gordon equations. We demonstrate that the scalar interactions can significantly affect the inspiral gravitational wave amplitude and phase, and the length of a potential hypermassive phase shortly after merger. If a black hole is formed after merger, we find its spin angular momentum to be consistent with similar binary black hole and binary neutron star merger remnants. Furthermore, we formulate a mapping that approximately predicts the remnant properties of any given binary boson star merger. Guided by this mapping, we use numerical evolutions to explicitly demonstrate, for the first time, that rotating boson stars can form as remnants from the merger of two non-spinning boson stars. We characterize this new formation mechanism and discuss its robustness. Finally, we comment on the implications for rotating Proca stars.

A heavy meta-stable field dominates the energy density of the universe after inflation. The dissipation of this field continuously sources high-energy particles. In general, the dissipation rate of this meta-stable field can have a non-trivial time dependence. We study the impact of this time-dependent dissipation rate on the thermalization of the high-energy decay products of the meta-stable field. These energetic particles can contribute substantially to dark matter production in addition to the usual production from the thermal bath particles during reheating. We investigate the impact of this generalized dissipation on dark matter production in a model-independent way. We illustrate the parameter space that explains the observed dark matter relic abundance in various cosmological scenarios. We observed that dark matter having a mass larger than the maximum temperature attained by the thermal bath can be produced from the collision of the high-energy particles which are not yet thermalized.

It is intriguing to ask whether the existence of primordial black holes (PBHs) in the early universe could significantly reduce the abundance of certain stable massive particles (SMP) via gravitational capture, after which the PBHs evaporate before BBN to avoid conflict with stringent bounds. For example, this mechanism is relevant to an alternative solution of the monopole problem proposed by Stojkovic and Freese, in which magnetic monopoles produced in the early universe are captured by PBHs, thus freeing inflation from having to occur during or after the corresponding phase transitions that produced the monopoles. In this work, we reanalyze the solution by modelling the capture process in the same way as the coexisting monopole annihilation, which exhibits typical features of a diffusive capture. A monochromatic PBH mass function and a radiation-dominated era before PBH evaporation are assumed. We found that for Pati-Salam monopoles corresponding to a symmetry breaking scale between $10^{10}\,\text{GeV}$ and $10^{15}\,\text{GeV}$, the capture rate is many orders of magnitude below what is needed to cause a significant reduction of the monopole density. The difference with respect to previous literature can be attributed to both the modelling of the capture process and also the assumption on the PBH mass function. Within our assumptions, we also found that the magnetic charge that is large enough to make an extremal magnetic black hole cosmologically stable cannot be obtained from magnetic charge fluctuation via monopole capture. The large magnetic charged required by cosmological stability can nevertheless be obtained from magnetic charge fluctuation at PBH formation, and if later the monopole abundance can be reduced significantly by some non-inflationary mechanism, long-lived near-extremal magnetic black holes of observational relevance might result.

In this work, we analyse the late-time evolution of the universe for a particular $f(R)$ gravity model built from an exponential function of the scalar curvature. Following the literature, we write the field equations in terms of a suited statefinder function ($y_H(z)$) and considering well motivated physical initial conditions, the resulting equations are solved numerically. Also, the cosmological parameters $w_{\rm{DE}}$, $w_{\rm{eff}}$, $\Omega_{\rm{DE}}$ and $H(z)$ and the statefinder quantities $q$, $j$, $s$ and $Om(z)$ are explicitly expressed in terms of $y_H(z)$ and its derivatives. Furthermore, setting an appropriate set of values for the model parameters, the cosmological parameters as well as the statefinder quantities are plotted, and their present values (at $z=0$), are shown to be compatible with Planck 2018 observations and the $\Lambda$CDM-model values. Considering updated measurements from the dynamics of the expansion of the universe, $H(z)$, we perform an statistical analysis to constrain the free parameters of the model, finding a particular set of values that fit the data well and predict acceptable values for the cosmological and statefinder parameters at present time. Therefore, the $f(R)$ gravity model is found to be consistent with the considered observational data, and a viable alternative to explain the late-time acceleration of the universe.