Papers for Wednesday, Oct 13 2021

We present a list of candidate tailed radio galaxies using the TIFR GMRT Sky Survey Alternative Data Release 1 (TGSS ADR1) at 150 MHz. We have visually examined 5336 image fields and found 268 candidates for the tailed galaxy. Tailed radio galaxies are classified as Wide Angle Tailed (WAT) galaxies and Narrow-Angle Tailed (NAT) galaxies, based on the angle between the two jets of the galaxy. We found a sample of tailed radio galaxies which includes 189 `WAT' and 79 `NAT' type sources. These newly identified tailed sources are significant inclusion to the list of known tailed radio galaxies. The source morphology, luminosity feature of the different candidate galaxies and their optical identifications are presented in the paper. Other radio properties and general features of the sources are also discussed.

We present millimeter (80-230 GHz), radio (6-45 GHz), and X-ray (0.2-10 keV) observations of ZTF20acigmel (AT2020xnd), a short-duration luminous optical transient at $z=0.2433$. The 100 GHz peak luminosity is similar to that of long-duration gamma-ray bursts ($2\times10^{30}$ erg sec$^{-1}$ Hz$^{-1}$) but the light curve rises on a much longer timescale (one month). In the standard framework of synchrotron self-absorption of electrons in a power-law energy distribution, the data imply a fast ($v\approx0.2c$) shock with large energy ($U\gtrsim 10^{49}$ erg) propagating in a medium with a steep ($n_e \propto r^{-3}$) density profile. The forward-shock properties are similar to those of the fast-luminous transient AT2018cow, and in both cases the model for the late-time ($\Delta t>70$ d) low-frequency ($\nu<40$ GHz) data is not consistent with the early-time ($\Delta t<40$d) high-frequency ($\nu>70$ GHz) emission. Motivated by the observation of a steep spectral index ($f_\nu \propto \nu^{-2}$) across the millimeter bands, we favor a thermal electron population (relativistic Maxwellian) for the synchrotron emission, the first such inference for a cosmic explosion. We find that the X-ray luminosity of $L_X\approx10^{43}$erg sec$^{-1}$ exceeds simple predictions from the radio and UVOIR luminosity and likely has a separate physical origin, such as a central engine. Our work suggests that luminous millimeter, radio, and X-ray emission are a generic feature of transients with fast ($\approx3$ days) and luminous ($M\approx-21$ mag) optical light curves. We estimate the rate at which transients like AT2018cow and AT2020xnd will be detected by future wide-field millimeter transient surveys like CMB-S4, and conclude that energetic explosions in dense environments may represent a significant population of extragalactic transients in the 100 GHz sky.

The next generation of wide-field cosmic microwave background (CMB) surveys are uniquely poised to open a new window for time-domain astronomy in the millimeter band. Here we explore the discovery phase space for extragalactic transients with near-term and future CMB experiments to characterize the expected population. We use existing millimeter-band light curves of known transients (gamma-ray bursts, tidal disruption events, fast blue optical transients, neutron star mergers) and theoretical models, in conjunction with known and estimated volumetric rates. Using Monte Carlo simulations of various CMB survey designs (area, cadence, depth, duration) we estimate the detection rates and the resulting light curve characteristics. We find that existing and near-term surveys will find tens to hundreds of long-duration gamma-ray bursts (LGRBs), driven primarily by detections of the reverse shock emission, and including off-axis LGRBs. Next-generation experiments (CMB-S4, CMB-HD) will find tens of fast blue optical transients (FBOTs) in the nearby universe and will detect a few tidal disruption events. CMB-HD will additionally detect a small number of short gamma-ray bursts, where these will be discovered within the detection volume of next generation gravitational wave experiments like Cosmic Explorer.

The globular cluster 47 Tucanae (47 Tuc) is one of the most massive star clusters in the Milky Way and is exceptionally rich in exotic stellar populations. For several decades it has been a favorite target of observers, and yet it is computationally very challenging to model because of its large number of stars ($N\gtrsim 10^6$) and high density. Here we present detailed and self-consistent 47 Tuc models computed with the \texttt{Cluster Monte Carlo} code (\texttt{CMC}). The models include all relevant dynamical interactions coupled to stellar and binary evolution, and reproduce various observations, including the surface brightness and velocity dispersion profiles, pulsar accelerations, and numbers of compact objects. We show that the present properties of 47 Tuc are best reproduced by adopting an initial stellar mass function that is both bottom-heavy and top-light relative to standard assumptions \citep[as in, e.g.,][]{Kroupa2001}, and an initial Elson profile \citep{Elson1987} that is overfilling the cluster's tidal radius. We include new prescriptions in \texttt{CMC} for the formation of binaries through giant star collisions and tidal captures, and we show that these mechanisms play a crucial role in the formation of neutron star binaries and millisecond pulsars in 47 Tuc; our best-fit model contains $\sim 50$ millisecond pulsars, $80\%$ of which are formed through giant collisions and tidal captures. Our models also suggest that 47 Tuc presently contains up to $\sim 200$ stellar-mass black holes, $\sim 5$ binary black holes, $\sim 15$ low-mass X-ray binaries, and $\sim 300$ cataclysmic variables.

