Papers for Monday, Apr 19 2021

We describe the processing of the PHANGS-ALMA survey and present the PHANGS-ALMA pipeline, a public software package that processes calibrated interferometric and total power data into science-ready data products. PHANGS-ALMA is a large, high-resolution survey of CO J=2-1 emission from nearby galaxies. The observations combine ALMA's main 12-m array, the 7-m array, and total power observations and use mosaics of dozens to hundreds of individual pointings. We describe the processing of the u-v data, imaging and deconvolution, linear mosaicking, combining interferometer and total power data, noise estimation, masking, data product creation, and quality assurance. Our pipeline has a general design and can also be applied to VLA and ALMA observations of other spectral lines and continuum emission. We highlight our recipe for deconvolution of complex spectral line observations, which combines multiscale clean, single scale clean, and automatic mask generation in a way that appears robust and effective. We also emphasize our two-track approach to masking and data product creation. We construct one set of "broadly masked" data products, which have high completeness but significant contamination by noise, and another set of "strictly masked" data products, which have high confidence but exclude faint, low signal-to-noise emission. Our quality assurance tests, supported by simulations, demonstrate that 12-m+7-m deconvolved data recover a total flux that is significantly closer to the total power flux than the 7-m deconvolved data alone. In the appendices, we measure the stability of the ALMA total power calibration in PHANGS--ALMA and test the performance of popular short-spacing correction algorithms.

The Constellation Project's planned return to the moon requires numerous landings at the same site. Since the top few centimeters are loosely packed regolith, plume impingement from the Lander ejects the granular material at high velocities. Much work is needed to understand the physics of plume impingement during landing in order to protect hardware surrounding the landing sites. While mostly qualitative in nature, the Apollo Lunar Module landing videos can provide a wealth of quantitative information using modern photogrammetry techniques. The authors have used the digitized videos to quantify plume impingement effects of the landing exhaust on the lunar surface. The dust ejection angle from the plume is estimated at 1-3 degrees. The lofted particle density is estimated at 10^8 - 10^13 particles/m^3. Additionally, evidence for ejection of large 10-15 cm sized objects and a dependence of ejection angle on thrust are presented. Further work is ongoing to continue quantitative analysis of the landing videos.

We present a scenario for the formation and the morphology of the arrow-shaped dwarf irregular galaxy NGC 1427A in the Fornax Cluster. This galaxy shows intriguing stellar and gaseous tails pointing in different directions for which alternative but not conclusive formation scenarios have been proposed in the literature. We performed N-body/SPH simulations of dwarf galaxies falling into a model of the Fornax cluster, exhibiting a jellyfish-like appearance while undergoing ram-pressure stripping. We noted that some of our models show interesting tail morphologies similar to that of NGC 1427A. In this way, the peculiar NGC 1427A structure can be studied using models whose stellar and neutral gas photometry and kinematics are in good agreement with the observed ones, without the need of invoking an interaction with a nearby galaxy. Thanks to the tails, we can identify the requirements for a galaxy to expose such a structure and assess the possible position and velocity of the galaxy in the cluster. This puts constraints on the orbit of the galaxy, its position in the cluster and the time since its pericentre passage. From the statistics of identified snapshots following our modelling, we found that the most likely position of the galaxy is around 200 kpc in front of the cluster centre, travelling towards the cluster with a velocity angle with respect to the line-of-sight direction of around 50 deg. This analysis can be useful in future observations of similar galaxies in clusters to characterise their position and velocity in the cluster and their formation.

The nature of dark matter (DM) is unknown. One compelling possibility is DM being composed of primordial black holes (PBHs), given the tight limits on some types of elementary particles as DM. Here, we show that Neptune's orbital eccentricity places constraints on PBHs in the sub-solar mass window. We also derive new limits on the local abundance of rogue interstellar planets.

All hydrodynamical simulations of turbulent astrophysical phenomena require sub-grid scale models to properly treat energy dissipation and metal mixing. We present the first implementation and application of an anisotropic eddy viscosity and metal mixing model in Lagrangian astrophysical simulations, including a dynamic procedure for the model parameter. We compare these two models directly to the common Smagorinsky and dynamic variant. Using the mesh-free finite mass method as an example, we show that the anisotropic model is best able to reproduce the proper Kolmogorov inertial range scaling in homogeneous, isotropic turbulence. Additionally, we provide a method to calibrate the metal mixing rate that ensures numerical convergence. In our first application to cosmological simulations, we find that all models strongly impact the early evolution of galaxies leading to differences in enrichment and thermodynamic histories. The anisotropic model has the strongest impact, with little difference between the dynamic and the constant-coefficient variant. We also find that the metal distribution functions in the circumgalactic gas are significantly tighter at all redshifts, with the anisotropic model providing the tightest distributions. This is contrary to a recent study that found metal mixing to be relatively unimportant on cosmological scales. In all of our experiments the constant-coefficient Smagorinsky and anisotropic models rivaled their dynamic counterparts, suggesting that the computationally inexpensive constant-coefficient models are viable alternatives in cosmological contexts.

Despite evidence for the existence of dark matter (DM) from very high and low redshifts, a moderate amount of DM particle decay remains a valid possibility. This includes both models with very long-lived yet unstable particles or mixed scenarios where only a small fraction of dark matter is allowed to decay. In this paper, we investigate how DM particles decaying into radiation affect non-linear structure formation. We look at the power spectrum and its redshift evolution, varying both the decay lifetime ($\tau$) and the fraction of decaying-to-total dark matter ($f$), and we propose a fitting function that reaches sub-percent precision below $k\sim10$ h/Mpc. Based on this fit, we perform a forecast analysis for a Euclid-like weak lensing (WL) survey, including both massive neutrino and baryonic feedback parameters. We find that with WL observations alone, it is possible to rule out decay lifetimes smaller than $\tau=75$ Gyr (at 95 percent CL) for the case that all DM is unstable. This constraint improves to $\tau=182$ Gyr if the WL data is combined with CMB priors from the Planck satellite and to $\tau=275$ Gyr if we further assume baryonic feedback to be fully constrained by upcoming Sunyaev-Zeldovich or X-ray data. The latter shows a factor of 3.2 improvement compared to constraints from CMB data alone. Regarding the scenario of a strongly decaying sub-component of dark matter with $\tau\sim 30$ Gyr or lower, it will be possible to rule out a decaying-to-total fraction of $f>0.49$, $f>0.21$, and $f>0.13$ (at the 95 percent CL) for the same three scenarios. We conclude that the upcoming stage-IV WL surveys will allow us to significantly improve current constraints on the stability of the dark matter sector.

We present the first study of the large-scale clustering of post-starburst (PSB) galaxies in the high redshift Universe ($0.5<z<3.0$). We select $\sim4000$ PSB galaxies photometrically, the largest high-redshift sample of this kind, from two deep large-scale near-infrared surveys: the UKIDSS Ultra Deep Survey (UDS) DR11 and the Cosmic Evolution Survey (COSMOS). Using angular cross-correlation techniques, we estimate the halo masses for this large sample of PSB galaxies and compare them with quiescent and star-forming galaxies selected in the same fields. We find that low-mass, low-redshift ($0.5<z<1.0$) PSB galaxies preferentially reside in very high-mass dark matter halos (M$_{\text{halo}}>10^{14}$M$_{\odot}$), suggesting they are likely to be infalling satellite galaxies in cluster-like environments. High-mass PSB galaxies are more weakly clustered at low redshifts, but they reside in higher mass haloes with increasing look-back time, suggesting strong redshift-dependent halo downsizing. These key results are consistent with previous results suggesting that two main channels are responsible for the rapid quenching of galaxies. While high-redshift ($z>1$) galaxies appear to be quenched by secular feedback mechanisms, processes associated with dense environments are likely to be the key driver of rapid quenching in the low-redshift Universe ($z<1$). Finally, we show that the clustering of photometrically selected PSBs are consistent with them being direct descendants of highly dust-enshrouded sub-millimetre galaxies (SMGs), providing tantalising evidence for the oft-speculated evolutionary pathway from starburst to quiescence.

We present the analysis of X-ray and optical observations of gas filaments observed in the radio source 3CR 318.1, associated with NGC 5920, the Brightest Cluster Galaxy (BCG) of MKW 3s, a nearby cool core galaxy cluster. This work is one of the first X-ray and optical analyses of filaments in cool core clusters carried out using MUSE observations. We aim at identifying the main excitation processes responsible for the emission arising from these filaments. We complemented the optical VLT/MUSE observations, tracing the colder gas phase, with X-ray $\textit{Chandra}$ observations of the hotter highly ionized gas phase. Using the MUSE observations, we studied the emission line intensity ratios along the filaments to constrain the physical processes driving the excitation, and, using the $\textit{Chandra}$ observations, we carried out a spectral analysis of the gas along these filaments. We found a spatial association between the X-ray and optical morphology of these filaments, which are colder and have lower metal abundance than the surrounding intra-cluster medium (ICM), as already seen in other BCGs. Comparing with previous results from the literature for other BCGs, we propose that the excitation process that is most likely responsible for these filaments emission is a combination of star formation and shocks, with a likely contribution from self-ionizing, cooling ICM. Additionally, we conclude that the filaments most likely originated from AGN-driven outflows in the direction of the radio jet.

We investigate the formation and growth of massive black hole (BH) seeds in dusty star-forming galaxies, relying and extending the framework proposed by Boco et al. 2020. Specifically, the latter envisages the migration of stellar compact remnants (neutron stars and stellar-mass black holes) via gaseous dynamical friction towards the galaxy nuclear region, and their subsequent merging to grow a massive central BH seed. In this paper we add two relevant ingredients: (i) we include primordial BHs, that could constitute a fraction $f_{\rm pBH}$ of the dark matter, as an additional component participating in the seed growth; (ii) we predict the stochastic gravitational wave background originated during the seed growth, both from stellar compact remnant and from primordial BH mergers. We find that the latter events contribute most to the initial growth of the central seed during a timescale of $10^6-10^7\,\rm yr$, before stellar compact remnant mergers and gas accretion take over. In addition, if the fraction of primordial BHs $f_{\rm pBH}$ is large enough, gravitational waves emitted by their mergers in the nuclear galactic regions could be detected by future interferometers like Einsten Telescope, DECIGO and LISA. As for the associated stochastic gravitational wave background, we predict that it extends over the wide frequency band $10^{-6}\lesssim f [{\rm Hz}]\lesssim 10$, which is very different from the typical range originated by mergers of isolated binary compact objects. On the one hand, the detection of such a background could be a smoking gun to test the proposed seed growth mechanism; on the other hand, it constitutes a relevant contaminant from astrophysical sources to be characterized and subtracted, in the challenging search for a primordial background of cosmological origin.

