Papers for Thursday, May 19 2022

Using spatially resolved H-alpha emission line maps of star-forming galaxies, we study the evolution of gradients in galaxy assembly over a wide range in redshift ($0.5<z<1.7$). Our $z\sim0.5$ measurements come from deep Hubble Space Telescope WFC3 G102 grism spectroscopy obtained as part of the CANDELS Lyman-alpha Emission at Reionization (CLEAR) Experiment. For star-forming galaxies with Log$(M_{*}/\mathrm{M}_{\odot})\geqslant8.96$, the mean H-alpha effective radius is $1.2\pm0.1$ times larger than that of the stellar continuum, implying inside-out growth via star formation. This measurement agrees within $1\sigma$ with those measured at $z\sim1$ and $z\sim1.7$ from the 3D-HST and KMOS-3D surveys respectively, implying no redshift evolution. However, we observe redshift evolution in the stellar mass surface density within 1 kiloparsec ($\Sigma_\mathrm{1kpc}$). Star-forming galaxies at $z\sim0.5$ with a stellar mass of Log$(M_{*}/\mathrm{M}_{\odot})=9.5$ have a ratio of $\Sigma_\mathrm{1kpc}$ in H-alpha relative to their stellar continuum that is lower by $(19\pm2)\%$ compared to $z\sim1$ galaxies. $\Sigma_{1\mathrm{kpc, H}\alpha}$/$\Sigma_{1\mathrm{kpc,Cont}}$ decreases towards higher stellar masses. The majority of the redshift evolution in $\Sigma_{1\mathrm{kpc,H}\alpha}$/$\Sigma_{1\mathrm{kpc,Cont}}$ versus stellar mass stems from the fact that Log($\Sigma_{1\mathrm{kpc, H}\alpha}$) declines twice as much as Log($\Sigma_{1\mathrm{kpc, Cont}}$) from $z\sim 1$ to 0.5 (at a fixed stellar mass of Log$(M_{*}/\mathrm{M}_{\odot})=9.5$). By comparing our results to the TNG50 cosmological magneto-hydrodynamical simulation, we rule out dust as the driver of this evolution. Our results are consistent with inside-out quenching following in the wake of inside-out growth, the former of which drives the significant drop in $\Sigma_{1\mathrm{kpc, H}\alpha}$ from $z\sim1$ to $z\sim0.5$.

TeV halos are regions of enhanced photon emissivity surrounding pulsars. While multiple sources have been discovered, a self-consistent explanation of their radial profile and spherically-symmetric morphology remains elusive due to the difficulty in confining high-energy electrons and positrons within ~20 pc regions of the interstellar medium. One proposed solution utilizes anisotropic diffusion to confine the electron population within a "tube" that is auspiciously oriented along the line of sight. In this work, we show that while such models may explain a unique source such as Geminga, the phase space of such solutions is very small and they are unable to simultaneously explain the size and approximate radial symmetry of the TeV halo population.

The population of binary black hole mergers identified through gravitational waves has uncovered unexpected features in the intrinsic properties of black holes in the Universe. One particularly surprising and exciting result is the possible existence of black holes in the pair-instability mass gap, $\sim50-120~M_\odot$. Dense stellar environments can populate this region of mass space through hierarchical mergers, with the retention efficiency of black hole merger products strongly dependent on the escape velocity of the host environment. We use simple toy models to represent hierarchical merger scenarios in various dynamical environments. We find that hierarchical mergers in environments with high escape velocities ($\gtrsim300$ km/s) are efficiently retained, leading to an abundance of high-mass mergers that is potentially incompatible with the empirical mass spectrum from the current catalog of binary black hole mergers. Models which efficiently generate hierarchical mergers must therefore be tuned to avoid over-producing binary black hole mergers within and above the pair-instability mass gap. If hierarchical formation is indeed an important feature of dynamical formation channels with high escape velocities, it must be inhibited by some physical mechanism to avoid a cluster catastrophe that produces too high a rate of enormous stellar-mass black hole mergers.

Reservoirs of dense atomic gas (primarily hydrogen), contain approximately 90 percent of the neutral gas at a redshift of 3, and contribute to 2-3 percent of the total baryons in the Universe. These damped Lyman-${\alpha}$ systems (so called because they absorb Lyman-${\alpha}$ photons from within and from background sources) have been studied for decades, but only through absorption lines present in the spectra of background quasars and gamma-ray bursts. Such pencil beams do not constrain the physical extent of the systems. Here, we report integral-field spectroscopy of a bright, gravitationally lensed galaxy at a redshift of 2.7 with two foreground damped Lyman-${\alpha}$ systems. These systems are $>$ 238 $kpc^2$ in extent, with column densities of neutral hydrogen varying by more than an order of magnitude on $<$ 3 kpc-scales. The mean column densities are $10^{20.46}$ - $10^{20.84} cm^{-2}$ and the total masses are $> 5.5 \times 10^{8}$ - $1.4 \times 10^{9} M_{\odot}$, showing that they contain the necessary fuel for the next generation of star formation, consistent with relatively massive, low-luminosity primeval galaxies at redshifts $>$ 2.

White dwarf (WD) stars are often associated with the central stars of planetary nebulae (CSPNe) on their way to the cooling track. A large number of WD star candidates have been identified thanks to optical large-scale surveys such as Gaia DR2 and EDR3. However, hot-WD/CSPNe stars are quite elusive in optical bands due to their high temperatures and low optical luminosities. The Galaxy Evolution Explorer (GALEX) matched with the INT Galactic Plane Survey (IGAPS) allowed us to identify hot-WD candidates by combining the GALEX far-UV (FUV) and near-UV (NUV) with optical photometric bands (g, r, i and H$\alpha$). After accounting for source confusion and filtering bad photometric data, a total of 236485 sources were found in the GALEX and IGAPS footprint (GaPHAS). A preliminary selection of hot stellar sources was made using the GALEX colour cut on FUV-NUV>-0.53, yielding 74 hot-WD candidates. We analysed their spectral energy distribution (SED) by developing a fitting program for single- and two-body SED using an MCMC algorithm; 41 are probably binary systems (a binary fraction of ~55% was estimated). Additionally, we classified the WD star candidates using different infrared (IR) colours available for our sample obtaining similar results as in the SED analysis for the single and binary systems. This supports the strength of the fitting method and the advantages of the combination of GALEX UV with optical photometry. Ground-based time-series photometry and spectra are required in order to confirm the nature of the WD star candidates.

55 Cnc e is in a 0.73 day orbit transiting a Sun-like star. It has been observed that the occultation depth of this Super-Earth, with a mass of 8$M_{\bigoplus}$ and radius of 2$R_{\bigoplus}$, changes significantly over time at mid-infrared wavelengths. Observations with Spitzer measured a change in its day-side brightness temperature of 1200 K, possibly driven by volcanic activity, magnetic star-planet interaction, or the presence of a circumstellar torus of dust. Previous evidence for the variability in occultation was in the infrared range. Here we aim to explore if the variability exists also in the optical. TESS observed 55 Cnc during sectors 21, 44 and 46. We carefully detrend the data and fit a transit and occultation model for each sector in a Markov Chain Monte Carlo routine. In a later stage we use the Leave-One-Out Cross-Validation statistic to compare with a model of constant occultation for the complete set and a model with no occultation. We report an occultation depth of 8$\pm$2.5 ppm for the complete set of TESS observations. In particular, we measured a depth of 15$\pm$4 ppm for sector 21, while for sector 44 we detect no occultation. In sector 46 we measure a weak occultation of 8$\pm$5 ppm. The occultation depth varies from one sector to the next between 1.6 and 3.4 $\sigma$ significance. We derive the possible contribution on reflected light and thermal emission, setting an upper limit on the geometric albedo. Based on our model comparison the presence of an occultation is favoured considerably over no occultation, where the model with varying occultation across sectors takes most of the statistical weight. Our analysis confirms a detection of the occultation in TESS. Moreover, our results weakly lean towards a varying occultation depth between each sector, while the transit depth is constant across visits.

