Papers for Friday, Dec 03 2021

We quantify the spatial distribution of stars for two subclusters centred around the massive/intermediate mass stars S Mon and IRS1/2 in the NGC2264 star-forming region. We find that both subclusters are have neither a substructured, nor a centrally concentrated distribution according to the Q-parameter. Neither subcluster displays mass segregation according to the $\Lambda_{\rm MSR}$ ratio, but the most massive stars in IRS1/2 have higher relative surface densities according to the $\Sigma_{\rm LDR}$ ratio. We then compare these quantities to the results of N-body simulations to constrain the initial conditions of NGC2264, which are consistent with having been dense ($\tilde{\rho} \sim 10^4$M$_\odot$pc$^{-3}$), highly substructured and subvirial. These initial conditions were also derived from a separate analysis of the runaway and walkaway stars in the region, and indicate that star-forming regions within 1kpc of the Sun likely have a broad range of initial stellar densities. In the case of NGC2264, its initial stellar density could have been high enough to cause the destruction or truncation of protoplanetary discs and fledgling planetary systems due to dynamical encounters between stars in the early stages of its evolution.

Making use of the publicly available 1D photoionization hydrodynamics code ATES we set out to investigate the combined effects of planetary gravitational potential energy ($\phi_p\equiv GM_p/R_p$) and stellar XUV irradiation ($F_{\rm XUV}$) on the evaporation efficiency ($\eta$) of moderately-to-highly irradiated gaseous planets, from sub-Neptunes through hot Jupiters. We show that the (known) existence of a threshold potential above which energy-limited escape (i.e., $\eta\simeq 1$) is unattainable can be inferred analytically. For $\phi_p\gtrsim \phi_p^{\rm thr}\approx [7.7-14.4]\times 10^{12}$ erg g$^{-1}$, most of the energy absorption occurs where the average kinetic energy acquired by the ions through photo-electron collisions is insufficient for escape. This causes the evaporation efficiency to plummet with increasing $\phi_p$,. Whether or not planets with $\phi_p\lesssim \phi_p^{\rm thr}$ exhibit energy-limited outflows is regulated primarily by the stellar irradiation level. Specifically, the following analytical threshold is derived: $F_{\rm XUV}^{\rm thr}\approx [2-8]\times 10^3 (g/g_\oplus)$ erg cm$^{-2}$ s$^{-1}$. Above this value, Ly$\alpha$ losses overtake adiabatic cooling and the evaporation efficiency of low-gravity planets, too, drops below the energy-limited approximation, albeit remaining largely independent of $\phi_p$. Further, we show that whereas $\eta$ increases as $F_{\rm XUV}$ increases for planets above $\phi^{\rm thr}_p$, the opposite is true for low-gravity planets. This behaviour can be understood by examining the relative fractional contributions of adiabatic vs. Ly$\alpha$ cooling as a function of atmospheric temperature. This novel framework enables a reliable, physically motivated prediction of the expected evaporation efficiency for a given planetary system; an analytical approximation of the best-fitting $\eta$ is given in the Appendix.

The observed correlation between outer giant planets and inner super-Earths is emerging as an important constraint on planet formation theories. In this study we focus on Kepler-167, which is currently the only system known to contain both inner transiting super-Earths and a confirmed outer transiting gas giant companion beyond 1 au. Using long term radial velocity monitoring, we measure the mass of the gas giant Kepler-167e ($P=1071$ days) to be $1.01^{+0.16}_{-0.15}$ M$_{\rm J}$, thus confirming it as a Jupiter analog. We re-fit the $Kepler$ photometry to obtain updated radii for all four planets. Using a planetary structure model, we estimate that Kepler-167e contains $66\pm19$ M$_{\oplus}$ of solids and is significantly enriched in metals relative to its solar-metallicity host star. We use these new constraints to explore the broader question of how systems like Kepler-167 form in the pebble accretion framework for giant planet core formation. We utilize simple disk evolution models to demonstrate that more massive and metal-rich disks, which are the most favorable sites for giant planet formation, can also deliver enough solids to the inner disk to form systems of super-Earths. We use these same models to constrain the nature of Kepler-167's protoplanetary disk, and find that it likely contained $\gtrsim 300$ M$_{\oplus}$ of dust and was $\gtrsim 40$ au in size. These values overlap with the upper end of the observed dust mass and size distributions of Class 0 and I disks, and are also consistent with the observed occurrence rate of Jupiter analogs around sun-like stars.

Observations of structure at sub-galactic scales are crucial for probing the properties of dark matter, which is the dominant source of gravity in the universe. It will become increasingly important for future surveys to distinguish between line-of-sight halos and subhalos to avoid wrong inferences on the nature of dark matter. We reanalyze a sub-galactic structure (in lens JVAS B1938+666) that has been previously found using the gravitational imaging technique in galaxy-galaxy lensing systems. This structure has been assumed to be a satellite in the halo of the main lens galaxy. We fit the redshift of the perturber of the system as a free parameter, using the multi-plane thin-lens approximation, and find that the redshift of the perturber is $z_\mathrm{int} = 1.22\substack{+0.11 \\ -0.11}$ (with a main lens redshift of $z=0.881$). Our analysis indicates that this structure is more massive than the previous result by more than an order of magnitude. This constitutes the first dark perturber shown to be a line-of-sight halo with a gravitational lensing method.

We present AT2020mrf (SRGe J154754.2$+$443907), an extra-galactic ($z=0.1353$) fast-rising ($t_{g,\rm rise}=3.7$ days) luminous ($M_{g,\rm peak}=-20.0$) optical transient. Its optical spectrum around peak shows a broad ($v\sim0.1c$) emission feature on a blue continuum ($T\sim2\times10^4$ K), which bears a striking resemblance to AT2018cow. Its bright radio emission ($\nu L_\nu = 1.2\times 10^{39}\,{\rm erg\,s^{-1}}$; $\nu_{\rm rest}= 7.4$ GHz; 261 days) is similar to four other AT2018cow-like events, and can be explained by synchrotron radiation from the interaction between a sub-relativistic ($\gtrsim0.07$--$0.08c$) forward shock and a dense environment ($\dot M \lesssim 10^{-3}\,M_\odot \,{\rm yr^{-1}}$ for $v_{\rm w}=10^3\,{\rm km\,s^{-1}}$). AT2020mrf occurs in a galaxy with $M_\ast \sim 10^8\,M_\odot$ and specific star formation rate $\sim 10^{-10}\, {\rm yr^{-1}}$, supporting the idea that AT2018cow-like events are preferentially hosted by dwarf galaxies. The X-ray properties of AT2020mrf are unprecedented. At 35--37 days, SRG/eROSITA detected luminous ($L_{\rm X}\sim 2\times 10^{43}\,{\rm erg\,s^{-1}}$; 0.3--10 keV) X-ray emission. The X-ray spectral shape ($f_\nu \propto \nu^{-0.8}$) and erratic intraday variability are reminiscent of AT2018cow, but the luminosity is a factor of $\sim20$ greater than AT2018cow. At 328 days, Chandra detected it at $L_{\rm X}\sim10^{42}\,{\rm erg\,s^{-1}}$, which is $>200$ times more luminous than AT2018cow and CSS161010. At the same time, the X-ray emission remains variable on the timescale of $\sim1$ day. We show that a central engine, probably a millisecond magnetar or an accreting black hole, is required to power the explosion. We predict the rates at which events like AT2018cow and AT2020mrf will be detected by present and future X-ray all-sky surveys, including SRG and Einstein Probe.

We investigate the possibility of cosmic ray (CR) confinement by charged dust grains through resonant drag instabilities (RDIs). We perform magnetohydrodynamic particle-in-cell simulations of magnetized gas mixed with charged dust and cosmic rays, with the gyro-radii of dust and GeV CRs on $\sim\mathrm{AU}$ scales fully resolved. As a first study, we focus on one type of RDI wherein charged grains drift super-Alfv\'enically, with Lorentz forces strongly dominating over drag forces. Dust grains are unstable to the RDIs and form concentrated columns and sheets, whose scale grows until saturating at the simulation box size. Initially perfectly-streaming CRs are strongly scattered by RDI-excited Alfv\'en waves, with the growth rate of the CR perpendicular velocity components equaling the growth rate of magnetic field perturbations. These rates are well-predicted by analytic linear theory. CRs finally become isotropized and drift at least at $\sim v_\mathrm{A}$ by unidirectional Alfv\'{e}n waves excited by the RDIs, with a uniform distribution of the pitch angle cosine $\mu$ and a flat profile of the CR pitch angle diffusion coefficient $D_{\mu\mu}$ around $\mu = 0$, without the "$90$ degree pitch angle problem." With CR feedback on the gas included, $D_{\mu\mu}$ decreases by a factor of a few, indicating a lower CR scattering rate, which contradicts predictions of standard self-confinement theory. Our study demonstrates that the dust-induced CR confinement can be very important under certain conditions, e.g., the dusty circumgalactic medium around quasars or superluminous galaxies.

Context. The fraction of field binaries on retrograde orbits about the Milky Way is significantly lower compared to its prograde counterpart. Chemical and dynamical evidence suggests that the retrograde stellar population originates from $\omega$ Centauri, which is either the most massive globular cluster (GC) of the Milky Way or the putative core of a former dwarf galaxy. Aims. Star formation conditions required to produce the retrograde binary population are constrained assuming that the retrograde stellar population originates from $\omega$ Centauri's progenitor. Methods. We match the observed low binary fraction with dynamical population synthesis models, including a universal initial binary population and dynamical processing in star clusters, making use of the publicly available binary population synthesis tool BiPoS1. Results. It is found that either the GC progenitor of $\omega$ Centauri must have formed with a stellar density of $\approx 10^8 \; M_{sun} \; pc^{-3}$ or the $\omega$ Centauri dwarf galaxy's progenitor star cluster population must have formed in an extreme starburst with a star formation rate exceeding $1000 \; M_{sun} \; yr^{-1}$ and probably a top-heavy embedded-cluster mass function with suppressed low-mass cluster formation. The separation and mass ratio distribution for retrograde field binaries are predicted for comparison with future observations. Conclusions. A viable solution for the deficiency of binaries on retrograde orbits is presented, and star formation conditions for $\omega$ Centauri as well as orbital parameter distributions for the Milky Way's retrograde binary population are predicted. The dwarf galaxy origin for $\omega$ Centauri is tentatively preferred within the present context.

