Papers for Thursday, Apr 06 2023

We present the JWST Cycle 1 53.8hr medium program FRESCO, short for "First Reionization Epoch Spectroscopically Complete Observations". FRESCO covers 62 arcmin$^2$ in each of the two GOODS/CANDELS fields for a total area of 124 arcmin$^2$ exploiting JWST's powerful new grism spectroscopic capabilities at near-infrared wavelengths. By obtaining ~2 hr deep NIRCam/grism observations with the F444W filter, FRESCO yields unprecedented spectra at R~1600 covering 3.8 to 5.0 $\mu$m for most galaxies in the NIRCam field-of-view. This setup enables emission line measurements over most of cosmic history, from strong PAH lines at z~0.2-0.5, to Pa$\alpha$ and Pa$\beta$ at z~1-3, HeI and [SIII] at z~2.5-4.5, H$\alpha$ and [NII] at z~5-6.5, up to [OIII] and H$\beta$ for z~7-9 galaxies, and possibly even [OII] at z~10-12. FRESCO's grism observations provide total line fluxes for accurately estimating galaxy stellar masses and calibrating slit-loss corrections of NIRSpec/MSA spectra in the same field. Additionally, FRESCO results in a mosaic of F182M, F210M, and F444W imaging in the same fields to a depth of ~28.2 mag (5 $\sigma$ in 0.32" diameter apertures). Together with this publication, the v1 imaging mosaics are released as high-level science products via MAST. Here, we describe the overall survey design and the key science goals that can be addressed with FRESCO. We also highlight several, early science results, including: spectroscopic redshifts of Lyman break galaxies that were identified almost 20 years ago, the discovery of broad-line active galactic nuclei at z>4, and resolved Pa$\alpha$ maps of galaxies at z~1.4. These results demonstrate the enormous power for serendipitous discovery of NIRCam/grism observations. Given the wealth of ancillary data available in these fields, the zero-proprietary time FRESCO data is poised to enable a large amount of legacy science by the community.

We study the nonlinear effects of minimally coupled, massless, cosmological scalar fields on the cosmic microwave background (CMB). These fields can exhibit post-recombination parametric resonance and subsequent nonlinear evolution leading to novel contributions to the gravitational potential. We compute the resulting contributions to the CMB temperature anisotropies through the time-variation of the gravitational potential (i.e., the integrated Sachs-Wolfe (ISW) effect). We find that fields that constitute 5% of the total energy density and become dynamical at $z_c \simeq 10^{4}$ can produce marginally observable ISW signals at multipoles $\ell \simeq 2000$. Fields that become dynamical at earlier times and/or have initial displacements at a flatter part of their potential, produce ISW contributions that are significantly larger and at higher multipoles. We calculate these dynamics and the resulting evolution of gravitational perturbations using analytic estimates alongside detailed nonlinear lattice simulations, which couple scalar fields and cosmological fluids to a perturbed metric. Finally, we discuss the possibility of detecting these features with future high-resolution CMB observations.

Most stars, binaries, and higher multiplicity systems are thought to form in stellar clusters and associations, which later dissociate. Very wide binaries can be easily disrupted in clusters due to dynamical evaporation (soft binaries) and/or due to tidal disruption by the gravitational potential of the cluster. Nevertheless, wide binaries are quite frequent in the field, where they can sometimes play a key role in the formation of compact binaries, and serve as tools to study key physical processes. Here we use analytic tools to study the dynamical formation of soft binaries in clusters, and their survival as field binaries following cluster dispersion. We derive the expected properties of very wide binaries both in clusters and in the field. We analytically derive their detailed distributions, including wide-binary fraction as a function of mass in different cluster environments, binaries mass functions and mass ratios, and the distribution of their orbital properties. We show that our calculations agree well on most aspects with the results of N-body simulations, but show some different binary-fraction dependence on the cluster mass. We find that the overall fraction of wide binaries scales as $\propto N_\star^{-1}$ where $N_\star$ is the size of the cluster, even for non-equal mass stars. More massive stars are more likely to capture wide companions, with most stars above five solar mass likely to capture at least one stellar companion, and triples formation is found to be frequent.

Stellar streams are potentially a very sensitive observational probe of galactic astrophysics, as well as the dark matter population in the Milky Way. On the other hand, performing a detailed, high-fidelity statistical analysis of these objects is challenging for a number of key reasons. Firstly, the modelling of streams across their (potentially billions of years old) dynamical age is complex and computationally costly. Secondly, their detection and classification in large surveys such as Gaia renders a robust statistical description regarding e.g., the stellar membership probabilities, challenging. As a result, the majority of current analyses must resort to simplified models that use only subsets or summaries of the high quality data. In this work, we develop a new analysis framework that takes advantage of advances in simulation-based inference techniques to perform complete analysis on complex stream models. To facilitate this, we develop a new, modular dynamical modelling code sstrax for stellar streams that is highly accelerated using jax. We test our analysis pipeline on a mock observation that resembles the GD1 stream, and demonstrate that we can perform robust inference on all relevant parts of the stream model simultaneously. Finally, we present some outlook as to how this approach can be developed further to perform more complete and accurate statistical analyses of current and future data.

We present a catalogue of dynamical properties for 2368 late-type galaxies from the MaNGA survey. The latter complements the catalogue of photometric properties for the same sample based on deep optical DESI photometry processed with AutoProf. Rotation curves (RCs), extracted by model fitting H$\alpha$ velocity maps from the MaNGA Data Analysis Pipeline, extend out to 1.4 (1.9) R$_{e}$ for the primary (secondary) MaNGA samples. The RCs and ancillary MaNGA Pipe3D data products were used to construct various fundamental galaxy scaling relations that are also compared uniformly with similar relations from NIHAO zoom-in simulations. Simulated NIHAO galaxies were found to broadly reproduce the observed MaNGA galaxy population for $\log (M_*/{\rm M_{\odot}) > 8.5}$. Some discrepancies remain, such as those pertaining to central stellar densities and the diversity of RCs due to strong feedback schemes. Also presented are spatially-resolved scatters for the velocity-size-stellar mass (VRM$_*$) structural relations using MaNGA and NIHAO samples. The scatter for these relations in the galaxian interiors is a consequence of the diversity of inner RC shapes, while scatter in the outskirts is dictated by the large range of stellar surface densities which itself is driven by sporadic star formation. The detailed spatially-resolved scatter analysis highlights the complex interplay between local and global astrophysical processes and provides a strong constraint to numerical simulations.

This is the fifth paper in a series of investigations of the clustering properties of luminous, broad-emission line active galactic nuclei (AGN) identified in the ROSAT All-Sky Survey (RASS) and Sloan Digital Sky Survey (SDSS). In this work we measure the cross-correlation function (CCF) between RASS/SDSS DR14 AGN with the SDSS CMASS galaxy sample at $0.44<z<0.64$. We apply halo occupation distribution (HOD) modeling to the CCF along with the auto-correlation function of the CMASS galaxies. We find that X-ray and optically-selected AGN at $0.44<z<0.64$ reside in statistically identical halos with a typical dark matter halo mass of $M_{\rm DMH}^{\rm typ,AGN} \sim 10^{12.7}\,h^{-1}\,\rm{M}_\odot$. The acceptable HOD parameter space for these two broad-line AGN samples have only statistically marginal differences caused by small deviations of the CCFs in the 1-halo dominated regime on small scales. In contrast to optically-selected AGN, the X-ray AGN sample may contain a larger population of satellites at $M_{\rm DMH} \sim 10^{13}\,h^{-1}\,\rm{M}_\odot$. We compare our measurements in this work with our earlier studies at lower independent redshift ranges, spanning a look-back time of 6 Gyr. The comparison over this wider redshift range of $0.07<z<0.64$ reveals: (i) no significant difference between the typical DMH masses of X-ray and optically-selected AGN, (ii) weak positive clustering dependencies of $M_{\rm DMH}^{\rm typ,AGN}$ with $L_X$ and $M_{\rm BH}$, (iii) no significant dependence of $M_{\rm DMH}^{\rm typ,AGN}$ on Eddington ratio, and (iv) the same DMH masses host more massive accreting black holes at high redshift than at low redshifts.

We present new HI observations of the nearby massive spiral galaxy M83, taken with the VLA at $21^{\prime\prime}$ angular resolution ($\approx500$ pc) of an extended ($\sim$1.5 deg$^2$) 10-point mosaic combined with GBT single dish data. We study the super-extended HI disk of M83 (${\sim}$50 kpc in radius), in particular disc kinematics, rotation and the turbulent nature of the atomic interstellar medium. We define distinct regions in the outer disk ($r_{\rm gal}>$central optical disk), including ring, southern area, and southern and northern arm. We examine HI gas surface density, velocity dispersion and non-circular motions in the outskirts, which we compare to the inner optical disk. We find an increase of velocity dispersion ($\sigma_v$) towards the pronounced HI ring, indicative of more turbulent HI gas. Additionally, we report over a large galactocentric radius range (until $r_{\rm gal}{\sim}$50 kpc) that $\sigma_v$ is slightly larger than thermal (i.e. $>8$km s$^{-1}$ ). We find that a higher star formation rate (as traced by FUV emission) is not always necessarily associated with a higher HI velocity dispersion, suggesting that radial transport could be a dominant driver for the enhanced velocity dispersion. We further find a possible branch that connects the extended HI disk to the dwarf irregular galaxy UGCA365, that deviates from the general direction of the northern arm. Lastly, we compare mass flow rate profiles (based on 2D and 3D tilted ring models) and find evidence for outflowing gas at r$_{\rm gal}$ $\sim$2 kpc, inflowing gas at r$_{\rm gal}$ $\sim$5.5 kpc and outflowing gas at r$_{\rm gal}$ $\sim$14 kpc. We caution that mass flow rates are highly sensitive to the assumed kinematic disk parameters, in particular, to the inclination.

We measure the mean free path ($\lambda_{\rm mfp,HI}$), photo-ionization rate ($\langle \Gamma_{\rm HI} \rangle$) and neutral fraction ($\langle f_{\rm HI} \rangle$) of hydrogen in 12 redshift bins at $4.85<z<6.05$ from a large sample of moderate resolution XShooter and ESI QSO absorption spectra. The fluctuations in ionizing radiation field are modeled by post-processing simulations from the Sherwood suite using our new code ''EXtended reionization based on the Code for Ionization and Temperature Evolution'' (EX-CITE). EX-CITE uses efficient Octree summation for computing intergalactic medium attenuation and can generate large number of high resolution $\Gamma_{\rm HI}$ fluctuation models. Our simulation with EX-CITE shows remarkable agreement with simulations performed with the radiative transfer code Aton and can recover the simulated parameters within $1\sigma$ uncertainty. We measure the three parameters by forward-modeling the Ly$\alpha$ forest and comparing the effective optical depth ($\tau_{\rm eff, HI}$) distribution in simulations and observations. The final uncertainties in our measured parameters account for the uncertainties due to thermal parameters, modeling parameters, observational systematics and cosmic variance. Our best fit parameters show significant evolution with redshift such that $\lambda_{\rm mfp,HI}$ and $\langle f_{\rm HI} \rangle$ decreases and increases by a factor $\sim 6$ and $\sim 10^{4}$, respectively from $z \sim 5$ to $z \sim 6$. By comparing our $\lambda_{\rm mfp,HI}$, $\langle \Gamma_{\rm HI} \rangle$ and $\langle f_{\rm HI} \rangle$ evolution with that in state-of-the-art Aton radiative transfer simulations and the Thesan and CoDa-III simulations, we find that our best fit parameter evolution is consistent with a model in which reionization completes by $z \sim 5.2$.

We introduce grlic, a publicly available Python tool for generating glass-like point distributions with a radial density profile $n(r)$ as it is observed in large-scale surveys of galaxy distributions on the past light cone. Utilising these glass-like catalogues, we assess the bias and variance of the Landy-Szalay (LS) estimator of the first three two-point correlation function (2PCF) multipoles in halo and particle catalogues created with the cosmological N-body code gevolution. Our results demonstrate that the LS estimator calculated with the glass catalogues is biased by less than $10^{-4}$ with respect to the estimate derived from Poisson-sampled random catalogues, for all multipoles considered and on all but the smallest scales. Additionally, the estimates derived from glass-like catalogues exhibit significantly smaller standard deviation than estimates based on commonly used Poisson-sampled random catalogues of comparable size. The standard deviation $\sigma$ of the estimate depends on a power of the number of objects $N_R$ in the random catalogue; we find a power law $\sigma \propto \alpha^{-0.9}$ for glass-like random catalogues as opposed to $\sigma \propto N_R^{-0.48}$ using Poisson-sampled random catalogues. Given a required precision, this allows for a much reduced number of objects in the glass-like random catalogues used for the LS estimate of the 2PCF multipoles, significantly reducing the computational costs of each estimate.

Due to their aerodynamical coupling with gas, pebbles in protoplanetary discs can drift over large distances to support planet growth in the inner disc. In the past decade, this pebble accretion has been studied extensively for aerodynamically small pebbles (Stokes number St < 1). However, accretion can also operate in the St > 1 mode, e.g., when planetesimals collisionally fragment to smaller bodies or when the primordial gas disc disperses. This work aims to extend the study of pebble accretion to these aerodynamically loosely coupled particles. We integrate the pebble's equation of motion, accounting for gas drag, stellar and planetary gravity, in the midplane of a laminar disc. The accretion probability ($\epsilon$) is calculated as function of Stokes number, disc pressure gradient index, planet mass and eccentricity. We find that for Stokes number above unity $\epsilon$(St) first rises, due to lower drift and aided by a large atmospheric capture radius, until it reaches a plateau where the efficiency approaches 100 per cent. At high St the plateau region terminates as particles become trapped in resonance. These results are well described by a semi-analytical "kick-and-drift" model and we also provide fully analytical prescriptions for $\epsilon$. We apply our model to the accretion of $\sim 30 \mu$m dust particles in a dispersing protoplanetary and secondary (CO-rich) debris disc. It shows that physically small particles are mainly accreted as aerodynamically large Stokes number pebbles during the debris disc phase. Earth-mass planets may obtain $\sim 25$ per cent of their heavy elements through this late accretion phase.

The formation of a cold Jupiter (CJ) is expected to quench the influx of pebbles and the migration of cores interior to its orbit, thus limiting the efficiency of rocky planet formation either by pebble accretion and/or orbital migration. Observations, however, show that the presence of outer CJs ( >1 au and >0.3 Jupiter masses) correlates with the presence of inner Super Earths (at <1 au). This observation may simply be a result of an enhanced initial reservoir of solids in the nebula required to form a CJ or a yet-to-be-determined mechanism assisted by the presence of the CJ. In this work, we focus on the latter alternative and study the orbital transport of planetesimals interior to a CJ subject to the gravity and drag from a viscously-evolving gaseous disk. We find that a secular resonance sweeping inwards through the disk gradually transports rings of planetesimals when their drag-assisted orbital decay is faster than the speed of the resonance scanning. This snowplow-like process leads to large concentration (boosted by a factor of ~10-100) of size-segregated planetesimal rings with aligned apsidal lines, making their expected collisions less destructive due to their reduced velocity dispersion. This process is efficient for a wide range of alpha-disk models and Jovian masses, peaking for ~1-5 Jupiter masses, typical of observed CJs in radial velocity surveys. Overall, our work highlights the major role that the disk's gravity may have on the orbital redistribution of planetesimals, depicting a novel avenue by which CJs may enhance the formation of inner planetary systems, including super-Earths and perhaps even warm and hot Jupiters.

We present long-term (2012-2022) optical monitoring of the candidate black hole X-ray binary Swift J1910.2-0546 with the Faulkes Telescopes and Las Cumbres Observatory (LCO) network. Following its initial bright 2012 outburst, we find that the source displayed a series of at least 7 quasi-periodic, high amplitude (~3 mags) optical reflares in 2013, with a recurrence time increasing from ~42 days to ~49 days. In 2014, the source experienced a mini-outburst with two peaks in the optical. We also study the recent 2022 outburst of the source at optical wavelengths, and perform a comparative analysis with the earlier rebrightenings. A single X-ray detection and only two radio detections were obtained during the 2013 reflaring period, and only optical detections were acquired in 2014. During the reflaring in both 2013 and 2014, the source showed bluer-when-brighter behavior, having optical colors consistent with a blackbody heating and cooling between 4500 and 9500 K, i.e. the temperature range in which hydrogen starts to ionize. Finally, we compare the flaring behavior of the source to re-brightening events in other X-ray binaries. We show that the repeated reflarings of Swift J1910.2-0546 are highly unusual, and propose that they arise from a sequence of repetitive heating and cooling front reflections travelling through the accretion disk.

Utilizing the James Webb Space Telescope Early Release NIRCam Observations, we perform a weak-lensing analysis of the massive galaxy cluster SMACS J0723.3-7327 ($z=0.39$). We investigate the spatial variation of the PSF from the stars in the mosaic image. Our measurements show that the PSF for both modules has very small spatial and temporal variation with average complex ellipticity components of $e_1=0.007\pm0.001$ and $e_2=0.029\pm0.001$ in the observed north-up reference frame. We create PSF models through a principal component analysis of the stars and show that they properly account for the ellipticity of the PSF with residual shapes of $e_1=(0.3\pm3.5)\times10^{-4}$ and $e_2=(1.8\pm4.0)\times10^{-4}$. We select background galaxies by their photometric redshift and measure galaxy shapes by model fitting. Our weak-lensing source catalog achieves 215 galaxies arcmin$^{-2}$. We map the projected mass density of SMACSJ0723 and detect the cluster with a peak significance of $12.2\sigma$. The mass distribution is found to elongate in the east-west direction with an extension to the northeast edge of the field of view where a candidate substructure is found in the Chandra X-ray imaging. We fit the tangential shear with a Navarro-Frenk-White model and estimate the mass of the cluster to be $M_{500}=7.9\pm1.1\times10^{14}$ M$_{\odot}$ ($M_{200}=11.4\pm1.5\times10^{14}$ M$_\odot$ ), which agrees with existing mass estimates. Combining the multiwavelength evidence from literature with our weak-lensing analysis, we hypothesize that SMACSJ0723 is observed near first pericenter passage and we identify candidate radio relics.

Previous studies analyzing the evanescent nature of acoustic waves in the lower solar atmosphere, up to 300\,km above the photosphere, have shown an unexpected phase shift of an order of 1\,s between different heights. Those studies investigated the spectral line \ion{Fe}{1} 6173.3\,\AA, commonly used for helioseismic measurements. Such phase-shifts can contribute to a misinterpretation of the measured travel times in local helioseismology, complicating inferences of, e.g., the deep meridional flow. In this study, we carry out phase-shift computations using a simulated, fully radiative, and convective atmosphere from which the \ion{Fe}{1} 6173.3\,\AA\ line is synthesized. The resulting phase-shifts as functions of frequency across multiple heights show non-zero values in evanescent waves, similar to what was found in observational data. Comparing the Doppler-velocities estimated from the synthesized absorption line with the true velocities directly obtained from the simulated plasma motions, we find substantial differences in phase-shifts between the two. This leads us to hypothesize that the non-adiabaticity of the solar atmosphere yields extra phase-shift contributions to Doppler velocities. Finally, computing phase-differences for different viewing angles reveals a systematic center-to-limb variation, similar to what is present in observations. Overall, this study helps to improve our understanding of the physical cause of the helioseismic center-to-limb effect.

The opening angles of some protostellar outflows appear too narrow to match the expected core-star mass efficiency SFE = 0.3-0.5 if outflow cavity volume traces outflow mass, with a conical shape and a maximum opening angle near 90 deg. However, outflow cavities with paraboloidal shape and wider angles are more consistent with observed estimates of the SFE. This paper presents a model of infall and outflow evolution based on these properties. The initial state is a truncated singular isothermal sphere which has mass $\approx$1 $M_\odot$, free fall time $\approx$80 kyr, and small fractions of magnetic, rotational, and turbulent energy. The core collapses pressure-free as its protostar and disk launch a paraboloidal wide-angle wind. The cavity walls expand radially and entrain envelope gas into the outflow. The model matches SFE values when the outflow mass increases faster than the protostar mass by a factor 1 - 2, yielding protostar masses typical of the IMF. It matches observed outflow angles if the outflow mass increases at nearly the same rate as the cavity volume. The predicted outflow angles are then typically $\sim$50 deg as they increase rapidly through the stage 0 duration of $\sim$40 kyr. They increase more slowly up to $\sim$110 deg during their stage I duration of $\sim$70 kyr. With these outflow rates and shapes, model predictions appear consistent with observational estimates of typical stellar masses, SFEs, stage durations, and outflow angles, with no need for external mechanisms of core dispersal.

We present the results of the first fully cosmological hydrodynamical simulations studying the merger-driven model for massive black hole (BH) seed formation via direct collapse. Using the zoom-in technique as well as particle splitting, we achieve a final spatial resolution of $2$ pc. We show that the major merger of two massive galaxies at redshift $z \sim 8$ results in the formation of a nuclear supermassive disk (SMD) of only $4$ pc in radius, owing to a prodigious gas inflow sustained at $100$-$1000$ $M_{\odot}$ yr$^{-1}$. The core of the merger remnant is metal-rich, well above solar abundance, and the SMD reaches a gaseous mass of $3 \times 10^8$ $M_{\odot}$ in less than a million years after the merger, despite a concurrent prominent nuclear starburst. Dynamical heating as gas falls into the deepest part of the potential well, and heating and stirring by supernova blastwaves, generate a turbulent multi-phase interstellar medium, with a gas velocity dispersion exceeding 100 km s$^{-1}$. As a result, only moderate fragmentation occurs in the inner $10$-$20$ pc despite the temperature falls below $1000$ K. The SMD is Jeans-unstable as well as bar-unstable and will collapse further adiabatically, becoming warm and ionized. We show that the SMD, following inevitable contraction, will become general relativistic unstable and directly form a supermassive BH of mass in the range $10^6$-$10^8$ $M_{\odot}$, essentially skipping the stage of BH seed formation. These results confirm that mergers between the most massive galaxies at $z \sim 8$-$10$ can naturally explain the rapid emergence of bright high-redshift quasars.

Next-generation large segmented mirror telescopes are expected to perform direct imaging and characterization of Earth-like rocky planets, which requires contrast limits of $10^{-7}$ to $10^{-8}$ at wavelengths from I to J band. One critical aspect affecting the raw on-sky contrast are polarization aberrations arising from the reflection from the telescope's mirror surfaces and instrument optics. We simulate the polarization aberrations and estimate their effect on the achievable contrast for three next-generation ground-based large segmented mirror telescopes. We performed ray-tracing in Zemax and computed the polarization aberrations and Jones pupil maps using the polarization ray-tracing algorithm. The impact of these aberrations on the contrast is estimated by propagating the Jones pupil maps through a set of idealized coronagraphs using hcipy, a physical optics-based simulation framework. The optical modeling of the giant segmented mirror telescopes (GSMTs) shows that polarization aberrations create significant leakage through a coronagraphic system. The dominant aberration is retardance defocus, which originates from the steep angles on the primary and secondary mirrors. The retardance defocus limits the contrast to $10^{-5}$ to $10^{-4}$ at 1 $\lambda/D$ at visible wavelengths, and $10^{-5}$ to $10^{-6}$ at infrared wavelengths. The simulations also show that the coating plays a major role in determining the strength of the aberrations. Polarization aberrations will need to be considered during the design of high-contrast imaging instruments for the next generation of extremely large telescopes. This can be achieved either through compensation optics, robust coronagraphs, specialized coatings, calibration, and data analysis approaches or by incorporating polarimetry with high-contrast imaging to measure these effects.

The census of obscured quasar populations is incomplete, and remains a major unsolved problem, especially at higher redshifts, where we expect a greater density of galaxy formation and quasar activity. We present Gemini GNIRS near-infrared spectroscopy of 24 luminous obscured quasar candidates from the Sloan Digital Sky Survey's Stripe 82 region. The targets were photometrically selected using a WISE/W4 selection technique that is optimized to identify IR-bright and heavily-reddened/optically-obscured targets at $z>1$. We detect emission lines of ${\rm H\alpha}$, ${\rm H\beta}$, and/or ${\rm[ O~III]}$ in 23 sources allowing us to measure spectroscopic redshifts in the range $1<z<3$ with bolometric luminosities spanning $L=10^{46.3}-10^{47.3}$ erg s$^{-1}$. We observe broad $10^3-10^4$ km s$^{-1}$ Balmer emissions with large ${\rm H\alpha}/{\rm H\beta}$ ratios, and we directly observe a heavily reddened rest-frame optical continuum in several sources, suggesting high extinction ($A_V\sim7-20$ mag). Our observations demonstrate that such optical/infrared photometric selection successfully recovers high-redshift obscured quasars. The successful identification of previously undetected red, obscured high-redshift quasar candidates suggests that there are more obscured quasars yet to be discovered.

Various models of modified gravity invoke ``screening'' mechanisms that are sensitive to the value of the local gravitational potential. This could have observable consequences for galaxies. These consequences might be seen by comparing two proxies for galaxy mass -- their luminosity and their internal kinematics -- as a function of local galaxy density. Motivated by this prospect, we have compared the observed properties of luminous red galaxies (LRGs) inside and outside of voids in the cosmic large scale structure. We used archival measurements of line widths, luminosities, redshifts, colors, and positions of galaxies in conjunction with recent void catalogs to construct comparison LRG samples inside and outside of voids. We fitted these two samples to the well-established fundamental plane of elliptical galaxies to constrain any differences between the inferred value of the Newtonian gravitational constant G for the two samples. We obtained a null result, with an upper limit on any fractional difference in G within and outside of cosmological voids to be $\alpha =\delta$$ G/G \sim$ 40\%. This upper bound is dominated by the small-number statistics of our N $\sim $ 100 within-void LRG sample. With the caveat that environmental effects could influence various parameters such as star formation, we estimate that a 1\% statistical limit on $\alpha$ could be attained with data from 10${^5}$ elliptical galaxies within voids. This is within the reach of future photometric and spectroscopic surveys, both of which are required to pursue this method.

As LIGO-Virgo-KAGRA enters its fourth observing run, a new opportunity to search for electromagnetic counterparts of compact object mergers will also begin. The light curves and spectra from the first "kilonova" associated with a binary neutron star binary (NSM) suggests that these sites are hosts of the rapid neutron capture ("$r$") process. However, it is unknown just how robust elemental production can be in mergers. Identifying signposts of the production of particular nuclei is critical for fully understanding merger-driven heavy-element synthesis. In this study, we investigate the properties of very neutron rich nuclei for which superheavy elements ($Z\geq 104$) can be produced in NSMs and whether they can similarly imprint a unique signature on kilonova light-curve evolution. A superheavy-element signature in kilonovae represents a route to establishing a lower limit on heavy-element production in NSMs as well as possibly being the first evidence of superheavy element synthesis in nature. Favorable NSMs conditions yield a mass fraction of superheavy elements is $X_{Z\geq 104}\approx 3\times 10^{-2}$ at 7.5 hours post-merger. With this mass fraction of superheavy elements, we find that kilonova light curves may appear similar to those arising from lanthanide-poor ejecta. Therefore, photometric characterizations of superheavy-element rich kilonova may possibly misidentify them as lanthanide-poor events.

The most reliable single-epoch supermassive black hole mass ($M_{\rm BH}$) estimates in quasars are obtained by using the velocity widths of low-ionization emission lines, typically the H$\beta$ $\lambda4861$ line. Unfortunately, this line is redshifted out of the optical band at $z\approx1$, leaving $M_{\rm BH}$ estimates to rely on proxy rest-frame ultraviolet (UV) emission lines, such as C IV $\lambda1549$ or Mg II $\lambda2800$, which contain intrinsic challenges when measuring, resulting in uncertain $M_{\rm BH}$ estimates. In this work, we aim at correcting $M_{\rm BH}$ estimates derived from the C IV and Mg II emission lines based on estimates derived from the H$\beta$ emission line. We find that employing the equivalent width of C IV in deriving $M_{\rm BH}$ estimates based on Mg II and C IV provides values that are closest to those obtained from H$\beta$. We also provide prescriptions to estimate $M_{\rm BH}$ values when only C IV, only Mg II, and both C IV and Mg II are measurable. We find that utilizing both emission lines, where available, reduces the scatter of UV-based $M_{\rm BH}$ estimates by $\sim15\%$ when compared to previous studies. Lastly, we discuss the potential of our prescriptions to provide more accurate and precise estimates of $M_{\rm BH}$ given a much larger sample of quasars at $3.20 \lesssim z \lesssim 3.50$, where both Mg II and H$\beta$ can be measured in the same near-infrared spectrum.

Initial orbit determination (IOD) from line-of-sight (i.e., bearing) measurements is a classical problem in astrodynamics. Indeed, there are many well-established methods for performing the IOD task when given three line-of-sight observations at known times. Interestingly, and in contrast to these existing methods, concepts from algebraic geometry may be used to produce a purely geometric solution. This idea is based on the fact that bearings from observers in general position may be used to directly recover the shape and orientation of a three-dimensional conic (e.g., a Keplerian orbit) without any need for knowledge of time. In general, it is shown that five bearings at unknown times are sufficient to recover the orbit -- without the use of any type of initial guess and without the need to propagate the orbit. Three bearings are sufficient for purely geometric IOD if the orbit is known to be (approximately) circular. The method has been tested over different scenarios, including one where extra observations make the system of equations over-determined.

Using time-distance local helioseismology flow maps within 1 Mm of the solar photosphere, we detect inflows toward activity belts that contribute to solar cycle scale variations in near-surface meridional flow. These inflows stretch out as far as 30 degrees away from active region centroids. If active region neighborhoods are excluded, the solar cycle scale variation in background meridional flow diminishes to below 2~m~s$^{-1}$, but still shows systematic variations in the absence of active regions between Sunspot Cycles 24 and 25. We, therefore, propose that the near-surface meridional flow is a three component flow made up of: a constant baseline flow profile that can be derived from quiet Sun regions, variations due to inflows around active regions, and solar cycle scale variation of the order of 2~m~s$^{-1}$. Torsional oscillation, on the other hand, is found to be a global phenomenon i.e. exclusion of active region neighborhoods does not affect its magnitude or phase significantly. This non-variation of torsional oscillation with distance away from active regions and the three-component breakdown of the near-surface meridional flow serve as vital constraints for solar dynamo models and surface flux transport simulations.

We have conducted numerical simulations to reproduce the observed optical energy profile of the 15 October 2021 (UT) impact flash on Jupiter, which was the largest and the most well-observed flash event detected by ground-based movie observations. The observed long-duration ($\sim 5.5~{\rm s}$) optical emission can be reproduced by an impact of an object with an exceptionally small angle of entry relative to the horizontal. The apparent lack of the impact debris feature despite the large impact object was possibly due to the shallower angle of entry ($\le 12^\circ$), which resulted in the lower ablation per unit volume at altitudes higher than $50 \, {\rm km}$, and the volume densities of the ablated materials were too low to allow the debris particulates to coagulate. The absence of temporal methane absorption change in the observed flash spectrum is consistent with the best-fit results. The model better fits the observed optical energy profile for weaker material (cometary and stony) cases than for metallic ones. Based on the simulation results, prospects for future observations of impact flashes are discussed.

The China Space Station Telescope (CSST) will enter a low Earth orbit around 2024 and operate for 10 years, with seven of those years devoted to surveying the area of the median-to-high Galactic latitude and median-to-high Ecliptic latitude of the sky. To maximize the scientific output of CSST, it is important to optimize the survey schedule. We aim to evaluate the astrometric capability of CSST for a given survey schedule and to provide independent suggestions for the optimization of the survey strategy. For this purpose, we first construct the astrometric model and then conduct simulated observations based on the given survey schedule. The astrometric solution is obtained by analyzing the simulated observation data. And then we evaluate the astrometric capability of CSST by analyzing the properties of the astrometric solution. We find that the accuracy of parallax and proper motion of CSST is better than 1 mas( yr1) for the sources of 18-22 mag in g band, and about 1-10 mas( yr1) for the sources of 22-26 mag in g band, respectively. The results from real survey could be worse since the assumptions are optimistic and simple. We find that optimizing the survey schedule can improve the astrometric accuracy of CSST. In the future, we will improve the astrometric capability of CSST by continuously iterating and optimizing the survey schedule.

Most stars in our Galaxy form in stellar aggregates, which can become long-lived structures called open clusters (OCs). Along their dynamical evolution, their gradual depletion leave some imprints on their structure. In this work, we employed astrometric, photometric and spectroscopic data from the \textit{Gaia} DR3 catalogue to uniformly characterize a sample of 60 OCs. Structural parameters (tidal, core and half-light radii, respectively, $r_t$, $r_c$ and $r_h$), age, mass ($M_{\textrm{clu}}$), distance, reddening, besides Jacobi radius ($R_J$) and half-light relaxation time ($t_{rh}$), are derived from radial density profiles and astrometrically decontaminated colour-magnitude diagrams. Ages and Galactocentric distances ($R_G$) range from 7.2$\,\lesssim\,$log($t.$yr$^{-1}$)$\,\lesssim\,$9.8 and 6$\,\lesssim\,R_G$(kpc)$\,\lesssim\,$12. Analytical expressions derived from $N$-body simulations, taken from the literature, are also employed to estimate the OC initial mass ($M_{\textrm{ini}}$) and mass loss due to exclusively dynamical effects. Both $r_c$ and the tidal filling ratio, $r_h/R_J$, tend to decrease with the dynamical age (=$t/t_{rh}$), indicating the shrinking of the OCs' internal structure as consequence of internal dynamical relaxation. This dependence seems differentially affected by the external tidal field, since OCs at smaller $R_G$ tend to be dynamically older and have smaller $M_{\textrm{clu}}/M_{\textrm{ini}}$ ratios. In this sense, for $R_G\lesssim8\,$kpc, the $r_h/R_J$ ratio presents a slight positive correlation with $R_G$. Beyond this limit, there is a dichotomy in which more massive OCs tend to be more compact and therefore less subject to tidal stripping in comparison to those less massive and looser OCs at similar $R_G$. Besides, the $r_t/R_J$ ratio also tends to correlate positively with $R_G$.

Disintegrating multiple systems have been previously discovered from kinematic studies of the $\it Hipparcos$ catalogue. They are presumably the result of dynamical encounters taking place in the Galactic disk between single/multiple systems. In this paper, we aim to expand the search for such systems, to study their properties, as well as to characterize possible low-mass ejecta (i.e. brown dwarfs and planets). We have assembled a list of 15 candidate systems using astrometry from the Tycho-Gaia astrometric solution (later upgraded with $\it Gaia$ DR3), and here we present the discovery and follow-up of 5 of them. We have obtained DECam imaging for all 5 systems and by combining near-infrared photometry and proper motion, we searched for ultra-cool ejected components. We find that the system consisting of TYC 7731-1951-1, TYC 7731-2128 AB, and TYC 7731-1995-1ABC?, contains one very promising ultra-cool dwarf candidate. Using additional data from the literature, we have found that 3 out of 5 disintegrating system candidates are likely to be true disintegrating systems.

As the most abundant element in the universe after hydrogen and helium, oxygen plays a key role in planetary, stellar, and galactic astrophysics. Its abundance is especially influential on stellar structure and evolution, and as the dominant opacity contributor at the base of the Sun's convection zone it is central to the discussion around the solar modelling problem. However, abundance analyses require complete and reliable sets of atomic data. We present extensive atomic data for O I, by using the multiconfiguration Dirac-Hartree-Fock and relativistic configuration interaction methods. Lifetimes and transition probabilities for radiative electric dipole transitions are given and compared with results from previous calculations and available measurements. The accuracy of the computed transition rates is evaluated by the differences between the transition rates in Babushkin and Coulomb gauges, as well as by a cancellation factor analysis. Out of the 989 computed transitions in this work, 205 are assigned to the accuracy classes AA-B, that is, with uncertainties less than 10%, following the criteria defined by the National Institute of Standards and Technology Atomic Spectra Database. We discuss the influence of the new log(gf) values on the solar oxygen abundance and ultimately advocate $\log\epsilon_{\mathrm{O}}=8.70\pm0.04$.

Aims. We estimate the dynamical evolution of the Globular Clusters interaction with the Galactic centre that dynamically changed in the past. Methods. We simulated the orbits of 147 globular clusters over 10 Gyr lookback time using the parallel N-body code phi-GPU. For each globular cluster, we generated 1000 sets of initial data with random proper motions and radial velocities based on the observed values. To distinguish globular clusters interacting with the galactic centre, we used the criterion of a relative distance of less than 100 pc. We used four external potentials from the IllustrisTNG-100 database, which were selected for their similarity to the present-day Milky Way, to simulate the structure of the Galaxy at different times. Results. We obtained 3-4 globular cluster interactions per Gyr at distances of less than 50 pc and 5-6 interactions per Gyr at distances of less than 80 pc among the studied 147 globular clusters that had close passages near the Galactic centre. We selected 10 of them for detailed study and found almost 100% probability of interaction with the Galactic centre for six of them. Conclusions. According to our results, the maximum interaction frequency of globular clusters with the Galactic centre in the Milky Way is likely to be a few dozens of passages per Gyr within a central zone of 100 pc. This low frequency may not be sufficient to fully explain the relatively high mass (of order 10^7 Msol) of the nuclear star cluster in the Milky Way, if we consider only the periodic capture of stars from globular clusters during close encounters. Therefore, we must also consider the possibility that some early globular clusters were completely tidally disrupted during interactions with the forming nuclear star cluster and the Galactic centre.

Cherenkov techniques are widely used in astroparticle experiments. This article reviews the various detection principles and the corresponding experiments, including some of the physics breakthroughs. In particular, it traces the development since the mid of the 1990s, a period when the field took a particularly dynamic development.

ERIS, the Enhanced Resolution Imager and Spectrograph, is an instrument that both extends and enhances the fundamental diffraction limited imaging and spectroscopy capability for the VLT. It replaces two instruments that were being maintained beyond their operational lifetimes, combines their functionality on a single focus, provides a new wavefront sensing module for natural and laser guide stars that makes use of the Adaptive Optics Facility, and considerably improves on their performance. The observational modes ERIS provides are integral field spectroscopy at 1-2.5 {\mu}m, imaging at 1-5 {\mu}m with several options for high contrast imaging, and longslit spectroscopy at 3-4 {\mu}m, The instrument is installed at the Cassegrain focus of UT4 at the VLT and, following its commissioning during 2022, has been made available to the community.

In wind-fed X-ray binaries, the radiatively driven wind of the primary star can be suppressed by the EUV irradiation of the compact secondary star, leading to an increased accretion rate. This causes feedback between the released accretion power and the luminosity of the compact star. We investigate the feedback process between the released accretion power and the X-ray luminosity of the compact star in the unique high-mass X-ray binary Cygnus X-3. We assume that a part of the wind-fed power experiences a small amplitude variability around the source luminosity. We propose a simple heuristic model to couple the influence of EUV irradiation on the stellar wind (from the Wolf-Rayet companion star) with the X-ray source itself. The resulting time profile of luminosity mimics that of the input variability, albeit with a larger amplitude. The most important property of the input variability are turnover times when it changes its sign and starts to have either positive or negative feedback. The bolometric luminosity derived by spectral modeling is the time average of the resulting feedback luminosity. We demonstrate that the erratic behavior of the X-ray light curve of Cygnus X-3 may have its origin in the small amplitude variability of the X-ray source and feedback with the companion wind. This variability could arise in the accretion flow and/or due to the loss of kinetic energy in a jet or an accretion disk wind. In order to produce similar properties of the simulated light curve as observed, we have to restrict the largest accretion radius to a changing level, and assume variable timescales for the rise and decline phases of the light curve.

Nearly all current simulations predict that outcomes of the star formation process, such as the fraction of stars that form in bound clusters (Gamma), depend on the intensity of star formation activity (SigmaSFR) in the host galaxy. The exact shape and strength of the predicted correlations, however, vary from simulation to simulation. Observational results also remain unclear at this time, because most works have mixed estimates made from very young clusters for galaxies with higher SigmaSFR with those from older clusters for galaxies with lower SigmaSFR. The three blue compact dwarf (BCD) galaxies ESO185-IG13, ESO338-IG04, and Haro11 have played a central role on the observational side because they have some of the highest known SigmaSFR and published values of Gamma. We present new estimates of Gamma for these BCDs in three age intervals (1-10 Myr, 10-100 Myr, 100-400 Myr), based on age-dating which includes Halpha photometry to better discriminate between clusters younger and older than ~10 Myr. We find significantly lower values for Gamma (1-10 Myr) than published previously. The likely reason for the discrepancy is that previous estimates appear to be based on age-reddening results that underestimated ages and overestimated reddening for many clusters, artificially boosting Gamma (1-10 Myr). We also find that fewer stars remain in clusters over time, with ~15-39% in 1-10 Myr, ~5-7% in 10-100 Myr, and ~1-2% in 100-400 Myr clusters. We find no evidence that Gamma increases with SigmaSFR. These results imply that cluster formation efficiency does not vary with star formation intensity in the host galaxy. If confirmed, our results will help guide future assumptions in galaxy-scale simulations of cluster formation and evolution.

We present observations with the Space Telescope Imaging Spectrograph (STIS) onboard the Hubble Space Telescope of eleven Lyman continuum (LyC) leaking galaxies at redshifts, z, in the range 0.29-0.43, with oxygen abundances 12+log(O/H)=7.64-8.16, stellar masses Mstar~10^7.8-10^9.8 Msun and O32=[OIII]5007/[OII]3727 of ~5-20, aiming to detect CIII]1908 emission line. We combine these observations with the optical Sloan Digital Sky Survey (SDSS) spectra for the determination of carbon, nitrogen and oxygen abundances. Our sample was supplemented by thirty one galaxies from the literature, for which carbon, nitrogen and oxygen abundances can be derived from the HST and SDSS spectra. These additional galaxies, however, do not have LyC observations. We find that log(C/O) for the entire sample at 12+log(O/H)<8.1 does not depend on metallicity, with a small dispersion of ~0.13 dex around the average value of ~ -0.75 dex. On the other hand, the log(N/O) in galaxies at z>0.1, including LyC leakers, is systematically higher compared to the rest of the sample with lower metallicity. We find that log(C/O) slightly decreases with increasing Mstar from ~ -0.65 at Mstar=10^6 Msun to ~ -0.80 at Mstar=10^9-10^10 Msun, whereas log(N/O) is considerably enhanced at Mstar>10^8 Msun. The origin of these trends remains basically unknown. One of the possible solutions is to assume that the upper mass limit of the stellar initial mass function (IMF) in more massive galaxies is higher. This would result in higher production of oxygen and larger fraction of massive stars with stellar wind polluting interstellar medium with nitrogen.