An energetic $\rm \gamma$-ray burst (GRB), GRB 211211A, was observed on 2021 December 11 by the Neil Gehrels Swift Observatory. Despite its long duration, typically associated with bursts produced by the collapse of massive stars, the discovery of an optical-infrared kilonova and a quasi-periodic oscillation during a gamma-ray precursor points to a compact object binary merger origin. The complete understanding of this nearby ($\sim$ 1 billion light-years) burst will significantly impact our knowledge of GRB progenitors and the physical processes that lead to electromagnetic emission in compact binary mergers. Here, we report the discovery of a significant ($\rm >5 \sigma$) transient-like emission in the high-energy $\rm \gamma$-rays (HE; E$>0.1$ GeV) observed by Fermi/LAT starting at $10^3$ s after the burst. After an initial phase with a roughly constant flux ($\rm \sim 5\times 10^{-10}\ erg\ s^{-1}\ cm^{-2}$) lasting $\sim 2\times 10^4$ s, the flux started decreasing and soon went undetected. The multi-wavelength afterglow emission observed at such late times is usually in good agreement with synchrotron emission from a relativistic shock wave that arises as the GRB jet decelerates in the interstellar medium. However, our detailed modelling of a rich dataset comprising public and dedicated multi-wavelength observations demonstrates that GeV emission from GRB 211211A is in excess with respect to the expectation of this scenario. We explore the possibility that the GeV excess is inverse Compton emission due to the interaction of a long-lived, low-power jet with an external source of photons. We discover that the kilonova emission can provide the necessary seed photons for GeV emission in binary neutron star mergers.

Recent years have seen unprecedentedly fast-growing prosperity in the commercial space industry. Several privately funded aerospace manufacturers, such as Space Exploration Technologies Corporation (SpaceX) and Blue Origin have innovated what we used to know about this capital-intense industry and gradually reshaped the future of human civilization. As private spaceflight and multi-planetary immigration gradually become realities from science fiction (sci-fi) and theory, both opportunities and challenges are presented. In this article, a review of the progress in space exploration and the underlying space technologies is firstly provided. For the next, a revisit and a prediction are paid and made to the K-Pg extinction event, the Chelyabinsk event, extra-terrestrialization, terraforming, planetary defense, including the emerging near-Earth object (NEO) observation and NEO impact avoidance technologies and strategies. Furthermore, a framework of the Solar Communication and Defense Networks (SCADN) with advanced algorithms and high efficacy is proposed to enable an internet of distributed deep-space sensing, communications, and defense to cope with disastrous incidents such as asteroid/comet impacts. Furthermore, the perspectives on the legislation, management, and supervision of founding the proposed SCADN are also discussed in depth.

The spin properties of merging black holes observed with gravitational waves can offer novel information about the origin of these systems. The magnitude and orientations of black hole spins offer a record of binaries' evolutionary history, encoding information about massive stellar evolution and the astrophysical environments in which binary black holes are assembled. Recent analyses of the binary black hole population have yielded conflicting portraits of the black hole spin distribution. Some work suggests that black hole spins are small but non-zero and exhibit a wide range of misalignment angles relative to binaries' orbital angular momenta. Other work concludes that the majority of black holes are non-spinning while the remainder are rapidly rotating and primarily aligned with their orbits. We revisit these conflicting conclusions, employing a variety of complementary methods to measure the distribution of spin magnitudes and orientations among binary black hole mergers. We find that the existence of a sub-population of black hole with vanishing spins is not required by current data. Should such a sub-population exist, we find that it must contain $\lesssim 60\%$ of binaries. Additionally, we find evidence for significant spin-orbit misalignment among the binary black hole population, with some systems exhibiting misalignment angles greater than $90^{\circ}$, and see no evidence for an approximately spin-aligned sub-population.

Particle acceleration during magnetic reconnection is a long-standing topic in space, solar and astrophysical plasmas. Recent 3D particle-in-cell simulations of magnetic reconnection show that particles can leave flux ropes due to 3D field-line chaos, allowing particles to access additional acceleration sites, gain more energy through Fermi acceleration, and develop a power-law energy distribution. This 3D effect does not exist in traditional 2D simulations, where particles are artificially confined to magnetic islands due to their restricted motions across field lines. Full 3D simulations, however, are prohibitively expensive for most studies. Here, we attempt to reproduce 3D results in 2D simulations by introducing ad hoc pitch-angle scattering to a small fraction of the particles. We show that scattered particles are able to transport out of 2D islands and achieve more efficient Fermi acceleration, leading to a significant increase of energetic particle flux. We also study how the scattering frequency influences the nonthermal particle spectra. This study helps achieve a complete picture of particle acceleration in magnetic reconnection.

The dynamo driven by the magnetorotational instability (MRI) is believed to play an important role in the dynamics of accretion discs and may also explain the origin of the extreme magnetic fields present in magnetars. Its saturation level is an important open question known to be particularly sensitive to the diffusive processes through the magnetic Prandtl number Pm (the ratio of viscosity to resistivity). Despite its relevance to proto-neutron stars and neutron star merger remnants, the numerically challenging regime of high Pm is still largely unknown. Using zero-net flux shearing box simulations in the incompressible approximation, we studied MRI-driven dynamos at unprecedentedly high values of Pm reaching 256. The simulations show that the stress and turbulent energies are proportional to Pm up to moderately high values ($\mathrm{Pm} \sim 50$). At higher Pm, they transition to a new regime consistent with a plateau independent of Pm for $\rm Pm \gtrsim 100$. This trend is independent of the Reynolds number, which may suggest an asymptotic regime where the energy injection and dissipation are independent of the diffusive processes. Interestingly, large values of Pm not only lead to intense small-scale magnetic fields but also to a more efficient dynamo at the largest scales of the box.

Understanding magnetic activity on the surface of stars other than the Sun is important for exoplanet analyses to properly characterize an exoplanet's atmosphere and to further characterize stellar activity on a wide range of stars. Modeling stellar surface features of a variety of spectral types and rotation rates are key to understanding of the magnetic activity of these stars. Using data from Kepler, we use the starspot modeling program STarSPot (STSP) to measure the position and size of spots for KOI-340 which is an eclipsing binary consisting of a subgiant star ($\rm T_{\rm eff} = 5593 \pm 27 K; R_{*} = 1.98 \pm 0.05 R_{\odot}$) with an M-dwarf companion ($M_{*} = 0.214 \pm 0.006 M_{\odot}$). STSP uses a novel technique to measure the spot positions and radii by using the transiting secondary to study and model individual active regions on the stellar surface using high-precision photometry. We find the average size of spot features on KOI-340's primary is $\sim$10% the radius of the star, i.e. two times larger than the mean size of Solar-maximum sunspots. The spots on KOI-340 are present at every longitude and show possible signs of differential rotation. The minimum fractional spotted area of KOI-340's primary is $2\substack{+12\\ -2} \%$ while the spotted area of the Sun is at most 0.2%. One transit of KOI-340 shows a signal in the transit consistent with a plage; this plage occurs right before a dark spot indicating the plage and spot might be co-located on the surface of the star.

The accurate localization of gamma-ray bursts remains a crucial task. While historically, improved localization have led to the discovery of afterglow emission and the realization of their cosmological distribution via redshift measurements, a more recent requirement comes with the potential of studying the kilonovae of neutron star mergers. Gravitational wave detectors are expected to provide locations to not better than 10 square degrees over the next decade. With their increasing horizon for merger detections also the intensity of the gamma-ray and kilonova emission drops, making their identification in large error boxes a challenge. Thus, a localization via the gamma-ray emission seems to be the best chance to mitigate this problem. Here we propose to equip some of the second generation Galileo satellites with dedicated GRB detectors. This saves costs for launches and satellites for a dedicated GRB network, the large orbital radius is beneficial for triangulation, and perfect positional and timing accuracy come for free. We present simulations of the triangulation accuracy, demonstrating that short GRBs as faint as GRB 170817A can be localized to 1 degree radius (1 sigma).