Context: Modelling satellite galaxy abundance $N_s$ in Galaxy Clusters (GCs) is a key element in modelling the Halo Occupation Distribution (HOD), which itself is a powerful tool to connect observational studies with numerical simulations. Aims: To study the impact of cosmological parameters on satellite abundance both in cosmological simulations and in mock observations. Methods: We build an emulator (HODEmu, \url{https://github.com/aragagnin/HODEmu/}) of satellite abundance based on cosmological parameters $\Omega_m, \Omega_b, \sigma_8, h_0$ and redshift $z.$ We train our emulator using \magneticum hydrodynamic simulations that span 15 different cosmologies, each over $4$ redshift slices between $0<z<0.5,$ and for each setup we fit normalisation $A$, log-slope $\beta$ and Gaussian fractional-scatter $\sigma$ of the $N_s-M$ relation. The emulator is based on multi-variate output Gaussian Process Regression (GPR). Results: We find that $A$ and $\beta$ depend on cosmological parameters, even if weakly, especially on $\Omega_m,$ $\Omega_b.$ This dependency can explain some discrepancies found in literature between satellite HOD of different cosmological simulations (Magneticum, Illustris, BAHAMAS). We also show that satellite abundance cosmology dependency differs between full-physics (FP) simulations, dark-matter only (DMO), and non-radiative simulations. Conclusions: This work provides a preliminary calibration of the cosmological dependency of the satellite abundance of high mass halos, and we showed that modelling HOD with cosmological parameters is necessary to interpret satellite abundance, and we showed the importance of using FP simulations in modelling this dependency.

We use N-body cosmological simulations and empirical galaxy models to study the merger history of dwarf-mass galaxies (with M_halo~10^10 M_Sun). Our input galaxy models describe the stellar mass-halo mass relation, and the galaxy occupation fraction. The number of major and minor mergers depends on the type of dark matter; in particular, minor mergers are greatly suppressed in warm dark matter models. In addition, the number of mergers that bring in stars is strongly dependent on the galaxy occupation model. For example, minor mergers are negligible for stellar halo growth in models with a high mass threshold for galaxy formation (i.e. 10^9.3 M_Sun at z=0). Moreover, this threshold for galaxy formation can also determine the relative difference (if any) between the stellar haloes of satellite and field dwarfs. Using isolated simulations of dwarf-dwarf mergers, we show that the relative frequency of major and minor mergers predict very different stellar haloes: typically, "intermediate" dark matter merger ratios (~1:5) maximise the growth of distant stellar haloes. We discuss the observability of dwarf stellar haloes and find that the surface brightness of these features are incredibly faint. However, when several dwarfs are stacked together models that form particularly rich stellar haloes could be detectable. Finally, we show that stellar streams in the Galactic halo overlapping in phase-space with known dwarf satellites are likely remnants of their stripped stellar haloes. The mere existence of dwarf stellar haloes can already put constraints on some small-scale models, and thus observational probes should be a high priority.

Gravitational microlensing is a powerful tool to probe the inner structure of strongly lensed quasars and to constrain parameters of the stellar mass function of lens galaxies. This is done by analysing microlensing light curves between the multiple images of strongly lensed quasars, under the influence of three main variable components: 1- the continuum flux of the source, 2- microlensing by stars in the lens galaxy and 3- reverberation of the continuum by the Broad Line Region (BLR). The latter, ignored by state-of-the-art microlensing techniques, can introduce high-frequency variations which we show carry information on the BLR size. We present a new method which includes all these components simultaneously and fits the power spectrum of the data in the Fourier space, rather than the observed light curve itself. In this new framework, we analyse COSMOGRAIL light curves of the two-image system QJ0158-4325 known to display high-frequency variations. Using exclusively the low frequency part of the power spectrum our constraint on the accretion disk radius agrees with the thin disk model estimate and previous work that fit the microlensing light curves in real space. However, if we also take into account the high-frequency variations, the data favour significantly smaller disk sizes than previous microlensing measurements. In this case, our results are in agreement with the thin disk model prediction only if we assume very low mean masses for the microlens population, i.e. <M> = 0.01 $M_\odot$. Eventually, including the differentially microlensed continuum reverberation by the BLR successfully explains the high frequencies without requiring such low mass microlenses. This allows us to measure, for the first time, the size of the BLR using single-band photometric monitoring, $R_{BLR}$ = $1.6^{+1.5}_{-0.8}\times 10^{17}$cm, in agreement with estimates using the BLR size-luminosity relation.

We present the analysis of the properties of large samples of protostellar discs formed in four radiation hydrodynamical simulations of star cluster formation. The four calculations have metallicities of 0.01, 0.1, 1 and 3 times solar metallicity. The calculations treat dust and gas temperatures separately and include a thermochemical model of the diffuse interstellar medium. We find the radii of discs of bound protostellar systems tend to decrease with decreasing metallicity, with the median characteristic radius of discs in the 0.01 and 3 times solar metallicity calculations being $\approx20$ and $\approx65$ au, respectively. Disc masses and radii of isolated protostars also tend to decrease with decreasing metallicity. We find that the circumstellar discs and orbits of bound protostellar pairs, and the two spins of the two protostars are all less well aligned with each other with lower metallicity than with higher metallicity. These variations with metallicity are due to increased small scale fragmentation due to lower opacities and greater cooling rates with lower metallicity, which increase the stellar multiplicity and increase dynamical interactions. We compare the disc masses and radii of protostellar systems from the solar metallicity calculation with recent surveys of discs around Class 0 and I objects in the Orion and Perseus star-forming regions. The masses and radii of the simulated discs have similar distributions to the observed Class 0 and I discs.

We combine multiple campaigns of K2 photometry with precision radial velocity measurements from Keck-HIRES to measure the masses of three sub-Neptune-size planets. We confirm the planetary nature of the massive sub-Neptune K2-182 b ($P_\mathrm{b}= 4.7$ days, $R_\mathrm{b} = 2.69$ $R_\oplus$) and derive refined parameters for K2-199 b and c ($P_\mathrm{b} = 3.2$ days, $R_\mathrm{b} = 1.73$ $R_\oplus$, and $P_\mathrm{c} = 7.4$ days, $R_\mathrm{c} = 2.85$ $R_\oplus$). These planets provide valuable data points in the mass-radius plane, especially as TESS continues to reveal an increasingly diverse sample of sub-Neptunes. K2-182 (EPIC 211359660) is a moderately bright ($V = 12.0$ mag) early-K dwarf observed during K2 campaigns 5 and 18. We find K2-182 b is potentially one of the densest sub-Neptunes known to date ($20 \pm 5$ $M_\oplus$ and $5.6 \pm 1.4$ g cm$^{-3}$). K2-199 (EPIC 212779596; $V = 12.3$ mag) is a K5V dwarf observed in K2 campaigns 6 and 17 which hosts two recently-confirmed planets. We refine the orbital and planetary parameters for K2-199 b and c by modeling both campaigns of K2 photometry and adding 12 Keck-HIRES measurements to the existing radial velocity data set ($N$ = 33). We find K2-199 b is likely rocky, at $6.9 \pm 1.8$ $M_\oplus$ and $7.2^{+2.1}_{-2.0}$ g cm$^{-3}$. K2-199 c has an intermediate density at $12.4 \pm 2.3$ $M_\oplus$ and $2.9^{+0.7}_{-0.6}$ g cm$^{-3}$. We contextualize these planets on the mass-radius plane, discuss a small but intriguing population of "super-dense" sub-Neptunes ($R_\mathrm{p} < 3$ $R_\oplus$, $M_\mathrm{p} > 20$ $M_\oplus$), and consider our prospects for the planets' atmospheric characterization.

Writing is a vital component of a modern career in astronomical research. Very few researchers, however, receive any training in how to produce high-quality written work in an efficient manner. We present a step-by-step guide to writing in astronomy. We concentrate on how to write scientific papers, and address various aspects including how to crystallise the ideas that underlie the research project, and how the paper is constructed considering the audience and the chosen journal. We also describe a number of grammar and spelling issues that often cause trouble to writers, including some that are particularly hard to master for non-native English speakers. This paper is aimed primarily at Master's and PhD level students who are presented with the daunting task of writing their first scientific paper, but more senior researchers or writing instructors may well find the ideas presented here useful.

Population studies of stellar-mass black-hole binaries have become major players in gravitational-wave astronomy. The underlying assumptions are that the targeted source parameters refer to the same quantities for all events in the catalog and are included when modeling selection effects. Both these points have so far been neglected when estimating the orientations of the black-hole spins. In particular, the detector-frame gravitational-wave frequency used to define frequency-dependent quantities (e.g., 20 Hz) introduces an inconsistent reference between events at the population level. We solve both issues by modeling binary black-hole populations and selection effects at past time infinity, corresponding to the well-defined reference frequency of 0 Hz. We show that, while current gravitational-wave measurement uncertainties obfuscate the influence of reference frequency in population inference, ignoring spins when estimating selection effects leads to an over-prediction of spin alignment in the underlying astrophysical distribution of merging black holes.

High cadence, high quality observations of active galactic nuclei (AGN) clearly show continuum variations with lags, relative to the shortest observed variable UV continuum, that increase with wavelength ("lag spectra"). These have been attributed to the irradiation and heating of the central accretion disk by the central X-ray emitting corona. An alternative explanation, connecting the observed lag-spectra to line and continuum emission from gas in the broad line region (BLR), has also been proposed. In this paper I show the clear spectral signature of the time-dependent diffuse gas emission in the lag-spectrum of 6 AGN. I also show a new lag-luminosity relationship for 9 objects which is a scaled down version of the well known lag(H-beta)-L(5100A) relationship in AGN. The shape of the lag-spectrum, and its normalization, are entirely consisted with diffuse emission from radiation pressure supported clouds in a BLR with a covering factor of about 0.2. While some contributions to the continuum lag from the irradiated disk cannot be excluded, there is no need for this explanation.

We present deep X-ray and radio observations of the Fast Blue Optical Transient (FBOT) AT2020xnd/ZTF20acigmel at $z=0.2433$ from $13$d to $269$d after explosion. AT2020xnd belongs to the category of optically luminous FBOTs with similarities to the archetypal event AT2018cow. AT2020xnd shows luminous radio emission reaching $L_{\nu}\approx8\times10^{29}$ergs$^{-1}$Hz$^{-1}$ at 20GHz and $75$d post explosion, accompanied by luminous and rapidly fading soft X-ray emission peaking at $L_{X}\approx6\times10^{42}$ergs$^{-1}$. Interpreting the radio emission in the context of synchrotron radiation from the explosion's shock interaction with the environment we find that AT2020xnd launched a high-velocity outflow ($v\sim$0.1-0.2$c$) propagating into a dense circumstellar medium (effective $\dot M\approx10^{-3}M_{\rm{sol}}$yr$^{-1}$ for an assumed wind velocity of $v_w=1000$kms$^{-1}$). Similar to AT2018cow, the detected X-ray emission is in excess compared to the extrapolated synchrotron spectrum and constitutes a different emission component, possibly powered by accretion onto a newly formed black hole or neutron star. These properties make AT2020xnd a high-redshift analog to AT2018cow, and establish AT2020xnd as the fourth member of the class of optically-luminous FBOTs with luminous multi-wavelength counterparts.

We present a new method for consistent, joint analysis of the pre- and post-reconstruction two-point functions of the BOSS survey. The post-reconstruction correlation function is used to accurately measure the distance-redshift relation and expansion history, while the pre-reconstruction power spectrum multipoles constrain the broad-band shape and the rate-of-growth of large-scale structure. Our technique uses Lagrangian perturbation theory to self-consistently work at the level of two-point functions, i.e.\ directly with the measured data, without approximating the constraints with summary statistics normalized by the drag scale. Combining galaxies across the full redshift range and both hemispheres we constrain $\Omega_m=0.304 \pm 0.0084$, $H_0=69.22 \pm 0.79$ and $\sigma_8=0.783 \pm 0.047$ within the context of $\Lambda$CDM. These constraints are in good agreement both with the Planck primary CMB anisotropy data and recent cosmic shear surveys.

The presence of planetary material in white dwarf atmospheres, thought to be accreted from a dusty debris disc produced via the tidal disruption of a planetesimal, is common. Approximately five per cent of these discs host a co-orbital gaseous component detectable via emission from atomic transitions - usually the 8600 Angstrom CaII triplet. These emission profiles can be highly variable in both morphology and strength. Furthermore, the morphological variations in a few systems have been shown to be periodic, likely produced by an apsidally precessing asymmetric disc. Of the known gaseous debris discs, that around HE1349-2305 has the most rapidly evolving emission line morphology, and we present updated spectroscopy of the CaII triplet of this system. The additional observations show that the emission line morphologies vary periodically and consistently, and we constrain the period to two aliases of 459$\pm$3d and 502$\pm$3d. We produce images of the CaII triplet emission from the disc in velocity space using Doppler tomography - only the second such imaging of a white dwarf debris disc. We suggest that the asymmetric nature of these velocity images is generated by gas moving on eccentric orbits with radially-dependent excitation conditions via photo-ionisation from the white dwarf. We also obtained short-cadence (~ 4 min) spectroscopy to search for variability on the time-scale of the disc's orbital period (~ hours) due to the presence of a planetesimal, and rule out variability at a level of ~ 1.4 per cent.

Based on observations with the IRAM 30-m and Yebes 40-m telescopes, we report evidence of the detection of Milky Way-like, low-excitation molecular gas, up to the transition CO($J=5-4$), in a distant, dusty star-forming galaxy at $z_{CO}=1.60454$. WISE J122651.0+214958.8 (alias SDSSJ1226, the Cosmic Seahorse), is strongly lensed by a foreground galaxy cluster at $z=0.44$ with a source magnification of $\mu=9.5\pm0.7$. This galaxy was selected by cross-correlating near-to-mid infrared colours within the full-sky AllWISE survey, originally aiming to discover rare analogs of the archetypical strongly lensed submillimeter galaxy SMM J2135-0102, the Cosmic Eyelash. We derive an apparent (i.e. not corrected for lensing magnification) rest-frame 8-1000 $\mu$m infrared luminosity of $\mu L_\mathrm{IR}=1.66^{+0.04}_{-0.04}\times 10^{13}$ L$_\odot$ and apparent star-formation rate $\mu\mathrm{SFR}_\mathrm{IR}=2960\pm70$ M$_\odot$ yr$^{-1}$. SDSSJ1226 is ultra-bright at $S_{350\mu m}\simeq170$ mJy and shows similarly bright low-$J$ CO line intensities as SMM J2135-0102, however, with exceptionally small CO($J=5-4$) intensity. We consider different scenarios to reconcile our observations with typical findings of high-$z$ starbursts, and speculate about the presence of a previously unseen star-formation mechanism in cosmic noon submillimeter galaxies. In conclusion, the remarkable low line luminosity ratio $r_{5,2}=0.11\pm0.02$ is best explained by an extended, main-sequence star-formation mode -- representing a missing link between starbursts to low-luminosity systems during the epoch of peak star-formation history.

We investigate whether the recently observed 2.6 $M_\odot$ compact object in the gravitational-wave event GW190814 can be a bosonic dark matter admixed compact star. By considering the three constraints in mass, radius and stability of such an object, we find that if the dark matter is made of QCD axions, their particle mass $m$ is constrained to a range that has already been ruled out by the independent constraint imposed by the stellar-mass black hole superradiance process. The 2.6 $M_\odot$ object can still be a neutron star admixed with at least 2.0 $M_\odot$ of dark matter made of axion-like particles (or even a pure axion-like particle star) if $2 \times 10^{-11}$ eV $\leq m \leq 2.4 \times 10^{-11}$ eV ($2.9 \times 10^{-11}$ eV $\leq m \leq 3.2 \times 10^{-11}$ eV) and with decay constant $f \geq 8 \times 10^{17}$ GeV.

The Gravitational-wave Optical Transient Observer (GOTO) is an array of wide-field optical telescopes, designed to exploit new discoveries from the next generation of gravitational wave detectors (LIGO, Virgo, KAGRA), study rapidly evolving transients, and exploit multi-messenger opportunities arising from neutrino and very high energy gamma-ray triggers. In addition to a rapid response mode, the array will also perform a sensitive, all-sky transient survey with few day cadence. The facility features a novel, modular design with multiple 40-cm wide-field reflectors on a single mount. In June 2017 the GOTO collaboration deployed the initial project prototype, with 4 telescope units, at the Roque de los Muchachos Observatory (ORM), La Palma, Canary Islands. Here we describe the deployment, commissioning, and performance of the prototype hardware, and discuss the impact of these findings on the final GOTO design. We also offer an initial assessment of the science prospects for the full GOTO facility that employs 32 telescope units across two sites.

We present a magnitude-limited set of lightcurves for stars observed over the TESS Extended Mission, as extracted from full-frame images (FFIs) by MIT's Quick-Look Pipeline (QLP). QLP uses multi-aperture photometry to produce lightcurves for ~1 million stars each 27.4-day sector, which are then searched for exoplanet transits. The per-sector lightcurves for 9.1 million unique targets observed over the first year of the Extended Mission (Sectors 27 - 39) are available as High-Level Science Products (HLSP) on the Mikulski Archive for Space Telescopes (MAST). As in our TESS Primary Mission QLP HLSP delivery (Huang et al. 2020), our available data products include both raw and detrended flux time series for all observed stars brighter than TESS magnitude T = 13.5, providing the community with one of the largest sources of FFI-extracted lightcurves to date.

The nuclear equation of state is an important component in the evolution of core collapse supernovae. In this paper we make a survey of various equations of state in the literature and analyze their effect on spherical core-collapse models in which the effects of three-dimensional turbulence is modeled by a general relativistic formulation of Supernova Turbulence in Reduced dimensionality (STIR). We show that the viability of the explosion is quite EOS dependent and that it best correlates with the early-time interior entropy density of the proto-neutron star. We check that this result is not progenitor dependent, although low-mass progenitors show different explosion properties, due to the different pre-collapse nuclear composition. Larger central entropies also induce more vigorous proto-neutron-star convection in our one-dimensional turbulence model, as well as a wider convective layer.

Over the course of several years, stars trace helical trajectories as they traverse across the sky due to the combined effects of proper motion and parallax. It is well known that the gravitational pull of an unseen companion can cause deviations to these tracks. Several studies have pointed out that the astrometric mission Gaia will be able to identify a slew of new exoplanets, stellar binaries, and compact object companions with orbital periods as short as tens of days to as long as Gaia's lifetime. Here, we use mock astrometric observations to demonstrate that Gaia can identify and characterize black hole companions to luminous stars with orbital periods longer than Gaia's lifetime. Such astrometric binaries have orbital periods too long to exhibit complete orbits, and instead are identified through curvature in their characteristic helical paths. By simultaneously measuring the radius of this curvature and the orbital velocity, constraints can be placed on the underlying orbit. We quantify the precision with which Gaia can measure orbital accelerations and apply that to model predictions for the population of black holes orbiting stars in the stellar neighborhood. Although orbital degeneracies imply that many of the accelerations induced by hidden black holes could also be explained by faint low-mass stars, we discuss how the nature of certain putative black hole companions can be confirmed with high confidence using Gaia data alone.

Context: As a result of Titan's migration and Saturn's probable capture in secular spin-orbit resonance, recent works show that Saturn's obliquity could be steadily increasing today and may reach large values in the next billions of years. Satellites around high-obliquity planets are known to be unstable near their Laplace radius, but the approximations used so far are invalidated in this regime. Aims: We aim to investigate the behaviour of a planet and its satellite when the satellite crosses its Laplace radius while the planet is locked in secular spin-orbit resonance. Methods: We expand on previous works and revisit the concept of Laplace surface. We use it to build an averaged analytical model that couples the planetary spin-axis and satellite dynamics. Results: We show that the dynamics is organised around a critical point, S1, at which the phase-space structure is singular, located at 90{\deg} obliquity and near the Laplace radius. If the spin-axis precession rate of the planet is maintained fixed by a resonance while the satellite migrates outwards or inwards, then S1 acts as an attractor towards which the system is forced to evolve. When it reaches the vicinity of S1, the entire system breaks down, either because the planet is expelled from the resonance or because the satellite is ejected or collides into the planet. Conclusions: Provided that Titan's migration is not halted in the future, Titan and Saturn may reach instability between a few gigayears and several tens of gigayears from now, depending on Titan's migration rate. The evolution would destabilise Titan and drive Saturn towards an obliquity of 90{\deg}. Our findings may have important consequences for Uranus. They also provide a straightforward mechanism for producing transiting exoplanets with a face-on massive ring, a configuration that is often put forward to explain some super-puff exoplanets.

Cosmic variance from large-scale structure will be a major source of uncertainty for galaxy surveys at z > 6, but that same structure will also provide an opportunity to identify and study dense environments in the early Universe. Using a robust model for galaxy clustering, we directly incorporate large-scale densities into an inference framework that simultaneously measures the high-z (z > 6) UV luminosity function and the average matter density of each distinct volume in a survey. Through this framework, we forecast the performance of several major upcoming James Webb Space Telescope (JWST) galaxy surveys. We find that they can constrain field matter densities down to the theoretical limit imposed by Poisson noise and unambiguously identify over-dense (and under-dense) regions on transverse scales of tens of comoving Mpc. We also predict JWST will measure the luminosity function with a precision at z = 12 comparable to existing Hubble Space Telescope's constraints at z = 8 (and even better for the faint-end slope). We also find that wide-field surveys are especially important in distinguishing luminosity function models.

The advent of high-resolution spectroscopy as a method for exoplanet atmospheric characterization has expanded our capability to study non-transiting planets, increasing the number of planets accessible for observation. Many of the most favorable targets for atmospheric characterization are hot Jupiters, where we expect large spatial variation in physical conditions such as temperature, wind speed, and cloud coverage, making viewing geometry important. Three-dimensional models have generally simulated observational properties of hot Jupiters assuming edge-on viewing, which neglects planets without near edge-on orbits. As the first investigation of how orbital inclination manifests in high-resolution emission spectra, we use a General Circulation Model to simulate the atmospheric structure of Upsilon Andromedae b, a non-transiting hot Jupiter. In order to accurately capture scattering from clouds, we implement a generalized two-stream radiative transfer routine for inhomogeneous multiple scattering atmospheres. We compare models with and without clouds, as cloud coverage intensifies spatial variations. Cloud coverage increases both the net Doppler shifts and the variation of the continuum flux amplitude over the course of the planet's orbit. As orbital inclination decreases, four key features also decrease in both the clear and cloudy models: 1) the average continuum flux level, 2) the amplitude of the variation in continuum with orbital phase, 3) net Doppler shifts of spectral lines, and 4) Doppler broadening in the spectra. Models capable of treating inhomogeneous cloud coverage and different viewing geometries are critical in understanding high-resolution emission spectra, enabling an additional avenue to investigate these extreme atmospheres.

We present a pilot library of synthetic NUV, U, B, V, and I photometry of star clusters with stochastically sampled IMFs and ionized gas for initial masses, $M_i=10^3$, $10^4$, and $10^5$ M$_{\odot}$; $t=1$, 3, 4, and 8 Myr; $Z=0.014$ and $Z=0.002$; and log(U$_{\rm S}$) =-2 and -3. We compare the library with predictions from deterministic models and observations of isolated low-mass ($<10^4$ M$_{\odot}$) star clusters with co-spatial compact H\2 regions. The clusters are located in NGC 7793, one of the nearest galaxies observed as part of the \hst LEGUS and H$\alpha$-LEGUS surveys. 1) For model magnitudes that only account for the stars: a) the residual |deterministic mag - median stochastic mag| can be $\ge0.5$ mag, even for $M_i=10^5$ M$_{\odot}$; and b) the largest spread in stochastic magnitudes occurs when Wolf-Rayet stars are present. 2) For $M_i=10^5$ M$_{\odot}$: a) the median stochastic mag with gas can be $>$1.0 mag more luminous than the median stochastic magnitude without gas; and b) nebular emission lines can contribute with $>50\%$ and $>30\%$ to the total emission in the V and I bands, respectively. 3) Age-dating OB-star clusters via deterministic tracks in the U-B vs. V-I plane is highly uncertain at $Z=0.014$ for $M_i\sim10^3$ M$_{\odot}$ and $Z=0.002$ for $M_i\sim10^3-10^5$ M$_{\odot}$. 4) For low-mass clusters, the V-band extinction derived with stochastic models significantly depends on the value of log(U$_{\rm S}$). 5) The youngest clusters tend to have higher extinction. 6) The majority of clusters have multi-peaked age PDFs. 7) Finally, we discuss the importance of characterising the true variance in the number of stars per mass bin in nature.

