Papers for Thursday, Jan 27 2022

We present a search for extreme emission line galaxies (EELGs) at z<1 in the COSMOS and North Ecliptic Pole (NEP) fields with imaging from Subaru/Hyper Suprime-Cam (HSC) and a combination of new and existing spectroscopy. We select EELGs on the basis of substantial excess flux in the z broadband, which is sensitive to H$\alpha$ at 0.3<z<0.42 and [OIII]$\lambda$5007 at 0.7<z<0.86. We identify 10,470 galaxies with z excesses in the COSMOS dataset and 91,385 in the NEP field. We cross-reference the COSMOS EELG sample with the zCOSMOS and DEIMOS 10k spectral catalogs, finding 1395 spectroscopic matches. We made an additional 71 (46 unique) spectroscopic measurements with Y<23 using the HYDRA multi-object spectrograph on the WIYN 3.5m telescope, and 204 spectroscopic measurements from the DEIMOS spectrograph on the Keck II telescope, providing a total of 1441/10,470 spectroscopic redshifts for the EELG sample in COSMOS (~14%). We confirm that 1418 (~98%) are H$\alpha$ or [OIII]$\lambda$5007 emitters in the above stated redshift ranges. We also identify 240 redshifted H$\alpha$ and [OIII]$\lambda$5007 emitters in the NEP using spectra taken with WIYN/HYDRA and Keck/DEIMOS. Using broadband selection techniques in g-r-i color space, we distinguish between H$\alpha$ and [OIII]$\lambda$5007 emitters with 98.6% accuracy. We test our EELG selection by constructing H$\alpha$ and [OIII]$\lambda$5007 luminosity functions and comparing to recent literature results. We conclude that broadband magnitudes from HSC, the Vera C. Rubin Observatory, and other deep optical multi-band surveys can be used to select EELGs in a straightforward manner.

We introduce a new model for understanding AGN continuum variability. We start from a Shakura--Sunyaev thin accretion disc with a steady-state radial temperature profile $T(R)$ and assume that the variable flux is due to axisymmetric temperature perturbations $\delta T(R,t)$. After linearizing the equations, we fit UV-optical AGN lightcurves to determine $\delta T(R,t)$ for a sample of seven AGNs. We see a diversity of $|\delta T/T| \sim 0.1$ fluctuation patterns which are not dominated by outgoing waves traveling at the speed of light as expected for the "lamppost" model used to interpret disc reverberation mapping studies. Rather, the most common pattern resembles slow ($v \ll c$) ingoing waves. An explanation for our findings is that these ingoing waves trigger central temperature fluctuations that act as a lamppost, producing lower amplitude temperature fluctuations moving outwards at the speed of light. The lightcurves are dominated by the lamppost signal -- even though the temperature fluctuations are dominated by other structures with similar variability time-scales -- because the discs exponentially smooth the contributions from the slower ($v \ll c$) moving fluctuations to the observed lightcurves. This leads to lightcurves that closely resemble the expectations for a lamppost model but with the slow variability time-scales of the ingoing waves. This also implies that longer time-scale variability signals will increasingly diverge from lamppost models because the smoothing of slower-moving waves steadily decreases as their period or spatial wavelength increases.

Subsurface convection zones are ubiquitous in early-type stars. Driven by narrow opacity peaks, these thin convective regions transport little heat but play an important role in setting the magnetic properties and surface variability of stars. Here we demonstrate that these convection zones are \emph{not} present in as wide a range of stars as previously believed. In particular, there are regions which 1D stellar evolution models report to be convectively unstable but which fall below the critical Rayleigh number for onset of convection. For sub-solar metallicity this opens up a \emph{stability window} in which there are no subsurface convection zones. For LMC metallicity this surface stability region extends roughly between $8M_\odot$ and $16M_\odot$, increasing to $8M_\odot$ -- $35M_\odot$ for SMC metallicity. Such windows are then an excellent target for probing the relative influence of subsurface convection and other sources of photometric variability in massive stars.

(Abridged) In this paper, we use the Tuffs et al. attenuation - inclination models in ultraviolet (UV), optical, and near-infrared (NIR) bands to investigate the average global dust properties in galaxies as a function of stellar mass$M_{*}$, stellar mass surface density $\mu_{*}$, star-formation rate $SFR$, specific star-formation rate $sSFR$, star-formation main-sequence offset $dMS$, and star-formation rate surface density $\Sigma_{SFR}$ at redshifts $z \sim 0$ and $z \sim 0.7$. We use star-forming galaxies from SDSS ($\sim$ 20000) and GAMA ($\sim$ 2000) to form our low-z sample at $0.04 < z < 0.1$ and star-forming galaxies from COSMOS ($\sim$ 2000) for the sample at $0.6 <z < 0.8$. We find that galaxies at $z \sim 0.7$ have higher optical depth $\tau_{B}^{f}$ and clumpiness $F$ than galaxies at $z \sim 0$. The increase in $F$ hints that the stars of $z \sim 0.7$ galaxies are less likely to escape their birth cloud, which might indicate that the birth clouds are larger. We also found that $\tau_{B}^{f}$ increases with $M_{*}$ and $\mu_{*}$independent of sample and therefore redshift. We found no clear trends in $\tau_{B}^{f}$ or $F$ with $SFR$, which could imply that the dust mass distribution is independent of $SFR$. In turn, this would imply that the balance of dust formation and destruction is independent of the $SFR$. Based on an analysis of the inclination-dependence of the Balmer decrement, we find that reproducing the Balmer line emission requires not only a completely optically thick dust component associated with star forming regions, as in the standard Tuffs et al. model, but an extra component of optically thin dust within the birth clouds. This new component implies the existence of dust inside HII regions that attenuates the Balmer emission before it escapes through gaps in the birth cloud and we find it is more important in high-mass galaxies.

