Papers for Monday, Feb 22 2021

White dwarfs with metal-polluted atmospheres have been studied widely in the context of the accretion of rocky debris from evolved planetary systems. One open question is the geometry of accretion and how material arrives and mixes in the white dwarf surface layers. Using the 3D radiation-hydrodynamics code CO$^5$BOLD, we present the first transport coefficients in degenerate star atmospheres which describe the advection-diffusion of a passive scalar across the surface-plane. We couple newly derived horizontal diffusion coefficients with previously published vertical diffusion coefficients to provide theoretical constraints on surface spreading of metals in white dwarfs. Our grid of 3D simulations probes the vast majority of the parameter space of convective white dwarfs, with pure-hydrogen atmospheres in the effective temperature range 6000-18000 K and pure-helium atmospheres in the range 12000-34000 K. Our results suggest that warm hydrogen-rich atmospheres (DA; $\gtrsim$13000 K) and helium-rich atmospheres (DB, DBA; $\gtrsim$30000 K) are unable to efficiently spread the accreted metals across their surface, regardless of the time dependence of accretion. This result may be at odds with the current non-detection of surface abundance variations at white dwarfs with debris discs. For cooler hydrogen- and helium-rich atmospheres, we predict a largely homogeneous distribution of metals across the surface within a vertical diffusion timescale. This is typically less than 0.1 per cent of disc lifetime estimates, a quantity which is revisited in this paper using the overshoot results. These results have relevance for studies of the bulk composition of evolved planetary systems and models of accretion disc physics.

Massive black holes often exist within dwarf galaxies, and both simulations and observations have shown that a substantial fraction of these may be off-center with respect to their hosts. We trace the evolution of off-center massive black holes (MBHs) in dwarf galaxies using cosmological hydrodynamical simulations, and show that the reason for off-center locations is mainly due to galaxy-galaxy mergers. We calculate dynamical timescales and show that off-center MBHs are unlikely to sink to their galaxys' centers within a Hubble time, due to the shape of the hosts' potential wells and low stellar densities. These wandering MBHs are unlikely to be detected electromagnetically, nor is there a measurable dynamical effect on the galaxy's stellar population. We conclude that off-center MBHs may be common in dwarfs, especially if the mass of the MBH is small or the stellar mass of the host galaxy is large. However detecting them is extremely challenging, because their accretion luminosities are very low and they do not measurably alter the dynamics of their host galaxies.

We use the data from Gaia Early Data Release 3 (EDR3) to study the kinematic properties of Milky Way globular clusters. We measure the mean parallaxes and proper motions (PM) for 170 clusters, determine the PM dispersion profiles for more than 100 clusters, uncover rotation signatures in more than 20 objects, and find evidence for radial or tangential PM anisotropy in a dozen richest clusters. At the same time, we use the selection of cluster members to explore the reliability and limitations of the Gaia catalogue itself. We find that the formal uncertainties on parallax and PM are underestimated by 10-20% in dense central regions even for stars that pass numerous quality filters. We explore the the spatial covariance function of systematic errors, and determine a lower limit on the uncertainty of average parallaxes and PM at the level 0.01 mas and 0.025 mas/yr, respectively. Finally, a comparison of mean parallaxes of clusters with distances from various literature sources suggests that the parallaxes (after applying the zero-point correction suggested by Lindegren et al. 2020) are overestimated by 0.006-0.009 mas. Despite these caveats, the quality of Gaia astrometry has been significantly improved in EDR3 and provides valuable insights into the properties of star clusters.

A Circumnuclear Disk (CND) of molecular gas occupies the central few parsecs of the Galactic Center. It is likely subject to turbulent disruptions from violent events in its surrounding environment, but the effect of such perturbations has not yet been investigated in detail. Here we perform 3D, N-body/smoothed particle hydrodynamic (SPH) simulations with an adapted general turbulence driving method to investigate the CND's structural evolution, in particular its reaction to varied scales of injected turbulence. We find that, because of shear flow in the disk, transient arcs of gas (streams) naturally arise when turbulence is driven on large scales (up to $\sim4$~pc), as might occur when a supernova blast wave encounters the CND. Because energetic events arise naturally and often in the central parsecs of our Galaxy, this result suggests that the transient structures that characterize the CND do not imply that the CND itself is a transient structure. We also note that features similar to the density concentrations, or ``clumps'', detailed in literature emerge when we account for the observed orientation of the disk and for the spatial resolution of observations. As such, clumps could be an artifact of observational limitations.

The evolution of a single raindrop falling below a cloud is governed by fluid dynamics and thermodynamics fundamentally transferable to planetary atmospheres beyond modern Earth's. Here, we show how three properties that characterize falling raindrops -- raindrop shape, terminal velocity, and evaporation rate -- can be calculated as a function of raindrop size in any planetary atmosphere. We demonstrate that these simple, interrelated characteristics tightly bound the possible size range of raindrops in a given atmosphere, independently of poorly understood growth mechanisms. Starting from the equations governing raindrop falling and evaporation, we demonstrate that raindrop ability to vertically transport latent heat and condensible mass can be well captured by a new dimensionless number. Our results have implications for precipitation efficiency, convective storm dynamics, and rainfall rates, which are properties of interest for understanding planetary radiative balance and (in the case of terrestrial planets) rainfall-driven surface erosion.

We present comprehensive, multi-wavelength observations of AT 2016jbu, an interacting transient. High cadence photometric coverage reveals that AT 2016jbu underwent significant photometric variability followed by two luminous events, the latter of which reached an absolute magnitude of M$_V\sim$-18.5 mag. This is similar to the transient SN 2009ip whose nature is still debated. Spectra are dominated by narrow emission lines and show a blue continuum during the peak of the second event. AT 2016jbu shows signatures of a complex, non-homogeneous circumstellar material (CSM). We see slowly evolving asymmetric hydrogen line profiles, with velocities of 500 km $s^{-1}$ seen in narrow emission features from a slow moving CSM, and up to 10,000 km $s^{-1}$ seen in broad absorption from some high velocity material. Similar velocities are seen in other SN 2009ip-like transients. Late-time spectra ($\sim$+1 year) show a lack of forbidden emission lines expected from a core-collapse supernova during the nebular phase and are dominated by strong emission from H, He I and Ca II. Strong asymmetric emission features, a bumpy lightcurve, and continually evolving spectra suggest late time CSM interaction is inhibiting the emergence of any nebular features. We compare the evolution of H$\alpha$ among SN 2009ip-like transients and find possible evidence for orientation angle effects. The light-curve evolution of AT 2016jbu suggests similar, but not identical, circumstellar environments to other SN 2009ip-like transients. In Paper II we continue the discussion of AT 2016jbu and SN 2009ip-like transients and using the data presented here, we focus on the local environment, the progenitor, and on modelling the transient itself.