Spirals in galaxies have long been thought to be caused by gravitational instability in the stellar component of the disk, but the precise mechanism had proved elusive. Tidal interactions, and perhaps bars, may provoke some spiral responses, but a self-excitation mechanism is still required for many galaxies. We survey the relevant observational data and aspects of disk dynamical theory. The origin of the recurring spiral patterns in simulations of isolated disk galaxies has recently become clear and it seems likely that the mechanism is the same in real galaxies, although evidence to confirm this supposition is hard to obtain. As transient spiral activity increases random motion, the patterns must fade over time unless the disk also contains a dissipative gas component. Continuing spiral activity alters the structure of the disks in other ways: reducing metallicity gradients and flattening rotation curves are two of the most significant. The overwhelming majority of spirals in galaxies have two- or three-fold rotational symmetry, indicating that the cool, thin disk component is massive. Spirals in simulations of halo-dominated disks instead manifest many arms, and consequently do not capture the expected full spiral-driven evolution. We conclude by identifying areas where further work is needed.

[abridged] We present a 1+1D global magnetically-driven disk accretion model that captures the essence of the MRI-driven accretion, without resorting to 3D global non-ideal MHD simulations. The gas dynamics is assumed to be solely controlled by the MRI and hydrodynamic instabilities. For given stellar and disk parameters, the Shakura-Sunyaev viscosity parameter $\alpha$ is determined self-consistently under the framework of viscously-driven accretion from detailed considerations of the MRI with non-ideal MHD effects (Ohmic resistivity and ambipolar diffusion), accounting for disk heating by stellar irradiation, non-thermal sources of ionization, and dust effects on the ionization chemistry. Additionally, the magnetic field strength is constrained and adopted to maximize the MRI activity. We demonstrate the use of our framework by investigating steady-state MRI-driven accretion in a fiducial protoplanetary disk model around a solar-type star. We find that the equilibrium solution displays no pressure maximum at the dead zone outer edge, except if a sufficient amount of dust particles have accumulated there before the disk reaches a steady-state accretion regime. Furthermore, the steady-state accretion solution describes a disk that displays a spatially extended long-lived inner disk gas reservoir (the dead zone) accreting a few $10^{-9}\, M_{\odot}.\rm{yr}^{-1}$. By conducting a detailed parameter study, we find that the extend to which the MRI can drive efficient accretion is primarily determined by the total disk gas mass, the representative grain size, the vertically-integrated dust-to-gas mass ratio, and the stellar X-ray luminosity. A self-consistent time-dependent coupling between gas, dust, stellar evolution models and our general framework on million-year timescales is required to fully understand the formation of dead zones and their potential to trap dust particles.