Despite traditional thinking, an appreciable population of massive black holes may be lurking in dwarf galaxies. Prior to the last decade, nearly all massive black holes were found in the nuclei of giant galaxies and the existence of massive black holes in dwarf galaxies was highly controversial. The field has now been transformed with a growing community of researchers working on a variety of observational studies and theoretical models of dwarf galaxies hosting massive black holes. Work in this area is not only important for a holistic understanding of dwarf galaxy evolution and feedback, but it may just tell us how the first "seeds" of massive black holes formed in the early Universe. In this Perspective, I discuss the current state of the field as well as future prospects. I also present new insights on the demographics of nearby dwarf galaxies, which can be used to help constrain the black hole occupation/active fraction as a function of mass and dwarf galaxy type.

The production and distribution of metals in the diffuse intergalactic medium (IGM) have implications for galaxy formation models and the baryon (re)cycling process. The relative abundance of metals in high versus low-ionization states has also been argued to be sensitive to the Universe's reionization history. However, measurements of the background metallicity of the IGM at z~4 are sparse and in poor agreement with one another, and reduced sensitivity in the near-IR renders detecting individual metal absorbers nearly impossible. We present a new clustering-based technique that enables the detection of these weak IGM absorbers by statistically averaging over all spectral pixels, here applied to the C IV forest. We simulate the z=4.5 IGM with different models of inhomogeneous metal distributions and investigate its two-point correlation function (2PCF) using mock skewers of the C IV forest. The 2PCF demonstrates a clear peak at the doublet separation of the C IV line. The peak amplitude scales quadratically with metallicity, while enrichment morphology affects both the shape and amplitude of the 2PCF. The effect of enrichment topology can also be framed in terms of the metal mass- and volume-filling factors, and we show their trends as a function of the enrichment topology. For models consistent with the distribution of metals at z~3, we find that we can constrain [C/H] to within 0.2 dex, log$\,M_{\rm{min}}$ to within 0.4 dex, and $R$ to within 15%. We show that strong absorbers arising from the circumgalactic medium of galaxies can be easily identified and masked, allowing one to recover the underlying IGM signal. The auto-correlation of the metal-line forest presents a new and compelling avenue to simultaneously constrain IGM metallicity and enrichment topology with high precision at z>4, thereby pushing such measurements into the Epoch of Reionization.

We study the relationships between stellar mass, size and age within the quiescent population, using two mass-complete spectroscopic samples with $\mathrm{log_{10}}(M_{\star}/\mathrm{M_{\odot}})>10.3$, taken from VANDELS at $1.0<z<1.3$, and LEGA-C at $0.6<z<0.8$. Using robust D$_{n}$4000 values, we demonstrate that the well-known 'downsizing' signature is already in place by $z\simeq1.1$, with D$_{n}$4000 increasing by $\simeq0.1$ across a $\simeq$ 1 dex mass interval for both VANDELS and LEGA-C. We then proceed to investigate the evolution of the quiescent galaxy stellar mass-size relation from $z\simeq1.1$ to $z\simeq0.7$. We find the median size increases by a factor of $1.9\pm{0.1}$ at $\mathrm{log_{10}}(M_{\star}/\mathrm{M_{\odot}})=10.5$, and see tentative evidence for flattening of the relation, finding slopes of $\alpha=0.72\pm0.06$ and $\alpha=$ $0.56\pm0.04$ for VANDELS and LEGA-C respectively. We finally split our sample into galaxies above and below our fitted mass-size relations, to investigate how size and D$_{n}$4000 correlate. For LEGA-C, we see a clear difference, with larger galaxies found to have smaller D$_{n}$4000 at fixed stellar mass. Due to the faintness and smaller numbers of the VANDELS sample, we cannot confirm whether a similar relation exists at $z\simeq1.1$. We consider whether differences in stellar age or metallicity are most likely to drive this size-D$_{n}$4000 relation, finding that any metallicity differences are unlikely to fully explain the observed offset, meaning smaller galaxies must be older than their larger counterparts. We find the observed evolution in size, mass and D$_{n}$4000 across the $\simeq2$ Gyr from $z\sim1.1$ to $z\sim0.7$ can be explained by a simple toy model in which VANDELS galaxies evolve passively, whilst experiencing a series of minor mergers.

Possible M-type contact binaries were investigated by selecting W UMa-type variables with orbital periods below the 0.22-day cutoff. Gaia parallaxes were combined with Gaia and 2MASS photometry to obtain G and J-band absolute magnitudes of 674 red variable stars catalogued as contact binaries in the VSX database. The absolute magnitude of main sequence cool dwarfs varies strongly with spectral type, which was used to create a selection of 218 systems with primaries potentially of spectral types M0-M3. Lightcurves of the 46 systems with lowest near-infrared luminosities were inspected individually to confirm or reject the W UMa classification where possible. The extinction limits and amplitudes were combined with the calculated absolute magnitudes in order to set upper limits for the masses of the primary components. This is achieved via the consideration that the luminosity of the primary is no less than that of a main sequence star of the same mass. For 26 possible contact binaries the mass of the primary was limited to less than 0.5 solar masses with some systems having upper limits for the primary as low as 0.35 solar masses. These rare systems are intriguing targets for more detailed observations since they are on the low-mass end of the contact binary distribution and their formation and evolution are still unclear.

Measured spectral shifts due to intrinsic stellar variability (e.g., pulsations, granulation) and activity (e.g., spots, plages) are the largest source of error for extreme precision radial velocity (EPRV) exoplanet detection. Several methods are designed to disentangle stellar signals from true center-of-mass shifts due to planets. The EXPRES Stellar Signals Project (ESSP) presents a self-consistent comparison of 22 different methods tested on the same extreme-precision spectroscopic data from EXPRES. Methods derived new activity indicators, constructed models for mapping an indicator to the needed RV correction, or separated out shape- and shift-driven RV components. Since no ground truth is known when using real data, relative method performance is assessed using the total and nightly scatter of returned RVs and agreement between the results of different methods. Nearly all submitted methods return a lower RV RMS than classic linear decorrelation, but no method is yet consistently reducing the RV RMS to sub-meter-per-second levels. There is a concerning lack of agreement between the RVs returned by different methods. These results suggest that continued progress in this field necessitates increased interpretability of methods, high-cadence data to capture stellar signals at all timescales, and continued tests like the ESSP using consistent data sets with more advanced metrics for method performance. Future comparisons should make use of various well-characterized data sets -- such as solar data or data with known injected planetary and/or stellar signals -- to better understand method performance and whether planetary signals are preserved.