This is the fourth paper of a series investigating the AGN fuelling/feedback processes in a sample of eleven nearby low-excitation radio galaxies (LERGs). In this paper we present follow-up Atacama Large Millimeter/submillimeter Array (ALMA) observations of one source, NGC3100, targeting the $^{12}$CO(1-0), $^{12}$CO(3-2), HCO$^{+}$(4-3), SiO(3-2) and HNCO(6-5) molecular transitions. $^{12}$CO(1-0) and $^{12}$CO(3-2) lines are nicely detected and complement our previous $^{12}$CO(2-1) data. By comparing the relative strength of these three CO transitions, we find extreme gas excitation conditions (i.e. $T_{\rm ex}\gtrsim50$ K) in regions that are spatially correlated with the radio lobes, supporting the case for a jet-ISM interaction. An accurate study of the CO kinematics demonstrates that, although the bulk of the gas is regularly rotating, two distinct non-rotational kinematic components can be identified in the inner gas regions: one can be associated to inflow/outflow streaming motions induced by a two-armed spiral perturbation; the second one is consistent with a jet-induced outflow with $v_{\rm max}\approx 200$ km s$^{-1}$ and $\dot{M}\lesssim 0.12$ M$_{\odot}$ yr$^{-1}$. These values indicate that the jet-CO coupling ongoing in NGC3100 is only mildly affecting the gas kinematics, as opposed to what expected from existing simulations and other observational studies of (sub-)kpc scale jet-cold gas interactions. HCO$^{+}$(4-3) emission is tentatively detected in a small area adjacent to the base of the northern radio lobe, possibly tracing a region of jet-induced gas compression. The SiO(3-2) and HNCO(6-5) shock tracers are undetected: this - along with the tentative HCO$^{+}$(4-3) detection - may be consistent with a deficiency of very dense (i.e. $n_{\rm crit} > 10^{6}$ cm$^{-3}$) cold gas in the central regions of NGC3100.

We present a case study of stellar clumps in G04-1, a clumpy, turbulent disk galaxy located at $z$ = 0.13 from the DYNAMO sample, using adaptive optics enabled K-band imaging ($\sim2.25$ kpc/arcsec) with Keck/NIRC2. We identify 15 stellar clumps in G04-1 with a range of masses from $3.6 \times 10^{6}$ to $2.7 \times 10^{8}\ \rm M_{\odot}$, and with a median mass of $\sim2.9 \times 10^{7}\ \rm M_{\odot}$. Note that these masses decrease by about one-half when we apply a light correction for the underlying stellar disk. A majority (12 of 15) of clumps observed in the $K_{P}$-band imaging have associated components in H$\alpha$ maps ($\sim2.75$ kpc/arcsec; $<$R$_{clump}> \sim$500 pc) and appear co-located ($\overline{\Delta x} \sim 0.1$ arcsec). Using Hubble Space Telescope WFC/ACS observations with the F336W and F467M filters, we also find evidence of radial trends in clump stellar properties: clumps closer to the centre of G04-1 are more massive (consistent with observations at high-$z$) and appear more red, suggesting they may be more evolved. Using our high-resolution data, we construct a star forming main sequence for G04-1 in terms of spatially-resolved quantities and find that all regions (both clump and intra-clump) within the galaxy are experiencing an enhanced mode of star formation routinely observed in galaxies at high-$z$. In comparison to recent simulations, our observation of a number of clumps with masses $10^{7}-10^{8}\ \rm M_{\odot}$ is not consistent with strong radiative feedback in this galaxy.

The velocity-space distribution of the solar neighbourhood stars shows complex substructures. Most of the previous studies use static potentials to investigate their origins. Instead we use a self-consistent $N$-body model of the Milky Way, whose potential is asymmetric and evolves with time. In this paper, we quantitatively evaluate the similarities of the velocity-space distributions in the $N$-body model and that of the solar neighbourhood, using Kullback-Leibler divergence (KLD). The KLD analysis shows the time evolution and spatial variation of the velocity-space distribution. The KLD fluctuates with time, which indicates the velocity-space distribution at a fixed position is not always similar to that of the solar neighbourhood. Some positions show velocity-space distributions with small KLDs (high similarities) more frequently than others. One of them locates at $(R,\phi)=(8.2\;\mathrm{kpc}, 30^{\circ})$, where $R$ and $\phi$ are the distance from the galactic centre and the angle with respect to the bar's major axis, respectively. The detection frequency is higher in the inter-arm regions than in the arm regions. In the velocity maps with small KLDs, we identify the velocity-space substructures, which consist of particles trapped in bar resonances. The bar resonances have significant impact on the stellar velocity-space distribution even though the galactic potential is not static.

We investigate $r$-process nucleosynthesis and kilonova emission resulting from binary neutron star (BNS) mergers based on a three-dimensional (3D) general-relativistic magnetohydrodynamic (GRMHD) simulation of a hypermassive neutron star (HMNS) remnant. The simulation includes a microphysical finite-temperature equation of state (EOS) and neutrino emission and absorption effects via a leakage scheme. We track the thermodynamic properties of the ejecta using Lagrangian tracer particles and determine its composition using the nuclear reaction network $\texttt{SkyNet}$. We investigate the impact of neutrinos on the nucleosynthetic yields by varying the neutrino luminosities during post-processing. The ejecta show a broad distribution with respect to their electron fraction $Y_e$, peaking between $\sim$0.25-0.4 depending on the neutrino luminosity employed. We find that the resulting $r$-process abundance patterns differ from solar, with no significant production of material beyond the second $r$-process peak when using luminosities recorded by the tracer particles. We also map the HMNS outflows to the radiation hydrodynamics code $\texttt{SNEC}$ and predict the evolution of the bolometric luminosity as well as broadband light curves of the kilonova. The bolometric light curve peaks on the timescale of a day and the brightest emission is seen in the infrared bands. This is the first direct calculation of the $r$-process yields and kilonova signal expected from HMNS winds based on 3D GRMHD simulations. For longer-lived remnants, these winds may be the dominant ejecta component producing the kilonova emission.

The dynamics of the intracluster medium (ICM) is affected by turbulence driven by several processes, such as mergers, accretion and feedback from active galactic nuclei. X-ray surface brightness fluctuations have been used to constrain turbulence in galaxy clusters. Here, we use simulations to further investigate the relation between gas density and turbulent velocity fluctuations, with a focus on the effect of the stratification of the ICM. In this work, we studied the turbulence driven by hierarchical accretion by analysing a sample of galaxy clusters simulated with the cosmological code ENZO. We used a fixed scale filtering approach to disentangle laminar from turbulent flows. In dynamically perturbed galaxy clusters, we found a relation between the root mean square of density and velocity fluctuations, albeit with a different slope than previously reported. The Richardson number is a parameter that represents the ratio between turbulence and buoyancy, and we found that this variable has a strong dependence on the filtering scale. However, we could not detect any strong relation between the Richardson number and the logarithmic density fluctuations, in contrast to results by recent and more idealised simulations. In particular, we find a strong effect from radial accretion, which appears to be the main driver for the gas fluctuations. The ubiquitous radial bias in the dynamics of the ICM suggests that homogeneity and isotropy are not always valid assumptions, even if the turbulent spectra follow Kolmogorov's scaling. Finally, we find that the slope of the velocity and density spectra are independent of cluster-centric radii.

All of the known circumbinary planets are large (> 3 Earth radii). Whilst observational biases may account for this dearth of small planets, in this paper we propose a theoretical explanation. Most of the known planets are near the stability limit, interspersed between potentially unstable 5 : 1, 6 : 1 and 7 : 1 mean motion resonances with the binary. It is believed that these planets did not form in situ, but rather migrated from farther out in the disc, and hence passed through these resonances. Planets are expected to migrate at a speed proportional to their mass, and a slower rate makes resonant capture and subsequent ejection more likely. Therefore, whilst large planets may be able to successfully "run the gauntlet", small planets may be imperiled. This hypothesis is tested using N-body integrations of migration in a truncated and turbulent disc. We discover that surprisingly none of the known planets exist interior to a fully unstable resonance. We demonstrate that resonant ejection of migrating planets may occur in nature, and that it does indeed disproportionately affect small planets, but it requires a highly turbulent disc and its efficiency is likely too low to fully explain a dearth of < 3 Earth radii planets.

This program is part of QUEST (Quasar/ULIRG Evolutionary Study) and seeks to examine the gaseous environments of z < 0.3 quasars and ULIRGs as a function of host galaxy properties and age across the merger sequence from ULIRGs to quasars. This first paper in the series focuses on 33 quasars from the QUEST sample and on the kinematics of the highly ionized gas phase traced by the N V 1238, 1243 A and O VI 1032, 1038 A absorption lines in high-quality Hubble Space Telescope (HST) Cosmic Origins Spectrograph (COS) data. N V and O VI outflows are present in about 60% of the QUEST quasars and span a broad range of properties, both in terms of equivalent widths (from 20 mA to 25 A) and kinematics (outflow velocities from a few x 100 km/s up to ~10,000 km/s). The rate of incidence and equivalent widths of the highly ionized outflows are higher among X-ray weak or absorbed sources. The weighted outflow velocity dispersions are highest among the X-ray weakest sources. No significant trends are found between the weighted outflow velocities and the properties of the quasars and host galaxies although this may be due to the limited dynamic range of properties of the current sample. These results will be re-examined in an upcoming paper where the sample is expanded to include the QUEST ULIRGs. Finally, a lower limit of ~0.1% on the ratio of time-averaged kinetic power to bolometric luminosity is estimated in the 2-4 objects with blueshifted P V 1117, 1128 absorption features.

We investigate the temporal and colour variability of 897 blazars, comprising 455 BL Lacertae objects (BL Lacs) and 442 Flat Spectrum Radio Quasars (FSRQs), selected from the Roma-BZCAT catalogue, using the multi-band light curves from the Zwicky Transient Facility (ZTF DR6) survey. Assessing the colour variability characteristics over ~2 year timescales, we found that 18.5 per cent (84 out of 455) BL Lacs showed a stronger bluer when brighter (BWB) trend, whereas 9.0 per cent (41 out of 455) showed a redder when brighter (RWB) trend. The majority (70 per cent) of the BL Lacs showing RWB are host galaxy dominated. For the FSRQ subclass, 10.2 per cent (45 out of 442) objects showed a strong BWB trend and 17.6 per cent (78 out of 442) showed a strong RWB trend. Hence we find that BL Lacs more commonly follow a BWB trend than do FSRQs. This can be attributed to the more dominant jet emission in the case of BL Lacs and the contribution of thermal emission from the accretion disc for FSRQs. In analysing the colour behaviour on shorter time windows, we find many blazars evince shorter partial trends of BWB or RWB nature (or occasionally both). Some of such complex colour behaviours observed in the colour-magnitude diagrams of the blazars may result from transitions between the jet-dominated state to the disc-dominated state and vice versa.