Determining the level of chemical complexity within dense starless and gravitationally bound prestellar cores is crucial for constructing chemical models, which subsequently constrain the initial chemical conditions of star formation. We have searched for complex organic molecules (COMs) in the young starless core L1521E, and report the first clear detection of dimethyl ether (CH$_3$OCH$_3$), methyl formate (HCOOCH$_3$), and vinyl cyanide (CH$_2$CHCN). Eight transitions of acetaldehyde (CH$_3$CHO) were also detected, five of which (A states) were used to determine an excitation temperature to then calculate column densities for the other oxygen-bearing COMs. If source size was not taken into account (i.e., if filling fraction was assumed to be one), column density was underestimated, and thus we stress the need for higher resolution mapping data. We calculated L1521E COM abundances and compared them to other stages of low-mass star formation, also finding similarities to other starless/prestellar cores, suggesting related chemical evolution. The scenario that assumes formation of COMs in gas-phase reactions between precursors formed on grains and then ejected to the cold gas via reactive desorption was tested and was unable to reproduce observed COM abundances, with the exception of CH$_3$CHO. These results suggest that COMs observed in cold gas are formed not by gas-phase reactions alone, but also through surface reactions on interstellar grains. Our observations present a new, unique challenge for existing theoretical astrochemical models.

Many black holes (BHs) detected by the Laser Interferometer Gravitational-wave Observatory (LIGO) and the Virgo detectors are multiple times more massive than those in X-ray binaries. One possibility is that some BBHs merge within a few Schwarzschild radii of a supermassive black hole (SMBH), such that the gravitational waves (GWs) are highly redshifted, causing the mass inferred from GW signals to appear higher than the real mass. The difficulty of this scenario lies in the delivery of BBH to such a small distance to a SMBH. Here we revisit the theoretical models for the migration of compact objects (COs) in the accretion discs of active galactic nuclei (AGNs). We find that when the accretion rate is high so that the disc is best described by the slim disc model, the COs in the disc could migrate to a radius close to the innermost stable circular orbit (ISCO) and be trapped there for the remaining lifetime of the AGN. The exact trapping radius coincides with the transition region between the sub- and super-Keplerian rotation of the slim disc. We call this region "the last migration trap" because inside it COs can no longer be trapped for a long time. We pinpoint the parameter space which could induce such a trap and we estimate that the last migration trap contributes a few per cent of the LIGO/Virgo events. Our result implies that a couple of BBHs discovered by LIGO/Virgo could have smaller intrinsic masses.

Astronomy has always been propelled by the discovery of new phenomena lacking precedent, often followed by new theories to explain their existence and properties. In the modern era of large surveys tiling the sky at ever high precision and sampling rates, these serendipitous discoveries look set to continue, with recent examples including Boyajian's Star, Fast Radio Bursts and `Oumuamua. Accordingly, we here look ahead and aim to provide a statistical framework for interpreting such events and providing guidance to future observations, under the basic premise that the phenomenon in question stochastically repeat at some unknown, constant rate, $\lambda$. Specifically, expressions are derived for 1) the a-posteriori distribution for $\lambda$, 2) the a-posteriori distribution for the recurrence time, and, 3) the benefit-to-cost ratio of further observations relative to that of the inaugural event. Some rule-of-thumb results for each of these are found to be 1) $\lambda < \{0.7, 2.3, 4.6\}\,t_1^{-1}$ to $\{50, 90, 95\}\%$ confidence (where $t_1=$ time to obtain the first detection), 2) the recurrence time is $t_2 < \{1, 9, 99\}\,t_1$ to $\{50, 90, 95\}\%$ confidence, with a lack of repetition by time $t_2$ yielding a $p$-value of $1/[1+(t_2/t_1)]$, and, 3) follow-up for $\lesssim 10\,t_1$ is expected to be scientifically worthwhile under an array of differing assumptions about the object's intrinsic scientific value. We apply these methods to the Breakthrough Listen Candidate 1 signal and tidal disruption events observed by TESS.

We study the contributions to the relativistic Fe $K_{\alpha}$ line profile from higher order images (HOIs) produced by strongly deflected rays from the disk which cross the plunging region, located between the innermost stable circular orbit (ISCO) radius and the event horizon of a Kerr black hole. We investigate the characteristics features imprinted by the HOIs in the line profile for different black hole spins, disk emissivity laws and inclinations. We find that they extend from the red wing of the profile up to energies slightly lower than those of the blue peak, adding $\sim 0.4 - 1.3$\% to the total line flux. The contribution to the specific flux is often in the $\sim 1$\% to 7\% range, with the highest values attained for low and negative spin ($a\lesssim 0.3$) black holes surrounded by intermediate inclination angle ($i\sim40^{\circ}$) disks. We simulate future observations of a black hole X-ray binary system with the Large Area Detector of the planned X-ray astronomy \emph{enhanced X-ray Timing and Polarimetry Mission} (eXTP) and find that the \fekal\ of systems accreting at $\lesssim 1 $\% the Eddington rate are affected by the HOI features for a range of parameters. This would provide evidence of the extreme gravitational lensing of HOI rays. Our simulations show also that not accounting for HOI contributions to the Fe $K_{\alpha}$ line profile may systematically bias measurements of the black hole spin parameter towards values higher by up to $\sim 0.3$ than the inputted ones.

LS 5039 is an enigmatic high-mass gamma-ray binary which hosts a powerful O6.5V companion, while the nature of the compact object is still to be established by multi-wavelength observations analysis. Phase-resolved multi-instrument spectra of nonthermal emission from LS 5039 are analyzed in order to produce reliable spectral models, which can be further employed to select between various scenarios and theoretical models of the binary. The combined phase-resolved hard X-ray and MeV-range gamma-ray spectra obtained with XMM-Newton, Suzaku, NuSTAR, INTEGRAL, and Comptel indicate a meaningful spectral hardening above 50 keV. The spectral break observed in both major phases of the binary may indicate the presence of two populations of accelerated leptons in the emitting region, which could originate from the collision of the O6.5V companion star wind with the flow from yet unidentified compact object.

We present PHANGS-ALMA, the first survey to map CO J=2-1 line emission at ~1" ~ 100pc spatial resolution from a representative sample of 90 nearby (d<~20 Mpc) galaxies that lie on or near the z=0 "main sequence" of star-forming galaxies. CO line emission traces the bulk distribution of molecular gas, which is the cold, star-forming phase of the interstellar medium. At the resolution achieved by PHANGS-ALMA, each beam reaches the size of a typical individual giant molecular cloud (GMC), so that these data can be used to measure the demographics, life-cycle, and physical state of molecular clouds across the population of galaxies where the majority of stars form at z=0. This paper describes the scientific motivation and background for the survey, sample selection, global properties of the targets, ALMA observations, and characteristics of the delivered ALMA data and derived data products. As the ALMA sample serves as the parent sample for parallel surveys with VLT/MUSE, HST, AstroSat, VLA, and other facilities, we include a detailed discussion of the sample selection. We detail the estimation of galaxy mass, size, star formation rate, CO luminosity, and other properties, compare estimates using different systems and provide best-estimate integrated measurements for each target. We also report the design and execution of the ALMA observations, which combine a Cycle~5 Large Program, a series of smaller programs, and archival observations. Finally, we present the first 1" resolution atlas of CO emission from nearby galaxies and describe the properties and contents of the first PHANGS-ALMA public data release.

This series of papers aims at understanding the formation and evolution of non-barred disc galaxies. We use the new spectro-photometric decomposition code, C2D, to separate the spectral information of bulges and discs of a statistically representative sample of galaxies from the CALIFA survey. Then, we study their stellar population properties analising the structure-independent datacubes with the Pipe3D algorithm. We find a correlation between the bulge-to-total ($B/T$) luminosity (and mass) ratio and galaxy stellar mass. The $B/T$ mass ratio has only a mild evolution with redshift, but the bulge-to-disc ($B/D$) mass ratio shows a clear increase of the disc component since redshift $z < 1$ for massive galaxies. The mass-size relation for both bulges and discs describes an upturn at high galaxy stellar masses (log{(M_{\star}/M_{\sun})} > 10.5). The relation holds for bulges but not for discs when using their individual stellar masses. We find a negligible evolution of the mass-size relation for both the most massive (log{(M_{\star \rm ,b,d}/M_{\sun})} > 10) bulges and discs. For lower masses, discs show a larger variation than bulges. We also find a correlation between the S\'ersic index of bulges and both galaxy and bulge stellar mass, which does not hold for the disc mass. Our results support an inside-out formation of nearby non-barred galaxies, and they suggest that i) bulges formed early-on and ii) they have not evolved much through cosmic time. However, we find that the early properties of bulges drive the future evolution of the galaxy as a whole, and particularly the properties of the discs that eventually form around them.

The plasma contributing to emission from the Sun between the cool chromosphere ($\le 10^4$K) and hot corona ($\ge 10^6$K) has been subjected to many different interpretations. Here we look at the magnetic structure of this transition region (TR) plasma, based upon the implications of CLASP2 data of an active region recently published by Ishikawa et al., and earlier IRIS and SDO data of quiet regions. Ishikawa et al. found that large areas of sunspot plages are magnetically unipolar as measured in the cores of \ion{Mg}{2} resonance lines, formed in the lower transition region under low plasma-$\beta$ conditions. Here we show that IRIS images in the line cores have fibrils which well aligned with the overlying coronal loop segments seen in the 171 \AA{} channel of SDO. When the TR emission in active regions arise from plasma magnetically and thermally connected to the corona, then the line cores can provide the first credible magnetic boundary conditions for force-free calculations extended to the corona. We also re-examine IRIS images of dynamic TR cool loops previously reported as a major contributor to transition region emission from the quiet Sun. Dynamic cool loops contribute only a small fraction of the total TR emission from the quiet Sun.

