Papers for Tuesday, Mar 22 2022

We use the TNG100 simulation of the IllustrisTNG project to investigate the stellar specific angular momenta ($j_{\ast}$) of $\sim$12~000 central galaxies at $z=0$ in a full cosmological context, with stellar masses ($M_{\ast}$) ranging from $10^{9}$ to $10^{12} \, {\rm M}_{\odot}$. We find that the $j_{\ast}$--$M_{\ast}$ relations for early-type and late-type galaxies in IllustrisTNG are in good overall agreement with observations, and that these galaxy types typically `retain' $\sim$10--20 and $\sim$50--60 per cent of their host haloes' specific angular momenta, respectively, with some dependence on the methodology used to measure galaxy morphology. We present results for \textit{kinematic} as well as \textit{visual-like} morphological measurements of the simulated galaxies. Next, we explore the scatter in the $j_{\ast}$--$M_{\ast}$ relation with respect to the spin of the dark matter halo and the mass of the supermassive black hole (BH) at the galactic centre. We find that galaxies residing in faster-spinning haloes, as well as those hosting less massive BHs, tend to have a higher specific angular momentum. We also find that, at fixed galaxy or halo mass, halo spin and BH mass are anticorrelated with each other, probably as a consequence of more efficient gas flow toward the galactic centre in slowly rotating systems. Finally, we show that halo spin plays an important role in determining galaxy sizes -- larger discs form at the centres of faster-rotating haloes -- although the trend breaks down for massive galaxies with $M_{\ast} \gtrsim 10^{11} \, {\rm M}_{\odot}$, roughly the mass scale at which a galaxy's stellar mass becomes dominated by accreted stars.

We explore the properties of Milky Way subhalos in self-interacting dark matter models for moderate cross sections of 1 to 5 cm$^2$g$^{-1}$ using high-resolution zoom-in N-body simulations. We include the gravitational potential of a baryonic disk and bulge matched to the Milky Way, which is critical for getting accurate predictions. The predicted number and distribution of subhalos within the host halo are similar for 1 and 5 cm$^2$g$^{-1}$ models, and they agree with observations of Milky Way satellite galaxies only if subhalos with peak circular velocity over all time > 4.5 km/s are able to form galaxies. We do not find distinctive signatures in the pericenter distribution of the subhalos that could help distinguish the models. Using an analytic model to extend the simulation results, we are able to show that subhalos in models with cross sections between 1 and 5 cm$^2$g$^{-1}$ are not dense enough to match the densest ultra-faint and classical dwarf spheroidal galaxies in the Milky Way. This motivates velocity-dependent cross sections with values larger than 5 cm$^2$g$^{-1}$ at the velocities relevant for the satellites such that core collapse would occur in some of the ultra-faint and classical dwarf spheroidals.

The ratio of baryonic to dark matter in present-day galaxies constrains galaxy formation theories and can be determined empirically via the baryonic Tully-Fisher relation (BTFR), which compares a galaxy's baryonic mass (M$_{bary}$) to its maximum rotation velocity (V$_{max}$). The BTFR is well-determined at M$_{bary}>10^8$ M$_{\odot}$, but poorly constrained at lower masses due to small samples and the challenges of measuring rotation velocities in this regime. For 25 galaxies with high-quality data and M$_{bary}<\sim10^8$ M$_{\odot}$, we estimate M$_{bary}$ from infrared, optical, and HI observations and Vmax from the HI gas rotation. Many of the V$_{max}$ values are lower limits because the velocities are still rising at the edge of the detected HI disks; consequently, most of our sample has lower velocities than expected from extrapolations of the BTFR at higher masses. To estimate V$_{max}$, we map each galaxy to a dark matter halo assuming density profiles with and without cores, and find that the cored profiles match the data better. When we compare the V$_{max}$ values derived from the cored density profiles to our M$_{bary}$ measurements, we find a turndown of the BTFR at low masses that is consistent with CDM predictions and implying baryon fractions of 1-10% of the cosmic value. Although we are limited by the sample size and assumptions inherent in mapping measured rotational velocities to theoretical rotation curves, our results suggest that the galaxy formation efficiency drops at masses below M$_{bary}\sim10^8$ M$_{\odot}$, corresponding to M$_{200}\sim10^{10}$ M$_{\odot}$.

The habitability of planets around M dwarfs ($\lesssim 0.5 M_\odot$) can be affected by the XUV (X rays + extreme UV) emission of these stars, with flares occasionally increasing the XUV flux by more than 2 orders of magnitude above quiescent levels. This wavelength range can warm and ionize terrestrial planets' upper atmospheres, which expands the planetary radius and promotes atmospheric loss. In this work, we study the contribution of the XUV flux due to flares on the atmospheric escape of Earth-like planets orbiting M dwarfs through numerical simulations. We considered the first Gyr of planets with initial surface water abundances between 1 and 10 terrestrial oceans (TO), a small primordial hydrogen envelope ($\le$ $10^{-3}$~$M_{\oplus}$), and with host star masses between 0.2 and 0.6 $M_{\odot}$. In this parameter range, we find that flares can remove up to two TO more than nonflaring stars, which, in some cases, translates to a doubling of the total water loss. We also find that flaring can increase atmospheric oxygen partial pressures by hundreds of bars in some cases. These results were obtained by adding a new module for flares to the \vplanet software package and upgrading its atmospheric escape module to account for Roche lobe overflow and radiation/recombination-limited escape.

We study a large sample of dwarf galaxies using archival Chandra X-ray observations, with the aim of detecting accreting intermediate-mass black holes (IMBHs). IMBHs are expected to inhabit dwarf galaxies and to produce specific signatures in terms of luminosity and X-ray spectra. We report the discovery of an X-ray source associated with an Abell 85 dwarf galaxy that fits the IMBH description. The stellar mass of the host galaxy is estimated to be 2 $\times$ 10$^8$ $M_\odot$, which makes it one of the least massive galaxies to potentially host an accreting black hole. The source is detected in the soft band, under 1 keV, while undetected at higher energies. The X-ray luminosity is $\approx$ 10$^{41}$ erg s$^{-1}$, making it almost three orders of magnitude more luminous than the most luminous stellar-mass supersoft emitters. From the galaxy stellar mass vs. black hole mass relation, we estimate the mass to be within the intermediate regime. Another method that resulted in an intermediate mass relies on the fact that supersoft emission is expected to be associated with high accretion rates, approaching the Eddington limit. We suggest that the observed offset of the X-ray source from the galactic center ($\approx$ 1.8 kpc) is due to galaxy interactions, and we present evidence from the literature that supports the relation between black hole activity and galaxy interactions.

Photonic spectrographs offer a highly miniaturized, flexible, and stable on-chip solution for astronomical spectroscopy and can be used for various science cases such as determining the atmospheric composition of exoplanets to understand their habitability, formation, and evolution. Arrayed Waveguide Gratings (AWGs) have shown the best promise to be used as an astrophotonic spectrograph. We developed a publically-available tool to conduct a preliminary examination of the capability of the AWGs in spectrally resolving exoplanet atmospheres. We derived the Line-Spread-Function (LSF) as a function of wavelength and the Full-Width-at-Half-Maximum (FWHM) of the LSF as a function of spectral line width to evaluate the response of a discretely- and continuously-sampled low-resolution AWG (R $\sim$ 1000). We observed that the LSF has minimal wavelength dependence ($\sim$5\%), irrespective of the offset with respect to the center-wavelengths of the AWG channels, contrary to the previous assumptions. We further confirmed that the observed FWHM scales linearly with the emission line width. Finally, we present simulated extraction of a sample molecular absorption spectrum with the discretely- and continuously-sampled low-resolution AWGs. From this, we show that while the discrete AWG matches its expected resolving power, the continuous AWG spectrograph can, in principle, achieve an effective resolution significantly greater ($\sim$ 2x) than the discrete AWG. This detailed examination of the AWGs will be foundational for future deployment of AWG spectrographs for astronomical science cases such as exoplanet atmospheres.

Stars are formed by gravitational collapse, spontaneously or, in some cases under the constructive influence of nearby massive stars, out of molecular cloud cores. Here we present an observational diagnosis of such triggered formation processes in the prominent \ion{H}{2} region Sh\,2-142, which is associated with the young star cluster NGC\,7380, and with some bright-rimmed clouds as the signpost of photoionization of molecular cloud surfaces. Using near- (2MASS) and mid-infrared (WISE) colors, we identified candidate young stars at different evolutionary stages, including embedded infrared sources having spectral energy distributions indicative of active accretion. We have also used data from our optical observations to be used in SEDs, and from Gaia EDR3 to study the kinematics of young objects. With this young stellar sample, together with the latest CO line emission data (spectral resolution $\sim 0.16$~km~s$^{-1}$, sensitivity $\sim 0.5$~K), a positional and ageing sequence relative to the neighboring cloud complex, and to the bright-rimmed clouds, is inferred. The propagating stellar birth may be responsible, at least partially, for the formation of the cluster a few million years ago, and for the ongoing activity now witnessed in the cloud complex.

Relativistic jets are believed to have a substantial impact on the gas dynamics and evolution of the interstellar medium (ISM) of their host galaxies. In this paper, we aim to draw a link between the simulations and the observable signatures of jet-ISM interactions by analyzing the emission morphology and gas kinematics resulting from jet-induced shocks in simulated disc and spherical systems. We find that the jet-induced laterally expanding forward shock of the energy bubble sweeping through the ISM causes large-scale outflows, creating shocked emission and high-velocity dispersion in the entire nuclear regions ($\sim2$ kpcs) of their hosts. The jetted systems exhibit larger velocity widths (> 800 km/s), broader Position-Velocity maps and distorted symmetry in the disc's projected velocities than systems without a jet. We also investigate the above quantities at different inclination angles of the observer with respect to the galaxy. Jets inclined to the gas disc of its host are found to be confined for longer times, and consequently couple more strongly with the disc gas. This results in prominent shocked emission and high-velocity widths, not only along the jet's path, but also in the regions perpendicular to them. Strong interaction of the jet with a gas disc can also distort its morphology. However, after the jets escape their initial confinement, the jet-disc coupling is weakened, thereby lowering the shocked emission and velocity widths.

We present Atacama Large Millimeter/submillimeter Array observations of multiple CO(1-0), $^{13}$CO(1-0), and C$^{18}$O(1-0) lines and 2.9 mm and 1.3 mm continuum emission toward the nearby interacting luminous infrared galaxy NGC 3110, supplemented with similar spatial resolution H$\alpha$, 1.4GHz continuum, and $K$-band data. We estimate the typical CO-to-H$_2$ conversion factor of 1.7 $M_{\odot}$ (K km s$^{-1}$ pc$^2$)$^{-1}$ within the disk using LTE-based and dust-based H$_2$ column densities, and measure the 1-kpc scale surface densities of star formation rate ($\Sigma_{\rm SFR}$), super star clusters ($\Sigma_{\rm SSC}$), molecular gas mass, and star formation efficiency (SFE) toward the entire gas disk. These parameters show a peak at the southern part of the southern spiral arm (SFE $\sim$ 10$^{-8.2}$ yr$^{-1}$, $\Sigma_{\rm SFR}$ $\sim$ 10$^{-0.6}$ $M_{\odot}$ kpc$^{-2}$ yr$^{-1}$, $\Sigma_{\rm SSC}$ $\sim$ 6.0 kpc$^{-2}$), which is likely attributed to the on-going tidal interaction with the companion galaxy MCG-01-26-013, as well as toward the circumnuclear region. We also find that thermal free-free emission contributes to a significant fraction of the millimeter continuum emission at the southern peak position. Those measurements imply that the peak of the southern arm is an active and young star-forming region, whereas the central part of NGC 3110 is a site of long-continued star formation. We suggest that, during the early stage of the galaxy-galaxy interaction with large mass ratio that in NGC 3110, fragmentation along the main galaxy's arms is an important driver of merger-induced star formation and massive gas inflow results in dusty nuclear starbursts.

The surface lithium abundance A(Li) of warm metal-poor dwarf stars exhibits a narrow plateau down to [Fe/H]~-2.8 dex, while at lower metallicities the average value drops by 0.3 dex with a significant star-by-star scatter (called lithium meltdown). This behaviour is in conflict with predictions of standard stellar evolution models calculated with the initial A(Li) provided by the standard Big Bang nucleosynthesis. The lower red giant branch (LRGB) stars provide a complementary tool to understand the initial A(Li) distribution in metal-poor stars. We have collected a sample of high-resolution spectra of 58 LRGB stars spanning a range of [Fe/H] between ~ -7.0 dex and ~ -1.3 dex. The LRGB stars display an A(Li) distribution clearly different from that of the dwarfs, without signatures of a meltdown and with two distinct components: (a) a thin A(Li) plateau with an average A(Li)=1.09+-0.01 dex (sigma=0.07 dex), and (b) a small fraction of Li-poor stars with A(Li) lower than ~0.7 dex. The A(Li) distribution observed in LRGB stars can be reconciled with an initial abundance close to the cosmological value, by including an additional chemical element transport in stellar evolution models. The required efficiency of this transport allows us to match also the Spite plateau lithium abundance measured in the dwarfs. The emerging scenario is that all metal-poor stars formed with the same initial A(Li) but those that are likely the product of coalescence or that experienced binary mass transfer and show lower A(Li) . We conclude that A(Li) in LRGB stars is qualitatively compatible with the cosmological A(Li) value and that the meltdown observed in dwarf stars does not reflect a real drop of the abundance at birth.

Differential Interferometry allows to obtain the differential visibility and phase, in addition to the spectrum. The differential phase contains important information about the structure and motion of stellar photosphere such as stellar spots and non-radial pulsations, and particularly the rotation. Thus, this interferometric observable strongly helps to constrain the stellar fundamental parameters of fast rotators. The spectro-astrometry mainly uses the photocentre displacements, which is a first approximation of the differential phase, and is applicable only for unresolved or marginally objects. We study here the sensitivity of relevant stellar parameters to the simulated photocentres using the SCIROCCO code: a semi-analytical algorithm dedicated to fast rotators, applied to two theoretical modeling stars based on Achernar and Regulus, in order to classify the importance of these parameters and their impact on the modeling. We compare our simulations with published VLTI/AMBER data. This current work sets the limits of application of photocentre displacements to fast rotators, and under which conditions we can use the photocentres and/or the differential phase, through a pre-established physical criterion. To validate our theoretical study, we apply our method of analysis on observed data of the edge-on fast rotator Regulus. For unresolved targets, with a visibility $V\sim 1$, the photocentre can constrain the main stellar fundamental parameters of fast rotators, whereas from marginally resolved objects ($0.8 \leq V < 1$), mainly the rotation axis position angle ($\rm PA_{\rm rot}$) can be directly deduced from the vectorial photocentre displacement, which is very important for young cluster studies.

We generalize the symmetron screening mechanism by allowing for an explicit symmetry breaking of the symmetron $\phi^4$ potential. A coupling to matter of the form $A(\phi)=1+\frac{\phi^2}{M^2}$ leads to an explicitly broken symmetry with effective potential $V_{eff}(\phi)=-\mu^2 (1-\frac{\rho}{\mu^2 M^2})\phi^2 +\frac{\lambda}{2}\phi^4 + 2 \varepsilon \phi^3+\frac{\lambda}{2}\eta^4$. Due to the explicit symmetry breaking induced by the cubic term we call this field the 'asymmetron'. For large matter density $\rho>\rho_*\equiv \mu^2M^2+\frac{9}{4}\varepsilon\eta M^2$ the effective potential has a single minimum at $\phi=0$ leading to restoration of General Relativity as in the usual symmetron screening mechanism. For low matter density however, there is a false vacuum and a single true vacuum due to the explicit symmetry breaking. This is expected to lead to an unstable network of domain walls with slightly different value of the gravitational constant $G$ on each side of the wall. This network would be in constant interaction with matter overdensities and would lead to interesting observational signatures which could be detected as gravitational and expansion rate transitions in redshift space. Such a gravitational transition has been recently proposed for the resolution of the Hubble tension.

Context. Interstellar scintillation (ISS) of pulsar emission can be used both as a probe of the ionised interstellar medium (IISM) and cause corruptions in pulsar timing experiments. Of particular interest are so-called scintillation arcs which can be used to measure time-variable interstellar scattering delays directly, potentially allowing high-precision improvements to timing precision. Aims. The primary aim of this study is to carry out the first sizeable and self-consistent census of diffractive pulsar scintillation and scintillation-arc detectability at low frequencies, as a primer for larger-scale IISM studies and pulsar-timing related propagation studies with the LOw-Frequency ARray (LOFAR) High Band Antennae (HBA). Results. In this initial set of 31 sources, 15 allow full determination of the scintillation properties; nine of these show detectable scintillation arcs at 120-180 MHz. Eight of the observed sources show unresolved scintillation; and the final eight don't display diffractive scintillation. Some correlation between scintillation detectability and pulsar brightness and dispersion measure is apparent, although no clear cut-off values can be determined. Our measurements across a large fractional bandwidth allow a meaningful test of the frequency scaling of scintillation parameters, uncorrupted by influences from refractive scintillation variations. Conclusions. Our results indicate the powerful advantage and great potential of ISS studies at low frequencies and the complex dependence of scintillation detectability on parameters like pulsar brightness and interstellar dispersion. This work provides the first installment of a larger-scale census and longer-term monitoring of interstellar scintillation effects at low frequencies.

Global 21-cm experiments require exquisitely precise calibration of the measurement systems in order to separate the weak 21-cm signal from Galactic and extragalactic foregrounds as well as instrumental systematics. Hitherto, experiments aiming to make this measurement have concentrated on measuring this signal using the single element approach. However, an alternative approach based on interferometers with short baselines is expected to alleviate some of the difficulties associated with a single element approach such as precision modelling of the receiver noise spectrum. Short spacing Interferometer Telescope probing cosmic dAwn and epoch of ReionisAtion (SITARA) is a short spacing interferometer deployed at the Murchison Radio-astronomy Observatory (MRO). It is intended to be a prototype or a test-bed to gain a better understanding of interferometry at short baselines, and develop tools to perform observations and data calibration. In this paper, we provide a description of the SITARA system and its deployment at the MRO, and discuss strategies developed to calibrate SITARA. We touch upon certain systematics seen in SITARA data and their modelling. We find that SITARA has sensitivity to all sky signals as well as non-negligible noise coupling between the antennas. It is seen that the coupled receiver noise has a spectral shape that broadly matches the theoretical calculations reported in prior works. We also find that when appropriately modified antenna radiation patterns taking into account the effects of mutual coupling are used, the measured data are well modelled by the standard visibility equation.

Exposed scarps images and ice-penetrating radar measurements in the North Polar Layered Deposits (NPLD) of Mars show alternating layers that provide an archive of past climate oscillations, that are thought to be linked to orbital variations, akin to Milankovitch cycles on Earth. We use the Laboratoire de Meteorologie Dynamique (LMD) Martian Global Climate Model (GCM) to study paleoclimate states to enable a better interpretation of the NPLD physical and chemical stratigraphy. When a tropical ice reservoir is present, water vapor transport from the tropics to the poles at low obliquity is modulated by the intensity of summer. At times of low and relatively constant obliquity, the flux still varies due to other orbital elements, promoting polar layer formation. Ice migrates from the tropics towards the poles in two stages. First, when surface ice is present in the tropics, and second, when the equatorial deposit is exhausted, from ice that was previously deposited in mid-high latitudes. The polar accumulation rate is significantly higher when tropical ice is available, forming thicker layers per orbital cycle. However, the majority of the NPLD is sourced from ice that temporary resided in the mid-high latitudes and the layers become thinner as the source location moves poleward. The migration stages imprint different D/H ratios in different sections in the PLDs. The NPLD is isotopically depleted compared to the SPLD in all simulations. Thus we predict the D/H ratio of the atmosphere in contact with NPLD upper layers is biased relative to the average global ice reservoirs.

By compiling a comprehensive census of literature studies, we investigate the evolution of the Main Sequence (MS) of star-forming galaxies (SFGs) in the widest range of redshift (0 < z < 6) and stellar mass ($10^{8.5}-10^{11.5}$ $M_{\odot}$) ever probed. We convert all observations to a common calibration and find a remarkable consensus on the variation of the MS shape and normalization across cosmic time. The relation exhibits a curvature towards the high stellar masses at all redshifts. The best functional form is governed by two parameters: the evolution of the normalization and the turn-over mass ($M_0(t)$), which both evolve as a power law of the Universe age. The turn-over mass determines the MS shape. It marginally evolves with time, making the MS slightly steeper towards $z\sim4-6$. At stellar masses below $M_0(t)$, SFGs have a constant specific SFR (sSFR), while above $M_0(t)$ the sSFR is suppressed. We find that the MS is dominated by central galaxies. This allows transforming $M_0(t)$ into the corresponding host halo mass. This evolves as the halo mass threshold between cold and hot accretion regimes, as predicted by the theory of accretion, where the central galaxy is fed or starved of cold gas supply, respectively. We, thus, argue that the progressive MS bending as a function of the Universe age is caused by the lower availability of cold gas in halos entering the hot accretion phase, in addition to black hole feedback. We also find qualitatively the same trend in the largest sample of star-forming galaxies provided by the IllustrisTNG simulation, although we note still a large discrepancy with respect to observations.

Through its magnetic activity, the Sun governs the conditions in Earth's vicinity, creating space weather events, which have drastic effects on our space- and ground-based technology. One of the most important solar magnetic features creating the space weather is the solar wind, that originates from the coronal holes (CHs). The identification of the CHs on the Sun as one of the source regions of the solar wind is therefore crucial to achieve predictive capabilities. In this study, we used an unsupervised machine learning method, $k$-means, to pixel-wise cluster the passband images of the Sun taken by the Atmospheric Imaging Assembly on {\it the Solar Dynamics Observatory} (AIA/SDO) in 171 \AA, 193 \AA\,, and 211 \AA\,in different combinations. Our results show that the pixel-wise $k$-means clustering together with systematic pre- and post-processing steps provides compatible results with those from complex methods, such as CNNs. More importantly, our study shows that there is a need for a CH database that a consensus about the CH boundaries are reached by observers independently. This database then can be used as the "ground truth", when using a supervised method or just to evaluate the goodness of the models.

Planets with non-zero obliquity and/or orbital eccentricity experience seasonal variations of stellar irradiation at local latitudes. The extent of the atmospheric response can be crudely estimated by the ratio between the orbital timescale and the atmospheric radiative timescale. Given a set of atmospheric parameters, we show that this ratio depends mostly on the stellar properties and is independent of orbital distance and planetary equilibrium temperature. For Jupiter-like atmospheres, this ratio is $\ll1$ for planets around very-low-mass M dwarfs and $\gtrsim1$ when the stellar mass is greater than about 0.6 solar mass. Complications can arise from various factors, including varying atmospheric metallicity, clouds, and atmospheric dynamics. Given the eccentricity and obliquity, the seasonal response is expected to be systematically weaker for gaseous exoplanets around low-mass stars and stronger for those around more massive stars. The amplitude and phase lag of atmospheric seasonal variations as a function of host stellar mass are quantified by idealized analytic models. At the infrared emission level in the photosphere, the relative amplitudes of thermal flux and temperature perturbations are negligible, and their phase lags are closed to $-90^{\circ}$ for Jupiter-like planets around very-low-mass stars. The relative amplitudes and phase lags increase gradually with increasing stellar mass. With a particular stellar mass, the relative amplitude and phase lag decrease from low to high infrared optical depth. We also present numerical calculations for a better illustration of the seasonal behaviors. Lastly, we discuss implications for the atmospheric circulation and future atmospheric characterization of exoplanets in systems with different stellar masses.

Thermal tides in the Martian atmosphere are analyzed using temperature profiles retrieved from nadir observations obtained by the TIRVIM Fourier-spectrometer, part of the Atmospheric Chemistry Suite (ACS) onboard the ExoMars Trace Gas Orbiter (TGO). The data is selected near the northern summer solstice at solar longitude (LS) 75{\deg}-105{\deg} of Martian Year (MY) 35. The observations have a full local time coverage, which enables analyses of daily temperature anomalies. The observed zonal mean temperature is lower by 4-6K at ~100Pa, but higher towards the summer pole, compared to the LMD Mars General Circulation Model (GCM). Wave mode decomposition shows dominant diurnal tide and important semi-diurnal tide and diurnal Kelvin wave, with maximal amplitudes of 5K, 3K, and 2.5K, respectively, from tens to hundreds of Pa. The results generally agree well with the LMD Mars GCM, but with noticeable earlier phases of diurnal (~1h) and semi-diurnal (~3h) tides.

Observations of brown dwarfs and relatively isolated young extrasolar giant planets have provided unprecedented details to probe atmospheric dynamics in a new regime. Questions about mechanisms governing global circulation remain to be addressed. Previous studies have shown that small-scale, randomly varying thermal perturbations resulting from interactions between convection and the overlying stratified layers can drive zonal jet streams, waves, and turbulence. Here, we improve upon our previous general circulation model by using a two-stream grey radiative transfer scheme to represent more realistic heating and cooling rates. We examine the formation of zonal jets and their time evolution, and vertical mixing of passive tracers including clouds and chemical species. Under relatively weak radiative and frictional dissipation, robust zonal jets with speeds up to a few hundred $\rm m\;s^{-1}$ are typical outcomes. The off-equatorial jets tend to be pressure-independent while the equatorial jets exhibit significant vertical wind shear. Models with strong dissipation inhibit jet formation and have isotropic turbulence in off-equatorial regions. Quasi-periodic oscillations of the equatorial flow with periods ranging from tens of days to months are prevalent at relatively low atmospheric temperatures. Sub-micron cloud particles can be transported to several scale heights above the condensation level, while larger particles form thinner layers. Cloud decks are inhomogeneous near their cloud tops. Chemical tracers with chemical timescales $>10^5$ s can be driven out of equilibrium. The equivalent vertical diffusion coefficients, $K_{\mathrm{zz}}$, for the global-mean tracer, are diagnosed from our models and are typically on the order of $1\sim10^2\rm m^2\;s^{-1}$. Finally, we derive an analytic estimation of $K_{\mathrm{zz}}$ for different types of tracers under relevant conditions.

We present a detailed study of the highly obscured active galaxy NGC 4507, performed using four Nuclear Spectroscopic Telescope Array (NuSTAR) observations carried out between May and August in 2015 (~ 130 ks in total). Using various phenomenological and physically motivated torus models, we explore the properties of the X-ray source and those of the obscuring material. The primary X-ray emission is found to be non-variable, indicating a stable accretion during the period of the observations. We find the equatorial column density of the obscuring materials to be ~ 2 x 10^24 cm^-2 while the line of sight column density to be ~ 7 - 8 x 10^23 cm^-22. The source is found to be deeply buried with the torus covering factor ~ 0.85. We observe variability in the line-of-sight column density on a timescale of < 35 days. The covering factor of the Compton-Thick material is found to be ~ 0.35, in agreement with the results of recent X-ray surveys. From the variability of the line-of-sight column density, we estimate that the variable absorbing material is likely located either in the BLR or in the torus.

Emerging high redshift cosmological probes, in particular quasars (QSOs), show a preference for larger matter densities, $\Omega_{m} \approx 1$, within the flat $\Lambda$CDM framework. Here, using the Risaliti-Lusso relation for standardizable QSOs, we demonstrate that the QSOs recover the \textit{same} Planck-$\Lambda$CDM Universe as Type Ia supernovae (SN), $\Omega_m \approx 0.3$ at lower redshifts $ 0 < z \lesssim 0.7$, before transitioning to an Einstein-de Sitter Universe ($\Omega_m =1$) at higher redshifts $z \gtrsim 1$. We illustrate the same trend, namely increasing $\Omega_{m}$ and decreasing $H_0$ with redshift, in SN but poor statistics prevent a definitive statement. We explain physically why the trend is expected in the flat $\Lambda$CDM cosmology, illustrate the intrinsic bias and non-Gaussian tails with mock Pantheon data, and identify a similar trend in BAO below $z=1$. Our results highlight an intrinsic bias in the flat $\Lambda$CDM Universe, whereby $\Omega_m$ increases, $H_0$ decreases and $S_8$ increases with effective redshift, thus providing a new perspective on $\Lambda$CDM tensions; even in a Planck-$\Lambda$CDM Universe the current tensions might have been expected.

We continue our series of papers on intergalactic medium (IGM) tracers using quasi-stellar objects (QSOs), having examined gamma-ray bursts (GRBs) and blazars in earlier studies. We have estimated the IGM properties of hydrogen column density (Nhxigm), temperature and metallicity using XMM-Newton QSO spectra over a large redshift range, with a collisional ionisation equilibrium (CIE) model for the ionised plasma. The Nhxigm parameter results were robust with respect to intrinsic power laws, spectral counts, reflection hump and soft excess features. There is scope for a luminosity bias given both luminosity and Nhxigm scale with redshift, but we find this unlikely given the consistent IGM parameter results across the other tracer types reviewed. The impact of intervening high column density absorbers was found to be minimal. The Nhxigm from the QSO sample scales as (1 + z)^1.5+/-0.2. The mean hydrogen density at z = 0 is n0 = (2.8 +/- 0.3) x 10^-7 cm^-3, the mean IGM temperature over the full redshift range is log(T/K) = 6.5+/-0.1, and the mean metallicity is [X/H] = -1.3+/-0.1(Z 0.05). Aggregating with our previous GRB and blazar tracers, we conclude that we have provided evidence of the IGM contributing substantially and consistently to the total X-ray absorption seen in the spectra. These results are based on the necessarily simplistic slab model used for the IGM, due to the inability of current X-ray data to constrain the IGM redshift distribution

We present a photometric analysis of the TESS light curve of contact binary system DY Cet and the behavior of its orbital period variation. The light curve and published radial velocity data analysis was performed using Wilson-Devinney code. As a result of simultaneous analysis of light curve with radial velocity data, the masses and radii of the system's components were determined as M1=1.55(2)Ms, M2=0.55(1)Ms and R1=1.51(2)Rs, R2=0.95(2)Rs, respectively. The degree of contact (f) and mass ratio (q) of the system were determined as 23% and 0.355(12), respectively. Orbital period analysis of DY Cet was conducted for the first time in this study. It was observed that the orbital period has a sinus-like change with decreasing parabola. To explain the orbital period change, mass transfer between components is proposed with the assumption of conservative mass, and the transfer rate was calculated to be dM/dt=1.1x10^(-7) Ms/yr. A possible third component is suggested for explaining the sinus-like change, and the mass of the unseen component was determined as 0.13 Ms. The age of the DY Cet system was estimated as 3.77 Gyr.

We analyze the white dwarf population in miniJPAS, the first square degree observed with 56 medium-band, 145 A in width optical filters by the Javalambre Physics of the accelerating Universe Astrophysical Survey (J-PAS), to provide a data-based forecast for the white dwarf science with low-resolution (R ~ 50) photo-spectra. We define the sample of the bluest point-like sources in miniJPAS with r < 21.5 mag, point-like probability larger than 0.5, (u-r) < 0.80 mag, and (g-i) < 0.25 mag. This sample comprises 33 sources with spectroscopic information, 11 white dwarfs and 22 QSOs. We estimate the effective temperature (Teff), the surface gravity, and the composition of the white dwarf population by a Bayesian fitting to the observed photo-spectra. The miniJPAS data permit the classification of the observed white dwarfs into H-dominated and He-dominated with 99% confidence, and the detection of calcium absorption and polluting metals down to r ~ 21.5 mag at least for sources with 7000 < Teff < 22000 K, the temperature range covered by the white dwarfs in miniJPAS. The effective temperature is estimated with a 2% uncertainty, close to the 1% from spectroscopy. A precise estimation of the surface gravity depends on the available parallax information. In addition, the white dwarf population at Teff > 7000 K can be segregated from the bluest extragalactic QSOs, providing a clean sample based on optical photometry alone. The J-PAS low-resolution photo-spectra provide precise and accurate effective temperatures and atmospheric compositions for white dwarfs, complementing the data from Gaia. J-PAS will also detect and characterize new white dwarfs beyond the Gaia magnitude limit, providing faint candidates for spectroscopic follow up.

We examine the influence of the quadrupole moment of a slowly rotating neutron star on the oscillations of non-slender accretion tori. We apply previously developed methods to perform analytical calculations of frequencies of the radial epicyclic mode of a torus in the specific case of the Hartle-Thorne geometry. We present here our preliminary results and provide a brief comparison between the calculated frequencies and the frequencies previously obtained assuming both standard and linearized Kerr geometry. Finally, we shortly discuss the consequences for models of high-frequency quasi-periodic oscillations observed in low-mass X-ray binaries.

The correlation between the mass of supermassive black holes (SMBHs; $\mathcal{M}_{\rm BH}$) and their host galaxies ($\mathcal{M}_\star$) in the reionization epoch provides valuable constraints on how their growth are jointly assembled in the early universe. High-redshift quasars typically have a $\mathcal{M}_{\rm BH}$/$\mathcal{M}_\star$ ratio significantly elevated in comparison to the local value. However, it is unclear by how much this apparent offset is driven by observational biases for the most distant quasars. To address this issue, we model the sample selection and measurement biases for a compilation of 20 quasars at $z\sim6$ with host properties based on ALMA observations. We find that the observed distribution of quasars in the $\mathcal{M}_{\rm BH} - \mathcal{M}_\star$ plane could be fully reproduced by assuming that the underlying SMBH populations at $z\sim6$ follow the relationship in the local universe. However, a positive or even a negative evolution in $\mathcal{M}_{\rm BH}$/$\mathcal{M}_\star$ can also explain the data, depending on how the intrinsic scatter evolves and the strength of various systematic uncertainties. This imposes critical requirements on improving the accuracy of mass measurements and expanding the current sample to lower $\mathcal{M}_{\rm BH}$ limits to break degeneracies. Interestingly, quasars that are significant outliers in $\mathcal{M}_{\rm BH}$/$\mathcal{M}_\star$ tend to move towards the local relation given their instantaneous BH accretion rate and star formation rate. This may provide tentative evidence that a self-regulated SMBH-galaxy coevolution scenario is already in place at $z\sim6$, with AGN feedback being a possible driver.

Odd Radio Circles (ORCs) are recently-discovered faint diffuse circles of radio emission, of unknown cause, surrounding galaxies at moderate redshift ($z ~ 0.2-0.6). Here we present detailed new MeerKAT radio images at 1284 MHz of the first ORC, originally discovered with the Australian Square Kilometre Array Pathfinder, with higher resolution (6 arcsec) and sensitivity (~ 2.4 uJy/bm). In addition to the new images, which reveal a complex internal structure consisting of multiple arcs, we also present polarisation and spectral index maps. Based on these new data, we consider potential mechanisms that may generate the ORCs.

Short-lived radioactive isotopes (SLRs) with half-lives between 0.1 to 100 Myr can be used to probe the origin of the Solar System. In this work, we examine the core-collapse supernovae production of the 15 SLRs produced: $^{26}$Al, $^{36}$Cl, $^{41}$Ca, $^{53}$Mn, $^{60}$Fe, $^{92}$Nb, $^{97}$Tc, $^{98}$Tc, $^{107}$Pd, $^{126}$Sn, $^{129}$I, $^{135}$Cs, $^{146}$Sm, $^{182}$Hf, and $^{205}$Pb. We probe the impact of the uncertainties of the core-collapse explosion mechanism by examining a collection of 62 core-collapse models with initial masses of 15, 20, and 25M$_{\odot}$, explosion energies between 3.4$\times$10$^{50}$ and 1.8$\times$10$^{52}$ ergs and compact remnant masses between 1.5M$_{\odot}$and 4.89M$_{\odot}$. We identify the impact of both explosion energy and remnant mass on the final yields of the SLRs. Isotopes produced within the innermost regions of the star, such as $^{92}$Nb and $^{97}$Tc, are the most affected by the remnant mass, $^{92}$Nb varying by five orders of magnitude. Isotopes synthesised primarily in explosive C-burning and explosive He-burning, such as $^{60}$Fe, are most affected by explosion energies. $^{60}$Fe increases by two orders of magnitude from the lowest to the highest explosion energy in the 15M$_{\odot}$model. The final yield of each examined SLR is used to compare to literature models.

We present an overview of the "FAIR Guiding Principles for scientific data management and stewardship", first published in 2016, and how they relate to astronomical data management. In particular, we discuss the connection between the FAIR principles and IVOA standards, and how data management systems with these standards implemented are close to compliance. We then look at what extra components are required to make astronomical data FAIR. Finally, we give a case study of the All-Sky Virtual Observatory (Australia's node of the VO) and their implementation of the FAIR principles.

Computational heliophysics has shed light on the fundamental physical processes inside the Sun, such as the differential rotation, meridional circulation, and dynamo-generation of magnetic fields. However, despite the substantial advances, the current results of 3D MHD simulations are still far from reproducing helioseismic inferences and surface observations. The reason is the multi-scale nature of the solar dynamics, covering a vast range of scales, which cannot be solved with the current computational resources. In such a situation, significant progress has been achieved by the mean-field approach, based on the separation of small-scale turbulence and large-scale dynamics. The mean-field simulations can reproduce solar observations, qualitatively and quantitatively, and uncover new phenomena. However, they do not reveal the complex physics of large-scale convection, solar magnetic cycles, and the magnetic self-organization that causes sunspots and solar eruptions. Thus, developing a synergy of these approaches seems to be a necessary but very challenging task.

A birefringent universe could show itself through a rotation of the plane of polarisation of the cosmic microwave background photons. This is usually investigated using polarisation $B$ modes, which is degenerate with miscalibration of the orientation of the polarimeters. Here we point out an independent method for extracting the birefringence angle using only temperature and $E$-mode signals. We forecast that, with an ideal cosmic-variance-limited experiment, we could constrain a birefringence angle of $0.3^\circ$ with $3\,\sigma$ statistical significance, which is close to the current constraints using $B$ modes. We explore how this method is affected by the systematic errors introduced by the polarisation efficiency. In the future, this could provide an additional way of checking any claimed $B$-mode derived birefringence signature.

Rings and asymmetries in protoplanetary disks are considered as signposts of ongoing planet formation. In this work, we conduct three-dimensional radiative transfer simulations to model the intriguing disk around HD 143006 that has three dust rings and a bright arc. A complex geometric configuration, with a misaligned inner disk, is assumed to account for the asymmetric structures. The two-dimensional surface density is constructed by iteratively fitting the ALMA data. We find that the dust temperature displays a notable discontinuity at the boundary of the misalignment. The ring masses range from 0.6 to 16Mearth that are systematically lower than those inferred in the younger HL Tau disk. The arc occupies nearly 20% of the total dust mass. Such a high mass fraction of dust grains concentrated in a local region is consistent with the mechanism of dust trapping into vortices. Assuming a gas-to-dust mass ratio of 30 that is constant throughout the disk, the dense and cold arc is close to the threshold of being gravitationally unstable, with the Toomre parameter Q~1.3. Nevertheless, our estimate of Q relies on the assumption for the unknown gas-to-dust mass ratio. Adopting a lower gas-to-dust mass ratio would increase the inferred Q value. Follow-up high resolution observations of dust and gas lines are needed to clarify the origin of the substructures.

We investigate the beaming of 11 Io-Jupiter decametric (Io-DAM) emissions observed by Juno/Waves, the Nan\c cay Decameter Array and NenuFAR. Using an up-to-date magnetic field model and three methods to position the active Io Flux Tube (IFT), we accurately locate the radiosources and determine their emission angle $\theta$ from the local magnetic field vector. These methods use (i) updated models of the IFT equatorial lead angle, (ii) ultraviolet (UV) images of Jupiter's aurorae and (iii) multi-point radio measurements. The kinetic energy $E_{e-}$ of source electrons is then inferred from $\theta$ in the framework of the Cyclotron Maser Instability. The precise position of the active IFT achieved from methods (ii,iii) can be used to test the effective torus plasma density. Simultaneous radio/UV observations reveal that multiple Io-DAM arcs are associated with multiple UV spots and provide the first direct evidence of an Io-DAM arc associated with a trans-hemispheric beam UV spot. Multi-point radio observations probe the Io-DAM sources at various altitudes, times and hemispheres. Overall, $\theta$ varies a function of frequency (altitude), by decreasing from $75^\circ-80^\circ$ to $70^\circ-75^\circ$ over $10-40$ MHz with slightly larger values in the northern hemisphere, and independently varies as a function of time (or longitude of Io). Its uncertainty of a few degrees is dominated by the error on the longitude of the active IFT. The inferred values of $E_{e-}$ also vary as a function of altitude and time. For the 11 investigated cases, they range from 3 to 16 keV, with a $6.6\pm2.7$ keV average.

We use updated Hubble parameter and baryon acoustic oscillation data, as well as other lower-redshift Type Ia supernova, Mg II reverberation-measured quasar, quasar angular size, H II starburst galaxy, and Amati-correlated gamma-ray burst data, to jointly constrain cosmological parameters in six cosmological models. The joint analysis provides model-independent determinations of the Hubble constant, $H_0=69.7\pm1.2$ $\rm{km \ s^{-1} \ Mpc^{-1}}$, and the current non-relativistic matter density parameter, $\Omega_{m0}=0.295\pm0.017$. These error bars are factors of 2.2 and 2.3 larger than the corresponding error bars in the flat $\Lambda$CDM model from Planck TT,TE,EE+lowE+lensing cosmic microwave background anisotropy data. Based on the deviance information criterion (DIC), the flat $\Lambda$CDM model is most favored but mild dark energy dynamics and a little spatial curvature are not ruled out.

Massive stars moving at supersonic peculiar velocities through the interstellar medium (ISM) can create bow shocks, arc-like structures at the interface between the stellar wind and the ISM. Many such bow shocks have been detected and catalogued at IR wavelengths, but detections in other wavebands remain rare. Strikingly, while electrons are expected to be accelerated in the bow shock and their non-thermal emission may include synchrotron emission at low frequencies, only two massive runaway stellar bow shocks have to date been detected in the radio band. Here, we examine a sample of fifty IR-detected bow shocks from the E-BOSS catalogues in recently released radio images from the Rapid ASKAP Continuum Survey (RACS). We identify three confident and three likely counterparts, as well as three inconclusive candidates requiring confirmation via follow-up observations. These detections significantly increase the number of known radio massive stellar bow shocks and highlight the advantage of dedicated searches with current and next-generation radio telescopes. We investigate the underlying radio emission mechanism for these radio sources, finding a mix of free-free-dominated and synchrotron-dominated systems. We also discuss the non-detected targets by putting constraints on their emission properties and investigating their detectability with future observations. Finally, we propose several future avenues of research to advance the study and understanding of bow shocks at radio frequencies.

Extremely Metal-Poor (EMP) stars provide a valuable probe of early chemical enrichment in the Milky Way. Here we leverage a large sample of $\sim600,000$ high-resolution stellar spectra from the GALAH survey plus a machine learning algorithm to find 54 candidates with estimated [Fe/H]~$\leq$~-3.0, 6 of which have [Fe/H]~$\leq$~-3.5. Our sample includes $\sim 20 \%$ main sequence EMP candidates, unusually high for \emp surveys. We find the magnitude-limited metallicity distribution function of our sample is consistent with previous work that used more complex selection criteria. The method we present has significant potential for application to the next generation of massive stellar spectroscopic surveys, which will expand the available spectroscopic data well into the millions of stars.

X-rays from High-Mass X-ray Binaries (HMXBs) are likely the main source of heating of the intergalactic medium (IGM) during Cosmic Dawn (CD), before the completion of reionization. This Epoch of Heating (EoH; $z\sim 10-15$) should soon be detected via the redshifted 21-cm line from neutral hydrogen, allowing us to indirectly study the properties of HMXBs in the unseen, first galaxies. Low-redshift observations, as well as theoretical models, imply that the integrated X-ray luminosity to star formation rate of HMXBs ($L_{\rm X}/{\rm SFR}$) should increase in metal-poor environments, typical of early galaxies. Here we study the impact of the metallicity ($Z$) dependence of $L_{\rm X}/{\rm SFR}$ during the EoH. For our fiducial models, galaxies with star formation rates of order $10^{-3} - 10^{-1}$ $M_\odot$ yr$^{-1}$ and metallicities of order $10^{-3} - 10^{-2}$ $Z_\odot$ are the dominant contributors to the X-ray background (XRB) during this period. Different $L_{\rm X}/{\rm SFR}$-$Z$ relations result in factors of $\sim$ 3 differences in these ranges, as well as in the mean IGM temperature and the large-scale 21-cm power, at a given redshift. We compute mock 21-cm observations adopting as a baseline a 1000h integration with the upcoming Square Kilometer Array (SKA), for two different $L_{\rm X}/{\rm SFR}$-$Z$ relations. We perform inference on these mock observations using the common simplification of a constant $L_{\rm X}/{\rm SFR}$, finding that constant $L_{\rm X}/{\rm SFR}$ models can recover the IGM evolution of the more complicated $L_{\rm X}/{\rm SFR}$-$Z$ simulations only during the EoH. At $z<10$, where the typical galaxies are more polluted, constant $L_{\rm X}/{\rm SFR}$ models over-predict the XRB and its relative contribution to the early stages of the reionization.

Context: Periodicity in 6.7 GHz methanol maser sources is a rare phenomenon that was discovered during long-term monitoring programmes. Understanding the underlying processes that lead to periodic variability might provide insights into the physical processes in high-mass star-forming regions. Aims: We aim to identify and describe new periodic methanol masers. Methods: The observations were obtained with the Torun 32m antenna. Time series analysis was conducted using well-proven statistical methods. Additionally, NEOWISE data were used to search for a correlation between infrared and maser fluxes. Results: We found two new periodic sources, G45.804-0.356 and G49.043-1.079, with periods of 416.9 and 469.3 days, respectively. For G49.043-1.079, infrared variability is simultaneous with methanol flares. Conclusions: A most likely cause of the periodicity in G49.043-1.079 is modulated accretion. For G45.804-0.356, the periodicity cannot be explained with the available data, and further research is needed.

Remarkable development of cosmology is benefited from the increasingly improved measurements of cosmic distances including absolute distances and relative distances. In recent years, however, the emerged cosmological tensions motivate us to explore the independent and precise late-universe probes. The two observational effects of strong gravitational lensing (SGL), the velocity dispersions of lens galaxies and the time delays between multiple images, can provide measurements of relative and absolute distances respectively, and their combination is possible to break the degeneracies between cosmological parameters and enable tight constraints on cosmological parameters. In this paper, we combine the observed 130 SGL systems with velocity-dispersion measurements and 7 SGL systems with time-delay measurements to constrain dark-energy cosmological models. It is found that the combination of the two effects does not significantly break the degeneracies between cosmological parameters as expected. However, with the simulations of 8000 SGL systems with well-measured velocity dispersions and 55 SGL systems with well-measured time delays based on the forthcoming LSST survey, we find that the combination of two effects can significantly break the parameter degeneracies, and make the constraint precision of cosmological parameters meet the standard of precision cosmology. We conclude that the observations of SGL will become a useful late-universe probe for precisely measuring cosmological parameters.

Isotopic ratios provide a powerful tool for understanding the origins of materials, including the volatile and refractory matter within solar system bodies. Recent high sensitivity observations of molecular isotopologues, in particular with ALMA, have brought us new information on isotopic ratios of hydrogen, carbon, nitrogen and oxygen in star and planet forming regions as well as the solar system objects. Solar system exploration missions, such as Rosetta and Cassini, have given us further new insights. Meanwhile, the recent development of sophisticated models for isotope chemistry including detailed gas-phase and grain surface reaction network has made it possible to discuss how isotope fractionation in star and planet forming regions is imprinted into the icy mantles of dust grains, preserving a record of the initial isotopic state of solar system materials. This chapter reviews recent progress in observations of molecular isotopologues in extra-solar planet forming regions, prestellar/protostellar cores and protoplanetary disks, as well as objects in our solar system -- comets, meteorites, and planetary/satellite atmospheres -- and discusses their connection by means of isotope chemical models.

The physical properties responsible for the formation and evolution of the corona and heliosphere are still not completely understood. 3D MHD global modeling is a powerful tool to investigate all the possible candidate processes. To fully understand the role of each of them, we need a validation process where the output from the simulations is quantitatively compared to the observational data. In this work, we present the results from our validation process applied to the wave turbulence driven 3D MHD corona-wind model WindPredict-AW. At this stage of the model development, we focus the work to the coronal regime in quiescent condition. We analyze three simulations results, which differ by the boundary values. We use the 3D distributions of density and temperature, output from the simulations at the time of around the First Parker Solar Probe perihelion (during minimum of the solar activity), to synthesize both extreme ultraviolet (EUV) and white light polarized (WL pB) images to reproduce the observed solar corona. For these tests, we selected AIA 193 A, 211 A and 171 A EUV emissions, MLSO K-Cor and LASCO C2 pB images obtained the 6 and 7 November 2018. We then make quantitative comparisons of the disk and off limb corona. We show that our model is able to produce synthetic images comparable to those of the observed corona.

We analyse measurements of the evolving stellar mass (M0) at which the bending of the star-forming main sequence (MS) occurs over 0<z<4. We find M0~10^{10}Msun over 0<z<1, then M0 rises up to ~10^{11}Msun at z=2, and then stays flat or slowly increases towards higher redshifts. When converting M0 values into hosting dark matter halo masses, we show that this behaviour is remarkably consistent with the evolving cold- to hot-accretion transition mass, as predicted by theory and defined by the redshift-independent Mshock at z<1.4 and by the rising Mstream at z>1.4 (for which we propose a revision in agreement with latest simulations). We hence argue that the MS bending is primarily due to the lessening of cold-accretion causing a reduction in available cold gas in galaxies and supports predictions of gas feeding theory. In particular, the rapidly rising M0 with redshift at z>1 is confirming evidence for the cold-streams scenario. In this picture, a progressive fueling reduction rather than its sudden suppression in halos more massive than Mshock/Mstream produces a nearly constant star-formation rate in galaxies with stellar masses larger than M0, and not their quenching, for which other physical processes are thus required. Compared to the knee M* in the stellar mass function of galaxies, M0 is significantly lower at z<1.5, and higher at z>2, suggesting that the imprint of gas deprivation on the distribution of galaxy masses happened at early times (z>1.5-2). The typical mass at which galaxies inside the MS become bulge-dominated evolves differently from M0, consistent with the idea that bulge-formation is a distinct process from the phasing-out of cold-accretion.

The photometric transit method has been the most effective method to detect and characterize exoplanets as several ground-based as well as space-based survey missions have discovered thousands of exoplanets using this method. With the advent of the upcoming next generation large telescopes, the detection of exomoons in a few of these exoplanetary systems is very plausible. In this paper, we present a comprehensive analytical formalism in order to model the transit light curves for such moon hosting exoplanets. In order to achieve analytical formalism, we have considered circular orbit of the exomoon around the host planet, which is indeed the case for tidally locked moons. The formalisms uses the radius and orbital properties of both the host planet and its moon as model parameters. The co-alignement or non-coalignment of the orbits of the planet and the moon is parametrized by two angles of separation and thus can be used to model all the possible orbital alignments for a star-planet-moon system. This formalism also provides unique and direct solutions to every possible star-planet-moon three circular body alignments. Using the formula derived, a few representative light curves are also presented.

We use the COSMOS2020 catalogue to measure the stellar-to-halo mass relation (SHMR) divided by central and satellite galaxies from $z=0.2$ to $z = 5.5$. Starting from accurate photometric redshifts we measure the near-infrared selected two-point angular correlation and stellar mass functions in ten redshift bins and fit them with an HOD-based model. At each redshift, we measure the ratio of stellar mass to halo mass, $M_*/M_h$, which shows the characteristic strong dependence of halo mass with a peak at $M_h^{\rm peak} \sim 2$. Our results are in accordance with the scenario in which the peak of star-formation efficiency moves towards more massive halos at higher redshifts. We also measure the fraction of satellites as a function of stellar mass and redshift. For all stellar mass thresholds the satellite fraction decreases at higher redshifts. At a given redshift there is a higher fraction of low-mass satellites. The satellite contribution to the total stellar mass budget in halos becomes more important than centrals at halo masses of about $M_h > 10^{13} \, M_{\odot}$ and always stays below by peak, indicating that quenching mechanisms are present in massive halos that keep the star-formation efficiency low. Finally, we compare our results with three hydrodynamical simulations Horizon-AGN, Illustris-TNG-100 and EAGLE. We find that the most significant discrepancy is at the high mass end, where the simulations generally show that satellites have a higher contribution to the total stellar mass budget than the observations. This, together with the finding that the fraction of satellites is higher in the simulations, indicates that the feedback mechanisms acting in group-and cluster-scale halos appear to be less efficient in quenching the mass assembly of satellites, and/or that quenching occurs much later in the simulations.

The intuition suggested by the Drake equation implies that technology should be less prevalent than biology in the galaxy. However, it has been appreciated for decades in the SETI community that technosignatures could be more abundant, longer-lived, more detectable, and less ambiguous than biosignatures. We collect the arguments for and against technosignatures' ubiquity and discuss the implications of some properties of technological life that fundamentally differ from nontechnological life in the context of modern astrobiology: It can spread among the stars to many sites, it can be more easily detected at large distances, and it can produce signs that are unambiguously technological. As an illustration in terms of the Drake equation, we consider two Drake-like equations, for technosignatures (calculating N(tech)) and biosignatures (calculating N(bio)). We argue that Earth and humanity may be poor guides to the longevity term L and that its maximum value could be very large, in that technology can outlive its creators and even its host star. We conclude that while the Drake equation implies that N(bio)>>N(tech), it is also plausible that N(tech)>>N(bio). As a consequence, as we seek possible indicators of extraterrestrial life, for instance, via characterization of the atmospheres of habitable exoplanets, we should search for both biosignatures and technosignatures. This exercise also illustrates ways in which biosignature and technosignature searches can complement and supplement each other and how methods of technosignature search, including old ideas from SETI, can inform the search for biosignatures and life generally.

Magnetic interactions between stars and close-in planets may lead to a detectable signal on the stellar disk. HD 189733 is one of the key exosystems thought to harbor magnetic interactions, which may have been detected in August 2013. We present a set of twelve wind models at that period, covering the possible coronal states and coronal topologies of HD 189733 at that time. We assess the power available for the magnetic interaction and predict its temporal modulation. By comparing the predicted signal with the observed signal, we find that some models could be compatible with an interpretation based on star-planet magnetic interactions. We also find that the observed signal can be explained only with a stretch-and-break interaction mechanism, while that the Alfv\'en wings scenario cannot deliver enough power. We finally demonstrate that the past observational cadence of HD 189733 leads to a detection rate of only between 12 to 23%, which could explain why star-planet interactions have been hard to detect in past campaigns. We conclude that the firm confirmation of their detection will require dedicated spectroscopic observations covering densely the orbital and rotation period, combined with scarcer spectropolarimetric observations to assess the concomitant large-scale magnetic topology of the star.

After NANOGrav, the IPTA collaboration also reports a strong evidence of a stochastic gravitation wave background. This hint has very important implications for fundamental physics. With the recent IPTA data release two, we attempt to search signals of light new physics. and give new constraints on the audible axion, domain walls and cosmic strings models. We find that the best fit point corresponding to a decay constant $F\approx5\times10^{17}$ GeV and an axion mass $m_a\approx2\times10^{-13}$ eV from NANOGrav data is ruled out by IPTA at beyond $2\sigma$ confidence level. Fixing the coupling strength $\lambda=1$, we obtain a $2\sigma$ lower bound on the breaking scale of $Z_2$ symmetry $\eta>135$ TeV. Interestingly, we give a very strong restriction on the cosmic-string tension $\mathrm{log}_{10}\,G\mu=-8.93_{-0.06}^{+0.12}$ at $1\sigma$ confidence level. Employing the rule of Bayes factor, we find that IPTA data has a moderate, strong and inconclusive preference of an uncorrelated common power-law (CPL) model over audible axion, domain walls and cosmic strings, respectively. This means that it is hard to distinguish CPL from cosmic strings with current observations and more pulsar timing data with high precision are required to give new clues of underlying physics.

The Sino-French space mission SVOM is mainly designed to detect, localize and follow-up Gamma-Ray Bursts and other high-energy transients. The satellite, to be launched mid 2023, embarks two wide-field gamma-ray instruments and two narrow-field telescopes operating at X-ray and optical wavelengths. It is complemented by a dedicated ground segment encompassing a set of wide-field optical cameras and two 1-meter class follow-up telescopes. In this contribution, we describe the main characteristics of the mission and discuss its scientific rationale and some original GRB studies that it will enable.

Taking advantage of the now available Gaia EDR3 parallaxes, we carry out an archival {\it Hubble Space Telescope} (HST) far ultraviolet spectroscopic analysis of 10 cataclysmic variable systems, including 5 carefully selected eclipsing systems. We obtain accurate white dwarf (WD) masses and temperatures, in excellent agreement with the masses for 4 of the eclipsing systems. For three systems in our sample, BD Pav, HS 2214, and TT Crt, we report the first robust masses for their WDs. We modeled the absorption lines to derive the WD chemical abundances and rotational velocities for each of the ten systems. As expected, for five higher inclination ($i \gtrsim 75^{\circ}$) systems, the model fits are improved with the inclusion of a cold absorbing slab (an iron curtain masking the WD) with $N_{\rm H} \approx 10^{20}-10^{22}$cm$^{-2}$. Modeling of the metal lines in the HST spectra reveals that 7 of the 10 systems have significant subsolar carbon abundance, and six have subsolar silicon abundance, thereby providing further evidence that CV WDs exhibit subsolar abundances of carbon and silicon. We suggest that strong aluminum absorption lines (and iron absorption features) in the spectra of some CV WDs (such as IR Com) may be due to the presence of a {\it thin} iron curtain ($N_{\rm H}\approx 10^{19}$cm$^{-2}$) rather than to suprasolar aluminum and iron abundances in the WD photosphere. The derived WD (projected) rotational velocities all fall in the range $\approx 100-400$~km/s, all sub-Keplerian similar to the values obtained in earlier studies.

Type II radio bursts are thought to be produced by shock waves in the solar atmosphere. However, what magnetic conditions are needed for the generation of type II radio bursts is still a puzzling issue. Here, we quantify the magnetic structure of a coronal shock associated with a type II radio burst. Based on the multi-perspective extreme-ultraviolet observations, we reconstruct the three-dimensional (3D) shock surface. By using a magnetic field extrapolation model, we then derive the orientation of the magnetic field relative to the normal of the shock front ($\theta_{\rm Bn}$) and Alfv\'{e}n Mach number ($M_A$) on the shock front. Combining the radio observations from Nancay Radio Heliograph, we obtain the source region of the type II radio burst on the shock front. It is found that the radio burst is generated by a shock with $M_A \gtrsim 1.5$ and a bimodal distribution of $\theta_{Bn}$. We also use the Rankine-Hugoniot relations to quantify the properties of the shock downstream. Our results provide a quantitative 3D magnetic structure condition of a coronal shock that produces a type II radio burst.

We report the first detection of large-scale matter flows around cosmic voids at a median redshift z = 2.49. Voids are identified within a tomographic map of the large-scale matter density built from eBOSS Lyman-$\alpha$ (Lya) forests in SDSS Stripe 82. We measure the imprint of flows around voids, known as redshift-space distortions (RSD), with a statistical significance of 10 $\sigma$. The observed quadrupole of the void-forest cross-correlation is described by a linear RSD model. The derived RSD parameter is $\beta = 0.52 \pm 0.05$. Our model accounts for the tomographic effect induced by the Lya data being located along parallel quasar lines of sight. This work paves the way towards growth-rate measurements at redshifts currently inaccessible to galaxy surveys.

The Tibet AS$\gamma$ collaboration has reported a diffuse $\gamma$-ray emission signal from the Galactic Plane. We consider the possibility that the diffuse emission from the outer Galactic Plane at the highest energies is produced by cosmic rays spreading from a single supernova-type source either in the Local or Perseus arm of the Milky Way. We show that anisotropic diffusion of multi-PeV cosmic rays along the Galactic magnetic field can produce an extended source spanning ten(s) of degrees on the sky, with a flux-per-unit-solid-angle consistent with Tibet AS$\gamma$ measurements. Observations of this new type of very extended sources, and measurements of their morphologies, can be used to characterize the anisotropic diffusion of PeV cosmic rays in the Galactic magnetic field, and to constrain the locations and properties of past PeVatrons.

Following the detection of a long GRB 190919B by INTEGRAL (INTErnational Gamma-Ray Astrophysics Laboratory), we obtained an optical photometric sequence of its optical counterpart. The light curve of the optical emission exhibits an unusually steep rise ~100 s after the initial trigger. This behaviour is not expected from a 'canonical' GRB optical afterglow. As an explanation, we propose a scenario consisting of two superimposed flares: an optical flare originating from the inner engine activity followed by the hydrodynamic peak of an external shock. The inner-engine nature of the first pulse is supported by a marginal detection of flux in hard X-rays. The second pulse eventually concludes in a slow constant decay, which, as we show, follows the closure relations for a slow cooling plasma expanding into the constant interstellar medium and can be seen as an optical afterglow sensu stricto.

SubHalo Abundance Matching (SHAM) is an empirical method for constructing galaxy catalogues based on high-resolution $N$-body simulations. We apply SHAM on the UNIT simulation to simulate SDSS BOSS/eBOSS Luminous Red Galaxies (LRGs) within a wide redshift range of $0.2 < z < 1.0$. Besides the typical SHAM scatter parameter $\sigma$, we include $v_{\rm smear}$ and $V_{\rm ceil}$ to take into account the redshift uncertainty and the galaxy incompleteness respectively. These two additional parameters are critical for reproducing the observed 2PCF multipoles on 5--25$\,h^{-1}\,{\rm Mpc}$. The redshift uncertainties obtained from the best-fitting $v_{\rm smear}$ agree with those measured from repeat observations for all SDSS LRGs except for the LOWZ sample. We explore several potential systematics but none of them can explain the discrepancy found in LOWZ. Our explanation is that the LOWZ galaxies might contain another type of galaxies which needs to be treated differently. The evolution of the measured $\sigma$ and $V_{\rm ceil}$ also reveals that the incompleteness of eBOSS galaxies decreases with the redshift. This is the consequence of the magnitude lower limit applied in eBOSS LRG target selection. Our SHAM also set upper limits for the intrinsic scatter of the galaxy--halo relation given a complete galaxy sample: $\sigma_{\rm int}<0.31$ for LOWZ at $0.2<z<0.33$, $\sigma_{\rm int}<0.36$ for LOWZ at $0.33<z<0.43$, and $\sigma_{\rm int}<0.46$ for CMASS at $0.43<z<0.51$. The projected 2PCFs of our SHAM galaxies also agree with the observational ones on the 2PCF fitting range.

The TESS space mission provides us with high-precision photometric observations of bright stars over more than 70% of the entire sky, allowing us to revisit and characterise well-known stars. We aim to conduct an asteroseismic analysis of the gamma Doradus star HD112429 using both the available ground-based spectroscopy and TESS photometry, and assess the conditions required to measure the near-core rotation rate and buoyancy travel time. We collect and reduce the available five sectors of short-cadence TESS photometry of this star, as well as 672 legacy observations from six medium- to high-resolution ground-based spectrographs. We determine the stellar pulsation frequencies from both data sets using iterative prewhitening, do asymptotic g mode modelling of the star and investigate the corresponding spectral line profile variations using the pixel-by-pixel method. We validate the pulsation frequencies from the TESS data up to $S/N \geq 5.6$, confirming recent reports in the literature that the classical criterion $S/N \geq 4$ does not suffice for space-based observations. We identify the pulsations as prograde dipole g modes and r-mode pulsations, and measure a near-core rotation rate of $1.536(3) d^{-1}$ and a buoyancy travel time $\Pi_0$ of 4190(50) s. These results are in agreement with the observed spectral line profile variations, which were qualitatively evaluated using a newly developed toy model. We establish a set of conditions that have to be fulfilled for an asymptotic asteroseismic analysis of g-mode pulsators. In the case of HD112429, two TESS sectors of space photometry suffice. Although a detailed asteroseismic modelling analysis is not viable for g-mode pulsators with only short or sparse light curves of space photometry, we find that it is possible to determine global asteroseismic quantities for a subset of these stars. (abbreviated.)

The massive outburst of the comet 17P/Holmes in October 2007 is the largest known outburst by a comet thus far. We present a new comprehensive model describing the evolution of the dust trails produced in this phenomenon. The model comprises of multiparticle Monte Carlo approach including the solar radiation pressure effects, gravitational disturbance caused by Venus, Earth and Moon, Mars, Jupiter and Saturn, and gravitational interaction of the dust particles with the parent comet itself. Good accuracy of computations is achieved by its implementation in Orekit, which executes Dormad-Prince numerical integration methods with higher precision. We demonstrate performance of the model by simulating particle populations with sizes from 0.001 mm to 1 mm with corresponding spherically symmetric ejection speed distribution, and towards the Sun outburst modelling. The model is supplemented with and validated against the observations of the dust trail in common nodes for 0.5 and 1 revolutions. In all cases, the predicted trail position showed a good match to the observations. Additionally, the hourglass pattern of the trail was observed for the first time within this work. By using variations of the outburst model in our simulations, we determine that the assumption of the spherical symmetry of the ejected particles leads to the scenario compatible with the observed hourglass pattern. Using these data, we make predictions for the two-revolution dust trail behavior near the outburst point that should be detectable by using ground-based telescopes in 2022.

We have used data from the UKIRT Hemisphere Survey (UHS) to search for substellar members of the Hyades cluster. Our search recovered several known substellar Hyades members, and two known brown dwarfs that we suggest may be members based on a new kinematic analysis. We uncovered thirteen new substellar Hyades candidates, and obtained near-infrared follow-up spectroscopy of each with IRTF/SpeX. Six candidates with spectral types between M7 and L0 are ruled out as potential members based on their photometric distances ($\gtrsim$100 pc). The remaining seven candidates, with spectral types between L5 and T4, are all potential Hyades members, with five showing strong membership probabilities based on BANYAN $\Sigma$ and a convergent point analysis. Distances and radial velocities are still needed to confirm Hyades membership. If confirmed, these would be some of the lowest mass free-floating members of the Hyades yet known, with masses as low as $\sim$30 $M_{\rm Jup}$. An analysis of all known substellar Hyades candidates shows evidence that the full extent of the Hyades has yet to be probed for low-mass members, and more would likely be recovered with deeper photometric and astrometric investigations.

The light-cone effect breaks the periodicity and statistical homogeneity (ergodicity) along the line-of-sight direction of cosmological emission/absorption line surveys. The spherically averaged power spectrum (SAPS), which by definition assumes ergodicity and periodicity in all directions, can only quantify some of the second-order statistical information in the 3D light-cone signals and therefore gives a biased estimate of the true statistics. The multi-frequency angular power spectrum (MAPS), by extracting more information from the data, does not rely on these assumptions. It is therefore better aligned with the properties of the signal. We have compared the performance of the MAPS and SAPS metrics for parameter estimation for a mock 3D light-cone observation of the 21-cm signal from the Epoch of Reionization. Our investigation is based on a simplified 3-parameter 21cmFAST model. We find that the MAPS produces parameter constraints which are a factor of $\sim 2$ more stringent than when the SAPS is used. The significance of this result does not change much even in the presence of instrumental noise expected for 128 hours of SKA-Low observations. Our results therefore suggest that a parameter estimation framework based on the MAPS metric would yield superior results over one using the SAPS metric.

We compare and validate COLA (COmoving Lagrangian Acceleration) simulations against existing emulators in the literature, namely Bacco and Euclid Emulator 2. Our analysis focuses on the non-linear response function, i.e., the ratio between the non-linear dark matter power spectrum in a given cosmology with respect to a pre-defined reference cosmology, which is chosen to be the Euclid Emulator 2 reference cosmology in this paper. We vary three cosmological parameters, the total matter density, the amplitude of the primordial scalar perturbations and the spectral index. By comparing the COLA non-linear response function with those computed from each emulator in the redshift range $0 \leq z \leq 3$, we find that the COLA method is in excellent agreement with the two emulators for scales up to $k \sim 1 \ h$/Mpc as long as the deviations of the matter power spectrum from the reference cosmology are not too large. We validate the implementation of massive neutrinos in our COLA simulations by varying the sum of neutrino masses to three different values, $0.0$ eV, $0.058$ eV and $0.15$ eV. We show that all three non-linear prescriptions used in this work agree at the $1\%$ level at $k \leq 1 \ h$/Mpc. We then introduce the Effective Field Theory of Dark Energy in our COLA simulations using the $N$-body gauge method. We consider two different modified gravity models in which the growth of structure is enhanced or suppressed at small scales, and show that the response function with respect to the change of modified gravity parameters depends weakly on cosmological parameters in these models.

The next generation of cosmological surveys will have unprecedented measurement precision, hence they hold the power to put theoretical ideas to the most stringent tests yet. However, in order to realise the full potential of these measurements, we need to ensure that we apply the most effective analytical tools. We need to identify which cosmological observables are the best cosmological probes. Two commonly used cosmological observables are the galaxy number counts, and the cosmic magnification. The galaxy number counts are a diagnostic of the growth of structure, hence are able to probe the nature of dark energy (DE). On the other hand, the cosmic magnification is a diagnostic of cosmic distances and sizes -- consequently, a diagnostic of the geometry of the large scale structure -- hence is also able to probe the nature of DE. Both of these observables have been investigated extensively in cosmological analyses, but only separately. In this paper, by using the angular power spectrum, we investigate both the galaxy number counts and the cosmic magnification, on ultra-large scales. Our results suggest that measuring relativistic effects in the cosmic magnification angular power spectrum will be relatively better than in the number-count angular power spectrum, at all redshifts z. On the other hand, we found that without relativistic effects, the number-count angular power spectrum will be relatively better in probing the imprint of interacting DE (IDE), at all z. At low z (up to around z = 0.1), relativistic effects enables the cosmic magnification angular power spectrum to be a relatively better probe of the IDE imprint; while at higher z (up to z < 3), the number-count angular power spectrum becomes the better probe of IDE imprint. However, at z = 3 and higher, our results suggest that either the number-count or the magnification angular power spectrum, will suffice to probe IDE.

The solar nebula contained a number of short-lived radionuclides (SLRs) with half-lives of tens of Myr or less, comparable to the timescales for formation of protostars and protoplanetary disks. Therefore, determining the origins of SLRs would provide insights into star formation and the Sun's astrophysical birth environment. In this chapter, we review how isotopic studies of meteorites reveal the existence and abundances of these now-extinct radionuclides; and the evidence that the SLR ${}^{10}{\rm Be}$, which uniquely among the SLRs is not produced during typical stellar nucleosynthesis, was distributed homogeneously in the solar nebula. We review the evidence that the SLRs ${}^{26}{\rm Al}$, ${}^{53}{\rm Mn}$, and ${}^{182}{\rm Hf}$, and other radionuclides, were also homogeneously distributed and can be used to date events during the Solar System's planet-forming epoch. The homogeneity of the SLRs, especially ${}^{10}{\rm Be}$, strongly suggests they were all inherited from the Sun's molecular cloud, and that production by irradiation within the solar nebula was very limited, except for ${}^{36}{\rm Cl}$. We review astrophysical models for the origin of ${}^{10}{\rm Be}$, showing that it requires that the Sun formed in a spiral arm of the Galaxy with higher star formation rate than the Galaxy-wide average. Likewise, we review the astrophysical models for the origins of the other SLRs and show that they likely arose from contamination of the Sun's molecular cloud by massive stars over tens of Myr, most likely dominated by ejecta from Wolf-Rayet stars. The other SLRs also demand formation of the Sun in a spiral arm of the Galaxy with a star formation rate as high as demanded by the Solar System initial ${}^{10}{\rm Be}$ abundance. We discuss the astrophysical implications, and suggest further tests of these models and future directions for the field.

In this second installment of SETI in 20xx, we very briefly and subjectively review developments in SETI in 2021. Our primary focus is 93 papers and books published or made public in 2021, which we sort into six broad categories: results from actual searches, new search methods and instrumentation, target and frequency selection, the development of technosignatures, theory of ETIs, and social aspects of SETI.

In this chapter we review recent advances in understanding the roles that magnetic fields play throughout the star formation process, gained through observations and simulations of molecular clouds, the dense, star-forming phase of the magnetised, turbulent interstellar medium (ISM). Recent results broadly support a picture in which the magnetic fields of molecular clouds transition from being gravitationally sub-critical and near equipartition with turbulence in low-density cloud envelopes, to being energetically sub-dominant in dense, gravitationally unstable star-forming cores. Magnetic fields appear to play an important role in the formation of cloud substructure by setting preferred directions for large-scale gas flows in molecular clouds, and can direct the accretion of material onto star-forming filaments and hubs. Low-mass star formation may proceed in environments close to magnetic criticality; high-mass star formation remains less well-understood, but may proceed in more supercritical environments. The interaction between magnetic fields and (proto)stellar feedback may be particularly important in setting star formation efficiency. We also review a range of widely-used techniques for quantifying the dynamic importance of magnetic fields, concluding that better-calibrated diagnostics are required in order to use the spectacular range of forthcoming observations and simulations to quantify our emerging understanding of how magnetic fields influence the outcome of the star formation process.

We report on Fe I in the day-side atmosphere of the ultra-hot Jupiter WASP-33b, providing evidence for a thermal inversion in the presence of an atomic species. We also introduce a new way to constrain the planet's brightness variation throughout its orbit, including its day-night contrast and peak phase offset, using high-resolution Doppler spectroscopy alone. We do so by analyzing high-resolution optical spectra of six arcs of the planet's phase curve, using ESPaDOnS on the Canada-France-Hawaii telescope and HDS on the Subaru telescope. By employing a likelihood mapping technique, we explore the marginalized distributions of parameterized atmospheric models, and detect Fe I emission at high significance ($>10.4\sigma$) in our combined data sets, located at $K_{\rm p}=222.1\pm0.4$ km/s and $v_{\rm sys}=-6.5\pm0.3$ km/s. Our values agree with previous reports. By accounting for WASP-33b's brightness variation, we find evidence that its night-side flux is $<10\%$ of the day-side flux and the emission peak is shifted westward of the substellar point, assuming the spectrum is dominated by Fe I. Our ESPaDOnS data, which cover phases before and after the secondary eclipse more evenly, weakly constrain the phase offset to $+22\pm12$ degrees. We caution that the derived volume-mixing-ratio depends on our choice of temperature-pressure profile, but note it does not significantly influence our constraints on day-night contrast or phase offset. Finally, we use simulations to illustrate how observations with increased phase coverage and higher signal-to-noise ratios can improve these constraints, showcasing the expanding capabilities of high-resolution Doppler spectroscopy.

We investigate a concrete scenario of a light scalar with a mass around 1 MeV which can be connected to the origin of neutrino masses and simultaneously survive current bounds on relativistic degrees of freedom in the early universe. In particular we show that a feeble coupling to the Standard Model neutrinos can relax the stringent bounds on the decays to photons inferred from the measured value of $N_{\rm eff}$. Interestingly, we find that such a scalar whose diphoton coupling is generated by a tree-level coupling to the muon of similar strength as that of the Standard Model Higgs boson can simultaneously explain the long-standing discrepancy in the measured value of the muon magnetic moment. We present a possible ultraviolet completion of this scenario providing a link between new physics in the early universe and the generation of neutrino masses.

Any new vector boson with non-zero mass (a `dark photon' or `Proca boson') that is present during inflation is automatically produced at this time from vacuum fluctuations and can comprise all or a substantial fraction of the observed dark matter density, as shown by Graham, Mardon, and Rajendran. We demonstrate, utilising both analytic and numerical studies, that such a scenario implies an extremely rich dark matter substructure arising purely from the interplay of gravitational interactions and quantum effects. Due to a remarkable parametric coincidence between the size of the primordial density perturbations and the scale at which quantum pressure is relevant, a substantial fraction of the dark matter inevitably collapses into gravitationally bound solitons, which are fully quantum coherent objects. The central densities of these `dark photon star', or `Proca star', solitons are typically a factor $10^6$ larger than the local background dark matter density, and they have characteristic masses of $10^{-16} M_\odot (10^{-5}{\rm eV}/m)^{3/2}$, where $m$ is the mass of the vector. During and post soliton production a comparable fraction of the energy density is initially stored in, and subsequently radiated from, long-lived quasi-normal modes. Furthermore, the solitons are surrounded by characteristic `fuzzy' dark matter halos in which quantum wave-like properties are also enhanced relative to the usual virialized dark matter expectations. Lower density compact halos, with masses a factor of $\sim 10^5$ greater than the solitons, form at much larger scales. We argue that, at minimum, the solitons are likely to survive to the present day without being tidally disrupted. This rich substructure, which we anticipate also arises from other dark photon dark matter production mechanisms, opens up a wide range of new direct and indirect detection possibilities, as we discuss in a companion paper.

