Papers for Friday, Aug 05 2022

We present the first detailed study comparing the populations of stellar streams in cosmological simulations to observed Milky Way dwarf galaxy streams. In particular, we compare streams identified around Milky Way analogs in the FIRE-2 simulations to stellar streams observed by the Southern Stellar Stream Spectroscopic Survey (S5). For an accurate comparison between the stream populations, we produce mock Dark Energy Survey (DES) observations of the FIRE streams and estimate the detectability of their tidal tails and progenitors. The number and stellar mass distributions of detectable stellar streams is consistent between observations and simulations. However, there are discrepancies in the distributions of pericenters and apocenters, with the detectable FIRE streams, on average, forming at larger pericenters (out to > 110 kpc) and surviving only at larger apocenters (> 40 kpc) than those observed in the Milky Way. We find that the population of high-stellar mass dwarf galaxy streams in the Milky Way is incomplete. Interestingly, a large fraction of the FIRE streams would only be detected as satellites in DES-like observations, since their tidal tails are too low-surface brightness to be detectable. We thus predict a population of yet-undetected tidal tails around Milky Way satellites, as well as a population of fully undetected low surface brightness stellar streams, and estimate their detectability with the Rubin Observatory. Finally, we discuss the causes and implications of the discrepancies between the stream populations in FIRE and the Milky Way, and explore future avenues for tests of satellite disruption in cosmological simulations.

The formation of close binaries has been an open question for decades. A large fraction of close binaries are in triple systems, suggesting that their formation may be associated with the Kozai-Lidov mechanism. However, this picture remains under debate because the configurations of many observed triples are unlikely to trigger the Kozai-Lidov mechanism. In this paper, we use the close binary samples, including eclipsing, spectroscopic, and astrometric binaries, from Gaia Data Release 3 to investigate the mysterious connection between inner binaries and their wide tertiaries. We show that the wide tertiary (at $10^3$-$10^4$ AU) fraction increases with decreasing orbital periods of the inner binaries. Compared to the field wide binary fraction, the wide tertiary fraction is $2.28\pm0.10$ times higher for eclipsing binaries (a median orbital period of $0.44$ day) and $0.65\pm0.03$ times lower for astrometric binaries (a median orbital period of $537$ days). The separation distribution of wide tertiaries is similar to wide binaries, with a tentative excess at $\sim10^4$ AU for tertiaries of eclipsing binaries. Using the $v$-$r$ angle distributions, we show that the wide tertiaries are consistent with isotropic orientations with respect to the inner binaries. The inferred eccentricity distribution of wide tertiaries is close to thermal ($f(e)\propto e$), similar to wide binaries at similar separations. The dynamical unfolding scenario is disfavored because it predicts highly eccentric wide tertiaries, which is inconsistent with our findings. For the Kozai-Lidov mechanism to be effective for wide tertiaries at $>10^3$ AU, the initial separations of the inner binaries need to be $>3$ AU. Future theoretical investigations are needed to explore the parameter space at these large initial separations and large tertiary separations.

The recent discovery of an unambiguous quiescent BH and main sequence O star companion in VFTS 243 opens the door to new constraints on theoretical stellar evolution and population models looking to reproduce the progenitors of black hole - black hole binaries. Here we show that the Binary Population and Spectral Synthesis fiducial models (BPASSv2.2.1) natively predict VFTS 243-like systems: We find that VFTS 243 likely originated from a binary system in a \about 15 day orbit with primary mass ranging from 40 to 50 \msol\, and secondary star with initial mass 24--25\msol. %BPASS systems with initial parameters similar to the ones inferred in the discovery paper result in a final system with an O star 10\msol more massive than indicated by the observations. Additionally we find that the death of the primary star must have resulted in a low energy explosion $E<10^{50}$ ergs. With a uniform prior we find that the kick velocity of the new-born black hole was $<33$ \kms (90 percent credible interval). The very low eccentricity reported for VFTS~243 and the subsequent conclusion by the authors that the SN kick must have been very small is in line with the peak in the posterior distribution between 0 and 5 \kms. Finally, the reduced Hobbs kick distribution commonly used in black hole population synthesis is strongly disfavoured, whereas the Bray kick with the most recent parameter calibration predicts 2 $\pm$ 3.5 \kms, which is very consistent with the posterior velocity distributions obtained for our matching VFTS 243-like models using a uniform kick prior.

The transformation of cold neutral intergalactic hydrogen into a highly ionized warm plasma marks the end of the cosmic dark ages and the beginning of the age of galaxies. The details of this process reflect the nature of the early sources of radiation and heat, the statistical characteristics of the large-scale structure of the Universe, the thermodynamics and chemistry of cosmic baryons, and the histories of star formation and black hole accretion. A number of massive data sets from new ground- and space-based instruments and facilities over the next decade are poised to revolutionize our understanding of primeval galaxies, the reionization photon budget, the physics of the intergalactic medium (IGM), and the fine-grained properties of hydrogen gas in the "cosmic web". In this review we survey the physics and key aspects of reionization-era modeling and describe the diverse range of computational techniques and tools currently available in this field.

Subsonic, compressive turbulence transfers energy to cosmic rays (CRs), a process known as non-resonant reacceleration. It is often invoked to explain observed ratios of primary to secondary CRs at $\sim \rm GeV$ energies, assuming wholly diffusive CR transport. However, such estimates ignore the impact of CR self-confinement and streaming. We study these issues in stirring box magnetohydrodynamic (MHD) simulations using Athena++, with field-aligned diffusive and streaming CR transport. For diffusion only, we find CR reacceleration rates in good agreement with analytic predictions. When streaming is included, reacceleration rates depend on plasma $\beta$. Due to streaming-modified phase shifts between CR and gas variables, they are slower than canonical reacceleration rates in low-$\beta$ environments like the interstellar medium (ISM) but remain unchanged in high-$\beta$ environments like the intracluster medium (ICM). We also quantify the streaming energy loss rate in our simulations. For sub-Alfv\'{e}nic turbulence, it is resolution-dependent (hence unconverged in large scale simulations) and heavily suppressed -- by an order of magnitude -- compared to the isotropic loss rate $v_{A} \cdot \nabla P_{\rm CR} / P_{\rm CR} \sim v_{A}/L_{0}$, due to misalignment between the mean field and isotropic CR gradients. Counterintuitively, and unlike acceleration efficiencies, CR losses are almost independent of magnetic field strength over $\beta \sim 1-100$ and are, therefore, not the primary factor behind lower acceleration rates when streaming is included. While this paper is primarily concerned with how turbulence affects CRs, in a follow-up paper (Bustard and Oh, in prep), we consider how CRs affect turbulence by diverting energy from the MHD cascade, altering the pathway to gas heating and steepening the turbulent power spectrum.

Laboratory testbeds are an integral part of conducting research and developing technology for high-contrast imaging and extreme adaptive optics. There are a number of laboratory groups around the world that use and develop resources that are imminently required for their operations, such as software and hardware controls. The CAOTIC (Community of Adaptive OpTics and hIgh Contrast testbeds) project is aimed to be a platform for this community to connect, share information, and exchange resources in order to conduct more efficient research in astronomical instrumentation, while also encouraging best practices and strengthening cross-team connections. In these proceedings, we present the goals of the CAOTIC project, our new website, and we focus in particular on a new approach to teaching version control to scientists, which is a cornerstone of successful collaborations in astronomical instrumentation.

Due to their low surface brightness, dwarf galaxies are particularly susceptible to tidal forces. The expected degree of disturbance depends on the assumed gravity law and whether they have a dominant dark halo. This makes dwarf galaxies useful for testing different gravity models. In this project, we use the Fornax Deep Survey (FDS) dwarf galaxy catalogue to compare the properties of dwarf galaxies in the Fornax Cluster with those predicted by the Lambda cold dark matter ($\Lambda$CDM) standard model of cosmology and Milgromian dynamics (MOND). We construct a test particle simulation of the Fornax system. We then use the MCMC method to fit this to the FDS distribution of tidal susceptibility $\eta$ (half-mass radius divided by theoretical tidal radius), the fraction of dwarfs that visually appear disturbed as a function of $\eta$, and the distribution of projected separation from the cluster centre. This allows us to constrain the $\eta$ value at which dwarfs should get destroyed by tides. Accounting for an $r'$-band surface brightness limit of 27.8 magnitudes per square arcsec, the required stability threshold is $\eta_{\textrm{destr}} = 0.25^{+0.07}_{-0.03}$ in $\Lambda$CDM and $ 1.88^{+0.85}_{-0.53}$ in MOND. The $\Lambda$CDM value is in tension with previous $\textit{N}$-body dwarf galaxy simulations, which indicate that $\eta_{\textrm{destr}} \approx 1$. Our MOND $\textit{N}$-body simulations indicate that $\eta_{\textrm{destr}} = 1.70 \pm 0.30$, which agrees well with our MCMC analysis of the FDS. We therefore conclude that the observed deformations of dwarf galaxies in the Fornax Cluster and the lack of low surface brightness dwarfs towards its centre are incompatible with $\Lambda$CDM expectations but well consistent with MOND.

We investigate the build-up of the stellar and dark matter haloes of Milky Way-like galaxies in cosmological hydrodynamics simulations. We show that the stellar halo is made up primarily of stars stripped from a small number of massive dwarfs, most of which are disrupted by the present day. The dark matter halo, on the other hand, is made up primarily of small unresolved subhaloes ($\lesssim 10^6$ M$_{\odot}$) and a ``smooth'' component consisting of particles which were never bound to a subhalo. Despite these differences, the massive dwarfs that make up the majority of the stellar halo also contribute a significant fraction of the dark matter. The stars and dark matter stripped from these dwarfs are related through their kinematics and this leaves imprints in the phase-space structure of the haloes. We examine the relation between the location of features, such as caustics, in the phase space of the stars and dark halo properties. We show that the ``edge'' of the stellar halo is a probe of dark matter halo mass and assembly history. The edges of Milky Way-mass galaxies should be visible at a surface brightness of 31-36 mag arcsec$^{-2}$.