We describe a new polarized imaging pipeline implemented in the FHD software package. The pipeline is based on the optimal mapmaking imaging approach and performs horizon-to-horizon image reconstruction in all polarization modes. We discuss the formalism behind the pipeline's polarized analysis, describing equivalent representations of the polarized beam response, or Jones matrix. We show that, for arrays where antennas have uniform polarization alignments, defining a non-orthogonal instrumental polarization basis enables accurate and efficient image reconstruction. Finally, we present a new calibration approach that leverages widefield effects to perform fully-polarized calibration. This analysis pipeline underlies the analysis of Murchison Widefield Array (MWA) data in Byrne et al. 2022.

Recently, many different pulsar timing array (PTA) collaborations have reported strong evidence for a common stochastic process in their data sets. The reported amplitudes are in tension with previously computed upper limits. In this paper, we investigate how using a subset of a set of pulsars biases Bayesian upper limit recovery. We generate 500 simulated PTA data sets based on the NANOGrav 11-year data set with an injected stochastic gravitational wave background (GWB). We then compute upper limits by sampling individual pulsar likelihoods, and combine them through a factorized version of the PTA likelihood to obtain upper limits on the GWB amplitude using different numbers of pulsars. We find that it is possible to recover an upper limit below the injected value, and that it is significantly more likely for this to occur when using a subset of pulsars to compute the upper limit. When picking pulsars to induce the maximum possible bias, we find that the upper limit recovered is below the injected value in 10.6% (53 of 500) of realizations. Further, we find that if we choose a subset of pulsars in order to obtain a lower upper limit than when using the full set of pulsars, the distribution of upper limits obtained from these 500 realizations is shifted to lower amplitude values.

We discuss a hypothetical existential threat from a 10 km diameter comet discovered 6 months prior to impact. We show that an extension of our work on bolide fragmentation using an array of penetrators, but modified with small nuclear explosive devices (NED) in the penetrators, combined with soon-to-be-realized heavy lift launch assets with positive $C_3$ such as NASA SLS or SpaceX Starship (with in-orbit refueling) is sufficient to mitigate this existential threat. A threat of this magnitude hitting the Earth at a closing speed of 40 km/s would have an impact energy of roughly 300 Teratons TNT, or about 40 thousand times larger than the current combined nuclear arsenal of the entire world. This is similar in energy to the KT extinction event that killed the dinosaurs some 66 million years ago. Such an event, if not mitigated, would be an existential threat to humanity. We show that mitigation is conceivable using existing technology, even with the short time scale of 6 months warning, but that the efficient coupling of the NED energy is critical.

Discoveries of gaps in data have been important in astrophysics. For example, there are kinematic gaps opened by resonances in dynamical systems, or exoplanets of a certain radius that are empirically rare. A gap in a data set is a kind of anomaly, but in an unusual sense: Instead of being a single outlier data point, situated far from other data points, it is a region of the space, or a set of points, that is anomalous compared to its surroundings. Gaps are both interesting and hard to find and characterize, especially when they have non-trivial shapes. We present methods to address this problem. First, we present a methodological approach to identify critical points, a criterion to select the most relevant ones and use those to trace the `valleys' in the density field. We then build on the observed properties of critical points to propose a novel gappiness criterion that can be computed at any point in the data space. This allows us to identify a broader variety of gaps, either by highlighting regions of the data-space that are `gappy' or by selecting data points that lie in local under densities. We also explore methodological ways to make the detected gaps robust to changes in the density estimation and noise in the data. We illustrate our methods on the velocity distribution of nearby stars in the Milky Way disk plane, which exhibits gaps that could originate from different processes. Identifying and characterizing those gaps could help determine their origins.

Using data from a year-long dedicated campaign to observe bright stars, we study the crosstalk channels present in the GPC1 camera. By analyzing these data, we construct a dataset that checks source stars on almost every CCD of every chip within the camera against all possible crosstalk destinations. We use a clustering algorithm to find potential crosstalk occurrences, and then also check all possible combinations (driven by the hardware layout) by eye. This results in a total of 640 rules, with a flux attenuation factor ranging from 2.5x10$^{2}$ for the bright end to 2.5$\times$10$^{4}$ at the faint end. The average value of m$_{cross}$-m$_{src}\approx$-10.25 corresponds to an attenuating factor of 1.25x10$^{4}$, which produces crosstalk ghosts with an average signal-to-noise ratio of 0.64$\pm$0.1 on the bright images. We find no evidence of crosstalk signals between CCDs not connected in the hardware setup. The distribution of attenuation factors is also found to be dependent on crosstalk movement. A clear dependence on cell column offsets is found, consistent with the idea that the source star charge is progressively attenuated during the traversal of cell readout lines. While we can see the trends, the uncertainties on aperture magnitude measurements are large at this stage.

We present CluSTAR-ND, a fast hierarchical galaxy/(sub)halo finder that produces Clustering Structure via Transformative Aggregation and Rejection in N-Dimensions. It is designed to improve upon Halo-OPTICS -- an algorithm that automatically detects and extracts significant astrophysical clusters from the 3D spatial positions of simulation particles -- by decreasing run-times, possessing the capability for metric adaptivity, and being readily applicable to data with any number of features. We directly compare these algorithms and find that not only does CluSTAR-ND produce a similarly robust clustering structure, it does so in a run-time that is at least $3$ orders of magnitude faster. In optimising CluSTAR-ND's clustering performance, we have also effectively made CluSTAR-ND a parameterless clustering algorithm since its parameters will be chosen automatically and optimally based on the input data. We conclude that CluSTAR-ND is a robust astrophysical clustering algorithm that can be leveraged to find stellar satellite groups on large synthetic or observational data sets.

