Papers for Friday, Jul 02 2021

Discoveries of extended rotation curves have suggested the presence of dark matter in spiral galaxy haloes. It has led to many studies that estimated the galaxy total mass, mostly by using the Navarro Frenk and White (NFW) density profile.We aim at verifying how the choice of the dark-matter profile may affect the predicted values of extrapolated total masses. We have considered the recent Milky Way rotation curve, firstly because of its unprecedented accuracy, and secondly because the Galactic disk is amongst the least affected by past major mergers having fully reshaped the initial disk. We find that the use of NFW profile (or its generalized form, gNFW) for calculating the dark-matter contribution to the Milky Way rotation curve generates apparently inconsistent results, e.g., an increase of the baryonic mass leads to increase of the dark matter mass. Furthermore we find that NFW and gNFW profile narrow the total mass range, leading to a possible methodological bias particularly against small MW masses. By using the Einasto profile that is more appropriate to represent cold dark matter haloes, we finally find that the Milky Way slightly decreasing rotation curve favors total mass that can be as small as 2.6 $\times 10^{11}$ $M_{\odot}$, disregarding any other dynamical tracers further out in the Milky Way. It is inconsistent with values larger than 18 $\times 10^{11}$ $M_{\odot}$ for any kind of CDM dark-matter halo profiles, under the assumption that stars and gas do not influence the predicted dark matter distribution in the Milky Way. This methodological paper encourages the use of the Einasto profile for characterizing rotation curves with the aim of evaluating their total masses.

Simulations of exoplanet albedo profiles are key to planning and interpreting future direct imaging observations. In this paper we demonstrate the use of the Planetary Spectrum Generator to produce simulations of reflected light exoplanet spectra. We use PSG to examine multiple issues relevant to all models of directly imaged exoplanet spectra and to produce sample spectra of the bright, nearby exoplanet $\upsilon$ Andromedae d, a potential direct imaging target for next-generation facilities. We introduce a new, fast, and accurate subsampling technique that enables calculations of disk-integrated spectra one order of magnitude faster than Chebyshev-Gauss sampling for moderate- to high-resolution sampling. Using this method and a first-principles-derived atmosphere for $\upsilon$ And d, we simulate phase-dependent spectra for a variety of different potential atmospheric configurations. The simulated spectra for $\upsilon$ And d include versions with different haze and cloud properties. Based on our combined analysis of this planet's orbital parameters, phase-and illumination-appropriate model spectra, and realistic instrument noise parameters, we find that $\upsilon$ And d is a potentially favorable direct imaging and spectroscopy target for the Coronagraph Instrument (CGI) on the Nancy Grace Roman Space Telescope. When a noise model corresponding to the Roman CGI SPC spectroscopy mode is included, PSG predicts the time required to reach a signal-to-noise ratio of 10 of the simulated spectra in both the central wavelength bin of the Roman CGI SPC spectroscopy mode (R=50 spectrum) and of the Band 1 HLC imaging mode is approximately 400 and less than 40 hr, respectively. We also discuss potential pathways to extricating information about the planet and its atmosphere with future observations and find that Roman observations may be able to bound the interior temperature of the planet.

Stellar streams from globular clusters (GCs) offer constraints on the nature of dark matter and have been used to explore the dark matter halo structure and substructure of our Galaxy. Detection of GC streams in other galaxies would broaden this endeavor to a cosmological context, yet no such streams have been detected to date. To enable such exploration, we develop the Hough Stream Spotter (HSS), and apply it to the Pan-Andromeda Archaeological Survey (PAndAS) photometric data of resolved stars in M31's stellar halo. We first demonstrate that our code can re-discover known dwarf streams in M31. We then use the HSS to blindly identify 27 linear GC stream-like structures in the PAndAS data. For each HSS GC stream candidate, we investigate the morphologies of the streams and the colors and magnitudes of all stars in the candidate streams. We find that the five most significant detections show a stronger signal along the red giant branch in color-magnitude diagrams (CMDs) than spurious non-stream detections. Lastly, we demonstrate that the HSS will easily detect globular cluster streams in future Nancy Grace Roman Space Telescope data of nearby galaxies. This has the potential to open up a new discovery space for GC stream studies, GC stream gap searches, and for GC stream-based constraints on the nature of dark matter.

Imaging the cosmic 21 cm signal will map out the first billion years of our Universe. The resulting 3D lightcone (LC) will encode the properties of the unseen first galaxies and physical cosmology. Unfortunately, there is no obvious summary statistic to use when interpreting this non-Gaussian image, and the commonly-used power spectrum may waste valuable information. Here we build on previous work using Convolutional Neural Networks (CNNs) to infer astrophysical parameters directly from 21 cm LC images. Guided by the properties of LCs, we combine recurrent layers characterizing evolution along the redshift axis with 2D convolutional layers characterizing local correlations in the sky-plane. Such Recursive Neural Networks (RNNs) are known for efficiently learning temporal correlations in sequential data. Using a large database of simulated cosmic 21 cm LCs, we confirm that RNNs outperform previously-used CNNs in recovering UV and X-ray galaxy properties, reducing the mean squared parameter estimation error by factors of $\sim 2 - 8$. We also corrupt the cosmic signal by adding noise expected from a 1000 h integration with the Square Kilometre Array, as well as excising a foreground-contaminated ''horizon wedge''. Parameter prediction errors increase when the NNs are trained on these contaminated LC images, though recovery is still good even in the most pessimistic case (with $R^2 \ge 0.5 - 0.95$). However, we find no notable differences in performance between network architectures on the contaminated images. We argue this is due to the size of our dataset, highlighting the need for larger datasets and/or better data augmentation in order to maximize the potential of NNs in 21 cm parameter estimation.

The extreme extinction ($A_V\sim30$) and its variation on arc-second scales towards the Galactic centre hamper the study of its stars. Their analysis is restricted to the near infrared (NIR) regime, where the extinction curve can be approximated by a broken power law for the $JHK_s$ bands. Therefore, correcting for extinction is fundamental to analyse the structure and stellar population of the central regions of our Galaxy. We aim at: (1) Discussing different strategies to de-redden the photometry and checking the usefulness of extinction maps to deal with variable stars. (2) Building extinction maps for the NIR bands $JHK_s$. (3) Creating a de-reddened GALACTICNUCLEUS (GNS) survey, identifying foreground stars. (4) Performing a preliminary analysis of the de-reddened $K_s$ luminosity functions (KLFs). We use photometry from the GNS survey to create extinction maps for the whole catalogue. We take red clump (RC) and red giant stars of similar brightness as a reference to build the maps, and de-redden the GNS photometry. We discuss the limitations of the process and analyse non-linear effects. We create high resolution ($\sim3''$) extinction maps with low statistical and systematics uncertainties ($\lesssim5$\,\%), and compute average extinctions for each of the regions covered by the GNS. We check that our maps effectively correct the differential extinction reducing the spread of the RC features by a factor of $\sim2$. We assess the validity of the broken power law approach computing two equivalent extinction maps $A_H$ using either $JH$ and $HK_s$ photometry for the same reference stars, and obtain compatible average extinctions within the uncertainties. Finally, we analyse de-reddened KLFs for different line-of-sights and find that the regions belonging to the NSD contain a homogeneous stellar population that is significantly different from the one in the innermost bulge regions.

We present a new, model-independent measurement of the clustering amplitude of galaxies and the growth of cosmic large-scale structures from the Baryon Oscillation Spectroscopic Survey (BOSS) 12th data release (DR12). This is achieved by generalising harmonic-space power spectra for galaxy clustering to measure separately the magnitudes of the density and of the redshift-space distortion terms, which are respectively related to the clustering amplitude, $b\sigma_8(z)$, and the growth, $f\sigma_8(z)$. We adopt a tomographic approach with 15 redshift bins in the range $z\in[0.15,0.67]$. We restrict our analysis to strictly linear scales, implementing a redshift-dependent maximum multipole for each of the tomographic bins. Thus, we obtain 30 data points in total, 15 for each of the quantities $b\sigma_8(z)$ and $f\sigma_8(z)$. The measurements do not appear to suffer from any apparent systematic effect and show excellent agreement with the theoretical prediction from a concordance cosmology as from the Planck satellite. Our results also agree with previous analyses by the BOSS collaboration. Although each single datum has, in general, a larger error bar than that obtained in configuration- or Fourier-space analyses, our study provides the community with a larger number of tomographic data points that allow for a complementary tracking in redshift of the evolution of fundamental cosmological quantities.

We report the detection of a single transit-like signal in the Kepler data of the slightly evolved F star KIC4918810. The transit duration is ~45 hours, and while the orbital period ($P\sim10$ years) is not well constrained, it is one of the longest among companions known to transit. We calculate the size of the transiting object to be $R_P = 0.910$ $R_J$. Objects of this size vary by orders of magnitude in their densities, encompassing masses between that of Saturn ($0.3$ $M_J$) and stars above the hydrogen-burning limit (~80 $M_J$). Radial-velocity observations reveal that the companion is unlikely to be a star. The mass posterior is bimodal, indicating a mass of either ~0.24 $M_J$ or ~26 $M_J$. Continued spectroscopic monitoring should either constrain the mass to be planetary or detect the orbital motion, the latter of which would yield a benchmark long-period brown dwarf with a measured mass, radius, and age.