Large-scale magnetic fields play a vital role in determining the angular momentum transport and in generating jets/outflows in the accreting systems, yet their origin remains poorly understood. We focus on radiatively inefficient accretion flows (RIAF) around the black holes, and conduct three-dimensional general-relativistic magnetohydrodynamic (GRMHD) simulations using the Athena++ code. We first re-confirm that the dynamo action alone cannot provide sufficient magnetic flux required to produce a strong jet. We next investigate the other possibility, where the large-scale magnetic fields are advected inward from external sources (e.g. the companion star in X-ray binaries, magnetized ambient medium in AGNs). Although the actual configuration of the external fields could be complex and uncertain, they are likely to be closed. As a first study, we treat them as closed field loops of different sizes, shapes and field strengths. Unlike earlier studies of flux transport, where magnetic flux is injected in the initial laminar flow, we injected the magnetic field loops in the quasi-stationary turbulent RIAF in inflow equilibrium and followed their evolution. We found that a substantial fraction ($\sim15\%-40\%$) of the flux injected at the large radii reaches the black hole with a weak dependence on the loop parameters except when the loops are injected at high latitudes, away from the mid-plane. Relatively high efficiency of flux transport observed in our study hints that a magnetically dominated RIAF, potentially a magnetically-arrested disk, might be formed relatively easily close to the black hole, provided that a source of the large-scale field exists at the larger radii.

We study the E and B mode polarisation of the cosmic microwave background (CMB) originating from the transverse peculiar velocity of free electrons, at second order in perturbation theory, during the reionisation and post-reionisation eras. Interestingly, the spectrum of this polarised kinetic Sunyaev-Zeldovich (SZ) effect can be decomposed into a blackbody part and a y-type distortion. The y-distortion part is distinguishable from the primary E and B modes and also the lensing B modes. Furthermore, it is also differentiable from the other y-type signals, such as the thermal SZ effect, which are unpolarised. We show that this signal is sensitive to the reionisation history, in particular to how fast reionisation happens. The E and B modes of y-type distortion provide a way to beat the cosmic variance of primary CMB anisotropies and are an independent probe of the cosmological parameters. The blackbody component of the pkSZ effect would be an important foreground for the primordial tensor modes for tensor to scalar ratio $r \lesssim 3\times10^{-5}$.

We introduce a model for the explicit evolution of interstellar dust in a cosmological galaxy formation simulation. We post-process a simulation from the Cosmic Reionization on Computers project (CROC, Gnedin 2014), integrating an ordinary differential equation for the evolution of the dust-to-gas ratio along pathlines in the simulation sampled with a tracer particle technique. This model incorporates the effects of dust grain production in asymptotic giant branch (AGB) star winds and supernovae (SN), grain growth due to the accretion of heavy elements from the gas phase of the interstellar medium (ISM), and grain destruction due to thermal sputtering in the high temperature gas of supernova remnants (SNRs). A main conclusion of our analysis is the importance of a carefully chosen dust destruction model, for which different reasonable parameterizations can predict very different values at the $\sim 100$ pc resolution of the ISM in our simulations. We run this dust model on the single most massive galaxy in a 10$h^{-1}$ co-moving Mpc box, which attains a stellar mass of $\sim 2\times10^9 M_{\odot}$ by $z=5$. We find that the model is capable of reproducing dust masses and dust-sensitive observable quantities broadly consistent with existing data from high-redshift galaxies. The total dust mass in the simulated galaxy is somewhat sensitive to parameter choices for the dust model, especially the timescale for grain growth due to accretion in the ISM. Consequently, observable quantities that can constrain galaxy dust masses at these epochs are potentially useful for placing constraints on dust physics.

We present the first Bayesian method for tomographic decomposition of the plane-of-sky orientation of the magnetic field with the use of stellar polarimetry and distance. This standalone tomographic inversion method presents an important step forward in reconstructing the magnetized interstellar medium (ISM) in 3D within dusty regions. We develop a model in which the polarization signal from the magnetized and dusty ISM is described by thin layers at various distances. Our modeling makes it possible to infer the mean polarization (amplitude and orientation) induced by individual dusty clouds and to account for the turbulence-induced scatter in a generic way. We present a likelihood function that explicitly accounts for uncertainties in polarization and parallax. We develop a framework for reconstructing the magnetized ISM through the maximization of the log-likelihood using a nested sampling method. We test our Bayesian inversion method on mock data taking into account realistic uncertainties from $Gaia$ and as expected for the optical polarization survey PASIPHAE according to the currently planned observing strategy. We demonstrate that our method is effective in recovering the cloud properties as soon as the polarization induced by a cloud to its background stars is higher than $\sim 0.1\%$, for the adopted survey exposure time and level of systematic uncertainty. Our method makes it possible to recover not only the mean polarization properties but also to characterize the intrinsic scatter, thus opening ways to characterize ISM turbulence and the magnetic field strength. Finally, we apply our method to an existing dataset of starlight polarization with known line-of-sight decomposition, demonstrating agreement with previous results and an improved quantification of uncertainties in cloud properties.

We present the conceptual design of the modular detector and readout system for the Cosmic Microwave Background Stage 4 (CMB-S4) ground-based survey experiment. CMB-S4 will map the cosmic microwave background (CMB) and the millimeter-wave sky to unprecedented sensitivity, using 500,000 superconducting detectors observing from Chile and Antarctica to map over 60 percent of the sky. The fundamental building block of the detector and readout system is a detector module package operated at 100 mK, which is connected to a readout and amplification chain that carries signals out to room temperature. It uses arrays of feedhorn-coupled orthomode transducers (OMT) that collect optical power from the sky onto dc-voltage-biased transition-edge sensor (TES) bolometers. The resulting current signal in the TESs is then amplified by a two-stage cryogenic Superconducting Quantum Interference Device (SQUID) system with a time-division multiplexer to reduce wire count, and matching room-temperature electronics to condition and transmit signals to the data acquisition system. Sensitivity and systematics requirements are being developed for the detector and readout system over a wide range of observing bands (20 to 300 GHz) and optical powers to accomplish CMB-S4's science goals. While the design incorporates the successes of previous generations of CMB instruments, CMB-S4 requires an order of magnitude more detectors than any prior experiment. This requires fabrication of complex superconducting circuits on over 10 square meters of silicon, as well as significant amounts of precision wiring, assembly and cryogenic testing.

Interstellar chemistry is important for galaxy formation, as it determines the rate at which gas can cool, and enables us to make predictions for observable spectroscopic lines from ions and molecules. We explore two central aspects of modelling the chemistry of the interstellar medium (ISM): (1) the effects of local stellar radiation, which ionises and heats the gas, and (2) the depletion of metals onto dust grains, which reduces the abundance of metals in the gas phase. We run high-resolution (400 M$_\odot$ per baryonic particle) simulations of isolated disc galaxies, from dwarfs to Milky Way-mass, using the FIRE galaxy formation models together with the CHIMES non-equilibrium chemistry and cooling module. In our fiducial model, we couple the chemistry to the stellar fluxes calculated from star particles using an approximate radiative transfer scheme, and we implement an empirical density-dependent prescription for metal depletion. For comparison, we also run simulations with a spatially uniform radiation field, and without metal depletion. Our fiducial model broadly reproduces observed trends in HI and H2 mass with stellar mass, and in line luminosity versus star formation rate for [CII] 158$\mu$m, [OI] 63$\mu$m, [OIII] 88$\mu$m, [NII] 122$\mu$m and H$\alpha$ 6563A. Our simulations with a uniform radiation field predict fainter luminosities, by up to an order of magnitude for [OIII] 88$\mu$m and H$\alpha$ 6563A, while ignoring metal depletion increases the luminosity of carbon and oxygen lines by a factor $\approx$2. However, the overall evolution of the galaxy is not strongly affected by local stellar fluxes or metal depletion, except in dwarf galaxies where the inclusion of local fluxes leads to weaker outflows and hence higher gas fractions.

We present the design, fabrication, and measured performance of metamaterial Anti-Reflection Cuttings (ARCs) for large-format alumina filters operating over more than an octave of bandwidth to be deployed on the Simons Observatory (SO). The ARC consists of sub-wavelength features diced into the optic's surface using a custom dicing saw with near-micron accuracy. The designs achieve percent-level control over reflections at angles of incidence up to 20$^\circ$. The ARCs were demonstrated on four 42 cm diameter filters covering the 75-170 GHz band and a 50 mm diameter prototype covering the 200-300 GHz band. The reflection and transmission of these samples were measured using a broadband coherent source that covers frequencies from 20 GHz to 1.2 THz. These measurements demonstrate percent-level control over reflectance across the targeted pass-bands and a rapid reduction in transmission as the wavelength approaches the length scale of the metamaterial structure where scattering dominates the optical response. The latter behavior enables the use of the metamaterial ARC as a scattering filter in this limit.