Gravitational wave observations of large mass ratio compact binary mergers like GW190814 highlight the need for reliable, high-accuracy waveform templates for such systems. We present NRHybSur2dq15, a new surrogate model trained on hybridized numerical relativity (NR) waveforms with mass ratios $q\leq15$, and aligned spins $|\chi_{1z}|\leq0.5$ and $\chi_{2z}=0$. We target the parameter space of GW190814-like events as large mass ratio NR simulations are very expensive. The model includes the (2,2), (2,1), (3,3), (4,4), and (5,5) spin-weighted spherical harmonic modes, and spans the entire LIGO bandwidth (with $f_{\mathrm{low}}=20$ Hz) for total masses $M \gtrsim 9.5 \, M_{\odot}$. NRHybSur2dq15 accurately reproduces the hybrid waveforms, with mismatches below $\sim 2 \times 10^{-3}$ for total masses $10 \, M_{\odot} \leq M \leq 300 \, M_{\odot}$. This is at least an order of magnitude improvement over existing semi-analytical models for GW190814-like systems. Finally, we reanalyze GW190814 with the new model and obtain source parameter constraints consistent with previous work.

In this paper we propose a new mathematical model capable of merging Darwinian Evolution, Human History and SETI into a single mathematical scheme: 1) Darwinian Evolution over the last 3.5 billion years is defined as one particular realization of a certain stochastic process called Geometric Brownian Motion (GBM). This GBM yields the fluctuations in time of the number of species living on Earth. Its mean value curve is an increasing exponential curve, i.e. the exponential growth of Evolution. 2) In 2008 this author provided the statistical generalization of the Drake equation yielding the number N of communicating ET civilizations in the Galaxy. N was shown to follow the lognormal probability distribution. 3) We call "b-lognormals" those lognormals starting at any positive time b ("birth") larger than zero. Then the exponential growth curve becomes the geometric locus of the peaks of a one-parameter family of b-lognormals: this is our way to re-define Cladistics. 4) b-lognormals may be also be interpreted as the lifespan of any living being (a cell, or an animal, a plant, a human, or even the historic lifetime of any civilization). Applying this new mathematical apparatus to Human History, leads to the discovery of the exponential progress between Ancient Greece and the current USA as the envelope of all b-lognormals of Western Civilizations over a period of 2500 years. 5) We then invoke Shannon's Information Theory. The b-lognormals' entropy turns out to be the index of "development level" reached by each historic civilization. We thus get a numerical estimate of the entropy difference between any two civilizations, like the Aztec-Spaniard difference in 1519. 6) In conclusion, we have derived a mathematical scheme capable of estimating how much more advanced than Humans an Alien Civilization will be when the SETI scientists will detect the first hints about ETs.