Recent gamma-ray and radio studies have obtained some stringent constraints on annihilating dark matter properties. However, only a few studies have focussed on using X-ray data to constrain annihilating dark matter. In this article, we perform the X-ray analysis of annihilating dark matter using the data of the Omega Centauri cluster. If dark matter is the correct interpretation of the non-luminous mass component derived in the Omega Centauri cluster, the conservative lower limits of thermal dark matter mass annihilating via the $\tau^+\tau^-$, $b\bar{b}$ and $W^+W^-$ channels can be significantly improved to 104(43) GeV, 650(167) GeV and 480(137) GeV respectively, assuming the diffusion coefficient $D_0 \le 10^{26}(10^{27})$ cm$^2$/s. These constraints can safely rule out the recent claims of dark matter interpretation of the gamma-ray excess and anti-proton excess seen in our Galaxy. Generally speaking, the conservative lower limits obtained for non-leptophilic annihilation channels are much more stringent than that obtained by gamma-ray analysis of nearby dwarf galaxies. We anticipate that this would open a new window for constraining annihilating dark matter.

Triple stars and compact objects are ubiquitously observed in nature. Their long-term evolution is complex; in particular, the von-Zeipel-Lidov-Kozai (ZLK) mechanism can potentially lead to highly eccentric encounters of the inner binary. Such encounters can lead to a plethora of interacting binary phenomena, as well as stellar and compact-object mergers. Here we find explicit analytical formulae for the maximal eccentricity, $e_{\rm max}$, of the inner binary undergoing ZLK oscillations, where both the test particle limit (parametrised by the inner-to-outer angular momentum ratio $\eta$) and the double-averaging approximation (parametrised by the period ratio, $\epsilon_{\rm SA}$) are relaxed, for circular outer orbits. We recover known results in both limiting cases (either $\eta$ or $\epsilon_{\rm SA} \to 0$) and verify the validity of our model using numerical simulations. We test our results with two accurate numerical N-body codes, $\texttt{Rebound}$ for Newtonian dynamics and $\texttt{Tsunami}$ for general-relativistic (GR) dynamics, and find excellent correspondence. We discuss the implications of our results for stellar triples and both stellar and supermassive triple black hole mergers.