Using high-resolution spectra from the Tillinghast Reflector Echelle Spectrograph (TRES) and photometry from sector 56 of the Transiting Exoplanet Survey Satellite (TESS), we report that the nearby M dwarf G 68-34 is a double-lined eclipsing binary. The pair is spin-orbit synchronized with a period of 0.655 days. The light curve shows significant spot modulation with a larger photometric amplitude than that of the grazing eclipses. We perform a joint fit to the spectroscopic and photometric data, obtaining masses of $0.3280\pm 0.0034$M$_\odot$ and $0.3207\pm 0.0036$M$_\odot$ and radii of $0.345\pm 0.014$R$_\odot$ and $0.342\pm 0.014$R$_\odot$ after marginalizing over unknowns in the starspot distribution. This system adds to the small but growing population of fully convective M dwarfs with precisely measured masses and radii that can be used to test models of stellar structure. The pair also has a white dwarf primary at 9" separation, with the system known to be older than 5 Gyr from the white-dwarf cooling age. The binarity of G 68-34 confirms our hypothesis from Pass et al. (2022): in that work, we noted that G 68-34 was both rapidly rotating and old, highly unusual given our understanding of the spindown of M dwarfs, and that a close binary companion may be responsible.

We have conducted a search for $z\simeq7$ Lyman break galaxies over 8.2 square degrees of near-infrared imaging from the VISTA Deep Extragalactic Observations (VIDEO) survey in the XMM-Newton - Large Scale Structure (XMM-LSS) and the Extended Chandra Deep Field South (ECDF-S) fields. Candidate galaxies were selected from a full photometric redshift analysis down to a $Y+J$ depth of 25.3 ($5\sigma$), utilizing deep auxiliary optical and Spitzer/IRAC data to remove brown dwarf and red interloper galaxy contaminants. Our final sample consists of 28 candidate galaxies at $6.5\le z \le7.5$ with $-23.5 \le M_{\mathrm{UV}} \le -21.6$. We derive stellar masses of $9.1 \le \mathrm{log}_{10}(M/M_{\odot}) \le 10.9$ for the sample, suggesting that these candidates represent some of the most massive galaxies known at this epoch. We measure the rest-frame UV luminosity function (LF) at $z\simeq7$, confirming previous findings of a gradual decline in number density at the bright-end ($M_{\mathrm{UV}} < -22$) that is well described by a double-power law (DPL). We show that quasar contamination in this magnitude range is expected to be minimal, in contrast to conclusions from recent pure-parallel Hubble studies. Our results are up to a factor of ten lower than previous determinations from optical-only ground-based studies at $M_{\rm UV} \lesssim - 23$. We find that the inclusion of $YJHK_{s}$ photometry is vital for removing brown-dwarf contaminants, and $z \simeq 7$ samples based on red-optical data alone could be highly contaminated ($\gtrsim 50$ per cent). In comparison with other robust $z > 5$ samples, our results further support little evolution in the very bright-end of the rest-frame UV LF from $z = 5-10$, potentially signalling a lack of mass quenching and/or dust obscuration in the most massive galaxies in the first Gyr.