The motion of a satellite can experience secular resonances between the precession frequencies of its orbit and the mean motion of the host planet around the star. Some of these resonances can significantly modify the eccentricity (evection resonance) and the inclination (eviction resonance) of the satellite. In this paper, we study in detail the secular resonances that can disturb the orbit of a satellite, in particular the eviction-like ones. Although the inclination is always disturbed while crossing one eviction-like resonance, capture can only occur when the semi-major axis is decreasing. This is, for instance, the case of Phobos, the largest satellite of Mars, that will cross some of these resonances in the future because its orbit is shrinking owing to tidal effects. We estimate the impact of resonance crossing in the orbit of the satellite, including the capture probabilities, as a function of several parameters, such as the eccentricity and the inclination of the satellite, and the obliquity of the planet. Finally, we use the method of the frequency map analysis to study the resonant dynamics based on stability maps, and we show that some of the secular resonances may overlap, which leads to chaotic motion for the inclination of the satellite.

We present spectroscopic and high-precision photometric observations, spanning the optical UV to the far red, before, during, and after the NASA Deep Impact event of July 4, 2005. The inner 2000 km of the pre and post-impact coma was about 0.3 magnitude redder in B-R than in the outer coma. The pre-impact spectrum was a faint reflected solar spectrum dominated by molecular emissions extending > 40000 km from the nucleus. The post-impact light curve in R and I showed a rapid rise consistent with an expanding optically thick cloud during the first 18 minutes after impact. During the next 8 minutes the cloud became optically thin. Sixty minutes after impact the impact R-band flux reached a plateau at 7.5 x 10-15 erg cm-2 s-1 {\AA}-1, the comet brightening by a factor of ~4.3 above its pre-impact value observed in a 15" aperture. The mean expansion velocity of the grains during the first 49 minutes was 229 {\pm} 49 m s-1. The spectrum became dominated by scattered sunlight during the first hour after impact. The volume scattering function (VSF) observed 32 minutes after impact shows strong reddening. At 49 minutes, however, the VSF shows an additional two-fold increase in the blue but only a 20 per-cent increase at 5500{\AA}. Post-impact spectra and R-I photometry showed rapid reddening. The particle size distribution, dominated by 1 to 2.5 micron particles shortly after impact, changed dramatically during the first hour due to sublimation of water-ice particles of this size. On the night following impact the comet was still substantially brighter than before impact, but R-I had returned to its pre-impact value. B-R remained significantly redder. The ejecta 25 hours after impact was fan-shaped subtending ~180{\deg} roughly symmetrical about position angle 225{\deg}. The mean expansion velocity 90{\deg} from the direction to the Sun was 185{\deg} 12 m s-1.

A limiting systematic effect in 21-cm interferometric experiments is the chromaticity due to the coupling between the sky and the instrument. This coupling is sourced by the instrument primary beam; therefore it is important to know the beam to extremely high precision. Here we demonstrate how known beam uncertainties can be characterized using databases of beam models. In this introductory work, we focus on beam errors arising from physically offset and/or broken antennas within a station. We use the public code OSKAR to generate an "ideal" SKA beam formed from 256 antennas regularly-spaced in a 35-m circle, as well as a large database of "perturbed" beams sampling distributions of broken/offset antennas. We decompose the beam errors ("ideal" minus "perturbed") using Principal Component Analysis (PCA) and Kernel PCA (KPCA). Using 20 components, we find that PCA/KPCA can reduce the residual of the beam in our datasets by 60-90% compared with the assumption of an ideal beam. Using a simulated observation of the cosmic signal plus foregrounds, we find that assuming the ideal beam can result in 1% error in the EoR window and 10% in the wedge of the 2D power spectrum. When PCA/KPCA is used to characterize the beam uncertainties, the error in the power spectrum shrinks to below 0.01% in the EoR window and <1% in the wedge. Our framework can be used to characterize and then marginalize over uncertainties in the beam for robust next-generation 21-cm parameter estimation.

We discuss a scenario in which TeV neutrinos are produced during explosions of Novae. It is argued that hadrons are accelerated to very high energies in the inner part of a Nova wind, as a result of reconnection of the strong magnetic field of a White Dwarf. Hadrons are expected to interact efficiently with a dense matter of the wind, either already during the acceleration process or during their advection with the equatorial wind. We calculate the neutrino spectra, and estimate the muon neutrino event rates in the IceCube telescope, in the case of a few Novae. In general, those event rates are unlikely to be detected with the present neutrino detectors. However, for favourable location of the observer, some neutrino events might be detected not only from the class of Novae recently detected in the GeV $\gamma$-rays by the {\it Fermi}-LAT telescope but also from novae not detected in $\gamma$-rays. The GeV $\gamma$-ray emission observed from Novae cannot originate in terms of the model discussed here, since protons are accelerated within a few stellar radii of the White dwarf, i.e. in the region in which GeV $\gamma$-rays are expected to be severely absorbed in the interactions with the radiation field and the matter of the wind.

Testing the Weak Equivalence Principle (WEP) to a precision of $10^{-15}$ requires a quantity of data that give enough confidence on the final result: ideally, the longer the measurement the better the rejection of thestatistical noise. The science sessions had a duration of 120 orbits maximum and were regularly repeated and spaced out to accommodate operational constraints but also in order to repeat the experiment in different conditions and to allow time to calibrate the instrument. Several science sessions were performed over the 2.5 year duration of the experiment. This paper aims to describe how the data have been produced on the basis of a mission scenario and a data flow process, driven by a tradeoff between the science objectives and the operational constraints. The mission was led by the Centre National d'Etudes Spatiales (CNES) which provided the satellite, the launch and the ground operations. The ground segment was distributed between CNES and Office National d'Etudes et de Recherches A\'erospatiales (ONERA). CNES provided the raw data through the Centre d'Expertise de Compensation de Tra\^{i}n\'{e}e (CECT: Drag-free expertise centre). The science was led by the Observatoire de la C\^ote d{'}Azur (OCA) and ONERA was in charge of the data process. The latter also provided the instrument and the Science Mission Centre of MICROSCOPE (CMSM).