The discovery of new exoplanets makes us wonder where each new exoplanet stands along its way to develop life as we know it on Earth. Our Evo-SETI Theory is a mathematical way to face this problem. We describe cladistics and evolution by virtue of a few statistical equations based on lognormal probability density functions (pdf) in the time. We call b-lognormal a lognormal pdf starting at instant b (birth). Then, the lifetime of any living being becomes a suitable b-lognormal in the time. Next, our "Peak-Locus Theorem" translates cladistics: each species created by evolution is a b-lognormal whose peak lies on the exponentially growing number of living species. This exponential is the mean value of a stochastic process called "Geometric Brownian Motion" (GBM). Past mass extinctions were all-lows of this GBM. In addition, the Shannon Entropy (with a reversed sign) of each b-lognormal is the measure of how evolved that species is, and we call it EvoEntropy. The "molecular clock" is re-interpreted as the EvoEntropy straight line in the time whenever the mean value is exactly the GBM exponential. We were also able to extend the Peak-Locus Theorem to any mean value other than the exponential. For example, we derive in this paper for the first time the EvoEntropy corresponding to the Markov-Korotayev (2007) "cubic" evolution: a curve of logarithmic increase.

In this paper, we implement the Adaptive Moving Mesh method (AMM) to the solution of initial value problems involving the Schr\"odinger equation, and more specifically the Schr\"odinger-Poisson system of equations. This method is based on the solution of the problem on a discrete domain, whose resolution is coordinate and time-dependent, and allows to dynamically assign numerical resolution in terms of desired refinement criteria. We apply the method to solve various test problems involving stationary solutions of the SP system, and toy scenarios related to the disruption of subhalo s made of ultralight bosonic dark matter traveling on top of host galaxies.