This paper details speckle observations of binary stars taken at the Lowell Discovery Telescope, the WIYN Telescope, and the Gemini telescopes between 2016 January and 2019 September. The observations taken at Gemini and Lowell were done with the Differential Speckle Survey Instrument (DSSI), and those done at WIYN were taken with the successor instrument to DSSI at that site, the NN-EXPLORE Exoplanet Star and Speckle Imager (NESSI). In total, we present 378 observations of 178 systems and we show that the uncertainty in the measurement precision for the combined data set is ~2 mas in separation, ~1-2 degrees in position angle depending on the separation, and $\sim$0.1 magnitudes in magnitude difference. Together with data already in the literature, these new results permit 25 visual orbits and one spectroscopic-visual orbit to be calculated for the first time. In the case of the spectroscopic-visual analysis, which is done on the trinary star HD 173093, we calculate masses with precision of better than 1% for all three stars in that system. Twenty-one of the visual orbits calculated have a K dwarf as the primary star; we add these to the known orbits of K dwarf primary stars and discuss the basic orbital properties of these stars at this stage. Although incomplete, the data that exist so far indicate that binaries with K dwarf primaries tend not to have low-eccentricity orbits at separations of one to a few tens of AU, that is, on solar-system scales.

Various remote sensing observations have been used so far to probe the turbulent properties of the solar wind. Using the recently reported density modulation indices that are derived using angular broadening observations of Crab Nebula during 1952 - 2013, we measured the solar wind proton heating using the kinetic $\rm Alfv\acute{e}n$ wave dispersion equation. The estimated heating rates vary from $\approx 1.58 \times 10^{-14}$ to $1.01 \times 10^{-8} ~\rm erg~ cm^{-3}~ s^{-1}$ in the heliocentric distance range 5 - 45 $\rm R_{\odot}$. Further, we found that heating rates vary with the solar cycle in correlation with density modulation indices. The models derived using in-situ measurements (for example, electron/proton density, temperature, and magnetic field) that the recently launched Parker Solar Probe observes (planned closest perihelia $\rm 9.86~ R_{\odot}$ from the center of the Sun) are useful in the estimation of the turbulent heating rate precisely. Further, we compared our heating rate estimates with the one derived using previously reported remote sensing and in-situ observations.

We present 0\farcs{14}-resolution Atacama Large Millimeter/submillimeter Array (ALMA) CO(2$-$1) observations of the circumnuclear gas disk in UGC 2698, a local compact galaxy. The disk exhibits regular rotation with projected velocities rising to 450 km s$^{-1}$ near the galaxy center. We fit gas-dynamical models to the ALMA data cube, assuming the CO emission originates from a dynamically cold, thin disk, and measured the mass of the supermassive black hole (BH) in UGC 2698 to be $M_{\mathrm{BH}} = (2.46 \pm{0.07}$ [$1\sigma$ stat] $^{+0.70}_{-0.78}$ [sys])$\times 10^9$ $M_\odot$. UGC 2698 is part of a sample of nearby early-type galaxies that are plausible $z\sim2$ red nugget relics. Previous stellar-dynamical modeling for three galaxies in the sample found BH masses consistent with the BH mass$-$stellar velocity dispersion ($M_{\mathrm{BH}}-\sigma_\star$) relation but over-massive relative to the BH mass$-$bulge luminosity ($M_{\mathrm{BH}}-L_{\mathrm{bul}}$) correlation, suggesting that BHs may gain the majority of their mass before their host galaxies. However, UGC 2698 is consistent with both $M_{\mathrm{BH}}-\sigma_\star$ and $M_{\mathrm{BH}}-L_{\mathrm{bul}}$. As UGC 2698 has the largest stellar mass and effective radius in the local compact galaxy sample, it may have undergone more recent mergers that brought it in line with the BH scaling relations. Alternatively, given that the three previously-measured compact galaxies are outliers from $M_{\mathrm{BH}}-L_{\mathrm{bul}}$, while UGC 2698 is not, there may be significant scatter at the poorly sampled high-mass end of the relation. Additional gas-dynamical $M_{\mathrm{BH}}$ measurements for the compact galaxy sample will improve our understanding of BH$-$galaxy co-evolution.

Stellar evolution models predict the existence of a gap in the black hole mass spectrum from $\sim55 M_\odot - 120 M_\odot$ due to pair-instability supernovae (PISNe). We investigate the possible existence of such an "upper" mass gap in the second gravitational wave transient catalog (GWTC-2) by hierarchically modeling the astrophysical distribution of black hole masses. We extend the Truncated and Powerlaw+Peak mass distribution families to allow for an explicit gap in the mass distribution, and apply the extended models to GWTC-2. We find that with the Truncated model there is mild evidence favoring an upper mass gap with log Bayes Factor $\ln \mathcal{B} = 2.79$, inferring the lower and upper bounds at $56.12_{-4.38}^{+7.54} M_\odot$, and $103.74_{-6.32}^{+17.01} M_\odot$ respectively. When using the Powerlaw+Peak model, we find no preference for the gap. When imposing tighter priors on the gap bounds centered on the expected PISN gap bounds, the log Bayes factors in favor of a gap mildly increase. These results are however contingent on the parameter inference for the most massive binary, GW190521, for which follow up analyses showed the source may be an intermediate mass ratio merger that has component masses straddling the gap. Using the GW190521 posterior samples from the analysis in Nitz & Capano (2021), we find an increase in Bayes factors in favor of the gap. However, the overall conclusions are unchanged: There is no preference for a gap when using the Powerlaw+Peak model. This work paves the way for constraining the physics of pair-instability and pulsational pair-instability supernovae and high-mass black hole formation.

Cross correlation analyses of high resolution spectroscopic data have recently shown great success in directly detecting planetary signals and enabling the characterization of their atmospheres. One such technique aims to observe a system at multiple epochs and combine the measured planetary radial velocities from each epoch into a measurement of the planetary Keplerian orbital velocity $K_p$, constituting a direct detection of the planetary signal. Recent work has shown that in few-epoch ($\sim$5) data sets, unintended structure can arise at a high level, obscuring the planetary detection. In this work, we look to simulations to examine whether there are ways to reduce this structured noise in few-epoch data sets by careful planning of observations. The choice of observation date allows observers to select the primary (stellar) velocity - through a set systemic velocity and chosen barycentric velocity - and the planetary orbital phase, and so we focus on the effects of these two parameters. We find that epochs taken when the primary velocity is near zero, and the stellar lines remain relatively fixed to the telluric rest-frame, greatly reduce the level of structured noise and allow for much stronger planetary detections, on average more than twice the significance of detections made with epochs using randomly selected primary velocities. Following these results, we recommend that observers looking to build up high-resolution multi-epoch data sets target nights when their system has a near-zero primary velocity.

A spectral-energy distribution (SED) model for Type Ia supernovae (SNe Ia) is a critical tool for measuring precise and accurate distances across a large redshift range and constraining cosmological parameters. We present an improved model framework, SALT3, which has several advantages over current models including the leading SALT2 model (SALT2.4). While SALT3 has a similar philosophy, it differs from SALT2 by having improved estimation of uncertainties, better separation of color and light-curve stretch, and a publicly available training code. We present the application of our training method on a cross-calibrated compilation of 1083 SNe with 1207 spectra. Our compilation is $2.5\times$ larger than the SALT2 training sample and has greatly reduced calibration uncertainties. The resulting trained SALT3.K21 model has an extended wavelength range $2000$-$11000$ angstroms (1800 angstroms redder) and reduced uncertainties compared to SALT2, enabling accurate use of low-$z$ $I$ and $iz$ photometric bands. Including these previously discarded bands, SALT3.K21 reduces the Hubble scatter of the low-z Foundation and CfA3 samples by 15% and 10%, respectively. To check for potential systematic uncertainties we compare distances of low ($0.01<z<0.2$) and high ($0.4<z<0.6$) redshift SNe in the training compilation, finding an insignificant $2\pm14$ mmag shift between SALT2.4 and SALT3.K21. While the SALT3.K21 model was trained on optical data, our method can be used to build a model for rest-frame NIR samples from the Roman Space Telescope. Our open-source training code, public training data, model, and documentation are available at https://saltshaker.readthedocs.io/en/latest/, and the model is integrated into the sncosmo and SNANA software packages.

Ultralight axions and other bosons are dark matter candidates present in many high energy physics theories beyond the Standard Model. In particular, the string axiverse postulates the existence of up to $\mathcal{O}(100)$ light scalar bosons constituting the dark sector. Considering a mixture of axions and cold dark matter, we obtain upper bounds for the axion relic density $\Omega_a h^2 < 0.004$ for axions of mass $10^{-31}\;\mathrm{eV}\leq m_a \leq 10^{-26}\;\mathrm{eV}$ at 95% confidence. We also improve existing constraints by a factor of over 4.5 and 2.1 for axion masses of $10^{-25}$ eV and $10^{-32}$ eV, respectively. We use the Fourier-space galaxy clustering statistics from the Baryon Oscillation Spectroscopic Survey (BOSS) and demonstrate how galaxy surveys break important degeneracies in the axion parameter space compared to the cosmic microwave background (CMB). We test the validity of the effective field theory of large-scale structure approach to mixed ultralight axion dark matter by making our own mock galaxy catalogs and find an anisotropic ultralight axion signature in the galaxy quadrupole. We also observe an enhancement of the linear galaxy bias from 1.8 to 2.4 when allowing for 5% of the dark matter to be composed of a $10^{-28}$ eV axion in our simulations. Finally, we develop an augmented interpolation scheme allowing a fast computation of the axion contribution to the linear matter power spectrum leading to a 70% reduction of the computational cost for the full Monte Carlo Markov chains analysis.

We present a pixelized source reconstruction method applied on Integral Field Spectroscopic (IFS) observations of gravitationally lensed galaxies. We demonstrate the effectiveness of this method in a case study on the clumpy morphology of a $z \sim 2$ lensed galaxy behind a group-scale lens. We use a Bayesian forward source modelling approach to reconstruct the surface brightness distribution of the source galaxy on a uniformly pixelized grid while accounting for the image point spread function (PSF). The pixelated approach is sensitive to clump sizes down to 100 pc and resolves smaller clump sizes with an improvement in the signal to noise ratio (SNR) by almost a factor of ten compared with more traditional ray-tracing approaches.