We report the discovery and spectroscopic identification of the bright doubly lensed quasar eRASS1 J050129.5-073309 at redshift $z=2.47$, selected from the first all-sky survey of the ${\it Spectrum\; Roentgen\; Gamma\; (SRG)}$ eROSITA telescope and the ${\it Gaia}$ EDR3 catalog. We systematically search for extragalactic sources with eROSITA X-ray positions having multiple ${\it Gaia}$ counterparts and have started spectroscopic follow-up of the most promising candidates using long-slit spectroscopy with NTT/EFOSC2 to confirm the lens nature. The two images are separated by $2.7''$ and their average ${\it Gaia}$ ${\it g}$-band magnitudes are 16.95 and 17.33. Legacy Survey DR10 imaging and image modeling reveal both the lensing galaxy and tentatively the lensed image of the quasar host galaxy. Archival optical light curves show evidence of a variability time delay with the fainter component lagging the brighter by about 100 days. The fainter image has also decreased its brightness by about 1 magnitude since 2019. This dimming was still obvious at the time of the spectroscopic observations and is probably caused by microlensing. The optical spectroscopic follow-up obtained from NTT/EFOSC2 and the evidence provided by the imaging and timing analysis allow us to confirm the lensed nature of eRASS1 J050129.5-073309.

The Lyapunov Characteristic Exponents are a useful indicator of chaos in astronomical dynamical systems. They are usually computed through a standard, very efficient and neat algorithm published in 1980. However, for Hamiltonian systems the expected result of pairs of opposite exponents is not always obtained with enough precision. We find here why in these cases the initial order of the deviation vectors matters, and how to sort them in order to obtain a correct result.