Binary stars in which oscillations can be studied in either or both components can provide powerful constraints on our understanding of stellar physics. The bright binary 12 Bo\"otis (12 Boo) is a particularly promising system because the primary is roughly 60 per cent brighter than the secondary despite being only a few per cent more massive. Both stars have substantial surface convection zones and are therefore, presumably, solar-like oscillators. We report here the first detection of solar-like oscillations and ellipsoidal variations in the TESS light curve of 12 Boo. Though the solar-like oscillations are not clear enough to unambiguously measure individual mode frequencies, we combine global asteroseismic parameters and a precise fit to the spectral energy distribution (SED) to provide new constraints on the properties of the system that are several times more precise than values in the literature. The SED fit alone provides new effective temperatures, luminosities and radii of $6115\pm45\,\mathrm{K}$, $7.531\pm0.110\,\mathrm{L}_\odot$ and $2.450\pm0.045\,\mathrm{R}_\odot$ for 12 Boo A and $6200\pm60\,\mathrm{K}$, $4.692\pm0.095\,\mathrm{L}_\odot$ and $1.901\pm0.045\,\mathrm{R}_\odot$ for 12 Boo B. When combined with our asteroseismic constraints on 12 Boo A, we obtain an age of $2.67^{+0.12}_{-0.16}\,\mathrm{Gyr}$, which is consistent with that of 12 Boo B.

With strong evidence of a common-spectrum stochastic process in the most recent datasets from the NANOGrav Collaboration, the European Pulsar Timing Array (PTA), Parkes PTA, and the International PTA, it is crucial to assess the effects of the several astrophysical and cosmological sources that could contribute to the stochastic gravitational wave background (GWB). Using the same dataset creation and injection techniques as in Pol et al. (2021), we assess the separability of multiple GWBs by creating single and multiple GWB source datasets. We search for these injected sources using Bayesian PTA analysis techniques to assess recovery and separability of multiple astrophysical and cosmological backgrounds. For a GWB due to supermassive black hole binaries and an underlying weaker background due to primordial gravitational waves with a GW energy density ratio of $\Omega_{\mathrm{PGW}}/\Omega_{\mathrm{SMBHB}} = 0.5$, the Bayes' factor for a second process exceeds unity at 17 years, and increases with additional data. At 20 years of data, we are able to constrain the spectral index and amplitude of the weaker GWB at this density ratio to a fractional uncertainty of 64% and 110%, respectively, using current PTA methods and techniques. Using these methods and findings, we outline a basic protocol to search for multiple backgrounds in future PTA datasets.

Does life exist outside our Solar System? A first step towards searching for life outside our Solar System is detecting life on Earth by using remote sensing applications. One powerful and unambiguous biosignature is the circular polarization resulting from the homochirality of biotic molecules and systems. We aim to investigate the possibility of identifying and characterizing life on Earth by using airborne spectropolarimetric observations from a hot air balloon during our field campaign in Switzerland, May 2022. In this work we present the optical-setup and the data obtained from aerial circular spectropolarimetric measurements of farmland, forests, lakes and urban sites. We make use of the well-calibrated FlyPol instrument that measures the fractionally induced circular polarization ($V/I$) of (reflected) light with a sensitivity of $<10^{-4}$. The instrument operates in the visible spectrum, ranging from 400 to 900 nm. We demonstrate the possibility to distinguish biotic from abiotic features using circular polarization spectra and additional broadband linear polarization information. We review the performance of our optical-setup and discuss potential improvements. This sets the requirements on how to perform future airborne spectropolarimetric measurements of the Earth's surface features from several elevations.

We use oxygen and argon abundances for planetary nebulae (PNe) with low internal extinction (progenitor ages of (>4.5 Gyr) and high extinction (progenitor ages <2.5 Gyr), as well as those of the H II regions, to constrain the chemical enrichment and star formation efficiency in the thin and thicker discs of M31. The argon element is produced in larger fraction by Type Ia supernovae (SNe) than oxygen. We find that the mean log(O/Ar) values of PNe as a function of their argon abundances, 12 + log(Ar/H), trace the inter-stellar matter (ISM) conditions at the time of birth of the M 31 disc PN progenitors. Thus the chemical enrichment and star formation efficiency information encoded in the [alpha/Fe] vs. [Fe/H] distribution of stars is also imprinted in the oxygen-to-argon abundance ratio log(O/Ar) vs. argon abundance for the nebular emissions of the different stellar evolution phases. We propose to use the log(O/Ar) vs. (12 + log(Ar/H)) distribution of PNe with different ages to constrain the star-formation histories of the parent stellar populations in the thin and thicker M31 discs. For the inner M31 disc (R_{GC} < 14 kpc), the chemical evolution model that reproduces the mean log(O/Ar) values as function of argon abundance for the high- and low-extinction PNe requires a second infall of metal poorer gas during a gas-rich (wet) satellite merger. In M31, the thin disc is younger and less radially extended, formed stars at a higher star formation efficiency, and had a faster chemical enrichment timescale than the more extended, thicker disc. Both the thin and thicker disc in M31 reach similar high argon abundances ( 12 + log(Ar/H) ) ~ 6.7. The chemical and structural properties of the thin/thicker discs in M31 are thus remarkably different from those determined for the Milky Way thin and thick discs.

In this paper, we investigate the feasibility of using subspace system identification techniques for estimating transient Structural-Thermal-Optical Performance (STOP) models of reflective optics. As a test case, we use a Newtonian telescope structure. This work is motivated by the need for the development of model-based data-driven techniques for prediction, estimation, and control of thermal effects and thermally-induced wavefront aberrations in optical systems, such as ground and space telescopes, optical instruments operating in harsh environments, optical lithography machines, and optical components of high-power laser systems. We estimate and validate a state-space model of a transient STOP dynamics. First, we model the system in COMSOL Multiphysics. Then, we use LiveLink for MATLAB software module to export the wavefront aberrations data from COMSOL to MATLAB. This data is used to test the subspace identification method that is implemented in Python. One of the main challenges in modeling and estimation of STOP models is that they are inherently large-dimensional. The large-scale nature of STOP models originates from the coupling of optical, thermal, and structural phenomena and physical processes. Our results show that large-dimensional STOP dynamics of the considered optical system can be accurately estimated by low-dimensional state-space models. Due to their low-dimensional nature and state-space forms, these models can effectively be used for the prediction, estimation, and control of thermally-induced wavefront aberrations. The developed MATLAB, COMSOL, and Python codes are available online.

We present a thorough search for signatures that would be suggestive of super-heavy $X$ particles decaying in the Galactic halo, in the data of the Pierre Auger Observatory. From the lack of signal, we derive upper limits for different energy thresholds above ${\gtrsim}10^8$\,GeV on the expected secondary by-product fluxes from $X$-particle decay. Assuming that the energy density of these super-heavy particles matches that of dark matter observed today, we translate the upper bounds on the particle fluxes into tight constraints on the couplings governing the decay process as a function of the particle mass. We show that instanton-induced decay processes allow us to derive a bound on the reduced coupling constant of gauge interactions in the dark sector: $\alpha_X \alt 0.09$, for $10^{9} \alt M_X/\text{GeV} < 10^{19}$. This upper limit on $\alpha_X$ is complementary to the non-observation of tensor modes in the cosmic microwave background in the context of Planckian-interacting massive particles for dark matter produced during the reheating epoch. Viable regions for this scenario to explain dark matter are delineated in several planes of the multidimensional parameter space that involves, in addition to $M_X$ and $\alpha_X$, the Hubble rate at the end of inflation, the reheating efficiency, and the non-minimal coupling of the Higgs with curvature.

Calcium monohydroxide radical (CaOH) is receiving an increasing amount of attention from the astrophysics community as it is expected to be present in the atmospheres of hot rocky super-Earth exoplanets as well as interstellar and circumstellar environments. Here, we report the high-resolution laboratory absorption spectroscopy on low-$J$ transitions in $\tilde{A}^2\Pi(0,0,0)-\tilde{X}^2\Sigma^+(0,0,0)$ band of buffer-gas-cooled CaOH. In total, 40 transitions out of the low-$J$ states were assigned, including 27 transitions which have not been reported in previous literature. The determined rotational constants for both ground and excited states are in excellent agreement with previous literature, and the measurement uncertainty for the absolute transition frequencies was improved by more than a factor of three. This will aid future interstellar, circumstellar, and atmospheric identifications of CaOH. The buffer-gas-cooling method employed here is a particularly powerful method to probe low-$J$ transitions and is easily applicable to other astrophysical molecules.

The Miniature X-ray Solar Spectrometer (MinXSS-1) CubeSat observed solar X-rays between 0.5 and 10 keV. A two-temperature, two-emission measure model is fit to each daily averaged spectrum. These daily average temperatures and emission measures are plotted against the corresponding daily solar 10.7 cm radio flux (F10.7) value and a linear correlation is found between each that we call the Schwab Woods Mason (SWM) model. The linear trends show that one can estimate the solar spectrum between 0.5 keV and 10 keV based on the F10.7 measurement alone. The cooler temperature component of this model represents the quiescent sun contribution to the spectra and is essentially independent of solar activity, meaning the daily average quiescent sun is accurately described by a single temperature (1.70 MK) regardless of solar intensity and only the emission measure corresponding to this temperature needs to be adjusted for higher or lower solar intensity. The warmer temperature component is shown to represent active region contributions to the spectra and varies between 5 MK to 6 MK. GOES XRS-B data between 1-8 Angstroms is used to validate this model and it is found that the ratio between the SWM model irradiance and the GOES XRS-B irradiance is close to unity on average. MinXSS-1 spectra during quiescent solar conditions have very low counts beyond around 3 keV. The SWM model can generate MinXSS-1 or DAXSS spectra at very high spectral resolution and with extended energy ranges to fill in gaps between measurements and extend predictions back to 1947.

We describe the design and initial results from a visible-light Lyot coronagraph for SCExAO/VAMPIRES. The coronagraph is comprised of four hard-edged, partially transmissive focal plane masks with inner working angles of 36 mas, 55 mas, 92 mas, and 129 mas, respectively. The Lyot stop is a reflective, undersized design with a geometric throughput of 65.7%. Our preliminary on-sky contrast is 1e-2 at 0.1" to 1e-4 at 0.75" for all mask sizes. The coronagraph was deployed in early 2022 and is available for open use.

We present an analysis of UV spectra of quasars at intermediate redshifts believed to belong to the extreme Population A (xA), aimed to estimate the chemical abundances of the broad line emitting gas. We follow the approach described in a previous work extending the sample to 42 sources. Our aim is to test the robustness of the analysis carried out previously, as well as to confirm the two most intriguing results of this investigation: evidence of very high solar metallicities, and deviation of the relative abundance of elements with respect to solar values. The basis of our analysis is multi-component fits in three regions of the spectra centered at 1900, 1550 and 1400 A in order to deblend the broad components of AlIII1860, CIII]1909, CIV1549, HeII1640, and SiIV1397+OIV]1402 and their blue excess. By comparing the observed flux ratios of these components with the same ratios predicted by photoionization code CLOUDY we found that the virialized gas (broad components) presents a metallicity Z higher than 10Z$_\odot$. For non-virialized clouds we derive a lower limit to the metallicity around $\sim$ 5Z$_\odot$ under the assumption of chemical composition proportional to the solar one, confirming the previous results. We especially rely on the ratios between metal lines and HeII1640. This allowed us to confirm systematic differences in the solar-scaled metallicity derived from the lines of Aluminium and Silicon, and of Carbon, with the first being a factor 2 higher. For luminous quasars accreting at high rates, high Z values are likely, but that Z scaled values are affected by the possible pollution due to highly-enriched gas associated with the circumnuclear star formation. The high-Z values suggest a complex process involving nuclear and circumnuclear star formation, interaction between nuclear compact objects and accretion disk, possibly with the formation of accretion-modified stars.