Dust evolution in protoplanetary disks from small dust grains to pebbles is key to the planet formation process. The gas in protoplanetary disks should influence the vertical distribution of small dust grains ($\sim$1 $\mu m$) in the disk.Utilizing archival near-infrared polarized light and millimeter observations, we can measure the scale height and the flare parameter $\beta$ of the small dust grain scattering surface and $^{12}$CO gas emission surface for three protoplanetary disks IM Lup, HD 163296, and HD 97048 (CU Cha). For two systems, IM Lup and HD 163296, the $^{12}$CO gas and small dust grains at small radii from the star have similar heights but at larger radii ($>$100 au) the dust grain scattering surface height is lower than the $^{12}$CO gas emission surface height. In the case of HD 97048, the small dust grain scattering surface has similar heights to the $^{12}$CO gas emission surface at all radii. We ran a protoplanetary disk radiative transfer model of a generic protoplanetary disk with TORUS and showed that there is no difference between the observed scattering surface and $^{12}$CO emission surface. We also performed analytical modeling of the system and found that gas-to-dust ratios larger than 100 could explain the observed difference in IM Lup and HD 163296. This is the first direct comparison of observations of gas and small dust grain heights distribution in protoplanetary disks. Future observations of gas emission and near-infrared scattered light instruments are needed to look for similar trends in other protoplanetary disks.

In this work we consider strange stars formed by quark matter in the color-flavor-locked (CFL) phase of color superconductivity. The CFL phase is described by a Nambu-Jona-Lasinio model with four-fermion vector and diquark interaction channels. The effect of the color superconducting medium on the gluons are incorporated into the model by including the gluon self-energy in the thermodynamic potential. We construct parametrizations of the model by varying the vector coupling $G_V$ and comparing the results to the data on tidal deformability from the GW170817 event, the observational data on maximum masses from massive pulsars such as the MSP J0740+6620, and the mass/radius fits to NICER data for PSR J003+0451. Our results points out to windows for the $G_V$ parameter space of the model, with and without gluon effects included, that are compatible with all these astrophysical constraints, namely, $0.21<G_V/G_S<0.4$, and $0.02<G_V/G_S<0.1$, respectively. We also observe a strong correlation between the tidal deformabilites of the GW170817 event and $G_V$. Our results indicate that strange stars cannot be ruled out in collisions of compact binaries from the structural point of view.

In this paper, we study the galactic cosmic ray (GCR) variations over the solar cycles 23 and 24, with measurements from the NASA's ACE/CRIS instrument and the ground-based neutron monitors (NMs). The results show that the maximum GCR intensities of heavy nuclei (nuclear charge 5-28, 50-500 MeV/nuc) at 1 AU during the solar minimum in 2019-2020 break their previous records, exceeding those recorded in 1997 and 2009 by ~25% and ~6%, respectively, and are at the highest levels since the space age. However, the peak NM count rates are lower than those in late 2009. The difference between GCR intensities and NM count rates still remains to be explained. Furthermore, we find that the GCR modulation environment during the solar minimum P24/25 are significantly different from previous solar minima in several aspects, including remarkably low sunspot numbers, extremely low inclination of the heliospheric current sheet, rare coronal mass ejections, weak interplanetary magnetic field and turbulence. These changes are conducive to reduce the level of solar modulation, providing a plausible explanation for the record-breaking GCR intensities in interplanetary space.

We present an analytical model to identify thin discs in galaxies, and apply this model to a sample of SDSS MaNGA galaxies. This model fits the velocity and velocity dispersion fields of galaxies with regular kinematics. By introducing two parameters $\zeta$ related to the comparison of the model's asymmetric drift correction to the observed gas kinematics and $\eta$ related to the dominant component of a galaxy, we classify the galaxies in the sample as "disc-dominated", "non-disc-dominated", or "disc-free" indicating galaxies with a dominating thin disc, a non-dominating thin disc, or no thin disc detection with our method, respectively. The dynamical mass resulting from our model correlates with stellar mass, and we investigate discrepancies by including gas mass and variation of the initial mass function. As expected, most spiral galaxies in the sample are disc-dominated, while ellipticals are predominantly disc-free. Lenticular galaxies show a dichotomy in their kinematic classification, which is related to their different star formation rates and gas fractions. We propose two possible scenarios to explain these results. In the first scenario, disc-free lenticulars formed in more violent processes than disc-dominated ones, while in the second scenario, the quenching processes in lenticulars lead to a change in their kinematic structures as disc-dominated lenticulars evolve to disc-free ones.

We present a novel, iterative method using an empirical Bayesian approach for modeling the limb darkened WASP-121b transit from the TESS light curve. Our method is motivated by the need to improve $R_{p}/R_{\ast}$ estimates for exoplanet atmosphere modeling, and is particularly effective with the limb darkening (LD) quadratic law requiring no prior central value from stellar atmospheric models. With the non-linear LD law, the method has all the advantages of not needing atmospheric models but does not converge. The iterative method gives a different $R_{p}/R_{\ast}$ for WASP-121b at a significance level of 1$\sigma$ when compared with existing non-iterative methods. To assess the origins and implications of this difference, we generate and analyze light curves with known values of the limb darkening coefficients (LDCs). We find that non-iterative modeling with LDC priors from stellar atmospheric models results in an inconsistent $R_{p}/R_{\ast}$ at 1.5$\sigma$ level when the known LDC values are as those previously found when modeling real data by the iterative method. In contrast, the LDC values from the iterative modeling yields the correct value of $R_{p}/R_{\ast}$ to within 0.25$\sigma$. For more general cases with different known inputs, Monte Carlo simulations show that the iterative method obtains unbiased LDCs and correct $R_{p}/R_{\ast}$ to within a significance level of 0.3$\sigma$. Biased LDC priors can cause biased LDC posteriors and lead to bias in the $R_{p}/R_{\ast}$ of up to 0.82$\%$, 2.5$\sigma$ for the quadratic law and 0.32$\%$, 1.0$\sigma$ for the non-linear law. Our improvement in $R_{p}/R_{\ast}$ estimation is important when analyzing exoplanet atmospheres.

Optical depth variations in the Galactic neutral interstellar medium (ISM) with spatial scales from hundreds to thousands of astronomical units have been observed through HI absorption against pulsars and continuum sources, while extremely small structures with spatial scales of tens of astronomical units remain largely unexplored. The nature and formation of such tiny-scale atomic structures (TSAS) need to be better understood. We report a tentative detection of TSAS with a signal-to-noise ratio of 3.2 toward PSR B1557$-$50 in the second epoch of two Parkes sessions just 0.36 yr apart, which are the closest-spaced spectral observations toward this pulsar. One absorption component showing marginal variations has been identified. Based on the pulsar's proper motion of 14 mas $\rm yr^{-1}$ and the component's kinematic distance of 3.3 kpc, the corresponding characteristic spatial scale is 17 au, which is among the smallest sizes of known TSAS. Assuming a similar line-of-sight (LOS) depth, the tentative TSAS cloud detected here is overdense and overpressured relative to the cold neutral medium (CNM), and can radiatively cool fast enough to be in thermal equilibrium with the ambient environment. We find that turbulence is not sufficient to confine the overpressured TSAS. We explore the LOS elongation that would be required for the tentative TSAS to be at the canonical CNM pressure, and find that it is $\sim5000$ -- much larger than filaments observed in the ISM. We see some evidence of line width and temperature variations in the CNM components observed at the two epochs, as predicted by models of TSAS-like cloud formation colliding warm neutral medium flows.

KMT-2016-BLG-2605, with planet-host mass ratio $q=0.012\pm 0.001$, has the shortest Einstein timescale, $t_\e = 3.41\pm 0.13\,$days, of any planetary microlensing event to date. This prompts us to examine the full sample of 7 short ($t_\e<7\,$day) planetary events with good $q$ measurements. We find that six have clustered Einstein radii $\theta_\e = 115\pm 20\,\muas$ and lens-source relative proper motions $\mu_\rel\simeq 9.5\pm 2.5\,\masyr$. For the seventh, these two quantities could not be measured. These distributions are consistent with a Galactic-bulge population of very low-mass (VLM) hosts near the hydrogen-burning limit. This conjecture could be verified by imaging at first adaptive-optics light on next-generation (30m) telescopes. Based on a preliminary assessment of the sample, "planetary" companions (i.e., below the deuterium-burning limit) are divided into "genuine planets", formed in their disks by core accretion, and very low-mass brown dwarfs, which form like stars. We discuss techniques for expanding the sample, which include taking account of the peculiar "anomaly dominated" morphology of the KMT-2016-BLG-2605 light curve.

The discovery of a large putative impact crater buried beneath Hiawatha Glacier along the margin of the northwestern Greenland Ice Sheet has reinvigorated interest into the nature of large impacts into thick ice masses. This circular structure is relatively shallow and exhibits a small central uplift, whereas a peak-ring morphology is expected. This discrepancy may be due to long-term and ongoing subglacial erosion but may also be explained by a relatively recent impact through the Greenland Ice Sheet, which is expected to alter the final crater morphology. Here we model crater formation using hydrocode simulations, varying pre-impact ice thickness and impactor composition over crystalline target rock. We find that an ice-sheet thickness of 1.5 or 2 km results in a crater morphology that is consistent with the present morphology of this structure. Further, an ice sheet that thick substantially inhibits ejection of rocky material, which might explain the absence of rocky ejecta in most existing Greenland deep ice cores if the impact occurred during the late Pleistocene. From the present morphology of the putative Hiawatha impact crater alone, we cannot distinguish between an older crater formed by a pre-Pleistocene impact into ice-free bedrock or a younger, Pleistocene impact into locally thick ice, but based on our modeling we conclude that latter scenario is possible.