Experiments deploying large arrays of transition-edge sensors (TESs) often require a robust method to monitor gain variations with minimal loss of observing time. We propose a sensitive and non-intrusive method for monitoring variations in TES responsivity using small square waves applied to the TES bias. We construct an estimator for a TES's small-signal power response from its electrical response that is exact in the limit of strong electrothermal feedback. We discuss the application and validation of this method using flight data from SPIDER, a balloon-borne telescope that observes the polarization of the cosmic microwave background with more than 2000 TESs. This method may prove useful for future balloon- and space-based instruments, where observing time and ground control bandwidth are limited.

The unknown physical nature of the Dark Energy motivates in cosmology the study of modifications of the gravity theory at large distances. One of these types of modifications are the theories known as $f(R)$ gravity. In this paper we use observational data to both constraint and test the Starobinsky $f(R)$ model \cite{Starobinsky2007}, using updated measurements from the dynamics of the expansion of the universe, $H(z)$, and the growth rate of cosmic structures, $[f\sigma_8](z)$, where the distinction between the concordance $\Lambda $CDM model and modified gravity models, $f(R)$, become clearer. We use MCMC likelihood analyses to explore the parameters space of the $f(R)$ model using $H(z)$ and $[f\sigma_8](z)$ data, both individually and jointly, and further check which of the models better fit the joint data. To further test the Starobinsky model, we use a method proposed by Linder \cite{Linder2017}, where the data from the observables is jointly binned in redshift space. This allows to further explore the model's parameter that better fits the data in comparison to the $\Lambda$CDM model. Our analyses show that the result from the MCMC run using $[f\sigma_8](z)$ data alone fits the analyzed model better than using the joint data. At the end, we reaffirm that this joint analysis is able to break the degenerescence between modified gravity models through the $f\sigma_8$ observable, as the original author proposed; our results show that the $f(R)$ Starobinsky model although not favored by the available data, cannot be discarded and is a possible alternative to the $\Lambda$CDM model.

Here we present an ensemble study of spin-orbit alignment in 51 close double star systems. We determine spin-orbit angles, obliquities, in 39 of these systems making use of recently improved apsidal motion rate measurements and apsidal motion constants. In the remaining 12 systems researchers have constrained spin-orbit alignment by different combinations of measurements of apsidal motion rates, projected obliquities and stellar inclinations. Of the 51 systems 48 are consistent with alignment albeit with some measurements having large uncertainties. A Fisher distribution with mean zero and a concentration factor $\kappa = 6.7$ represents this ensemble well. Indeed employing a bootstrapping resampling technique we find our data on these 48 systems is consistent with perfect alignment. We also confirm significant misalignment in two systems which travel on eccentric orbits and where misalignments have been reported on before, namely DI Her and AS Cam. The third misaligned system CV Vel orbits on a circular orbit. So while there are some glaring exceptions, the majority of close double star systems for which data are available appear to be well aligned.

The Laser Interferometer Space Antenna (LISA) is a future space-based gravitational wave (GW) detector designed to be sensitive to sources radiating in the low frequency regime (0.1 mHz to 1 Hz). LISA's interferometer signals will be dominated by laser frequency noise which has to be suppressed by about 7 orders of magnitude using an algorithm called Time-Delay Interferometry (TDI). Arm locking has been proposed to reduce the laser frequency noise by a few orders of magnitude to reduce the potential risks associated with TDI. In this paper, we present an updated performance model for arm locking for the new LISA mission using 2.5 Gm arm lengths, the currently assumed clock noise, spacecraft motion, and shot noise. We also update the Doppler frequency pulling estimates during lock acquisition.

We present a catalog of unbound stellar pairs, within 100 pc of the Sun, that are undergoing close, hyperbolic, encounters. The data are drawn from the GAIA EDR3 catalogue, and the limiting factors are errors in the radial distance and unknown velocities along the line of sight. Such stellar pairs have been suggested (Hansen & Zuckerman 2021) to be possible events associated with the migration of technological civilisations between stars. As such, this sample may represent a finite set of targets for a SETI search based on this hypothesis. Our catalog contains a total of 132 close passage events, featuring stars from across the entire main sequence, with 16 pairs featuring at least one main sequence star of spectral type between K1 and F3. Many of these stars are also in binaries, so that we isolate eight single stars as the most likely candidates to search for an ongoing migration event -- HD 87978, HD 92577, HD 50669, HD 44006, HD 80790, LSPM J2126+5338, LSPM J0646+1829 and HD 192486. Amongst host stars of known planets, the stars GJ 433 and HR 858 are the best candidates.

We present baseband radio observations of the millisecond pulsar J1909$-$3744, the most precisely timed pulsar, using the MeerKAT telescope as part of the MeerTime pulsar timing array campaign. During a particularly bright scintillation event the pulsar showed strong evidence of pulse mode changing, among the first millisecond pulsars and the shortest duty cycle millisecond pulsar to do so. Two modes appear to be present, with the weak (lower signal-to-noise ratio) mode arriving $9.26$ $\pm 3.94$ $\mu$s earlier than the strong counterpart. Further, we present a new value of the jitter noise for this pulsar of $8.20 \pm 0.14$ ns in one hour, finding it to be consistent with previous measurements taken with the MeerKAT ($9 \pm 3$ ns) and Parkes ($8.6 \pm 0.8$ ns) telescopes, but inconsistent with the previously most precise measurement taken with the Green Bank telescope ($14 \pm 0.5$ ns). Timing analysis on the individual modes is carried out for this pulsar, and we find an approximate $10\%$ improvement in the timing precision is achievable through timing the strong mode only as opposed to the full sample of pulses. By forming a model of the average pulse from templates of the two modes, we time them simultaneously and demonstrate that this timing improvement can also be achieved in regular timing observations. We discuss the impact an improvement of this degree on this pulsar would have on searches for the stochastic gravitational wave background, as well as the impact of a similar improvement on all MeerTime PTA pulsars.

Stations of dipole antennas for SKA1-Low will comprise 256 elements spread over an area with a diameter of 38~m. We consider the effect of residual unsubtracted sources well outside of the main beam for differing numbers of unique station configurations, in the Epoch of Reionisation (EoR) and the Cosmic Dawn (CD), both in simulation and with theoretical considerations. We find that beam sidelobes imprint power that renders the cosmological signal unobservable over a range of scales unless compact sources are subtracted beyond \theta_Z=30 degrees from zenith, and that station apodization will likely be required to control far sidelobes. An array with N_b=4 unique station configurations is sufficient to reduce the contamination, with an increase to N_b=8 showing little improvement. Comparison with an achromatic Airy disk beam model shows that beam sidelobe level is the main contributor to excess power in the EoR window, and beam chromaticity is less relevant. In the EoR, z=8.5, subtracting sources above 200~mJy out to \theta_Z=45 degrees, will be required to access relevant modes of the power spectrum.

We present the Engineering Development Array 2, which is one of two instruments built as a second generation prototype station for the future Square Kilometre Array Low Frequency Array. The array is comprised of 256 dual-polarization dipole antennas that can work as a phased array or as a standalone interferometer. We describe the design of the array and the details of design changes from previous generation instruments, as well as the motivation for the changes. Using the array as an imaging interferometer, we measure the sensitivity of the array at five frequencies ranging from 70 to 320 MHz.

Recently, FRB 190520B with the largest extragalactic dispersion measure (DM), was discovered by the Five-hundred-meter Aperture Spherical radio Telescope (FAST). The DM excess over the intergalactic medium and Galactic contributions is estimated as $\sim 900$ pc cm$^{-3}$, which is nearly ten times higher than other fast radio bursts (FRBs) host galaxies. The DM decreases with the rate $\sim0.1$ pc cm$^{-3}$ per day. It is the second FRB associated with a compact persistent radio source (PRS). The rotation measure (RM) is found to be larger than $1.8 \times 10^{5} \mathrm{rad} ~\mathrm{m}^{-2}$. In this letter, we argue that FRB 190520B is powered by a young magentar formed by core-collapse of massive stars, embedded in a composite of magnetar wind nebula (MWN) and supernova remnant (SNR). The energy injection of the magnetar drives the MWN and SN ejecta to evolve together, and the PRS is generated by the synchrotron radiation of the MWN. The magnetar has the interior magnetic field $B_{\text{int}}\sim (2-4)\times 10^{16}$ G and the age $t_{\text{age}}\sim 14-22$ yr. The dense SN ejecta and the shocked shell contribute a large fraction of the observed DM and RM. Our model can naturally explain the luminous PRS, decreasing DM and extreme RM of FRB 190520B simultaneously.

This invited white paper, submitted to the National Science Foundation in January of 2020, discusses the current challenges faced by the United States astronomical instrumentation community in the era of extremely large telescopes. Some details may have changed since submission, but the basic tenets are still very much valid. The paper summarizes the technical, funding, and personnel challenges the US community faces, provides an informal census of current instrumentation groups in the US, and compares the state-of-affairs in the US with that of the European community, which builds astronomical instruments from consortia of large hard-money funded instrument centers in a coordinated fashion. With the recent release of the Decadal Survey on Astronomy and Astrophysics 2020 (Astro2020), it is clear that strong community support exists for this next generation of large telescopes in the US. Is the US ready? Is there sufficient talent, facilities, and resources in the community today to meet the challenge of developing the complex suite of instruments envisioned for two US ELTs? These questions are addressed, along with thoughts on how the National Science Foundation can help build a more viable and stable instrumentation program in the US. These thoughts are intended to serve as a starting point for a broader discussion, with the end goal being a plan that puts the US astronomical instrumentation community on solid footing and poised to take on the challenges presented by the ambitious goals we have set in the era of ELTs.

The growing catalog of circumbinary planets strengthens the notion that planets form in a diverse range of conditions across the cosmos. Transiting circumbinary planets yield especially important insights and many examples are now known, in broadly coplanar obits with respect to their binary. Studies of circumbinary disks suggest misaligned transiting examples could also plausibly exist, but their existence would exacerbate the already challenging feat of automatic detection. In this work, we synthesize populations of such planets and consider the number of transits per epoch they produce, forming integer sequences. For isotropic distributions, such sequences will appear foreign to conventional expectation, rarely (~3%) producing the signature double-transits we have come to expect for circumbinaries, instead producing sparse sequences dominated by zero-transit epochs (~70%). Despite their strangeness, we demonstrate that these sequences will be non-random and that the two preceding epochs predict the next to high accuracy. Additionally, we show that even when clustering the transits into grouped epochs, they often appear unphysical if erroneously assuming a single star, due to the missing epochs. Crucially, missing epochs mean highly isotropic populations can trick the observer into assigning the wrong period in up to a quarter of cases, adding further confusion. Finally, we show that the transit sequences encode the inclination distribution and demonstrate a simple inference method that successfully matches the injected truth. Our work highlights how the simple act of flagging transits can be used to provide an initial, vetting-level analysis of misaligned transiting circumbinary planets.