We discuss two-stage dilaton-axion inflation models [1] and describe $\alpha$-attractor models with either exponential or polynomial approach to the plateau. We implement one of the models of primordial black hole production proposed in [2] in the $\alpha$-attractor context, and develop its supergravity version. The predictions of this model following from its polynomial attractor properties are: $n_s$ and $r$ are $\alpha$-independent, $r$ depends on the mass parameter $\mu$ defining the approach to the plateau. The tachyonic instability at the transition point between the two stages of inflation is proportional to the negative curvature of the hyperbolic space $\mathcal{R}_K=-2/3\alpha$. Therefore the masses of primordial black holes (PBHs) and the frequencies of small-scale gravitational waves (GWs) in this model show significant dependence on $\alpha$.

We present a novel approach to modeling the ground state mass of atomic nuclei based directly on a probabilistic neural network constrained by relevant physics. Our Physically Interpretable Machine Learning (PIML) approach incorporates knowledge of physics by using a physically motivated feature space in addition to a soft physics constraint that is implemented as a penalty to the loss function. We train our PIML model on a random set of $\sim$20\% of the Atomic Mass Evaluation (AME) and predict the remaining $\sim$80\%. The success of our methodology is exhibited by the unprecedented $\sigma_\textrm{RMS}\sim186$ keV match to data for the training set and $\sigma_\textrm{RMS}\sim316$ keV for the entire AME with $Z \geq 20$. We show that our general methodology can be interpreted using feature importance.