We report on simultaneous optical and X-ray observations of the dwarf nova SS Cyg with Tomo-e Gozen/1.05 m Kiso Schmidt and Neutron star Interior Composition ExploreR (it NICER) / International Space Station (ISS). A total of four observations were carried out in the quiescent state and highly correlated light variations between the two wavelengths were detected. We have extracted local brightness peaks in the light curves with a binning interval of 1 sec, called `shots', and have evaluated time lags between the optical and X-ray variations by using a cross-correlation function. Some shots exhibit significant optical lags to X-ray variations and most of them are positive ranging from $+$0.26 to 3.11 sec, which have never been detected. They may be ascribable to X-ray reprocessing in the accretion disk and/or the secondary star. This analysis is possible thanks to the high timing accuracy and the high throughput of NICER and the matching capability of Tomo-e Gozen. Also, it is confirmed that the correlation between the optical and X-ray variations changed in the middle of one of our observation windows and the X-ray spectrum softer than 2 keV varied in accordance with the correlation.

Radio frequency interference (RFI) is a significant challenge faced by today's radio astronomers. While most past efforts were devoted to cleaning the RFI from the data, we develop a novel method for categorizing and cataloguing RFI for forensic purpose. We present a classifier that categorizes RFI into different types based on features extracted using Principal Component Analysis (PCA) and Fourier analysis. The classifier can identify narrowband non-periodic RFI above 2 sigma, narrowband periodic RFI above 3 sigma, and wideband impulsive RFI above 5 sigma with F1 scores between 0.87 and 0.91 in simulation. This classifier could be used to identify the sources of RFI as well as to clean RFI contamination (particularly in pulsar search). In the long-term analysis of the categorized RFI, we found a special type of drifting periodic RFI that is detrimental to pulsar search. We also found evidences of an increased rate of impulsive RFI when the telescope is pointing toward the cities. These results demonstrate this classifier's potential as a forensic tool for RFI environment monitoring of radio telescopes.

We explore the formation, energetics and geometry of the jets and variability of their central engine by means of 3D general relativistic magneto-hydrodynamical simulations of the magnetically arrested disk accretion in Kerr geometry. We study both fast and slowly rotating black holes, and address our simulations to both active galaxy centers and Gamma Ray Burst engines. The structured jets are postulated to account for emission properties of high energy sources across the mass scale, launched from stellar mass black holes in GRBs and from supermassive black holes in AGN. Their active cores contain magnetized accretion disks and rotation of the Kerr black hole provides mechanism for launching relativistic jets. This process works most effectively if the mode of accretion turns out to be magnetically arrested. In this mode, the modulation of jets launched from the engine is related to internal instabilities in the accretion flow, that work on smallest time and spatial scales. As these scales are related to the light crossing time, and the black hole gravitational radius, the universal model of jet-disk connection should scale with the mass of the black hole. The system evolution is governed by the physical parameters of the engine, such as the black hole spin, and disk size, as well as disk magnetisation. We find that stronger magnetic fields may lead to jet quenching. The effect is supposed to be important mainly for Gamma Ray Burst jets and may be related to the magnetically driven winds from their engines.

Spiral galaxies can be classified into the Grand-designs and Flocculents, based on the nature of their spiral arms. The Grand-designs exhibit almost continuous, high contrast spiral arms and are believed to be driven by density waves; the Flocculents, on the other hand, have patchy or discontinuous spiral features and are mostly stochastic in origin. We construct a convolutional neural network (CNN) model that classifies spirals into Grand-designs and Flocculents, with a testing accuracy of $\mathrm{97.2\%}$. We then use the above model for classifying $\mathrm{1,220}$ new spirals from the SDSS. Out of these, $\mathrm{721}$ are identified as Flocculents, the rest being Grand-designs. The mean asymptotic rotational velocity of our sample Grand-designs and Flocculents are $\mathrm{218 \; km \; s^{-1}}$ and $\mathrm{145 \; km \; s^{-1}}$ respectively, while their respective de Vaucouleur morphological type indices are $\mathrm{2.6}$ and $\mathrm{4.7}$. This possibly indicates that Grand-designs are mostly ordinary high surface brightness galaxies like our Milky Way, while Flocculents are intermediate-mass low surface brightness galaxies.

We report on a very bright, long-duration gamma-ray burst (GRB), GRB~220426A, observed by \emph{Fermi} satellite. GRB~220426A with total duration of $T_{90}=6$~s is composed with two main pulses and some sub-peaks. The spectral analysis of this burst with Band function reveals that both the time-integrated and the time-resolved spectra are very narrow with high $\alpha \gtrsim 0.2$ and low $\beta\lesssim -3.1$. It is strong reminiscent of GRB~090902B, a special GRB with identification of the photospheric emission. Then, we perform the spectral analysis of this burst based on a non-dissipated photospheric emission, which can be well modelled as the multicolor-blackbody with a cutoff power-law distribution of the thermal temperature. The spectral fittings reveal that the photospheric emission can well describe the radiation spectrum of this burst. We conclude that this burst would be a second burst in the class of GRB~090902B observed by \emph{Fermi} satellite. We also discuss the physics of photosphere and the origin of the high-energy component in GRB~220426A .

Context: Class II methanol masers at 6.7 GHz originate close to high-mass young stellar objects (HMYSOs). Their high sensitivity to local condition variations makes them a useful marker of the activity of the emerging massive stars. Aims: We aim to closely examine the variability of the 6.7 GHz methanol masers in Cep A HW2 using the new and archival single-dish and high-resolution Very-Long-Baseline Interferometry (VLBI) datasets. Methods: We monitored 6.7 GHz methanol masers towards the target between 2009 and 2021 using the Torun 32 m radio telescope, and analysed nine datasets of VLBI observations taken over 16 yr. Results: Faint, extremely redshifted maser emission located close to the presumed central star position and coincident with the radio jet shows flaring activity with a period of ~ 5 yr. Flares have an asymmetric profile with a rise-to-decay time ratio of 0.18 and relative amplitude higher than 10. The velocity of the flaring cloudlets drifts at a rate of 3-4 * 10^-5 km/s/d for about 11.5 yr of the monitoring. The time-lag between the peaks of the two flaring features implies a propagation speed of the exciting factor of ~ 1000 km/s. Synchronised and anticorrelated variations of the flux density of blue- and redshifted features begin ~ 2.5 yr after flares of the extremely redshifted emission.

We present observations of a solar plage in the millimeter-continuum with the ALMA and in the Ca 8542 and Na 5896 spectral lines with the Interferometric BIdimensional Spectrometer (IBIS). Our goal is to compare the measurement of local gas temperatures provided by ALMA with the temperature diagnostics provided by non-LTE inversions using STIC. In performing these inversions, we find that using column mass as the reference height scale, rather than optical depth, provides more reliable atmospheric profiles above the temperature minimum and that the treatment of non- LTE hydrogen ionization brings the inferred chromospheric temperatures into better agreement with the ALMA measurements. The Band 3 brightness temperatures are higher but well correlated with the inversion-derived temperatures at the height of formation of the Ca 8542 line core. The Band 6 temperatures instead do not show good correlations with the temperatures at any specific layer in the inverted atmospheres. We then performed inversions that included the millimeter continuum intensities as an additional constraint. Incorporating Band 3 generally resulted in atmospheres showing a strong temperature rise in the upper atmosphere, while including Band 6 led to significant regions of anomalously low temperatures at chromospheric heights. This is consistent with the idea that the Band 6 emission can come from a range of heights. The poor constraints on the chromospheric electron density with existing inversion codes introduces difficulties in determining the height(s) of formation of the millimeter continuum as well as uncertainties in the temperatures derived from the spectral lines.

Young multiple systems accrete most of their final mass in the first few Myr of their lifetime, during the protostellar and protoplanetary phases. Previous studies showed that in binary systems the majority of the accreted mass falls onto the lighter star, thus evolving to mass equalisation. However, young stellar systems often comprise more than two stars, which are expected to be in hierarchical configurations. Despite its astrophysical relevance, differential accretion in hierarchical systems remains to be understood. In this work, we investigate whether the accretion trends expected in binaries are valid for higher order multiples. We performed a set of 3D Smoothed Particle Hydrodynamics simulations of binaries and of hierarchical triples (HTs) embedded in an accretion disc, with the code Phantom. We identify for the first time accretion trends in HTs and their deviations compared to binaries. These deviations, due to the interaction of the small binary with the infalling material from the circum-triple disc, can be described with a semi-analytical prescription. Generally, the smaller binary of a HT accretes more mass than a single star of the same mass as the smaller binary. We found that in a HT, if the small binary is heavier than the third body, the standard differential accretion scenario (whereby the secondary accretes more of the mass) is hampered. Reciprocally, if the small binary is lighter than the third body, the standard differential accretion scenario is enhanced. The peculiar differential accretion mechanism we find in HTs is expected to affect their mass ratio distribution.

We use sparse regression methods (SRM) to build accurate and explainable models that predict the stellar mass of central and satellite galaxies as a function of properties of their host dark matter halos. SRM are machine learning algorithms that provide a framework for modelling the governing equations of a system from data. In contrast with other machine learning algorithms, the solutions of SRM methods are simple and depend on a relatively small set of adjustable parameters. We collect data from 35,459 galaxies from the EAGLE simulation using 19 redshift slices between $z=0$ and $z=4$ to parameterize the mass evolution of the host halos. Using an appropriate formulation of input parameters, our methodology can model satellite and central halos using a single predictive model that achieves the same accuracy as when predicted separately. This allows us to remove the somewhat arbitrary distinction between those two galaxy types and model them based only on their halo growth history. Our models can accurately reproduce the total galaxy stellar mass function and the stellar mass-dependent galaxy correlation functions ($\xi(r)$) of EAGLE. We show that our SRM model predictions of $\xi(r)$ is competitive with those from sub-halo abundance matching and performs better than results from extremely randomized trees. We suggest SRM as an encouraging approach for populating the halos of dark matter only simulations with galaxies and for generating mock catalogues that can be used to explore galaxy evolution or analyse forthcoming large-scale structure surveys.

KIC 5359678 is a 6.231-d period F-type eclipsing binary system whose component stars both show starspot activity. It was observed by the Kepler satellite in long cadence for the full four-year duration of the mission. Wang et al (2021) obtained radial velocity measurements of the two stars and analysed these plus the Kepler data to study their spot activity and measure their physical properties, but left several questions unanswered. We have performed an independent analysis and determined the masses (1.252 +/- 0.018 and 1.065 +/- 0.013 Msun) and radii (1.449 +/- 0.012 and 1.048 +/- 0.017 Rsun) of the stars to high precision. The distance we find to the system is slightly shorter than that from Gaia EDR3 for unknown reason(s). We also investigated the precision of the numerical integration applied to the model light curve to match the 1765-s sampling cadence of the Kepler observations. We found that ignoring this temporal smearing leads to biased radius measurements for the stars: that for the primary is too small by 4 sigma and that for the secondary is too large by 10 sigma. Doubling the sampling rate of the model light curve is sufficient to remove most of this bias, but for precise results a minimum of five samples per observed datapoint is required.

V1388 Ori is an early-B type detached eclipsing binary whose physical properties have previously been measured from dedicated spectroscopy and a ground-based survey light curve. We reconsider the properties of the system using newly-available light curves from the Transiting Exoplanet Survey Satellite (TESS). We discover two frequencies in the system, at 2.99 d$^{-1}$ and 4.00 d$^{-1}$ which are probably due to beta Cephei or slowly-pulsating B-star pulsations. A large number of additional significant frequencies exist at multiples of the orbital frequency, 0.4572 d$^{-1}$. We are not able to find a fully satisfactory model of the eclipses, but the best attempts show highly consistent values for the fitted parameters. We find masses of 7.24 +/- 0.08 Msun and 5.03 +/- 0.04 Msun, and radii of 5.30 +/- 0.07 Rsun and 3.14 +/- 0.06 Rsun. The properties of the system are in good agreement with the predictions of theoretical stellar evolutionary models and the Gaia EDR3 parallax if the published temperature estimates are revised downwards by 1500 K, to 19000 K for the larger and more massive star and 17000 K for its companion.

Thermal conductivity is one of the important mechanisms of heat transfer in the solar corona. In the limit of strongly magnetized plasma, it is typically modeled by Spitzer's expression where the heat flux is aligned with the magnetic field. This paper describes the implementation of the heat conduction into the code MANCHA3D with an aim of extending single-fluid MHD simulations from the upper convection zone into the solar corona. Two different schemes to model heat conduction are implemented: (1) a standard scheme where a parabolic term is added to the energy equation, and (2) a scheme where the hyperbolic heat flux equation is solved. The first scheme limits the time step due to the explicit integration of a parabolic term, which makes the simulations computationally expensive. The second scheme solves the limitations on the time step by artificially limiting the heat conduction speed to computationally manageable values. The validation of both schemes is carried out with standard tests in one, two, and three spatial dimensions. Furthermore, we implement the model for heat flux derived by Braginskii (1965) in its most general form, when the expression for the heat flux depends on the ratio of the collisional to cyclotron frequencies of the plasma, and, therefore on the magnetic field strength. Additionally, our implementation takes into account the heat conduction in parallel, perpendicular, and transverse directions, and provides the contributions from ions and electrons separately. The model also transitions smoothly between field-aligned conductivity and isotropic conductivity for regions with a low or null magnetic field. Finally, we present a two-dimensional test for heat conduction using realistic values of the solar atmosphere where we prove the robustness of the two schemes implemented.

Through several articles we have developed a model for dark matter as consisting of bubbles of a new (speculated) type of vacuum, starting from cm-sized pearls or balls down to atomic size ones and now we believe they have nanometer sizes. In the latest development of our model we have the bubbles of the new vacuum imbedded in dust grains very similar to the grains present in interstellar and intergalactic space anyway, although the presence of the bubble with a very large homolumo gap in its single electron spectrum influences the dust grain material so as to become denser and harder. We have earlier explained how our dark matter particles get stopped in the shielding, so that normally expected nucleonic collisions are not observable. The signal of the dark matter in the underground experiments rather becomes decays of excited particles actually with the energy of the homolumo gap, which is also equal to the photon energy of the X-ray line presumably observed astronomically from galaxy clusters etc. A new calculation here is a fitting of the velocity dependence of the dark matter self-interaction as estimated by Correa [15], using deviations from the only gravitationally interacting dark matter in dwarf galaxies. Let us stress that apart from the speculated new vacuum we have no new physics, and if the couplings in the Standard Model were adjusted to make degenerate vacua as speculated according to our Multiple Point Principle (MPP) we would only need the Standard Model, so dark matter would not require new physics.

Displacement sensors have a variety of applications within gravitational wave detectors. The seismic isolation chain of the LIGO core optics utilises optical shadow sensors for their stabilisation. Future upgrades, such as LIGO Voyager, plan to operate at cryogenic temperatures to reduce their thermal noise and will require cryogenic displacement sensors. We present the results of simulations and experimental tests of the shadow sensors embedded in the Birmingham Optical Sensors and Electromagnetic Motors (BOSEMs). We determine that the devices can reliably viability operate at 100 K. We also show that the performance of the BOSEM sensors improves at cryogenic temperatures.

Active galactic nuclei (AGN) jets play an important role as a feedback mechanism that quenches the growth of massive galaxies. These jets have so far been simulated almost exclusively with grid-based codes. In this work we present results from hydrodynamical tests of AGN jets simulated with the SWIFT code and its implementation of smoothed particle hydrodynamics (SPH). We reach numerical resolutions of more than a million particles per jet, unprecedented for an SPH code. In all cases we find broad agreement between our jets and theoretical predictions for the jet and lobe lengths and widths. At very low resolutions, typical of low-resolution cosmological simulations ($m_\mathrm{gas}\approx10^7$ $\mathrm{M}_\odot$), the shape of the lobes is close to self-similar ones at a level of $15\%$ accuracy. This indicates that the basics of jet lobe physics can be captured even at such resolutions ($\approx500$ particles per jet). The jets first evolve ballistically, and then transition to a self-similar phase, during which the kinetic and thermal energies in the lobes and the shocked ambient medium are constant fractions of the total injected energy. In our standard simulation, two thirds of the initially injected energy is transferred to the ambient medium by the point the jets are turned off, mainly through a bow shock. Of that, three quarters is in thermal form, indicating that the bow shock thermalizes efficiently. We find that jet launching velocity and power, in addition to determining the temperature, mass and overall length of the jets, also play an important role in the resolution of the jets.

We aim to explore spectral signatures of the predicted multi-ion UFOs in the broadband X-ray spectra of active galactic nuclei (AGNs) by exploiting an accretion disk wind model in the context of a simple magnetohydrodynamic (MHD) framework. We are focused primarily on examining the spectral dependences on a number of key properties; (1) ionizing luminosity ratio $\lambda_{\rm ion}$, (2) line-of-sight wind density slope $p$, (3) optical/UV-to-X-ray strength $\alpha_{\rm OX}$, (4) inclination $\theta$, (5) X-ray photon index $\Gamma$ and (6) wind density factor $f_D$. With an emphasis on radio-quiet Seyferts in sub-Eddington regime, multi-ion UFO spectra are systematically calculated as a function of these parameters to show that MHD-driven UFOs imprint a unique asymmetric absorption line profile with a pronounced blue tail structure on average. Such a characteristic line signature is generic to the simplified MHD disk-wind models presented in this work due to their specific kinematics and density structure. The properties of these absorption line profiles could be utilized as a diagnostics to distinguish between different wind driving mechanisms or even the specific values of a given MHD wind parameters. We also present high fidelity microcalorimeter simulations in anticipation of the upcoming {\it XRISM}/Resolve and {\it Athena}/X-IFU instruments to demonstrate that such a "tell-tale" sign may be immune to a spectral contamination by the presence of additional warm absorbers and partially covering gas.

We report on the first detection of linearly polarized x-ray emission from an ultra-magnetized neutron star with the Imaging X-ray Polarimetry Explorer (IXPE). The IXPE 35 observations of the anomalous x-ray pulsar 4U 0142+61 reveal a linear polarization degree of $(12\pm 1)\%$ throughout the IXPE 2--8 keV band. We detect a substantial variation of the polarization with energy: the degree is $(14\pm 1)\%$ at 2--4 keV and $(41\pm 7)\%$ at 5.5--8 keV, while it drops below the instrumental sensitivity around 4--5 keV, where the polarization angle swings by $\sim 90^\circ$. The IXPE observations give us completely new information about the properties of the neutron star surface and magnetosphere and lend further support to the presence of the quantum mechanical effect of vacuum birefringence.

We simulate the space environment around AU Microscopii b and the interaction between the magnetized stellar wind with a planetary atmospheric outflow for ambient stellar wind conditions and Coronal Mass Ejection (CME) conditions. We also calculate synthetic Ly$\alpha$ absorption due to neutral hydrogen in the ambient and the escaping planetary atmosphere affected by this interaction. We find that the Ly$\alpha$ absorption is highly variable due to the highly-varying stellar wind conditions. A strong Doppler blue-shift component is observed in the Ly$\alpha$ profile, in contradiction to the actual escape velocity observed in the simulations themselves. This result suggest that the strong Doppler blue-shift is likely attributed to the stellar wind, not the escaping neutral atmosphere, either through its advection of neutral planetary gas, or through the creation of a fast neutral flow via charge exchange between the stellar wind ions and the planetary neutrals. Indeed, our CME simulations indicate a strong stripping of magnetospheric material from the planet, including some of the neutral escaping atmosphere. Our simulations show that the pressure around close-in exoplanets is not much lower, and may be even higher, than the pressure at the top of the planetary atmosphere. Thus, the neutral atmosphere is hydrodynamically escaping with a very small velocity ($<15~km~s^{-1}$). Moreover, our simulations show that an MHD treatment is essential in order to properly capture the coupled magnetized stellar wind and the escaping atmosphere, despite of the atmosphere being neutral. This coupling should be considered when interpreting Ly$\alpha$observations in the context of exoplanets atmospheric escape.

Precise Gaia measurements of positions, parallaxes, and proper motions provide an opportunity to calculate 3D positions and 2D velocities (i.e. 5D phase-space) of Milky Way stars. Where available, spectroscopic radial velocity (RV) measurements provide full 6D phase-space information, however there are now and will remain many stars without RV measurements. Without an RV it is not possible to directly calculate 3D stellar velocities, however one can infer 3D stellar velocities by marginalizing over the missing RV dimension. In this paper, we infer the 3D velocities of stars in the Kepler field in Cartesian Galactocentric coordinates (vx, vy, vz). We directly calculate velocities for around a quarter of all Kepler targets, using RV measurements available from the Gaia, LAMOST and APOGEE spectroscopic surveys. Using the velocity distributions of these stars as our prior, we infer velocities for the remaining three-quarters of the sample by marginalizing over the RV dimension. The median uncertainties on our inferred vx, vy, and vz velocities are around 4, 18, and 4 km/s, respectively. We provide 3D velocities for a total of 148,590 stars in the Kepler field. These 3D velocities could enable kinematic age-dating, Milky Way stellar population studies, and other scientific studies using the benchmark sample of well-studied Kepler stars. Although the methodology used here is broadly applicable to targets across the sky, our prior is specifically constructed from and for the Kepler field. Care should be taken to use a suitable prior when extending this method to other parts of the Galaxy.

We present [CII] synthetic observations of smoothed particle hydrodynamics (SPH) simulations of a dwarf galaxy merger. The merging process varies the star-formation rate by more than three orders of magnitude. Several star clusters are formed, the feedback of which disperses and unbinds the dense gas through expanding HII regions and supernova (SN) explosions. For galaxies with properties similar to the modelled ones, we find that the [CII] emission remains optically thin throughout the merging process. We identify the Warm Neutral Medium ($3<\log T_{\rm gas}<4$ with $\chi_{\rm HI}>2\chi_{\rm H2}$) to be the primary source of [CII] emission ($\sim58\%$ contribution), although at stages when the HII regions are young and dense (during star cluster formation or SNe in the form of ionized bubbles) they can contribute $\gtrsim50\%$ to the total [CII] emission. We find that the [CII]/FIR ratio decreases due to thermal saturation of the [CII] emission caused by strong FUV radiation fields emitted by the massive star clusters, leading to a [CII]-deficit medium. We investigate the [CII]-SFR relation and find an approximately linear correlation which agrees well with observations, particularly those from the Dwarf Galaxy Survey. Our simulation reproduces the observed trends of [CII]/FIR versus $\Sigma_{\rm SFR}$ and $\Sigma_{\rm FIR}$, and it agrees well with the Kennicutt relation of SFR-FIR luminosity. We propose that local peaks of [CII] in resolved observations may provide evidence for ongoing massive cluster formation.

We study the merger history of primordial black holes (PBHs) in a scenario where they represent the dominant dark matter component of a typical dwarf galaxies' core. We investigate the possibility of a sequence of collisions resulting in a hierarchical merger of black holes and look at the final mass spectrum in such {\it clusters}, which initially present a monochromatic (single-mass) PBH population. Our study shows that the merging process results in the transfer of about $40\%$ of the total mass of the core to the merger products regardless of the initial mass of PBHs, with about $5\%$ of energy radiated out in the form of gravitational waves. We find that, in the lighter mass limit, black holes up to eight times more massive than the original population can be formed within a Hubble time.

The search for primordial $B$-mode polarization of the CMB is limited by the sample variance of $B$-modes produced at later times by gravitational lensing. Constraints can be improved by `delensing': using some proxy of the matter distribution to partially remove the lensing-induced $B$-modes. Current and soon-upcoming experiments will infer a matter map -- at least in part -- from the temperature anisotropies of the CMB. These reconstructions are contaminated by extragalactic foregrounds: radio-emitting galaxies, the cosmic infrared background, or the Sunyaev--Zel'dovich effects. Using the Websky simulations, we show that the foregrounds add spurious power to the angular auto-spectrum of delensed $B$-modes via non-Gaussian higher-point functions, biasing constraints on the tensor-to-scalar ratio, $r$. We consider an idealized experiment similar to the Simons Observatory, with no Galactic or atmospheric foregrounds. After removing point sources detectable at 143 GHz and reconstructing lensing from CMB temperature modes $l<3500$ using a Hu-Okamoto quadratic estimator (QE), we infer a value of $r$ that is $1.5\,\sigma$ higher than the true $r=0$. Reconstructing instead from a minimum-variance ILC map only exacerbates the problem, bringing the bias above $3\,\sigma$. When the $TT$ estimator is co-added with other QEs or with external matter tracers, new couplings ensue which partially cancel the diluted bias from $TT$. We provide a simple and effective prescription to model these effects. In addition, we demonstrate that the point-source-hardened or shear-only QEs can not only mitigate the biases to acceptable levels, but also lead to lower power than the Hu-Okamoto QE after delensing. Thus, temperature-based reconstructions remain powerful tools in the quest to measure $r$.