We revisit the system consisting of a neutron star that harbors a small, possibly primordial, black hole at its center, focusing on a nonspinning black hole embedded in a nonrotating neutron star. Extending earlier treatments, we provide an analytical treatment describing the rate of secular accretion of the neutron star matter onto the black hole, adopting the relativistic Bondi accretion formalism for stiff equations of state that we presented elsewhere. We use these accretion rates to sketch the evolution of the system analytically until the neutron star is completely consumed. We also perform numerical simulations in full general relativity for black holes with masses up to nine orders of magnitude smaller than the neutron star mass, including a simulation of the entire evolution through collapse for the largest black hole mass. We construct relativistic initial data for these simulations by generalizing the black hole puncture method to allow for the presence of matter, and evolve these data with a code that is optimally designed to resolve the vastly different length scales present in this problem. We compare our analytic and numerical results, and provide expressions for the lifetime of neutron stars harboring such endoparasitic black holes.

In the second of two papers on the peculiar interacting transient AT 2016jbu, we present the bolometric lightcurve, identification and analysis of the progenitor candidate, as well as preliminary modelling to help elucidate the nature of this event. We identify the progenitor candidate for AT 2016jbu in quiescence, and find it to be consistent with a $\sim$20 M$_{\odot}$ yellow hypergiant surrounded by a dusty circumstellar shell. We see evidence for significant photometric variability in the progenitor, as well as strong H$\alpha$ emission consistent with pre-existing circumstellar material. The age of the resolved stellar population surrounding AT 2016jbu, as well as integral-field unit spectra of the region support a progenitor age of >16 Myr, again consistent with a progenitor mass of $\sim$20 M$_{\odot}$. Through a joint analysis of the velocity evolution of AT 2016jbu, and the photospheric radius inferred from the bolometric lightcurve, we find that the transient is consistent with two successive outbursts or explosions. The first outburst ejected a shell of material with velocity 650 km $s^{-1}$, while the second more energetic event ejected material at 4500 km $s^{-1}$. Whether the latter is the core-collapse of the progenitor remains uncertain, as the required ejecta mass is relatively low (few tenths of M$_{\odot}$). We also place a restrictive upper limit on the ejected $^{56}$Ni mass of <0.016 M$_{\odot}$. Using the BPASS code, we explore a wide range of possible progenitor systems, and find that the majority of these are in binaries, some of which are undergoing mass transfer or common envelope evolution immediately prior to explosion. Finally, we use the SNEC code to demonstrate that the low-energy explosion of some of these systems together with sufficient CSM can reproduce the overall morphology of the lightcurve of AT 2016jbu.

Connecting the composition of planet-forming disks with that of gas giant exoplanet atmospheres, in particular through C/O ratios, is one of the key goals of disk chemistry. Small hydrocarbons like $\rm C_2H$ and $\rm C_3H_2$ have been identified as tracers of C/O, as they form abundantly under high C/O conditions. We present resolved $\rm C_3H_2$ observations from the TW Hya Rosetta Stone Project, a program designed to map the chemistry of common molecules at $15-20$ au resolution in the TW Hya disk. Augmented by archival data, these observations comprise the most extensive multi-line set for disks of both ortho and para spin isomers spanning a wide range of energies, $E_u=29-97$ K. We find the ortho-to-para ratio of $\rm C_3H_2$ is consistent with 3 throughout extent of the emission, and the total abundance of both $\rm C_3H_2$ isomers is $(7.5-10)\times10^{-11}$ per H atom, or $1-10$% of the previously published $\rm C_2H$ abundance in the same source. We find $\rm C_3H_2$ comes from a layer near the surface that extends no deeper than $z/r=0.25$. Our observations are consistent with substantial radial variation in gas-phase C/O in TW Hya, with a sharp increase outside $\sim30$ au. Even if we are not directly tracing the midplane, if planets accrete from the surface via, e.g., meridonial flows, then such a change should be imprinted on forming planets. Perhaps interestingly, the HR 8799 planetary system also shows an increasing gradient in its giant planets' atmospheric C/O ratios. While these stars are quite different, hydrocarbon rings in disks are common, and therefore our results are consistent with the young planets of HR 8799 still bearing the imprint of their parent disk's volatile chemistry.

We present the first set of cosmological baryonic zoom-in simulations of galaxies including dissipative self-interacting dark matter (dSIDM). These simulations utilize the Feedback In Realistic Environments (FIRE-2) galaxy formation physics, but allow the dark matter to have dissipative self-interactions analogous to Standard Model forces, parameterized by the self-interaction cross-section per unit mass, $(\sigma/m)$, and the dimensionless degree of dissipation, $0<f_{\rm diss}<1$. We survey this parameter space, including constant and velocity-dependent cross-sections, and focus on structural and kinematic properties of dwarf galaxies with $M_{\rm halo} \simeq 10^{10-11} {\rm M}_{\odot}$. Central density profiles of simulated dwarfs become cuspy when $(\sigma/m)_{\rm eff} \gtrsim 0.1\,{\rm cm^{2}\,g^{-1}}$ (and $f_{\rm diss}=0.5$ as fiducial). The power-law slopes asymptote to $\alpha \approx -1.5$ in low-mass dwarfs independent of cross-section, which arises from a dark matter ``cooling flow''. Through comparisons with dark matter only simulations, we find the profile in this regime is insensitive to the inclusion of baryons. However, when $(\sigma/m)_{\rm eff} \ll 0.1\,{\rm cm^{2}\,g^{-1}}$, baryonic effects can produce cored density profiles comparable to non-dissipative cold dark matter (CDM) runs but at smaller radii. Simulated galaxies with $(\sigma/m) \gtrsim 10\,{\rm cm^{2}\,g^{-1}}$ develop significant coherent rotation of dark matter, accompanied by halo deformation, but this is unlike the well-defined thin ``dark disks'' often attributed to baryon-like dSIDM. The density profiles in this high cross-section model exhibit lower normalizations given the onset of halo deformation. For our surveyed dSIDM parameters, halo masses and galaxy stellar masses do not show appreciable difference from CDM, but dark matter kinematics and halo concentrations/shapes can differ.

Using radio observations with the Green Bank Telescope, evidence has now been found for a second five-membered ring in the dense cloud Taurus Molecular Cloud-1 (TMC-1). Based on additional observations of an ongoing, large-scale, high-sensitivity spectral line survey (GOTHAM) at centimeter wavelengths toward this source, we have used a combination of spectral stacking, Markov chain Monte Carlo (MCMC), and matched filtering techniques to detect 2-cyanocyclopentadiene, a low-lying isomer of 1-cyanocyclopentadiene, which was recently discovered there by the same methods. The new observational data also yields a considerably improved detection significance for the more stable isomer and evidence for several individual transitions between 23 - 32 GHz. Through our MCMC analysis, we derive total column densities of $8.3\times10^{11}$ and $1.9\times10^{11}$ cm$^{-2}$ for 1- and 2-cyanocyclopentadiene respectively, corresponding to a ratio of 4.4(6) favoring the former. The derived abundance ratios point towards a common formation pathway - most likely being cyanation of cyclopentadiene by analogy to benzonitrile.