Main-sequence intermediate-mass stars present a radiative envelope that supports internal gravity waves (IGWs). Excited at the boundary with the convective core, IGWs propagate towards the stellar surface and are suspected to impact physical processes such as rotation and chemical mixing. Using the fully compressible time-implicit code MUSIC, we study IGWs in two-dimensional simulations of a zero-age-main-sequence 5 solar mass star model up to 91\% of the stellar radius with different luminosity and radiative diffusivity enhancements. Our results show that low frequency waves excited by core convection are strongly impacted by radiative effects as they propagate. This impact depends on the radial profile of radiative diffusivity which increases by almost 5 orders of magnitude between the centre of the star and the top of the simulation domain. In the upper layers of the simulation domain, we observe an increase of the temperature. Our study suggests that this is due to heat added in these layers by IGWs damped by radiative diffusion. We show that non-linear effects linked to large amplitude IGWs may be relevant just above the convective core. Both these effects are intensified by the artificial enhancement of the luminosity and radiative diffusivity, with enhancement factors up to $10^4$ times the realistic values. Our results also highlight that direct comparison between numerical simulations with enhanced luminosity and observations must be made with caution. Finally, our work suggests that thermal effects linked to the damping of IGWs could have a non-negligible impact on stellar structure.

We use low-redshift background cosmology data to place quantitative constraints on three separate modified gravity models, each of which aims to explain the low-redshift acceleration through a different physical mechanism. The Lifshitz cosmology is effectively a parametric extension of the canonical $\Lambda$CDM model, where a time-dependent cosmological constant originates from vacuum energy. The Infinite Statistics model is also a parametric extension of $\Lambda$CDM, where the dark energy is dynamic and originates from the curvature of a dual space-time. We show that the data restricts the additional parameters in these models to be consistent with their $\Lambda$CDM values, and in particular that it implies that the theoretically predicted value for a dimensionless coupling parameter in the Lifshitz model is ruled out at more than six standard deviations. In the Regge-Teitelboim model, gravity is described by embedding the usual space-time manifold in a fixed higher-dimensional background, and there is no parametric $\Lambda$CDM limit. We study several separate realizations of the model, respectively introduced by Davidson, by Fabi \textit{et al.}, and by Stern \& Xu, and show that the first two are ruled out by the low-redshift data we use, while the latter is consistent with this data but requires a non-standard value of the matter density. Overall, our analysis highlights the tight constraints imposed by current data on the allowed low-redshift deviations from the standard $\Lambda$CDM background evolution.