The formation of the interstellar complex organic molecules (iCOMs) is a hot topic in astrochemistry. One of the main paradigms trying to reproduce the observations postulates that iCOMs are formed on the ice mantles covering the interstellar dust grains as a result of radical--radical coupling reactions. We investigate iCOMs formation on the icy surfaces by means of computational quantum mechanical methods. In particular, we study the coupling and direct hydrogen abstraction reactions involving the CH$_3$ + X systems (X = NH$_2$, CH$_3$, HCO, CH$_3$O, CH$_2$OH) and HCO + Y (Y = HCO, CH$_3$O, CH$_2$OH), plus the CH$_2$OH + CH$_2$OH and CH$_3$O + CH$_3$O systems. We computed the activation energy barriers of these reactions as well as the binding energies of all the studied radicals, by means of density functional theory (DFT) calculations on two ice water models, made of 33 and 18 water molecules. Then, we estimated the efficiency of each reaction using the reaction activation, desorption and diffusion energies and derived kinetics with the Eyring equations. We find that radical--radical chemistry on surfaces is not as straightforward as usually assumed. In some cases, direct H abstraction reactions can compete with radical--radical couplings, while in others they may contain large activation energies. Specifically, we found that (i) ethane, methylamine and ethylene glycol are the only possible products of the relevant radical--radical reactions; (ii) glyoxal, methyl formate, glycolaldehyde, formamide, dimethyl ether and ethanol formation is likely in competition with the respective H-abstraction products, and (iii) acetaldehyde and dimethyl peroxide do not seem a likely grain surface products.