The Large-aperture Experiment to Detect the Dark Age (LEDA) was designed to measure brightness temperature fluctuations in the Cosmic Microwave Background due to 21-cm absorption/emission by neutral hydrogen in the early Universe. Measurements of 21-cm fluctuations will provide constraints on astrophysical processes during Cosmic Dawn ($z \approx 15-30$). LEDA uses observations made with the Owens Valley Long Wavelength Array (OVRO-LWA) in California, a compact imaging interferometer. We present a power spectrum of observations at 48 MHz ($z \approx 28$), and an analysis of the sensitivity of the OVRO-LWA telescope. Recently, Eastwood et al. published a power spectrum of OVRO-LWA observations generated from m-mode analysis, using statistical techniques to remove foregrounds. We use a method that isolates foreground power in a "wedge" in the power spectrum; outside of the wedge, in the "Cosmic Dawn" window, 21-cm fluctuations may be detected. Using 4 hrs of incoherently integrated OVRO-LWA observations, we measure a power of $\Delta^2(k) \approx 2 \times 10^{12}$ mK$^2$ in the Cosmic Dawn window. Given that a value of $\Delta^2(k) \approx 100$ mK$^2$ is predicted from theory, we conclude that our power spectrum is dominated by telescope thermal noise and systematic effects. By modelling the OVRO-LWA thermal noise, we show that OVRO-LWA has sufficient sensitivity for a 21-cm detection if 3000 hrs of observations are integrated using a mix of incoherent and coherent averaging. We show that OVRO-LWA will then be able to detect theoretically predicted Ly-$\alpha$ and X-ray peaks in 21-cm power spectra at Cosmic Dawn, according to theoretical models.

We present ALMA $^{12}$CO (2-1), $^{13}$CO (2-1), C$^{18}$O (2-1) molecular line observations of a very young proto-brown dwarf system, ISO-OPH 200. We have conducted physical+chemical modelling of the complex internal structure for this system using the core collapse simulations for brown dwarf formation. The model at an age of $\sim$6000 yr can provide a good fit to the observed kinematics, spectra, and reproduce the complex structures seen in the moment maps. Results from modelling indicate that $^{12}$CO emission is tracing an extended ($\sim$1000 au) molecular outflow and a bright shock knot, $^{13}$CO is tracing the outer ($\sim$1000 au) envelope/pseudo-disc, and C$^{18}$O is tracing the inner ($\sim$500 au) pseudo-disc. The source size of $\sim$8.6 au measured in the 873$\mu$m image is comparable to the inner Keplerian disc size predicted by the model. A 3D model structure of ISO-OPH 200 suggests that this system is viewed partially through a wide outflow cavity resulting in a direct view of the outflow and a partial view of the envelope/pseudo-disc. We have argued that ISO-OPH 200 has been mis-classified as a Class Flat object due to the unusual orientation. The various signatures of this system, notably, the young $\sim$616 yr outflow dynamical age and high outflow rate ($\sim$1$\times$10$^{-7}$ M$_{\odot}$ yr$^{-1}$), silicate absorption in the 10$\mu$m mid-infrared spectrum, pristine ISM-like dust in the envelope/disc, comparable sizes of the extended envelope and outflow, indicate that ISO-OPH 200 is an early Class 0 stage system formed in a star-like mechanism via gravitational collapse of a very low-mass core.

Six binary-merger progenitors of Supernova 1987A (SN 1987A), whose properties are close to those of the blue supergiant Sanduleak -69 202, are exploded by neutrino heating in three-dimensions (3D), and light-curve calculations are performed in spherical symmetry, starting long after shock breakout. Our results confirm basic findings of previous works using single-star progenitors: (1) 3D neutrino-driven explosions with SN 1987A-like energies synthesize an amount of Ni-56 that is consistent with the radioactive tail of the light curve; (2) hydrodynamic models mix hydrogen inward to minimum velocities below 40 km/s as required by spectral observations of SN 1987A; and (3) for given explosion energy the efficiency of outward radioactive Ni-56 mixing depends mainly on high growth factors of Rayleigh-Taylor instabilities at the (C+O)/He and He/H composition interfaces and a weak interaction of fast plumes with the reverse shock occurring below the He/H interface. All binary-merger models possess presupernova radii consistent with the photometric radius of Sanduleak -69 202 and a structure of the outer layers allowing them to reproduce the observed initial luminosity peak in the first about 7 days. The light curve shape of models that mix about 0.5 Msun of hydrogen into the helium shell and that exhibit strong outward mixing of Ni-56 with maximum velocities exceeding the 3000 km/s observed for the bulk of ejected Ni-56 is in very good agreement with the dome of the SN 1987A light curve. A comparative analysis of light-curve models of SN 1987A based on 3D neutrino-driven explosions of single-star and binary-merger progenitors shows that only one binary model matches all observational constraints with one exception.

We propose a novel method to reconstruct high-resolution three-dimensional mass maps using data from photometric weak-lensing surveys. We apply an adaptive LASSO algorithm to perform a sparsity-based reconstruction on the assumption that the underlying cosmic density field is represented by a sum of Navarro-Frenk-White halos. We generate realistic mock galaxy shape catalogues by considering the shear distortions from isolated halos for the configurations matched to Subaru Hyper Suprime-Cam Survey with its photometric redshift estimates. We show that the adaptive method significantly reduces line-of-sight smearing that is caused by the correlation between the lensing kernels at different redshifts. Lensing clusters with lower mass limits of $10^{14.0} h^{-1}M_{\odot}$, $10^{14.7} h^{-1}M_{\odot}$, $10^{15.0} h^{-1}M_{\odot}$ can be detected with 1.5-$\sigma$ confidence at the low ($z<0.3$), median ($0.3\leq z< 0.6$) and high ($0.6\leq z< 0.85$) redshifts, respectively, with an average false detection rate of 0.022 deg$^{-2}$. The estimated redshifts of the detected clusters are systematically lower than the true values by $\Delta z \sim 0.03$ for halos at $z\leq 0.4$, but the relative redshift bias is below $0.5\%$ for clusters at $0.4<z\leq 0.85$. The standard deviation of the redshift estimation is $0.092$. Our method enables direct three-dimensional cluster detection with accurate redshift estimates.

In addition to the well-known gas phase mass-metallicity relation (MZR), recent spatially-resolved observations have shown that local galaxies also obey a mass-metallicity gradient relation (MZGR) whereby metallicity gradients can vary systematically with galaxy mass. In this work, we use our recently-developed analytic model for metallicity distributions in galactic discs, which includes a wide range of physical processes -- radial advection, metal diffusion, cosmological accretion, and metal-enriched outflows -- to simultaneously analyse the MZR and MZGR. We show that the same physical principles govern the shape of both: centrally-peaked metal production favours steeper gradients, and this steepening is diluted by the addition of metal-poor gas, which is supplied by inward advection for low-mass galaxies and by cosmological accretion for massive galaxies. The MZR and the MZGR both bend at galaxy stellar mass $\sim 10^{10} - 10^{10.5}\,\rm{M_{\odot}}$, and we show that this feature corresponds to the transition of galaxies from the advection-dominated to the accretion-dominated regime. We also find that both the MZR and MZGR strongly suggest that low-mass galaxies preferentially lose metals entrained in their galactic winds. While this metal-enrichment of the galactic outflows is crucial for reproducing both the MZR and the MZGR at the low-mass end, we show that the flattening of gradients in massive galaxies is expected regardless of the nature of their winds.