In the broad subclass of Horndeski theories with a luminal speed of gravitational waves, we derive gravitational waveforms emitted from a compact binary by considering the wave propagation on a spatially flat cosmological background. A scalar field nonminimally coupled to gravity gives rise to hairy neutron star (NS) solutions with a nonvanishing scalar charge, whereas black holes (BHs) do not have scalar hairs in such theories. A binary system containing at least one hairy neutron star modifies the gravitational waveforms in comparison to those of the BH-BH binary. Using the tensor gravitational waveforms, we forecast the constraints on a parameter characterizing the difference of scalar charges of NS-BH or NS-NS binaries for Advanced LIGO and Einstein Telescope. We illustrate how these constraints depend on redshift and signal-to-noise ratio, and on different possible priors. We show that in any case it is possible to constrain the scalar charge precisely, so that some scalarized NS solutions known in the literature can be excluded.

The majority of massive stars are born in binaries, and most unbind upon the first supernova. With precise proper motion surveys such as Gaia, it is possible to trace back the motion of stars in the vicinity of young remnants to search for ejected companions. Establishing the fraction of remnants with an ejected companion, and the photometric and kinematic properties of these stars, offers unique insight into supernova progenitor systems. In this paper, we employ binary population synthesis to produce kinematic and photometric predictions for ejected secondary stars. We demonstrate that the unbound neutron star velocity distribution from supernovae in binaries closely traces the input kicks. Therefore, the observed distribution of neutron star velocities should be representative of their natal kicks. We evaluate the probability for any given filter, magnitude limit, minimum measurable proper motion (as a function of magnitude), temporal baseline, distance and extinction that an unbound companion can be associated with a remnant. We compare our predictions with results from previous companion searches, and demonstrate that the current sample of stars ejected by the supernova of their companion can be increased by a factor of 5-10 with Gaia data release 3. Further progress in this area is achievable by leveraging the absolute astrometric precision of Gaia, and by obtaining multiple epochs of deep, high resolution near-infrared imaging with the Hubble Space Telescope, JWST and next-generation wide-field near-infrared observatories such as Euclid or the Nancy Grace Roman Space Telescope.