We present a new machine learning model for estimating photometric redshifts with improved accuracy for galaxies in Pan-STARRS1 data release 1. Depending on the estimation range of redshifts, this model based on neural networks can handle the difficulty for inferring photometric redshifts. Moreover, to reduce bias induced by the new model's ability to deal with estimation difficulty, it exploits the power of ensemble learning. We extensively examine the mapping between input features and target redshift spaces to which the model is validly applicable to discover the strength and weaknesses of trained model. Because our trained model is well calibrated, our model produces reliable confidence information about objects with non-catastrophic estimation. While our model is highly accurate for most test examples residing in the input space, where training samples are densely populated, its accuracy quickly diminishes for sparse samples and unobserved objects (i.e., unseen samples) in training. We report that out-of-distribution (OOD) samples for our model contain both physically OOD objects (i.e., stars and quasars) and galaxies with observed properties not represented by training data. The code for our model is available at https://github.com/GooLee0123/MBRNN for other uses of the model and retraining the model with different data.

The simplest $\alpha$-attractor $T$ model is given by the potential $V=V_0\tanh^2(\lambda\phi/M_{pl})$. However its generalization to the class of models of the type $V = V_0 \tanh^p (\lambda \phi/M_{pl})$ is difficult to interpret as a model of inflation for most values of $p$. Keeping the basic model, we propose a new generalization, where the final potential is of the form $V = V_0 (1-\sech^p (\lambda \phi/M_{pl}))$, which does not present any of the problems that plague the original generalization, allowing a successful interpretation as a model of inflation for any value of $p$ and, at the same time, providing the potential with a region where reheating can occur for any $p$ (including odd and fractional values) without difficulty. In the cases $p = 1, 2, 4$ we obtain the solutions $r(n_s, N_{ke})$ where $r$ is the tensor-to-scalar ratio, $n_s$ the spectral index and $N_{ke}$ the number of $e$-folds during inflation. We also show how these solutions connect to the $\phi^2$ monomial.

The detection of disk-integrated polarization from Luhman 16A and B in H-band, and subsequent modeling, has been interpreted in the framework of zonal cloud bands on these bodies. Recently, Tan and Showman (2021) investigated three-dimensional atmospheric circulation and cloud structures of brown dwarfs with general circulation models (GCMs), and their simulations yield complex cloud distributions showing some aspects of zonal jets, but also complex vortices that cannot be captured by a simple model. Here we use these 3D GCMs specific to Luhman 16A and B, along with the three-dimensional Monte Carlo radiative transfer code ARTES, to calculate their polarization signals. We adopt the 3D temperature-pressure and cloud profiles from the GCMs as our input atmospheric structures. Our polarization calculations at 1.6 $\mu$m agree well with the measured degree of linear polarization from both Luhman 16 A and B. Our calculations reproduce the measured polarization for both the objects with cloud particle sizes between 0.5-1 \,$\mu$m for Luhman 16 A and 5 \,$\mu$m for Luhman 16 B. We find that the degree of linear polarization can vary on hour-long timescales over the course of a rotation period. We also show that models with azimuthally symmetric band-like cloud geometries, typically used for interpreting polarimetry observations of brown dwarfs, over-predict the polarization signal if the cloud patterns do not include complex vortices within these bands. This exploratory work shows that GCMs are promising for modeling and interpreting polarization signals of brown dwarfs.

We study the \textit{Roman} sensitivity to exoplanets in the Habitable Zone (HZ). The \textit{Roman}~efficiency for detecting habitable planets is maximized for three classes of planetary microlensing events with close caustic topologies. (a) The events with the lens distances of $D_{\rm l} \gtrsim 7$ kpc, the host lens masses of $M_{\rm h}\gtrsim 0.6M_{\odot}$. By assuming Jupiter-mass planets in the HZs, these events have $q\lesssim 0.001$ and $d\gtrsim 0.17$ ($q$ is their mass ratio and $d$ is the projected planet-host distance on the sky plane normalized to the Einstein radius). The events with primary lenses, $M_{\rm h} \lesssim 0.1 M_{\odot}$, while their lens systems are either (b) close to the observer with $D_{\rm l}\lesssim 1$ kpc or (c) close to the Galactic bulge, $D_{\rm l}\gtrsim 7$ kpc. For Jupiter-mass planets in the HZs of the primary lenses, the events in these two classes have $q\gtrsim 0.01$, $d\lesssim 0.04$. The events in the class (a) make larger caustics. By simulating planetary microlensing events detectable by \textit{Roman},~we conclude that the \textit{Roman}~efficiencies for detecting Earth- and Jupiter-mass planets in the Optimistic HZs (OHZs, which is the region between $[0.5,~2]$ AU around a Sun-like star) are $0.01\%$ and $5\%$, respectively. If we assume that one exoplanet orbits each microlens in microlensing events detectable by \textit{Roman}~( i.e., $\sim 27000$ ),~this telescope has the potential to detects $35$ exoplanets with the projected planet-host distances in the OHZs with only one having a mass $\lesssim 10M_{\oplus}$. According to the simulation, $27$ of these exoplanets are actually in the OHZs.

Line intensity mapping is emerging as a novel method that can measure the collective intensity fluctuations of atomic/molecular line emission from distant galaxies. Several observational programs with various wavelengths are ongoing and planned, but there remains a critical problem of line confusion; emission lines originating from galaxies at different distances are confused at an observed wavelength. We devise a generative adversarial network that extracts designated emission line signals from noisy three-dimensional data. Our novel network architecture allows two input data at different wavelengths so that it discerns the co-existence and the correlation of two targeted lines, $\rm H\alpha$ and [OIII]. After being trained with a large number of realistic mock catalogs, the network is able to reconstruct the three-dimensional distribution of emission-line galaxies at $z = 1.3-2.4$. Bright galaxies are identified with a precision of 82 percent, and the cross-correlation coefficients between the true and reconstructed intensity maps are as high as 0.8. Our deep-learning method can be readily applied to data from planned space-borne and ground-based experiments.

In this letter, we report a robust measurement of the morphology and color dependences of the stellar-halo mass relation (SHMR) at the high mass end ($10^{11.3}{\rm M_{\odot}}<M_{\star}<10^{11.7}{\rm M_{\odot}}$) at redshift $z_s\sim0.6$. Applying our method, Photometric objects Around Cosmic webs (PAC), developed in a previous work to CMASS and HSC-SSP observations, we measure the excess surface density ($\bar{n}_2w_p(r_p)$) of satellites around massive central galaxies with different morphologies indicated by S\'ersic index $n$. We find that more compact (larger $n$) central galaxies are surrounded by more satellites. With the abundance matching method, we estimate halo mass for the central galaxies, and find that halo mass is increased monotonically with $n$, solid evidence for a morphology dependence of SHMR. Specifically, our results show that the most compact galaxies ($n>6$) have the halo mass around 5.5 times larger than the disk galaxies ($n<2$). Similarly, using the rest-frame $u-r$ color, we find that red galaxies reside in halos $2.6$ times larger than those hosting blue galaxies.

Analytical tools are valuable to study gravitational collapse. However, solutions are hard to find due to the highly non-linear nature. Only a few simple but powerful tools exist so far. Two examples are the spherical collapse model (SCM) and stable clustering hypothesis (SCH). We present a new tool based on an elementary step of mass cascade, i.e. a two-body collapse model. TBCM plays the same role as harmonic oscillator in dynamics and can be fundamental to understand structure evolution. For convenience, TBCM is formulated for gravity with any exponent $n$ in a static background with fixed damping. The competition between gravity, expanding background (or damping), and angular momentum classifies two-body collapse into: 1) free fall collapse for weak angular momentum, where free fall time is greater if same system starts to collapse at earlier time; 2) equilibrium collapse for weak damping that persists longer in time, where perturbative solutions lead to power-law evolution of system energy and momentum. Two critical values $\beta_{s1}=1$ and $\beta_{s2}=1/3\pi$ are identified that quantifies the competition between damping and gravity. Value $\beta_{s2}$ only exists for discrete values of n=(2-6m)/(1+3m)= -1,-10/7...for integer m. Critical density ratio ($18\pi^2$) is obtained for $n$=-1 that is consistent with SCM. TBCM predicts angular velocity $\propto Hr^{-3/2}$ with r. The isothermal density is a result of infinitesimal halo lifetime. TBCM is able to demonstrate SCP, i.e. mean pairwise velocity (first moment) $\langle\Delta u\rangle=-Hr$. A generalized SCH is developed for higher order moments $\langle\Delta u^{2m+1}\rangle=-(2m+1)\langle\Delta u^{2m}\rangle Hr$. Energy evolution in TBCM is independent of mass and energy equipartition does not apply. TBCM can be considered as a non-radial SCM. Both models predict same critical ratio, while TBCM contains much richer information.

Jupiter's tenuous dust ring system is embedded in the planet's inner magnetosphere, and - among other structures - contains a very tenuous protrusion called the Thebe extension. In an attempt to explain the existence of this swath of particles beyond Thebe's orbit, Hamilton and Kr\"uger (2008) proposed that the dust particle motion is driven by a shadow resonance caused by variable dust charging on the day and night side of Jupiter. However, the model by Divine (1983) together with recent observations by the Juno spacecraft indicates a warm and rather dense inner magnetosphere of Jupiter which implies that the mechanism of the shadow resonance does not work. Instead we argue that dust grains ejected from Thebe due to micrometeoroid bombardment become the source of dust in the Thebe extension. We show that large (grain radii of a few micrometers up to multi-micrometers) charged dust grains having significant initial velocities oscillate in the Thebe extension. Smaller charged grains (with sub-micrometer radii) ejected from Thebe do not spend much time in the Thebe extension and migrate into the Thebe ring. At the same time, if such grains are ejected from larger dust grains in the Thebe extension due to fragmentation, they continue to oscillate within the Thebe extension for years. We argue that fragmentation of large dust grains in the Thebe extension could be the main source of sub-micrometer grains detected in the Thebe extension.