In this work, we explore the diversity of ionised gas kinematics (rotational velocity $v_{\phi}$ and velocity dispersion $\sigma_g$) and gas-phase metallicity gradients at $0.1 \leq z \leq 2.5$ using a compiled data set of 74 galaxies resolved with ground-based integral field spectroscopy. We find that galaxies with the highest and the lowest $\sigma_g$ have preferentially flat metallicity gradients, whereas those with intermediate values of $\sigma_g$ show a large scatter in the metallicity gradients. Additionally, steep negative gradients appear almost only in rotation-dominated galaxies ($v_{\phi}/\sigma_g > 1$), whereas most dispersion-dominated galaxies show flat gradients. We use our recently developed analytic model of metallicity gradients to provide a physical explanation for these observed trends. In the case of high $\sigma_g$, the inward radial advection of gas dominates over metal production and causes efficient metal mixing, thus giving rise to flat gradients. For low $\sigma_g$, it is the cosmic accretion of metal-poor gas diluting the metallicity that gives rise to flat gradients. Finally, the reason for intermediate $\sigma_g$ showing the steepest negative gradients is that both inward radial advection and cosmic accretion are weak as compared to metal production, which leads to the creation of steeper gradients. The larger scatter at intermediate $\sigma_g$ may be due in part to preferential ejection of metals in galactic winds, which can decrease the strength of the production term. Our analysis shows that gas kinematics plays a critical role in setting metallicity gradients in high-redshift galaxies.

HO Puppis (HO Pup) was considered as a Be-star candidate based on its gamma-Cassiopeiae-type light curve, but lacked spectroscopic confirmation. Using distance measured from Gaia Data Release 2 and the spectral-energy-distribution (SED) fit on broadband photometry, the Be-star nature of HO Pup is ruled out. Furthermore, based on the 28,700 photometric data points collected from various time-domain surveys and dedicated intensive-monitoring observations, the light curves of HO Pup closely resemble IW And-type stars (as pointed out in Kimura et al. 2020a), exhibiting characteristics such as quasi-standstill phase, brightening, and dips. The light curve of HO Pup displays various variability timescales, including brightening cycles ranging from 23 to 61 days, variations with periods between 3.9 days and 50 minutes during the quasi-standstill phase, and a semi-regular ~14-day period for the dip events. We have also collected time-series spectra (with various spectral resolutions), in which Balmer emission lines and other expected spectral lines for an IW And-type star were detected (even though some of these lines were also expected to be present for Be stars). We detect Bowen fluorescence near the brightening phase, and that can be used to discriminate between IW And-type stars and Be stars. Finally, despite only observing for four nights, the polarization variation was detected, indicating that HO Pup has significant intrinsic polarization.

Quasars behind the Galactic plane (GPQs) are important astrometric references and useful probes of Milky Way gas. However, the search for GPQs is difficult due to large extinctions and high source densities in the Galactic plane. Existing selection methods for quasars developed using high Galactic latitude (high-$b$) data cannot be applied to the Galactic plane directly because the photometric data obtained from high-$b$ regions and the Galactic plane follow different probability distributions. To alleviate this dataset shift problem for quasar candidate selection, we adopt a Transfer Learning Framework at both data and algorithm levels. At the data level, to make a training set in which dataset shift is modeled, we synthesize quasars and galaxies behind the Galactic plane based on SDSS sources and Galactic dust map. At the algorithm level, to reduce the effect of class imbalance, we transform the three-class classification problem for stars, galaxies, and quasars to two binary classification tasks. We apply XGBoost algorithm on Pan-STARRS1 (PS1) and AllWISE photometry for classification, and additional cut on Gaia proper motion to remove stellar contaminants. We obtain a reliable GPQ candidate catalog with 160,946 sources located at $|b|\leq 20^{\circ}$ in PS1-AllWISE footprint. Photometric redshifts of GPQ candidates achieved with XGBoost regression algorithm show that our selection method can identify quasars in a wide redshift range ($0<z\lesssim5$). This study extends the systematic searches for quasars to the dense stellar fields and shows the feasibility of using astronomical knowledge to improve data mining under complex conditions in the Big Data era.

The radio, optical, and $\gamma$-ray light curves of the blazar S5 1803+784, from the beginning of the {\it Fermi} Large Area Telescope (LAT) mission in August 2008 until December 2018, are presented. The aim of this work is to look for correlations among different wavelengths useful for further theoretical studies. We analyzed all the data collected by {\it Fermi} LAT for this source, taking into account the presence of nearby sources, and we collected optical data from our own observations and public archive data to build the most complete optical and $\gamma$-ray light curve possible. Several $\gamma$-ray flares ($\mathrm{F>2.3~10^{-7} ph(E>0.1 GeV)~cm^{-2}~s^{-1}}$) with optical coverage were detected, all but one with corresponding optical enhancement; we also found two optical flares without a $\gamma$-ray counterpart. We obtained two {\it Swift} Target of Opportunity observations during the strong flare of 2015. Radio observations performed with VLBA and EVN through our proposals in the years 2016-2020 were analyzed to search for morphological changes after the major flares. The optical/$\gamma$-ray flux ratio at the flare peak varied for each flare. Very minor optical V-I color changes were detected during the flares. The X-ray spectrum was well fitted by a power law with photon spectral index $\alpha$=1.5, nearly independent of the flux level: no clear correlation with the optical or the $\gamma$-ray emission was found. The $\gamma$-ray spectral shape was well fitted by a power law with average photon index $\alpha$= 2.2. These findings support an Inverse Compton origin for the high-energy emission of the source, nearly co-spatial with the optically emitting region. The radio maps showed two new components originating from the core and moving outwards, with ejection epochs compatible with the dates of the two largest $\gamma$-ray flares.

We present near-infrared (IR) spectra of two planetary nebula (PN) candidates in close lines of sight toward the Galactic center (GC) using GNIRS at Gemini North. High-resolution images from radio continuum and narrow-band IR observations reveal ring-like morphologies of these objects, and their mid-IR spectra from the Spitzer Space Telescope exhibit rich emission lines from highly excited species such as [S IV], [Ne III], [Ne V], and [O IV]. We also derive elemental abundances using the Cloudy synthetic models, and find an excess amount of the s-process element Krypton in both targets, which supports their nature as PNe. We estimate foreground extinction toward each object using near-IR hydrogen recombination lines, and find significant visual extinctions ($A_V > 20$). The distances inferred from the size versus surface brightness relation of other PNe are $9.0\pm1.6$ kpc and $7.6\pm1.6$ kpc for SSTGC 580183 and SSTGC 588220, respectively. These observed properties along with abundance patterns and their close proximity to Sgr A$^*$ (projected distances $<20$ pc) make it highly probable that these objects are the first confirmed PNe objects in the nuclear stellar disk. The apparent scarcity of such objects resembles the extremely low rate of PN formation in old stellar systems, but is in line with the current rate of the sustained star formation activity in the Central Molecular Zone.