Funded by a $3M Department of Defense (DoD) National Defense Education Program (NDEP) award, we are developing and deploying on a national scale a follow-up curriculum to "Our Place In Space!", or OPIS!, in which approx. 3,500 survey-level astronomy students are using our global network of "Skynet" robotic telescopes each year. The goal of this new curriculum, called "Astrophotography of the Multi-Wavelength Universe!", or MWU!, is to boost the number of these students who choose STEM majors. During Y1, our participating educators have developed MWU!'s (now renumbered) 2nd and 4th modules, and are in the process of developing its 3rd and 7th modules (out of 7). Solid progress has also been made on the software front, (1) where we have developed new graphing/analysis/modeling interfaces in support of Modules 2 and 4, and in response to feedback from the participating educators; and (2) where we are in the process of developing and adding astrophotography capabilities to Afterglow Access (AgA), our student-level, web-based, image processing and analysis application, in support of Modules 1 - 3 and 5 - 7. On the hardware front, development of our first four signal-processing units proceeds on schedule; these are key to Skynet's integration of a global network of radio telescopes, capable of exploring the invisible universe. Preparations have also been made on the evaluation and accessibility fronts, for when the first MWU! modules are deployed in Spring 2023.

Sun-as-a-star coronal plasma composition, derived from full-Sun spectra, and the F10.7 radio flux (2.8 GHz) have been shown to be highly correlated (r = 0.88) during solar cycle 24. However, this correlation becomes nonlinear during increased solar magnetic activity. Here, we use co-temporal, high spatial resolution, multi-wavelength images of the Sun to investigate the underlying causes of the non-linearity between coronal composition (FIP bias) and F10.7 solar index correlation. Using the Karl G. Jansky Very Large Array (JVLA), Hinode/EIS (EUV Imaging Spectrometer), and the Solar Dynamic Observatory (SDO), we observed a small active region, AR 12759, throughout the solar atmosphere from the photosphere to the corona. Results of this study show that the magnetic field strength (flux density) in active regions plays an important role in the variability of coronal abundances, and it is likely the main contributing factor to this non-linearity during increased solar activity. Coronal abundances above cool sunspots are lower than in dispersed magnetic plage regions. Strong magnetic concentrations are associated with stronger F10.7 cm gyroresonance emission. Considering that as the solar cycle moves from minimum to maximum, the size of sunspots and their field strength increase with gyroresonance component, the distinctly different tendencies of radio emission and coronal abundances in the vicinity of sunspots is the likely cause of saturation of Sun-as-a-star coronal abundances during solar maximum, while the F10.7 index remains well correlated with the sunspot number and other magnetic field proxies.

The upcoming Laser Interferometer Space Antenna (LISA) is expected to detect gravitational waves (GWs) from massive black hole binaries (MBHB). Finding the electromagnetic (EM) counterparts for these GW events will be crucial for understanding how and where MBHBs merge, measuring their redshifts, constraining the Hubble constant and the graviton mass, and for other novel science applications. However, due to poor GW sky localisation, multi-wavelength, time-dependent electromagnetic (EM) models are needed to identify the right host galaxy among many candidates. We studied merging MBHBs embedded in a circumbinary disc using high-resolution two-dimensional simulations, with a $\Gamma$-law equation of state, incorporating viscous heating, shock heating, and radiative cooling. We simulate the binary from large separation until after merger, allowing us to model the decoupling of the binary from the circumbinary disc (CBD). We compute the EM signatures and identify distinct features before, during, and after the merger. Our main result is a multi-band EM signature: we find that the MBHB produces strong thermal X-ray emission until 1-2 days prior to the merger. However, as the binary decouples from the CBD, the X-ray-bright minidiscs rapidly shrink in size, become disrupted, and the accretion rate drops precipitously. As a result, the thermal X-ray luminosity drops by orders of magnitude, and the source remains X-ray dark for several days after the merger, regardless of any post-merger effects such as GW recoil or mass loss. Looking for the abrupt spectral change where the thermal X-ray disappears is a tell-tale EM signature of LISA mergers that does not require extensive pre-merger monitoring.

Feedback from active galactic nuclei (AGN) has long been invoked to explain the correlation between black hole mass and stellar velocity dispersion (M-{\sigma}) discovered in low redshift galaxies. We describe the time evolution of AGN in the M-{\sigma} plane based on our gap model (Garofalo, Evans & Sambruna 2010) for black hole accretion and jet formation illustrating a fundamental difference between jetted and non-jetted AGN. While the latter tend to evolve diagonally upward with black hole mass increasing along with stellar dispersion, we show that jetted AGN tend on average to move initially more upwards because their effect on velocity dispersion is weaker than for non-jetted AGN. But this initial phase is followed by a shift in the nature of the feedback, from positive to negative, a transition that is more dramatic on average in denser cluster environments. The feedback gets its kick from tilted jets which shut down star formation but increase velocity dispersion values. As this change in the nature of the feedback takes tens of million to hundreds of millions of years, jetted AGN triggered in mergers will evolve mostly upwards for up to order 10^8 years, followed by an extremely long phase in which low excitation progressively slows black hole growth but dramatically affects stellar dispersion. As a result, powerful jetted AGN evolve for most of their lives almost horizontally on the M-{\sigma} plane. The prediction is that strongest AGN feedback on stellar dispersion is a late universe phenomenon with M87 being a case in point. We show how jetted and non-jetted AGN parallel the Sersic and core-Sersic galaxy paths in the M-{\sigma} plane found by Sahu et al (2019).

Deep surveys of the CMB polarization have more information on the lensing signal than the quadratic estimators (QE) can capture. We showed in a recent work that a CMB lensing power spectrum built from a single optimized CMB lensing mass map, working in close analogy to state-of-the-art QE techniques, can result in an essentially optimal spectrum estimator at reasonable numerical cost. We extend this analysis here to account for real-life non-idealities including masking and realistic instrumental noise maps. As in the QE case, it is necessary to include small corrections to account for the estimator response to these anisotropies, which we demonstrate can be estimated easily from simulations. The realization-dependent debiasing of the spectrum remains robust, allowing unbiased recovery of the band powers even in cases where the statistical model used for the lensing map reconstruction is grossly wrong. This allows now robust and at the same time optimal CMB lensing constraints from CMB data, on all scales relevant for the inference of the neutrino mass, or other parameters of our cosmological model.

We report two low-frequency measurements of the power-law index for the amplitudes of giant radio pulses from the Crab pulsar. The two observations were taken with the Arecibo and Green Bank radio telescopes at center frequencies of 327 MHz and 350 MHz, respectively. We find best-fit values for the differential power-law index $\beta$ (where $dN/dS \propto S^\beta$ and $S$ is pulse amplitude) of $-2.63 \pm 0.05$ and $-3.6 \pm 0.5$ from the Arecibo and Green Bank data sets, respectively. Both values are broadly consistent with other values previously measured for the Crab pulsar at low radio frequencies. These reported values may be useful in future giant pulse studies of the Crab pulsar.