In order to describe inflation in general relativity, scalar fields must inevitably be used, with all the setbacks of that description. On the other hand, $f(R)$ gravity and other modified gravity theories seem to provide a unified description of early and late-time dynamics without resorting to scalar or phantom theories. The question is, can modified gravity affect directly the mysterious radiation domination era? Addressing this question is the focus in this work, and we shall consider the case for which in the early stages of the radiation domination era, namely during the reheating era, the background equation of state parameter is different from $w=1/3$. As we show, in the context of $f(R)$ gravity, an abnormal reheating era can affect the primordial gravitational wave energy spectrum today. Since future interferometers will exactly probe this era, which consists of subhorizon modes that reentered the horizon during the early stages of the radiation domination era, the focus in this work is how a short abnormal reheating era that deviates from the standard perfect fluid pattern with $w\neq 1/3$, and generated by higher order curvature terms, can affect the primordial gravitational wave energy spectrum. Using a WKB approach, we calculate the effect of an $f(R)$ gravity generated abnormal reheating era, and as we show the primordial gravitational wave spectrum is significantly amplified, a result which is in contrast to the general relativistic case, where the effect is minor.

We investigate the dynamics of purely kinetic $k$-essence scalar field in three different scenarios. We begin with the case of a purely kinetic $k$-essence theory and try to formulate a late time accelerating universe model. Then we introduce a perfect fluid in presence of the purely kinetic $k$-essence field and calculate how the features of the previous model are modified. We discuss two cases primarily, in the first case the fluid and the scalar field does not exchange energy and momentum and in the later case the scalar field and the fluid interact non-minimally and exchange energy and momentum. Using dynamical analysis technique, we examine the stability of these systems in the background of a spatially flat Friedman-Lemaitre-Robertson-Walker (FLRW) universe for a particular chosen functional form of $k$-essence kinetic energy term. From the evolution of the effective equation of state parameter of the field-fluid system we conclude that a non-minimally coupled field-fluid scenario can produce physically interesting results. We also find that the $k$-essence field can effectively hide in a radiation background.