We developed an operational solar flare prediction model using deep neural networks, named Deep Flare Net (DeFN). DeFN can issue probabilistic forecasts of solar flares in two categories, such as >=M-class and <M-class events or >=C-class and <C-class events, occurring in the next 24 h after observations and the maximum class of flares occurring in the next 24 h. DeFN is set to run every 6 h and has been operated since January 2019. The input database of solar observation images taken by the Solar Dynamic Observatory (SDO) is downloaded from the data archive operated by the Joint Science Operations Center (JSOC) of Stanford University. Active regions are automatically detected from magnetograms, and 79 features are extracted from each region nearly in real time using multiwavelength observation data. Flare labels are attached to the feature database, and then, the database is standardized and input into DeFN for prediction. DeFN was pretrained using the datasets obtained from 2010 to 2015. The model was evaluated with the skill score of the true skill statistics (TSS) and achieved predictions with TSS = 0.80 for >=M-class flares and TSS = 0.63 for >=C-class flares. For comparison, we evaluated the operationally forecast results from January 2019 to June 2020. We found that operational DeFN forecasts achieved TSS = 0.70 (0.84) for >=C-class flares with the probability threshold of 50 (40)%, although there were very few M-class flares during this period and we should continue monitoring the results for a longer time. Here, we adopted a chronological split to divide the database into two for training and testing. The chronological split appears suitable for evaluating operational models. Furthermore, we proposed the use of time-series cross-validation. The procedure achieved TSS = 0.70 for >=M-class flares and 0.59 for >=C-class flares using the datasets obtained from 2010 to 2017.

The ARIANNA experiment is an Askaryan detector designed to record radio signals induced by neutrino interactions in the Antarctic ice. Because of the low neutrino flux at high energies ($E > 10^{16} $), the physics output is limited by statistics. Hence, an increase in sensitivity significantly improves the interpretation of data and offers the ability to probe new parameter spaces. The amplitudes of the trigger threshold are limited by the rate of triggering on unavoidable thermal noise fluctuations. We present a real-time thermal noise rejection algorithm that enables the trigger thresholds to be lowered, which increases the sensitivity to neutrinos by up to a factor of two (depending on energy) compared to the current ARIANNA capabilities. A deep learning discriminator, based on a Convolutional Neural Network (CNN), is implemented to identify and remove thermal events in real time. We describe a CNN trained on MC data that runs on the current ARIANNA microcomputer and retains 95 percent of the neutrino signal at a thermal noise rejection factor of $10^5$, compared to a template matching procedure which reaches only $10^2$ for the same signal efficiency. Then the results are verified in a lab measurement by feeding in generated neutrino-like signal pulses and thermal noise directly into the ARIANNA data acquisition system. Lastly, the same CNN is used to classify cosmic-rays events to make sure they are not rejected. The network classified 102 out of 104 cosmic-ray events as signal.

The abundance of volatile sulfur in dense clouds is an old problem in studies of the physics and chemistry of star forming regions. Sulfur is an important species because its low ionization potential makes it a potentially important charge carrier. Observed sulfur-bearing species in the gas-phase of dense clouds represent only a minor fraction of the cosmic sulfur abundance, which has been interpreted as a signature of sulfur depletion into ices at the surface of dust grains. Yet, atomic sulfur, which could be the main gas-phase carrier, can not be observed directly in cold cores. We present measurements of the NS radical toward four dense cores performed with the IRAM-30m telescope. Analytical chemical considerations and chemical models over a wide parameter space show that the NS:N2H+ abundance ratio provides a direct constraint on the abundance of gas-phase atomic sulfur. Toward early-type cores, we find that n(S)/nH is close, or even equal, to the cosmic abundance of sulfur, 14x10^-6 , demonstrating that sulfur is not depleted and is atomic, in agreement with chemical models. More chemically evolved cores show sulfur depletion by factors up to 100 in their densest parts. In L1544, atomic sulfur depletion is shown to increase with increasing density. Future observations are needed to discover the solid-phase carrier of sulfur. The initial steps of the collapse of pre-stellar cores in the high sulfur abundance regime will also need to be explored from chemical and dynamical perspectives.

Observation of redshifted 21-cm signal from the Epoch of Reionization (EoR) is challenging due to contamination from the bright foreground sources that exceed the signal by several orders of magnitude. The removal of this very high foreground relies on accurate calibration to keep the intrinsic property of the foreground with frequency. Commonly employed calibration techniques for these experiments are the sky model-based and the redundant baseline-based calibration approaches. However, the sky model-based and redundant baseline-based calibration methods could suffer from sky-modeling error and array redundancy imperfection issues, respectively. In this work, we introduce the hybrid correlation calibration ("CorrCal") scheme, which aims to bridge the gap between redundant and sky-based calibration by relaxing redundancy of the array and including sky information into the calibration formalisms. We demonstrate the slight improvement of power spectra, about $-6\%$ deviation at the bin right on the horizon limit of the foreground wedge-like structure, relative to the power spectra before the implementation of "CorrCal" to the data from the Precision Array for Probing the Epoch of Reionization (PAPER) experiment, which was otherwise calibrated using redundant baseline calibration. This small improvement of the foreground power spectra around the wedge limit could be suggestive of reduced spectral structure in the data after "CorrCal" calibration, which lays the foundation for future improvement of the calibration algorithm and implementation method.

OB star clusters originate from parsec-scale massive molecular clumps. We aim to understand the evolution of temperature and density structures on the intermediate-scale ($\lesssim$0.1-1 pc) extended gas of massive clumps. We performed $\sim$0.1 pc resolution observations (SMA+APEX) of multiple molecular line tracers (e.g., CH$_{3}$CCH, H$_{2}$CS, CH$_{3}$CN, CH$_{3}$OH) which cover a wide range of excitation conditions, towards a sample of eight massive clumps. Based on various radiative transfer models, we constrain the gas temperature and density structures and establish an evolutionary picture, aided by a spatially-dependent virial analysis and abundance ratios of multiple species. We determine temperature radial profiles varying between 30-200 K over a continuous scale, from the center of the clumps out to 0.3-0.4 pc radii. The clumps' radial gas density profiles, described by radial power-laws with slopes between -0.6 and $\sim$-1.5, are steeper for more evolved sources, as suggested by results based on both dust continuum, representing the bulk of the gas ($\sim$10$^{4}$ cm$^{-3}$), and CH$_{3}$OH lines probing the dense gas ($\gtrsim$10$^{6}$-10$^{8}$ cm$^{-3}$) regime. The density contrast between the dense gas and the bulk gas increases with evolution, and may be indicative of spatially and temporally varying star formation efficiencies. The radial profiles of the virial parameter show a global variation towards a sub-virial state as the clump evolves. The line-widths decline with increasing radius around the central core region and increase in the outer envelope, with a slope shallower than the case of the supersonic turbulence ($\,\propto\,$$r^{0.5}$) and the subsonic Kolmogorov scaling ($\,\propto\,$$r^{0.33}$). In the context of clump evolution, we also find that the abundance ratios of [CCH]/[CH$_{3}$OH] and [CH$_{3}$CN]/[CH$_{3}$OH] show correlations with clump $L/M$.

We study the impact of young clusters on the gravitational wave background from compact binary coalescence. We simulate a catalog of sources from population I/II isolated binary stars and stars born in young clusters, corresponding to one year of observations with second-generation (2G) detectors. Taking into account uncertainties on the fraction of dynamical binaries and star formation parameters, we find that the background is dominated by the population of binary black holes, and we obtain a value of $\Omega_{gw}(25 \rm{Hz}) = 1.2^{+1.38}_{-0.65} \times 10^{-9}$ for the energy density, in agreement with the actual upper limits derived from the latest observation run of LIGO--Virgo. We demonstrate that a large number of sources in a specific corrected mass range yields to a bump in the background. This background could be detected with 8 years of coincident data by a network of 2G detectors.

The signs of an expanding atmosphere of HD209458b have been observed with far-ultraviolet transmission spectroscopy and in the measurements of transit absorption by metastable HeI. These observations are interpreted using the hydrodynamic and Monte-Carlo numerical simulations of various degree of complexity and consistency. At the same time, no attempt has been made to model atmospheric escape of a magnetized HD209458b, to see how the planetary magnetic field might affect the measured transit absorption lines. This paper presents the global 3D MHD self-consistent simulations of the expanding upper atmosphere of HD209458b interacting with the stellar wind, and models the observed HI (Lya), OI (1306 A), CII (1337 A), and HeI (10830 A) transit absorption features. We find that the planetary dipole magnetic field with the equatorial surface value of Bp = 1 G profoundly changes the character of atmospheric material outflow and the related absorption. We also investigate the formation of planetary magnetosphere in the stellar wind and show that its size is more determined by the escaping atmosphere flow rather than by the strength of magnetic field. Fitting of the simulation results to observations enables constraining the stellar XUV flux and He abundance at ~10 erg cm2/s (at 1 a.u.) and He/H=0.02, respectively, as well as setting an upper limit for the dipole magnetic field of Bp<0.1 G on the planetary surface at the equator. This implies that the magnetic dipole moment of HD209458b should be less than 0.06 of the Jovian value.

Neural networks have provided powerful approaches to solve various scientific problems. Many of them are even difficult for human experts who are good at accessing the physical laws from experimental data. We investigate whether neural networks can assist us in exploring the fundamental laws of classical mechanics from data of planetary motion. Firstly, we predict the orbits of planets in the geocentric system using the gate recurrent unit, one of the common neural networks. We find that the precision of the prediction is obviously improved when the information of the Sun is included in the training set. This result implies that the Sun is particularly important in the geocentric system without any prior knowledge, which inspires us to gain Copernicus' heliocentric theory. Secondly, we turn to the heliocentric system and make successfully mutual predictions between the position and velocity of planets. We hold that the successful prediction is due to the existence of enough conserved quantities (such as conservations of mechanical energy and angular momentum) in the system. Our research provides a new way to explore the existence of conserved quantities in mechanics system based on neural networks.