We have examined the light curves of the planetary nebulae (PNe) in all Kepler/K2 campaigns (0 through 19) to identify central star (CS) variability that may indicate a binary CS. We found no variable CS among the three PNe in Campaign 0, but we did identify one likely variable among the four PNe in Campaign 2, three variables among the 12 PN candidates in Campaign 7, one possible variable in the single PN in Campaign 15, and one very likely binary CS in the only PN in Campaign 16 (Abell 30). Our primary effort, though, was focused on Campaign 11 which targeted a Galactic bulge field having approximately 183 PNe in which we identified 21 strong candidate variable CS, plus another 9 possible, but less convincing, variables. The periods of our variables range from 2.3h to 30d. Most of the variables are likely to be binary stars. From our sample of 204 target PNe, we find the fraction of PNe having a binary (or, at least variable) CS to be 20-25 percent, depending on the selected subsample (e.g., all PN candidates or only `true' PNe) and details of the incompleteness correction. We believe that these fractions are lower limits for various reasons, primarily due to reduced sensitivity to detecting variability caused by dilution and noise from the nebula and neighbouring stars. This degradation is especially severe for longer period, highly separated systems.

We use the 21 cm emission line data from the DINGO-VLA project to study the atomic hydrogen gas H\,{\textsc i} of the Universe at redshifts $z<0.1$. Results are obtained using a stacking analysis, combining the H\,{\textsc i} signals from 3622 galaxies extracted from 267 VLA pointings in the G09 field of the Galaxy and Mass Assembly Survey (GAMA). Rather than using a traditional one-dimensional spectral stacking method, a three-dimensional cubelet stacking method is used to enable deconvolution and the accurate recovery of average galaxy fluxes from this high-resolution interferometric dataset. By probing down to galactic scales, this experiment also overcomes confusion corrections that have been necessary to include in previous single dish studies. After stacking and deconvolution, we obtain a $30\sigma$ H\,{\textsc i} mass measurement from the stacked spectrum, indicating an average H\,{\textsc i} mass of $M_{\rm H\,{\textsc i}}=(1.674\pm 0.183)\times 10^{9}~{\Msun}$. The corresponding cosmic density of neutral atomic hydrogen is $\Omega_{\rm H\,{\textsc i}}=(0.377\pm 0.042)\times 10^{-3}$ at redshift of $z=0.051$. These values are in good agreement with earlier results, implying there is no significant evolution of $\Omega_{\rm H\,{\textsc i}}$ at lower redshifts.

Simulating the long-term evolution of temperature and magnetic fields in neutron stars is a major effort in astrophysics, having significant impact in several topics. A detailed evolutionary model requires, at the same time, the numerical solution of the heat diffusion equation, the use of appropriate numerical methods to control non-linear terms in the induction equation, and the local calculation of realistic microphysics coefficients. Here we present the latest extension of the magneto-thermal 2D code in which we have coupled the crustal evolution to the core evolution, including ambipolar diffusion. It has also gained in modularity, accuracy, and efficiency. We revise the most suitable numerical methods to accurately simulate magnetar-like magnetic fields, reproducing the Hall-driven magnetic discontinuities. From the point of view of computational performance, most of the load falls on the calculation of microphysics coefficients. To a lesser extent, the thermal evolution part is also computationally expensive because it requires large matrix inversions due to the use of an implicit method. We show two representative case studies: (i) a non-trivial multipolar configuration confined to the crust, displaying long-lived small-scale structures and discontinuities; and (ii) a preliminary study of ambipolar diffusion in normal matter. The latter acts on timescales that are too long to have relevant effects on the timescales of interest but sets the stage for future works where superfluid and superconductivity need to be included.

Using the Large Binocular Telescope (LBT)/Multi-Object Dual Spectrograph (MODS), we have obtained optical spectroscopy of one of the most metal-poor dwarf star-forming galaxies (SFG) in the local Universe, J2229+2725. This galaxy with a redshift z=0.0762 was selected from the Data Release 16 (DR16) of the Sloan Digital Sky Survey (SDSS). Its properties derived from the LBT observations are most extreme among SFGs in several ways. Its oxygen abundance 12+logO/H=7.085+/-0.031 is among the lowest ever observed for a SFG. With its very low metallicity, an absolute magnitude Mg=-16.39 mag, a low stellar mass Mstar=9.1x10^6 Msun and a very low mass-to-light ratio Mstar/Lg~0.0166 (in solar units), J2229+2725 deviates strongly from the luminosity-metallicity relation defined by the bulk of the SFGs in the SDSS. J2229$+$2725 has a very high specific star-formation rate sSFR~75 Gyr^-1, indicating very active ongoing star formation. Three other features of J2229+2725 are most striking, being the most extreme among lowest-metallicity SFGs: 1) a ratio O32=I([OIII]5007)/I([OII]3727)~53, 2) an equivalent width of the Hbeta emission line EW(Hbeta) of 577A, and 3) an electron number density of ~1000 cm^-3. These properties imply that the starburst in J2229+2725 is very young. Using the extremely high O32 in J2229+2725, we have improved the strong-line calibration for the determination of oxygen abundances in the most metal-deficient galaxies, in the range 12 + logO/H<7.3.

The chemical compounds carrying the thiol group (-SH) have been considered essential in recent prebiotic studies regarding the polymerization of amino acids. We have searched for this kind of compounds toward the Galactic Centre quiescent cloud G+0.693-0.027. We report the first detection in the interstellar space of the trans-isomer of monothioformic acid (t-HC(O)SH) with an abundance of $\sim\,$1$\,\times\,$10$^{-10}$. Additionally, we provide a solid confirmation of the gauche isomer of ethyl mercaptan (g-C$_2$H$_5$SH) with an abundance of $\sim\,$3$\,\times\,$10$^{-10}$, and we also detect methyl mercaptan (CH$_3$SH) with an abundance of $\sim\,$5$\,\times\,$10$^{-9}$. Abundance ratios were calculated for the three SH-bearing species and their OH-analogues, revealing similar trends between alcohols and thiols with increasing complexity. Possible chemical routes for the interstellar synthesis of t-HC(O)SH, CH$_3$SH and C$_2$H$_5$SH are discussed, as well as the relevance of these compounds in the synthesis of prebiotic proteins in the primitive Earth.

We present a study of far and near-ultraviolet emission from the accretion disk in a powerful Seyfert 1 galaxy IC4329A using observations performed with the Ultraviolet Imaging Telescope (UVIT) onboard AstroSat. These data provide the highest spatial resolution and deepest images of IC4329A in the far and near UV bands acquired to date. The excellent spatial resolution of the UVIT data has allowed us to accurately separate the extended emission from the host galaxy and the AGN emission in the far and near UV bands. We derive the intrinsic AGN flux after correcting for the Galactic and internal reddening, as well as for the contribution of emission lines from the broad and narrow-line regions. The intrinsic UV continuum emission shows a marked deficit compared to that expected from the "standard" models of the accretion disk around an estimated black hole mass of 1-2x10^8Msun when the disk extends to the innermost stable circular orbit. We find that the intrinsic UV continuum is fully consistent with the standard disk models, but only if the disk emits from distances larger than 80-150 gravitational radii.

We study the morphology of the stellar periphery of the Magellanic Clouds in search of substructure using near-infrared imaging data from the VISTA Hemisphere Survey (VHS). Based on the selection of different stellar populations using the ($J-K_\mathrm{s}$, $K_\mathrm{s}$) colour-magnitude diagram, we confirm the presence of substructures related to the interaction history of the Clouds and find new substructures on the easter side of the LMC disc which may be owing to the influence of the Milky Way, and on the northern side of the SMC, which is probably associated to the ellipsoidal structure of the galaxy. We also study the luminosity function of red clump stars in the SMC and confirm the presence of a bi-modal distance distribution, in the form of a foreground population. We find that this bi-modality is still detectable in the eastern regions of the galaxy out to a 10 deg distance from its centre. Additionally, a background structure is detected in the North between 7 and 10 deg from the centre which might belong to the Counter Bridge, and a foreground structure is detected in the South between 6 and 8 deg from the centre which might be linked to the Old Bridge.

Self-gravitating astronomical objects often show a central plateau in the density profile (core) whose physical origin is hotly debated. Cores are theoretically expected in N-body systems of maximum entropy, however, they are not present in the canonical N-body numerical simulations of cold dark matter (CDM). Our work shows that despite this apparent contradiction between theory and numerical simulations, they are fully consistent. Simply put, cores are characteristic of systems in thermodynamic equilibrium, but thermalizing collisions are purposely suppressed in CDM simulations. When collisions are allowed, N-body numerical simulations develop cored density profiles, in perfect agreement with the theoretical expectation. We compare theory and two types of numerical simulations: (1) when DM particles are self-interacting (SIDM) with enough cross-section, then the effective two-body relaxation timescale becomes shorter than the Hubble time resulting in cored DM haloes. The haloes thus obtained, with masses from dwarf galaxies to galaxy clusters, collapse to a single shape after normalization, and this shape agrees with the polytropic density profile theoretically expected. (2) The inner radii in canonical N-body numerical simulations are always discarded because the use of finite-mass DM particles artificially increases the two-body collision rate. We show that the discarded radii develop cores that are larger than the employed numerical softening and have polytropic shapes independently of halo mass. Our work suggests that the presence of cores in simulated (or observed) density profiles can used as evidence for systems in thermodynamic equilibrium.

The absorption lines with high blue-shifted velocities are frequently found in the ultraviolet (UV) and X-ray spectra of luminous active galactic nuclei (AGNs). This implies that the high-velocity winds/outflows are common in AGNs. In order to study the formation of high-velocity winds, especially ultra-fast outflows (UFOs), we perform two-dimensional magnetohydrodynamic (MHD) simulations. Initially, a magnetic field is set to be weaker than the gas pressure at disk surface. In our simulations, line force operates on the region like filaments, because the X-ray radiation from corona is shielded by dense gas in the inner region at some angle. The location of filaments changes with time and then the line-driven winds are in exposure to X ray and become highly ionized. Line force at UV bands does not directly drive the highly-ionized winds. In the sense of time average, the properties of high-velocity winds meet the formation condition of UFOs. Compared with line force, the function of magnetic field is negligible in directly driving winds. In the MHD model, the region around the rotational axis becomes magnetic-pressure-dominated, which prevents gases from spreading to higher latitudes and then enhances the gas column density at middle and low latitudes (20$^{\rm o}$--70$^{\rm o}$). Higher column density is helpful to shield X-ray photons, which causes that line force is more effective in the MHD model than in the HD model. Higher-velocity winds with broader opening angle are produced in the MHD model.

We present evolutionary models for solar-like stars with an improved treatment of convection that results in a more accurate estimate of the radius and effective temperature. This is achieved by improving the calibration of the mixing-length parameter, which sets the length scale in the 1D convection model implemented in the stellar evolution code. Our calibration relies on the results of 2D and 3D radiation hydrodynamics simulations of convection to specify the value of the adiabatic specific entropy at the bottom of the convective envelope in stars as a function of their effective temperature, surface gravity and metallicity. For the first time, this calibration is fully integrated within the flow of a stellar evolution code, with the mixing-length parameter being continuously updated at run-time. This approach replaces the more common, but questionable, procedure of calibrating the length scale parameter on the Sun, and then applying the solar-calibrated value in modeling other stars, regardless of their mass, composition and evolutionary status. The internal consistency of our current implementation makes it suitable for application to evolved stars, in particular to red giants. We show that the entropy calibrated models yield a revised position of the red giant branch that is in better agreement with observational constraints than that of standard models.

We study the ionised ISM in NGC 7793 with MUSE/AO, at a spatial resolution of $\sim$ 10 pc. The data are complemented with young star clusters (YSCs), O stars and GMCs observed with HST and ALMA. Using a strong-line method, we find a median $12 + \log(O/H) \sim 8.37$ with a scatter of 0.25 dex, in agreement with previous estimates. The abundance map is rich in substructures, surrounding clusters and massive stars, although clear degeneracies with photoionisation are present. We determine the observed total amount of ionising photons, $Q(H^0)$, from the extinction corrected H$\alpha$ luminosity, and compare it to the expected $Q(H^0)$ obtained by summing the contributions of YSCs and massive stars, to obtain an escape fraction ($f_{esc}$). Overall, we find $f_{esc, HII} = 0.67_{-0.12}^{+0.08}$ for the population of HII regions. We also conclude that the sources of ionisation observed within the FoV are more than sufficient to explain the amount of diffuse ionised gas observed in this region of the galaxy. In general, we find that YSCs located in HII regions have a higher probability to be younger, less massive, and to emit a higher number of ionising photons than clusters in the rest of the field. Finally, we study the optical depth of the regions traced by [SII]/[OIII], finding no systematic trend between the resulting ionisation structure and $f_{esc}$. [abridged]

Magnetic fields influence the formation and evolution of stars and impact the observed stellar properties. Ap stars (magnetic A-type stars) are a prime example of this. Access to precise and accurate determinations of their stellar fundamental properties, such as masses and ages, is crucial to understand the origin and evolution of fossil magnetic fields. We propose using the radii and luminosities determined from interferometric measurements, in addition to seismic constraints when available, to infer fundamental properties of 14 Ap stars pr\'eviously characterised. We used a grid-based modelling approach, employing stellar models computed with the \textsc{cestam} stellar evolution code, and the parameter search performed with the \textsc{aims} optimisation method. The stellar model grid was built using a wide range of initial helium abundances and metallicities in order to avoid any bias originating from the initial chemical composition. The large frequency separations ($\Delta\nu$) of HR\,1217 (HD\,24712) and $\alpha$~Cir (HD\,128898), two rapidly oscillating Ap stars of the sample, were used as seismic constraints. We inferred the fundamental properties of the 14 stars in the sample. The overall results are consistent within $1\sigma$ with previous studies, however, the stellar masses inferred in this study are higher. This trend likely originates from the broader range of chemical compositions considered in this work. We show that the use of $\Delta\nu$ in the modelling significantly improves our inferences, allowing us to set reasonable constraints on the initial metallicity which is, otherwise, unconstrained. This gives an indication of the efficiency of atomic diffusion in the atmospheres of roAp stars and opens the possibility of characterising the transport of chemical elements in their interiors.

Despite being the best studied red supergiant star in our Galaxy, the physics behind the photometric variability and mass loss of Betelgeuse is poorly understood. Moreover, recently the star has experienced an unusual fading with its visual magnitude reaching a historical minimum. We investigate the nature of this event with the help of a recently developed tomographic technique. Tomography allows us to probe different depths in the stellar atmosphere and to recover the corresponding disk-averaged velocity field. We apply the tomographic method to a few-year time-series high-resolution spectroscopic observations of Betelgeuse in order to relate its atmospheric dynamics to the photometric variability. Our results show that a sudden increase of the molecular opacity in the cooler upper atmosphere of Betelgeuse is likely the reason of the observed unusual decrease of the star's brightness.

We provide the largest and most homogeneous sample of $\alpha$-element (Mg, Ca, Ti) and iron abundances for field RR Lyrae (RRLs, 162 variables) by using high-resolution spectra. The current measurements were complemented with similar abundances available in the literature for 46 field RRLs brought to our metallicity scale. We ended up with a sample of old (t$\ge$ 10 Gyr), low-mass stellar tracers (208 RRLs: 169 fundamental, 38 first overtone, 1 mixed mode) covering three dex in iron abundance (-3.00$\le$[Fe/H]$\le$0.24). We found that field RRLs are $\sim$0.3 dex more $\alpha$-poor than typical Halo tracers in the metal-rich regime, ([Fe/H]$\ge$-1.2) while in the metal-poor regime ([Fe/H]$\le$-2.2) they seem to be on average $\sim$0.1 dex more $\alpha$-enhanced. This is the first time that the depletion in $\alpha$-elements for solar iron abundances is detected on the basis of a large, homogeneous and coeval sample of old stellar tracers. Interestingly, we also detected a close similarity in the [$\alpha$/Fe] trend between $\alpha$-poor, metal-rich RRLs and red giants (RGs) in the Sagittarius dwarf galaxy as well as between $\alpha$-enhanced, metal-poor RRLs and RGs in ultra faint dwarf galaxies. These results are supported by similar elemental abundances for 46 field Horizontal Branch (HB) stars. These stars share with RRLs the same evolutionary phase and the same progenitors. This evidence further supports the key role that old stellar tracers play in constraining the early chemical enrichment of the Halo and, in particular, in investigating the impact that dwarf galaxies have had in the mass assembly of the Galaxy.

We discuss the spectral energy distributions and physical properties of six galaxies whose photometric redshifts suggest they lie beyond a redshift $z\simeq$9. Each was selected on account of a prominent excess seen in the Spitzer/IRAC 4.5$\mu$m band which, for a redshift above $z=9.0$, likely indicates the presence of a rest-frame Balmer break and a stellar component that formed earlier than a redshift $z\simeq10$. In addition to constraining the earlier star formation activity on the basis of fits using stellar population models with BAGPIPES, we have undertaken the necessary, but challenging, follow-up spectroscopy for each candidate using various combinations of Keck/MOSFIRE, VLT/X-shooter, Gemini/FLAMINGOS2 and ALMA. Based on either Lyman-$\alpha$ or [OIII] 88 $\mu$m emission, we determine a convincing redshift of $z$=8.78 for GN-z-10-3 and a likely redshift of $z$=9.28 for the lensed galaxy MACS0416-JD. For GN-z9-1, we conclude the case remains promising for a source beyond $z\simeq$9. Together with earlier spectroscopic data for MACS1149-JD1, our analysis of this enlarged sample provides further support for a cosmic star formation history extending beyond redshifts $z\simeq$10. We use our best-fit stellar population models to reconstruct the past rest-frame UV luminosities of our sources and discuss the implications for tracing earlier progenitors of such systems with the James Webb Space Telescope.

In the currently favored cosmological paradigm galaxies form hierarchically through the accretion of numerous satellite galaxies. Since the satellites are much less massive than the host halo, they occupy a small fraction of the volume in action space defined by the potential of the host halo. Since actions are conserved when the potential of the host halo changes adiabatically, stars from an accreted satellite are expected to remain clustered in action space as the host halo evolves. In this paper, we identify accreted satellites in three Milky Way like disk galaxies from the cosmological baryonic FIRE-2 simulations by tracking satellite galaxies through simulation snapshots. We then try to recover these satellites by applying the cluster analysis algorithm Enlink to the orbital actions of accreted star particles in the present-day snapshot. We define several metrics to quantify the success of the clustering algorithm and use these metrics to identify well-recovered and poorly-recovered satellites. We plot these satellites in the infall time-progenitor mass (or stellar mass) space, and determine the boundaries between the well-recovered and poorly-recovered satellites in these two spaces with classification tree method. The groups found by Enlink are more likely to correspond to a real satellite if they have high significance, a quantity assigned by Enlink. Since cosmological simulations predict that most stellar halos have a population of insitu stars, we test the ability of Enlink to recover satellites when the sample is contaminated by 10-50% of insitu star particles, and show that most of the satellites well-recovered by Enlink in the absence of insitu stars, stay well-recovered even with 50% contamination. We thus expect that, in the future, cluster analysis in action space will be useful in upcoming data sets (e.g. Gaia) for identifying accreted satellites in the Milky Way.

The non-thermal acceleration of particles in magnetohydrodynamic (MHD) turbulence plays a central role in a wide variety of astrophysical sites. This physics is addressed here in the context of a strong turbulence, composed of coherent structures rather than waves, outside the realm of quasilinear theory. The present description tracks the momentum of the particle along its history through a sequence of (non-inertial) frames in which the electric field vanishes, in the spirit of the original Fermi scenario. It thus connects the sources of energy gain (or loss) to the gradients of the velocity of the magnetic field lines, in particular the acceleration and the shear of their velocity flow projected along the field line direction, as well as their compression in the transverse plane. Momentum diffusion coefficients associated to these different contributions are expressed as a function of the power spectrum of those velocity fluctuations. Importantly, those velocity gradients are subject to intermittency: they are spatially localized and their strengths obey powerlaw distributions, as demonstrated through direct measurements in the incompressible MHD simulation of the Johns Hopkins University database. This intermittency impacts the acceleration process in a significant way, which opens up prospects for a rich phenomenology, in particular non-thermal spectra with extended powerlaw tails. Those findings generally agree with recent kinetic numerical simulations. Extensions to this description and possible avenues of exploration are discussed.

The large disparity in physical conditions from the diffuse interstellar medium (ISM) to denser clouds such as photon-dominated regions (PDRs) triggers an evolution of the dust properties (i.e. composition, size, and shape). The gas physics and chemistry are tightly connected to these dust properties and are therefore affected by dust evolution and especially the nano-grain depletion in the outer irradiated part of PDRs. We highlight the influence of nano-grain depletion on the gas physics and chemistry in the Horsehead nebula, a prototypical PDR. We used a model for atomic and molecular gas in PDRs, the Meudon PDR code, using diffuse ISM-like dust and Horsehead-like dust to study the influence of nano-grain depletion on the gas physics and chemistry, focusing on the impact on photoelectric heating and H2 formation and, therefore, on the H2 gas lines. We find that nano-grain depletion in the Horsehead strongly affects gas heating through the photoelectric effect and thus the gas temperature and the H2 formation, hence the H -> H2 position. Consequently, the first four pure rotational lines of H2 (e.g. 0-0 S(0), S(1), S(2), and S(3)) vary by a factor of 2 to 14. The 0-0 S(3) line that is often underestimated in models is underestimated even more when taking nano-grain depletion into account due to the decrease in gas heating through the photoelectric effect. This strongly suggests that our understanding of the excitation of H2 and/or of heating processes in the Horsehead, and more generally in PDRs, is still incomplete. Nano-grain depletion in the outer part of the Horsehead has a strong influence on several gas tracers that will be prominent in JWST observations of irradiated clouds. We therefore need to take this depletion into account in order to improve our understanding of the Horsehead, and more generally PDRs, and to contribute to the optimal scientific return of the mission.

In the modelling of blazars and other radiating plasma flows, correctly evolving in time a series of kinetic equation coupled for different particle species is a computationally challenging problem. In this paper we introduce the code Katu, capable of simultaneously evolving a large number of particle species including photons, leptons and hadrons and their interactions. For each particle, the numerical expressions for the kinetic equations are optimized for efficiency and the code is written to easily allow for extensions and modifications. Additionally, we tie Katu to two pieces of statistical software, emcee and pymultinest, which makes a suite ready to apply to blazars and extract relevant statistical information from their electromagnetic (and neutrino, if applicable) flux. As an example, we apply Katu to Mrk 421 and compare a leptonic and a leptohadronic model. In this first ever direct Bayesian comparison between full kinetic simulations of these models, we find substantial evidence (Bayes factor $\sim7$) for the leptohadronic model purely from the electromagnetic spectrum. The code is available online at https://github.com/hveerten/katu under the LGPL license.

We use both new and archival ALMA data of three energy lines each of CN and HCN to explore intensity ratios in dense gas in NGC 3256, NGC 7469, and IRAS 13120-5453. The HCN (3-2)/HCN (1-0) intensity ratio varies in NGC 3256 and NGC 7469, with superlinear trends of 1.53$\pm$0.07 and 1.55$\pm$0.05, respectively. We find an offset to higher HCN (3-2)/HCN (1-0) intensity ratios (~0.8) in IRAS 13120-5453 compared to NGC 3256 (~0.3-0.4) and NGC 7469 (~0.3-0.5). The HCN (4-3)/HCN (3-2) intensity ratio in NGC 7469 has a slope of 1.34$\pm$0.05. We attribute the variation within NGC 3256 to excitation associated with the northern and southern nuclei. In NGC 7469, the variations are localized to the region surrounding the active galactic nucleus. At our resolution (~700 pc), IRAS 13120-5453 shows little variation in the HCN intensity ratios. Individual galaxies show nearly constant CN (2-1)/CN (1-0) intensity ratios. We find an offset to lower CN (2-1)/CN (1-0) intensity ratios (~0.5) in NGC 3256 compared to the other two galaxies (~0.8). For the CN (3-2)/CN (2-1) intensity ratio, NGC 7469 has a superlinear trend of 1.55$\pm$0.04, with the peak localized toward the active galactic nucleus. We find high (~1.7) CN (1-0)/HCN (1-0) intensity ratios in IRAS 13120-5453 and in the northern nucleus of NGC 3256, compared to a more constant ratio (~1.1) in NGC 7469 and non-starbursting regions of NGC 3256.

The observing strategy of a galaxy survey influences the degree to which its resulting data can be used to accomplish any science goal. LSST is thus seeking metrics of observing strategies for multiple science cases in order to optimally choose a cadence. Photometric redshifts are essential for many extragalactic science applications of LSST's data, including but not limited to cosmology, but there are few metrics available, and they are not straightforwardly integrated with metrics of other cadence-dependent quantities that may influence any given use case. We propose a metric for observing strategy optimization based on the potentially recoverable mutual information about redshift from a photometric sample under the constraints of a realistic observing strategy. We demonstrate a tractable estimation of a variational lower bound of this mutual information implemented in a public code using conditional normalizing flows. By comparing the recoverable redshift information across observing strategies, we can distinguish between those that preclude robust redshift constraints and those whose data will preserve more redshift information, to be generically utilized in a downstream analysis. We recommend the use of this versatile metric to observing strategy optimization for redshift-dependent extragalactic use cases, including but not limited to cosmology, as well as any other science applications for which photometry may be modeled from true parameter values beyond redshift.

Silver (Ag) mirrors for astronomical telescopes consist of multiple metallic and dielectric thin films. Furthermore, the topmost surface of such Ag mirrors needs to be covered by a protection coating. While the protection coating is often deposited at room temperature and the entire mirrors are also handled at room temperature, various thin film deposition techniques offer protection coatings with improved characteristics when carried out at elevated temperatures. Thus, in this work, high-performance Ag mirrors were designed and fabricated with a new benchmark. The resulting Ag mirrors were annealed (i.e., post-fabrication annealing) at various temperatures to investigate the viability of introducing thermal processes during and/or after fabrication in improving overall optical performance and durability of protected silver mirrors. In our experiments, Ag mirror samples were deposited by electron-beam evaporation and subsequently annealed at various temperatures in the range from 60 {\deg}C to 300 {\deg}C, and then the mirror samples underwent an environmental stress test at 80 {\deg}C and 80% humidity for 10 days. While all the mirror samples annealed below 200 {\deg}C showed negligible corrosion after undergoing the stress testing, those annealed below 160 {\deg}C presented spectral reflectivity comparable to or higher than that of as-deposited reference samples. In contrast, the mirror samples annealed above 200 {\deg}C exhibited significant degradation after the stress testing. The comprehensive analysis indicated that delamination and voids caused by the growth of Ag grains during the annealing are the primary mechanisms of the degradation.

We analyze high-resolution (400pc) 220GHz continuum and CO(2-1) ALMA observations of a representative sample of 22 local (z<0.165) ULIRG systems (32 individual nuclei) as part of the "Physics of ULIRGs with MUSE and ALMA" (PUMA) project. The deconvolved half-light radii of the 220GHz continuum sources are between <60-350 pc (median 90pc). We associate these regions with the regions emitting the bulk of the infrared luminosity. The good agreement, within a factor of 2, between the 220GHz fluxes and the extrapolation of the infrared gray-body, and the small synchrotron and free-free contributions support this assumption. The cold molecular gas emission sizes are 60-700 pc and are on average 2.6 times larger than the continuum. We derive L_IR and cold molecular gas surface densities: log Sigma(L_IR)=11.5-14.3 Lsun/kpc^2 and log Sigma(H2)=2.9-4.2 Msun/pc^2. Assuming that the L_IR is produced by star-formation, this corresponds to median Sigma(SFR)=2500 Msun/yr/kpc^2 which would imply extremely short depletion times, <1-15 Myr, and unphysical SF efficiencies >1 for 70% of the sample. Therefore, this favors the presence of obscured AGN that could dominate the L_IR. We also classify the ULIRG nuclei in two groups: (a) compact nuclei (r<130 pc) with high mid-IR excess emission found in optically classified AGN; and (b) nuclei following a relation with decreasing mid-IR excess for decreasing r. The majority, 65%, of the interacting nuclei lie in the low end (<130 pc) of this relation, while only 25% of the advanced mergers do so. This suggests that in the early stages of the interaction, the activity occurs in more compact and dust-obscured regions than in advanced merger stages. Approximately two thirds of the nuclei are above the Eddington limit which is consistent with the detection of massive outflows in local ULIRGs and the potential role of radiation pressure in the launching process.

We investigate the 3D spin alignment of galaxies with respect to the large-scale filaments using the MaNGA survey. The cosmic web is reconstructed from the Sloan Digital Sky Survey using Disperse and the 3D spins of MaNGA galaxies are estimated using the thin disk approximation with integral field spectroscopy kinematics. Late-type spiral galaxies are found to have their spins parallel to the closest filament's axis. The alignment signal is found to be dominated by low-mass spirals. Spins of S0-type galaxies tend to be oriented preferentially in perpendicular direction with respect to the filament's axis. This orthogonal orientation is found to be dominated by S0s that show a notable misalignment between their kinematic components of stellar and ionised gas velocity fields and/or by low mass S0s with lower rotation support compared to their high mass counterparts. Qualitatively similar results are obtained when splitting galaxies based on the degree of ordered stellar rotation, such that galaxies with high spin magnitude have their spin aligned, and those with low spin magnitude in perpendicular direction to the filaments. In the context of conditional tidal torque theory, these findings suggest that galaxies' spins retain memory of their larger-scale environment. In agreement with measurements from hydrodynamical cosmological simulations, the measured signal at low redshift is weak, yet statistically significant. The dependence of the spin-filament orientation of galaxies on their stellar mass, morphology and kinematics highlights the importance of sample selection to detect the signal.

Starting from a Born-Oppenheimer decomposition of the Wheeler-DeWitt equation for the quantum cosmology of the matter-gravity system, we have performed a Wigner-Weyl transformation and obtained equations involving a Wigner function for the scale factor and its conjugate momentum. This has allowed us to study in more detail than previously the approach to the classical limit of gravitation and the way time emerges in such a limit. To lowest order we reproduce the Friedmann equation and the previously obtained equation for the evolution of matter. We also obtain expressions for higher order corrections to the semi-classical limit.

Gravitational-wave data can be used to test general relativity in the extreme, highly nonlinear, and strong-field regime. Modified gravity theories such as Einstein-dilation-Gauss-Bonnet and dynamical Chern-Simons gravity are supposed to be well checked with some specific signals. We analyze gravitational-wave data from the first half of the third observing run of Advanced LIGO/Virgo to place constraints on the parameters of these two theories. The dynamical Chern-Simons gravity remains unconstrained. While for the Einstein-dilation-Gauss-Bonnet gravity, we have $\sqrt{\alpha_{\rm EdGB}}\lesssim 0.40\,\rm km$ with the data of GW190814 and $\sqrt{\alpha_{\rm EdGB}}\lesssim 0.24\,\rm km$ if GW190425 is a neutron star$-$black hole merger event, as allowed by the current observations. Such constraints are improved by a factor of $\sim 10-20$ in comparison to that set by the GWTC-1 events. We also demonstrate that the Bayes approach is more effective than the reweight method to constrain the Einstein-dilation-Gauss-Bonnet model.

We develop an analytical approach, verified by the nice agreement with the numerical relativity simulation results, to calculate the quasinormal mode spectrum of the Kerr-Newman black hole. Then we analyze the gravitational wave data with the ringdown waveform model including both the fundamental mode and the overtone modes, and find that it can efficiently constrain the charge of the source. For GW150914, the charge-to-mass ratio of the remnant black hole is constrained to be $\leq 0.38$ at $90\%$ credibility. Our waveform model can be widely applied to other GW150914 like events. With the sole ringdown data, it is capable of constraining the electric, magnetic, or other $U(1)$ dark charges carried by black holes, as well as the deviation parameter ($\alpha$) of the scalar-tensor-vector gravity. Indeed, a constraint of $\alpha\leq 0.17$ is achieved with the ringdown data alone for the first time.

The most promising indirect search for the existence of axion dark matter uses radio telescopes to look for narrow spectral lines generated from the resonant conversion of axions in the magnetospheres of neutron stars. Unfortunately, a large list of theoretical uncertainties has prevented this search strategy from being fully accepted as robust. In this work we attempt to address major outstanding questions related to the role and impact of the plasma, including: $(i)$ does refraction and reflection of radio photons in the magnetosphere induce strong inhomogeneities in the flux, $(ii)$ can refraction induce premature axion-photon de-phasing, $(iii)$ to what extent do photon-plasma interactions induce a broadening of the spectral line, $(iv)$ does the flux have a strong time dependence, and $(v)$ can radio photons sourced by axions be absorbed by the plasma. We present an end-to-end analysis pipeline based on ray-tracing that exploits a state-of-the-art auto-differentiation algorithm to propagate photons from the conversion surface to asymptotically large distances. Adopting a charge symmetric Goldreich-Julian model for the magnetosphere, we show that for reasonable parameters one should expect a strong anisotropy of the signal, refraction induced axion-photon de-phasing, significant line-broadening, a variable time-dependence of the flux, and, for large enough magnetic fields, anisotropic absorption. Our simulation code is flexible enough to serve as the basis for follow-up studies with a large range of magnetosphere models.

An exact and analytical solution of four dimensional vacuum General Relativity representing a system of two static black holes at equilibrium is presented. The metric is completely regular outside the event horizons, both from curvature and conical singularities. The balance between the two Schwarzschild sources is granted by an external gravitational field, without the need of extra matter fields besides gravity, nor strings or struts. The geometry of the solution is analysed. The Smarr law, the first and the second law of black hole thermodynamics are discussed.

In this work, we apply the tools of the dynamical system theory in order to revisit and uncover the structure of a nongravitational interaction between pressureless dark matter and dark energy described by a scalar field, which has been previously investigated in the literature. For a coupling function $Q = -(\alpha \dot{\rho}_m + \beta \dot{\rho}_{\phi} )$, we have found that it can be rewritten in the form $Q = 3H (\alpha \rho_m + \beta \dot{\phi}^2)/(1-\alpha +\beta)$, so that its dependence on the dark matter density and on the kinetic term of the scalar field is linear and proportional to the Hubble parameter. We analyze the following scenarios $\alpha=0$, $\alpha = \beta$ and $\alpha = -\beta$, separately and in order to describe the cosmological evolution for each solution we have calculated various observables. We find that there are not any new stable late-time solutions apart from those found of standard quintessence, nevertheless, the stability conditions are severely altered. A notable result found with respect to previous works is that in our case, with the exception of the matter dominated solution, the remaining critical points behave as scaling although the stiff matter solution and the dark energy dominated state can be recovered in the limit $\beta \rightarrow 0$ and $\beta \rightarrow 1$, respectively. Moreover, it is shown that for $\alpha = \beta $ and $\alpha = - \beta$ (in general for $\alpha \neq 0$), a separatrix arises modifying prominently the structure of the phase space. This represents a novel feature no mentioned before in the literature.

We calculate the exact solutions to the equations of motion that govern the light ray trajectories as they travel in a Kerr black hole's exterior that is considered to be filled with an inhomogeneous and anisotropic plasmic medium. This is approached by characterizing the plasma through conceiving a radial and an angular structure function, which are let to be constant. The description of the motion is carried out by using the Hamilton-Jacobi method, that allows defining two effective potentials, characterizing the evolution of the polar coordinates. The (hyper-)elliptic integrals of motion are then solved analytically, and the evolution of coordinates is expressed in terms of the Mino time. This way, the three-dimensional demonstrations of the light ray trajectories are given respectively.

A new formalism for the nonlinear Alfv\'enic states sustainable in Hall Magnetohydrodynamics is developed in a complete basis provided by the circularly polarized Beltrami Vectors, the eigenstates of linear HMHD. Nonlinear HMHD is, then, reduced to a rather simple looking set of scalar equations from which a model problem of three interacting Beltrami modes is formulated and analytically solved. The triplet interactions span a variety of familiar nonlinear processes leading to a redistribution as well as periodic exchange of energy. The energy exchange processes (whose strength is measured by an energy exchange /depletion time) will, perhaps, play a dominant role in determining the spectral content of an eventual Alfv\'enic state. All nonlinearities (sensitive functions of the interacting wave vectors) operate at par, and none is dominant over any substantial region of k-space; their intricate interplay prevents a "universal" picture from emerging; few generalizations on the processes that may, for instance, lead to a turbulent state, are possible. However, the theory can definitely claim: 1) the energy tends to flow from lower to higher $k$, and 2) the higher $k_z$ (in the direction of the ambient magnetic field) components of a mode with a given $k$ are depleted/oscillate faster -in some cases much faster. It is noteworthy that the mode coupling is the strongest (with the shortest depletion time) when the participating wave vectors are nearly perpendicular; perhaps, an expected consequence of the curl (cross product) nonlinearities. Numerical simulations will be necessary to help create a fully reliable picture.

Bayesian parameter estimation of gravitational waves from compact binary coalescence (CBC) typically requires more than millions of evaluations of computationally expensive template waveforms. We propose a technique to reduce the cost of waveform generation by exploiting the chirping behavior of CBC signal. Our technique does not require waveforms at all frequencies in the frequency range used in the analysis, and does not suffer from the fixed cost due to the upsampling of waveforms. Our technique speeds up the parameter estimation of typical binary neutron star signal by a factor of $\mathcal{O}(10)$ for the low-frequency cutoff of $20\,\mathrm{Hz}$, and $\mathcal{O}(10^2)$ for $5\,\mathrm{Hz}$. It does not require any offline preparations or accurate estimates of source parameters provided by detection pipelines.

We argue that in quantum gravity there is no Born rule. The quantum-gravity regime, described by a non-normalisable Wheeler-DeWitt wave functional $\Psi$, must be in quantum nonequilibrium with a probability distribution $P\neq\left\vert \Psi\right\vert ^{2}$ (initially and always). A Born rule can emerge only in the semiclassical regime of quantum systems on a classical spacetime background, with normalisable Schr\"{o}dinger wave functions $\psi$. Conditioning on the underlying quantum-gravitational ensemble yields a nonequilibrium distribution $\rho\neq\left\vert \psi\right\vert ^{2}$ at the beginning of the semiclassical regime, with quantum relaxation $\rho\rightarrow\left\vert \psi\right\vert ^{2}$ taking place only afterwards. Quantum gravity naturally creates an early nonequilibrium universe. We also show how small corrections to the Schr\"{o}dinger equation yield an intermediate regime in which the Born rule is unstable: an initial distribution $\rho=\left\vert \psi\right\vert ^{2}$ can evolve to a final distribution $\rho\neq\left\vert \psi\right\vert ^{2}$. These results arise naturally in the de Broglie-Bohm pilot-wave formulation of quantum gravity. We show that quantum instability during inflation generates a large-scale deficit $\sim1/k^{3}$ in the primordial power spectrum at wavenumber $k$, though the effect is too small to observe. Similarly we find an unobservably large timescale for quantum instability in a radiation-dominated universe. Quantum instability may be important in black-hole evaporation, with a final burst of Hawking radiation that violates the Born rule. Deviations from the Born rule can also be generated for atomic systems in the gravitational field of the earth, though the effects are unlikely to be observable. The most promising scenario for the detection of Born-rule violations appears to be in radiation from exploding primordial black holes.

In review of the White Papers from the Voyage 2050 process and after the public presentation of a number of these papers in October 2019 in Madrid, we as White Paper lead authors have identified a coherent science theme that transcends the divisions around which the Topical Teams are structured. This note aims to highlight this synergistic science theme and to make the Topical Teams and the Voyage 2050 Senior Committee aware of the wide importance of these topics and the broad support that they have across the worldwide science community.

In this work we present the foundations of generalized scalar-tensor theories arising from vector bundle constructions, and we study the kinematic, dynamical and cosmological consequences. In particular, over a pseudo-Riemannian space-time base manifold, we define a fiber structure with two scalar fields. The resulting space is a 6-dimensional vector bundle endowed with a non-linear connection. We provide the form of the geodesics and the Raychaudhuri and general field equations, both in Palatini and metrical method. When applied at a cosmological framework, this novel geometrical structure induces extra terms in the modified Friedmann equations, leading to the appearance of an effective dark energy sector, as well as of an interaction of the dark mater sector with the metric. We show that we can obtain the standard thermal history of the universe, with the sequence of matter and dark-energy epochs, and furthermore the effective dark-energy equation-of-state parameter can lie in the quintessence or phantom regimes, or exhibit the phantom-divide crossing.

Magnetic field amplification by relativistic streaming plasma instabilities is central to a wide variety of high-energy astrophysical environments as well as to laboratory scenarios associated with intense lasers and electron beams. We report on a new secondary nonlinear instability which arises for relativistic dilute electron beams after the saturation of the linear Weibel instability. This instability grows due to the transverse magnetic pressure associated with the beam current filaments, which cannot be quickly neutralized due to the inertia of background ions. We show that it can amplify the magnetic field strength and spatial scale by orders of magnitude, leading to large-scale plasma cavities with strong magnetic field and to very efficient conversion of the beam kinetic energy into magnetic energy. The instability growth rate, saturation level, and scale length are derived analytically and shown to be in good agreement with fully-kinetic simulations.