The characteristics of Alfv\'en waves propagating in a direction oblique to the ambient magnetic field in a stellar wind environment are discussed. A kinetic formulation for a magnetized dusty plasma is adopted considering Maxwellian distributions of electrons and ions, and immobile dust particles electrically charged by absorption of plasma particles and by photoionization. The dispersion relation is numerically solved and the results are compared with situations previously studied where dust particles were not charged by photoionization, which is an important process in a stellar wind of a relatively hot star. We show that the presence of dust causes the shear Alfv\'en waves to present a region of wavenumber values with zero frequency and that the minimum wavelength for which the mode becomes dispersive again is roughly proportional to the radiation intensity to which the dust grains are exposed. The damping rates of both shear and compressional Alfv\'en waves are observed to decrease with increasing radiation flux, for the parameters considered. For the particular case where both modes present a region with null real frequency when the radiation flux is absent or weak, it is shown that when the radiation flux is sufficiently strong, the photoionization mechanism may cause this region to get smaller or even to vanish, for compressional Alfv\'en waves. In that case, the compressional Alfv\'en waves present non zero frequency for all wavenumber values, while the shear Alfv\'en waves still present null frequency in a certain interval of wavenumber values, which gets smaller with the presence of radiation.

Nonlinear electrodynamics, which acts as a source of gravity Einstein field equations, leads to emergent cosmology, an alternative solution which can avoid Big Bang singularity. In this paper, we explore the emerging universe in models of non-linear electrodynamics (described by dimensional parameter $\beta$) by using the equation of state parameter $\omega$ and see how the parameter $\beta$ helps the universe to cause a transition from a quasi-static Minkowski phase to the inflationary phase of expansion through the point of emergence and subsequently to the phase of normal thermal expansion. We predict the spectral index parameter $n_s = 0.97467$ (scalar spectral index), $r =0.10133$ (tensor to scalar ratio) and $n_T = -0.01267$ (tensor spectral index) of the inflationary perturbation in emergent cosmology of nonlinear electrodynamics corresponding to $\beta$ = 0.1 and $B_0=10^{-10}$G.