We report on the detection, for the first time in space, of the radical HCCCO and of pentacarbon monoxide, C5O. The derived column densities are (1.6+/-0.2)e11 cm-2 and (1.5+/-0.2)e10 cm-2, respectively. We have also analysed the data for all the molecular species of the families HCnO and CnO within our QUIJOTE line survey. Upper limits are obtained for HC4O, HC6O, C4O, and C6O. We report a robust detection of HC5O and HC7O based on 14 and 12 rotational lines detected with a signal-to-noise ratio >30 and >5,respectively. The derived N(HC3O)/N(HC5O) abundance ratio is 0.09+/-0.03, while N(C3O)/N(C5O) is 80+/-2, and N(HC5O)/N(HC7O) is 2.2+/-0.3. As opposed to the cyanopolyyne family, HC2n+1N, which shows a continuous decrease in the abundances with increasing n, the CnO and HCnO species show a clear abundance maximum for n=3 and 5, respectively. They also show an odd and even abundance alternation, with odd values of n being the most abundant, which is reminiscent of the behaviour of CnH radicals, where in that case species with even values of n are more abundant. We explored the formation of these species through two mechanisms previously proposed, which are based on radiative associations between CnHm+ ions with CO and reactions of Cn- and CnH- anions with O atoms, and we found that several species, such as C5O, HC4O, and HC6O, are significantly overestimated. Our understanding of how these species are formed is incomplete as of yet. Other routes based on neutral-neutral reactions such as those of Cn and CnH carbon chains with O, OH, or HCO, could be behind the formation of these species.

X-ray emission from counterparts of historical classical novae (CNe) in our Galaxy is studied. To this end, we use data from three SRG/eROSITA sky surveys in the hemisphere analyzed by the RU eROSITA consortium. Out of 309 historical CNe, X-ray emission has been detected from 52 sources with 0.3-2.3 keV luminosities in the $\rm L_X\approx 10^{30}\sim 10^{34}$ erg/s range. Among them, two sources have supersoft spectra and are associated with the post-nova supersoft X-ray emission. Hardness of X-ray spectra of some of the remaining sources tentatively suggests that magnetized white dwarfs (WDs) may account for some fraction of CN counterparts detected in X-rays. This hypothesis will be further tested in the course of the following SRG/eROSITA sky surveys. The CN counterparts represent a bona fide sample of accreting WDs with unstable hydrogen burning on their surface, while their X-ray luminosity in quiescence is a reasonable proxy for the accretion rate in the binary system. Using this fact, we have constructed the accretion rate distribution of WDs with unstable hydrogen burning and compared it with the accretion rate distribution of known steady supersoft X-ray sources in our Galaxy and nearby external galaxies. There is a pronounced dichotomy between these two distributions with the CN counterparts and the steady supersoft sources occupying different domains along the mass accretion rate axis, in accordance with the predictions of the theory of hydrogen burning on the WD surface.

Microturbulence (\xi) is a key parameter introduced in stellar spectroscopy to explain the strength of saturated lines by formally incorporating an additional thermal broadening term in the line opacity profile. Although our Sun can serve as an important testing bench to check the usual assumption of constant \xi, the detailed behavior of how \xi varies from the disk center through the limb seems to have never been investigated so far. In order to fill this gap, local \xi values on the solar disk were determined from the equivalent widths of 46 Fe I lines at 32 points from the center to the limb by requiring the consistency between the abundances derived from lines of various strengths. The run of \xi with \theta (angle between line of sight and the surface normal) was found to be only gradual from ~1.0km/s (at sin\theta = 0: disk center) to ~1.3km/s (at sin\theta ~ 0.7: two-thirds of radial distance); but thereafter increasing more steeply up to ~2km/s (at sin\theta = 0.97: limb). This result further suggests that the microturbulence derived from the flux spectrum of the disk-integrated Sun is by ~20% larger than that of the disk-center value, which is almost consistent with the prediction from 3D hydrodynamical model atmospheres.

An inconsistency between the theoretical analysis and numerical calculations of the relativistic $r$-modes puzzles the neutron star community since the Kojima's finding of the continuous part in the $r$-mode oscillation spectrum in 1997. In this paper, after a brief review of the Newtonian $r$-mode theory and of the literature devoted to the continuous spectrum of $r$-modes, we apply our original approach to the study of relativistic oscillation equations. Working within the Cowling approximation, we derive the general equations, governing the dynamics of discrete relativistic $r$-modes for both barotropic (isentropic) and nonbarotropic stars. A detailed analysis of the obtained equations in the limit of extremely slow stellar rotation rate reveals that, because of the effect of inertial reference frame-dragging, the relativistic $r$-mode eigenfunctions and eigenfrequencies become {\it non-analytic} functions of the stellar angular velocity, $\Omega$. We also derive the explicit expressions for the $r$-mode eigenfunctions and eigenfrequencies for very small values of $\Omega$. These expressions explain the asymptotic behavior of the numerically calculated eigenfrequencies and eigenfunctions in the limit $\Omega\to 0$. All the obtained $r$-mode eigenfrequencies take discrete values in the frequency range, usually associated with the continuous part of the spectrum. No indications of the continuous spectrum, at least in the vicinity of the Newtonian $l=m=2$ $r$-mode frequency $\sigma=-4/3 \ \Omega$, are found.

Different properties of quasars may be observed and analysed through the many ranges of the electromagnetic spectrum. Pioneering studies showed that an "H-R diagram" for quasars was needed to organize these data, and that more than two dimensions were necessary: a four dimensional Eigenvector (4DE1) parameter space was proposed. The 4DE1 makes use of independent observational properties obtained from the optical and UV emission lines, as well as from the soft-X rays. The 4DE1 "optical plane", also known as the quasar Main Sequence (MS), identifies different spectral types in order to describe a consistent picture of QSOs. In this work we present a spectroscopic analysis focused on the comparison between two sources, one radio-loud (PKS2000-330, z = 3.7899) and one radio-quiet (Q1410+096, z = 3.3240), both showing Population A quasar spectral properties. Optical spectra were observed in the infrared with VLT/ESO, and the additional measures in UV were obtained through the fitting of archive spectra. The analysis was performed through a non-linear multi-component decomposition of the emission line profiles. Results are shown in order to highlight the effects of the radio-loudness on their emission line properties. The two quasars share similar optical spectroscopic properties and are very close on the MS classification while presenting significant differences on the UV data. Both sources show significant blueshifts in the UV lines but important differences in their UV general behaviour. While the radio-quiet source Q1410+096 shows a typical Pop A UV spectrum with similar intensities and shapes on both CIV1549 and SiIV1392, the UV spectrum of the strong radio-loud PKS2000-330 closely resembles the one of population B of quasars.