A system of 5,020 robotic fiber positioners was installed in 2019 on the Mayall Telescope, at Kitt Peak National Observatory. The robots automatically re-target their optical fibers every 10 - 20 minutes, each to a precision of several microns, with a reconfiguration time less than 2 minutes. Over the next five years, they will enable the newly-constructed Dark Energy Spectroscopic Instrument (DESI) to measure the spectra of 35 million galaxies and quasars. DESI will produce the largest 3D map of the universe to date and measure the expansion history of the cosmos. In addition to the 5,020 robotic positioners and optical fibers, DESI's Focal Plane System includes 6 guide cameras, 4 wavefront cameras, 123 fiducial point sources, and a metrology camera mounted at the primary mirror. The system also includes associated structural, thermal, and electrical systems. In all, it contains over 675,000 individual parts. We discuss the design, construction, quality control, and integration of all these components. We include a summary of the key requirements, the review and acceptance process, on-sky validations of requirements, and lessons learned for future multi-object, fiber-fed spectrographs.

The structures in the Universe are distributed in the cosmic web. The structures' distribution, statistics, and evolution can be considered probes of cosmological models. We investigate the number density of voids and dark matter halo-in-voids in the Excursion Set Theory (EST). We study both the spherical and ellipsoidal collapse models in the Markov and non-Markov frameworks of the EST. Then, we compare the number density of voids and halo-in-voids in standard $\Lambda$CDM and reconstructed model. The reconstructed model is a model-independent expansion history reconstruction from background observations. In this work, the effect of the collapse model barrier in exclusion of the trajectories in EST is explored. This exclusion affects the statistics of voids and the halo-in-voids. We show that the statistic of the voids in the non-Markov extension of EST is almost independent of the collapse model. Finally, we show that the number density of halo-in-voids can distinguish these two models and address the anomalies such as void phenomena.