Timing a pulsar in a close orbit around the supermassive black hole SgrA* at the center of the Milky Way would open the window for an accurate determination of the black hole parameters and for new tests of General Relativity and alternative modified gravity theories. An important relativistic effect which has to be taken into account in the timing model is the propagation delay of the pulses in the gravitational field of the black hole. Due to the extreme mass ratio of the pulsar and the supermassive back hole we use the test particle limit to derive an exact analytical formula for the propagation delay in a Kerr spacetime and deduce a relativistic formula for the frame dragging effect on the arrival time. As an illustration, we treat an edge-on orbit in which the frame dragging effect is expected to be maximal. We compare our formula for the propagation time delay with Post-Newtonian approaches, and in particular with the frame dragging terms derived in previous works by Wex & Kopeikin and Rafikov & Lai. Our approach correctly identifies the asymmetry of the frame dragging delay with respect to superior conjunction, avoids singularities in the time delay, and indicates that in the Post-Newtonian approach frame dragging effects are generally slightly overestimated.

We construct an effective four dimensional string-corrected black hole (4D SCBH) by rescaling the string coupling parameter in a $D$-dimensional Callan-Myers-Perry black hole. From the theoretical point of view, the most interesting findings are that the string corrections coincide with the so-called generalized uncertainty principle (GUP) corrections to black hole solutions, Bekenstein-Hawking entropy acquires logarithmic corrections, and that there exists a critical value of the coupling parameter for which the black hole temperature vanishes. We also find that, due to the string corrections the nature of the central singularity may be altered from spacelike to timelike singularity. In addition, we study the possibility of testing such a black hole with astrophysical observations. Since the dilaton field does not decouple from the metric it is not a priori clear that the resulting 4D SCBH offers only small corrections to the Schwarzschild black hole. We used motion of the S2 star around the black hole at the center of our galaxy to constrain the parameters (the string coupling parameter and ADM mass) of the 4D SCBH. To test the weak gravity regime we calculate the deflection angle in this geometry and apply it to gravitational lensing. To test the strong field regime, we calculate the black hole shadow radius. While we find that the observables change as we change the string coupling parameter, the magnitude of the change is too small to distinguish it from the Schwarzschild black hole. With the current precision, to the leading order terms, the 4D SCBH cannot be distinguished from the Schwarzschild black hole.