Large format arrays of transition edge sensor (TES) are crucial for the next generation of X-ray space observatories. Such arrays are required to achieve an energy resolution of $\mathrm{\Delta}E<$3 eV full-width-half-maximum (FWHM) in the soft X-ray energy range. We are currently developing X-ray microcalorimeter arrays as a backup option for the X-IFU instrument on board of ATHENA space telescope, led by ESA and foreseen to be launched in 2031. In this contribution, we report on the development and the characterization of a uniform 32$\times$32 pixel array with (length$\times $width) 140$\times$30 $\mu$m$^2$ TiAu TESs, which have \textcolor{black}{a 2.3 $\mu$m} thick Au absorber for X-ray photons. The pixels have a typical normal resistance $R_\mathrm{n}$ = 121 m$\Omega$ and a critical temperature $T_\mathrm{c}\sim$ 90 mK. We performed extensive measurements on 60 pixels out of the array in order to show the uniformity of the array. We obtained an energy resolutions between 2.4 and 2.6 eV (FWHM) at 5.9 keV, measured in a single-pixel mode at AC bias frequencies ranging from 1 to 5 MHz, with a frequency domain multiplexing (FDM) readout system, which is developed at SRON/VTT. We also present the detector energy resolution at X-ray with different photon energies generated by a modulated external X-ray source from 1.45 keV up to 8.9 keV. Multiplexing readout across several pixels has also been performed to evaluate the impact of the thermal crosstalk to the instrument's energy resolution budget requirement. This value results in a derived requirement, for the first neighbour, that is less than 1$\times$10$^{-3}$ when considering the ratio between the amplitude of the crosstalk signal to an X-ray pulse (for example at 5.9 keV)

In this study, we investigated the energy partition of four confined circular-ribbon flares (CRFs) near the solar disk center, which are observed simultaneously by SDO, GOES, and RHESSI. We calculated different energy components, including the radiative outputs in 1$-$8, 1$-$70, and 70$-$370 {\AA}, total radiative loss, peak thermal energy derived from GOES and RHESSI, nonthermal energy in flare-accelerated electrons, and magnetic free energy before flares. It is found that the energy components increase systematically with the flare class, indicating that more energies are involved in larger flares. The magnetic free energies are larger than the nonthermal energies and radiative outputs of flares, which is consistent with the magnetic nature of flares. The ratio $\frac{E_{nth}}{E_{mag}}$ of the four flares, being 0.70$-$0.76, is considerably higher than that of eruptive flares. Hence, this ratio may serve as an important factor that discriminates confined and eruptive flares. The nonthermal energies are sufficient to provide the heating requirements including the peak thermal energy and radiative loss. Our findings impose constraint on theoretical models of confined CRFs and have potential implication for the space weather forecast.

We present the first BVR photometry, period variation, and photometric light-curve analysis of two poorly studied eclipsing binaries V1321 Cyg and CR Tau. Observations were carried out from November 2017 to January 2020 at the observatory of Uzhhorod National University. Period variations were studied using all available early published as well as our minima times. We have used newly developed ELISa code for the light curve analysis and determination of photometric parameters of both systems. We found that V1321 Cyg is a close detached eclipsing system with a low photometric mass ratio of $q=0.28$ which suggests that the binary is a post mass transfer system. No significant period changes in this system are detected. CR Tau is, on the other hand, a semi-detached system where the secondary component almost fills its Roche lobe. We detected a long-term period increase at a rate of $1.49 \times 10^{-7} d/y$, which support mass transfer from lower mass secondary component to the more massive primary.

The chemistry of the diffuse interstellar medium is driven by the combined influences of cosmic rays, ultraviolet (UV) radiation, and turbulence. Previously detected at the outer edges of photodissociation regions (PDRs) and formed from the reaction of C+ and OH, CO+ is the main chemical precursor of HCO+ and CO in a thermal, cosmic-ray, and UV-driven chemistry. Our aim was to test whether the thermal cosmic-ray and UV-driven chemistry is producing CO in diffuse interstellar molecular gas through the intermediate formation of CO+ We searched for CO+ absorption with the Atacama Large Millimeter Array (ALMA) toward two quasars with known Galactic foreground absorption from diffuse interstellar gas, J1717-3342 and J1744-3116, targeting the two strongest hyperfine components of the J=2-1 transition near 236 GHz. We could not detect CO+ but obtained sensitive upper limits toward both targets. The derived upper limits on the CO+ column densities represent about 4% of the HCO+ column densities. The corresponding upper limit on the CO+ abundance relative to H2 is <1.2 x 10^{-10}. The non-detection of CO$^+$ confirms that HCO+ is mainly produced in the reaction between oxygen and carbon hydrides, CH2+ or CH3+ , induced by suprathermal processes, while CO+ and HOC+ result from reactions of C+ with OH and H2O. The densities required to form CO molecules at low extinction are consistent with this scheme.

Abundances and partitioning of ices and gases produced by gas-grain chemistry are governed by adsorption and desorption on grains. Understanding astrophysical observations rely on laboratory measurements of adsorption and desorption rates on dust grains analogs. On flat surfaces, gas adsorption probabilities (or sticking coefficients) have been found close to unity for most gases. Here we report a strong decrease of the sticking coefficients of H2O and CO2 on substrates more akin to cosmic dust, such as submicrometer-sized particles of carbon and olivine, bare or covered with ice. This effect results from the local curvature of the grains, and then extends to larger grains made of aggregated small particles, such as fluffy or porous dust in more evolved media (e.g. circumstellar disks). The main astrophysical implication is that accretion rates of gases are reduced accordingly, slowing the growth of cosmic ices. Furthermore, volatile species that are not adsorbed on a grain at their freeze-out temperature will pertain in the gas phase, which will impact gas-ice partitions. We also found that thermal desorption of H2O is not modified by grains size, and thus the snowlines temperature should be independent on the dust size distribution.

Large-scale sky surveys have played a transformative role in our understanding of astrophysical transients, only made possible by increasingly powerful machine learning-based filtering to accurately sift through the vast quantities of incoming data generated. In this paper, we present a new real-bogus classifier based on a Bayesian convolutional neural network that provides nuanced, uncertainty-aware classification of transient candidates in difference imaging, and demonstrate its application to the datastream from the GOTO wide-field optical survey. Not only are candidates assigned a well-calibrated probability of being real, but also an associated confidence that can be used to prioritise human vetting efforts and inform future model optimisation via active learning. To fully realise the potential of this architecture, we present a fully-automated training set generation method which requires no human labelling, incorporating a novel data-driven augmentation method to significantly improve the recovery of faint and nuclear transient sources. We achieve competitive classification accuracy (FPR and FNR both below 1%) compared against classifiers trained with fully human-labelled datasets, whilst being significantly quicker and less labour-intensive to build. This data-driven approach is uniquely scalable to the upcoming challenges and data needs of next-generation transient surveys. We make our data generation and model training codes available to the community.