The 4m International Liquid Mirror Telescope (ILMT) installation activities have recently been completed at the Devasthal observatory (Uttarakhand, India). The ILMT will perform continuous observation of a narrow strip of the sky ($\sim$27$'$) passing over the zenith in the SDSS $g'$, $r'$ and $i'$ bands. In combination with a highly efficient 4k $\times$ 4k CCD camera and an optical corrector, the images will be secured at the prime focus of the telescope using the Time Delayed Integration technique. The ILMT will reach $\sim$22.5 mag ($g'$-band) in a single scan and this limiting magnitude can be further improved by co-adding the nightly images. The uniqueness of the one-day cadence and deeper imaging with the ILMT will make it possible to discover and study various galactic and extra-galactic sources, especially variable ones. Here, we present the latest updates of the ILMT facility and discuss the preparation for the first light, which is expected during early 2022. We also briefly explain different steps involved in the ILMT data reduction pipeline.

The most spiral galaxies have a flat rotational velocity curve, according to the different observational techniques used in several wavelengths domain. In this work, we show that non-linear terms are able to balance the dispersive effect of the wave, thus reviving the observed rotational curve profiles without inclusion of any other but baryonic matter concentrated in the bulge and disk. In order to prove that the considered model is able to restore a flat rotational curve, Milky Way has been chosen as the best mapped galaxy to apply on. Using the gravitational N-body simulations with up to $10^7$ particles, we test this dynamical model in the case of the Milky Way with two different approaches. Within the direct approach, as an input condition in the simulation runs we set the spiral surface density distribution which is previously obtained as an explicit solution to non-linear Schr\"{o}dinger equation (instead of a widely used exponential disk approximation). In the evolutionary approach, we initialize the runs with different initial mass and rotational velocity distributions, in order to capture the natural formation of spiral arms, and to determine their role in the disk evolution. In both cases we are able to reproduce the stable and non-expanding disk structures at the simulation end times of $\sim10^9$ years, with no halo inclusion. Although the given model doesn't take into account the velocity dispersion of stars and finite disk thickness, the results presented here still imply that non-linear effects can significantly alter the amount of dark matter which is required to keep the galactic disk in stable configuration.

The large gap between a galactic dark matter subhalo's velocity and its own gravitational binding velocity creates the situation that dark matter soft-scattering on baryons to evaporate the subhalo, if kinetic energy transfer is efficient by low momentum exchange. Small subhalos can evaporate before dark matter thermalize with baryons due to the low binding velocity. In case dark matter acquires an electromagnetic dipole moment, the survival of low-mass subhalos requires stringent limits on the photon-mediated soft scattering. We calculate the subhalo evaporation rate via soft collision by ionization gas and accelerated cosmic rays, and show the stability of subhalos lighter than $10^{-5}M_{\odot}$ in the gaseous inner galactic region is sensitive to dark matter's effective electric and magnetic dipole moments below current direct detection limits.

LINC-NIRVANA (LN) is one of the instruments on-board the Large Binocular Telescope (LBT). LN is a high-resolution, near-infrared imager equipped with an advanced adaptive optics module. LN implements layer-oriented Multi-Conjugate Adaptive Optics (MCAO) approach using two independent wavefront sensors per side of the binocular telescope measuring the turbulence volume above the telescope. The capability of acquiring up to 20 Natural Guide Stars simultaneously from two distinct fields of view, and using them for wavefront sensing with 20 separate pyramids per side of the telescope makes the LN MCAO system one of a kind. Commissioning of the left MCAO channel is almost complete, while that of the right arm is on-going. The Science Verification on the left side is expected to start soon after the MCAO performance is optimised for faint guide stars. In this article, we put together the lessons learned during the commissioning of the LN MCAO module. We hope and believe that this article will help the future MCAO instrument commissioning teams.

We study the thermodynamics of galaxy clusters in a modified Newtonian potential motivated by a general solution to Newton's "sphere-point" equivalence theorem. We obtain the $N$ particle partition function by evaluating the configurational integral while accounting for the extended nature of galaxies (via the inclusion of the softening parameter $\epsilon$ into the potential energy function). This softening parameter takes care of the galaxy-halos whose effect on structuring the shape of the galactic disc has been found recently. The spatial distribution of the particles (galaxies) is also studied in this framework. A comparison of the new clustering parameter $b_+$ to the original clustering parameters is presented in order to visualize the effect of the modified gravity. We also discuss the possibility of system symmetry breaking via the behavior of the specific heat as a function of temperature.

We seek to differentiate dynamical and morphological attributes between globular clusters that were formed inside their own dark matter mini-halo, and those who were not. For that, we employ high resolution $N$-body simulations of globular clusters with (and without) an enveloping dark matter mini-halo, orbiting a host galaxy. We set the same prescriptions of the Fornax dwarf spheroidal galaxy and its main five globular clusters, and use $N$-body particles for all components (i.e., stars and dark matter, for both Fornax and its clusters). For clusters embedded in dark matter, we observe that the increment of mass from the extra dark component triggers a tidal radius growth that allows the mini-halo to work as a protective shield against tidal stripping, being itself stripped beforehand the stars. Consequently, tidal effects such as inflation of the stellar velocity dispersion, development of prominent tidal tails, ellipticity increase and diffusion of the stellar distribution profile are generally much milder in clusters originally embedded in dark matter. However, this shielding effect becomes negligible after an important amount of dark matter has been stripped, which happens faster for clusters having simultaneously short orbital periods, low typical orbital radii and relatively high eccentricities. Finally, we notice that even for clusters that retain a large amount of dark matter at redshift zero, their inner regions are still predominantly composed of stars, with the typical density ratio of dark matter to cluster stars being of the order of $1\%$ up to roughly $10~$pc away from the clusters' centre.

Scheduled to launch in late 2021, the Imaging X-ray Polarimetry Explorer (IXPE) is a NASA Small Explorer Mission in collaboration with the Italian Space Agency (ASI). The mission will open a new window of investigation - imaging X-ray polarimetry. The observatory features 3 identical telescopes each consisting of a mirror module assembly with a polarization-sensitive imaging X-ray detector at the focus. A coilable boom, deployed on orbit, provides the necessary 4-m focal length. The observatory utilizes a 3-axis-stabilized spacecraft which provides services such as power, attitude determination and control, commanding, and telemetry to the ground. During its 2-year baseline mission, IXPE will conduct precise polarimetry for samples of multiple categories of X-ray sources, with follow-on observations of selected targets.

Elongation of the point spread function due to atmospheric dispersion becomes a severe problem for high resolution imaging instruments, if an atmospheric dispersion corrector is not present. This work reports on a novel technique to measure this elongation, corrected or uncorrected, from imaging data. By employing a simple diffraction mask it is possible to magnify the chromatic elongation caused by the atmosphere and thus make it easier to measure. We discuss the theory and design of such a mask and report on two proof of concept observations using the 40 cm Gratama telescope of the University of Groningen. We evaluate the acquired images using a geometric approach, a forward modelling approach and from a direct measurement of the length of the point spread function. For the first two methods we report measurements consistent with atmospheric dispersion models to within 0.5 arcsec. Direct measurements of the elongation do not prove suitable for the characterisation of atmospheric dispersion. We conclude that the addition of this type of diffraction mask can be valuable for measurements of atmospheric dispersion enabling high precision correction of this atmospheric effect on future instruments.

Extracting information from stochastic fields or textures is a ubiquitous task in science, from exploratory data analysis to classification and parameter estimation. From physics to biology, it tends to be done either through a power spectrum analysis, which is often too limited, or the use of convolutional neural networks (CNNs), which require large training sets and lack interpretability. In this paper, we advocate for the use of the scattering transform (Mallat 2012), a powerful statistic which borrows mathematical ideas from CNNs but does not require any training, and is interpretable. We show that it provides a relatively compact set of summary statistics with visual interpretation and which carries most of the relevant information in a wide range of scientific applications. We present a non-technical introduction to this estimator and we argue that it can benefit data analysis, comparison to models and parameter inference in many fields of science. Interestingly, understanding the core operations of the scattering transform allows one to decipher many key aspects of the inner workings of CNNs.

Plasma loops or plumes rooted in sunspot umbrae often harbor downflows with speeds of 100 km/s. These downflows are supersonic at transition region temperatures of 0.1 MK. The source of these flows is not well understood. We aim to investigate the source of sunspot supersonic downflows (SSDs) in AR 12740 using simultaneous spectroscopic and imaging observations. We identified SSD events from multiple raster scans of a sunspot by the Interface Region Imaging Spectrograph, and calculated the electron densities, mass fluxes and velocities of these SSDs. The EUV images provided by the AIA onboard the SDO and the EUVI onboard the STEREO were employed to investigate the origin of these SSDs and their associated coronal rain. Almost all the identified SSDs appear at the footpoints of sunspot plumes and are temporally associated with appearance of chromospheric bright dots inside the sunspot umbra. Dual-perspective EUV imaging observations reveal a large-scale closed magnetic loop system spanning the sunspot region and a remote region. We observed that the SSDs are caused by repeated coronal rain that forms and flows along these closed magnetic loops toward the sunspot. One episode of coronal rain clearly indicates that reconnection near a coronal X-shaped structure first leads to the formation of a magnetic dip. Subsequently, hot coronal plasma catastrophically cools from 2 MK in the dip region via thermal instability. This results in the formation of a transient prominence in the dip, from which the cool gas mostly slides into the sunspot along inclined magnetic fields under the gravity. This drainage process manifests as a continuous rain flow, which lasts for around 2 hrs and concurrently results in a nearly steady SSD event. Our results demonstrate that coronal condensation in magnetic dips can result in the quasi-steady sunspot supersonic downflows.

The use of Bayesian neural networks is a novel approach for the classification of gamma-ray sources. We focus on the classification of Fermi-LAT blazar candidates, which can be divided into BL Lacertae objects and Flat Spectrum Radio Quasars. In contrast to conventional dense networks, Bayesian neural networks provide a reliable estimate of the uncertainty of the network predictions. We explore the correspondence between conventional and Bayesian neural networks and the effect of data augmentation. We find that Bayesian neural networks provide a robust classifier with reliable uncertainty estimates and are particularly well suited for classification problems that are based on comparatively small and imbalanced data sets. The results of our blazar candidate classification are valuable input for population studies aimed at constraining the blazar luminosity function and to guide future observational campaigns.

We study the anisotropy of phase transition gravitational wave (PTGW) and find it could be a new approach to probe the primordial density perturbation, especially at small scale. Generated long before recombination and free from Silk damping, PTGW could directly provide the information of density perturbation seeded from inflation or alternatives. The power spectrum of anisotropic PTGW even goes up at high $\ell$ due to the early integrated Sachs Wolfe effect, making it a powerful probe into these small scale perturbations. Moreover, we find anisotropic PTGW around nano Hertz could be detected by pulsar timing array projects like the Square Kilometre Array.

The analysis of debris disk observations is often based on the assumption of a dust phase composed of compact spherical grains consisting of astronomical silicate. Instead, observations indicate the existence of water ice in debris disks. We quantify the impact of water ice as a potential grain constituent in debris disks on the disk parameter values estimated from photometric and spatially resolved observations in the mid- and far-infrared. We simulated photometric measurements and radial profiles of debris disks containing water ice and analyzed them by applying a disk model purely consisting of astronomical silicate. Subsequently, we quantified the deviations between the derived and the true parameter values. As stars in central positions we discuss a $\beta$ Pic sibling and main-sequence stars with spectral types ranging from A0 to K5. To simulate observable quantities we employed selected observational scenarios regarding the choice of wavelengths and instrument characteristics. For the $\beta$ Pic stellar model and ice fractions $\geq 50\ \%$ the derived inner disk radius is biased by ice sublimation toward higher values. However, the derived slope of the radial density profile is mostly unaffected. Along with an increasing ice fraction, the slope of the grain size distribution is overestimated by up to a median factor of $\sim 1.2$ for an ice fraction of $90\ \%$ while the total disk mass is underestimated by a factor of $\sim 0.4$. The reliability of the derived minimum grain size strongly depends on the spectral type of the central star. For an A0-type star the minimum grain size can be underestimated by a factor of $\sim 0.2$, while for solar-like stars it is overestimated by up to a factor of $\sim 4 - 5$. Neglecting radial profile measurements and using solely photometric measurements, the factor of overestimation increases for solar-like stars up to $\sim 7 - 14$.

We compare some predictions of Wagoner & Tandon (2021) with the results of the hydrodynamic and magnetohydrodynamic (MHD) simulations of Reynolds & Miller (2009). It appears that the MHD simulations were not run for long enough and the numerical damping was not small enough to produce the observed high-frequency QPOs (and the g-mode seen in the hydro simulations).

There are currently many large-field surveys operational and planned including the powerful Vera Rubin Observatory Legacy Survey for Space and Time (LSST). These surveys will increase the number and diversity of transients dramatically. However, for some transients, like supernovae (SNe), we can gain more understanding by directed observations than by simply increasing the sample size. For example, the initial emission from these transients can be a powerful probe of these explosions. Upcoming ground-based detectors are not ideally suited to observe the initial emission (shock emergence) of these transients. These observations require a large field-of-view X-ray mission with a UV follow up within the first hour of the event. The emission in the first one hour to even one day provides strong constraints on the stellar radius and asymmetries in the outer layers of stars, the properties of the circumstellar medium (e.g. inhomogeneities in the wind for core-collapse SNe, accreting companion in thermonuclear SNe), and the transition region between these two. This paper describes a simulation for the number of SNe that could be seen by a large field of view lobster eye X-ray and UV observatory.

We use telemetry data from the Gemini North ALTAIR adaptive optics system to investigate how well the commands to the wavefront correctors (Tip/Tilt and high-order turbulence) can be forecasted in order to reduce lag error and improve delivered image quality. We show that a high level of reduction ($\sim$ 5 for Tip-Tilt and $\sim$ 2 for high-order modes) can be achieved by using a "forecasting module" based on a linear auto-regressive models with only a few coefficients ($\sim$ 30 for Tip-Tilt and $\sim$ 5 for high-order modes) to complement the existing integral servo-controller. Updating this module to adapt to evolving observing conditions is computationally inexpensive and requires less than 10 seconds worth of telemetry data. We also use several machine learning models to evaluate whether further improvements could be achieved with a more sophisticated non-linear model. Our attempts showed no perceptible improvements over linear auto-regressive predictions, even for large lags where residuals from the linear models are high, suggesting that non-linear wavefront distortions for ALTAIR at the Gemini North telescope may not be forecasted with the current setup.

Faraday tomography through broadband polarimetry can provide crucial information on magnetized astronomical objects, such as quasars, galaxies, or galaxy clusters. However, the limited wavelength coverage of the instruments requires that we solve an ill-posed inverse problem when we want to obtain the Faraday dispersion function (FDF), a tomographic distribution of the magnetoionic media along the line of sight. This paper explores the use of wavelet transforms and the sparsity of the transformed FDFs in the form of wavelet shrinkage (WS) for finding better solutions to the inverse problem. We recently proposed the Constraining and Restoring iterative Algorithm for Faraday Tomography (CRAFT; Cooray et al. 2021), a new flexible algorithm that showed significant improvements over the popular methods such as Rotation Measure Synthesis. In this work, we introduce CRAFT+WS, a new version of CRAFT incorporating the ideas of wavelets and sparsity. CRAFT+WS exhibit significant improvements over the original CRAFT when tested for a complex FDF of realistic Galactic model. Reconstructions of FDFs demonstrate super-resolution in Faraday depth, uncovering previously unseen Faraday complexities in observations. The proposed approach will be necessary for effective cosmic magnetism studies using the Square Kilometre Array and its precursors.

We present wind models of ten young Solar-type stars in the Hercules-Lyra association and the Coma Berenices cluster aged around 0.26 Gyr and 0.58 Gyr respectively. Combined with five previously modelled stars in the Hyades cluster, aged 0.63 Gyr, we obtain a large atlas of fifteen observationally based wind models. We find varied geometries, multi-armed structures in the equatorial plane, and a greater spread in quantities such as the angular momentum loss. In our models we infer variation of a factor of ~6 in wind angular momentum loss $\dot J$ and a factor of ~2 in wind mass loss $\dot M$ based on magnetic field geometry differences when adjusting for the unsigned surface magnetic flux. We observe a large variation factor of ~4 in wind pressure for an Earth-like planet; we attribute this to variations in the 'magnetic inclination' of the magnetic dipole axis with respect to the stellar axis of rotation. Within our models, we observe a tight correlation between unsigned open magnetic flux and angular momentum loss. To account for possible underreporting of the observed magnetic field strength we investigate a second series of wind models where the magnetic field has been scaled by a factor of 5. This gives $\dot M \propto B^{0.4}$ and $\dot J \propto B^{1.0}$ as a result of pure magnetic scaling.

We conduct three-dimensional hydrodynamical simulations of common envelope jets supernova (CEJSN) events where we assume that a neutron star (NS) launches jets as it orbits inside the outer zones of a red supergiant (RSG) envelope, and find the negative jet feedback coefficient to be 0.1-0.2. This coefficient is the factor by which the jets reduce the mass accretion rate onto the NS as they remove mass from the envelope and inflate bubbles (cocoons). Our results suggest that in most CEJSN events the NS-RSG binary system experiences the grazing envelope evolution (GEE) before it enters a full common envelope evolution (CEE). We also find that the jets induce up and down flows in the RSG envelope. These flows together with the strong convection of RSG stars might imply that energy transport by convection in CEJSNe is very important. Because of limited numerical resources we do not include in the simulations the gravity of the NS, nor the accretion process, nor the jets launching process, and nor the gravity of the deformed envelope. Future numerical simulations of CEE with a NS/BH companion should include the accretion process onto the NS and vary the jets power accordingly, the full gravitational interaction of the NS with the RSG, and energy transport by the strong convection.

High-quality stellar spectra are in great demand now - they are the most important ingredient in the stellar population synthesis to study galaxies and star clusters. Here we describe the procedures to increase the quality of flux calibration of stellar spectra. We use examples of NIR intermediate-resolution Echelle spectra collected with the Folded InfraRed Echellete (R~6500, Magellan Baade) and high-resolution UV-optical spectra observed with UVES (R~80000, ESO VLT). By using these procedures, we achieved the quality of the global spectrophotometric calibration as good as 1-2%, which fulfills the requirements for the quality of stellar spectra intended to be used in the stellar population synthesis

We examine the redshifts of a comprehensive set of published Type Ia supernovae, and provide a combined, improved catalog with updated redshifts. We improve on the original catalogs by using the most up-to-date heliocentric redshift data available; ensuring all redshifts have uncertainty estimates; using the exact formulae to convert heliocentric redshifts into the Cosmic Microwave Background (CMB) frame; and utilizing an improved peculiar velocity model that calculates local motions in redshift-space and more realistically accounts for the external bulk flow at high-redshifts. In total we reviewed 2821 supernova redshifts; 534 are comprised of repeat-observations of the same supernovae and 1764 pass the cosmology sample quality cuts. We found 5 cases of missing or incorrect heliocentric corrections, 44 incorrect or missing supernova coordinates, 230 missing heliocentric or CMB frame redshifts, and 1200 missing redshift uncertainties. Of the 2287 unique Type Ia supernovae in our sample (1594 of which satisfy cosmology-sample cuts) we updated 990 heliocentric redshifts. The absolute corrections range between $10^{-8} \leq \Delta z \leq 0.038$, and RMS$(\Delta z) \sim 3\times 10^{-3}$. The sign of the correction was essentially random, so the mean and median corrections are small: $4\times 10^{-4}$ and $4\times 10^{-6}$ respectively. We examine the impact of these improvements for $H_0$ and the dark energy equation of state $w$ and find that the cosmological results change by $\Delta H_0 = -0.11$ km s$^{-1}$ Mpc$^{-1}$ and $\Delta w = -0.001$, both significantly smaller than previously reported uncertainties for $H_0$ of 1.4 km s$^{-1}$ Mpc$^{-1}$ and $w$ of 0.04 respectively.

Recently, Active Galactic Nuclei (AGNs) have been proposed as standardizable candles, thanks to an observed non-linear relation between their X-ray and optical-ultraviolet (UV) luminosities, which provides an independent measurement of their distances. In this paper, we use these observables for the first time to estimate the parameters of f(R) gravity models (specifically the Hu-Sawicki and the exponential models) together with the cosmological parameters. The importance of this type of modified gravity theories lies in the fact that they can explain the late time accelerated expansion of the universe without the inclusion of a dark energy component. We have also included other observable data to the analyses such as estimates of the Hubble parameter H(z) from Cosmic Chronometers, the Pantheon Type Ia supernovae compilation, and Baryon Acoustic Oscillations measurements. Our results show that the allowed space parameter is restricted when both AGN and BAO data are added to CC and SnIa data, being the BAO data set the most restrictive one. We can also conclude that even though our results are consistent with the ones from the LCDM model, small deviations from General Relativity, than can be successfully described by the f(R) models studied in this paper, are also allowed by the considered data sets.

The ARIANNA detector is designed to detect neutrinos with energies above $10^{17}$eV. Due to the similarities in generated radio signals, cosmic rays are often used as test beams for neutrino detectors. Some ARIANNA detector stations are equipped with antennas capable of detecting air showers. Since the radio emission properties of air showers are well understood, and the polarization of the radio signal can be predicted from the arrival direction, cosmic rays can be used as a proxy to assess the reconstruction capabilities of the ARIANNA neutrino detector. We report on dedicated efforts of reconstructing the polarization of cosmic-ray radio pulses. After correcting for difference in hardware, the two stations used in this study showed similar performance in terms of event rate and agreed with simulation. Subselecting high quality cosmic rays, the polarizations of these cosmic rays were reconstructed with a resolution of $2.5^{\circ}$ (68% containment), which agrees with the expected value obtained from simulation. A large fraction of this resolution originates from uncertainties in the predicted polarization because of the contribution of the subdominant Askaryan effect in addition to the dominant geomagnetic emission. Subselecting events with a zenith angle greater than $70^{\circ}$ removes most influence of the Askaryan emission, and, with limited statistics, we found the polarization uncertainty is reduced to $1.3^{\circ}$ (68% containment).

We study the evolution of light spectator fields in an asymmetric polynomial potential. During inflation, stochastic fluctuations displace the spectator field from the global minimum of its potential, populating the false vacuum state and thereby allowing for the formation of false vacuum bubbles. By using a lattice simulation, we show that these bubbles begin to contract once they re-enter the horizon and, if sufficiently large, collapse into black holes. This process generally results in the formation of primordial black holes, which, due to the specific shape of their mass function, are constrained to yield at most 1% of the total dark matter abundance. However, the resulting population can source gravitational wave signals observable at the LIGO-Virgo experiments, provide seeds for supermassive black holes or cause a transient matter-dominated phase in the early Universe.

We present version 1.0 of the cosmological simulation code $\scriptstyle{\rm CO}N{\rm CEPT}$, designed for simulations of large-scale structure formation. $\scriptstyle{\rm CO}N{\rm CEPT}\, 1.0$ contains a P$^3$M gravity solver, with the short-range part implemented using a novel (sub)tiling strategy, coupled with individual and adaptive particle time-stepping. In addition to this, $\scriptstyle{\rm CO}N{\rm CEPT}$ contains a (linear or non-linear) fluid solver to treat non-baryonic components which are not easily treatable using the particle implementation. This allows e.g. for the inclusion of non-linear massive neutrinos (which may be relativistic) and for simulations that are consistent with general relativistic perturbation theory. Decaying dark matter scenarios are also fully supported. A primary objective of $\scriptstyle{\rm CO}N{\rm CEPT}$ is ease of use. To this end, it has built-in initial condition generation and can produce output in the form of snapshots, power spectra and direct visualisations. It also comes with a robust installation script and thorough documentation. $\scriptstyle{\rm CO}N{\rm CEPT}$ is the first massively parallel cosmological simulation code written in Python. Despite of this, excellent performance is obtained, even comparing favourably to other codes such as $\scriptstyle{\rm GADGET}$ at similar precision, in the case of low to moderate clustering. We find extraordinary good agreement between results produced using $\scriptstyle{\rm CO}N{\rm CEPT}\, 1.0$ and $\scriptstyle{\rm GADGET}$, at all scales and times. The $\scriptstyle{\rm CO}N{\rm CEPT}$ code itself along with documentation is openly released at https://github.com/jmd-dk/concept .

Perturbative and non-perturbative terms of the cross sections of ultraperipheral production of lepton pairs in ion collisions are taken into account. It is shown that production of low-mass $e^+e^-$ pairs is strongly enhanced (compared to perturbative estimates) due to the non-perturbative Sommerfeld-Gamow-Sakharov (SGS) factor. Coulomb attraction of the non-relativistic components of those pairs leads to the finite value of their mass distribution at lowest relative velocities. Their annihilation can result in the increased intensity of 511 keV photons. It can be recorded at the NICA collider and is especially crucial in astrophysical implications regarding the 511 keV line emitted from the Galactic center. The analogous effect can be observed in lepton pairs production at LHC. Energy spectra of lepton pairs created in ultraperipheral nuclear collisions and their transverse momenta are calculated.

Sterile neutrinos with keV-scale masses are popular candidates for warm dark matter. In the most straightforward case they are produced via oscillations with active neutrinos. We introduce effective self-interactions of active neutrinos and investigate the effect on the parameter space of sterile neutrino mass and mixing. Our focus is on mixing with electron neutrinos, which is subject to constraints from several upcoming or running experiments like TRISTAN, ECHo, BeEST and HUNTER. Depending on the size of the self-interaction, the parameter space moves closer to, or further away from, the one testable by those future experiments. In particular, we show that phase 3 of the HUNTER experiment would test a larger amount of parameter space in the presence of self-interactions than without them. We also investigate the effect of the self-interactions on the free-streaming length of the sterile neutrino dark matter, which is important for structure formation observables.

We propose that a first-order phase transition of a sub-eV scalar field in the dark sector, triggered by the decreasing temperature as the universe expands, can simultaneously relieve the Hubble tension and explain neutrino masses. Here, the supercooled vacuum of the scalar field gives rise to a sizable fraction of an early dark energy component that boosts the expansion before recombination and subsequently decays. The neutrino masses are generated through the inverse seesaw mechanism by making a set of sterile Majorana fermions massive when the scalar field picks up its vacuum expectation value. We embed this low-energy theory in a larger gauge group that is partially broken above the TeV scale. This novel theory, which could even be motivated independently of the Hubble tension, completes the high-energy corner of the inverse seesaw mechanism and explains the mass of a dark matter candidate that can be produced through gravitational interactions at high energies. An approximate global lepton symmetry that is spontaneously broken during the low-energy phase transition protects the neutrino masses against loop corrections.

We consider gravitational wave (GW) production during preheating in hybrid inflation models where an axion-like waterfall field couples to Abelian gauge fields. Based on a linear analysis, we find that the GW signal from such models can be within the reach of a variety of foreseeable GW experiments such as LISA, AEDGE, ET and CE, and is close to that of LIGO A+, both in terms of frequency range and signal strength. Furthermore, the resultant GW signal is helically polarized and thus may distinguish itself from other sources of stochastic GW background. Finally, such models can produce primordial black holes that can compose dark matter and lead to merger events detectable by GW detectors.

New early dark energy (NEDE) makes the cosmic microwave background consistent with a higher value of the Hubble constant inferred from supernovae observations. It is an improvement over the old early dark energy model (EDE) because it explains naturally the decay of the extra energy component in terms of a vacuum first-order phase transition that is triggered by a subdominant scalar field at zero temperature. With hot NEDE, we introduce a new mechanism to trigger the phase transition. It relies on thermal corrections that subside as a subdominant radiation fluid in a dark gauge sector cools. We explore the phenomenology of hot NEDE and identify the strong supercooled regime as the scenario favored by phenomenology. In a second step, we propose different microscopic embeddings of hot NEDE. This includes the (non-)Abelian dark matter model, which has the potential to also resolve the LSS tension through interactions with the dark radiation fluid. We also address the coincidence problem generically present in EDE models by relating NEDE to the mass generation of neutrinos via the inverse seesaw mechanism. We finally propose a more complete dark sector model, which embeds the NEDE field in a larger symmetry group and discuss the possibility that the hot NEDE field is central for spontaneously breaking lepton number symmetry.

We consider single graviton loop corrections to the effective field equation of a massless, minimally coupled scalar on de Sitter background in the simplest gauge. We find a large temporal logarithm in the approach to freeze-in at late times, but no correction to the feeze-in amplitude. We also find a large spatial logarithm (at large distances) in the scalar potential generated by a point source, which can be explained using the renormalization group with one of the higher derivative counterterms regarded as a curvature-dependent field strength renormalization. We discuss how these results set the stage for a project to purge gauge dependence by including quantum gravitational corrections to the source which disturbs the effective field and to the observer who measures it.

Heavy atomic nuclei have an excess of neutrons over protons. This leads to the formation of a neutron skin whose thickness, $R_\mathrm{skin}$, is sensitive to details of the nuclear force -- linking atomic nuclei to properties of neutron stars, thereby relating objects that differ in size by 18 orders of magnitude [1, 2]. ${}^{208}$Pb is of particular interest here because it exhibits a simple structure and is accessible to experiment. However, computing such a heavy nucleus has been out of reach for ab initio theory. By combining advances in quantum many-body methods, statistical tools, and emulator technology, we make quantitative predictions for the properties of ${}^{208}$Pb starting from nuclear forces that are consistent with symmetries of low-energy quantum chromodynamics. We explore $10^9$ different nuclear-force parameterisations via history matching, confront them with data in select light nuclei, and arrive at an importance-weighted ensemble of interactions. We accurately reproduce bulk properties of ${}^{208}$Pb and find $R_\mathrm{skin}({}^{208}\mathrm{Pb}) = 0.14-0.20$ fm which is smaller than a recent extraction from parity-violating electron scattering [3] but in agreement with other experimental probes. The allowable range of $R_\mathrm{skin}({}^{208}\mathrm{Pb})$ is significantly constrained by nucleon-nucleon scattering data, ruling out very thick skins. This work demonstrates that nuclear forces constrained to light systems extrapolate reliably to even the heaviest nuclei, and that we can make quantitative predictions across the nuclear landscape.

Magnetic helicity, and more broadly magnetic field line topology, imposes constraints on the plasma dynamics. Helically interlocked magnetic rings are harder to bring into a topologically non-trivial state than two rings that are not linked. This particular restriction has the consequence that helical plasmas exhibit increased stability in laboratory devices, in the Sun and in the intergalactic medium. Here we discuss how a magnetic field is stabilizing the plasma and preventing it from disruption by the presence of magnetic helicity. We present observational results, numerical experiments and analytical results that illustrate how helical magnetic fields strongly contribute to the long-term stability of some plasmas. We discuss several cases, such as that of solar corona, tokamaks, the galactic and extragalactic medium, with a special emphasis on extragalactic bubbles.

We present an analytic solution for accretion of a gaseous medium with adiabatic equation of state onto a charged dilaton black hole which moves at a constant velocity. We determine the four-velocity of accreted flow and find that it possesses axial symmetry. We obtain the accretion rate which depends on the mass, the magnetic charge, and the dilation of black hole, meaning that these parameters take important roles in the process of accretion. The results may help us to get deeper understanding of the behavior of accreted flow near the event horizon of black hole.

We conduct a linear analysis of axisymmetric magnetorotational instability (MRI) in a magnetized cylindrical Taylor-Couette (TC) flow for its standard version (SMRI) with a purely axial background magnetic field and two further types -- helically modified SMRI (H-SMRI) and helical MRI (HMRI) -- in the presence of combined axial and azimuthal magnetic fields. This study is intended as preparatory for upcoming large-scale liquid sodium MRI experiments planned within the DRESDYN project at Helmholtz-Zentrum Dresden-Rossendorf, so we explore these instability types for typical values of the main parameters: the magnetic Reynolds number, the Lundquist number and the ratio of the angular velocities of the cylinders, which are attainable in these experiments. In contrast to previous attempts of detecting MRI in the lab, our results demonstrate that SMRI and its helically modified version can in principle be detectable in the DRESDYN-TC device for the range of the above parameters, including the astrophysically most important Keplerian rotation, despite the extremely small magnetic Prandtl number of liquid sodium. Since in the experiments we plan to approach (H-)SMRI from the previously studied HMRI regime, we characterize the continuous and monotonous transition between the both regimes. We show that H-SMRI, like HMRI, represents an overstability (travelling wave) with non-zero frequency linearly increasing with azimuthal field. Because of its relevance to finite size flow systems in experiments, we also analyze the absolute form of H-SMRI and compare its growth rate and onset criterion with the convective one.

Inspired by Gounaris-Sakurai and Lee-Zumino, we postulate that the weak vector and axial vector currents are dominated by $J^{PC} = 1^{--}$ and $1^{++}$ resonances respectively in the appropriate channels of $\nu + \bar \nu$ annihilation into quark-antiquark pairs when an ultrahigh-energy incoming $\nu \ (\bar \nu)$ strikes a relic $\bar \nu \ (\nu)$. Despite this and some other ideas, it appears the detection of relic neutrinos with just the Standard Model interactions seems extremely difficult at existing or future neutrino telescopes. Thus any positive signal would be due to some non-standard interactions of neutrinos.

We explore the impact of non-perturbative effects, namely Sommerfeld enhancement and bound state formation, on the cosmological production of non-thermal dark matter. For this purpose, we focus on a class of simplified models with t-channel mediators. These naturally combine the requirements for large corrections in the early Universe, i.e. beyond the Standard Model states with long range interactions, with a sizable new physics production cross section at the LHC. We find that the dark matter yield of the superWIMP mechanism is suppressed considerably due to the non-perturbative effects under consideration. This leads to a significant shift in the cosmologically preferred parameter space of non-thermal dark matter in these models. We also revisit the implications of LHC bounds on long-lived particles associated with non-thermal dark matter and find that testing this scenario at the LHC is a bigger challenge than previously anticipated.

Bound state formation can have a large impact on the dynamics of dark matter freeze-out in the early Universe, in particular for colored coannihilators. We present a general formalism to include an arbitrary number of excited bound states in terms of an effective annihilation cross section, taking bound state formation, decay as well as transitions into account, and derive analytic approximations in the limiting cases of no or efficient transitions. Furthermore, we provide explicit expressions for radiative bound state formation rates for states with arbitrary principal and angular quantum numbers $n,\ell$ for a mediator in the fundamental representation of $SU(3)_c$, as well as electromagnetic transition rates among them in the Coulomb approximation. We then assess the impact of bound states within a model with Majorana dark matter and a colored scalar $t$-channel mediator. We consider the regime of coannihilation as well as conversion-driven freeze-out (or coscattering), where the relic abundance is set by the freeze-out of conversion processes. We find that the region in parameter space where the latter occurs is considerably enhanced into the multi-TeV regime. For conversion-driven freeze-out dark matter is very weakly coupled, evading direct and indirect detection constraints, but leading to prominent signatures of long-lived particles that provide great prospects to be probed by dedicated searches at the upcoming LHC runs.