We present a new model for the probability that the Disturbance storm time (Dst) index exceeds -100 nT, with a lead time between 1 and 3 days. $Dst$ provides essential information about the strength of the ring current around the Earth caused by the protons and electrons from the solar wind, and it is routinely used as a proxy for geomagnetic storms. The model is developed using an ensemble of Convolutional Neural Networks (CNNs) that are trained using SoHO images (MDI, EIT and LASCO). The relationship between the SoHO images and the solar wind has been investigated by many researchers, but these studies have not explicitly considered using SoHO images to predict the $Dst$ index. This work presents a novel methodology to train the individual models and to learn the optimal ensemble weights iteratively, by using a customized class-balanced mean square error (CB-MSE) loss function tied to a least-squares (LS) based ensemble. The proposed model can predict the probability that Dst<-100 nT 24 hours ahead with a True Skill Statistic (TSS) of 0.62 and Matthews Correlation Coefficient (MCC) of 0.37. The weighted TSS and MCC from Guastavino et al. (2021) is 0.68 and 0.47, respectively. An additional validation during non-Earth-directed CME periods is also conducted which yields a good TSS and MCC score.

We introduce a simulator of charge transport in fully-depleted, thick CCDs that include Coulomb repulsion between carriers. The calculation of this long-range interaction is highly intensive computationally, and only a few thousands of carriers can be simulated in reasonable times using regular CPUs. G-CoReCCD takes advantage of the high number of multiprocessors available in a graphical processing unit (GPU) to parallelize the operations and thus achieve a massive speedup. We can simulate the path inside the CCD bulk for up to hundreds of thousands of carriers in only a few hours using modern GPUs.