We report on the discovery and characterisation of three planets orbiting the F8 star HD~28109, which sits comfortably in \tess's continuous viewing zone. The two outer planets have periods of $\rm 56.0067 \pm 0.0003~days$ and $\rm 84.2597_{-0.0008}^{+0.0010}~days$, which implies a period ratio very close to that of the first-order 3:2 mean motion resonance, exciting transit timing variations (TTVs) of up to $\rm 60\,mins$. These two planets were first identified by \tess, and we identified a third planet in the \textcolor{black}{\tess photometry} with a period of $\rm 22.8911 \pm 0.0004~days$. We confirm the planetary nature of all three planetary candidates using ground-based photometry from Hazelwood, ASTEP and LCO, including a full detection of the $\rm \sim9\,h$ transit of HD~28109 c from Antarctica. The radii of the three planets are \textcolor{black}{$\rm R_b=2.199_{-0.10}^{+0.098} ~R_{\oplus}$, $\rm R_c=4.23\pm0.11~ R_{\oplus}$ and $\rm R_d=3.25\pm0.11 ~R_{\oplus}$}; we characterise their masses using TTVs and precise radial velocities from ESPRESSO and HARPS, and find them to be $\rm M_b=18.5_{-7.6}^{+9.1}~M_{\oplus}$, $\rm M_c=7.9_{-3.0}^{+4.2}~M_{\oplus}$ and $\rm M_d=5.7_{-2.1}^{+2.7}~M_{\oplus}$, making planet b a dense, massive planet while c and d are both under-dense. We also demonstrate that the two outer planets are ripe for atmospheric characterisation using transmission spectroscopy, especially given their position in the CVZ of JWST. The data obtained to date are consistent with resonant (librating) and non-resonant (circulating) solutions; additional observations will show whether the pair is actually locked in resonance or just near-resonant.

The observational high-resolution data on the $HI$ gas disk rotation curve of the spiral galaxy $NGC\;6946$ at large radii and the recent increased amount and quality of data on linearly polarized radio continuum emission for this galaxy provide us to re-investigate the role of the regular magnetic fields in the rotation of gas and to test whether disk magnetic fields should be considered as a non-negligible dynamical ingredient. The spiral galaxy $NGC\;6946$ hosts two symmetric bright magnetic spiral arms which have been revealed in the radio polarization map. By taking into account the dynamical effect of the regular magnetic fields caused by two main magnetic arms, on the circular gas rotation and considering two dark matter mass density models, ISO and the universal NFW profile, the shape of the $HI$ gas rotation curve of this galaxy is fitted better, especially in the outer part. The contribution of the regular magnetic fields in the rotation velocity of the gas has a positive value and shows an ascending curve with a typical amplitude of about $ 7 - 14 \; km s^{-1}$ in the outer gaseous disc of the galaxy $NGC\;6946$. We also generate the map of the modeled regular magnetic field strength for the spiral galaxy $NGC\;6946$, which clearly shows two main inner contours with the spiral pattern and the highest magnetic field strength, with the field directed towards the galaxy's center in both, almost compatible with two inner main magnetic arms detected in the polarized synchrotron intensity map.