Natal supernova kicks, the linear momentum compact remnants receive during their formation, are an essential part of binary population synthesis (BPS) models. Although these kicks are well-supported by evidence, their underlying distributions and incorporation into BPS models is uncertain. In this work, we investigate the nature of natal kicks using a previously proposed analytical prescription where the strength of the kick is linearly proportional to the ejecta-remnant mass ratio. We vary the free parameters over large ranges of possible values, comparing these synthetic populations simultaneously against four constraints: the merger rate of compact binary neutron star (BNS) systems, the period-eccentricity distribution of galactic BNSs, the velocity distribution of single-star pulsars, and the likelihood for low-ejecta mass supernovae to produce low-velocity kicks. We find that different samples of the parameter space satisfy each tests, and only 1 per cent of the models satisfy all four constraints simultaneously. Although we cannot identify a single best kick model, we report $\alpha = 80 \pm 30$ km s$^{-1}$, $\beta = 0 \pm 20$ km s$^{-1}$ as the center of the region of the parameter space that fulfils all of our constraints, and expect $\beta \geq 0$ km s$^{-1}$ as a further constraint. We also suggest further observations that will enable future refinement of the kick model. A sensitive test for the kick model will be the redshift evolution of the BNS merger rate since this is effectively a direct measure of the delay-time distribution for mergers. For our best fitting values, we find that the peak of the BNS merger rate is the present-day.

As a major approach to looking for life beyond the Earth, the search for extraterrestrial intelligence (SETI) is committed to detecting technosignatures such as engineered radio signals that are indicative of technologically capable life. In this paper, we report a targeted SETI campaign employing an observation strategy named multi-beam coincidence matching (MBCM) at the Five-hundred-meter Aperture Spherical radio Telescope (FAST) towards 33 known exoplanet systems, searching for ETI narrow-band drifting signals across 1.05-1.45 GHz in two orthogonal linear polarization directions separately. A signal at 1140.604 MHz detected from the observation towards Kepler-438 originally peaked our interest because its features are roughly consistent with assumed ETI technosignatures. However, evidences such as its polarization characteristics are almost able to eliminate the possibility of an extraterrestrial origin. Our observations achieve an unprecedented sensitivity since the minimum equivalent isotropic radiated power (EIRP) we are able to detect reaches 1.48 x10^9 W.

We report analytical and numerical investigations of sub-ion-scale turbulence in low-beta plasmas, focusing on the spectral properties of the fluctuations and electron heating. In the isothermal limit, the numerical results strongly support a description of the turbulence as a critically-balanced Kolmogorov-like cascade of kinetic Alfv\'en wave fluctuations, as amended by Boldyrev & Perez (Astrophys. J. Lett. 758, L44 (2012)) to include intermittent effects. When the constraint of isothermality is removed (i.e., with the inclusion of electron kinetic physics), the energy spectrum is found to steepen due to electron Landau damping, which is enabled by the local weakening of advective nonlinearities around current sheets, and yields significant energy dissipation via a velocity-space cascade. The use of a Hermite-polynomial representation to express the velocity-space dependence of the electron distribution function allows us to obtain an analytical, lowest-order solution for the Hermite moments of the distribution, which is borne out by numerical simulations.

This work aimed to find the relationship between quasars' optical variability and spectral features to reveal the regularity behind the random variation. It is known that quasar's FeII/Hbeta flux ratio and equivalent width of [OIII]5007 are negatively correlated, called Eigenvector 1. In this work, we visualized the relationship between the position on this Eigenvector 1 (EV1) plane and how they had changed their brightness after ~10 years. We conducted three analyses using different quasar samples each. The first analysis showed the relation between their distributions on the EV1 plane and how much they had changed brightness, using 13,438 Sloan Digital Sky Survey quasars. This result shows how brightness changes later are clearly related to the position on the EV1 plane. In the second analysis, we plotted the sources reported as Changing-Look(State) Quasars on the EV1 plane. This result shows that the position on the EV1 plane corresponds activity level of each source, the bright or dim state of them are distributed on the opposite sides divided by the typical quasar distribution. In the third analysis, we examined the transition vectors on the EV1 plane using sources with multiple-epoch spectra. This result shows that the brightening and dimming sources move on the similar path and they turn into the position corresponding to the opposite activity level. We also found this trend is opposite to the empirical rule that RFeII positively correlated with the Eddington ratio, which has been proposed based on the trends of a large number of quasars. From all these analyses, it is indicated that quasars tend to oscillate between both sides of the distribution ridge on the EV1 plane; each of them corresponds to a dim state and a bright state. This trend in optical variation suggests that significant brightness changes, such as Changing-Look quasars, are expected to repeat.

A detailed comparison of HI and CO line cube data of the Galactic Center (GC) region from the archives is obtained. The central molecular zone (CMZ) is shown to be embedded in the HI disc (central HI zone, CHZ) of radius $\sim 320$ pc and vertical scale height $\sim 70$ pc. A radio continuum belt is shown to run parallel to molecular Arms I and II. The belt draws a double-infinity shape on the sky, connecting Sgr E ($l\sim -1^\circ.2)$, C, B1, B2 and Sgr D ($+1^\circ.2$), and is interpreted as a warping star-forming ring. The molecular Arms are closely associated with the HI arms on the longitude-velocity diagram (LVD), showing coherent rigid-body ridges. Due to the close relationship between HI and CO, the HI line absorption can be used to determine the Arms' position relative to Sgr A, B1, B2 and C.Combining the trigonometric data of proper motions of Sgr A$^*$ and maser sources of Sgr B2 as well as radial velocities, the 3D velocity vector of Sgr B2 is determined. From these analyses, the molecular Arm I with Sgr B2 is shown to be located in the near side of Sgr A$^*$, and Arm II with Sgr C in the other side, both composing a pair of symmetrical Arms around the GC. We present a possible 3D view of Sgr A through E and Arms I and II along with a parameter list.

Surrounding the Milky Way (MW) is the circumgalactic medium (CGM), an extended reservoir of hot gas that has significant implications for the evolution of the MW. We used the HaloSat all-sky survey to study the CGM's soft X-ray emission in order to better define its distribution and structure. We extend a previous HaloSat study of the southern CGM (Galactic latitude b < -30 deg) to include the northern CGM (b > 30 deg) and find evidence that at least two hot gas model components at different temperatures are required to produce the observed emission. The cooler component has a typical temperature of kT ~ 0.18 keV, while the hotter component has a typical temperature of kT ~ 0.7 keV. The emission measure in both the warm and hot components has a wide range (~ 0.005 - 0.03, ~ 0.0005 - 0.004 cm-6 pc respectively), indicating that the CGM is clumpy. A patch of relatively consistent CGM was found in the north, allowing for the CGM spectrum to be studied in finer detail using a stacked spectrum. The stacked spectrum is well described with a model including two hot gas components at temperatures of kT = 0.166 +/- 0.005 keV and kT = 0.69 +0.04 -0.05 keV. As an alternative to adding a hot component, a neon-enhanced single-temperature model of the CGM was also tested and found to have worse fit statistics and poor residuals.

The spectrum of the Wolf-Rayet (WR) star WR 63 contains spectral lines of two different O stars that show regular radial velocity (RV) variations with amplitudes of ~160 and ~225 km/s on a ~4.03 d period. The light-curve shows two narrow eclipses that are 0.2 mag deep on the same period as the RV changes. On the other hand, our data show no significant RV variations for the WR spectral lines. Those findings are compatible with WR 63 being a triple system composed of two non-interacting late O stars orbiting a WR star on a period larger than 1000 days. The amplitude of the WR spectral line-profile variability reaches 7-8% of the line intensity and seems related to a 0.04 mag periodic photometric variation. Large wind density structures are a possible origin of this variability, but our data are not sufficient to verify this. Our analysis shows that, should the three stars be bound, they would be coeval with an age of about 5.9+/-1.4 Myrs. The distance to the O stars is estimated to be 3.4+/-0.5 kpc. Their dynamical masses are 14.3+/-0.1 and 10.3+/-0.1 M_sol. Using rotating, single star evolutionary tracks, we estimate their initial masses to be 18+/-2 and 16+/-2 M_sol for the primary and the secondary, respectively. Regular spectral monitoring is required in the future to detect RV variations of the WR star that would prove that it is gravitationally bound to the close O+OB system and to determine its mass.