The CoPhyLab (Cometary Physics Laboratory) project is designed to study the physics of comets through a series of earth-based experiments. For these experiments, a dust analogue was created with physical properties comparable to those of the non-volatile dust found on comets. This "CoPhyLab dust" is planned to be mixed with water and CO$_2$ ice and placed under cometary conditions in vacuum chambers to study the physical processes taking place on the nuclei of comets. In order to develop this dust analogue, we mixed two components representative for the non-volatile materials present in cometary nuclei. We chose silica dust as representative for the mineral phase and charcoal for the organic phase, which also acts as a darkening agent. In this paper, we provide an overview of known cometary analogues before presenting measurements of eight physical properties of different mixtures of the two materials and a comparison of these measurements with known cometary values. The physical properties of interest are: particle size, density, gas permeability, spectrophotometry, mechanical, thermal and electrical properties. We found that the analogue dust that matches the highest number of physical properties of cometary materials consists of a mixture of either 60\%/40\% or 70\%/30\% of silica dust/charcoal by mass. These best-fit dust analogue will be used in future CoPhyLab experiments.

The repeating fast radio burst FRB 20200120E is located in a globular cluster belonging to the nearby M81 galaxy. Its small distance (3.6 Mpc) and accurate localization make it an interesting target to search for bursting activity at high energies. From November 2003 to September 2021, the INTEGRAL satellite has obtained an exposure time of 18 Ms on the M81 sky region. We used these data to search for hard X-ray bursts from FRB 20200120E using the IBIS/ISGRI instrument, without finding any significant candidate, down to an average fluence limit of $\sim10^{-8}$ erg cm$^{-2}$ (20-200 keV). The corresponding limit on the isotropic luminosity for a burst of duration $\Delta t$ is $\sim10^{45} \left ( \frac{10~ms}{\Delta t} \right )$ erg s$^{-1}$, the deepest limit obtained for an extragalactic FRB in the hard X-ray range. This rules out the emission of powerful flares at a rate higher than 0.1 yr$^{-1}$ that could be expected in models invoking young hyper-active magnetars.

The observed atmospheric spectrum of exoplanets and brown dwarfs depends critically on the presence and distribution of atmospheric condensates. The Ackerman & Marley (2001) methodology for predicting the vertical distribution of condensate particles is widely used to study cloudy atmospheres and has recently been implemented in an open-source python package virga. The model relies upon input parameter $f_{\text{sed}}$, the sedimentation efficiency, which until now has been held constant. The relative simplicity of this model renders it useful for retrieval studies due to its rapidly attainable solutions. However, comparisons with more complex microphysical models such as CARMA have highlighted inconsistencies between the two approaches, namely that the cloud parameters needed for radiative transfer produced by virga are dissimilar to those produced by CARMA. To address these discrepancies, we have extended the original Ackerman and Marley methodology in virga to allow for non-constant $f_{\text{sed}}$ values, in particular those that vary with altitude. We discuss one such parameterization and compare the cloud mass mixing ratio produced by virga with constant and variable $f_{\text{sed}}$ profiles to that produced by CARMA. We find that the variable $f_{\text{sed}}$ formulation better captures the profile produced by CARMA with heterogeneous nucleation, yet performs comparatively to constant $f_{\text{sed}}$ for homogeneous nucleation. In general, virga has the capacity to handle any $f_{\text{sed}}$ with an explicit anti-derivative, permitting a plethora of alternative cloud profiles that are otherwise unattainable by constant $f_{\text{sed}}$ values. The ensuing flexibility has the potential to better agree with increasingly complex models and observed data.

DHOST inflation models where deviations from a pure de Sitter background are induced by an axion-like potential can lead to large non-Gaussianities. We investigate the nature of non-Gaussianities in these models and compare to the results given by the Planck experiment. The overlap between the DHOST non-Gaussianities and the equilateral, orthogonal and local templates can be rendered arbitrarily small. On the other hand, this does not preclude DHOST models from showing large non-Gaussianities as exemplified by their reduced bispectrum. As a result, they could be probed by future experiments and also by a more thorough analysis of the existing Planck data.

We present results of the analysis of light and radial velocity (RV) curves of eight detached eclipsing binaries observed by the All-Sky Automated Survey, which we have followed up with high-resolution spectroscopy, and were later observed by the $Kepler$ satellite as part of the $K2$ mission. The RV measurements came from spectra obtained with OAO-188/HIDES, MPG-2.2m/FEROS, SMARTS 1.5m/CHIRON, Euler/CORALIE, ESO-3.6m/HARPS, and OHP-1.93/ELODIE instruments. The $K2$ time-series photometry was analyzed with the JKTEBOP code, with out-of-eclipse modulations of different origin taken into account. Individual component spectra were retrieved with the FD3 code, and analyzed with the code iSpec in order to determine effective temperatures and metallicities. Absolute values of masses, radii, and other stellar parameters are calculated, as well as ages, found through isochrone fitting. For five systems such analysis has been done for the first time. The presented sample consists of a variety of stars, from low-mass dwarfs, through G and F-type Main Sequence objects, to evolved active sub-giants, one of which is found to be crossing the Hertzsprung gap. One target may contain a $\gamma$ Dor-type pulsator, two more are parts of higher-order multiples, and spectra of their tertiaries were also retrieved and used to constrain the properties of these systems.

In this paper, we present an application of the Shannon entropy in the case of the planar (non-restricted) four-body problem. Specifically, the Kepler-60 extrasolar system is being investigated with a primary interest in the resonant configuration of the planets that exhibit a chain of mean-motion commensurabilities with the ratios 5:4:3. In the dynamical maps provided, the Shannon entropy is utilized to explore the general structure of the phase space, while, based on the time evolution of the entropy, we determine also the extent and rate of the chaotic diffusion as well as the characteristic times of stability for the planets. Two cases are considered: (i) the pure Laplace resonance when the critical angles of the 2-body resonances circulate and that of the 3-body resonance librates; and (ii) the chain of two 2-body resonances when all the critical angles librate. Our results suggest that case (ii) is the more favourable configuration but we state too that, in either case, the relevant resonance plays an important role to stabilize the system. The derived stability times are no shorter than $10^8$ yrs in the central parts of the resonances.

Detection of black holes (BHs) with detached luminous companions (LCs) can be instrumental in connecting the BH properties with their progenitors' since the latter can be inferred from the observable properties of the LC. Past studies showed the promise of Gaia astrometry in detecting BH-LC binaries. We build upon these studies by: 1) initialising the zero-age binary properties based on realistic, metallicity-dependent star-formation history in the Milky Way (MW), 2) evolving these binaries to current epoch to generate realistic MW populations of BH-LC binaries, 3) distributing these binaries in the MW preserving the complex age-metallicity-Galactic position correlations, 4) accounting for extinction and reddening using three-dimensional dust maps, 5) examining the extended Gaia mission's ability to resolve BH-LC binaries. We restrict ourselves to detached BH-LC binaries with orbital period <10 yr such that Gaia can observe at least one full orbit. We find: 1) the extended Gaia mission can astrometrically resolve 30-300 detached BH-LC binaries depending on our assumptions of supernova physics and astrometric detection threshold; 2) Gaia's astrometry alone can indicate BH candidates for 10-100 BH-LC binaries by constraining the dark primary mass >3 Msun; 3) distributions of observables including orbital periods, eccentricities, and component masses are sensitive to the adopted binary evolution model, hence can directly inform binary evolution models. Finally, we comment on the potential to further characterise these BH binaries through radial velocity measurements and observation of X-ray counterparts.

A detailed investigation of the magnetic properties of young Sun-like stars can provide valuable information on our Sun's magnetic past and its impact on the early Earth. We determine the properties of the moderately rotating young Sun-like star kappa Ceti's magnetic and activity cycles using 50 years of chromospheric activity data and six epochs of spectropolarimetric observations. The chromospheric activity was determined by measuring the flux in the Ca II H and K lines. A generalised Lomb-Scargle periodogram and a wavelet decomposition were used on the chromospheric activity data to establish the associated periodicities. The vector magnetic field of the star was reconstructed using the technique of Zeeman Doppler imaging on the spectropolarimetric observations. Our period analysis algorithms detect a 3.1 year chromospheric cycle in addition to the star's well-known ~6 year cycle period. Although the two cycle periods have an approximate 1:2 ratio, they exhibit an unusual temporal evolution. Additionally, the spectropolarimetric data analysis shows polarity reversals of the star's large-scale magnetic field, suggesting a ~10 year magnetic or Hale cycle. The unusual evolution of the star's chromospheric cycles and their lack of a direct correlation with the magnetic cycle establishes kappa Ceti as a curious young Sun. Such complex evolution of magnetic activity could be synonymous with moderately active young Suns, which is an evolutionary path that our own Sun could have taken.

There are different environments in the interstellar medium (ISM), depending on the density, temperature and chemical composition. Among them, molecular clouds, often referred to as the cradle of stars, are paradigmatic environments relative to the chemical diversity and complexity in space. Indeed, there, radio to far-infrared observations revealed the presence of several molecules in the gas phase, while near-infrared spectroscopy detected the existence of submicron sized dust grains covered by H2O -dominated ice mantles. The interaction between gas-phase species and the surfaces of water ices is measured by the binding energy (BE), a crucial parameter in astrochemical modelling. In this work, the BEs of a set of sulphur-containing species on water ice mantles have been computed by adopting a periodic ab initio approach using a crystalline surface model. The Density Functional Theory (DFT)-based B3LYP-D3(BJ) functional was used for the prediction of the structures and energetics. DFT BEs were refined by adopting an ONIOM-like procedure to estimate them at CCSD(T) level toward complete basis set extrapolation, in which a very good correlation between values has been found. Moreover, we show that geometry optimization with the computationally cheaper HF-3c method followed by single point energy calculations at DFT to compute the BEs is a suitable cost-effective recipe to arrive at BE values of the same quality as those computed at full DFT level. Finally, computed data were compared with the available literature data.

NectarCAM is a camera for the medium-sized telescopes of the Cherenkov Telescope Array (CTA), which covers the energy range of 100 GeV to 30 TeV. The camera is equipped with 265 focal plane modules (FPMs). Each FPM comprises 7 pixels, each consisting of a photo-multiplier tube, a preamplifier, an independently controlled power supply, and a common control system. We developed a dedicated test bench to validate and qualify the industrial FPM production and to measure the performance of each FPM in a dark room before its integration in the camera. We report the measured performance of 61 FPM prototypes obtained with our experimental setup. We demonstrate that the gains of the photo multiplier tubes are stable and that pulse widths, transit time spreads, afterpulse rates and charge resolutions are within the specifications for NectarCAM.