3D maps of the extinction density in the Galaxy can be built through the inversion of catalogues of distance-extinction pairs for individual target stars. The spatial resolution of the maps that can be achieved increases with the spatial density of the targets, and subsequently with the combination of catalogues. However, this requires their careful inter-calibration. Our aim is to develop methods of inter-calibration of two different catalogues. We used as reference a spectrophotometric catalogue. A principal component analysis was performed in G,GB,GR,J,H,K multi-colour space for the second catalog. The subspace constituted by the two first components was split into cells in which we estimated deviations from the reference. The deviations were computed using all targets from the reference located at a short spatial distance of each secondary target. Corrections and filtering were deduced for each cell in the multi-colour space. We applied the technique to two different datasets: on the one hand, the spectrophotometric catalogue, and, on the other hand, a catalogue of extinctions based on photometry of Gaia eDR3 and 2MASS. After calibration, the dispersion of the extinction among neighbouring points in the second catalogue is reduced, regardless of whether reference targets are present locally. Weak structures are then more apparent. The extinction of high Galactic latitude targets is more tightly correlated with the dust emission, a property acquired from the first catalogue. A hierarchical inversion technique was applied to the two merged inter-calibrated catalogues to produce 3D extinction density maps corresponding to different volumes and maximum spatial resolution. The maximum resolution is 10pc for a 3000x3000x800pc3 volume around the Sun, and the maximum size of the maps is 10x10x0.8 kpc3 for a resolution of 50pc. Maps can be downloaded or used by means of on-line tools.

The forthcoming solar cycle (SC) 25 was beleived to be rather low when using the sunspot number (SN) as a measurement of the level of activity. The most popular prediction was made by the panel of NASA in 2019, including works based on extrapolations of dynamo-type models. We however discovered that using different observations to measure the level of polar regions activity several years before the start of SC25 and also after the start of the SC25 in 2020, the height of the SC25 could be high. The polar regions activity we considered seems related to the polar coronal holes (CH) activity and it is found significantly higher before the SC25 that it was before the SC24 and accordingly, we suggest that the SN cycle could indeed be much higher than during the SC24 that was a low SN height cycle.

Optical interferometers may not require a phase-stable optical link between the stations if instead sources of quantum-mechanically entangled pairs could be provided to them, enabling long baselines. We developed a new variation of this idea, proposing that photons from two different astronomical sources could be interfered at two decoupled stations. Interference products can then be calculated in post-processing or requiring only a slow, classical connection between stations. In this work, we investigated practical feasibility of this approach. We developed a Bayesian analysis method for the earth rotation fringe scanning technique and showed that in the limit of high signal-to-noise ratio it reproduced the results from a simple Fisher matrix analysis. We identify candidate stair pairs in the northern hemisphere, where this technique could be applied. With two telescopes with an effective collecting area of $\sim 2$ m$^2$, we could detect fringing and measure the astrometric separation of the sources at $\sim 10\,\mu$as precision in a few hours of observations, in agreement with previous estimates.

The linear point (LP) standard ruler was identified as the basis of a purely geometric method for exploiting the Baryon Acoustic Oscillations (BAO). The LP exploits the BAO feature imprinted in the galaxy two-point correlation function to measure cosmological distances independent of any specific cosmological model. We forecast the expected precision of future and ongoing spectroscopic galaxy surveys to measure distances leveraging the linear point. We investigate the cosmological implications of our forecasted results. We focus in particular on a relevant working example: the detection of the late-time cosmic acceleration independent of other cosmological probes. Our findings show that, even within the LCDM standard cosmological paradigm, estimated distances need to be reliable over a very wide parameter range in order to realize their maximum utility. This is particularly relevant if we aim to properly characterize cosmological tensions. The LP is a promising candidate approach to achieve this reliability. In contrast, widely employed procedures in BAO analysis estimate distances keeping fixed cosmological parameters to fiducial values close to cosmic-microwave-background constraints in flat-LCDM. It is unclear whether they are purely geometric methods. Moreover, they rely on untested extrapolations to explore the parameter space away from those fiducial flat-LCDM values. We recommend that all BAO methodologies be validated across the full range of models and parameters over which their results are quoted, first by means of linear predictions and then N-body simulations.

The detection of the Stochastic Gravitational Wave Background (SGWB) is essential for understanding black hole populations, especially for supermassive black hole binaries. The recent promising results from various Pulsar Timing Array (PTA) collaborations allude to an imminent detection. In this paper, we investigate the relative astrometric gravitational wave detection method, which can contribute to SGWB studies in the microhertz range. We consider the Roman Space Telescope and Gaia as candidates and quantitatively discuss the survey sensitivity in both the frequency and spatial domains. We emphasize the importance of survey specific constraints on performance estimates by considering mean field of view (FoV) signal subtraction and angular power spectrum binning. We conclude that if the SGWB is at a similar level as in PTA estimates, both Roman and Gaia have the potential to detect this frequency-domain power excess. However, both Roman and Gaia are subject to FoV limitations, and are unlikely to be sensitive to the spatial pattern of the SGWB.

Inspirals of stellar-mass compact objects into massive black holes, known as extreme mass ratio inspirals (EMRIs), are one of the key targets for upcoming space-based gravitational-wave detectors. In this paper we take the first steps needed to systematically incorporate the effect of external gravitating matter on EMRIs. We model the inspiral as taking place in the field of a Schwarzschild black hole perturbed by the gravitational field of a far axisymmetric distribution of mass enclosing the system. We take into account the redshift, frame-dragging, and quadrupolar tide caused by the enclosing matter, thus incorporating all effects to inverse third order of the characteristic distance of the enclosing mass. Then, we use canonical perturbation theory to obtain the action-angle coordinates and Hamiltonian for mildly eccentric precessing test-particle orbits in this background. Finally, we use this to efficiently compute mildly eccentric inspirals in this field and document their properties. This work shows the advantages of canonical perturbation theory for the modeling EMRIs, especially in the cases when the background deviates from the standard black hole fields.

Black hole spectroscopy with gravitational waves is an important tool to measure the mass and spin of astrophysical black holes and to test their Kerr nature. Next-generation ground- and space-based detectors will observe binary black hole mergers with large signal-to-noise ratios and perform spectroscopy routinely. It was recently shown that small perturbations due, e.g., to environmental effects (the "flea") to the effective potential governing gravitational-wave generation and propagation in black hole exteriors (the "elephant") can lead to arbitrarily large changes in the black hole's quasinormal spectrum, including the fundamental mode, which is expected to dominate the observed signal. This raises an important question: is the black hole spectroscopy program robust against perturbations? We clarify the physical behavior of time-domain signals under small perturbations in the potential, and we show that changes in the amplitude of the fundamental mode in the prompt ringdown signal are parametrically small. This implies that the fundamental quasinormal mode extracted from the observable time-domain signal is stable against small perturbations. The stability of overtones deserves further investigation.

Binaries of relatively massive black holes like GW190521 have been proposed to form in dense gas environments, such as the disks of Active Galactic Nuclei (AGNs), and they might be associated with transient electromagnetic counterparts. The interactions of this putative environment with the binary could leave a significant imprint at the low gravitational wave frequencies observable with the Laser Interferometer Space Antenna (LISA). We show that LISA will be able to detect up to ten GW190521-like black hole binaries, with sky position errors $\lesssim1$ deg$^2$. Moreover, it will measure directly various effects due to the orbital motion around the supermassive black hole at the center of the AGN, especially the Doppler modulation and the Shapiro time delay. Thanks to a careful treatment of their frequency domain signal, we were able to perform the full parameter estimation of Doppler and Shapiro-modulated binaries as seen by LISA. We find that the Doppler and Shapiro effects will allow for measuring the AGN parameters (radius and inclination of the orbit around the AGN, central black hole mass) with up to percent-level precision. Properly modeling these low-frequency environmental effects is crucial to determine the binary formation history, as well as to avoid biases in the reconstruction of the source parameters and in tests of general relativity with gravitational waves.