Quasi-periodic pulsations (QPPs) with various periods that originate in the underlying magnetohydrodynamic processes of the flaring structures are detected repeatedly in the solar flare emissions. We apply a 2D cellular automaton (CA) avalanche model to simulate QPPs as a result of a repetitive load/unload mechanism. We show that the frequent occurrence of magnetic reconnections in a flaring loop could induce quasi-periodic patterns in the detected emissions. We obtain that among 21070 simulated flares, 813 events endure over 50 seconds, scaled with the temporal resolution of the Yohkoh Hard X-ray Telescope, and about 70 percent of these rather long-lasting events exhibit QPPs. We also illustrate that the applied CA model provides a wide range of periodicities for QPPs. Furthermore, we observe the presence of multiple periods in nearly 50 percent of the cases applying the Lomb-Scargle periodogram. A lognormal distribution is fitted to the unimodal distribution of the periods as a manifestation of an underlying multiplicative mechanism that typifies the effect of the system's independent varying parameters. The global maximum of the periods' lognormal distribution is located at 29.29 seconds. We compare statistics of the simulated QPPs with parameters of the host flares and discuss the impacts of flare properties on QPPs' periods. Considering the intrinsic characteristic of CA models, namely the repetitive load/unload mechanism, and the obtained pieces of evidence, we suggest that CA models may generate QPPs. We also examine the applicability of the autoregressive integrated moving average models to describe the simulated and observational QPPs.

We investigate the effects of magnetic field configurations on the ionization rate by cosmic rays in protoplanetary disks. First, we consider cosmic-ray propagation from the interstellar medium (ISM) to the protoplanetary disks and showed that the cosmic-ray density around the disk should be 4 times lower than the ISM value. Then, we compute the attenuation of cosmic rays in protoplanetary disks. The magnetic fields in the disk are stretched to the azimuthal directions, and cosmic rays need to detour while propagating to the midplane. We show that the detouring effectively enhances the column density by about two orders of magnitudes. In the case of the disk around IM lup, this increases the ionization rate over an order of magnitude for $r\gtrsim\,100$ au. On the other hand, for $r\lesssim\,100$ au, the cosmic rays are shielded at the disk midplane while the ionization rate is also enhanced at $z\sim\,2H$. Our results are consistent with the recent ALMA observation that indicates the radial gradient in the cosmic-ray ionization rate. The elevated ionization rate in the outer radii of disks may activate the magnetorotational instability that was thought to be suppressed due to ambipolar diffusion. These results will have a strong influence on the dynamical and chemical evolutions of protoplanetary disks.

The Doppler shift of a radio signal is caused by the relative motion between the transmitter and receiver. The change in frequency of the signal over time is called drift rate. In the studies of radio SETI (Search for Extraterrestrial Intelligence), extraterrestrial narrow-band signals are expected to appear ``chirped'' since both the exoplanet and the Earth are moving. Such planet rotation and orbital revolution around the central star can cause a non-zero drift rate. Other relative motions between the transmitter and receiver, such as the gravitational redshift and galactic potential, are negligible. In this paper, we mainly consider the common cases that the drift rate is contributed by the rotations and orbits of the Earth and exoplanet in celestial mechanics perspective, and briefly discuss other cases different from the Earth-exoplanet one. We can obtain the expected pseudo-sinusoidal drifting result with long-term observations, shorter orbital periods of exoplanets. Exoplanets with higher orbital eccentricities can cause asymmetric drifting. The expected result should be intermittent pseudo-sinusoidal curves in long-term observations. The characteristics of pseudo-sinusoidal curves, as another new criterion for extraterrestrial signals, can be applied to long-term SETI re-observations in future research.

Companions at sub-arcsecond separation from young stars are difficult to image. However their presence can be inferred from the perturbations they create in the dust and gas of protoplanetary disks. Here we present a new interpretation of SPHERE polarised observations that reveal the previously detected inner spiral in the disk of HD 100546. The spiral coincides with a newly detected 12CO inner spiral and the previously reported CO emission Doppler-flip, which has been interpreted as the signature of an embedded protoplanet. Comparisons with hydrodynamical models indicate that this Doppler-flip is instead the kinematic counterpart of the spiral, which is likely generated by an inner companion inside the disk cavity.

Hierarchical mergers in a dense environment are one of the primary formation channels of intermediate-mass black hole (IMBH) binary system. We expect that the resulting massive binary system will exhibit mass asymmetry. The emitted gravitational-wave (GW) carry significant contribution from higher-order modes and hence complex waveform morphology due to superposition of different modes. Further, IMBH binaries exhibit lower merger frequency and shorter signal duration in the LIGO detector which increases the risk of them being misclassified as short-duration noisy glitches. Deep learning algorithms can be trained to discriminate noisy glitches from short GW transients. We present the $\mathtt{THAMES}$ -- a deep-learning-based end-to-end signal detection algorithm for GW signals from quasi-circular nearly edge-on, mass asymmetric IMBH binaries in advanced GW detectors. Our study shows that it outperforms matched-filter based $\mathtt{PyCBC}$ searches for higher mass asymmetric, nearly edge-on IMBH binaries. The maximum gain in the sensitive volume-time product for mass ratio $q \in (5, 10)$ is by a factor of 5.24 (2.92) against $\mathtt{PyCBC-IMBH}$ ($\mathtt{PyCBC-HM}$) search at a false alarm rate of 1 in 100 years. Compared to the broad $\mathtt{PyCBC}$ search this factor is $\sim100$ for the $q \in (10,18)$. One of the reasons for this leap in volumetric sensitivity is its ability to discriminate between signals with complex waveform morphology and noisy transients, clearly demonstrating the potential of deep learning algorithms in probing into complex signal morphology in the field of gravitational wave astronomy. With the current training set, $\mathtt{THAMES}$ slightly underperforms with respect to $\mathtt{PyCBC}$-based searches targeting intermediate-mass black hole binaries with mass ratio $q \in (5, 10)$ and detector frame total mass $M_T(1+z) \in (100,200)~M_\odot$.

Observational evidence is mounting regarding the population demographics of Massive Black Holes (MBHs), from the most massive cluster galaxies down to the dwarf galaxy regime. However, the progenitor pathways from which these central MBHs formed remain unclear. Here we report a potentially powerful observational signature of MBH formation in dwarf galaxies. We argue that a continuum in the mass spectrum of MBHs in (fossil) dwarf galaxies would be a unique signature of a heavy seed formation pathway. Under the robust assumption of initial fragmentation of the parent gas cloud resulting in a burst of heavy seed production, a significant fraction of these seeds will survive to the present day as off-nuclear MBHs with masses less than that of the central object. Motivated by the recent discovery of a MBH in the relatively low central density Leo I galaxy, we show that such a continuum in MBH seed masses should persist from the lightest black hole masses up to the mass of the central MBH in contrast to the light seeding scenario where no such continuum should exist. The detection of off-centered MBHs and a central MBH would represent strong evidence of a heavy seeding pathway.

We analyze the rest-frame near-UV and optical nebular spectra of three z > 7 galaxies from the Early Release Observations taken with the Near-Infrared Spectrograph (NIRSpec) on the James Webb Space Telescope (JWST). These three high-z galaxies show the detection of several strong-emission nebular lines, including the temperature-sensitive [O III] $\lambda$4363 line, allowing us to directly determine the nebular conditions and gas-phase abundances for O/H, C/O, Ne/O, and Fe/O. We derive O/H abundances and ionization parameters that are generally consistent with other recent analyses. The lowest-mass galaxy has a large O/H uncertainty, which as a significant effect on anchoring the mass-metallicity relationship (i.e., slope) and tests of its redshift evolution. We also detect the C III] $\lambda$$\lambda$1907,1909 emission in a z > 8 galaxy from which we determine the most distant C/O abundance to date. This valuable detection provides the first test of C/O redshift evolution out to high-redshift. For neon, we use the high-ionization [Ne III] $\lambda$3869 line to measure the first Ne/O abundances at z>7, finding no evolution in this $\alpha$-element ratio. To investigate the Fe abundance, we explore the tentative detection of weak [Fe II] and [Fe III] lines in a z>8 galaxy, which would indicate a rapid build up of metals. Importantly, we demonstrate that properly flux-calibrated and higher S/N spectra are crucial to robustly determine the abundance pattern in z>7 galaxies with NIRSpec/JWST.

Wolf-Rayet (WR) stars comprise a class of stars whose spectra are dominated by strong, broad emission lines that are associated with copious mass loss. In the massive-star regime, roughly 90% of the known WR stars are thought to have evolved off the main sequence. Dubbed classical WR (cWR) stars, these hydrogen-depleted objects represent a crucial evolutionary phase preceding core collapse into black holes, and offer a unique window into hot-star wind physics. Their formation is thought to be rooted in either intrinsic mass-loss or binary interactions. Results obtained from analyses using contemporary model atmospheres still fail to reconcile the derived properties of WR stars with predictions from stellar evolution. Importantly, stellar evolution models cannot reproduce the the bulk of cWR stars, a problem that becomes especially severe at subsolar metallicity. Next-generation model atmospheres and upcoming observational campaigns to hunt for undetected companions promise a venue for progress.