Stars originate from the dense interstellar medium, which exhibits filamentary structure to scales of $\sim 1$ kpc in galaxies like our Milky Way. We explore quantitatively how much resulting large-scale correlation there is among different stellar clusters and associations in $\textit{orbit phase space}$, characterized here by actions and angles. As a starting point, we identified 55 prominent stellar overdensities in the 6D space of orbit (actions) and orbital phase (angles), among the $\sim$ 2.8 million stars with radial velocities from Gaia EDR3 and $d \leq 800$ pc. We then explored the orbital $\textit{phase}$ distribution of all sample stars in the same $\textit{orbit}$ patch as any one of these 55 overdensities. We find that very commonly numerous other distinct orbital phase overdensities exist along these same orbits, like pearls on a string. These `pearls' range from known stellar clusters to loose, unrecognized associations. Among orbit patches defined by one initial orbit-phase overdensity 50% contain at least 8 additional orbital-phase pearls of 10 cataloged members; 20% of them contain 20 additional pearls. This is in contrast to matching orbit patches sampled from a smooth mock catalog, or random nearby orbit patches, where there are only 2 (or 5, respectively) comparable pearls. Our findings quantify for the first time how common it is for star clusters and associations to form at distinct orbital phases of nearly the same orbit. This may eventually offer a new way to probe the 6D orbit structure of the filamentary interstellar medium.

We leverage state-of-the-art machine learning methods and a decade's worth of archival data from the Canada-France-Hawaii Telescope (CFHT) to predict observatory image quality (IQ) from environmental conditions and observatory operating parameters. Specifically, we develop accurate and interpretable models of the complex dependence between data features and observed IQ for CFHT's wide field camera, MegaCam. Our contributions are several-fold. First, we collect, collate and reprocess several disparate data sets gathered by CFHT scientists. Second, we predict probability distribution functions (PDFs) of IQ, and achieve a mean absolute error of $\sim0.07''$ for the predicted medians. Third, we explore data-driven actuation of the 12 dome ``vents'', installed in 2013-14 to accelerate the flushing of hot air from the dome. We leverage epistemic and aleatoric uncertainties in conjunction with probabilistic generative modeling to identify candidate vent adjustments that are in-distribution (ID) and, for the optimal configuration for each ID sample, we predict the reduction in required observing time to achieve a fixed SNR. On average, the reduction is $\sim15\%$. Finally, we rank sensor data features by Shapley values to identify the most predictive variables for each observation. Our long-term goal is to construct reliable and real-time models that can forecast optimal observatory operating parameters for optimization of IQ. Such forecasts can then be fed into scheduling protocols and predictive maintenance routines. We anticipate that such approaches will become standard in automating observatory operations and maintenance by the time CFHT's successor, the Maunakea Spectroscopic Explorer (MSE), is installed in the next decade.

Based on new observations during 2015-2020 and published data, the unusual eruptive variables PV Cep and V350 Cep are examined. It is shown that PV Cep underwent a regular outburst followed by a drop in brightness that lasted overall from 2011 to 2019 and is still in a deep minimum. The outburst was accompanied by substantial changes in the intensity and profiles of a number of lines, including Ha, [SII], and [OI]. The forbidden lines generally have negative radial velocities and can be divided into four components, with variable velocities and relative intensities. V350 Cep essentially is at a maximum brightness level over the entire time and its spectrum is practically unaltered. The available data suggest that the pronounced P Cyg profile of the Ha line in the spectrum of V350 Cep appeared several years after the luminosity rise, in 1986. The luminosities of the stars in the current state are estimated to be 20 L(sun) and 3.3 L(sun), respectively. It is concluded that both stars may represent a so-called intermediate objects between the FUor and EXor classes.

Our project MOMO (Multiwavelength observations and modelling of OJ 287) consists of dedicated, dense, long-term flux and spectroscopic monitoring and deep follow-up observations of the blazar OJ 287 at >13 frequencies from the radio to the X-ray band since late 2015. In particular, we are using Swift to obtain optical-UV-X-ray spectral energy distributions (SEDs) and the Effelsberg telescope to obtain radio measurements between 2 and 40 GHz. MOMO is the densest long-term monitoring of OJ 287 involving X-rays and broad-band SEDs. The theoretical part of the project aims at understanding jet and accretion physics of the blazar central engine in general and the supermassive binary black hole scenario in particular. Results are presented in a sequence of publications and so far included: detection and detailed analysis of the bright 2016/17 and 2020 outbursts and the long-term light curve; Swift, XMM and NuSTAR spectroscopy of the 2020 outburst around maximum; and interpretation of selected events in the context of the binary black hole scenario of OJ 287 (papers I-IV). Here, we provide a description of the project MOMO, a summary of previous results, the latest results, and we discuss future prospects.

We consider the possibility that the Milky Way's dark matter halo possesses a non vanishing equation of state. Consequently, we evaluate the contribution due to the speed of sound, assuming that the dark matter content of the galaxy behaves like a fluid with pressure. In particular, we model the dark matter distribution via an exponential sphere profile in the galactic core, and inner parts of the galaxy whereas we compare the exponential sphere with three widely-used profiles for the halo, i.e. the Einasto, Burkert and Isothermal profile. For the galactic core we also compare the effects due to a dark matter distribution without black hole with the case of a supermassive black hole in vacuum and show that present observations are unable to distinguish them. Finally we investigate the expected experimental signature provided by gravitational lensing due to the presence of dark matter in the core.

We present the analysis of the microlensing event KMT-2018-BLG-1743. The light curve of the event, with a peak magnification $A_{\rm peak}\sim 800$, exhibits two anomaly features, one around the peak and the other on the falling side of the light curve. An interpretation with a binary lens and a single source (2L1S) cannot describe the anomalies. By conducting additional modeling that includes an extra lens (3L1S) or an extra source (2L2S) relative to a 2L1S interpretation, we find that 2L2S interpretations with a planetary lens system and a binary source best explain the observed light curve with $\Delta\chi^2\sim 188$ and $\sim 91$ over the 2L1S and 3L1S solutions, respectively. Assuming that these $\Delta\chi^2$ values are adequate for distinguishing the models, the event is the fourth 2L2S event and the second 2L2S planetary event. The 2L2S interpretations are subject to a degeneracy, resulting in two solutions with $s>1.0$ (wide solution) and $s<1.0$ (close solution). The masses of the lens components and the distance to the lens are $(M_{\rm host}/M_\odot, M_{\rm planet}/M_{\rm J}, D_{\rm L}/{\rm kpc}) \sim (0.19^{+0.27}_{-0.111}, 0.25^{+0.34}_{-0.14}, 6.48^{+0.94}_{-1.03})$ and $\sim (0.42^{+0.34}_{-0.25}, 1.61^{+1.30}_{-0.97}, 6.04^{+0.93}_{-1.27})$ according to the wide and close solutions, respectively. The source is a binary composed of an early G dwarf and a mid M dwarf. The values of the relative lens-source proper motion expected from the two degenerate solutions, $\mu_{\rm wide}\sim 2.3 $mas yr$^{-1}$ and $\mu_{\rm close} \sim 4.1 $mas yr$^{-1}$, are substantially different, and thus the degeneracy can be broken by resolving the lens and source from future high-resolution imaging observations.

We study differential rotation in late-stage shell convection in a 3D hydrodynamic simulation of a rapidly rotating $16M_\odot$ helium star with a particular focus on the convective oxygen shell. We find that the oxygen shell develops a quasi-stationary pattern of differential rotation that is neither described by uniform angular velocity as assumed in current stellar evolution models of supernova progenitors, nor by uniform specific angular momentum. Instead, the oxygen shell develops a positive angular velocity gradient with faster rotation at the equator than at the pole by tens of percent. We show that the angular momentum transport inside the convection zone is not adequately captured by a diffusive mixing-length flux proportional to the angular velocity or angular momentum gradient. Zonal flow averages reveal stable large-scale meridional flow and an entropy deficit near the equator that mirrors the patterns in the angular velocity. The structure of the flow is reminiscent of simulations of stellar surface convection zones and the differential rotation of the Sun, suggesting that similar effects are involved; future simulations will need to address in more detail how the interplay of buoyancy, inertial forces, and turbulent stresses shapes differential rotation during late-stage convection in massive stars. Our findings may have implications for neutron star birth spin periods and supernova explosion scenarios that involve rapid core rotation. If convective regions develop positive angular velocity gradients, angular momentum could be shuffled out of the core region more efficiently, potentially making the formation of millisecond magnetars more difficult.

MeV gamma-ray is a unique window for the direct measurement of line emissions from radioisotopes, but there is no significant progress in the observation after COMPTEL/{\it CGRO}. Hence, for observing celestial objects in this band, we are developing an electron-tracking Compton camera (ETCC), which enables us to perform true imaging spectroscopy similar to X-ray or GeV telescopes. Therefore, we can obtain the energy spectrum of the observation target by a simple ON-OFF method using the correctly defined a proper point-spread function. For validating the performance of celestial object observation using an ETCC, the second balloon SMILE-2+, which had an ETCC based on a gaseous electron tracker with a volume of 30$\times$30$\times$30~cm$^3$, was launched at Alice Springs, Australia on April 7, 2018. SMILE-2+ observed the southern sky including the Crab nebula with a live time of 5.1 h at the zenith angle of $\sim$50 degrees and detected gamma-rays from the Crab nebula with a significance of 4.0$\sigma$ at the energy range of 0.15--2.1~MeV. Additionally, an enhancement of gamma-ray events due to the Galactic center region was clearly observed in the light curve. The realized detection sensitivity agrees well with the sensitivity estimated before launching based on the total background of extragalactic diffuse, atmospheric gamma-rays, and a small number of instrumental gamma-rays suppressed to one-third of the total background. We have succeeded to overcome the most difficult and serious problem of huge background for the stagnation of MeV gamma-ray astronomy for the first time in the world, and thus demonstrate that an ETCC can pioneer a deeper survey than COMPTEL in MeV gamma-ray astronomy.

The effects in the local U-V velocity field due to orbital trapping by bar resonances have been studied computing fifteen resonant families in a non-axisymmetric Galactic potential, considering the bar's angular velocity between 35 and 57.5 ${\rm\,km\,s^{-1}{kpc}^{-1}}$. Only cases in the low, 37.5, 40 ${\rm\,km\,s^{-1}{kpc}^{-1}}$, and high, 55, 57.5 ${\rm\,km\,s^{-1}{kpc}^{-1}}$, velocity ranges give trapping structures that have some similarity with observed features in the velocity distribution. The resulting structures in the local U-V plane form resonant bands appearing at various levels in velocity V. Cases with angular velocity 40 and 55 ${\rm\,km\,s^{-1}{kpc}^{-1}}$ show the greatest similarity with observed branches. Our best approximation to the local velocity field by orbital trapping is obtained with a bar angular velocity of 40 ${\rm\,km\,s^{-1}{kpc}^{-1}}$ and a bar angle of 40${^\circ}$. With this solution, three main observed features can be approximated: i) the Hercules branch at V=$-50$ ${\rm\,km\,s^{-1}}$ produced by the resonance 8/1 outside corotation, and the close features produced by resonances 5/1 and 6/1, ii) the newly detected low-density arch at V $\simeq$ 40 ${\rm\,km\,s^{-1}}$ produced approximately by the resonance 4/3, iii) the inclined structure below the Hercules branch, also observed in the $\textit{Gaia}$ DR2 data, produced by tube orbits around Lagrange point $L_5$ at corotation. Some predicted contributions due to orbital trapping in regions of the U-V plane corresponding to the Galactic halo are given, which could help to further restrict the value of the angular velocity of the Galactic bar. No support by orbital trapping is found for the Arcturus stream at V $\approx$ $-100$ ${\rm\,km\,s^{-1}}$.

The study aims to examine the spectral dynamics of the low-frequency, optically thick gyrosynchrotron microwave emission in solar flares to determine the characteristics of the emitting source. We present the high-resolution spectra of a set of microwave bursts observed by the Expanded Owens Valley Solar Array (EOVSA) during its commissioning phase in the $2.5-18$ GHz frequency range with $1$ second time resolution. Out of the 12 events analyzed in this study, nine bursts exhibit a direct decrease with time in the optically thick spectral index $\alpha_l$, an indicator of source morphology. Particularly, five bursts display "flat" spectrum ($\alpha_l\leq1.0$) compared to that expected for a homogeneous/uniform source ($\alpha_l\approx2.9$). These flat spectra at the low-frequencies (<$10$ GHz) can be defined as the emission from a spatially inhomogeneous source with a large area and/or with multiple emission components. In a subset of six events with partial cross-correlation data, both the events with flat spectra show a source size of $\sim120$ arcsec at $2.6-3$ GHz. Modeling based on inhomogeneity supports the conclusion that multiple discrete sources can only reproduce a flat spectrum. We report that these flat spectra appear predominantly in the decay phase and typically grow flatter over the duration in most of the bursts, which indicates the increasing inhomogeneity and complexity of the emitting volume as the flare progresses. This large volume of flare emission filled with the trapped energetic particles is often invisible in other wavelengths, like hard X-rays, presumably due to the collisionless conditions in these regions of low ambient density and magnetic field strength.

Aims: We investigate the photometric properties of the M7.5 equal-mass binary VHS J1256-1257AB, which, combined with the late-L dwarf VHS J1256-1257 b, forms one of the few young triple systems of ultra-cool dwarfs currently known. Methods: We analyzed two-minute TESS and two-second Spitzer archival data with total durations of about 25 days and 36 hours, respectively. Typical precision in the data is $\pm$1.5 % for TESS and $\pm$0.1 % (in 1 minute) for Spitzer. Results: The optical and infrared light curves periodically exhibit epochs of quasi-sinusoidal modulation followed by epochs of stochastic variability, which resembles the beat pattern created by two waves of similar frequencies that interfere with each other. Our two-wave model for the TESS data shows that the components of VHS J1256-1257AB rotate with periods of $2.0782\pm0.0004$ h and $2.1342\pm0.0003$ h, which is also supported by the Spitzer observations. As a result, the fluxes of the equally bright VHS J1256-1257A and B alternate between states of phase and anti-phase, explaining the observed photometric variability in their combined light. The projected spectroscopic velocity of VHS J1256-1257AB is remarkably similar to those obtained by combining the measured rotation periods and the expected radii, which indicates that the spin axes of VHS J1256-1257A and B are likely inclined at nearly 90 deg, as previously reported for VHS J1256-1257 b.

We carried out an extensive analysis of the stability of the outer solar system, making use of the frequency analysis technique over short-term integrations of nearly a hundred thousand test particles, as well as a statistical analysis of 200, 1 Gyr long numerical simulations, which consider the mutual perturbations of the giant planets and the 34 largest trans-Neptunian objects (we have called all 34 objects ``dwarf planets'', DPs, even if probably only the largest of them are true DPs). From the frequency analysis we produced statistical diffusion maps for a wide region of the $a$-$e$ phase-space plane; we also present the average diffusion time for orbits as a function of perihelion. We later turned our attention to the 34 DPs making an individualized analysis for each of them and producing a first approximation of their future stability. From the 200 distinct realizations of the orbital evolution of the 34 DPs, we classified the sample into three categories including 17 Stable, 11 Unstable, and 6 Resonant objects each; we also found that statistically, 2 objects from the sample will leave the trans-Neptunian region within the next Gyr, most likely being ejected from the solar system, but with a non-negligible probability of going inside the orbit of Neptune, either to collide with a giant planet or even falling to the inner solar system, where our simulations are no longer able to resolve their continuous evolution.

We compare mid-IR and ALMA far-IR images of 12 nearby Seyferts selected from GATOS. The mid-IR unresolved emission contributes more than 60% of the nuclear emission in most galaxies. By contrast, the ALMA 870micron continuum emission is mostly resolved and typically along the torus equatorial direction (Paper I, Garcia-Burillo et al. 2021). The Eddington ratios and nuclear hydrogen column densities NH of half the sample are favorable to launching polar and/or equatorial dusty winds, according to simulations. Six show mid-IR extended emission in the polar direction as traced by the NLR and perpendicular to the ALMA emission. In a few, the nuclear NH might be too high to uplift large quantities of dusty material along the polar direction. Five galaxies have low NH and/or Eddington ratios and thus polar dusty winds are not likely. We generate new CAT3D-WIND disk-wind model images. At low wind-to-disk cloud ratios the far-IR model images have disk- and ring-like morphologies. The X-shape associated with dusty winds is seen better in the far-IR at intermediate-high inclinations for the extended-wind configurations. In most models, the mid-IR emission comes from the inner part of the disk/cone. Extended bi-conical and one-sided polar mid-IR emission is seen in extended-wind configurations and high wind-to-disk cloud ratios. When convolved to our resolution, the model images reproduce qualitative aspects of the observed morphologies. Low-intermediate wind-to-disk ratios are required to account for the large fractions of unresolved mid-IR emission. This work and Paper I provide observational support for the torus+wind scenario. The wind component is more relevant at high Eddington ratios and/or AGN luminosities, and polar dust emission is predicted at NH of up to $10^{24}$cm$^{-2}$. The torus/disk component, on the other hand, prevails at low luminosities and/or Eddington ratios. (Abridged)

We perform two-dimensional particle-in-cell simulations of reconnection in magnetically dominated electron-positron plasmas subject to strong Compton cooling. We vary the magnetization $\sigma\gg1$, defined as the ratio of magnetic tension to plasma inertia, and the strength of cooling losses. Magnetic reconnection under such conditions can operate in magnetically dominated coronae around accreting black holes, which produce hard X-rays through Comptonization of seed soft photons. We find that the particle energy spectrum is dominated by a peak at mildly relativistic energies, which results from bulk motions of cooled plasmoids. The peak has a quasi-Maxwellian shape with an effective temperature of $\sim 100$~keV, which depends only weakly on the flow magnetization and the strength of radiative cooling. The mean bulk energy of the reconnected plasma is roughly independent of $\sigma$, whereas the variance is larger for higher magnetizations. The spectra also display a high-energy tail, which receives $\sim 25$% of the dissipated reconnection power for $\sigma=10$ and $\sim 40$% for $\sigma=40$. We complement our particle-in-cell studies with a Monte-Carlo simulation of the transfer of seed soft photons through the reconnection layer, and find the escaping X-ray spectrum. The simulation demonstrates that Comptonization is dominated by the bulk motions in the chain of Compton-cooled plasmoids and, for $\sigma\sim 10$, yields a spectrum consistent with the typical hard state of accreting black holes.

The high interstellar abundances of polycyclic aromatic hydrocarbons (PAHs) and their size distribution are the result of complex chemical processes implying dust, UV radiation, and the main gaseous components (H, C+, and O). These processes must explain the high abundance of relatively small PAHs in the diffuse interstellar medium (ISM) and imply the continuous formation of some PAHs that are small enough (number of carbon atoms NC <~ 35-50) to be completely dehydrogenated by interstellar UV radiation. The carbon clusters Cn thus formed are constantly exposed to the absorption of ~10-13.6 eV UV photons, allowing isomerization and favoring the formation of the most stable isomers. They might tend to form irregular carbon cages. The frequent accretion of interstellar C+ ions could favor further cage isomerization, as is known in the laboratory for C60, possibly yielding most stable fullerenes, such as C40, C44, and C50. These fullerenes are expected to be very stable in the diffuse ISM because C2 ejection is not possible by single UV photon absorption, but could need rare two-photon absorption. It is possible that at least one of these fullerenes or its cation is as abundant as C60 or C60+ in the diffuse ISM, although this abundance is limited by the lack of observed matching features in observed mid-infrared spectra. B3LYP calculations of the visible spectrum for a number of fullerene isomers with 40 <~ NC <~ 50 show that they generally have a few spectral bands in the visible range, with f-values in the range of a few 10-2. This could make such fullerenes interesting candidates for the carriers of some diffuse interstellar bands.

Our recent work attempted to account for the anisotropy of galactic cosmic rays (GCRs) through the contribution of nearby Geminga SNR source, however the anisotropic phase of CRs is about $\sim R.A.= 5^h$ below 100 TeV, which deviates significantly from the experimental data. Recent observations of CR anisotropy indicate that the phase of the anisotropy of CRs below 100 TeV is coincident with local interstellar magnetic field $(l,b= 210.5^\circ,-57.1^\circ)$ observed by Interstellar Boundary Explorer (IBEX). In this work, we consider simultaneously the contributions of both local interstellar magnetic fields and nearby Geminga SNR source to anisotropy of CRs. We found that the anisotropic phase of CRs points to $\sim R.A.= 3^h$ that accord with local regular magnetic field below 100 TeV, which is due to effect of magnetic field deflection on CRs. We further analyze the influence of the ratio of perpendicular and parallel to the magnetic diffusion coefficient on the anisotropy of CRs. The results illustrate that with the decrease of ratio, the anisotropic phase of CRs changes from the direction of nearby source to local regular magnetic field below 100 TeV, meanwhile from the galactic center (GC) to opposite direction of local regular magnetic field above 100 TeV.

White dwarfs which exhibit transit signatures of planetary debris and accreted planetary material provide exceptional opportunities to probe the material composition and dynamical structure of planetary systems. Although previous theoretical work investigating the role of minor body disruption around white dwarfs has focussed on spherical bodies, Solar System asteroids can be more accurately modelled as triaxial ellipsoids. Here we present an analytical framework to identify the type of disruption (tidal fragmentation, total sublimation or direct impact) experienced by triaxial asteroids approaching white dwarfs on extremely eccentric (e \sim 1) orbits. This framework is then used to identify the outcomes for simplified Main belt analogues of 100 bodies across five different white dwarf temperatures. We also present an empirical relationship between cooling age and effective temperature for both DA and DB white dwarfs to identify the age of the white dwarfs considered here. We find that using a purely spherical shape model can underestimate the physical size and radial distance at which an asteroid is subjected to complete sublimation, and these differences increase with greater elongation of the body. Contrastingly, fragmentation always occurs in the largest semi-axis of a body and so can be modelled by a sphere of that radius. Both fragmentation and sublimation are greatly affected by the body's material composition, and hence by the composition of their progenitor asteroid belts. The white dwarf temperature, and hence cooling age, can affect the expected debris distribution: higher temperatures sublimate large elongated asteroids, and cooler temperatures accommodate more direct impacts.

The peculiar motion of the observer, if not accurately accounted for, is bound to induce a well-defined clustering signal in the distribution of galaxies. This signal is related to the Kaiser rocket effect. Here we examine the amplitude and form of this effect, both analytically and numerically, and discuss possible implications for the analysis and interpretation of forthcoming cosmological surveys. For an idealistic cosmic variance dominated full-sky survey with a Gaussian selection function peaked at $z\sim 1.5$ it is a $> 5\sigma$ effect and it can in principle bias very significantly the inference of cosmological parameters, especially for primordial non-Gaussianity. For forthcoming surveys, with realistic masks and selection functions, the Kaiser rocket is not a significant concern for cosmological parameter inference. However, it is a systematic effect whose origin, nature and imprint on galaxy maps are well known and thus should be subtracted or mitigated. We present several approaches to do so.

Planet formation models begin with proto-embryos and planetesimals already fully formed, missing out a crucial step, the formation of planetesimals/proto-embryos. In this work, we include prescriptions for planetesimal and proto-embryo formation arising from pebbles becoming trapped in short-lived pressure bumps, in thermally evolving viscous discs to examine the sizes and distributions of proto-embryos and planetesimals throughout the disc. We find that planetesimal sizes increase with orbital distance, from ~10 km close to the star to hundreds of kilometres further away. Proto-embryo masses are also found to increase with orbital radius, ranging from $10^{-6} M_{\rm \oplus}$ around the iceline, to $10^{-3} M_{\rm \oplus}$ near the orbit of Pluto. We include prescriptions for pebble and planetesimal accretion to examine the masses that proto-embryos can attain. Close to the star, planetesimal accretion is efficient due to small planetesimals, whilst pebble accretion is efficient where pebble sizes are fragmentation limited, but inefficient when drift dominated due to low accretion rates before the pebble supply diminishes. Exterior to the iceline, planetesimal accretion becomes inefficient due to increasing planetesimal eccentricities, whilst pebble accretion becomes more efficient as the initial proto-embryo masses increase, allowing them to significantly grow before the pebble supply is depleted. Combining both scenarios allows for more massive proto-embryos at larger distances, since the accretion of planetesimals allows pebble accretion to become more efficient, allowing giant planet cores to form at distances upto 10 au. By including more realistic initial proto-embryo and planetesimal sizes, as well as combined accretion scenarios, should allow for a more complete understanding in the beginning to end process of how planets and planetary systems form.

Machine learning techniques have been increasingly used in astronomical applications and have proven to successfully classify objects in image data with high accuracy. The current work uses archival data from the Faint Images of the Radio Sky at Twenty Centimeters (FIRST) to classify radio galaxies into four classes: Fanaroff-Riley Class I (FRI), Fanaroff-Riley Class II (FRII), Bent-Tailed (BENT), and Compact (COMPT). The model presented in this work is based on Convolutional Neural Networks (CNNs). The proposed architecture comprises three parallel blocks of convolutional layers combined and processed for final classification by two feed-forward layers. Our model classified selected classes of radio galaxy sources on an independent testing subset with an average of 96\% for precision, recall, and F1 score. The best selected augmentation techniques were rotations, horizontal or vertical flips, and increase of brightness. Shifts, zoom and decrease of brightness worsened the performance of the model. The current results show that model developed in this work is able to identify different morphological classes of radio galaxies with a high efficiency and performance

We derive a catalog of early-type emission-line stars including 30,048 spectra of 25,886 stars from LAMOST DR7, in which 3,922 have Simbad records. The sample is obtained using K-Nearest Neighbor and Random Forest methods and visually inspected. The spectra are classified into 3 morphological types (10 subtypes) based on H$\alpha$ emission line profiles. Some spectra contaminated by nebula emission lines such as from HII regions are flagged in the catalog. We also provide a specific sub-catalog of 101 stars with the stellar wind by calculating stellar wind or accretion flow velocities based on P Cygni or inverse P Cygni profiles, in which 74\% of them having velocities below 400km/s. More important, with two color-color diagrams (H-K, K- W1) and (H-K, J-H) of a collection of known Herbig Ae/Be stars (HAeBes) and classical Ae/Be stars (CAeBes), we propose an updated criterion to separate HAeBes from CAeBes. By the criterion, we select 201 HAeBes candidates and 5,547 CAeBes candidates from the sample. We confirmed 66 of the 201 HAeBes by using both WISE images and LAMOST spectra and present a specific HAeBe sub-catalog, in which 58 are newly identified. In addition, the WISE colors (W1-W2, W1- W3, and W1-W4) show the distribution consistency between our confirmed HAeBes and that of the known HAeBes. Most of the 66 confirmed HAeBes locate in the lower edge of the main sequence of hot end in the HR diagram, while the distances of about 77\% exceed 1Kpc, which enlarges the number of far distant known HAeBes.

Nuclear rings are excellent laboratories for studying intense star formation. We present results from a study of nuclear star-forming rings in five nearby normal galaxies from the Star Formation in Radio Survey (SFRS) and four local LIRGs from the Great Observatories All-sky LIRG Survey (GOALS) at sub-kpc resolutions using VLA high-frequency radio continuum observations. We find that nuclear ring star formation (NRSF) contributes 49 - 60\% of the total star formation of the LIRGs, compared to 7 - 40\% for the normal galaxies. We characterize a total of 58 individual star-forming regions in these rings, and find that with measured sizes of 10 - 200 pc, NRSF regions in the LIRGs have SFR and $\Sigma_\mathrm{SFR}$ up to 1.7 M$_\odot$yr$^{-1}$ and 402 M$_\odot$yr$^{-1}$kpc$^{-2}$, respectively, which are about 10 times higher than NRSF regions in the normal galaxies with similar sizes, and comparable to lensed high-$z$ star-forming regions. At $\sim 100 - 300$ pc scales, we estimate low contributions ($< 50\%$) of thermal free-free emission to total radio continuum emission at 33 GHz in the NRSF regions in the LIRGs, but large variations possibly exist at smaller physical scales. Finally, using archival sub-kpc resolution CO (J=1-0) data of nuclear rings in the normal galaxies and NGC 7469 (LIRG), we find a large scatter in gas depletion times at similar molecular gas surface densities, which tentatively points to a multi-modal star formation relation on sub-kpc scales.

The chemical and physical evolution of starless and pre-stellar cores are of paramount importance to understanding the process of star formation. The Taurus Molecular Cloud cores TMC 1-C and TMC 1-CP share similar initial conditions and provide an excellent opportunity to understand the evolution of the pre-stellar core phase. We investigated the evolutionary stage of starless cores based on observations towards the prototypical dark cores TMC 1-C and TMC 1-CP, mapping them in the CS $3\rightarrow 2$, C$^{34}$S $3\rightarrow 2$, $^{13}$CS $2\rightarrow 1$, DCN $1\rightarrow 0$, DCN $2\rightarrow 1$, DNC $1\rightarrow 0$, DNC $2\rightarrow 1$, DN$^{13}$C $1\rightarrow 0$, DN$^{13}$C $2\rightarrow 1$, N$_2$H$^+$ $1\rightarrow 0$, and N$_2$D$^+$ $1\rightarrow 0$ transitions. We performed a multi-transitional study of CS and its isotopologs, DCN, and DNC lines to characterize the physical and chemical properties of these cores. We studied their chemistry using the state-of-the-art gas-grain chemical code Nautilus and pseudo time-dependent models to determine their evolutionary stage. Observational diagnostics seem to indicate that TMC 1-C is in a later evolutionary stage than TMC 1-CP, with a chemical age $\sim$1 Myr. TMC 1-C shows signs of being an evolved core at the onset of star formation, while TMC 1-CP appears to be in an earlier evolutionary stage due to a more recent formation or, alternatively, a collapse slowed down by a magnetic support.

In this article a new, multi-functional, high-vacuum astrophysical ice setup, VIZSLA (Versatile Ice Zigzag Sublimation Setup for Laboratory Astrochemistry), is introduced. The instrument allows the investigation of astrophysical processes both in a low-temperature para-H2 matrix and in astrophysical analog ices. In para-H2 matrix the reaction of astrochemical molecules with H atoms and H+ ions can be studied very effectively. For the investigation of astrophysical analog ices the setup is equipped with different irradiation and particle sources: an electron gun, for modeling cosmic rays; an H atom beam source (HABS); a microwave H atom lamp, for generating H Lyman-alpha radiation, and a tunable (213 nm to 2800 nm) laser source. For analysis, an FT-IR (and a UV-Visible) spectrometer and a quadrupole mass analyzer are available. The setup has two cryostats, offering novel features for analysis. Upon the so-called temperature-programmed desorption (TPD) the molecules, desorbing from the first cryostat, can be mixed with Ar and can be deposited onto the substrate of the other cryostat. The well-resolved spectrum of the molecules isolated in an Ar matrix serves a unique opportunity to identify the desorbing products of a processed ice. Some examples are provided to show how the para-H2 matrix experiments and the TPD -- matrix-isolation recondensation experiments can help to understand astrophysically important chemical processes at a low temperature. It is also discussed, how these experiments can complement the studies carried out by similar astrophysical ice setups.

In order to study the ram-pressure interaction between radio galaxies and the intracluster medium, we analyse a sample of 208 highly-bent narrow-angle tail radio sources (NATs) in clusters, detected by the LOFAR Two-metre Sky Survey. For NATs within $7\,R_{500}$ of the cluster centre, we find that their tails are distributed anisotropically with a strong tendency to be bent radially away from the cluster, which suggests that they are predominantly on radially inbound orbits. Within $0.5\,R_{500}$, we also observe an excess of NATs with their jets bent towards the cluster core, indicating that these outbound sources fade away soon after passing pericentre. For the subset of NATs with spectroscopic redshifts, we find the radial bias in the jet angles exists even out to $10\,R_{500}$, far beyond the virial radius. The presence of NATs at such large radii implies that significant deceleration of the accompanying inflowing intergalactic medium must be occurring there to create the ram pressure that bends the jets, and potentially even triggers the radio source.

Two gravitational wave events, namely GW200105 and GW200115, were observed by the Advanced LIGO and Virgo detectors recently. In this work, we show that they can be explained by a scenario of primordial black hole binaries that are formed in the early Universe. The merger rate predicted by such a scenario could be consistent with the one estimated from LIGO and Virgo, even if primordial black holes constitute a fraction of cold dark matter. The required abundance of primordial black holes is compatible with the existing upper limits from microlensing, caustic crossing and cosmic microwave background observations.

The eclipsing binary epsilon Aur is unique in being a very long-period binary involving an evolved, variable F star and a suspected B main-sequence star enshrouded in an opaque circumstellar disk. The geometrical arrangement is that the disk is viewed almost perfectly edge-on, with alignment leading to a partial eclipse of the F star. Despite a global observing campaign for the 2009-11 eclipse, there remain outstanding questions about the nature of the binary, its components, the disk, and the evolutionary state of the system. We analyze optical-band polarimetry in conjunction with broad-band color variations to interpret brightness variations across the surface of the F star. We model this both during and after the 1982-84 eclipse for which an extensive and dense data set exists. We develop a model in terms of surface temperature variations characterized by a small global variation overlaid with a temperature variation described with low-order spherical harmonics. While not providing a detailed fit to the dataset, our modeling captures the overall characterization of the color and polarimetric variability. In particular, we are able to recover the gross behavior of the polarimetric excursion in the Q-U plane as observed during eclipse of the F star when compared to post-eclipse behavior.

In this paper we show how effects from small scales enter the angular-redshift power spectrum $C_\ell(z,z')$. In particular, we show that spectroscopic surveys with high redshift resolution are affected by small scales already on large angular scales, i.e. at low multipoles. Therefore, when considering the angular power spectrum with spectroscopic redshift resolution, it is important to account for non-linearities relevant on small scales even at low multipoles. This may also motivate the use of the correlation function instead of the angular power spectrum. These effects, which are very relevant for bin auto-correlations, but not so important for cross-correlations, are quantified in detail.

Systematic effects arising from cosmic rays have been shown to be a significant threat to space telescopes using high-sensitivity bolometers. The LiteBIRD space mission aims to measure the polarised Cosmic Microwave Background with unprecedented sensitivity, but its positioning in space will also render it susceptible to cosmic ray effects. We present an end-to-end simulator for evaluating the expected scale of cosmic ray effect on the LiteBIRD space mission, which we demonstrate on a subset of detectors on the 166 GHz band of the Low Frequency Telescope. The simulator couples the expected proton flux at L2 with a model of the thermal response of the LFT focal plane and the electrothermal response of its superconducting detectors, producing time-ordered data which is projected into simulated sky maps and subsequent angular power spectra.

Tools for computing detailed optically thick spectral line profiles out of local thermodynamic equilibrium have always been focused on speed, due to the large computational effort involved. With the Lightweaver framework, we have produced a more flexible, modular toolkit for building custom tools in a high-level language, Python, without sacrificing speed against the current state of the art. The goal of providing a more flexible method for constructing these complex simulations is to decrease the barrier to entry and allow more rapid exploration of the field. In this paper we present an overview of the theory of optically thick NLTE radiative transfer, the numerical methods implemented in Lightweaver including the problems of time-dependent populations and charge-conservation, as well as an overview of the components most users will interact with, to demonstrate their flexibility.

Asymptotic giant branch (AGB) stars are known to lose a significant amount of mass by a stellar wind, which controls the remainder of their stellar lifetime. High angular-resolution observations show that the winds of these cool stars typically exhibit mid- to small-scale density perturbations such as spirals and arcs, believed to be caused by the gravitational interaction with a (sub-)stellar companion. We aim to explore the effects of the wind-companion interaction on the 3D density and velocity distribution of the wind, as a function of three key parameters: wind velocity, binary separation and companion mass. For the first time, we compare the impact on the outflow of a planetary companion to that of a stellar companion. We intend to devise a morphology classification scheme based on a singular parameter. With our grid of models we cover the prominent morphology changes in a companion-perturbed AGB outflow: slow winds with a close, massive binary companion show a more complex morphology. Additionally, we prove that massive planets are able to significantly impact the density structure of an AGB wind. We find that the interaction with a companion affects the terminal velocity of the wind, which can be explained by the gravitational slingshot mechanism. We distinguish between two types of wind focussing to the orbital plane resulting from distinct mechanisms: global flattening of the outflow as a result of the AGB star's orbital motion and the formation of an EDE as a consequence of the companion's gravitational pull. We investigate different morphology classification schemes and uncover that the ratio of the gravitational potential energy density of the companion to the kinetic energy density of the AGB outflow yields a robust classification parameter for the models presented in this paper.

We present a Hubble Space Telescope/Wide-Field Camera 3 near infrared spectrum of the archetype Y dwarf WISEP 182831.08+265037.8. The spectrum covers the 0.9-1.7 um wavelength range at a resolving power of lambda/Delta lambda ~180 and is a significant improvement over the previously published spectrum because it covers a broader wavelength range and is uncontaminated by light from a background star. The spectrum is unique for a cool brown dwarf in that the flux peaks in the Y, J, and H band are of near equal intensity in units of f_lambda. We fail to detect any absorption bands of NH_3 in the spectrum, in contrast to the predictions of chemical equilibrium models, but tentatively identify CH_4 as the carrier of an unknown absorption feature centered at 1.015 um. Using previously published ground- and spaced-based photometry, and using a Rayleigh Jeans tail to account for flux emerging longward of 4.5 um, we compute a bolometric luminosity of log (L_bol/L_sun)=-6.50+-0.02 which is significantly lower than previously published results. Finally, we compare the spectrum and photometry to two sets of atmospheric models and find that best overall match to the observed properties of WISEP 182831.08+265037.8 is a ~1 Gyr old binary composed of two T_eff~325 K, ~5 M_Jup brown dwarfs with subsolar [C/O] ratios.

Using large-scale fully-kinetic two-dimensional particle-in-cell simulations, we investigate the effects of shock rippling on electron acceleration at low-Mach-number shocks propagating in high-$\beta$ plasmas, in application to merger shocks in galaxy clusters. We find that the electron acceleration rate increases considerably when the rippling modes appear. The main acceleration mechanism is stochastic shock-drift acceleration, in which electrons are confined at the shock by pitch-angle scattering off turbulence and gain energy from the motional electric field. The presence of multi-scale magnetic turbulence at the shock transition and the region immediately behind the main shock overshoot is essential for electron energization. Wide-energy non-thermal electron distributions are formed both upstream and downstream of the shock. The maximum energy of the electrons is sufficient for their injection into diffusive shock acceleration. We show for the first time that the downstream electron spectrum has a~power-law form with index $p\approx 2.5$, in agreement with observations.

ExoFit trials are field campaigns financed by ESA to test the Rosalind Franklin rover and to enhance collaboration practices between ExoMars working groups. During the first trial, a replicate of the ExoMars rover was remotely operated from Oxfordshire (United Kingdom) to perform a complex sequence of scientific operation at the Tabernas Desert (Spain). By following the ExoMars Reference Surface Mission (RSM), the rover investigated the Badlands subsoil and collected drill cores, whose analytical study was entrusted to the RLS (Raman Laser Spectrometer) team. The preliminary characterization of core samples was performed in-situ through the RLS Engineering and Qualification Model (EQM-2) and the Raman Demonstrator (RAD1), being this a new, portable emulator of the RLS. In-situ results where then complemented by laboratory analysis using the RLS ExoMars simulator and the commercial version of the Curiosity/CheMin XRD system. Raman data, obtained by closely simulating the operational constraints of the mission, successfully disclosed the mineralogical composition of the samples, reaching the detection of minor/trace phases that were not detected by XRD. More importantly, Raman analysis detected many organic functional groups, proving the presence of extremophile organisms in the arid sub-surface of the Tabernas Desert. In light of the forthcoming ExoMars mission, the results here presented proves that RLS could play a critical role in the characterization of Martian sub-surface environments and in the analytical detection of potential traces of live tracers.

In some short gamma-ray bursts, precursor flares occurring $\sim$ seconds prior to the main episode have been observed. These flares may then be associated with the last few cycles of the inspiral when the orbital frequency is a few hundred Hz. During these final cycles, tidal forces can resonantly excite quasi-normal modes in the inspiralling stars, leading to a rapid increase in their amplitude. It has been shown that these modes can exert sufficiently strong strains onto the neutron star crust to instigate yieldings. Due to the typical frequencies of $g$-modes being $\sim 100\text{ Hz}$, their resonances with the orbital frequency match the precursor timings and warrant further investigation. Adopting realistic equations of state and solving the general-relativistic pulsation equations, we study $g$-mode resonances in coalescing quasi-circular binaries, where we consider various stellar rotation rates, degrees of stratification, and magnetic field structures. We show that for some combination of stellar parameters, the resonantly excited $g_1$- and $g_2$-modes may lead to crustal failure and trigger precursor flares.

The existence of a kinematic morphology-density relation remains uncertain, and instead stellar mass appears the more dominant driver of galaxy kinematics. We investigate the dependence of the stellar spin parameter proxy $\lambda_{R_e}$ on environment using a marked cross-correlation method with data from the SAMI Galaxy Survey. Our sample contains 710 galaxies with spatially resolved stellar velocity and velocity dispersion measurements. By utilising the highly complete spectroscopic data from the GAMA survey, we calculate marked cross-correlation functions for SAMI galaxies using a pair count estimator and marks based on stellar mass and $\lambda_{R_e}$. We detect an anti-correlation of stellar kinematics with environment at the 3.2$\sigma$ level, such that galaxies with low $\lambda_{R_e}$ values are preferably located in denser galaxy environments. However, a significant correlation between stellar mass and environment is also found (correlation at 2.4$\sigma$), as found in previous works. We compare these results to mock-observations from the cosmological EAGLE simulations, where we find a similar significant $\lambda_{R_e}$ anti-correlation with environment, and a mass and environment correlation. We demonstrate that the environmental correlation of $\lambda_{R_e}$ is not caused by the mass-environment relation. The significant relationship between $\lambda_{R_e}$ and environment remains when we exclude slow rotators. The signals in SAMI and EAGLE are strongest on small scales (10-100 kpc) as expected from galaxy interactions and mergers. Our work demonstrates that the technique of marked correlation functions is an effective tool for detecting the relationship between $\lambda_{R_e}$ and environment.

We create a continuous series of daily and monthly hemispheric sunspot numbers (HSNs) from 1874 to 2020, which will be continuously expanded in the future with the HSNs provided by SILSO. Based on the available daily measurements of hemispheric sunspot areas from 1874 to 2016 from Greenwich Royal Observatory and NOAA, we derive the relative fractions of the northern and southern activity. These fractions are applied to the international sunspot number (ISN) to derive the HSNs. This method and obtained data are validated against published HSNs for the period 1945--2020. We provide a continuous data series and catalogue of daily, monthly mean, and 13-month smoothed monthly mean HSNs for the time range 1874--2020 that are consistent with the newly calibrated ISN. Validation of the reconstructed HSNs against the direct data available since 1945 reveals a high level of consistency, with a correlation of r=0.94 (0.97) for the daily (monthly) data. The cumulative hemispheric asymmetries for cycles 12-24 give a mean value of 16%, with no obvious pattern in north-south predominance over the cycle evolution. The strongest asymmetry occurs for cycle no. 19, in which the northern hemisphere shows a cumulated predominance of 42%. The phase shift between the peaks of solar activity in the two hemispheres may be up to 28 months, with a mean absolute value of 16.4 months. The phase shifts reveal an overall asymmetry of the northern hemisphere reaching its cycle maximum earlier (in 10 out of 13 cases). Relating the ISN and HSN peak growth rates during the cycle rise phase with the cycle amplitude reveals higher correlations when considering the two hemispheres individually, with r = 0.9. Our findings demonstrate that empirical solar cycle prediction methods can be improved by investigating the solar cycle dynamics in terms of the hemispheric sunspot numbers.

A novel 4-dimensional Einstein-Gauss-Bonnet (4D EGB) gravity has been proposed that asserts to bypass the Lovelock's theorem and to result in a non-trivial contribution to the gravitational dynamics in four dimensions. In this work, we study the integrated Sachs-Wolfe (ISW) effect in the 4D EGB model. In this regard, we calculate the evolution of the gravitational potential and linear growth factor as a function of redshift for the 4D EGB model and compare it with the corresponding result obtained from the $\Lambda$-cold dark matter ($\Lambda$CDM) model. We also calculate the ISW auto-power-spectrum and the ISW cross-power-spectrum as functions of cosmic microwave background (CMB) multipoles for the 4D EGB model and compare those with the one obtained from the $\Lambda$CDM one. To do this, we use the strongest constraint for the coupling constant that has been proposed for the 4D EGB model. To measure the ISW effect for the 4D EGB model, we employ three large-scale structure surveys from different wavelengths. The results display that the ISW effect in the 4D EGB model is higher than the one obtained from the $\Lambda$CDM model. Moreover, we show the 4D EGB model can amplify the ISW cross-power-spectrum that can be considered as a relative advantage of the 4D EGB model. Also, we indicate that the deviation from the $\Lambda$CDM model is proportional to the value of $\beta$.

Exploding granules have drawn renewed interest because of their interaction with the magnetic field. Especially the newly forming downflow lanes developing in their centre seem to be eligible candidates for the intensification of magnetic fields. We analyse spectroscopic data from two different instruments in order to study the intricate velocity pattern within the newly forming downflow lanes in detail. We aim to examine general properties of a number of exploding granules. To gain a better understanding of the formation process of the developing intergranular lane in exploding granules, we study the temporal evolution and height dependence of the line-of-sight velocities at their formation location. Additionally, we search for evidence that exploding granules act as acoustic sources. We investigated the evolution of several exploding granules using data taken with the Interferometric Bidimensional Spectrometer and the Imaging Magnetograph eXperiment. Velocities for different heights of the solar atmosphere were determined by computing bisectors of the Fe I 6173.0{\AA} and the Fe I 5250.2{\AA} lines. We performed a wavelet analysis to study the intensity and velocity oscillations within and around exploding granules. We also compared our findings with predictions of numerical simulations. We found that exploding granules have significantly longer lifetimes than regular granules. Exploding granules larger than 3.8 arcsec form an independent intergranular lane during their decay phase, while smaller granules usually fade away or disappear into the intergranular area. For all exploding granules that form a new intergranular downflow lane, we find a temporal height-dependent shift with respect to the maximum of the downflow velocity. Our suggestion that this results from a complex atmospheric structure within the newly forming downflow lane is supported by the simulations.

We study the interaction of three solar wind structures, two stream interaction regions and one interplanetary coronal mass ejection, with Mars' plasma environment during 20-27 November 2007. This period corresponds to the solar minimum between the solar cycles 23 and 24 which was characterized by very low values of the solar wind density and dynamic pressure and low IMF magnitude. During that time the Mars-Express orbit was in the terminator plane, while the Earth, Sun, and Mars were almost aligned, so we use the ACE and STEREO probes as solar wind monitors in order to identify and characterize the structures that later hit Mars. We find that the passage of these structures caused strong variations of in the bow shock location (between 2.2 and 3.0~R$_M$), compression of the magnetospheric cavity (up to 45~\%) and an increased transterminator flow below 2~R$_M$ (by a factor of $\leq$8). This study shows that during times of low solar activity, modest space weather phenomena may cause large variations of plasma flow at Mars.

In order to address the question of whether spiral disturbances in galaxy discs are gravitationally coupled to the halo, we conduct simulations of idealized models of disc galaxies. We compare growth rates of spiral instabilities in identical mass models in which the halo is held rigid or is represented by particles drawn from an equilibrium distribution function. We examine cases of radial and azimuthal bias in the halo velocity ellipsoid in one of our models, and an isotropic velocity distribution in both. We find at most marginal evidence for an enhanced growth rate of spiral modes caused by a halo supporting response. We also find evidence for very mild dynamical friction between the spiral disturbance and the halo. We offer an explanation to account for the different behaviour between spiral modes and bar modes, since earlier work had found that bar instabilities became significantly more vigorous when a responsive halo was substituted for an equivalent rigid mass distribution. The barely significant differences found here justify the usual simplifying approximation of a rigid halo made in studies of spiral instabilities in galaxies.

Spinning black holes create electromagnetic storms when immersed in ambient magnetic fields, illuminating the otherwise epically dark terrain. In an electromagnetic extension of the Penrose process, tremendous energy can be extracted, boosting the energy of radiating particles far more efficiently than the mechanical Penrose process. We locate the regions from which energy can be mined and demonstrate explicitly that they are no longer restricted to the ergosphere. We also show that there can be toroidal regions that trap negative energy particles in orbit around the black hole. We find that the effective charge coupling between the black hole and the super-radiant particles decreases as energy is extracted, much like the spin of a black hole decreases in the mechanical analogue. While the effective coupling decreases, the actual charge of the black hole increases in magnitude reaching the energetically-favored Wald value, at which point energy extraction is impeded. We demonstrate the array of orbits for products from the electromagnetic Penrose process.

Our understanding of quantum correlators in cosmological spacetimes, including those that we can observe in cosmological surveys, has improved qualitatively in the past few years. Now we know many constraints that these objects must satisfy as consequences of general physical principles, such as symmetries, unitarity and locality. Using this new understanding, we derive the most general scalar four-point correlator, i.e., the trispectrum, to all orders in derivatives for manifestly local contact interactions. To obtain this result we use techniques from commutative algebra to write down all possible scalar four-particle amplitudes without assuming invariance under Lorentz boosts. We then input these amplitudes into a contact reconstruction formula that generates a contact cosmological correlator in de Sitter spacetime from a contact scalar or graviton amplitude. We also show how the same procedure can be used to derive higher-point contact cosmological correlators. Our results further extend the reach of the boostless cosmological bootstrap and build a new connection between flat and curved spacetime physics.

Low-energy effective field theories containing a light scalar field are used extensively in cosmology, but often there is a tension between embedding such theories in a healthy UV completion and achieving a phenomenologically viable screening mechanism in the IR. Here, we identify the range of interaction couplings which allow for a smooth resummation of classical non-linearities (necessary for kinetic/Vainshtein-type screening), and compare this with the range allowed by unitarity, causality and locality in the underlying UV theory. The latter region is identified using positivity bounds on the $2\to2$ scattering amplitude, and in particular by considering scattering about a non-trivial background for the scalar we are able to place constraints on interactions at all orders in the field (beyond quartic order). We identify two classes of theories can both exhibit screening and satisfy existing positivity bounds, namely scalar-tensor theories of $P(X)$ or quartic Horndeski type in which the leading interaction contains an odd power of $X$. Finally, for the quartic DBI Galileon (equivalent to a disformally coupled scalar in the Einstein frame), the analogous resummation can be performed near two-body systems and imposing positivity constraints introduces a non-perturbative ambiguity in the screened scalar profile. These results will guide future searches for UV complete models which exhibit screening of fifth forces in the IR.

Hawking evaporation of black holes in the early Universe is expected to copiously produce all kinds of particles, regardless of their charges under the Standard Model gauge group. For this reason, any fundamental particle, known or otherwise, could be produced during the black hole lifetime. This certainly includes dark matter (DM) particles. This paper improves upon previous calculations of DM production from primordial black holes (PBH) by consistently including the greybody factors, and by meticulously tracking a system of coupled Boltzmann equations. We show that the initial PBH densities required to produce the observed relic abundance depend strongly on the DM spin, varying in about $\sim 2$ orders of magnitude between a spin-2 and a scalar DM in the case of non-rotating PBHs. For Kerr PBHs, we have found that the expected enhancement in the production of bosons reduces the initial fraction needed to explain the measurements. We further consider indirect production of DM by assuming the existence of additional and unstable degrees of freedom emitted by the evaporation, which later decay into the DM. For a minimal setup where there is only one heavy particle, we find that the final relic abundance can be increased by at most a factor of $\sim 4$ for a scalar heavy state and a Schwarzschild PBH, or by a factor of $\sim 4.3$ for a spin-2 particle in the case of a Kerr PBH.

In this paper, we study how the evaporation of PBHs can affect the production of dark matter (DM) particles through thermal processes. We consider fermionic DM interacting with Standard Model particles via a spin-1 mediator in the context of a Freeze-Out or Freeze-In mechanism. We show that when PBHs evaporate after dominating the Universe's energy density, PBHs act as a source of DM and continuously inject entropy into the visible sector that can affect the thermal production in three qualitatively different ways. We compute the annihilation cross-sections which account for the interactions between and within the PBH produced and thermally produced DM populations, and establish a set of Boltzmann equations which we solve to obtain the correct relic abundance in those different regimes and confront the results with a set of different cosmological constraints. We provide analytic formulae to calculate the relic abundance for the Freeze-Out and Freeze-In mechanism in a PBH dominated early Universe. We identify regions of the parameter space where the PBHs dilute the relic density and thermalization occurs. Furthermore, we have made our code that numerically solves the Boltzmann equations publicly available.

In this work we study static black holes in the regularized 4D Einstein-Gauss-Bonnet theory of gravity; a shift-symmetric scalar-tensor theory that belongs to the Horndeski class. This theory features a simple black hole solution that can be written in closed form, and which we show is the unique static, spherically-symmetric and asymptotically-flat black hole vacuum solution of the theory. We further show that no asymptotically-flat, time-dependent, spherically-symmetric perturbations to this geometry are allowed, which suggests that it may be the only spherically-symmetric vacuum solution that this theory admits (a result analogous to Birkhoff's theorem). Finally, we consider the thermodynamic properties of these black holes, and find that their final state after evaporation is a remnant with a size determined by the coupling constant of the theory. We speculate that remnants of this kind from primordial black holes could act as dark matter, and we constrain the parameter space for their formation mass, as well as the coupling constant of the theory.

We study a cross-shaped cavity filled with superfluid $^4$He as a prototype resonant-mass gravitational wave detector. Using a membrane and a re-entrant microwave cavity as a sensitive optomechanical transducer, we were able to observe the thermally excited high-$Q$ acoustic modes of the helium at 20 mK temperature and achieved a strain sensitivity of $8 \times 10^{-19}$ Hz$^{-1/2}$ to gravitational waves. To facilitate the broadband detection of continuous gravitational waves, we tune the kilohertz-scale mechanical resonance frequencies up to 173 Hz/bar by pressurizing the helium. With reasonable improvements, this architecture will enable the search for GWs in the 1-30 kHz range, relevant for a number of astrophysical sources both within and beyond the Standard Model.

The precise calibration of the strain readout of the LIGO gravitational wave observatories is paramount to the accurate interpretation of gravitational wave events. This calibration is traditionally done by imparting a known force on the test masses of the observatory via radiation pressure. Here we describe the implementation of an alternative calibration scheme: the Newtonian Calibrator. This system uses a rotor consisting of both quadrupole and hexapole mass distributions to apply a time-varying gravitational force on one of the observatory's test masses. The force produced by this rotor can be predicted to $<1\%$ relative uncertainty and is well-resolved in the readout of the observatory. This system currently acts as a cross-check of the existing absolute calibration system.

We present experimental results using a single-electron resolution Skipper-CCD running above ground level to demonstrate the potential of this technology for reactor neutrinos and other low-energy particle interactions experiments. Operating conditions and data selection criteria are provided to decouple most of the background rate at low energies. Our final results for events with energies as low as 5 ionized electron-hole pairs show that the exponentially increasing rate of events seen in other technologies is not present in our data.

We study inflationary universes with an SU(3) gauge field coupled to an inflaton through a gauge kinetic function. Although the SU(3) gauge field grows at the initial stage of inflation due to the interaction with the inflaton, nonlinear self-couplings in the kinetic term of the gauge field become significant and cause nontrivial dynamics after sufficient growth. We investigate the evolution of the SU(3) gauge field numerically and reveal attractor solutions in the Bianchi type I spacetime. In general cases where all the components of the SU(3) gauge field have the same magnitude initially, they all tend to decay eventually because of the nonlinear self-couplings. Therefore, the cosmic no-hair conjecture generically holds in a mathematical sense. Practically, however, the anisotropy can be generated transiently in the early universe, even for an isotropic initial condition. Moreover, we find particular cases for which several components of the SU(3) gauge field survive against the nonlinear self-couplings. It occurs due to flat directions in the potential of a gauge field for Lie groups whose rank is higher than one. Thus, an SU(2) gauge field has a specialty among general non-Abelian gauge fields.

The temporal coincidences of events detected in four neutrino detectors and two gravitational antennas still remains among the most puzzling phenomena associated with SN1987A. The coincidences form a six-hour signal approximately coincident in time with the well-known LSD signal at 2h52m UT on 23/02/1987. After 30 years of research, the characteristics and the shape of the six-hour signal have been studied quite well, but the mechanisms of its formation have not been fully understood as of yet. Here we suggest that data obtained from another technology, radioactive decays, might provide new insights into the origin of signals previously seen in neutrino detectors and gravity wave detectors. On August 17, 2017, at 12h41m UT, the GW170817 signal was detected by LIGO and Virgo. At the same time, an approximately 7-hour long signal coincident with GW170817 was detected in the Si/Cl experiment on precision measurement of the $^{32}$Si half-life. We show that the Si/Cl signal is unexpectedly similar to the six-hour signal from SN1987A. In addition, we establish that the sources of the coinciding events are similar to those of the Si/Cl signal. To explain the surprising similarities in both signals, we present a mechanism which could in principle account for this phenomenon in terms of a local increase in the density of axionic dark matter induced by a gravity wave.

In this paper, we investigate the four-dimensional Einstein-Gauss-Bonnet black hole. The thermodynamic variables and equations of state of black holes are obtained in terms of a new parameterization. We discuss a formulation of the van der Waals equation by studying the effects of the temperature on P-V isotherms. We show the influence of the Cauchy horizon on the thermodynamic parameters. We prove by different methods, that the black hole entropy obey area law (plus logarithmic term that depends on the Gauss-Bonnet coupling {\alpha}). We propose a physical meaning for the logarithmic correction to the area law. This work can be extended to the extremal EGB black hole, in that case, we study the relationship between compressibility factor, specific heat and the coupling {\alpha}.

In the current 3nu paradigm, flavor oscillations probe 3 mixing angles (theta_12, theta_23, theta_13), one CP phase delta, and two squared mass differences delta m^2>0 and Delta m^2, where sign(Delta m^2)=+ (-) for normal (inverted) ordering. Absolute nu masses can be probed by the effective m_beta in beta decay, by the total mass Sigma in cosmology and, if neutrinos are Majorana, by another effective m_{beta beta} in 0nu2beta decay. Within an updated global analysis of (non)oscillation data, we constrain these 3nu parameters, both separately and in selected pairs, and highlight the concordance or discordance among different constraints. Five oscillation parameters (delta m^2, Delta m^2, theta_12, theta_23, theta_13) are consistently measured, with an overall accuracy ranging from ~1% for Delta m^2 to ~6% for sin^2(theta_23) (due to its octant ambiguity). We find overall hints for normal ordering (at 2.5 sigma), as well as for theta_23<pi/4 and for sin(delta)<0 (both at 90% C.L.), and discuss some tensions among datasets. Concerning nonoscillation data, we include the recent KATRIN constraints on m_beta, and we combine the latest 76-Ge, 130-Te and 136-Xe bounds on m_{beta beta}, accounting for NME covariances. We also discuss some variants related to CMB anisotropy and lensing data, which may affect cosmological constraints on Sigma and hints on sign(Delta m^2). The default option, including all Planck results, irrespective of the lensing anomaly, sets upper bounds on Sigma at the level of ~10^-1 eV, and further favors normal ordering up to ~3 sigma. An alternative option, that includes recent ACT results + other independent results (from WMAP and selected Planck data) globally consistent with standard lensing, is insensitive to the ordering but prefers Sigma ~(few) x 10^-1 eV, with different implications for m_beta and m_{beta beta} searches. (Abridged)

We discuss that it is quite possible to realize the smooth transition of the universe between a matter/radiation dominated deceleration and a dark energy dominated acceleration, even with a variation of proton-to-electron mass ratio. The variation is incorporated into the theory of gravity using a cosmological Higgs scalar field with a non-trivial self-interaction potential, leading to a varying Higgs vacuum expectation value (VEV). This matches well with the data from molecular absorption spectra of a series of Quasars. In comparison with late-time cosmology, an observational consistency is reached using a Markov chain Monte Carlo simulation and JLA+OHD+BAO data sets. We find that the pattern of variation is embedded within the evolving Equation of State (EOS) of the scalar Dark Energy/Matter components, but leaves negligible trace on the effective EOS of the system. We discuss three cases of scalar extended theory of gravity, (a) a minimally coupled scalar, (b) a non-minimally coupled scalar and (c) a generalized Brans-Dicke setup. We also give a toy model of a unified cosmic history from inflation to the present era and discuss how the Higg VEV might have changed as a function of look back time.

The Italian institute for nuclear physics (INFN) has financed the SIMP project (2019-2021) in order to strengthen its skills and technologies in the field of meV detectors with the ultimate aim of developing a single microwave photon detector. This goal will be pursued by improving the sensitivity and the dark count rate of two types of photodetectors: current biased Josephson Junction (JJ) for the frequency range 10-50 GHz and Transition Edge Sensor (TES) for the frequency range 30-100 GHz. Preliminary results on materials and devices characterization are presented.

We report on an all-sky search for continuous gravitational waves in the frequency band 20-2000\,Hz and with a frequency time derivative in the range of $[-1.0, +0.1]\times10^{-8}$\,Hz/s. Such a signal could be produced by a nearby, spinning and slightly non-axisymmetric isolated neutron star in our galaxy. This search uses the LIGO data from the first six months of Advanced LIGO's and Advanced Virgo's third observational run, O3. No periodic gravitational wave signals are observed, and 95\%\ confidence-level (CL) frequentist upper limits are placed on their strengths. The lowest upper limits on worst-case (linearly polarized) strain amplitude $h_0$ are $~1.7\times10^{-25}$ near 200\,Hz. For a circularly polarized source (most favorable orientation), the lowest upper limits are $\sim6.3\times10^{-26}$. These strict frequentist upper limits refer to all sky locations and the entire range of frequency derivative values. For a population-averaged ensemble of sky locations and stellar orientations, the lowest 95\%\ CL upper limits on the strain amplitude are $\sim1.\times10^{-25}$. These upper limits improve upon our previously published all-sky results, with the greatest improvement (factor of $\sim$2) seen at higher frequencies, in part because quantum squeezing has dramatically improved the detector noise level relative to the second observational run, O2. These limits are the most constraining to date over most of the parameter space searched.