The binary systems consisting of a Be star and a white dwarf (BeWDs) are very interesting.They can originate from the binaries composed of a Be star and a subdwarf O or B star (BesdOBs), and they can merge into red giants via luminous red nova or can evolve into double WD potentially detected by $LISA$ mission. Using the method of population synthesis, we investigate the formation and the destiny of BeWDs,and discuss the effects of the metallicity ($Z$) and the common envelope evolution parameters. We find that BesdOBs are significant progenitors of BeWDs. About 30\% ($Z=0.0001$)-50\% ($Z=0.02$) of BeWDs come from BesdOBs. About 60\% ($Z=0.0001$) -70\% ($Z=0.02$) of BeWDs turn into red giants via a merger between a WD and a non-degenerated star. About 30\% ($Z=0.0001$) -40\% ($Z=0.02$) of BeWDs evolve into double WDs which are potential gravitational waves of $LISA$ mission at a frequency band between about $3\times10^{-3}$ and $3\times10^{-2}$ Hz. The common envelope evolution parameter introduces an uncertainty with a factor of about 1.3 on BeWD populations in our simulations.

Context: Solar flares are the result of the sudden release of magnetic energy in the corona. Much of this energy goes into accelerating charged particles to high velocity. These particles travel along the magnetic field and the energy is dissipated when the density gets high enough, primarily in the solar chromosphere. Modelling this region is difficult because the radiation energy balance is dominated by strong, optically thick spectral lines. Aims: Our aim is to provide the community with realistic simulations of a flaring loop with an emphasis on the detailed treatment of the chromospheric energy balance. This will enable a detailed comparison of existing and upcoming observations with synthetic observables from the simulations, thereby elucidating the complex interactions in a flaring chromosphere. Methods: We used the 1D radiation hydrodynamics code RADYN to perform simulations of the effect of a beam of electrons injected at the apex of a solar coronal loop. A grid of models was produced, varying the total energy input, the steepness, and low-energy cutoff of the beam energy spectrum. Results: The full simulation results for a grid of models are made available online. Some general properties of the simulations are discussed.

Geminga is an enigmatic radio-quiet gamma-ray pulsar located at a mere 250 pc distance from Earth. Extended very-high-energy gamma-ray emission around the pulsar was discovered by Milagro and later confirmed by HAWC, which are both water Cherenkov detector-based experiments. However, evidence for the Geminga pulsar wind nebula in gamma rays has long evaded detection by imaging atmospheric Cherenkov telescopes (IACTs) despite targeted observations. The detection of gamma-ray emission on angular scales > 2 deg poses a considerable challenge for the background estimation in IACT data analysis. With recent developments in understanding the complementary background estimation techniques of water Cherenkov and atmospheric Cherenkov instruments, the H.E.S.S. IACT array can now confirm the detection of highly extended gamma-ray emission around the Geminga pulsar with a radius of at least 3 deg in the energy range 0.5-40 TeV. We find no indications for statistically significant asymmetries or energy-dependent morphology. A flux normalisation of $(2.8\pm0.7)\times10^{-12}$ cm$^{-2}$s$^{-1}$TeV$^{-1}$ at 1 TeV is obtained within a 1 deg radius region around the pulsar. To investigate the particle transport within the halo of energetic leptons around the pulsar, we fitted an electron diffusion model to the data. The normalisation of the diffusion coefficient obtained of $D_0 = 7.6^{+1.5}_{-1.2} \times 10^{27}$ cm$^2$s$^{-1}$, at an electron energy of 100 TeV, is compatible with values previously reported for the pulsar halo around Geminga, which is considerably below the Galactic average.

There is evidence that exoplanet systems display intra-system uniformity in mass, radius, and orbital spacing (like "peas in a pod") when compared with the system-to-system variations of planetary systems. This has been interpreted as the outcome of the early stages of planet formation, indicative of a picture in which planets form at characteristic mass scales with uniform separations. In this paper, we argue instead that intra-system uniformity in planet sizes and orbital spacings likely arose from the long-term dynamical sculpting of initially-overly-packed planetary systems (in other words, the giant impact phase). With a suite of $N$-body simulations, we demonstrate that systems with random initial masses and compact planet spacings naturally develop intra-system uniformity, in quantitative agreement with observations, due to collisions between planets. Our results suggest that the pre-giant impact planet mass distribution is fairly wide and provide evidence for the prevalence of dynamical sculpting in shaping the observed population of exoplanets.

A network of cosmic strings (CS), if present, would continue emitting gravitational waves (GW) as it evolves throughout the history of the Universe. This results in a characteristic broad spectrum making it a perfect source to infer the expansion history. In particular, a short inflationary period caused by a supercooled phase transition would cause a drop in the spectrum at frequencies corresponding to that event. However, the impact on the spectrum is similar to the ones caused by an early matter-dominated era or from particle production, making it difficult to disentangle these different physical origins. We point out that, in the case of a short inflationary period, the GW spectrum receives an additional contribution from the phase transition itself. This leads to a characteristic imprint of a peak on top of a wide plateau both visible at future GW observatories.

We present a search for host galaxy associations for the third set of repeating fast radio burst (FRB) sources discovered by the CHIME/FRB Collaboration. Using the ~1 arcmin CHIME/FRB baseband localizations and probabilistic methods. We identify potential host galaxies of two FRBs, 20200223B and 20190110C at redshifts of 0.06024(2) and 0.12244(6), respectively. We also discuss the properties of a third marginal candidate host galaxy association for FRB 20191106C with a host redshift of 0.10775(1). The three putative host galaxies are all relatively massive, fall on the standard mass-metallicity relationship for nearby galaxies, and show evidence of ongoing star formation. They also all show signatures of being in a transitional regime, falling in the "green valley" which is between the bulk of star-forming and quiescent galaxies. The plausible host galaxies identified by our analysis are consistent with the overall population of repeating and non-repeating FRB hosts while increasing the fraction of massive and bright galaxies. Coupled with these previous host associations, we identify a possible excess of FRB repeaters whose host galaxies have M_u - M_r colors redder than the bulk of star-forming galaxies. Additional precise localizations are required to confirm this trend.

A successful measurement of the Stochastic Gravitational Wave Background (SGWB) in Pulsar Timing Arrays (PTAs) would open up a new window through which to test the predictions of General Relativity (GR). We consider how these measurements might reveal deviations from GR by studying the overlap reduction function -- the quantity that in GR is approximated by the Hellings-Downs curve -- in some sample modifications of gravity, focusing on the generic prediction of a modified dispersion relation for gravitational waves. We find a distinct signature of such modifications to GR -- a shift in the minimum angle of the angular distribution -- and demonstrate that this shift is quantitatively sensitive to any change in the phase velocity. In a given modification of gravity, this result can be used, in some regions of parameter space, to distinguish the effect of a modified dispersion relation from that due to the presence of extra polarization modes.

We compute the causality/positivity bounds on the Wilson coefficients of scalar-tensor effective field theories. Two-sided bounds are obtained by extracting IR information from UV physics via dispersion relations of scattering amplitudes, making use of the full crossing symmetry. The graviton $t$-channel pole is carefully treated in the numerical optimization, taking into account the constraints with fixed impact parameters. It is shown that the typical sizes of the Wilson coefficients can be estimated by simply inspecting the dispersion relations. We carve out sharp bounds on the leading coefficients, particularly, the scalar-Gauss-Bonnet couplings, and discuss how some bounds vary with the leading $(\partial\phi)^4$ coefficient and as well as phenomenological implications of the causality bounds.

We present a gravitoelectric quadrupolar dynamical tidal-interaction Hamiltonian for a compact binary system, that is valid to second order in the post-Newtonian expansion. Our derivation uses the diagrammatic effective field theory approach, and involves Feynman integrals up to two loops, evaluated with the dimensional regularization scheme. We also derive the effective Hamiltonian for adiabatic tides, obtained by taking the appropriate limit of the dynamical effective Hamiltonian, and we check its validity by verifying the complete Poincar\'e algebra. In the adiabatic limit, we also calculate two gauge-invariant observables, namely, the binding energy for a circular orbit and the scattering angle in a hyperbolic scattering. Our results are important for developing accurate gravitational waveform models for neutron-star binaries for present and future gravitational-wave observatories.

The current and upcoming generations of gravitational wave experiments represent an exciting step forward in terms of detector sensitivity and performance. For example, key upgrades at the LIGO, Virgo and KAGRA facilities will see the next observing run (O4) probe a spatial volume around four times larger than the previous run (O3), and design implementations for e.g. the Einstein Telescope, Cosmic Explorer and LISA experiments are taking shape to explore a wider frequency range and probe cosmic distances. In this context, however, a number of very real data analysis problems face the gravitational wave community. For example, it will be crucial to develop tools and strategies to analyse (amongst other scenarios) signals that arrive coincidentally in detectors, longer signals that are in the presence of non-stationary noise or other shorter transients, as well as noisy, potentially correlated, coherent stochastic backgrounds. With these challenges in mind, we develop peregrine, a new sequential simulation-based inference approach designed to study broad classes of gravitational wave signal. In this work, we describe the method and implementation, before demonstrating its accuracy and robustness through direct comparison with established likelihood-based methods. Specifically, we show that we are able to fully reconstruct the posterior distributions for every parameter of a spinning, precessing compact binary coalescence using one of the most physically detailed and computationally expensive waveform approximants (SEOBNRv4PHM). Crucially, we are able to do this using only 2\% of the waveform evaluations that are required in e.g. nested sampling approaches. Finally, we provide some outlook as to how this level of simulation efficiency and flexibility in the statistical analysis could allow peregrine to tackle these current and future gravitational wave data analysis problems.