Low-surface-brightness galaxies (LSBs) contribute to a significant fraction of all the galaxies in the Universe. Ultra-diffuse galaxies (UDGs) form a subclass of LSBs that has attracted a lot of attention in recent years (although its definition may vary between studies). Although UDGs are found in large numbers in galaxy clusters, groups, and in the field, their formation and evolution are still very much debated. Using a comprehensive set of multiwavelength data from the NGVS (optical), VESTIGE (H$\alpha$ narrowband), and GUViCS (UV) surveys, we studied a sample of 64 diffuse galaxies and UDGs in the Virgo cluster to investigate their formation history. We analyzed the photometric colors and surface-brightness profiles of these galaxies and then compared them to models of galaxy evolution, including ram-pressure stripping (RPS) events to infer any possible strong interactions with the hot cluster gas in the past. While our sample consists mainly of red LSBs, which is typical in cluster environments, we found evidence of a color variation with the cluster-centric distance. Blue, HI-bearing, star-forming diffuse galaxies are found at larger distances from the cluster center than the rest of the sample. The comparison of our models with multifrequency observations suggests that most of the galaxies of the sample might have undergone a strong RPS event in their lifetime, on average 1.6 Gyr ago (with a large dispersion, and RPS still ongoing for some of them). This process resulted in the transformation of initially gas-rich diffuse blue galaxies into gas-poor and red ones that form the dominant population now, the more extreme UDGs having undergone the process in a more distant past on average. The RPS in dense environments could be one of the major mechanisms for the formation of the large number of quiescent UDGs we observe in galaxy clusters.

The presence of a strong magnetic field is a feature common to a significant fraction of degenerate stars, yet little is understood about field origin and evolution. New observational constraints from volume-limited surveys point to a more complex situation than a single mechanism valid for all stars. We show that in high-mass white dwarfs, which are probably the results of mergers, magnetic fields are extremely common and very strong, and appear immediately in the cooling phase. These fields may have been generated by a dynamo active during the merging. Lower mass white dwarfs, which are often the product of single star evolution, are rarely detectably magnetic at birth, but fields appear very slowly, and very weakly, in about a quarter of them. What we may see is an internal field produced in an earlier evolutionary stage that gradually relaxes to the surface from the interior. The frequency and strength of magnetic fields continue to increase to eventually rival those of highly-massive stars, particularly after the stars cool past the start of core crystallisation, an effect that could be responsible for a dynamo mechanism similar to the one that is active in Earth's interior.

Laser Guide Star (LGS) Shack-Hartmann (SH) wavefront sensors for next generation Extremely Large Telescopes (ELTs) require low-noise, large format (about 1Mpx), fast detectors to match the need for a large number of subapertures and a good sampling of the very elongated spots. One path envisaged to fulfill this need has been the adoption of CMOS detectors with a rolling shutter read-out scheme, that allows low read-out noise and fast readout time at the cost of image distortion due to the detector rows exposed in different moments. In this work we analyze the impact of the rolling shutter read-out scheme when used for LGS SH wavefront sensing of the Multiconjugate adaptive Optic Relay For ELT Observations (MORFEO, formerly known as MAORY) for ESO ELT; in particular, we focus on the impact on the adaptive optics correction of the distortion-induced aberrations created by the rolling exposure in the case of fast varying aberrations, like the ones coming from the LGS tilt jitter due to the up-link propagation of laser beams. We show that the LGS jitter-induced aberration for MORFEO can be as large as 100nm rms and we discuss possible mitigation strategies.

The Multiconjugate adaptive Optic Relay For ELT Observations (MORFEO, formerly known as MAORY) is the adaptive optics (AO) module for the Extremely Large Telescope (ELT) aimed at providing a 1 arcmin corrected field to the Multi-AO Imaging CamerA for Deep Observations (MICADO) and to a future client instrument. It should provide resolution close to the diffraction limit on a large portion of the sky and in a wide range of atmospheric conditions. Its ability to provide a flat wavefront must face the known aspect of the atmospheric turbulence and telescope environment, but also the final characteristic of a telescope still to be fully developed and built. In this work we focused on issues related to the segmentation of the telescope pupil (like low wind effect, residual phasing error at handover and control related issues), that could limit the system performance. MORFEO currently does not foresee a dedicated sensor to measure the phase step between adjacent mirror segments: in this work we study the possibility to use the low order wavefront sensors designed to sense and correct tip-tilt and focus as phasing sensors thanks to the linearized focal-plane technique (LIFT).

The MCAO Assisted Visible Imager and Spectrograph (MAVIS) is a new visible instrument for ESO Very Large Telescope (VLT). Its Adaptive Optics Module (AOM) must provide extreme adaptive optics correction level at low galactic latitude and high sky coverage at the galactic pole on the FoV of 30arcsec of its 4k x 4k optical imager and on its monolithic Integral Field Unit, thanks to 3 deformable mirrors (DM), 8 Laser Guide Stars (LGS), up to 3 Natural Guide Stars (NGS) and 11 Wave Front Sensors (WFS). A careful performance estimation is required to drive the design of this module and to assess the fulfillment of the system and subsystems requirements. Here we present the work done on this topic during the last year: we updated the system parameters to account for the phase B design and for more realistic conditions, and we produced a set of results from analytical and end-to-end simulations that should give a as complete as possible view on the performance of the system.

Context. B-type supergiants are versatile tools to address various astrophysical topics, ranging from stellar atmospheres over stellar and galactic evolution to the cosmic distance scale. Aims. A hybrid non-LTE approach - line-blanketed model atmospheres computed under the assumption of local thermodynamic equilibrium (LTE) in combination with line formation calculations that account for deviations from LTE - is tested for quantitative analyses of B-type supergiants with masses $M<30 M_{\odot}$, characterising a sample of 14 Galactic objects. Methods. Hydrostatic plane-parallel atmospheric structures and synthetic spectra computed with Kurucz's Atlas12 code together with the non-LTE line-formation codes Detail/Surface are compared to results from full non-LTE calculations with Tlusty, and the effects of turbulent pressure on the models are investigated. High-resolution spectra are analysed for atmospheric parameters, using Stark-broadened hydrogen lines and multiple metal ionisation equilibria, and for elemental abundances. Fundamental stellar parameters are derived by considering stellar evolution tracks and Gaia EDR3 parallaxes. Interstellar reddening towards the target stars is determined by matching model spectral energy distributions to observed ones. Results. Our hybrid non-LTE approach turns out to be equivalent to hydrostatic full non-LTE modelling for the deeper photospheric layers of the B-type supergiants considered. Turbulent pressure can become relevant for microturbulent velocities larger than 10 km s$^{-1}$. High precision and accuracy is achieved for all derived parameters by bringing multiple indicators to agreement simultaneously. Abundances for chemical species (He, C, N, O, Ne, Mg, Al, Si, S, Ar, Fe) are derived with uncertainties of 0.05 to 0.10 dex. The derived ratios N/C vs. N/O tightly follow the predictions from Geneva stellar evolution models.

MagAO-X system is a new adaptive optics for the Magellan Clay 6.5m telescope. MagAO-X has been designed to provide extreme adaptive optics (ExAO) performance in the visible. VIS-X is an integral-field spectrograph specifically designed for MagAO-X, and it will cover the optical spectral range (450 - 900 nm) at high-spectral (R=15.000) and high-spatial resolution (7 mas spaxels) over a 0.525 arsecond field of view. VIS-X will be used to observe accreting protoplanets such as PDS70 b and c. End-to-end simulations show that the combination of MagAO-X with VIS-X is 100 times more sensitive to accreting protoplanets than any other instrument to date. VIS-X can resolve the planetary accretion lines, and therefore constrain the accretion process. The instrument is scheduled to have its first light in Fall 2021. We will show the lab measurements to characterize the spectrograph and its post-processing performance.

GRB 171205A is a low-luminosity, long-duration gamma-ray burst (GRB) associated with SN 2017iuk, a broad-line type Ic supernova (SN). It is consistent with being formed in the core-collapse of a single CO star, or in a widely separated binary, which we have called the Binary driven Hypernova (BdHN) of type III. The core-collapse of the CO star forms a newborn NS ($\nu$NS) and the SN explosion. Fallback accretion transfers mass and angular momentum to the $\nu$NS. The accretion energy injected into the expanding stellar layers powers the prompt emission. The multiwavelength power-law afterglow is explained by the synchrotron radiation of electrons in the SN ejecta, powered by energy injected by the spinning $\nu$NS. We calculate the amount of mass and angular momentum gained by the $\nu$NS, as well as the $\nu$NS rotational evolution. The $\nu$NS spins up to a period of $58$ ms, then releases its rotational energy powering the synchrotron emission of the afterglow. The paucity of the $\nu$NS spin explains the low-luminosity characteristic and that the optical emission of the SN from the nickel radioactive decay outshines the optical emission from the synchrotron radiation. From the $\nu$NS evolution, we infer that the SN explosion had to occur at most $7.36$ h before the GRB trigger. Therefore, for the first time, the analysis of the GRB data leads to the time of occurrence of the associated SN explosion, setting a stringent delay time between the neutrino emission associated with the SN and the electromagnetic emission of the GRB event.

Observations of the Cosmic Microwave Background rely on cryogenic instrumentation with cold detectors, readout, and optics providing the low noise performance and instrumental stability required to make more sensitive measurements. It is therefore critical to optimize all aspects of the cryogenic design to achieve the necessary performance, with low temperature components and acceptable system cooling requirements. In particular, we will focus on our use of thermal filters and cold optics, which reduce the thermal load passed along to the cryogenic stages. To test their performance, we have made a series of in situ measurements while integrating the third receiver for the BICEP Array telescope. In addition to characterizing the behavior of this receiver, these measurements continue to refine the models that are being used to inform design choices being made for future instruments.