This study explores the gravitational lensing effects of supermassive black holes (SMBHs) in galaxy clusters. While the presence of central SMBHs in galaxies is firmly established, recent work from high-resolution simulations predict the existence of an additional population of wandering SMBHs. Though the masses of these SMBHs are a minor perturbation on the larger scale and individual galaxy scale dark matter components in the cluster, they can impact statistical lensing properties and individual lensed image configurations. Probing for these potentially observable signatures, we find that SMBHs imprint detectable signatures in rare, higher-order strong lensing image configurations although they do not manifest any statistically significant detectable evidence in either the magnification distribution or the integrated shear profile. Investigating specific lensed image geometries, we report that a massive, near point-like, potential of an SMBH causes the following detectable effects: (ii) image splitting leading to the generation of extra images; (ii) positional and magnification asymmetries in multiply imaged systems; and (iii) the apparent disappearance of a lensed counter-image. Of these, image splitting inside the cluster tangential critical curve, is the most prevalent notable observational signature. We demonstrate these possibilities in two cases of observed giant arcs in $SGAS\,J003341.5+024217$ and $RX\,J1347.5-1145$, wherein specific image configurations seen can be reproduced with SMBHs. Future observations with high-resolution instrumentation (e.g. MAVIS-Very Large Telescope, MICADO-Extremely Large Telescope, and the upgraded ngVLA, along with data from the \textit{Euclid} \& \textit{Nancy Grace Roman} Space Telescopes and the Rubin LSST Observatory are likely to allow us to probe these unique yet rare SMBHs lensing signatures.

We present a study of the current state of knowledge concerning spacecraft operations and potential hazards while operating near a comet nucleus. Starting from simple calculations comparing the cometary coma environment to benign conditions on Earth, we progress to sophisticated engineering models of spacecraft behavior, and then confront these models with recent spacecraft proximity operations experience. Finally, we make recommendations from lessons learned for future spacecraft missions that enter into orbit around a comet for long-term operations. All of these considerations indicate that, with a proper spacecraft design and operations planning, the near-nucleus environment can be a relatively safe region in which to operate, even for an active short period comet near perihelion with gas production rates as high as 1e29 molecules/s. With gas densities similar to those found in good laboratory vacuums, dust densities similar to Class 100 cleanrooms, dust particle velocities of 10s of m/s, and microgravity forces that permit slow and deliberate operations, the conditions around a comet are generally more benign than a typical day on Mars. Even in strong dust jets near the nucleus surface, dust densities tend to be only a few grains/cm3, about the same as in a typical interior room on Earth. Stochastic forces on a modern spacecraft with tens of square meters of projected surface area can be accounted for using modern Attitude Control Systems to within tens of meters navigation error; surface contamination issues are only important for spacecraft spending months to years within a few kilometers of the nucleus surface; and the issues the Rosetta spacecraft faced, confusion of celestial star trackers by sunlit dust particles flying past the spacecraft, will be addressed using the next generation of star trackers implementing improved transient rejection algorithms.

The Giant Stellar Stream (GSS) is a prominent tidal feature in the halo of the Andromeda Galaxy (M31), representing the ongoing destruction of a satellite galaxy. In this paper, we investigate the formation of the GSS through detailed numerical simulations of the tidal disruption of a progenitor system. Assuming that the stream was created in a single merger event between M31 and a dwarf spheroidal galaxy with stellar mass of $10^{9}M_{\odot}$, we successfully reproduce the dynamical properties of the GSS. As the metallicity distribution along the stream has been well determined from the observations (PAndAS and AMIGA data sets), we use Monte Carlo simulations to reconstruct the original metallicity distribution of the dwarf progenitor. We find that a progenitor dwarf galaxy with a negative radial metallicity gradient, $\Delta$FeH = -$0.3 \pm 0.2$, reproduces the observed GSS properties at a time between 2.4 and 2.9 Gyrs into the merger. We also show that the observed double peak metallicity distribution along the stream is a transitory structure caused by unique merger circumstances where two groups of streaming stars are moving in opposite directions, intersecting to produce the peaks.

Many of the Milky Way's globular clusters are likely accreted from satellite galaxies that have long since merged with the Milky Way. When these globular clusters are susceptible to tidal disruption, this process likely starts already inside the parent satellite leading to an early stellar stream within the satellite. When the parent satellite merges with the Milky Way, the globular cluster and its pre-merger stellar stream are accreted in a somewhat chaotic process. Here, we investigate the properties of the accreted stream after the merger as we would see it today using a suite of simulations of accretion events. We find that the accretion process leads to a wide range of behaviors, but generally scatters the accreted stream over a wide, two-dimensional area of the sky. The behavior ranges from a set of a few or more well-defined "sub-streams" extending out from the post-merger thin stream by tens of degrees, to more widely dispersed debris over much of the sky, depending on how close to the center of the MilkyWay the merger happened. Using mock Gaia-like observations of the simulated streams, we demonstrate that an accreted-stream component can explain the off-track features observed in the GD-1 stream. Sub-streams can appear like thin tidal streams themselves that are seemingly unassociated with the post-merger stream, raising the possibility that some of the progenitor-less streams observed in the Milky Way are part of a single or a few accreted streams created in an ancient merger event.

NGC 1566 is a changing look AGN known to exhibit recurrent X-ray outbursts with each lasting for several years. The most recent X-ray outburst is observed on 2018, with a substantial increase of 2--10 keV flux by a factor of ~24 than the historical minimum. We re-analyze the XMM-Newton and NuSTAR observations covering the pre-outburst, outburst and post-outburst epochs, and confirm the discovery of the broad feature in the ~5--7 keV band during the period of outburst that could be interpreted as a relativistic Fe K_alpha emission line. Our analysis suggests that its flux has increased in tandem with the 2--10 keV continuum, making it the second changing look AGN in which the broad Fe K_alpha line responds to the X-ray continuum variability. This behavior strongly supports the idea that X-rays originates in a corona above the accretion disk, and disk reflection produces the relativistic Fe K_alpha line. In addition, we find the response of narrow Fe K_alpha emission line to the changes in the X-ray continuum on a time-scale as short as four months, allowing to put the location of line-emitting region at <0.1 pc, comparable to the size of optical BLR. By comparing to the changing look AGN NGC 2992, the Fe K_alpha variation rate (the ratio of Fe K_alpha variation to luminosity variation) in NGC 1566 appears greater, which could be possibly explained by larger amount of gas or Fe abundance responsible for producing the Fe K_alpha line for the latter. The strength of variable broad Fe K_alpha line as well as the soft X-ray excess emission appears to be correlated with the accretion rate, which could be explained as due to the state transition associated with the changing-look phenomenon.

We present the results of the reverberation monitoring aimed at MgII broad line and FeII pseudocontinuum for the luminous quasar CTS C30.10 (z = 0.90052) with the Southern African Large Telescope covering the years 2012-2021. We aimed at disentangling the MgII and UV FeII variability and the first measurement of UV FeII time delay for a distant quasar. We used several methods for time-delay measurements and determined both FeII and MgII time delays as well as performed a wavelength-resolved time delay study for a combination of MgII and FeII in the 2700 - 2900 \AA restframe wavelength range. We obtain the time delay for MgII of $275.5^{+12.4}_{-19.5}$ days in the rest frame, while for FeII we have two possible solutions of $270.0^{+13.8}_{-25.3}$ days and $180.3^{+26.6}_{-30.0}$ in the rest frame. Combining this result with the old measurement of FeII UV time delay for NGC 5548 we discuss for the first time the radius-luminosity relation for UV FeII with the slope consistent with $0.5$ within uncertainties. Since FeII time delay has a shorter time-delay component but lines are narrower than MgII, we propose that the line delay measurement is biased towards the BLR part facing the observer, with the bulk of the Fe II emission may arise from the more distant BLR region, one that is shielded from the observer.

We evaluate the cosmological coalescence and detection rates for massive black hole (MBH) binaries targeted by the gravitational wave observatory Laser Interferometer Space Antenna ({\it LISA}). Our calculation starts with a population of gravitationally unbound MBH pairs, drawn from the TNG50-3 cosmological simulation, and follows their orbital evolution from kpc scales all the way to coalescence using a semi-analytic model developed in our previous work. We find that for a majority of MBH pairs that coalesce within a Hubble time dynamical friction is the most important mechanism that determines their coalescence rate. Our model predicts a MBH coalescence rate $\lesssim 0.45$~yr$^{-1}$ and a {\it LISA} detection rate $\lesssim 0.34$~yr$^{-1}$. Most {\it LISA} detections should originate from $10^{\rm 6} - 10^{\rm 6.8}\,M_{\rm \odot}$ MBHs in gas-rich galaxies at redshifts $1.6 \leq z \leq 2.4$, and have a characteristic signal to noise ratio SNR $\sim 100$. We however find a dramatic reduction in the coalescence and detection rates, as well as the average SNR, if the effects of radiative feedback from accreting MBHs are taken into account. In this case, the MBH coalescence rate is reduced by $78\%$ (to $\lesssim 0.1$~yr$^{-1}$), and the \textit{LISA} detection rate is reduced by $94\%$ (to $0.02$~yr$^{-1}$), whereas the average SNR is $\sim 10$. We emphasize that our model provides a lower limit on the \textit{LISA} detection rate, consistent with other works in the literature that draw their MBH pairs from cosmological simulations.

How does a smooth cosmic distance ladder emerge from observations made from a single location in a lumpy Universe? Distances to the Type Ia supernova (SN1A) in the Hubble flow are anchored on local distance measurements to sources that are not in the Hubble flow. We described how this configuration could be built in a perturbed universe where lumpiness is described as small perturbations on top of a flat Friedmann-Lema{\i}tre Robertson-Walker (FLRW) spacetime. We show that there is a non-negligible modification (about 11\%) to the background FLRW area distance due to the presence of inhomogeneities in the immediate neighbourhood of an observer. We find that the modification is sourced by the electric part of the Weyl tensor indicating a tidal deformation of the local spacetime of the observer. We show in detail how it impacts the calibration of the SN1A absolute magnitude in the Hubble flow. We show that it resolves the SN1A absolute magnitude and Hubble tension simultaneously without the need for any of the dark stuff.

The structure of the Small Magellanic Cloud (SMC) is very complex, in particular in the periphery that suffers more from the interactions with the Large Magellanic Cloud (LMC). A wealth of observational evidence has been accumulated revealing tidal tails and bridges made up of gas, stars and star clusters. Nevertheless, a full picture of the SMC outskirts is only recently starting to emerge with a 6D phase-space map plus age and metallicity using star clusters as tracers. In this work, we continue our analysis of another outer region of the SMC, the so-called West Halo, and combined it with the previously analysed Northern Bridge. We use both structures to define the Bridge and Counter-bridge trailing and leading tidal tails. These two structures are moving away from each other, roughly in the SMC-LMC direction. The West Halo form a ring around the SMC inner regions that goes up to the background of the Northern Bridge shaping an extended layer of the Counter-bridge. Four old Bridge clusters were identified at distances larger than 8 kpc from the SMC centre moving towards the LMC, which is consistent with the SMC-LMC closest distance of 7.5 kpc when the Magellanic Bridge was formed about 150Myr ago; this shows that the Magellanic Bridge was not formed only by pulled gas, but it also removed older stars from the SMC during its formation. We also found age and metallicity radial gradients using projected distances on sky, which are vanished when we use the real 3D distances.

The $\ion{He}{I}$ 10833\,$\AA$ triplet is a powerful tool for characterising the upper atmosphere of exoplanets and tracing possible mass loss. Here, we analysed one transit of GJ\,1214\,b observed with the CARMENES high-resolution spectrograph to study its atmosphere via transmission spectroscopy around the $\ion{He}{I}$ triplet. Although previous studies using lower resolution instruments have reported non-detections of $\ion{He}{I}$ in the atmosphere of GJ\,1214\,b, we report here the first potential detection. We reconcile the conflicting results arguing that previous transit observations did not present good opportunities for the detection of $\ion{He}{I}$, due to telluric H$_2$O absorption and OH emission contamination.We simulated those earlier observations, and show evidence that the planetary signal was contaminated. From our single non-telluric-contaminated transit, we determined an excess absorption of 2.10$^{+0.45}_{-0.50}$\,\% (4.6\,$\sigma$) with a full width at half maximum (FWHM) of 1.30$^{+0.30}_{-0.25}$\,\AA. The detection of \ion{He}{I} is statistically significant at the 4.6\,$\sigma$ level, but repeatability of the detection could not be confirmed due to the availability of only one transit. By applying a hydrodynamical model and assuming an H/He composition of 98/2, we found that GJ\,1214\,b would undergo hydrodynamic escape in the photon-limited regime, losing its primary atmosphere with a mass-loss rate of (1.5--18)\,$\times$\,10$^{10}$\,g\,s$^{-1}$ and an outflow temperature in the range of 2900--4400\,K. Our $\ion{He}{I}$ excess absorption is the first tentative detection of a chemical species in the atmosphere of this benchmark sub-Neptune planet.

Recent observations have revealed remarkable insights into the gas reservoir in the circumgalactic medium (CGM) of galaxy haloes. In this paper, we characterise the gas in the vicinity of Milky Way and Andromeda analogues in the HESTIA (High resolution Environmental Simulations of The Immediate Area) suite of constrained Local Group (LG) simulations. The HESTIA suite comprise of a set of three high-resolution {\sc arepo}-based simulations of the LG, run using the Auriga galaxy formation model. For this paper, we focus only on the $z = 0$ simulation datasets and generate mock skymaps along with a power spectrum analysis to show that the distributions of ions tracing low-temperature gas (HI and SiIII) are more clumpy in comparison to warmer gas tracers (OVI, OVII and OVIII). We compare to the spectroscopic CGM observations of M31 and low-redshift galaxies. HESTIA under-produces the column densities of the M31 observations, but the simulations are consistent with the observations of low-redshift galaxies. A possible explanation for these findings is that the spectroscopic observations of M31 are contaminated by gas residing in the CGM of the Milky Way.

Dyson spheres are hypothetical megastructures built by advanced extraterrestrial civilizations to harvest radiation energy from stars. Here, we combine optical data from Gaia DR2 with mid-infrared data from AllWISE to set the strongest upper limits to date on the prevalence of partial Dyson spheres within the Milky Way, based on their expected waste-heat signatures. Conservative upper limits are presented on the fraction of stars at G $\leq$ 21 that may potentially host non-reflective Dyson spheres that absorb 1 - 90$\%$ of the bolometric luminosity of their host stars and emit thermal waste-heat in the 100 - 1000 K range. Based on a sample of $\approx$ $2.7\mathrm{e}\,5$ stars within 100 pc, we find that a fraction less than $\approx$ $2\mathrm{e}\,-5$ could potentially host $\sim$300 K Dyson spheres at 90$\%$ completion. These limits become progressively weaker for less complete Dyson spheres due to increased confusion with naturally occurring sources of strong mid-infrared radiation, and also at larger distances, due to the detection limits of WISE. For the $\sim2.9\mathrm{e}\,8$ stars within 5 kpc in our Milky Way sample, the corresponding upper limit on the fraction of stars that could potentially be $\sim$300 K Dyson spheres at 90$\%$ completion is $\leq$ $8\mathrm{e}\,-4$.

We study the interaction of upstream ultra-low frequency (ULF) waves with collisionless shocks by analyzing the outputs of eleven 2D local hybrid simulation runs. Our simulated shocks have Alfv\'enic Mach numbers between 4.29-7.42 and their $\theta_{BN}$ angles are 15$^\circ$, 30$^\circ$, 45$^\circ$ and 50$^\circ$. The ULF wave foreshocks develop upstream of all of them. The wavelength and the amplitude of the upstream waves exhibit a complex dependence on the shock's M$_A$ and $\theta_{BN}$. The wavelength positively correlates with both parameters, with the dependence on $\theta_{BN}$ being much stronger. The amplitude of the ULF waves is proportional to the product of the reflected beam velocity and density, which also depend on M$_A$ and $\theta_{BN}$. The interaction of the ULF waves with the shock causes large-scale (several tens of upstream ion inertial lengths) shock rippling. The properties of the shock ripples are related to the ULF wave properties, namely thier wavelength and amplitude. In turn, the ripples have a large impact on the ULF wave transmission across the shock because they change local shock properties ($\theta_{BN}$, strength), so that different sections of the same ULF wave front encounter shock with different characteristics. Downstream fluctuations do not resemble the upstream waves in terms the wavefront extension, orientation or their wavelength. However some features are conserved in the Fourier spectra of downstream compressive waves that present a bump or flattening at wavelengths approximately corresponding to those of the upstream ULF waves. In the transverse downstream spectra these features are weaker.

We study the system of interacting axions and magnetic fields in the early universe after the quantum chromodynamics phase transition, when axions acquire masses. Both axions and magnetic fields are supposed to be spatially inhomogeneous. We derive the equations for the spatial spectra of these fields, which depend on conformal time. In case of the magnetic field, we deal with the spectra of the energy density and the magnetic helicity density. The evolution equations are obtained in the closed form within the mean field approximation. We choose the parameters of the system and the initial condition which correspond to realistic primordial magnetic fields and axions. The system of equations for the spectra is solved numerically. We study the evolution of the fields for different seed spectra of axions. We find that the mutual interaction between axions and magnetic fields is suppressed for the chosen realistic parameters of the system, that is in agreement with previous results where only the zero mode of axions was accounted for.

Ground-based gravitational-wave detectors like Cosmic Explorer can be tuned to improve their sensitivity at high or low frequencies by tuning the response of the signal extraction cavity. Enhanced sensitivity above 2 kHz enables measurements of the post-merger gravitational-wave spectrum from binary neutron star mergers, which depends critically on the unknown equation of state of hot, ultra-dense matter. Improved sensitivity below 500 Hz favors precision tests of extreme gravity with black hole ringdown signals and improves the detection prospects while facilitating an improved measurement of source properties for compact binary inspirals at cosmological distances. At intermediate frequencies, a more sensitive detector can better measure the tidal properties of neutron stars. We present and characterize the performance of tuned Cosmic Explorer configurations that are designed to optimize detections across different astrophysical source populations. These tuning options give Cosmic Explorer the flexibility to target a diverse set of science goals with the same detector infrastructure. We find that a 40 km Cosmic Explorer detector outperforms a 20 km in all key science goals other than access to post-merger physics. This suggests that Cosmic Explorer should include at least one 40 km facility.

Computational studies that use block-structured adaptive mesh refinement (AMR) approaches suffer from unnecessarily high mesh resolution in regions adjacent to important solution features. This deficiency limits the performance of AMR codes. In this work a novel hybrid adaptive multiresolution (HAMR) approach to AMR-based calculations is introduced to address this issue. The multiresolution (MR) smoothness indicators are used to identify regions of smoothness on the mesh where the computational cost of individual physics solvers may be decreased by replacing direct calculations with interpolation. We suggest an approach to balance the errors due to the adaptive discretization and the interpolation of physics quantities such that the overall accuracy of the HAMR solution is consistent with that of the MR-driven AMR solution. The performance of the HAMR scheme is evaluated for a range of test problems, from pure hydrodynamics to turbulent combustion.

Long range scalar fields with a coupling to matter appear to violate known bounds on gravitation in the solar system and the laboratory. This is evaded thanks to screening mechanisms. In this short review, we shall present the various screening mechanisms from an effective field theory point of view. We then investigate how they can and will be tested in the laboratory and on astrophysical and cosmological scales.

The Blandford and Znajek (BZ) split-monopole serves as an important theoretical example of the mechanism that can drive the electromagnetic extraction of energy from Kerr black holes. It is constructed as a perturbative low spin solution of Force Free Electrodynamics (FFE). Recently, Armas $et~al.$ put this construction on a firmer footing by clearing up issues with apparent divergent asymptotics. This was accomplished by resolving the behavior around the outer light surface, a critical surface of the FFE equations. Building on this, we revisit the BZ perturbative expansion, and extend the perturbative approach to higher orders in the spin parameter of the Kerr black hole. We employ matched-asymptotic-expansions and semi-analytic techniques to extend the split-monopole solution to the sixth-order in perturbation theory. The expansion necessarily includes novel logarithmic contributions in the spin parameter. We show that these higher order terms result in non-analytic contributions to the power and angular momentum output. In particular, we compute for the first time the perturbative contributions to the energy extraction at seventh- and eighth-order in the spin parameter. The resulting formula for the energy extraction improves the agreement with numerical simulations at finite spin. Moreover, we present a novel numerical procedure for resolving the FFE equations across the outer light surface, resulting in significantly faster convergence and greater accuracy, and extend this to higher orders as well. Finally, we include a general discussion of light surfaces as critical surfaces of the FFE equations.

Recent simulations show that very large electric and magnetic fields near the kilo Tesla strength will likely be generated by ultra-intense lasers at existing facilities over distances of hundreds of microns in underdense plasmas. Stronger ones are even expected in the future although some technical dificulties must be overcome. In addition, it has been shown that vacuum exhibits a peculiar non-linear behaviour in presence of high magnetic and electric field strengths. In this work we are interested in the analysis of thermodynamical contributions of vacuum to the energy density and pressure when radiation interacts with it in the presence of an external magnetic field. Using the Euler-Heisenberg formalism in the regime of weak fields i.e. smaller than critical Quantum Electrodynamics field strength values, we evaluate these magnitudes and analyze the highly anisotropic behaviour we find. Our work has implications for photon-photon scattering with lasers and astrophysically magnetized underdense systems far outside their surface where matter effects are increasingly negligible.