In this work we shall study the late-time dynamics of several $F(R)$ gravity models. By appropriately expressing the field equations in terms of the redshift and of a statefinder function, we shall solve numerically the field equations using appropriate physical motivated initial conditions. We consider models which, by construction, are described by a nearly $R^2$-model at early epochs and we fine tune the parameters to achieve viability and compatibility with the latest Planck constraints at late times. Such models provide a unified description of inflation and dark energy era and notably a common feature of all the models is the presence of dark energy oscillations. Furthermore, we show that, in contrast to general relativistic fluids and scalar field descriptions, a large spectrum of different dark energy physics is generated by simple $F(R)$ gravity models, varying from phantom, to nearly de Sitter and to quintessential dark energy eras.

Recent advances with space navigation technologies developed by NASA in space-based atomic clocks and pulsar X-ray navigation combined with past successes in autonomous navigation using optical imaging, brings to the forefront the need to compare space navigation using optical, radiometric, and pulsar-based measurements using a common set of assumptions and techniques. This review article examines these navigation data types in two different ways. First, a simplified deep space orbit determination problem is posed that captures key features of the dynamics and geometry, and then each data type is characterized for its ability to solve for the orbit. The data types are compared and contrasted using a semi-analytical approach with geometric dilution of precision techniques. The results provide useful parametric insights into the strengths of each data type. In the second part of the paper, a high-fidelity, Monte Carlo simulation of a Mars cruise, approach, and entry navigation problem is studied. The results found complement the semi-analytic results in the first part, and illustrate specific issues such as each data type's quantitative impact on solution accuracy and their ability to support autonomous delivery to a planet.

Using the method we proposed previously, in which the spin-$1/2$ massless spinors (Dirac-Weyl fields) are regarded as basic units, we study the Weyl double copy for vacuum solutions with a cosmological constant. The result explicitly demonstrates that the single and zeroth copy satisfy conformally invariant field equations on conformally flat spacetime, based on the exact non-twisting vacuum type $N$ and vacuum type $D$ cases. Furthermore, irrespective of the presence of a cosmological constant, we show that the zeroth copy not only connects gravity fields with the single copy but also connects degenerate electromagnetic fields with Dirac-Weyl fields in the curved spacetime. Moreover, the study shows that the zeroth copy plays a vital role in time-dependent radiation spacetimes. In particular, for Robinson-Trautman ($\Lambda$) gravitational waves, as opposed to the single copy, we find that the zeroth copy carries extra information to indicate whether the sources of associated gravitational waves are time-like, null, or space-like.

Turbulence plays an integral role in astrophysical phenomena, including core-collapse supernovae (CCSN). Unfortunately, current simulations must resort to using subgrid models for turbulence treatment, as direct numerical simulations (DNS) are too expensive to run. However, subgrid models used in CCSN simulations lack accuracy compared to DNS results. Recently, Machine Learning (ML) has shown impressive prediction capability for turbulence closure. We have developed a physics-informed, deep convolutional neural network (CNN) to preserve the realizability condition of Reynolds stress that is necessary for accurate turbulent pressure prediction. The applicability of the ML model was tested for magnetohydrodynamic (MHD) turbulence subgrid modeling in both stationary and dynamic regimes. Our future goal is to utilize our ML methodology within the MHD CCSN framework to investigate the effects of accurately-modeled turbulence on the explosion rate of these events.

I define here a novel function on a modeled space of gravitational-wave signals, before studying its properties as a statistic for detection, as an objective function for identification, and as an effective likelihood function for inference. The main motivation behind this work is the open data-analysis problem for signals from extreme-mass-ratio inspirals, which is severely hindered by the presence of strong non-local parameter degeneracy in the signal space. I demonstrate the utility of the proposed function for the analysis of such signals, and suggest various possible directions for future research.

We study the parametric resonance excitation of the electromagnetic field by a gravitational wave. We show that there is narrow band resonance. For an electromagnetic field in the vacuum the resonance occurs only in the second band, and its strength is thus suppressed by two powers of the Planck mass. On the other hand, in the case of an electromagnetic field in a medium with the speed of light smaller than 1 (in natural units), there is a band of Fourier modes which undergo resonance in the first band.

We present the detection potential for the diffuse supernova neutrino background (DSNB) at the Jiangmen Underground Neutrino Observatory (JUNO), using the inverse-beta-decay (IBD) detection channel on free protons. We employ the latest information on the DSNB flux predictions, and investigate in detail the background and its reduction for the DSNB search at JUNO. The atmospheric neutrino induced neutral current (NC) background turns out to be the most critical background, whose uncertainty is carefully evaluated from both the spread of model predictions and an envisaged \textit{in situ} measurement. We also make a careful study on the background suppression with the pulse shape discrimination (PSD) and triple coincidence (TC) cuts. With latest DSNB signal predictions, more realistic background evaluation and PSD efficiency optimization, and additional TC cut, JUNO can reach the significance of 3$\sigma$ for 3 years of data taking, and achieve better than 5$\sigma$ after 10 years for a reference DSNB model. In the pessimistic scenario of non-observation, JUNO would strongly improve the limits and exclude a significant region of the model parameter space.

We propose a new four-parameter entropy function that generalizes the Tsallis, R\'{e}nyi, Barrow, Sharma-Mittal, Kaniadakis and Loop Quantum Gravity entropies for suitable limits of the parameters. Consequently, we address the early and late universe cosmology corresponding to the proposed four-parameter entropy function. As a result, it turns out that the entropic cosmology from the generalized entropy function can unify the early inflation to the late dark energy era of the universe. In such a unified scenario, we find that -- (1) the inflation era is described by a quasi de-Sitter evolution of the Hubble parameter, which has an exit at around 58 e-folding number, (2) the inflationary observable quantities like the spectral index for primordial scalar perturbation and the tensor-to-scalar ratio are simultaneously compatible with the recent Planck data, and (3) regarding the late time cosmology, the dark energy EoS parameter is found to be consistent with the Planck result for the same values of the entropy parameters that lead to the viable inflation during the early universe. Furthermore, we show that the entropic cosmology from the proposed entropy function is equivalent to holographic cosmology, where the respective holographic cut-offs are determined in terms of either particle horizon and its derivative or future horizon and its derivative.

Field distortion is widespread in imaging systems. If it cannot be measured and corrected well, it will affect the accuracy of photogrammetry. To this end, we proposed a general field distortion model based on Fredholm integration, which uses a reconstructed high-resolution reference point spread function (PSF) and two sets of 4-variable polynomials to describe an imaging system. The model includes the point-to-point positional distortion from the object space to the image space and the deformation of the PSF so that we can measure an actual field distortion with arbitrary accuracy. We also derived the formula required for correcting the sampling effect of the image sensor. Through numerical simulation, we verify the effectiveness of the model and reconstruction algorithm. This model will have potential application in high-precision image calibration, photogrammetry and astrometry.

We present ${\tt NRPMw}$, an analytical model of gravitational-waves from neutron star merger remnants informed using 618 numerical relativity (NR) simulations. ${\tt NRPMw}$ is designed in the frequency domain using a combination of complex Gaussian wavelets. The wavelet's parameters are calibrated to equations of state (EOS) insensitive relations from NR data. The NR simulations are computed with 21 EOS (7 of which are finite-temperature microphysical models, and 3 of which contain quark phase transitions or hyperonic degrees of freedom) and span total binary masses $M\in[2.4,3.4]~{\rm M}_\odot$, mass ratios up to $q=2$, and (nonprecessing) dimensionless spins magnitudes up to ${0.2}$. The theoretical uncertainties of the EOS-insensitive relations are incorporated in ${\tt NRPMw}$ using recalibration parameters that enhance the flexibility and accuracy of the model. ${\tt NRPMw}$ is NR-faithful with fitting factors ${\gtrsim}0.9$ computed on an independent validation set of 102 simulations.