We study how stellar rotation curves (RCs) of galaxies are correlated on average with morphology and stellar mass ($M_\mathrm{star}$) using the final release of Sloan Digital Sky Survey IV MaNGA data. We use the visually assigned $T$-types for the morphology indicator, and adopt a functional form for the RC that can model non-flat RCs at large radii. We discover that within the radial coverage of the MaNGA data, the popularly known flat rotation curve at large radii applies only to the particular classes of galaxies, i.e., massive late types ($T$-type $\geq1$, $M_\mathrm{star}\gtrsim10^{10.8}M_\odot$) and S0 types ($T$-type$=-1$ or $0$, $M_\mathrm{star}\gtrsim10^{10.0}M_\odot$). The RC of late-type galaxies at large radii rises more steeply as $M_\mathrm{star}$ decreases, and its slope increases to about $+9$ km s$^{-1}$kpc$^{-1}$ at $M_\mathrm{star}\approx10^{9.7}M_\odot$. By contrast, elliptical galaxies ($T$-type $\le-2$) have descending RCs at large radii. Their slope becomes more negative as $M_\mathrm{star}$ decreases, and reaches as negative as $-15$ km s$^{-1}$kpc$^{-1}$ at $M_\mathrm{star}\approx10^{10.2}M_\odot$. We also find that the inner slope of the RC is highest for elliptical galaxies with $M_\mathrm{star}\approx10^{10.5}M_\odot$, and decreases as $T$-type increases or $M_\mathrm{star}$ changes away from $10^{10.5}M_\odot$. The velocity at the turnover radius $R_t$ is higher for higher $M_\mathrm{star}$, and $R_t$ is larger for higher $M_\mathrm{star}$ and later $T$-types. We show that the inner slope of the RC is coupled with the central surface stellar mass density, which implies that the gravitational potential of central regions of galaxies is dominated by baryonic matter. With the aid of simple models for matter distribution, we discuss what determines the shapes of RCs.

Extremely powerful magnetic fields are contained inside neutron stars. Their effect is to deform the shape of the star, leading to the emission of continuous gravitational waves. The magnetic deformation of neutron stars depends on the details of their magnetic field, that is its geometry and strength. Moreover, it depends on their composition, described by the equation of state. Unfortunately, both the configuration of the magnetic field and the equation of state of neutron stars are unkown, and assessing the detectability of continuous gravitational waves from neutron stars suffers from these uncertainties. Using our recent results relating the magnetic deformation of a neutron star to its mass and radius, and considering the Galactic pulsar population, we assess the detectability of continuous gravitational waves from pulsars in the Galaxy - described by realistic equations of state currently allowed by observational and nuclear physics constraints - by gravitational waves detectors.

Nucleosynthetic Fe isotopic anomalies in meteorites may be used to reconstruct the early dynamical evolution of the solar system and to identify the origin and nature of the material that built planets. Using high-precision iron isotopic data of 23 iron meteorites from nine major chemical groups we show that all iron meteorites show the same fundamental dichotomy between non-carbonaceous (NC) and a carbonaceous (CC) meteorites previously observed for other elements. The Fe isotopic anomalies are predominantly produced by variation in 54Fe, where all CC iron meteorites are characterized by an excess in 54Fe relative to NC iron meteorites. This excess in 54Fe is accompanied by an excess in 58Ni observed in the same CC meteorite groups. Together, these overabundances of 54Fe and 58Ni are produced by nuclear statistical equilibrium either in type Ia supernovae or in the Si/S shell of core-collapse supernovae. The new Fe isotopic data reveal that Earth's mantle plots on or close to correlations defined by Fe, Mo, and Ru isotopic anomalies in iron meteorites, indicating that throughout Earth's accretion, the isotopic composition of its building blocks did not drastically change. While Earth's mantle has a similar Fe isotopic composition to CI chondrites, the latter are clearly distinct from Earth's mantle for other elements (e.g., Cr and Ni) whose delivery to Earth coincided with Fe. The fact that CI chondrites exhibit large Cr and Ni isotopic anomalies relative to Earth's mantle, therefore, demonstrates that CI chondrites are unlikely to have contributed significant Fe to Earth.

Observing gravitationally lensed objects in the time domain is difficult, and well-observed time-varying sources are rare. Lensed gamma-ray bursts (GRBs) offer improved timing precision to this class of objects complementing observations of quasars and supernovae. The rate of lensed GRBs is highly uncertain, approximately 1 in 1000. The Gamma-ray Burst Monitor (GBM) onboard the Fermi Gamma-ray Space Telescope has observed more than 3000 GRBs making it an ideal instrument to uncover lensed bursts. Here we present observations of GRB 210812A showing two emission episodes, separated by 33.3 s, and with flux ratio of about 4.5. An exhaustive temporal and spectral analysis shows that the two emission episodes have the same pulse and spectral shape, which poses challenges to GRB models. We report multiple lines of evidence for a gravitational lens origin. In particular, modeling the lightcurve using nested sampling we uncover strong evidence in favor of the lensing scenario. Assuming a point mass lens, the mass of the lensing object is about 1 million solar masses. High-resolution radio imaging is needed for future lens candidates to derive tighter constraints.

The morphology of the Milky Way is still a matter of debate. In order to shed light on uncertainties surrounding the structure of the Galaxy, in this paper, we study the imprint of spiral arms on the distribution and properties of its molecular gas. To do so, we take full advantage of the SEDIGISM survey that observed a large area of the inner Galaxy in the $^{13}$CO(2-1) line at an angular resolution of 28". We analyse the influences of the spiral arms by considering the features of the molecular gas emission as a whole across the longitude-velocity map built from the full survey. Additionally, we examine the properties of the molecular clouds in the spiral arms compared to the properties of their counterparts in the inter-arm regions. Through flux and luminosity probability distribution functions, we find that the molecular gas emission associated with the spiral arms does not differ significantly from the emission between the arms. On average, spiral arms show masses per unit length of $\sim10^5-10^6$ M$_{\odot} $kpc$^{-1}$. This is similar to values inferred from data sets in which emission distributions were segmented into molecular clouds. By examining the cloud distribution across the Galactic plane, we infer that the molecular mass in the spiral arms is a factor of 1.5 higher than that of the inter-arm medium, similar to what is found for other spiral galaxies in the local Universe. We observe that only the distributions of cloud mass surface densities and aspect ratio in the spiral arms show significant differences compared to those of the inter-arm medium; other observed differences appear instead to be driven by a distance bias. By comparing our results with simulations and observations of nearby galaxies, we conclude that the measured quantities would classify the Milky Way as a flocculent spiral galaxy, rather than as a grand-design one.

We present measurements of the Minkowski functionals extracted from the SDSS-III BOSS catalogs. After defining the Minkowski functionals, we describe how an unbiased reconstruction of these statistics can be obtained from a field with masked regions and survey boundaries, validating our methodology with Gaussian random fields and mock galaxy snapshot data. From the BOSS galaxy data we generate a set of four density fields in three dimensions corresponding to the northern and southern skies of LOWZ and CMASS catalogs, smoothing over large scales such that the field is perturbatively non-Gaussian. We extract the Minkowski functionals from each data set separately, and measure their shapes and amplitudes by fitting a Hermite polynomial expansion. For the shape parameter of the Minkowski functional curves $a_0$, that is related to the bispectrum of the field, we find that the LOWZ-South data presents a systematically lower value of $a_0 = -0.080 \pm 0.040$ than its northern sky counterpart $a_0 = 0.032 \pm 0.024$. Although the significance of this discrepancy is low, it potentially indicates some systematics in the data or that the matter density field exhibits anisotropy at low redshift. By assuming a standard isotropic flat $\Lambda$CDM cosmology, the amplitudes of Minkowski functionals from the combination of northern and southern sky data give the constraints $\Omega_{\rm c} h^2 n_{\rm s} = 0.110 \pm 0.006$ and $0.111 \pm 0.008$ for CMASS and LOWZ, respectively, which is in agreement with the Planck $\Lambda$CDM best-fit $\Omega_{\rm c}h^{2} n_{\rm s} = 0.116 \pm 0.001$.

Our understanding of the formation and evolution of the corona and the heliosphere is linked to our capability of properly interpreting the data from remote sensing and in-situ observations. In this respect, being able to correctly connect in-situ observations with their source regions on the Sun is the key for solving this problem. In this work we aim at testing a diagnostics method for this connectivity. This paper makes use of a coronal jet observed on 2010 August 2nd in active region 11092 as a test for our connectivity method. This combines solar EUV and in-situ data together with magnetic field extrapolation, large scale MHD modeling and FIP (First Ionization Potential) bias modeling to provide a global picture from the source region of the jet to its possible signatures at 1AU. Our data analysis reveals the presence of outflow areas near the jet which are within open magnetic flux regions and which present FIP bias consistent with the FIP model results. In our picture, one of these open areas is the candidate jet source. Using a back-mapping technique we identified the arrival time of this solar plasma at the ACE spacecraft. The in-situ data show signatures of changes in the plasma and magnetic field parameters, with FIP bias consistent with the possible passage of the jet material. Our results highlight the importance of remote sensing and in-situ coordinated observations as a key to solve the connectivity problem. We discuss our results in view of the recent Solar Orbiter launch which is currently providing such unique data.

The InSight mission, featuring continuous high-frequency high-sensitivity pressure measurements, is in ideal position to study the active atmospheric turbulence of Mars. Data acquired during 1.25 Martian year allows us to study the seasonal evolution of turbulence and its diurnal cycle. We investigate vortices (abrupt pressure drops), local turbulence (frequency range 0.01-2 Hz) and non-local turbulence often caused by convection cells and plumes (frequency range 0.002-0.01 Hz). Contrary to non-local turbulence, local turbulence is strongly sensitive at all local times and seasons to the ambient wind. We report many remarkable events with the arrival of northern autumn at the InSight landing site: a spectacular burst of daytime vortices, the appearance of nighttime vortices, and the development of nighttime local turbulence as intense as its daytime counterpart. Nighttime turbulence at this dusty season appears as a result of the combination of a stronger low-level jet, producing shear-driven turbulence, and a weaker stability.