Context. The xGASS and xCOLD GASS surveys have measured the atomic (HI) and molecular gas (H2) content of a large and representative sample of nearby galaxies (redshift range of 0.01 $\lt$ z $\lt$ 0.05). Aims. We present optical longslit spectra for a subset of the xGASS and xCOLD GASS galaxies to investigate the correlation between radial metallicity profiles and cold gas content. In addition to data from Moran et al. (2012), this paper presents new optical spectra for 27 galaxies in the stellar mass range of 9.0 $\leq$ log Mstar/Msun $\leq$ 10.0. Methods. The longslit spectra were taken along the major axis of the galaxies, allowing us to obtain radial profiles of the gas-phase oxygen abundance (12 + log(O/H)). The slope of a linear fit to these radial profiles is defined as the metallicity gradient. We investigated correlations between these gradients and global galaxy properties, such as star formation activity and gas content. In addition, we examined the correlation of local metallicity measurements and the global HI mass fraction. Results. We obtained two main results: (i) the local metallicity is correlated with the global HI mass fraction, which is in good agreement with previous results. A simple toy model suggests that this correlation points towards a 'local gas regulator model'; (ii) the primary driver of metallicity gradients appears to be stellar mass surface density (as a proxy for morphology). Conclusions. This work comprises one of the few systematic observational studies of the influence of the cold gas on the chemical evolution of star-forming galaxies, as considered via metallicity gradients and local measurements of the gas-phase oxygen abundance. Our results suggest that local density and local HI mass fraction are drivers of chemical evolution and the gas-phase metallicity.

The Laser Guide Star Facility (LGSF), as the most important part of the adaptive optics system of the large ground-based telescope, is aimed to generate multiple laser guide stars at the sodium layer. Laser Launch Telescope is employed to implement this requirement by projecting the Gauss beam to the sodium layer with a small beam size in LGSF system. As the diffraction and interference effects of laser's long-distance transmission, the conventional optical design based on the geometrical optics mechanism cannot achieve the expected laser propagation. In this paper, we propose a method to design optical system for laser launch telescope based on the physical optics theorem to generate an acceptable light spot at the sodium layer in the atmosphere. Besides, a tolerance analysis method based on physical optics propagation is also demonstrated to be necessitated to optimize the system's instrumentation performance. The numerical results show that the optical design considering physical optics propagation is highly rewarding and even necessitated in many occasions, especially for laser beam propagation systems.

Evolution of the cosmic star formation rate (SFR) and molecular mass density is expected to be matched by a similarly strong evolution of the fraction of atomic hydrogen (HI) in the cold neutral medium (CNM). We use results from a recent commissioning survey for intervening 21-cm absorbers with the Australian Square Kilometre Array Pathfinder (ASKAP) to construct a Bayesian statistical model of the $N_\mathrm{HI}$-weighted harmonic mean spin temperature ($T_\mathrm{s}$) at redshifts between $z = 0.37$ and 1.0. We find that $T_\mathrm{s} \leq 274$K with 95 per cent probability, suggesting that at these redshifts the typical HI gas in galaxies at equivalent DLA column densities may be colder than the Milky Way interstellar medium $(T_\mathrm{s,MW} \sim 300$K). This result is consistent with an evolving CNM fraction that mirrors the molecular gas towards the peak in SFR at $z \sim 2$. We expect that future surveys for HI 21-cm absorption with the current SKA pathfinder telescopes will be able to provide constraints on the CNM fraction that are an order of magnitude greater than presented here.

High throughput optical system is defined to possess the features of both large field of view (FOV) and high resolution. However, it is full of challenge to design such a telescope with the two conflicting specifications at the same time. In this paper, we propose a method to design a high throughput telescope based on the classical off-axis Three-Mirror Anastigmat (TMA) configuration by introducing a scanning mechanism. We derive the optimum initial design for the TMA system with no primary aberrations through characteristic ray tracing. During the design process, a real exit pupil is necessitated to accommodate the scanning mirror. By gradually increasing the system's FOV during the optimization procedure, we finally obtained a high throughput telescope design with an F-number of 6, a FOV of 60$^{\circ}$*1.5$^{\circ}$, and a long focal length of 876mm. In addition, the tolerance analysis is also conducted to demonstrate the instrumentation feasibility. We believe that this kind of large rectangle FOV telescope with high resolution has broad future applications in the optical remote sensing field.

It is known that hydrodynamic triangular jets are unstable to antisymmetric kink perturbations. The inclusion of magnetic field may lead to the stabilisation of the jets. Jets and complex magnetic fields are ubiquitous in the solar atmosphere, which suggests the possibility of the kink instability in certain cases. The aim of the paper is to study the kink instability of triangular jets sandwiched between magnetic tubes/slabs and its possible connection to observed properties of the jets in the solar atmosphere. A dispersion equation governing the kink perturbations is obtained through matching of analytical solutions at the jet boundaries. The equation is solved analytically and numerically for different parameters of jets and surrounding plasma. The analytical solution is accompanied by a numerical simulation of fully nonlinear MHD equations for a particular situation of solar type II spicules. MHD triangular jets are unstable to the dynamic kink instability depending on the Alfven Mach number (the ratio of flow to Alfven speeds) and the ratio of internal and external densities. When the jet has the same density as the surrounding plasma, then only super Alfvenic flows are unstable. However, denser jets are unstable also in sub Alfvenic regime. Jets with an angle to the ambient magnetic field have much lower thresholds of instability than field-aligned flows. Growth times of the kink instability are estimated as 6-15 min for type I spicules and 5-60 s for type II spicules matching with their observed life times. Numerical simulation of full nonlinear equations shows that the transverse kink pulse locally destroys the jet in less than a minute in the conditions of type II spicules. Dynamic kink instability may lead to full breakdown of MHD flows and consequently to observed disappearance of spicules in the solar atmosphere.

Temperatures of the outer planet thermospheres exceed those predicted by solar heating alone by several hundred degrees. Enough energy is deposited at auroral regions to heat the entire thermosphere, but models predict that equatorward distribution is inhibited by strong Coriolis forces and ion drag. A better understanding of auroral energy deposition and circulation are critical to solving this so-called energy crisis. Stellar occultations observed by the Ultraviolet Imaging Spectrograph instrument during the Cassini Grand Finale were designed to map the thermosphere from pole to pole. We analyze these observations, together with earlier observations from 2016 and 2017, to create a two-dimensional map of densities and temperatures in Saturns thermosphere as a function of latitude and depth. The observed temperatures at auroral latitudes are cooler and peak at higher altitudes and lower latitudes than predicted by models, leading to a shallower meridional pressure gradient. Under modified geostrophy, we infer slower westward zonal winds that extend to lower latitudes than predicted, supporting equatorward flow from approximately 70 to 30 degrees latitude in both hemispheres. We also show evidence of atmospheric waves in the data that can contribute to equatorward redistribution of energy through zonal drag.