The two-body problem is extensively studied in open systems and asymptotically flat spacetimes. However, there are many systems where radiation is trapped: they range from radiating charges in cavities to low-energy excitations of massive degrees of freedom, to anti-de Sitter spacetimes. Here, we study the problem of motion of a pointlike particle orbiting a massive compact object inside a cavity. We first show that - assuming circular motion - there are initial conditions for which the self-force vanishes and the binary is eternal. We then consider the evolution of the system under radiation reaction in a toy model which we argue captures the essentials of orbiting particles. We show that eternal circular binaries may exist. We also show that the presence of cavity modes leads to chaos in regimes of strong coupling or when the system is initialized close enough to a resonance. Our results have implications for physics in anti-de Sitter spacetimes and possibly for binaries evolving within dark matter haloes, if it consists on massive fundamental fields.

In this paper, we constrain the linear dark-matter-related parameter of a static spherically symmetric f (R) black hole spacetime regarding the observed angular diameters of M87* and Sgr A* from the EHT. We then investigate the light deflection angles inferred from direct analytical calculation of null geodesics and that obtained from the Gauss-Bonnet theorem. Assuming an optically thin accretion disk for the black hole and after discussing its properties, we conceive different emission profiles and investigate the shadow cast of this black hole when it is illuminated by the disk. Furthermore, we simulate the brightness of an infalling spherical accretion in the context of the silhouette imaging of the black hole. We find that, excluding some specific cases, the specific observed brightness of the accretion disk consists of the direct emission, rather than that for the lensing and photon rings. Furthermore, it is revealed that the linear dark parameter of the black hole has considerable effects on the size of the shadow and its brightness. The discussion is done both analytically and numerically, and ray-tracing methods are employed to generate proper visualizations

We describe the gravitational lens on a constant-curvature background by using an optical constant-curvature (OCC) approach that allows to explicitly take into account a global geometry of the background space. First, light rays are curves in the space described by an optical metric. The OCC approach focuses on the case that the optical metric for the background spacetime has a constant curvature, for which the exact lens equation on an OCC background [Phys. Rev. D 105, 084022 (2022)] can be used. As a concrete example, next we discuss the gravitational lens in Mannheim-Kazanas (MK) solution, which include Rindler and de Sitter terms. By fully taking into account a background dependence of the light deflection, the deflection angle of light consistent with the OCC approach is well defined at large distance. In the OCC approach, finally, we examine the global behavior of the deflection angle and the gravitational lens observables in the Weyl gravity.

Gravitational-wave black-hole spectroscopy provides a unique opportunity to test the strong-field regime of gravity and the nature of the final object formed in the aftermath of a merger. Here we investigate the prospects for black-hole spectroscopy with third-generation gravitational-wave detectors, in particular the Einstein Telescope in different configurations, possibly in combination with Cosmic Explorer. Using a state-of-the-art population model for stellar-origin binary black holes informed by LIGO-Virgo-KAGRA data, we compute the average number of expected events for precision black-hole spectroscopy using a Fisher-matrix analysis. We find that Einstein Telescope will measure two independent quasinormal modes within ${\cal O}(1)\%$ (resp. ${\cal O}(10)\%$) relative uncertainty for at least ${\cal O}(1)$ (resp. ${\cal O}(500)$) events per year, with similar performances in the case of a single triangular configuration or two L-shaped detectors with same arm length. A 15-km arm-length configuration would improve rates by roughly a factor of two relative to a 10-km arm-length configuration. When operating in synergy with Cosmic Explorer the rates will improve significantly, reaching few-percent accuracy for ${\cal O}(100)$ events per year.

The perspectives of detecting the general relativistic Lense-Thirring effect on the orbits of the Galilean moons of Jupiter induced by its angular momentum ${\boldsymbol{S}}$ are preliminarily investigated. Numerical integrations over one century show that the expected gravitomagnetic signatures of some observables such as the satellites' right ascension $\alpha$ and declination $\delta$ are as large as tens of arcseconds for Io, while for Callisto they drop to the $\simeq 0.2\,\mathrm{arcseconds}$ level; the shifts in the transverse component $T$ of the orbit range from 40 km for Io to 2 km for Callisto. Major competing effects due to the mismodeling in the zonal multipoles $J_\ell,\,\ell=2,\,3,\,4,\,\ldots$ of the Jovian non-spherically symmetric gravity field and in the Jupiter's spin axis ${\boldsymbol{\hat{k}}}$ should have a limited impact, especially in view of the future improvements in determing such parameters expected after the completion of the ongoing Juno mission in the next few years. Present-day accuracy in knowing the orbits of Io, Europa, Ganymede and Callisto is of the order of 10 milliarcseconds, to be likely further improved by the approved JUICE and Clipper missions. This suggests that the Lense-Thirring effect in the main Jovian system of moons might be detectable with dedicated data reductions in which the gravitomagnetic field is explicitly modeled and solved-for.

We study the properties of the stochastic gravitational wave background (SGWB) produced by domain walls (DWs) during inflation without forming a network. We numerically simulate the DW production caused by a second-order phase transition and calculate the SGWB spectrum using a $1000\times1000\times1000$ lattice. We show that the SGWB can be observed directly by future terrestrial and spatial gravitational wave detectors and through the B-mode spectrum in CMB. This signal can also explain the common noise process observed by pulsar timing array experiments. With numerical simulations, we derive an empirical formula for the strength and qualitative features of the SGWB spectrum. The details of the SGWB spectrum also contain information about the later evolution of the universe.

We formulate a version of the low-scale Majoron model equipped with an inverse seesaw mechanism featuring lepton-number preserving dimension-6 operators in the scalar potential. Contrary to its dimension-4 counterpart, we find that the model can simultaneously provide light and ultralight Majorons, neutrino masses and their mixing, while featuring strong first-order cosmological phase transitions associated to the spontaneous breaking of the lepton number and the electroweak symmetries in the early Universe. We show by a detailed numerical analysis that under certain conditions on the parameter space accounted for in collider physics, the model can be probed via the primordial gravitational wave spectrum potentially observable at LISA and other planned facilities.

We study consistently the effects of magnetic field on hot and dense matter. In particular, we look for differences that arise due to assumptions that reproduce the conditions produced in particle collisions or astrophysical scenarios, such as in the core of fully evolved neutron stars. We assume the magnetic field to be either constant or follow a profile extracted from general relativity calculations of magnetars and make use of two realistic models that can consistently describe chiral symmetry restoration and deconfinement to quark matter, {the CMF and the PNJL models}. We find that net isospin, strangeness, and {weak} chemical equilibrium with leptons can considerably change the effects of temperature and magnetic fields on particle content and deconfinement in dense matter. We finish by discussing the possibility of experimentally detecting quark deconfinement in dense and/or hot matter and the possible role played by magnetic fields.

This white paper outlines the plans of the History Philosophy Culture Working Group of the Next Generation Event Horizon Telescope Collaboration.

There are several unanswered questions regarding neutrinos which pave the way for physics beyond the standard model (SM) of particle physics. Generalized interactions of neutrinos provide a way to characterize these effects in a manner which is even more general than the oft-studied non-standard neutrino interactions. These interactions are described by higher dimensional operators maintaining the SM gauge symmetries. On the other hand cosmic neutrino background, although yet to be detected directly, is a robust prediction of the SM and the standard cosmology. We perform a global analysis of the relevant generalized neutrino interactions which are expressly relevant for the proposed cosmic neutrino detector PTOLEMY. The electron spectrum due to the capture of cosmic neutrinos on radioactive tritium gets modified due to the presence of these generalized interactions. We also show how the differential electron spectrum is sensitive to the finite experimental resolution, mass of the lightest neutrino eigenstate, the strength of these interactions and the ordering of neutrino mass.

We study stochastic gravitational waves from cosmic strings generated in an ultraviolet-complete model for pseudo-Nambu-Goldstone dark matter with a hidden $\mathrm{U(1)}$ gauge symmetry. The dark matter candidate in this model can naturally evade direct detection bounds and easily satisfy other phenomenological constraints. The bound on the dark matter lifetime implies an ultraviolet scale higher than $10^9~\mathrm{GeV}$. The spontaneous $\mathrm{U(1)}$ symmetry breaking at such a high scale would induce cosmic strings with high tension, resulting in a stochastic gravitational wave background with a high energy density. We investigate the constraints from current gravitational wave experiments as well as the future sensitivity. We find that most of the viable parameter points could be well studied in future gravitational wave experiments.

We make use of Borel resummation to extract the exact time dependence from the divergent series found in the context of stochastic inflation. Correlation functions of self-interacting scalar fields in de Sitter spacetime are known to develop secular IR divergences via loops, and the first terms of the divergent series have been consistently computed both with standard techniques for curved spacetime quantum field theory and within the framework of stochastic inflation. We show that Borel resummation can be used to interpret the divergent series and to correctly infer the time evolution of the correlation functions. In practice, we adopt a method called Borel--Pad\'{e} resummation where we approximate the Borel transformation by a Pad\'{e} approximant. We also discuss the singularity structures of Borel transformations and mention possible applications to cosmology.

We have studied the role of Hawking evaporation of primordial black hole in the production of the supersymmetric particles like sneutrinos - the super-partner of heavy right handed neutrinos. Considering lepton number violating decays of such particles and $CP$ violating phases due to soft supersymmetry breaking terms, we have obtained the baryonic asymmetry of the universe which depends on the mass of primordial black holes. Apart from CMB and BBN constraints on such mass, we have shown more stringent upper bound on this mass from the requirement of black hole evaporation temperature to be above the temperature required for almost first order phase transition so that leptogenesis could create the observed baryonic asymmetry. We have shown how the primordial black hole mass, heavy right handed neutrino mass and soft supersymmetry breaking parameters are related from the requirement of successful leptogenesis and the allowed parameter space to avoid gravitino problem. Considering experimental constraint on branching ratio of $\mu \rightarrow e \gamma $, we have shown the connection of right handed neutrino mass scale with the lower bound of the typical mass scale of supersymmetric particles.