We explore the possibility of detecting Baryon Acoustic Oscillations (BAO) solely from gravitational wave observations of binary neutron star mergers with third generation (3G) gravitational wave (GW) detectors like Cosmic Explorer and the Einstein Telescope. These measurements would provide a new independent probe of cosmology. The detection of the BAO peak with current generation GW detectors (solely from GW observations) is not possible because i) unlike galaxies, the GW mergers are poorly localized and ii) there are not enough merger events to probe the BAO length scale. With the 3G GW detector network, it is possible to observe $\sim \mathcal{O}(1000)$ binary neutron star mergers per year localized well within one square degree in the sky for redshift $z \leq 0.3$. We show that 3G observatories will enable precision measurements of the BAO feature in the large-scale two-point correlation function; the effect of BAO can be independently detected at different reshifts, with a log-evidence ratio of $\sim$ 23, 17, or 3 favouring a model with a BAO peak at redshift of 0.2, 0.25, or 0.3, respectively, using a redshift bin corresponding to a shell of thickness $~150 h^{-1}$ Mpc.

Apep is the brightest and most luminous non-thermal colliding-wind binary by over an order of magnitude. It has been suggested from infrared observations that one of the Wolf-Rayet stars in Apep is launching an anisotropic wind. Here we present radio observations of Apep from 0.2 to 20 GHz taken over 33 years. The spectrum reveals an extremely steep turnover in the flux density at low frequencies, where the flux density decreases by two orders of magnitude over only 325 MHz of bandwidth. This exponential decline is best described by free-free absorption, with a turnover frequency at 0.54 $\pm$ 0.01 GHz. Above the turnover, the spectrum is well described by a power-law and a high-frequency cut-off likely caused by inverse-Compton cooling. The lightcurve of Apep shows significant variation over the observing period, with Apep brightening by over 50 mJy in a span of 25 years at 1.4 GHz. Models that assume spherical winds do not replicate all of the structure evident in the radio lightcurve. We derived a model that allows one of the winds in the system to be anisotropic. This anisotropic model recovers most of the structure of the lightcurve and is a significantly better statistical fit to the data than the spherical wind model. We suggest such a result is independent support that one of the Wolf-Rayet stars in Apep is launching an anisotropic wind. If the anisotropic wind model is correct, we predict a ~25% decrease of the 1.4 GHz flux density of Apep over the next five years.

Plasma outflow from a gravitational potential well with cosmic rays and self-excited Alfv\'en waves with cooling and wave damping is studied in the hydrodynamics regime. We study outflows in the presence of cosmic ray and Alfv\'en waves including the effect of cooling and wave damping. We seek physically allowable steady-state subsonic-supersonic transonic solutions. We adopted a multi-fluid hydrodynamical model for the cosmic ray plasma system. Thermal plasma, cosmic rays, and self-excited Alfv\'en waves are treated as fluids. Interactions such as cosmic-ray streaming instability, cooling, and wave damping were fully taken into account. We considered one-dimensional geometry and explored steady-state solutions. The model is reduced to a set of ordinary differential equations, which we solved for subsonic-supersonic transonic solutions with given boundary conditions at the base of the gravitational potential well. We find that physically allowable subsonic-supersonic transonic solutions exist for a wide range of parameters. We studied the three-fluid system (considering only forward-propagating Alfv\'en waves) in detail. We examined the cases with and without cosmic ray diffusion separately. Comparisons of solutions with and without cooling and with and without wave damping for the same set of boundary conditions (on density, pressures of thermal gas, cosmic rays and waves) are presented. We also present the interesting case of a four-fluid system (both forward- and backward-propagating Alfv\'en waves are included), highlighting the intriguing relation between different components.

The Square Kilometre Array (SKA) is a planned radio interferometer of unprecedented scale that will revolutionize low-frequency radio astronomy when completed. In particular, one of its core science drivers is the systematic mapping of the Cosmic Dawn and Epoch of Reionization, which mark the birth of the first stars and galaxies in the Universe and their subsequent ionization of primordial intergalactic hydrogen, respectively. The SKA will offer the most sensitive view of these poorly understood epochs using the redshifted 21 cm hyperfine signal from intergalactic hydrogen. However, significant technical challenges stand in the way of realizing this scientific promise. These mainly involve the mitigation of systematics coming from astrophysical foregrounds, terrestrial radio interference, and the instrumental response. The Low Frequency Array, the Murchison Widefield Array and the Hydrogen Epoch of Reionization Array are SKA pathfinder experiments that have developed a variety of strategies for addressing these challenges, each with unique characteristics that stem largely from their different instrumental designs. We outline these various directions, highlighting key differences and synergies, and discuss how these relate to the future of low-frequency intensity mapping with the SKA. We also briefly summarize the challenges associated with modeling the 21 cm signal and discuss the methodologies being proposed for inferring constraints on astrophysical models.

We present GLADE+, an extended version of the GLADE galaxy catalogue introduced in our previous paper for multimessenger searches with advanced gravitational-wave detectors. GLADE+ combines data from six separate but not independent astronomical catalogues: the GWGC, 2MPZ, 2MASS XSC, HyperLEDA, and WISExSCOSPZ galaxy catalogues, and the SDSS-DR16Q quasar catalogue. To allow corrections of CMB-frame redshifts for peculiar motions, we calculated peculiar velocities along with their standard deviations of all galaxies having $B$-band magnitude data within redshift $z=0.05$ using the "Bayesian Origin Reconstruction from Galaxies" formalism. GLADE+ is complete up to luminosity distance $d_L=47^{+4}_{-2}$ Mpc in terms of the cumulative $B$-band luminosity of galaxies, and contains all of the brightest galaxies giving half of the total $B$-band luminosity up to $d_L\simeq 250$ Mpc. We include estimations of stellar masses and individual binary neutron star merger rates for galaxies with $W1$ magnitudes in GLADE+. These parameters can help in ranking galaxies in a given gravitational wave localization volume in terms of their likelihood of being hosts, thereby possibly reducing the number of pointings and total integration time needed to find the electromagnetic counterpart.

We report NanoSIMS Si and Mg-Al isotopic data (and C, N, and Ti isotopic data when available) for 85 submicron- to micron-sized presolar SiC grains from the CM2 Murchison meteorite, including 60 MS, 8 AB1, 8 X, 7 AB2, and 2 Y grains. The MS and Y grain data demonstrate that (1) C and N contamination mainly appears as surface contamination, and sufficient presputtering is needed to expose a clean grain surface for obtaining intrinsic C and N signals, and (2) Mg and Al contamination appears as adjacent grains and rims, and high-resolution imaging and the choice of small regions of interest during data reduction together are effective in suppressing the contamination. Our results strongly indicate that previous studies on presolar SiC grains could have sampled differing degrees of contamination for C, N, Mg, and Al. Compared to the literature data, our new MS and Y grains are in better agreement with carbon star observations for both the C and N isotopic ratios. By comparing our new, tighter distributions of 12C/13C, 14N/15N, and initial 26Al/27Al ratios for MS and Y grains with FRUITY AGB stellar models, we provide more stringent constraints on the occurrence of cool bottom processing and the production of 26Al in N-type carbon stars, classical asymptotic giant branch stars.

We present a new derivation of relativistic second-order dissipative hydrodynamics for quantum systems using Zubarev's non-equilibrium statistical-operator formalism. This is achieved by a systematic expansion of the energy-momentum tensor and the charge current to second order in deviations from equilibrium. As a concrete example, we obtain the relaxation equations for the shear-stress tensor, the bulk-viscous pressure, and the charge-diffusion currents required to close the set of equations of motion for relativistic second-order dissipative hydrodynamics. We also identify new transport coefficients which describe the relaxation of dissipative processes to second order and express them in terms of equilibrium correlation functions, thus establishing new Kubo-type formulas for second-order transport coefficients.

We study the possibility that, after inflation, the inflaton reaches thermal equilibrium with the Standard Model thermal bath and eventually freezes-out in the non-relativistic regime. When the inflaton decay is the sole source of (non-thermal) dark matter, its relic density is automatically suppressed. We delineate parameter space leading to the correct dark matter abundance. The model allows for a significant Higgs-inflaton coupling which may lead to invisible Higgs decay into inflaton pairs at the LHC.

It has been recently argued \cite{Henry_2014,Akshaya_2018,2019MNRAS.489.1120A} that there is a strong component of the diffuse far-ultraviolet (FUV) background which is hard to explain by conventional physics in terms of the dust-scattered starlight. We propose that this excess in FUV radiation might be result of the dark matter annihilation events within the so-called axion quark nugget (AQN) dark matter model, which was originally invented for completely different purpose to explain the observed similarity between the dark and the visible components in the Universe, i.e. $\Omega_{\rm DM}\sim \Omega_{\rm visible}$. We support this proposal by demonstrating that intensity and the spectral features of the AQN induced emissions are consistent with the corresponding characteristics of the observed excess of the FUV radiation.

We construct a relativistic model for bulk viscosity and heat conduction in a superfluid. Building on the principles of Unified Extended Irreversible Thermodynamics, the model is derived from Carter's multifluid approach for a theory with three currents, where the quasi-particle current is an independent hydrodynamic degree of freedom. For small deviations from local thermodynamic equilibrium, the model reduces to an extension of the Israel-Stewart theory to superfluid systems. It can, therefore, be made hyperbolic, causal and stable if the microscopic input is accurate. The non-dissipative limit of the model is the relativistic two-fluid model of Carter, Khalatnikov and Gusakov. The Newtonian limit of the model is an Extended-Irreversible-Thermodynamic extension of Landau's two-fluid model. The model predicts the existence of four bulk viscosity coefficients and accounts for their microscopic origin, providing their exact formulas in terms of the quasi-particle creation rate. Furthermore, when fast oscillations of small amplitude around the equilibrium are considered, the relaxation-time term in the telegraph-type equations for the bulk viscosities accounts directly for their expected dependence on the frequency.

A primordial black hole in the last stages of evaporation and located in the local neighborhood can produce a detectable signal in gamma ray and neutrino telescopes. We re-evaluate the expected gamma ray and neutrino fluxes from these transient point events and discuss the consequences for existing constraints. For gamma rays we improve the current bounds by a factor of few, while for neutrinos we obtain significantly different results than the existing literature. The capability and advantages of neutrino telescopes in the search for primordial black holes is discussed thoroughly. The correlations of gamma ray and neutrino energy and time profiles will be promoted as a powerful tool in identifying the primordial black holes, in case of detection.

The shell evolution of neutron-rich nuclei with temperature is studied in a beyond-mean-field framework rooted in the meson-nucleon Lagrangian. The temperature-dependent Dyson equation with the dynamical kernel taking into account the particle-vibration coupling (PVC) is solved for the fermionic propagators in the basis of the thermal relativistic mean-field Dirac spinors. The calculations are performed for $^{68-78}$Ni in a broad range of temperatures $0 \leq T \leq 4$ MeV. The special focus is put on the fragmentation pattern of the single-particle states, which is further investigated within toy models in strongly truncated model spaces. Such models allow for quantifying the sensitivity of the fragmentation to the phonon frequencies, the PVC strength and to the mean-field level density. The model studies provide insights into the temperature evolution of the PVC mechanism in real nuclear systems under the conditions which may occur in astrophysical environments.