We analyse the coronal elemental abundances during a small flare using Hinode/EIS observations. Compared to the pre-flare elemental abundances, we observed a strong increase in coronal abundance of Ca XIV 193.84 {\AA}, an emission line with low first ionisation potential (FIP < 10 eV), as quantified by the ratio Ca/Ar during the flare. This is in contrast to the unchanged abundance ratio observed using Si X 258.38 {\AA}/S X 264.23 {\AA}. We propose two different mechanisms to explain the different composition results. Firstly, the small flare-induced heating could have ionised S, but not the noble gas Ar, so that the flare-driven Alfv\'en waves brought up Si, S and Ca in tandem via the ponderomotive force which acts on ions. Secondly, the location of the flare in strong magnetic fields between two sunspots may suggest fractionation occurred in the low chromosphere, where the background gas is neutral H. In this region, high-FIP S could behave more like a low-FIP than a high-FIP element. The physical interpretations proposed generate new insights into the evolution of plasma abundances in the solar atmosphere during flaring, and suggests that current models must be updated to reflect dynamic rather than just static scenarios.

The Cosmic Microwave Background anisotropies are difficult to measure at large angular scales. In this paper, we present a new analysis of the \Planck\ High Frequency Instrument data that brings the cosmological part and its major foreground signal close to the detector noise. The solar dipole signal, induced by the motion of the solar system with respect to the CMB, is a very efficient tool to calibrate a detector or a set of detectors with high accuracy. In this work, the solar dipole signal is used to extract corrections of the frequency maps offsets reducing significantly uncertainties. The solar dipole parameters are refined together with the improvement of the high frequency foregrounds, and of the CMB large scales cosmological anisotropies. The stability of the solar dipole parameters is a powerful way to control the galactic foregrounds removal in the component separation process. It is used to build a model for Spectral Energy Distribution spatial variations of the interstellar dust emission. The knowledge of these variations will help future CMB analyses in intensity, and also in polarization to measure faint signal related to the optical reionization depth and the tensor-to-scalar ratio of the primordial anisotropies. The results of this work are: improved solar dipole parameters, a new interstellar dust model, and a large scale cosmological anisotropies map.

In their recent paper, Chan and Lee discuss an interesting possibility: radio continuum emission from a dwarf irregular galaxy may be used to constrain upper limits on the cross section of annihilating dark matter. They claim that the contributions from nonthermal and thermal emission can be estimated with such accuracy that one can place new upper limits on the annihilation cross section. We argue that the observations presented can be explained entirely with a standard spectrum and no contribution from dark matter. As a result, the estimated upper limits of Chan and Lee are atleast by a factor of 100 too low.

The analysis of 20 years of spectrophotometric data of the double shell planetary nebula PM\,1-188 is presented, aiming to determine the time evolution of the emission lines and the physical conditions of the nebula, as a consequence of the systematic fading of its [WC\,10] central star whose brightness has declined by about 10 mag in the past 40 years. Our main results include that the [\ion{O}{iii}], [\ion{O}{ii}], [\ion{N}{ii}] line intensities are increasing with time in the inner nebula as a consequence of an increase in electron temperature from 11,000 K in 2005 to more than 14,000 K in 2018, due to shocks. The intensity of the same lines are decreasing in the outer nebula, due to a decrease in temperature, from 13,000 K to 7,000 K, in the same period. The chemical composition of the inner and outer shells was derived and they are similar. Both nebulae present subsolar O, S and Ar abundances, while they are He, N and Ne rich. For the outer nebula the values are 12+log He/H= 11.13$\pm$0.05, 12+log O/H = 8.04$\pm$0.04, 12+log N/H= 7.87$\pm$0.06, 12+log S/H = 7.18$\pm$0.10 and 12+log Ar = 5.33$\pm$0.16. The O, S and Ar abundances are several times lower than the average values found in disc non-Type I PNe, and are reminiscent of some halo PNe. From high resolution spectra, an outflow in the N-S direction was found in the inner zone. Position-velocity diagrams show that the outflow expands at velocities in the $-$150 to 100 km s$^{-1}$ range, and both shells have expansion velocities of about 40 km s$^{-1}$.

The so-called aromatic infrared bands are attributed to emission of polycyclic aromatic hydrocarbons. The observed variations toward different regions in space are believed to be caused by contributions of different classes of PAH molecules, i.e. with respect to their size, structure, and charge state. Laboratory spectra of members of these classes are needed to compare them to observations and to benchmark quantum-chemically computed spectra of these species. In this paper we present the experimental infrared spectra of three different PAH dications, naphthalene$^{2+}$, anthracene$^{2+}$, and phenanthrene$^{2+}$, in the vibrational fingerprint region 500-1700~cm$^{-1}$. The dications were produced by electron impact ionization of the vapors with 70 eV electrons, and they remained stable against dissociation and Coulomb explosion. The vibrational spectra were obtained by IR predissociation of the PAH$^{2+}$ complexed with neon in a 22-pole cryogenic ion trap setup coupled to a free-electron infrared laser at the Free-Electron Lasers for Infrared eXperiments (FELIX) Laboratory. We performed anharmonic density-functional theory calculations for both singly and doubly charged states of the three molecules. The experimental band positions showed excellent agreement with the calculated band positions of the singlet electronic ground state for all three doubly charged species, indicating its higher stability over the triplet state. The presence of several strong combination bands and additional weaker features in the recorded spectra, especially in the 10-15~$\mu$m region of the mid-IR spectrum, required anharmonic calculations to understand their effects on the total integrated intensity for the different charge states. These measurements, in tandem with theoretical calculations, will help in the identification of this specific class of doubly-charged PAHs as carriers of AIBs.

The spherical Fourier-Bessel (SFB) decomposition is a natural choice for the radial/angular separation that allows optimal extraction of cosmological information from large volume galaxy surveys. In this paper we develop a SFB power spectrum estimator that allows the measurement of the largest angular and radial modes with the next generation of galaxy surveys. The code measures the pseudo-SFB power spectrum, and takes into account mask, selection function, pixel window, and shot noise. We show that the local average effect is significant only in the largest-scale mode, and we provide an analytical covariance matrix. By imposing boundary conditions at the minimum and maximum radius encompassing the survey volume, the estimator does not suffer from the numerical instabilities that have proven challenging in the past. The estimator is demonstrated on simplified Roman-like, SPHEREx-like, and Euclid-like mask and selection functions. For intuition and validation, we also explore the SFB power spectrum in the Limber approximation. We release the associated public code written in Julia.

Single phonon excitations are sensitive probes of light dark matter in the keV-GeV mass window. For anisotropic target materials, the signal depends on the direction of the incoming dark matter wind and exhibits a daily modulation. We discuss in detail the various sources of anisotropy, and carry out a comparative study of 26 crystal targets, focused on sub-MeV dark matter benchmarks. We compute the modulation reach for the most promising targets, corresponding to the cross section where the daily modulation can be observed for a given exposure, which allows us to combine the strength of DM-phonon couplings and the amplitude of daily modulation. We highlight Al$_2$O$_3$ (sapphire), CaWO$_4$ and h-BN (hexagonal boron nitride) as the best polar materials for recovering a daily modulation signal, which feature $\mathcal{O}(1 - 100)\%$ variations of detection rates throughout the day, depending on the dark matter mass and interaction. The directional nature of single phonon excitations offers a useful handle to mitigate backgrounds, which is crucial for fully realizing the discovery potential of near future experiments.