Context. Plasma inflows accompany a variety of processes in the solar atmosphere such as heating of coronal loops and formation of prominences. Aims. We model a stratified solar atmosphere, within which a simulated prominence thread experiences density accumulation via a plasma inflow designed to mimic the formation process. We aim to investigate the interaction of such a system with torsional perturbations, and the possible consequences. Methods. The linearised equations of motion and induction are integrated to analyse the spatial and temporal evolution of torsional perturbations that are randomly driven at the photospheric footpoints. Results. Our results demonstrate that magnetic threads will experience twist amplification. Different sources and sinks of energy and the corresponding amplification mechanisms are identified. Threads reaching chromospheric heights are most susceptible to magnetic twisting with the maximum twist occurring near their footpoints. The amplifying twists are associated with a standing wave behaviour along the simulated threads. Conclusions. Our work suggests that torsional perturbations may be amplified within prominence threads, with strong magnetic twists forming at the footpoints. The amplification process is facilitated by small length scales in the background magnetic field. On the other hand, a small length scale in the background density inhibits growth. Possible consequences of the amplified twists, including their role in supporting the dense plasma within a prominence structure are discussed.

We present photometric and spectroscopic data on three extragalactic luminous red novae (LRNe): AT2018bwo, AT2021afy, and A2021blu. AT2018bwo was discovered in NGC45 (at about 6.8 Mpc) a few weeks after the outburst onset. During the monitoring period, the transient reached a peak luminosity of 10^40 erg/s. AT2021afy, hosted by UGC10043 (~49.2 Mpc), showed a double-peaked light curve, with an early maximum reaching a luminosity of 2.1x10^41 erg/s. Finally, for AT2021blu in UGC5829 (~8.6 Mpc), the pre-outburst phase was well-monitored by several photometric surveys, and the object showed a slow luminosity rise before the outburst. The light curve of AT2021blu was sampled with an unprecedented cadence until the object disappeared behind the Sun, and was then recovered at late phases. The light curve of LRN AT2021blu shows a double peak, with a prominent early maximum reaching a luminosity of 6.5x10^40 erg/s, half that of AT2021afy. The spectra of AT2021afy and AT 2021blu display the expected evolution for LRNe: a blue continuum dominated by prominent Balmer lines in emission during the first peak, and a redder continuum consistent with that of a K-type star with narrow absorption metal lines during the second, broad maximum. The spectra of AT2018bwo are markedly different, with a very red continuum dominated by broad molecular features in absorption. As these spectra closely resemble those of LRNe after the second peak, probably AT 2018bwo was discovered at the very late evolutionary stages. This would explain its fast evolution and the spectral properties compatible with that of an M-type star. From the analysis of deep frames of the LRN sites years before the outburst, and considerations of the light curves, the quiescent progenitor systems of the three LRNe were likely massive, with primaries ranging from about 13Mo for AT2018bwo, to 13-18Mo for AT 2021blu, and over 40Mo for AT2021afy.

The current state of far-infrared astronomy drives the need to develop compact, sensitive spectrometers for future space and ground-based instruments. Here we present details of the $\rm \mu$-Spec spectrometers currently in development for the far-infrared balloon mission EXCLAIM. The spectrometers are designed to cover the $\rm 555 - 714\ \mu$m range with a resolution of $\rm R\ =\ \lambda / \Delta\lambda\ =\ 512$ at the $\rm 638\ \mu$m band center. The spectrometer design incorporates a Rowland grating spectrometer implemented in a parallel plate waveguide on a low-loss single-crystal Si chip, employing Nb microstrip planar transmission lines and thin-film Al kinetic inductance detectors (KIDs). The EXCLAIM $\rm \mu$-Spec design is an advancement upon a successful $\rm R = 64\ \mu$-Spec prototype, and can be considered a sub-mm superconducting photonic integrated circuit (PIC) that combines spectral dispersion and detection. The design operates in a single $M{=}2$ grating order, allowing one spectrometer to cover the full EXCLAIM band without requiring a multi-order focal plane. The EXCLAIM instrument will fly six spectrometers, which are fabricated on a single 150 mm diameter Si wafer. Fabrication involves a flip-wafer-bonding process with patterning of the superconducting layers on both sides of the Si dielectric. The spectrometers are designed to operate at 100 mK, and will include 355 Al KID detectors targeting a goal of NEP ${\sim}8\times10^{-19}$ $\rm W/\sqrt{Hz}$. We summarize the design, fabrication, and ongoing development of these $\rm \mu$-Spec spectrometers for EXCLAIM.

We compute 1-loop corrections to the redshift space galaxy power spectrum in cosmologies containing additional scales, and hence kernels different from Einstein-de Sitter (EdS). Specifically, our method is tailored for cosmologies in the presence of massive neutrinos and some modified gravity models; in this article we concentrate on the former case. The perturbative kernels have contributions that we notice appear either from the logarithmic growth factor $f(k,t)$, which is scale-dependent because of the neutrino free-streaming, or from the failure of the commonly used approximation $f^2=\Omega_m$. The latter contributions make the computation of loop corrections quite slow, precluding full-shape analyses for parameter estimation. However, we identify that the dominant pieces of the kernels come from the growth factor, allowing us to simplify the kernels but retaining the characteristic free-streaming scale introduced by the neutrinos' mass. Moreover, with this simplification one can exploit FFTLog methods to speed up the computations even more. We validate our analytical modeling and numerical method with halo catalogs extracted from the Quijote simulations finding good agreement with the, a priori, known cosmological parameters. We make public our Python code FOLPS$\nu$ to compute the redshift space power spectrum in a fraction of second. Code available at https://github.com/henoriega/FOLPS-nu.

$JWST$'s first glimpse of the $z>10$ Universe has yielded a surprising abundance of luminous galaxy candidates. Here we present the most extreme of these systems: CEERS-1749. Based on $0.6-5\mu$m photometry, this strikingly luminous ($\approx$26 mag) galaxy appears to lie at $z\approx17$. This would make it an $M_{\rm{UV}}\approx-22$, $M_{\rm{\star}}\approx5\times10^{9}M_{\rm{\odot}}$ system that formed a mere $\sim220$ Myrs after the Big Bang. The implied number density of this galaxy and its analogues challenges virtually every early galaxy evolution model that assumes $\Lambda$CDM cosmology. However, there is strong environmental evidence supporting a secondary redshift solution of $z\approx5$: all three of the galaxy's nearest neighbors at $<2.5$" have photometric redshifts of $z\approx5$. Further, we show that CEERS-1749 may lie in a $z\approx5$ protocluster that is $\gtrsim5\times$ overdense compared to the field. Intense line emission at $z\approx5$ from a quiescent galaxy harboring ionized gas, or from a dusty starburst, may provide satisfactory explanations for CEERS-1749's photometry. The emission lines at $z\approx5$ conspire to boost the $>2\mu$m photometry, producing an apparent blue slope as well as a strong break in the SED. Such a perfectly disguised contaminant is possible only in a narrow redshift window ($\Delta z\lesssim0.1$), implying that the permitted volume for such interlopers may not be a major concern for $z>10$ searches, particularly when medium-bands are deployed. If CEERS-1749 is confirmed to lie at $z\approx5$, it will be the highest-redshift quiescent galaxy, or one of the lowest mass dusty galaxies of the early Universe detected to-date. Both redshift solutions of this intriguing galaxy hold the potential to challenge existing models of early galaxy evolution, making spectroscopic follow-up of this source critical.

We present a set of coupled Boltzmann equations describing the intensity and polarisation Stokes parameters of the SGWB, including collision terms which account for gravitational Compton scattering off of massive objects. This set resembles that for the CMB Stokes parameters, but the different spin nature of the gravitational radiation and the physics involved in the scattering process determine crucial differences. In the case of gravitational Compton scattering, all the SGWB intensity multipoles $\ell$ (with $m=0$ in the case of scalar metric perturbations alone) are scattered out, therefore producing outgoing intensity anisotropies of any order $\ell$ if they are present in the incoming radiation. SGWB linear polarisation can be generated from unpolarised anisotropic radiation only with $m=\pm 4$, which requires at least an hexadecapole anisotropy ($\ell\ge 4$) in the incoming intensity. We confirm the contribution of the gravitational Compton scattering to the SGWB anisoptropies is extremely small for collisions with massive compact objects (BH and SMBH) in the frequency range of current and upcoming surveys. However, we stress that the system of coupled Boltzmann equations presented here provides an accurate estimate of the total amount of anisotropies generated by multiple SGWB scattering processes off of massive objects, as well as the interplay between polarisation and intensity, during the GW propagation across the LSS of the universe. These results will be useful for the full treatment of the astrophysical SWGB anisotropies in view of upcoming gravitational waves detectors.

Moduli potential loses its minima due to external energy sources of inflaton energy density or radiation produced at the end of inflation. But, the non-existence of minima does not necessarily mean destabilization of moduli. In fact, the destabilization of moduli is always dependent on the initial field values of the fields. In this work, we study carefully how the effects of reheating ease the problem of moduli destabilization. The associated time scale to produce the thermal bath allows a larger initial field range to stabilize the field. Contrary to the usual notion, the allowed initial field range is larger for higher temperatures when the effective potential is of a run-away nature. This eases the moduli destabilization problem for heavy mass moduli. For low mass moduli ($\lesssim$ 30 TeV), the allowed field range still causes the cosmological moduli problem by violating the BBN constraints unless its initial abundance is suppressed.

In this work, we invent the Anyon Cavity Resonator. The resonator is based on twisted hollow structures, which allow select resonant modes to exhibit non-zero helicity. Depending on the cross-section of the cavity, the modes have more general symmetry than what has been studied before. For example, with no twist, the mode is the form of a boson, while with a $180^{o}$ twist the symmetry is in the form of a fermion. We show that the generally twisted resonator is in the form of an anyon. The non-zero helicity couples the mode to axions, and we show in the upconversion limit the mode couples to ultra-light axions within the bandwidth of the resonator. The coupling adds amplitude modulated sidebands and allows a simple sensitive way to search for ultra-light axions using only a single mode within the resonator's bandwidth.