Typical materials for optical Kinetic Inductance Detetectors (KIDs) are metals with a natural absorption of 30-50% in the visible and near-infrared. To reach high absorption efficiencies (90-100%) the KID must be embedded in an optical stack. We show an optical stack design for a 60 nm TiN film. The optical stack is modeled as sections of transmission lines, where the parameters for each section are related to the optical properties of each layer. We derive the complex permittivity of the TiN film from a spectral ellipsometry measurement. The designed optical stack is optimised for broadband absorption and consists of, from top (illumination side) to bottom: 85 nm SiOx, 60 nm TiN, 23 nm of SiOx, and a 100 nm thick Al mirror. We show the modeled absorption and reflection of this stack, which has >80% absorption from 400 nm to 1550 nm and near-unity absorption for 500 nm to 800 nm. We measure transmission and reflection of this stack with a commercial spectrophotometer. The results are in good agreement with the model.

Axion-like particles (ALPs) are pseudoscalar bosons predicted by string theory. The ALPs have a shallower potential than a quadratic one, which induces the instability and can form the solitonic object called oscillon/I-ball. Although the lifetime of oscillons can be very long for some type of potentials, they finally decay until the present. We perform the numerical lattice simulations to investigate the decay process of oscillons and evaluate the averaged momentum of ALPs emitted from the oscillon decay. It is found that, if oscillons decay in the early universe, the free-streaming length of ALPs becomes too long to explain the small-scale observations of the matter power spectrum. We show that oscillons with long lifetimes can change the density fluctuations on small scales, which leads to stringent constraints on the ALP mass and the oscillon lifetime.

We report on the search for weakly interacting massive particle (WIMP) dark matter candidates in the galactic halo that interact with sodium and iodine nuclei in the COSINE-100 experiment and produce energetic electrons that accompany recoil nuclei via the the Migdal effect. The WIMP mass sensitivity of previous COSINE-100 searches that relied on the detection of ionization signals produced by target nuclei recoiling from elastic WIMP-nucleus scattering was restricted to WIMP masses above $\sim$5 GeV/$c^2$ by the detectors' 1 keVee energy-electron-equivalent threshold. The search reported here looks for recoil signals enhanced by the Migdal electrons that are ejected during the scattering process. This is particularly effective for the detection of low-mass WIMP scattering from the crystals' sodium nuclei in which a relatively larger fraction of the WIMP's energy is transferred to the nucleus recoil energy and the excitation of its orbital electrons. In this analysis, the low-mass WIMP search window of the COSINE-100 experiment is extended to WIMP mass down to 200\,MeV/$c^2$. The low-mass WIMP sensitivity will be further improved by lowering the analysis threshold based on a multivariable analysis technique. We consider the influence of these improvements and recent developments in detector performance to re-evaluate sensitivities for the future COSINE-200 experiment. With a 0.2 keVee analysis threshold and high light-yield NaI(Tl) detectors (22 photoelectrons/keVee), the COSINE-200 experiment can explore low-mass WIMPs down to 20 MeV/$c^2$ and probe previously unexplored regions of parameter space.

We examine in greater detail the recent proposal that time is the conjugate of the constants of nature. Fundamentally distinct times are associated with different constants and we should select the one related to the constant dominating the dynamics in each region or epoch. We show in detail how in regions dominated by a single constant the Hamiltonian constraint can be reframed as a Schrodinger equation in the corresponding time, solved in the connection representation by outgoing-only monochromatic plane waves moving in a "space" that generalizes the Chern-Simons functional. We pay special attention to the issues of unitarity and the measure employed for the inner product. Normalizable superpositions can be built, including solitons, "light-rays" and coherent/squeezed states saturating a Heisenberg uncertainty relation between constants and their times. A healthy classical limit is obtained for factorizable coherent states, with classical cosmology seen through the prism of the connection (the comoving Hubble length) rather than the more conventional expansion factor (metric). A brief discussion of the arrow of time within this framework is included. In this multi-time setting we show how to deal with transition regions, where one is passing on the baton from one time to to another, and investigate the fate of the subdominant clock. For this purpose minisuperspace is best seen as a dispersive medium, with packets moving with a group speed distinct from the phase speed. We show that the motion of the packets' peaks reproduces the classical limit even during the transition periods, and for subdominant clocks once the transition is over. Strong deviations from the coherent/semi-classical limit are expected in these cases, however. Could these be a "sign of the times", accessible notably in the transition period (from matter to Lambda domination) we live in?

The QSNET consortium is building a UK network of next-generation atomic and molecular clocks that will achieve unprecedented sensitivity in testing variations of the fine structure constant, $\alpha$, and the electron-to-proton mass ratio, $\mu$. This in turn will provide more stringent constraints on a wide range of fundamental and phenomenological theories beyond the Standard Model and on dark matter models.

This work investigates the formation of primordial black holes within a radiation fluid with an anisotropic pressure. We focus our attention on the initial conditions describing cosmological perturbations on the super horizon regime, using a covariant form of the equation of state in terms of pressure and energy density gradients. The effect of the anisotropy is to modify the initial shape of the cosmological perturbations with respect to the isotropic case. Using the dependence of the threshold $\delta_\mathrm{c}$ for primordial black holes on the shape of cosmological perturbations, derived in the isotropic limit, we estimate here how the threshold $\delta_\mathrm{c}$ is varying with respect to the amplitude of the anisotropies. If the anisotropy is large enough this could lead to a significant variation of the abundance of PBHs.

With a laser interferometric gravitational-wave detector in separate free flying spacecraft, the only way to achieve detection is to mitigate the dominant noise arising from the frequency fluctuations of the lasers via postprocessing. The noise can be effectively filtered out on the ground through a specific technique called time-delay interferometry (TDI), which relies on the measurements of time-delays between spacecraft and careful modeling of how laser noise enters the interferometric data. Recently, this technique has been recast into a matrix-based formalism by several authors, offering a different perspective on TDI, particularly by relating it to principal component analysis (PCA). In this work, we demonstrate that we can cancel laser frequency noise by directly applying PCA to a set of shifted data samples, without any prior knowledge of the relationship between single-link measurements and noise, nor time-delays. We show that this fully data-driven algorithm achieves a gravitational-wave sensitivity similar to classic TDI.

In earlier works, we studied the validity of Extended Effective Field Theory of Inflation (EEFToI) in the regime where initial conditions are set with dispersion relations $\omega^2 \propto k^6$. We had also evaluated and examined the power spectrum for some interesting corners of the parameter space. In this paper, we compute the bispectrum in the EEFToI, take a closer look at the strong coupling constraints and calculate the size of the non-Gaussianities in those regions of parameter space. We also investigate the shape of triangles that contribute to the enhancement of non-Gaussianities in this regime. We find that there are allowed parts of parameter spaces where EEFToI description with initial conditions set with $\omega^2 \propto k^6$ is sensible and interesting.

Non-minimal coupling between the Riemann curvature and the electromagnetic field appears as quantum corrections when gravity is coupled to the standard model of particle physics. The non-minimal coupling is expected to be dominant in the strong-gravity regions such as near black-holes or the early Universe. With better instruments planned shortly, electromagnetic fields are an important source of astrophysical observations to test general relativity in strong gravity regimes. However, to precisely test general relativity in strong gravity using electromagnetic fields, it is \emph{imperative} to obtain constraints on the non-minimal coupling parameter to electromagnetic fields. As a step in this direction, we calculate the deflection angle of non-minimally coupled electromagnetic fields in the vicinity of a dynamical, spherically symmetric black hole described by the Sultana-Dyer metric. We compare the deflection angle of the photon modes for the Sultana-Dyer black hole with the Schwarzschild black hole. We show that the difference in the deflection angle for the Schwarzschild black hole is always negative, while for Sultana-Dyer is always positive. Thus, our analysis points out that the two black holes provide a distinct signature irrespective of the black hole mass. We discuss the implications of the results for future astrophysical observations.

We report the results from a haloscope search for axion dark matter in the $3.3\text{-}4.2~{\mu}$eV mass range. This search excludes the axion-photon coupling predicted by one of the benchmark models of "invisible" axion dark matter, the KSVZ model. This sensitivity is achieved using a large-volume cavity, a superconducting magnet, an ultra low noise Josephson parametric amplifier, and sub-Kelvin temperatures. The validity of our detection procedure is ensured by injecting and detecting blind synthetic axion signals.

We present the validation of ADBSat, a novel implementation of the panel method including a fast pseudo-shading algorithm, that can quickly and accurately determine the forces and torques on satellites in free-molecular flow. Our main method of validation is comparing test cases between ADBSat, the current de facto standard of direct simulation Monte Carlo (DSMC), and published literature. ADBSat broadly performs well, except where deep concavities are present in the satellite models. The shading algorithm also experiences problems when a large proportion of the satellite surface area is oriented parallel to the flow, but this can be mitigated by examining the body at small angles to this configuration (${\pm}$ 0.1{\deg}). We determine the error interval on ADBSat outputs to be 1-3% whilst exhibiting a significantly shorter runtime than comparable methods. ADBSat can therefore be used as a viable alternative to DSMC for preliminary design studies involving a wide range of geometries and cases. It can also be used in a complementary manner to identify cases that warrant further investigation using methods such as DSMC. Thus, it is an ideal tool for determining the aerodynamic characteristics of future missions to VLEO.

We present new constraints on the merging rates of planetary-mass and asteroid-mass primordial black hole binaries using limits on continuous waves(quasi-monochromatic, quasi-infinite duration signals) derived from an all-sky search for isolated compact objects in the first six months of the third observing run (O3a) of LIGO/Virgo. We calculate the merging rates of these binaries in a model-independent way, and convert them to constraints on the primordial black hole abundance with minimal modelling assumptions. Our results show that we are sensitive to sources at most $\mathcal{O}(10$ pc) away for systems with chirp masses of $\mathcal{O}(10^{-5}M_\odot)$ at gravitational-wave frequencies around 30-40 Hz. These results also show that continuous-wave searches could in the future directly probe the existence of planetary-mass and asteroid-mass primordial black holes, especially those in binaries with asymmetric mass ratios. Furthermore, they demonstrate that new methods accounting for the full nonlinear gravitational-wave frequency evolution are needed to improve constraints on primordial black holes.