A problematic feature of low energy scale inflationary models, such as Starobinsky inflation, in a spatially closed universe is the occurrence of a recollapse and a big crunch singularity before inflation can even set in. In a recent work it was shown that this problem can be successfully resolved in loop quantum cosmology for a large class of initial conditions due to a non-singular cyclic evolution and a hysteresis-like phenomena. However, for certain highly unfavorable initial conditions the onset of inflation was still difficult to obtain. In this work, we explore the role of dissipative particle production, which is typical in warm inflation scenario, in the above setting. We find that entropy production sourced by such dissipative effects makes hysteresis-like phenomena stronger. As a result, the onset of inflation is quick in general including for highly unfavorable initial conditions where it fails or is significantly delayed in the absence of dissipative effects. We phenomenologically consider three warm inflation scenarios with distinct forms of dissipation coefficient, and from dynamical solutions and phase space portraits find that the phase space of favorable initial conditions turns out to be much larger than in cold inflation.

Gravitational nonlocality effects at large scale can be investigated using the cosmological structure formation. In particular, it is possible to show if and how gravitational nonlocality modifies the clustering properties of galaxies and of clusters of galaxies. The thermodynamics of such systems can be used to obtain important information about the effects of such lager scale nonlocality on clustering. We will compare the effects of such larger scale nonlocality with observational data. It will be demonstrated that the observations seem to point to a characteristic scale of nonlocality of the gravitational interactions at galactic scales. However, at larger scales such statistical inferences are much weaker.

I review the "Trans-Planckian Censorship Conjecture" (TCC) and its implications for cosmology, in particular for the inflationary universe scenario. Whereas the inflationary scenario is tightly constrained by the TCC, alternative early universe scenarios are not restricted.

We provide a first systematic and quantitative discussion of the phenomenology of the non-relativistic effective Hamiltonian describing the nuclear scattering process for a Weakly Interacting Massive Particle (WIMP) of arbitrary spin $j_\chi$. To this aim we obtain constraints from a representative sample of present direct detection experiments assuming the WIMP-nucleus scattering process to be driven by each one of the 44 effective couplings that arise for $j_\chi\le$ 2. We find that a high value of the multipolarity $s\le 2 j_\chi$ of the coupling, related to the power of the momentum transfer $q$ appearing in the scattering amplitude, leads to a suppression of the expected rates and pushes the expected differential spectra to large recoil energies $E_R$. For $s\le$ 4 the effective scales probed by direct detection experiments can be suppressed by up to 5 orders of magnitude compared to the case of a standard spin-independent interaction. For operators with large $s$ the expected differential spectra can be pushed to recoil energies in the MeV range, with the largest part of the signal concentrated at $E_R\gtrsim$ 100 keV and a peculiar structure of peaks and minima arising when both the nuclear target and the WIMP are heavy. As a consequence the present bounds on the effective operators can be significantly improved by extending the recoil energy intervals to higher recoil energies. Our analysis assumes effective interaction operators that are irreducible under the rotation group. Such operators drive the interactions of high-multipole dark matter candidates, i.e. states that possess only the highest multipole allowed by their spin. As a consequence our analysis represents also the first phenomenological study of the direct detection of quadrupolar, octupolar, and hexadecapolar dark matter.

We study dynamical screening effects of nuclear charge on big bang nucleosynthesis (BBN). A moving ion in plasma creates a distorted electric potential leading to a screening effect which is different from the standard static Salpeter formula. We consider the electric potential for a moving test charge, taking into account dielectric permittivity in the unmagnetized Maxwellian plasma during the BBN epoch. Based on the permittivity in a BBN plasma condition, we present the Coulomb potential for a moving nucleus, and show that enhancement factor for the screening of the potential increases the thermonuclear reaction rates by a factor order of 10^(-7). In the Gamow energy region for nuclear collisions, we find that the contribution of the dynamical screening is less than that of the static screening case, consequently which primordial abundances hardly change. Based on the effects of dynamical screening under various possible astrophysical conditions, we discuss related plasma properties required for possible changes of the thermal nuclear reactions.

Current cosmological data restricts Holographic Dark Energy (HDE) to a regime of parameter space with a turning point in the Hubble parameter. In particular, Cosmic Microwave Background (CMB) and baryon acoustic oscillation (BAO) data together favour a turning point in the observational regime, thereby alleviating Hubble tension. Nevertheless, including Pantheon supernovae (SNE) pushes $H_0$ back to a value consistent with Planck-$\Lambda$CDM. Noting that SNE are weighted to lower redshifts relative to BAO, this amounts to an evolution in $H_0$ with redshift within the HDE model. Since $H_0$ is an integration constant in any FLRW cosmology, this suggests that HDE may be at odds with the cosmological paradigm.

If a homogeneous field evolves within a medium, with the latter gradually picking up a temperature, then the friction felt by the field depends on how its evolution rate compares with medium time scales. We suggest a framework which permits to incorporate the contributions from all medium time scales. As an example, we illustrate how warm axion inflation can be described by inputting the retarded pseudoscalar correlator of a thermal Yang-Mills plasma. Adopting a semi-realistic model for the latter, and starting the evolution at almost vanishing temperature, we show how the system heats up and then enters the "weak" or "strong" regime of warm inflation. Previous approximate treatments are scrutinized.

The symmetry energy and its density dependence are crucial inputs for many nuclear physics and astrophysics applications, as they determine properties ranging from the neutron-skin thickness of nuclei to the crust thickness and the radius of neutron stars. Recently, PREX-II reported a value of $0.29\pm 0.07$ fm for the neutron-skin thickness of $^{208}$Pb, implying a slope parameter $L = 110\pm 37$ MeV, larger than most ranges obtained from microscopic calculations and other nuclear experiments. We use a nonparametric equation of state representation based on Gaussian processes to constrain the symmetry energy $S_0$, $L$, and $R_{\rm skin}^{^{208}\text{Pb}}$ directly from observations of neutron stars with minimal modeling assumptions. The resulting astrophysical constraints from heavy pulsar masses, LIGO/Virgo, and NICER clearly favor smaller values of the neutron skin and $L$, as well as negative symmetry incompressibilities. Combining astrophysical data with PREX-II and chiral effective field theory constraints yields $S_0 = 34^{+3}_{-3}$ MeV, $L=58^{+19}_{-19}$ MeV, and $R_{\rm skin}^{^{208}\text{Pb}} = 0.19^{+0.03}_{-0.04}$ fm.

The Hubble tension is shown to be solvable, without any free parameter, conceptually and quantitatively, within the approach of modified weak-field General Relativity involving the cosmological constant $\Lambda$. That approach enables one to describe in a unified picture both the dynamics of dark matter containing galaxies and the accelerated expansion of the Universe, thus defining a {\it local} Hubble constant of a local flow and the {\it global} one. The data on the dark matter content of peculiar galaxy samples are shown to be compatible to that unified picture. Future more refined surveys of galaxy distribution, hierarchical dynamics and flows within the vicinity of the Local group and the Virgo supercluster can be decisive in revealing the possible common nature of the dark sector.