We explore potential uses of physics formulated in Kleinian (i.e., $2+2$) signature spacetimes as a tool for understanding properties of physics in Lorentzian (i.e., $3+1$) signature. Much as Euclidean (i.e., $4+0$) signature quantities can be used to formally construct the ground state wavefunction of a Lorentzian signature quantum field theory, a similar analytic continuation to Kleinian signature constructs a state of low particle flux in the direction of analytic continuation. There is also a natural supersymmetry algebra available in $2+2$ signature, which serves to constrain the structure of correlation functions. Spontaneous breaking of Lorentz symmetry can produce various $\mathcal{N} = 1/2$ supersymmetry algebras that in $3 + 1$ signature correspond to non-supersymmetric systems. We speculate on the possible role of these structures in addressing the cosmological constant problem.

Gravitational-wave parameter estimation for compact binary signals typically relies on sequential estimation of the properties of the detector Gaussian noise and of the binary parameters. This procedure assumes that the noise variance, expressed through its power spectral density, is perfectly known in advance. We assess the impact of this approximation on the estimated parameters by means of an analysis that simultaneously estimates the noise and compact binary parameters, thus allowing us to marginalize over uncertainty in the noise properties. We compare the traditional sequential estimation method and the new full marginalization method using events from the GWTC-3 catalog. We find that the recovered signals and inferred parameters agree to within their statistical measurement uncertainty. At current detector sensitivities, uncertainty about the noise power spectral density is a subdominant effect compared to other sources of uncertainty.

Energy dissipation in collisionless plasmas is a longstanding fundamental physics problem. Although it is well known that magnetic reconnection and turbulence are coupled and transport energy from system-size scales to sub-proton scales, the details of the energy distribution and energy dissipation channels remain poorly understood. Especially, the energy transfer and transport associated with three dimensional (3D) small-scale reconnection that occurs as a consequence of a turbulent cascade is unknown. We use an explicit fully kinetic particle-in-cell code to simulate 3D small scale magnetic reconnection events forming in anisotropic and Alfv\'enic decaying turbulence. We identify a highly dynamic and asymmetric reconnection event that involves two reconnecting flux ropes. We use a two-fluid approach based on the Boltzmann equation to study the spatial energy transfer associated with the reconnection event and compare the power density terms in the two-fluid energy equations with standard energy-based damping, heating and dissipation proxies. Our findings suggest that the electron bulk flow transports thermal energy density more efficiently than kinetic energy density. Moreover, in our turbulent reconnection event, the energy-density transfer is dominated by plasma compression. This is consistent with turbulent current sheets and turbulent reconnection events, but not with laminar reconnection.

Detectable electromagnetic counterparts to gravitational waves from compact binary mergers can be produced by outflows from the black hole-accretion disk remnant during the first ten seconds after the merger. Two-dimensional axisymmetric simulations with effective viscosity remain an efficient and informative way to model this late-time post-merger evolution. In addition to the inherent approximations of axisymmetry and modeling turbulent angular momentum transport by a viscosity, previous simulations often make other simplifications related to the treatment of the equation of state and turbulent transport effects. In this paper, we test the effect of these modeling choices. By evolving with the same viscosity the exact post-merger initial configuration previously evolved in Newtonian viscous hydrodynamics, we find that the Newtonian treatment provides a good estimate of the disk ejecta mass but underestimates the outflow velocity. We find that the inclusion of heavy nuclei causes a notable increase in ejecta mass. An approximate inclusion of r-process effects has a comparatively smaller effect, except for its designed effect on the composition. Diffusion of composition and entropy, modeling turbulent transport effects, has the overall effect of reducing ejecta mass and giving it a speed with lower average and more tightly-peaked distribution. Also, we find significant acceleration of outflow even at distances beyond 10,000\,km, so that thermal wind velocities only asymptote beyond this radius and at somewhat higher values than previously reported.

We argue that Nash theory, a quadratic theory of Gravity, can describe a late-time cosmic acceleration without any exotic matter or cosmological constant. The observational viability of an exact cosmological solution of Nash theory is adjudged using a Markov chain Monte Carlo simulation and JLA+OHD+BAO data sets. Departures from standard {\Lambda}CDM cosmology are noted and analyzed. We prove that the Nash vacuum dynamics is equivalent to the dynamics of an Einstein vacuum plus a self-interacting Higgs scalar field, only if a mild evolution of the Higgs Vacuum Expectation Value is allowed. This leads to variations in the mass scales of fundamental particles and the fine structure constant. The variations are found to fit nicely with the analysis of molecular absorption spectra from a series of Quasars.

We have already shown how a modified Friedmann equation, originating from a model of the Universe built from a certain $W_3$ algebra, is able to explainthe difference between the Hubble constants extracted from CMB data and from local measurements. In this article we show that the same model also describes aspects of the large scale structure of the Universe well.

In this paper, we consider a gravitational action containing a combination of the Ricci scalar, $R$, and the topological Gauss--Bonnet term, $G$. Specifically, we study the cosmological features of a particular class of modified gravity theories selected by symmetry considerations, namely the $f(R,G)= R^n G^{1-n}$ model. In the context of a spatially flat, homogeneous and isotropic background, we show that the currently observed acceleration of the Universe can be addressed through geometry, hence avoiding \emph{de facto} the shortcomings of the cosmological constant. We thus present a strategy to numerically solve the Friedmann equations in presence of pressureless matter and obtain the redshift behavior of the Hubble expansion rate. Then, to check the viability of the model, we place constraints on the free parameters of the theory by means of a Bayesian Monte Carlo method applied to late-time cosmic observations. Our results show that the $f(R,G)$ model is capable of mimicking the low-redshift behavior of the standard $\Lambda$CDM model, though substantial differences emerge when going toward high redshifts, leading to the absence of a standard matter-dominated epoch. Finally, we investigate the energy conditions and show that, under suitable choices for the values of the cosmographic parameters, they are all violated when considering the mean value of $n$ obtained from our analysis, as occurs in the case of a dark fluid.

In this paper, an effective field theory for proton superconductor (SC) interacting with neutron superfluid (SF), both with scalar order parameters, is developed and applied to the surface energy (SE) of a magnetized SC body in neutron stars (NS). Essentially, the SE studied here differs from the nuclear SE: here, the proton SF density decays to zero while the total proton density is constant across the surface. Interactions between the condensates are parameterized phenomenologically and their effects determined from calculations of a planar SE as the ranges of parameters are varied. The critical Ginzburg-Landau (GL) parameter $\kappa_c$ which renders the SE equal to zero is found analytically by noting that in a system with vanishing SE the thermodynamic critical MF is equivalent to the upper critical MF. In the case of weak coupling, $\kappa_c$ is shown to be a linear function of SF-SF density-density coupling, in agreement with the earlier results based on asymptotic intervortex interactions. Numerical simulations corroborate our analytical predictions. Coupling due to the mixed term arising from a scalar product of gradients of the SF densities, which had been considered in the earlier literature, is seen to have practically no effect on the superconductivity type. However, this coupling does produce a frozen wave packet of the SF neutron density localized at the surface. It is shown that the leading contribution from the gradient coupling arises from a novel mixed quantum pressure term, but still does not affect the planar SE. The present calculations provide an initial map of superconductivity types in the phenomenological effective field theory and will serve as a landmark for future studies, which require microscopic calculations of the coupling parameters introduced here phenomenologically.

The quadrupole normal-mode oscillation frequency $f_{n}$ of multiquark star are computed for $n=1-5$. At the transition from low to high density multiquark in the core region, the first 2 modes jump to larger values, a distinctive signature of the presence of the high-density core. When the star oscillation couples with spacetime, gravitational waves~(GW) will be generated and the star will undergo damped oscillation. The quasinormal modes~(QNMs) of the oscillation are computed using two methods, direct scan and WKB, for QNMs with small and large imaginary parts respectively. The small imaginary QNMs have frequencies $1.5-2.6$ kHz and damping times $0.19-1.7$ secs for multiquark star with mass $M=0.6-2.1 M_{\odot}$~(solar mass). The WKB QNMs with large imaginary parts have frequencies $5.98-9.81$ kHz and damping times $0.13-0.46$ ms for $M\simeq 0.3-2.1 M_{\odot}$. They are found to be the fluid $f-$modes and spacetime curvature $w-$modes respectively.

We apply common gravitational wave inference procedures on binary black hole merger waveforms from a theory of gravity beyond general relativity. We consider dynamical Chern-Simons gravity, a modified theory of gravity with origins in string theory and loop quantum gravity. This theory introduces an additional parameter $\ell$, corresponding to the length-scale below which quantum gravity effects become important. We simulate data based on numerical relativity waveforms produced under this theory that differ from the predictions of general relativity in the strongly nonlinear merger regime. We consider a system with parameters similar to GW150914 with different values of $\ell$ and signal-to-noise ratios. We perform two analyses of the simulated data. The first is a template-based analysis that uses waveforms derived under general relativity and allows us to identify degeneracies between waveforms predicted by the two theories of gravity. The second is a morphology-independent analysis based on BayesWave that does not assume that the signal is consistent with general relativity. Under the BayesWave analysis, the simulated signals can be faithfully reconstructed. However, waveform models derived under general relativity are unable to fully mimic the simulated modified-gravity signals and such a deviation would be identifiable with existing inference tools. Depending on the magnitude of the deviation $\ell$, we find that the templated analysis can under perform the morphology-independent analysis in fully recovering simulated beyond-GR waveforms even for achievable signal-to-noise ratios $\gtrsim 20{-}30$.