Papers for Friday, Feb 19 2021

At redshift $z>5$ the far-infrared (FIR) continuum spectra of main-sequence galaxies are sparsely sampled, often with a single data point. The dust temperature $T_{\rm d, SED}$ thus has to be assumed in the FIR continuum fitting. This introduces large uncertainties regarding the derived dust mass ($M_{\rm d}$), FIR luminosity, and obscured fraction of the star formation rate. These are crucial quantities to quantify the effect of dust obscuration in high-$z$ galaxies. To overcome observations limitations, we introduce a new method that combines dust continuum information with the overlying $\hbox{[C $\scriptstyle\rm II $]} 158\mu$m line emission. By breaking the $M_{\rm d} - T_{\rm d, SED}$ degeneracy, with our method, we can reliably constrain the dust temperature with a single observation at $158\mu$m. This method can be applied to all ALMA and NOEMA $\hbox{[C $\scriptstyle\rm II $]}$ observations and exploited in ALMA Large Programs such as ALPINE and REBELS targeting $\hbox{[C $\scriptstyle\rm II $]}$ emitters at high-$z$. We also provide a physical interpretation of the empirical relation recently found between $molecular$ gas mass and $\hbox{[C $\scriptstyle\rm II $]}$ luminosity. We derive an analogous relation linking the $total$ gas surface density and $\hbox{[C $\scriptstyle\rm II $]}$ surface brightness. By combining the two, we predict the cosmic evolution of the surface density ratio $\Sigma_{\rm H_2} / \Sigma_{\rm gas}$. We find that $\Sigma_{\rm H_2} / \Sigma_{\rm gas}$ slowly increases with redshift, which is compatible with current observations at $0 < z < 4$.

Controversial studies on the jet collimation profile of BL Lacertae (BL Lac), the eponymous blazar of BL Lac objects class, complicate the scenario in this already puzzling class of objects. Understanding the jet geometry, in connection with the jet kinematics and the physical conditions in the surrounding medium, is fundamental to better constrain the formation, acceleration and collimation mechanisms in extragalactic jets. With the aim of investigating the jet geometry in the innermost regions of the jet of BL Lac, and solving the controversy, we explore the radio jet in this source, using high resolution millimeter-wave VLBI data. We collect 86GHz GMVA and 43GHz VLBA data to obtain stacked images that we use to infer the jet collimation profile by means of two comparable methods. We analyze the kinematics at 86GHz, and we discuss it in the context of the jet expansion. Finally we consider a possible implication of the Bondi sphere in shaping the different expanding region observed along the jet. We found that the jet in BL Lac expands with an overall conical geometry. A higher expanding rate region is observed between ~5 and 10 pc (de-projected) from the black hole. Such a region is associated with the decrease in brightness usually observed in high-frequency VLBI images of BL Lac. The jet retrieves the original jet expansion around 17 pc, where the presence of a recollimation shock is supported by both the jet profile and the 15GHz kinematics (MOJAVE survey). The change in the jet expansion profile occurring at ~5 pc could be associated with a change in the external pressure profile in correspondence of the Bondi radius (~3.3X10$^5$$R_s$).

We present a new method for inferring the gravitational potential of the Galactic disk, using the time-varying structure of a phase-space spiral in the $(z,w)$-plane (where $z$ and $w$ represent vertical position and vertical velocity). Our method of inference extracts information from the shape of the spiral, and disregards the bulk density distribution that is usually used to perform dynamical mass measurements. In this manner, it is complementary to traditional methods that are based on the assumption of a steady state. Our method consists of fitting an analytical model for the phase-space spiral to data, where the spiral is seen as a perturbation of the stellar number density in the $(z,w)$-plane. We tested our method on one-dimensional simulations, which were initiated in a steady state and then perturbed by an external force similar to that of a passing satellite. We were able to retrieve the true gravitational potentials of the simulations with high accuracy. The gravitational potential at 400-500 parsec distance from the disk mid-plane was inferred with an error of only a few percent. This is a first paper of a series in which we plan to test and refine our method on more complex simulations, as well as apply our method to Gaia data.

We investigate the infrared (IR) emission of high-redshift ($z\sim 6$), highly star-forming (${ {\rm SFR} > 100}$ $M_{\rm \odot} {\rm yr}^{-1}$) galaxies, with/without Active Galactic Nuclei (AGN), using a suite of cosmological simulations featuring dust radiative transfer. Synthetic Spectral Energy Distributions (SEDs) are used to quantify the relative contribution of stars/AGN to dust heating. In dusty (${M_{\rm d}\gtrsim 3\times 10^7 M_{\rm \odot}}$) galaxies, $\gtrsim 50-90 \%$ of the UV radiation is obscured by dust inhomogeneities on scales ${\gtrsim 100}$ pc. In runs with AGN, a clumpy, warm ($\approx 250$ K) dust component co-exists with a colder ($\approx 60$ K) and more diffuse one, heated by stars. Warm dust provides up to ${50 \%}$ of the total IR luminosity, but only $\lesssim 0.1 \%$ of the total mass content. The AGN boosts the MIR flux by ${10-100 \times}$ with respect to star forming galaxies, without significantly affecting the FIR. Our simulations successfully reproduce the observed SED of bright (${M_{\rm UV}\sim -26}$) ${z\sim 6}$ quasars, and show that these objects are part of complex, dust-rich merging systems, containing multiple sources (accreting BHs and/or star forming galaxies) in agreement with recent HST and ALMA observations. Our results show that the proposed ORIGINS missions will be able to investigate the MIR properties of dusty star forming galaxies and to obtain good quality spectra of bright quasars at $z\sim 6$. Finally, the MIR-to-FIR flux ratio of faint (${M_{\rm UV}\sim -24}$) AGN is ${>10\times}$ higher than for normal star forming galaxies. This implies that combined JWST/ORIGINS/ALMA observations will be crucial to identify faint and/or dust-obscured AGN in the distant Universe.

GW190814 was a compact object binary coalescence detected in gravitational waves by Advanced LIGO and Advanced Virgo that garnered exceptional community interest due to its excellent localization and the uncertain nature of the binary's lighter-mass component (either the heaviest known neutron star, or the lightest known black hole). Despite extensive follow up observations, no electromagnetic counterpart has been identified. Here we present new radio observations of 75 galaxies within the localization volume at $\Delta t\approx 35-266$ days post-merger. Our observations cover $\sim32$% of the total stellar luminosity in the final localization volume and extend to later timescales than previously-reported searches, allowing us to place the deepest constraints to date on the existence of a radio afterglow from a highly off-axis relativistic jet launched during the merger (assuming that the merger occurred within the observed area). For a viewing angle of $\sim46^{\circ}$ (the best-fit binary inclination derived from the gravitational wave signal) and assumed electron and magnetic field energy fractions of $\epsilon_e=0.1$ and $\epsilon_B=0.01$, we can rule out a typical short gamma-ray burst-like Gaussian jet with isotropic-equivalent kinetic energy $2\times10^{51}$ erg propagating into a constant density medium $n\gtrsim0.01$ cm$^{-3}$. These are the first limits resulting from a galaxy-targeted search for a radio counterpart to a gravitational wave event, and we discuss the challenges, and possible advantages, of applying similar search strategies to future events using current and upcoming radio facilities.

The abundance of dark matter haloes is one of the key probes of the growth of structure and expansion history of the Universe. Theoretical predictions for this quantity usually assume that, when expressed in a certain form, it depends only on the mass variance of the linear density field. However, cosmological simulations have revealed that this assumption breaks, leading to 10-20% systematic effects. In this paper we employ a specially-designed suite of simulations to further investigate this problem. Specifically, we carry out cosmological N-body simulations where we systematically vary growth history at a fixed linear density field, or vary the power spectrum shape at a fixed growth history. We show that the halo mass function generically depends on these quantities, thus showing a clear signal of non-universality. Most of this effect can be traced back to the way in which the same linear fluctuation grows differently into the nonlinear regime depending on details of its assembly history. With these results, we propose a parameterization with explicit dependence on the linear growth rate and power spectrum shape. Using an independent suite of simulations, we show that this fitting function accurately captures the mass function of haloes over cosmologies spanning a vast parameter space, including massive neutrinos and dynamical dark energy. Finally, we employ this tool to improve the accuracy of so-called cosmology-rescaling methods and show they can deliver 2% accurate predictions for the halo mass function over the whole range of currently viable cosmologies.

We studied the apparent galaxy pair NGC 1232 / NGC 1232A with Chandra, looking for evidence of interactions and collisions. We report that there is no cloud of diffuse emission in NGC 1232, contrary to previous claims in the literature. Instead, we find that the small "companion" galaxy NGC 1232A contains three ultraluminous X-ray sources (ULXs) with peak 0.3-10 keV luminosities above 10^{40} erg/s (assuming a cosmological distance of 93 Mpc for this galaxy). For its mass, morphology, metal abundance and bright ULX population, NGC 1232A is analogous to the more nearby late-type spiral NGC 1313.

The GeV break in spectra of the blazar 3C 454.3 is a special observation feature that has been discovered by the {\it Fermi}-LAT. The origin of the GeV break in the spectra is still under debate. In order to explore the possible source of GeV spectral break in 3C 454.3, a one-zone homogeneous leptonic jet model, as well as the {\it McFit} technique are utilized for fitting the quasi-simultaneous multi-waveband spectral energy distribution (SED) of 3C 454.3. The outside border of the broad-line region (BLR) and inner dust torus are chosen to contribute radiation in the model as external, seed photons to the external-Compton process, considering the observed $\gamma$-ray radiation. The combination of two components, namely the Compton-scattered BLR and dust torus radiation, assuming a broken power-law distribution of emitted particles, provides a proper fitting to the multi-waveband SED of 3C 454.3 detected 2008 Aug 3 - Sept 2 and explains the GeV spectral break. We propose that the spectral break of 3C 454.3 may originate from an inherent break in the energy distribution of the emitted particles and the Klein-Nishina effect. A comparison is performed between the energy density of the 'external' photon field for the whole BLR $U_{\rm BLR}$ achieved via model fitting and that constrained from the BLR data. The distance from the position of the $\gamma$-ray radiation area of 3C 454.3 to the central black hole could be constrained at $\sim 0.78$pc ($\sim 4.00 R_{\rm BLR}$, the size of the BLR).

Primordial supermassive stars (SMSs) formed in atomic-cooling halos at $z \sim$ 15 - 20 are leading candidates for the seeds of the first quasars. Past numerical studies of the evolution of SMSs have typically assumed constant accretion rates rather than the highly variable flows in which they form. We model the evolution of SMSs in the cosmological flows that create them using the Kepler stellar evolution and implicit hydrodynamics code. We find that they reach masses of $1 - 2 \times 10^5 M_{\odot}$ before undergoing direct-collapse to black holes (DCBHs) during or at the end of their main-sequence hydrogen burning, at 1 - 1.5 Myr, regardless of halo mass, spin, or merger history. Our models confirm that the accretion histories predicted for purely atomic-cooling halos impose a narrow spectrum of masses on the seeds of the first massive quasars. Our results also indicate that multiple SMSs at disparate stages of evolution can form in these halos, raising the possibility of SMS binaries and supermassive X-ray binaries (SMXBs), as well as DCBH mergers which could be detected by LISA.

We present high-resolution near-infrared spectra taken during eight transits of 55 Cancri e, a nearby low-density super-Earth with a short orbital period (< 18 hours). While this exoplanet's bulk density indicates a possible atmosphere, one has not been detected definitively. Our analysis relies on the Doppler cross-correlation technique, which takes advantage of the high spectral resolution and broad wavelength coverage of our data, to search for the thousands of absorption features from hydrogen-, carbon-, and nitrogen-rich molecular species in the planetary atmosphere. Although we are unable to detect an atmosphere around 55 Cancri e, we do place strong constraints on the levels of HCN, NH${}_3$, and C${}_2$H${}_2$ that may be present. In particular, at a mean molecular weight of 5 amu we can rule out the presence of HCN in the atmosphere down to a volume mixing ratio (VMR) of 0.02%, NH${}_3$ down to a VMR of 0.08%, and C${}_2$H${}_2$ down to a VMR of 1.0%. If the mean molecular weight is relaxed to 2 amu, we can rule out the presence of HCN, NH${}_3$, and C${}_2$H${}_2$ down to VMRs of 0.001%, 0.0025%, and 0.08% respectively. Our results reduce the parameter space of possible atmospheres consistent with the analysis of HST/WFC3 observations by Tsiaras et al. (2016), and indicate that if 55 Cancri e harbors an atmosphere, it must have a high mean molecular weight and/or clouds.

In this paper, we systematically study gravitational waves (GWs) produced by remote compact astrophysical sources. To describe such GWs properly, we introduce three scales, $\lambda, \; L_c$ and $L$, denoting, respectively, the typical wavelength of GWs, the scale of the cosmological perturbations, and the size of the observable universe. For GWs to be detected by the current and foreseeable detectors, the condition $\lambda \ll L_c \ll L$ holds, and such GWs can be well approximated as high-frequency GWs. In order for the backreaction of the GWs to the background to be negligible, we must assume that $\left|h_{\mu\nu}\right| \ll 1$, in addition to the condition $\epsilon \ll 1$, which are also the conditions for the linearized Einstein field equations for $h_{\mu\nu}$ to be valid, where $g_{\mu\nu} = \gamma_{\mu\nu} + \epsilon h_{\mu\nu}$, and $\gamma_{\mu\nu}$ denotes the background. To simplify the field equations, we show that the spatial, traceless, and Lorentz gauge conditions can be imposed simultaneously, even when the background is not vacuum, as long as the high-frequency GW approximation is valid. However, to develop the formulas that can be applicable to as many cases as possible, we first write down explicitly the linearized Einstein field equations by imposing only the spatial gauge. Applying the general formulas together with the geometrical optics approximation to such GWs, we find that they still move along null geodesics and its polarization bi-vector is parallel-transported, even when both the cosmological scalar and tensor perturbations are present. In addition, we also calculate the gravitational integrated Sachs-Wolfe effects, whereby the dependences of the amplitude, phase and luminosity distance of the GWs on these two kinds of perturbations are read out explicitly.

Analysis of self-similarity in scale-free $N$-body simulations reveals the spatial and temporal scales for which statistics measured in cosmological simulations are converged to the physical continuum limit. We examine how the range of scales in which the two-point correlation function is converged depends on the force softening length and whether it is held constant in comoving or proper coordinates. We find that a proper softening that reaches roughly 1/30th of the inter-particle spacing by the end of the simulation resolves the same spatial and temporal scales as a comoving softening of the same length while using a third fewer time steps, for a range of scale factors typical to $\Lambda$CDM simulations. We additionally infer an inherent resolution limit, set by the particle mass and scaling as $a^{-1/2}$, beyond which reducing the softening does not improve the resolution. We postulate a mapping of these results with spectral index $n=-2$ to $\Lambda$CDM simulations.

Context: The number of known strong gravitational lenses is expected to grow substantially in the next few years. The statistical combination of large samples of lenses has the potential of providing strong constraints on the inner structure of galaxies. Aims: We investigate to what extent we can calibrate stellar mass measurements and constrain the average dark matter density profile of galaxies by statistically combining strong lensing data from thousands of lenses. Methods: We generate mock samples of axisymmetric lenses. We assume that, for each lens, we have measurements of two image positions of a strongly lensed background source, as well as magnification information from full surface brightness modelling, and a stellar population synthesis-based estimate of the lens stellar mass. We then fit models describing the distribution of the stellar population synthesis mismatch parameter $\alpha_{sps}$ (the ratio between the true stellar mass and the stellar population synthesis-based estimate) and dark matter density profile of the population of lenses to an ensemble of 1000 mock lenses. Results: The average $\alpha_{sps}$, projected dark matter mass and dark matter density slope can be obtained with great precision and accuracy, compared with current constraints. A flexible model and the knowledge of the lens detection efficiency as a function of image configuration are required in order to avoid a biased inference. Conclusions: Statistical strong lensing inferences from upcoming surveys have the potential to calibrate stellar mass measurements and to constrain the inner dark matter density profile of massive galaxies.

$K2$ greatly extended $Kepler$'s ability to find new planets, but it was typically limited to identifying transiting planets with orbital periods below 40 days. While analyzing $K2$ data through the Exoplanet Explorers project, citizen scientists helped discover one super-Earth and four sub-Neptune sized planets in the relatively bright ($V=12.21$, $K=10.3$) K2-138 system, all which orbit near 3:2 mean motion resonances. The $K2$ light curve showed two additional transit events consistent with a sixth planet. Using $Spitzer$ photometry, we validate the sixth planet's orbital period of $41.966\pm0.006$ days and measure a radius of $3.44^{+0.32}_{-0.31}\,R_{\oplus}$, solidifying K2-138 as the $K2$ system with the most currently known planets. There is a sizeable gap between the outer two planets, since the fifth planet in the system, K2-138 f, orbits at 12.76 days. We explore the possibility of additional non-transiting planets in the gap between f and g. Due to the relative brightness of the K2-138 host star, and the near resonance of the inner planets, K2-138 could be a key benchmark system for both radial velocity and transit timing variation mass measurements, and indeed radial velocity masses for the inner four planets have already been obtained. With its five sub-Neptunes and one super-Earth, the K2-138 system provides a unique test bed for comparative atmospheric studies of warm to temperate planets of similar size, dynamical studies of near resonant planets, and models of planet formation and migration.

Brighter Type Ia supernovae (SNe Ia) prefer less massive host galaxies and regions of higher star formation, but the dependence on host property measurement technique on this bias is not known. Using the PISCO IFS sample of SN host galaxies alongside SDSS, GALEX, and 2MASS photometry, we compared host galaxy stellar mass and star formation rate (SFR) estimates using different observation methods and fitting techniques. Mass estimates were consistent regardless of fitting technique. For optical wavelengths, the magnitude correction size ("magnitude step") for mass estimates varied from -0.04 mag for FAST++ and STARLIGHT to -0.02 mag for ZPEG. Including UV information for mass estimates results in a mass-dependent shift in estimates of the individual galaxies but had no effect on measured mass step size or location. Regardless of technique, the measured magnitude step for global specific SFR (sSFR) was ~0.1 mag; UV-calibrated SFR estimates were consistent with those from H{\alpha}. Only 49 hosts were used to measure the sSFR step as opposed to 76 FAST++ mass estimates or 97 STARLIGHT mass estimates to measure the mass step. Local sSFR step model fits were inconclusive due to a diminished sample size of 38 SNe Ia. Stronger correlation between GALEX UV photometry global and local sSFR estimates relative H{\alpha}estimates due to a large GALEX PSF inhibits its use to measure local UV 1 kpc flux. Overall, observation method and fitting technique had little influence on host bias detection. Instead, of importance is the property used to measure the host bias.

Potential features in the primordial power spectrum, such as oscillatory patterns, have been searched for in galaxy surveys in recent years, since these features can assist in understanding the nature of inflation and distinguishing between different scenarios of inflation. The null detection to date suggests that any such features should be fairly weak, and next-generation galaxy surveys, with their unprecedented sizes and precisions, are in a position to place stronger constraints than before. However, even if such primordial features once existed in the early Universe, they would have been significantly weakened or even wiped out on small scales in the late Universe due to nonlinear structure formation, which makes them difficult to be directly detected in real observations. A potential way to tackle this challenge for probing the features is to undo the cosmological evolution, i.e., using reconstruction to obtain an approximate linear density field. By employing a suite of large N-body simulations, we show that a recently-proposed nonlinear reconstruction algorithm can effectively retrieve lost oscillatory features from the mock galaxy catalogues and improve the accuracy of the measurement of feature parameters (assuming such primordial features do exist). We do a Fisher analysis to forecast how reconstruction affects the constraining power, and find that it can lead to significantly more robust constraints on the oscillation amplitude for a DESI-like survey. In particular, we compare the application of reconstruction with other ways of improving constraints, such as increasing the survey volume and range of scales, and show that it can achieve what the latter do, but at a much lower cost.

Using a sample of 133 991 galaxies distributed in the sky region $100^{\circ} <\alpha<270^{\circ}$ and $7^{\circ}<\delta<65^{\circ}$, extracted from the SDSS NASA/AMES Value Added Galaxy Catalog (AMES-VAGC), we estimate the fractal dimension using two different methods. First, using an algorithm to estimate the correlation dimension. The second method, in a novel approach, creates a graph from the data and estimates the graph dimension purely from connectivity information. In both methods we found a dimension $D\approx 2$ in scales below 20 Mpc, which agrees with previous works. This result shows the non-homogeneity of galaxies distribution at certain scales.

The Galactic HII region luminosity function (LF) is an important metric for understanding global star formation properties of the Milky Way, but only a few studies have been done and all use relatively small numbers of HII regions. We use a sample of 797 first Galactic quadrant HII regions compiled from the WISE Catalog of Galactic HII Regions to examine the form of the LF at multiple infrared and radio wavelengths. Our sample is statistically complete for all regions powered by single stars of type O9.5V and earlier. We fit the LF at each wavelength with single and double power laws. Averaging the results from all wavelengths, the mean of the best-fit single power law index is $\langle\alpha\rangle=-1.75\,\pm\,0.01$. The mean best-fit double power law indices are $\langle\alpha_1\rangle=-1.40\,\pm\,0.03$ and $\langle\alpha_2\rangle=-2.33\,\pm\,0.04$. We conclude that neither a single nor a double power law is strongly favored over the other. The LFs show some variation when we separate the HII region sample into subsets by heliocentric distance, physical size, Galactocentric radius, and location relative to the spiral arms, but blending individual HII regions into larger complexes does not change the value of the power law indices of the best-fit LF models. The consistency of the power law indices across multiple wavelengths suggests that the LF is independent of wavelength. This implies that infrared and radio tracers can be employed in place of H$\alpha$.

In 1972, Zinn, Newell, & Gibson (ZNG) published a list of 156 candidate "UV-bright" stars they had found in 27 Galactic globular clusters (GCs), based on photographs in the U and V bands. UV-bright stars lie above the horizontal branch (HB) and blueward of the asymptotic-giant branch (AGB) and red-giant branch in the clusters' color-magnitude diagrams. They are in rapid evolutionary phases--if they are members and not unrelated bright foreground stars. The ZNG list has inspired numerous follow-up studies, aimed at understanding late stages of stellar evolution. However, the ZNG candidates were presented only in finding charts, and celestial coordinates were not given. Using my own collection of CCD frames in u and V, I have identified all of the ZNG objects, and have assembled their coordinates, parallaxes, and proper motions from the recent Gaia Early Data Release 3 (EDR3). Based on the Gaia astrometry, I have determined which objects are probable cluster members (45% of the sample). For the members, using photometry from EDR3, I have assigned the stars to various evolutionary stages, including luminous post-AGB stars, and stars above the HB. I point out several ZNG stars of special interest that have still, to my knowledge, never been studied in detail. This study is an adjunct to a forthcoming survey of the Galactic GCs in the uBVI photometric system, designed for detection of low-gravity stars with large Balmer discontinuities.

Earth Trojan Asteroids are an important but elusive population that co-orbit with Earth at the L4 and L5 Lagrange points. There is only one known, but a large population is theoretically stable and could provide insight into our solar system's past and present as well as planetary defense. In this paper, we present the results of an Earth Trojan survey that uses a novel shift-and-stack detection method on two nights of data from the Dark Energy Camera. We find no new Earth Trojan Asteroids. We calculate an upper limit on the population that is consistent with previous searches despite much less sky coverage. Additionally, we elaborate on previous upper limit calculations using current asteroid population statistics and an extensive asteroid simulation to provide the most up to date population constraints. We find an L4 Earth Trojan population of NET < 1 for H = 14, NET < 7 for H = 16, and NET < 642 for H = 22.

The Juno orbiter continues to collect data on Jupiter's gravity field with unprecedented precision since 2016, recently reporting a non-hydrostatic component in the tidal response of the planet. At the mid-mission perijove 17, Juno registered a Love number $k_2=0.565\pm0.006$ that is $-4\pm1\%$ ($1\sigma$) from the theoretical hydrostatic $k_2^{(hs)}=0.590$. Here we assess whether the aforementioned departure of tides from hydrostatic equilibrium represents the neglected gravitational contribution of dynamical tides. We employ perturbation theory and simple tidal models to calculate a fractional dynamical correction $\Delta k_2$ to the well-known hydrostatic $k_2$. Exploiting the analytical simplicity of a toy uniform-density model, we show how the Coriolis acceleration motivates the negative sign in the $\Delta k_2$ observed by Juno. By simplifying Jupiter's interior into a core-less, fully-convective, and chemically-homogeneous body, we calculate $\Delta k_2$ in a model following an $n=1$ polytrope equation of state. Our numerical results for the $n=1$ polytrope qualitatively follow the behaviour of the uniform-density model, mostly because the main component of the tidal flow is similar in each case. Our results indicate that the gravitational effect of the Io-induced dynamical tide leads to $\Delta k_2=-4\pm1\%$, in agreement with the non-hydrostatic component reported by Juno. Consequently, our results suggest that Juno obtained the first unambiguous detection of the gravitational effect of dynamical tides in a gas giant planet. These results facilitate a future interpretation of Juno tidal gravity data with the purpose of elucidating the existence of a dilute core in Jupiter.

We analyze observations of the mass profiles of 175 late-type galaxies in the Spitzer Photometry \& Accurate Rotation Curves (SPARC) database to construct the temperature profile of their dark matter (DM) halos by assuming that (1) DM in the halos obeys either the Fermi-Dirac or the Maxwell-Boltzmann distribution, and (2) the halos are in the virial state. We derive the dispersion velocity of DM at the center of the halos and show that its correlation with the halo's total mass is the same as the one estimated in N-body simulations. The correlation is also the same as the observed relation between the two variables for visible matter in galaxies. Taking the latter agreement as a validation of our analysis, we derive the mass to the temperature of DM at the edge of the halos and show that it is galaxy independent and is equal to $m/T_R\simeq 10^{10}$ in natural units. Since the analyzed galaxies are far away in the sky, we conclude that DM is a thermal relic, and $T_R$ in the above ratio can be expressed in terms of the temperature of the cosmic microwave background (CMB) at the time of DM decoupling. This result is used to study possible cosmological scenarios. We show that observations are at odds with (1) non-thermal DM, (2) hot DM, and (3) collision-less cold DM. Our findings are in favor of a warm DM with a mass of $\sim 1\,$MeV.

Recent observations of the high-mass X-ray binary Cygnus X-1 have shown that both the companion star (41 solar masses) and the black hole (21 solar masses) are more massive than previously estimated. Furthermore, the black hole appears to be nearly maximally spinning. Here we present a possible formation channel for the Cygnus X-1 system that matches the observed system properties. In this formation channel, we find that the orbital parameters of Cygnus X-1, combined with the observed metallicity of the companion, imply a significant reduction in mass loss through winds relative to commonly used prescriptions for stripped stars.

Cygnus X-1 is a well-studied persistent black hole X-ray binary. Recently, the three parameters needed to estimate the black hole spin of this system, namely the black hole mass $M$, the orbital inclination $i$ and the source distance $D$, have been updated. In this work we redetermine the spin parameter using the continuum-fitting technique for those updated parameter values. Based on the assumption that the spin axis of the black hole is aligned with the orbital plane, we fit the thermal disk component to a fully relativistic thin accretion disk model. The error in the spin estimate arising from the combined observational uncertainties is obtained via Monte Carlo (MC) simulations. We demonstrate that, without considering the counteracting torque effect, the new spin parameter is constrained to be a$_* > 0.9985$ (3$\sigma$), which confirms that the spin of the black hole in Cygnus X-1 is extreme.

Accretion is an essential physical process in black-hole X-ray binaries (BHXRBs) and active galactic nuclei. The properties of accretion flows and their radiation are generally considered to be uniquely determined by the mass accretion rate of the disk; however, the "hysteresis effect" observed during outbursts of nearly all BHXRBs seriously challenges this paradigm. The hysteresis effect is referred to that apparently similar spectral state transitions take place at very different luminosities during an outburst cycle. Phenomenologically, this effect is also represented as the so-called "q"-shaped hardness-intensity diagram (HID), which has been proposed as a unified scene for BHXRBs. However, even with the high cadence pointing observations, the distinctly important rapid-rise stage has been caught to date in only a few BHXRB outbursts, and it is still limited to narrow energy/wavelength bands. As a result, there is still lack of a quantitative theoretical interpretation and observational understanding on the "q"-diagram. Here, we present a detailed time-lag analysis on Insight-HXMT's intensive monitoring data of a newly found BHXRB, MAXI J1348-630, over a broad energy band (1--150 keV). We find the first observational evidence that the observed time-lag between radiations of the accretion disk and the corona leads naturally to the hysteresis effect and the "q"-diagram. Moreover, complemented by the quasi-simultaneous Swift data, we achieve a panorama of the accretion flow: the hard X-rays from the corona heat the outer disk and induce the corresponding optical brightening; thereafter, the enhanced accretion rate in the outer disk propagates inwards to the vicinity of the central black-hole, generating the soft X-rays in the inner disk region, at a viscous timescale of $\sim 8-12$ days.

We study molecular emission in a massive condensation at the border of the HII region RCW 120, paying particular attention to the Core 1 and Core 2 objects, the most massive fragments of the condensation found previously by ALMA. The latter fragment was previously suggested to host a high-mass analogue of Class 0 young stellar object. We present spectra of molecular emission in the 1 mm range made with the APEX telescope. We detect CH$_3$OH and C$^{34}$S lines in Core 1 and Core 2. The CH$_3$CN series and the SO$_2$ lines are only found in Core 2. We estimate gas physical parameters using methanol lines and obtain gas temperature less than 100 K in both regions. Molecular hydrogen number density in Core 2 is in the range of $10^5-10^7$ cm$^{-3}$ and is more uncertain in Core 1. However, the detection of the CH$_3$CN lines corresponding to highly excited transitions ($E_{\rm u}> 400$~K) in Core~2 indicates that the region contains hot gas, while the abundances of CH$_3$OH, CS, SO$_2$ and CH$_3$CN are quite low for a hot core stage. We propose that Core 2 is in the warm-up phase prior to the establishing of the hot gas chemistry. We suggest that Core 2 is in the beginning of the hot core stage. There are no detected CH$_3$CN lines in Core 1, therefore, it might be on an even less evolved evolutionary stage.

In circumstellar disks, the size of dust particles varies from submicron to several centimeters, while planetesimals have sizes of hundreds of kilometers. Therefore, various regimes for the aerodynamic drag between solid bodies and gas can be realized in these disks, depending on the grain sizes and velocities: Epstein, Stokes, and Newton, as well as transitional regimes between them. For small bodies moving in the Epstein regime, the time required to establish the constant relative velocity between the gas and bodies can be much less than the dynamical time scale for the problem - the time for the rotation of the disk about the central body. In addition, the dust may be concentrated in individual regions of the disk, making it necessary to take into account the transfer of momentum between the dust and gas. It is shown that, for a system of equations for gas and monodisperse dust, a semi-implicit first-order approximation scheme in time in which the interphase interaction is calculated implicitly, while other forces, such as the pressure gradient and gravity are calculated explicitly, is suitable for stiff problems with intense interphase interactions and for computations of the drag in non-linear regimes. The piece-wise drag coefficient widely used in astrophysical simulations has a discontinuity at some values of the Mach and Knudsen numbers that are realized in a circumstellar disk. A continuous drag coefficient is presented, which corresponds to experimental dependences obtained for various drag regimes.

Latest study reports that plasma emission can be generated by energetic electrons of DGH distribution via the electron cyclotron maser instability (ECMI) in plasmas characterized by a large ratio of plasma oscillation frequency to electron gyro-frequency ($\omega_{pe}/\Omega_{ce}$). In this study, on the basis of the ECMI-plasma emission mechanism, we examine the double plasma resonance (DPR) effect and the corresponding plasma emission at both harmonic (H) and fundamental (F) bands using PIC simulations with various $\omega_{pe}/\Omega_{ce}$. This allows us to directly simulate the feature of zebra pattern (ZP) observed in solar radio bursts for the first time. We find that (1) the simulations reproduce the DPR effect nicely for the upper hybrid (UH) and Z modes, as seen from their variation of intensity and linear growth rate with $\omega_{pe}/\Omega_{ce}$, (2) the intensity of the H emission is stronger than that of the F emission by $\sim$ 2 orders of magnitude and vary periodically with increasing $\omega_{pe}/\Omega_{ce}$, while the F emission is too weak to be significant, therefore we suggest that it is the H emission accounting for solar ZPs, (3) the peak-valley contrast of the total intensity of H is $\sim 4$, and the peak lies around integer values of $\omega_{pe}/\Omega_{ce}$ (= 10 and 11) for the present parameter setup. We also evaluate the effect of energy of energetic electrons on the characteristics of ECMI-excited waves and plasma radiation. The study provides novel insight on the physical origin of ZPs of solar radio bursts.

The scientific value of the next generation of large continuum surveys would be greatly increased if the redshifts of the newly detected sources could be rapidly and reliably estimated. Given the observational expense of obtaining spectroscopic redshifts for the large number of new detections expected, there has been substantial recent work on using machine learning techniques to obtain photometric redshifts. Here we compare the accuracy of the predicted photometric redshifts obtained from Deep Learning(DL) with the k-Nearest Neighbour (kNN) and the Decision Tree Regression (DTR) algorithms. We find using a combination of near-infrared, visible and ultraviolet magnitudes, trained upon a sample of SDSS QSOs, that the kNN and DL algorithms produce the best self-validation result with a standard deviation of {\sigma} = 0.24. Testing on various sub-samples, we find that the DL algorithm generally has lower values of {\sigma}, in addition to exhibiting a better performance in other measures. Our DL method, which uses an easy to implement off-the-shelf algorithm with no filtering nor removal of outliers, performs similarly to other, more complex, algorithms, resulting in an accuracy of {\Delta}z < 0.1$ up to z ~ 2.5. Applying the DL algorithm trained on our 70,000 strong sample to other independent (radio-selected) datasets, we find {\sigma} < 0.36 over a wide range of radio flux densities. This indicates much potential in using this method to determine photometric redshifts of quasars detected with the Square Kilometre Array.

Southern Spectroscopic Survey Telescope (SSST) is a wide-field spectroscopic survey telescope that China plans to build in Chile in the next few years. As an instrument for astronomical spectroscopic survey, the multi-object and fiber-fed spectrograph (MOFFS) is one of the most important scientific instruments for SSST. In this paper, we present a recommended optical design for the MOFFS system based on the Volume Phase Holographic Gratings (VPHG). The whole design philosophy and procedure, including the analytic method to determine the initial structure, optimization procedures of the VPHG and the camera groups, are demonstrated in detail. The numerical results of the final obtained spectrograph show a superior imaging quality and a relatively high transmittance for the whole working waveband and the field of view. The design method proposed in this paper can provide a reference for the design of MOFFS accommodated in spectroscopic survey telescopes.

We collect a sample of stars observed both in LAMOST and Gaia which have colors implying a temperature hotter than 7000 K. We train a machine learning algorithm on LAMOST spectroscopic data which has been tagged with stellar classifications and metallicities, and use this machine to construct a catalog of Blue Horizontal Branch stars (BHBs) with metallicity information. Another machine is trained using Gaia parallaxes to predict absolute magnitudes for these stars. The final catalog of 13,693 BHBs is thought to be about 86\% pure, with $\sigma_{[Fe/H]}\sim$0.35 dex and $\sigma_{G}\sim$0.31 mag. These values are confirmed via comparison to globular clusters, although a covariance error seems to affect our magnitude and abundance estimates. We analyze a subset of this catalog in the Galactic Halo. We find that BHB populations in the outer halo appear redder, which could imply a younger population, and that the metallicity gradient is relatively flat around [Fe/H] = -1.9 dex over our sample footprint. We find that our metal rich BHB stars are on more radial velocity dispersion dominated orbits ($\beta \sim 0.70$) at all radii than our metal poor BHB stars ($\beta \sim 0.62$).

We present rest-frame $B$ and $I$ imaging of 35 low-redshift ($z < 0.5$) Palomar-Green quasars using the Hubble Space Telescope Wide Field Camera 3. We perform multi-component two-dimensional image decomposition to separate the host galaxy from its bright active nucleus, characterize its morphology, and measure its photometric properties. Special care is devoted to quantifying the structural parameters of the galaxy bulge, determine its $B-I$ color, and estimate its stellar mass. Roughly half of the sample, comprising the less luminous ($L_{5100} \lesssim 10^{45}\,\mathrm{erg\,s^{-1}}$) but most high Eddington ratio quasars, reside in disk galaxies that are often barred and possess pseudo bulges. The large stellar masses, large effective radii, and faint surface brightnesses suggest that the host galaxies of the most luminous quasars are mostly ellipticals. Major mergers constitute only a minority ($\lesssim 20\%$) of our sample. Our quasar sample roughly obeys the scaling relations between black hole mass and host galaxy (bulge, core, total) stellar mass. Hosts with black holes more massive than $\sim 10^8\,M_\odot$ behave similarly to classical bulges and early-type galaxies, while those with less massive black holes, particular the narrow-line Seyfert 1s, are consistent with pseudo bulges in late-type galaxies. The host galaxy bulges, irrespective of whether they are classical or pseudo, follow the relatively tight inverse relation between effective radius and mean effective surface brightness of inactive classical bulges and ellipticals. We argue that pseudo bulges experience recent or ongoing nuclear star formation.

The quintessence field coupled to the cosmic neutrino background (CNB) has been widely discussed as an alternative mechanism to address the coincidence problem. As it is well known, it is possible to extend such models to obtain quintessential inflation, that is, to incorporate inflationary stage as well. Taking an alternative route, one can start from the well established inflationary models and obtain successful quintessence models at the expense of coupling with the CNB. To Follow this route, we use a slightly reformulated model addressed in PRD95, 123521 (2017). This particular model assumes $\mathcal{Z}_2$ symmetry for both scalar field potential and coupling term, which then breaks down in course of the cosmological evolution. For our discussion, however, the $\mathcal{Z}_2$ symmetry of the potential is not mandatory the model to work. The conventional mechanism of particle production by the oscillating inflaton field (and their subsequent thermalization) remains operative. It is plain to see that the proposed construction can be easily applied for many successful models of inflation to incorporate dark energy at the expense of coupling with the CNB. We address the issue of neutrino nuggets from the quantum field theory point of view. Namely, these nuggets are considered as bound states caused basically by the Yukawa force, which arises in the framework of linear perturbation theory due to exchange of virtual quanta of quintessence field between the neutrinos.

The origin of cosmic rays in our Galaxy remains a subject of active debate. While supernova remnant shocks are often invoked as the sites of acceleration, it is now widely accepted that the difficulties of such sources in reaching PeV energies are daunting and it seems likely that only a subclass of rare remnants can satisfy the necessary conditions. Moreover the spectra of cosmic rays escaping the remnants have a complex shape that is not obviously the same as the spectra observed at the Earth. Here we investigate the process of particle acceleration at the termination shock that develops in the bubble excavated by star clusters' winds in the interstellar medium. While the main limitation to the maximum energy in supernova remnants comes from the need for effective wave excitation upstream so as to confine particles in the near-shock region and speed up the acceleration process, at the termination shock of star clusters the confinement of particles upstream in guaranteed by the geometry of the problem. We develop a theory of diffusive shock acceleration at such shock and we find that the maximum energy may reach the PeV region for powerful clusters in the high end of the luminosity tail for these sources. A crucial role in this problem is played by the dissipation of energy in the wind to magnetic perturbations. Under reasonable conditions the spectrum of the accelerated particles has a power law shape with a slope $4\div 4.3$, in agreement with what is required based upon standard models of cosmic ray transport in the Galaxy.

We report the detection of [O I]145.5um in the BR 1202-0725 system, a compact group at z=4.7 consisting of a quasar (QSO), a submillimeter-bright galaxy (SMG), and three faint Lya emitters. By taking into account the previous detections and upper limits, the [O I]/[C II] line ratios of the now five known high-z galaxies are higher than or on the high-end of the observed values in local galaxies ([O I]/[C II]$\gtrsim$0.13). The high [O I]/[C II] ratios and the joint analysis with the previous detection of [N II] lines for both the QSO and the SMG suggest the presence of warm and dense neutral gas in these highly star-forming galaxies. This is further supported by new CO (12-11) line detections and a comparison with cosmological simulations. There is a possible positive correlation between the [NII]122/205 line ratio and the [O I]/[C II] ratio when all local and high-z sources are taken into account, indicating that the denser the ionized gas, the denser and warmer the neutral gas (or vice versa). The detection of the [O I] line in the BR1202-0725 system with a relatively short amount of integration with ALMA demonstrates the great potential of this line as a dense gas tracer for high-z galaxies.

LOFAR is the only radio telescope that is presently capable of high-sensitivity, high-resolution (<1 mJy/b and <15") observations at ultra-low frequencies (<100 MHz). To utilise these capabilities, the LOFAR Surveys Key Science Project is undertaking a large survey to cover the entire northern sky with Low Band Antenna (LBA) observations. The LOFAR LBA Sky Survey (LoLSS) aims to cover the entire northern sky with 3170 pointings in the frequency range 42-66 MHz, at a resolution of 15 arcsec and at a sensitivity of 1 mJy/b. Here we outline the survey strategy, the observational status, the current calibration techniques, and briefly describe several scientific motivations. We also describe the preliminary public data release. The preliminary images were produced using a fully automated pipeline that aims to correct all direction-independent effects in the data. Whilst the direction-dependent effects, such as those from the ionosphere, are not yet corrected, the images presented in this work are still 10 times more sensitive than previous surveys available at these low frequencies. The preliminary data release covers 740 sqdeg around the HETDEX spring field region at a resolution of 47" with a median noise level of 5 mJy/b. The images and the catalogue with 25,247 sources are publicly released. We demonstrate that the system is capable of reaching an rms noise of 1 mJy/b and the resolution of 15" once direction-dependent effects are corrected for. LoLSS will provide the ultra-low-frequency information for hundreds of thousands of radio sources, providing critical spectral information and producing a unique dataset that can be used for a wide range of science topics such as: the search for high redshift galaxies and quasars, the study of the magnetosphere of exoplanets, and the detection of the oldest populations of cosmic-rays in galaxies, clusters of galaxies, and from AGN activity.

Growing evidence indicates that the synchrotron radiation mechanism may be responsible for the prompt emission of gamma-ray bursts (GRBs). In the synchrotron radiation scenario, the electron energy spectrum of the prompt emission is diverse in theoretical works and has not been estimated from observations in a general way (i.e., without specifying a certain physical model for the electron spectrum). In this paper, we creatively propose a method to directly estimate the electron spectrum for the prompt emission, without specifying a certain physical model for the electron spectrum in the synchrotron radiation scenario. In this method, an empirical function (i.e., a four-order Bezier curve jointed with a linear function at high-energy) is applied to describe the electron spectrum in log-log coordinate. It is found that our empirical function can well mimic the electron spectra obtained in many numerical calculations or simulations. Then, our method can figure out the electron spectrum for the prompt emission without specifying a model. By employing our method on observations, taking GRB 180720B and GRB 160509A as examples, it is found that the obtained electron spectra are generally different from that in the standard fast-cooling scenario and even a broken power law. Moreover, the morphology of electron spectra in its low-energy regime varies with time in a burst and even in a pulse. Our proposed method provides a valuable way to confront the synchrotron radiation mechanism with observations.

To predict line emission in the solar atmosphere requires models which are fundamentally different depending on whether the emission is from the chromosphere or the corona. At some point between the two regions, there must be a change between the two modelling regimes. Recent extensions to the coronal modelling for carbon and oxygen lines in the solar transition region have shown improvements in the emission of singly- and doubly-charged ions, along with Li-like ions. However, discrepancies still remain, particularly for singly-charged ions and intercombination lines. The aim of this work is to explore additional atomic processes that could further alter the charge state distribution and the level populations within ions, in order to resolve some of the discrepancies. To this end, excitation and ionisation caused by both the radiation field and by atom-ion collisions have been included, along with recombination through charge transfer. The modelling is carried out using conditions which would be present in the quiet Sun, which allows an assessment of the part atomic processes play in changing coronal modelling, separately from dynamic and transient events taking place in the plasma. The effect the processes have on the fractional ion populations are presented, as well as the change in level populations brought about by the new excitation mechanisms. Contribution functions of selected lines from low charge states are also shown, to demonstrate the extent to which line emission in the lower atmosphere could be affected by the new modelling.

We present a concept of a millimeter wavefront sensor that allows real-time sensing of the surface of a ground-based millimeter/submillimeter telescope. It is becoming important for ground-based millimeter/submillimeter astronomy to make telescopes larger with keeping their surface accurate. To establish `millimetric adaptive optics (MAO)' that instantaneously corrects the wavefront degradation induced by deformation of telescope optics, our wavefront sensor based on radio interferometry measures changes in excess path lengths from characteristic positions on the primary mirror surface to the focal plane. This plays a fundamental role in planned 50-m class submillimeter telescopes such as LST and AtLAST.

Using particle-in-cell (PIC) numerical simulations with electron-positron pair plasma, we study how the efficiencies of magnetic dissipation and particle acceleration scale with the initial coherence length $\lambda_0$ in relation to the system size $L$ of the two-dimensional (2D) `Arnold-Beltrami-Childress' (ABC) magnetic field configurations. Topological constraints on the distribution of magnetic helicity in 2D systems, identified earlier in relativistic force-free (FF) simulations, that prevent the high-$(L/\lambda_0)$ configurations from reaching the Taylor state, limit the magnetic dissipation efficiency to about $\epsilon_{\rm diss} \simeq 60\%$. We find that the peak growth time scale of the electric energy $\tau_{\rm E,peak}$ scales with the characteristic value of initial Alfven velocity $\beta_{\rm A,ini}$ like $\tau_{\rm E,peak} \propto (\lambda_0/L)\beta_{\rm A,ini}^{-3}$. The particle energy change is decomposed into non-thermal and thermal parts, with non-thermal energy gain dominant only for high initial magnetisation. The most robust description of the non-thermal high-energy part of the particle distribution is that the power-law index is a linear function of the initial magnetic energy fraction.

Solar cycle 24 is one of the weakest solar cycles recorded, but surprisingly the declining phase of it had a slow CME which evolved without any low coronal signature and is classified as a stealth CME which was responsible for an intense geomagnetic storm at Earth (Dst = -176 nT). The impact of this space weather event on the terrestrial ionosphere has been reported. However,the propagation of this CME beyond 1 au and the impact of this CME on other planetary environments have not been studied so far. In this paper, we analyse the data from Sun-Earth L1 point as well as from the Martian orbit (near 1.5 au) to understand the characteristics of the stealth CME as observed beyond 1 au. The observations near Earth are using data from the Solar Dynamics Observatory (SDO) and the Advanced Composition Explorer (ACE) satellite located at L1 point whereas those near Mars are from the instruments for plasma and magnetic field measurements on board Mars Atmosphere and Volatile EvolutioN (MAVEN) mission. The observations show that the stealth CME has reached 1.5 au after 7 days of its initial observations at the Sun and caused depletion in the nightside topside ionosphere of Mars, as observed during the inbound phase measurements of the Langmuir Probe and Waves (LPW) instrument on board MAVEN. These observations have implications on the ion escape rates from the Martian upper atmosphere.

We present an updated version of the semi-analytic Just-Jahrei{\ss} (JJ) model of the Galactic disk and constrain its parameters in the Solar neighbourhood. The new features of the JJ model include a simple two-component gaseous disk, a star-formation rate (SFR) function of the thick disk that has been extended in time, and a correlation between the kinematics of molecular gas and thin-disk populations. Here, we study the vertical number density profiles and W-velocity distributions determined from ~2 million local stars of the Gaia DR2. We also investigate an apparent Hess diagram of the Gaia stars selected in a conic volume towards the Galactic poles. Using a stellar evolution library, we synthesise stellar populations with a four-slope broken power-law initial mass function (IMF), the SFR, and an age-metallicity relation (AMR). The latter is consistently derived with the observed metallicity distribution of the local Red Clump (RC) giants from the APOGEE. Working within a Bayesian approach, we sample the posterior probability distribution in a multidimensional parameter space using the MCMC method. We find that the spatial distribution and motion of the Gaia stars imply two recent SF bursts centered at ages of ~0.5 Gyr and ~3 Gyr and characterised by a ~30% and ~55% SF enhancement, respectively, relative to a monotonously declining SFR continuum. The stellar populations associated with this SF excess are found to be dynamically hot for their age: they have W-velocity dispersions of ~12.5 km/s and ~26 km/s. The new JJ model is able to reproduce the local star counts with an accuracy of ~5 %. Using Gaia DR2 data, we self-consistently constrained 22 parameters of the updated JJ model. Our optimised model predicts two SF bursts within the last ~4 Gyr, which may point to recent episodes of gas infall.

We performed direct N-body simulations of open clusters with four different metallicities. To investigate the effective spins of merging binary black holes (BBHs) originated from these open clusters, we calculated the spin evolution of Wolf-Rayet (WR) stars with close companion stars (BBH progenitors), taking into account stellar wind mass loss and tidal spin-up of the WR stars. We found that BBH progenitors with smaller semi-major axes evolve to merging BBHs with greater effective spins because of strong tidal forces. In the local Universe, about 16% of merging BBHs get effective spins larger than 0.1 even if BHs and their progenitors do not get spin angular momenta other than tidal forces exerted by their companion stars. If we assume that WR stars have flat and isotropic distribution of dimensionless spins just after common envelope phases, the effective spin distribution of merging BBHs is similar to that inferred from gravitational wave observations with LIGO and Virgo.

The progenitor systems accounting for explosions of type Ia supernovae (SNe Ia) is still under debate. Symbiotic channel is one of the possible progenitor scenarios, in which the WDs in these systems increase in mass through wind accretion from their red giant companions. The mass-loss processes of the giants in the symbiotic systems could produce amount of circumstellar medium (CSM), and the detections of the interaction signals between SN ejecta and CSM can be used as an ideal way to distinguish the different progenitor models. However, the density distribution and geometric structure of the CSM around the symbiotic systems remain highly uncertain. By assuming that the tidal torque from binary interaction can increase the mass-loss rate of the red giant and cause the stellar wind concentrate towards the equatorial plane, we provide a simplified method to estimate the density and the degree of deviation from spherical symmetry of the CSM. Based on the calculations of the binary evolutions of symbiotic systems using stellar evolution code MESA, we obtained the parameter space for producing SNe Ia. We found that SNe Ia could originate from symbiotic systems with massive carbon-oxygen white dwarfs (CO WDs), while the half-opening angle of the stellar wind from red giant towards the WD varies with the binary evolution, resulting in the formation of surrounding CSM with different geometric structures. The corresponding properties of ejecta-CSM interactions may be examined by the spectropolarimetry observations in the future, from which one may find additional relationship between circumstellar environment of SNe Ia and their progenitor systems.

Uniform large transition-edge sensor (TES) arrays are fundamental for the next generation of X-ray space observatories. These arrays are required to achieve an energy resolution $\Delta E$ < 3 eV full-width-half-maximum (FWHM) in the soft X-ray energy range. We are currently developing X-ray microcalorimeter arrays for use in future laboratory and space-based X-ray astrophysics experiments and ground-based spectrometers. In this contribution we report on the development and the characterization of a uniform 32$\times$32 pixel array with 140$\times$30 $\mu$m$^2$ Ti/Au TESs with Au X-ray absorber. We report upon extensive measurements on 60 pixels in order to show the uniformity of our large TES array. The averaged critical temperature is $T_\mathrm{c}$ = 89.5$\pm$0.5 mK and the variation across the array ($\sim$1 cm) is less than 1.5 mK. We found a large region of detector's bias points between 20\% and 40\% of the normal-state resistance where the energy resolution is constantly lower than 3 eV. In particular, results show a summed X-ray spectral resolution $\Delta E_\mathrm{FWHM}$ = 2.50$\pm$0.04 eV at a photon energy of 5.9 keV, measured in a single-pixel mode using a frequency domain multiplexing (FDM) readout system developed at SRON/VTT at bias frequencies ranging from 1 to 5 MHz. Moreover we compare the logarithmic resistance sensitivity with respect to temperature and current ($\alpha$ and $\beta$ respectively) and their correlation with the detector's noise parameter $M$, showing an homogeneous behaviour for all the measured pixels in the array.

Generation of magnetic field during inflation can explain its presence over a wide range of scales in the Universe. In Ref.\cite{sharma2017}, we proposed a model to generate these fields during inflation. These fields have nonzero anisotropic stress which lead to the generation of a stochastic background of gravitational waves (GW) in the early universe. Here we show that for a scenario of magnetogenesis where reheating takes place around QCD epoch, this stochastic GW background lies in the $95\%$ confidence region of the GW signal probed by NANOGrav collaboration. This is the case when the generated electromagnetic field (EM) energy density is $3-10\%$ of the background energy density at the end of reheating. For this case, the values of magnetic field strength $B_0 \sim (3.8-6.9) \times 10^{-11}$G and its coherence length $\sim 30$ kpc at the present epoch. These values are for the models in which EM fields are of nonhelical nature. For the helical nature of the fields, these values are $B_0 \sim (1.1-1.9) \times 10^{-9}$G and its coherence length $\sim 0.8$ Mpc.

We present the theoretical framework and the numerical setting of JASMINE2, a code designed to efficiently solve the Jeans equations for multi-component axisymmetric stellar systems. The models may include an arbitrary number of stellar distributions, a dark matter halo, and a central supermassive black hole; each stellar distribution is implicitly described by a two-integral distribution function, and the stellar components can have different structural (density profile, flattening, mass,scale-length), dynamical (rotation, velocity dispersion anisotropy), and population (age, metallicity, initial mass function, mass-to-light ratio) properties. In order to determine the ordered rotational velocity and the azimuthal velocity dispersion fields of each component, we introduce a decomposition that can be used when the commonly adopted Satoh decomposition cannot be applied. The numerical implementation of JASMINE2, and of the post-processing procedures (including projection), are optimised to fully exploit the scalings allowed by the Poisson and the Jeans equations. For illustrative purposes, we present three multi-component galaxy models with a central black hole and a dark matter halo; one of the models is also used to test JASMINE2 against available analytical solutions.

M-dwarfs have proven to be ideal targets for planetary radial velocity (RV) searches due to their higher planet-star mass contrast. The HADES and CARMENES programs aim to carry out extensive searches of exoplanetary systems around this type of stars in the northern hemisphere, allowing us to address statistically the properties of the planets orbiting these objects. In this work, we perform a spectroscopic and photometric study of one of the program stars (GJ 740), which exhibits a short-period RV signal compatible with a planetary companion. We carried out a spectroscopic analysis based on 129 HARPS-N spectra taken over a time-span of 6 yr combined with 57 HARPS spectra taken over 4 yr, as well as 32 CARMENES spectra taken during more than 1 yr, resulting in a dataset with a time coverage of 10 yr. We also relied on 459 measurements from the public ASAS survey with a time-coverage of 8 yr along with 5 yr of photometric magnitudes from the EXORAP project taken in the $V$, $B$, $R$, and $I$ filters to carry out a photometric study. Both analyses were made using Markov Chain Monte Carlo (MCMC) simulations and Gaussian Process regression to model the activity of the star. We present the discovery of a short-period super-Earth with an orbital period of 2.37756$^{+0.00013}_{-0.00011}$ d and a minimum mass of 2.96$^{+0.50}_{-0.48}$ M$_{\oplus}$. We offer an update to the previously reported characterization of the magnetic cycle and rotation period of the star, obtaining values of $P_{\rm rot}$=35.563$\pm$0.071 d and $P_{\rm cycle}$=2800$\pm$150 d. Furthermore, the RV time-series exhibits a possibly periodic long-term signal which might be related to a Saturn-mass planet of $\sim$ 100 M$_{\oplus}$.

We present spectroscopy of the P~Cygni profile of the 1.083-micron He I line in the WC7 + O5 colliding-wind binary (CWB) WR 140 (HD 193793), observed in 2008, before its periastron passage in 2009, and in 2016-17, spanning the subsequent periastron passage. Both absorption and emission components showed strong variations. The variation of the absorption component as the O5 star was occulted by the wind-collision region (WCR) sets a tight constraint on its geometry. While the sightline to the O5 star traversed the WCR, the strength and breadth of the absorption component varied significantly on time-scales of days. An emission sub-peak was observed on all our profiles. The variation of its radial velocity with orbital phase was shown to be consistent with formation in the WCR as it swung round the stars in their orbit. Modelling the profile gave a measure of the extent of the sub-peak forming region. In the phase range 0.93-0.99, the flux in the sub-peak increased steadily, approximately inversely proportionally to the stellar separation, indicating that the shocked gas in the WCR where the line was formed was adiabatic. After periastron, the sub-peak flux was anomalously strong and varied rapidly, suggesting formation in clumps down-stream in the WCR. For most of the time, its flux exceeded the 2-10-keV X-ray emission, showing it to be a significant coolant of the shocked wind.

Various properties of Jovian trojan asteroids such as composition, rotation periods, and photometric amplitudes, or the rate of binarity in the population can provide information and constraints on the evolution of the group and of the Solar System itself. Here we present new photometric properties of 45 Jovian trojans from the K2 mission of the Kepler space telescope, and present phase-folded light curves for 44 targets, including (11351) Leucus, one of the targets of the Lucy mission. We extend our sample to 101 asteroids with previous K2 Trojan measurements, then compare their combined amplitude- and frequency distributions to other ground-based and space data. We show that there is a dichotomy in the periods of Trojans with a separation at $\sim 100$ hr. We find that 25% of the sample are slow rotators (P$\geq$30 hr), which excess can be attributed to binary objects. We also show that 32 systems can be classified as potential detached binary systems. Finally, we calculate density and rotation constraints for the asteroids. Both the spin barrier and fits to strengthless ellipsoid models indicate low densities and thus compositions similar to cometary and TNO populations throughout the sample. This supports the scenario of outer Solar System origin for Jovian trojans.

How the solar electromagnetic energy entering the Earth's atmosphere varied since pre-industrial times is an important consideration in the climate change debate. Detrimental to this debate, estimates of the change in total solar irradiance (TSI) since the Maunder minimum, an extended period of weak solar activity preceding the industrial revolution, differ markedly, ranging from a drop of 0.75 Wm-2 to a rise of 6.3 Wm-2. Consequently, the exact contribution by solar forcing to the rise in global temperatures over the past centuries remains inconclusive. Adopting a novel approach based on state-of-the-art solar imagery and numerical simulations, we establish the TSI level of the Sun when it is in its least-active state to be 2.0 +/- 0.7 Wm-2 below the 2019 level. This means TSI could not have risen since the Maunder minimum by more than this amount, thus restricting the possible role of solar forcing in global warming.

The WALOP-South instrument will be mounted on the 1 m SAAO telescope in South Africa as part of the PASIPHAE program to carry out a linear imaging polarization survey of the Galactic polar regions in the optical band. Designed to achieve polarimetric sensitivity of $0.05~\%$ across a $35\times35$ arcminute field of view, it will be capable of measuring the Stokes parameters I, q and u in a single exposure in the SDSS-r broadband and narrowband filters between $0.5~{\mu}m - 0.7~{\mu}m$. For each measurement, four images of the full field corresponding to linear polarization angles of 0 deg, 45 deg, 90 deg and 135 deg in the instrument coordinate system will be created on four detectors from which the Stokes parameters can be found using differential photometry. In designing the optical system, major challenges included correcting for the dispersion introduced by large split angle Wollaston Prisms used as analysers as well as other aberrations from the entire field to obtain imaging quality PSF at the detector. We present the optical design of the WALOP-South instrument which overcomes these challenges and delivers near seeing limited PSFs for the entire field of view.

Recently a new kind of quintessential inflation coming from the Lorentzian distribution has been introduced in [1,2]. The model leads to a very simple potential, which basically depends on two parameters, belonging to the class of $\alpha$-attractors and depicting correctly the early and late time accelerations of our universe. The potential emphasizes a {\it cosmological see-saw mechanism} (CSSM) that produces a large inflationary vacuum energy in one side of the potential and a very small value of dark energy on the right hand side of the potential. {Here we show that the model agrees with the recent observations and with the reheating constraints. Therefore the model gives a reasonable scenario beyond the standard $\Lambda$CDM that includes the inflationary epoch.}

WALOP (Wide-Area Linear Optical Polarimeter)-South, to be mounted on the 1m SAAO telescope in South Africa, is first of the two WALOP instruments currently under development for carrying out the PASIPHAE survey. Scheduled for commissioning in the year 2021, the WALOP instruments will be used to measure the linear polarization of around $10^{6}$ stars in the SDSS-r broadband with $0.1~\%$ polarimetric accuracy, covering 4000 square degrees in the Galactic polar regions. The combined capabilities of one-shot linear polarimetry, high polarimetric accuracy ($< 0.1~\%$) and polarimetric sensitivity ($< 0.05~\%$), and a large field of view (FOV) of $35\times35~arcminutes$ make WALOP-South a unique astronomical instrument. In a single exposure, it is designed to measure the Stokes parameters $I$, $q$ and $u$ in the SDSS-r broadband and narrowband filters between $500-700~nm$. During each measurement, four images of the full field corresponding to the polarization angles of $0^{\circ}$, $45^{\circ}$, $90^{\circ}$ and $135^{\circ}$ will be imaged on four detectors and carrying out differential photometry on these images will yield the Stokes parameters. Major challenges in designing WALOP-South instrument include- (a) in the optical design, correcting for the spectral dispersion introduced by large split angle Wollaston Prisms used as polarization analyzers as well as aberrations from the wide field, and (b) making an optomechanical design adherent to the tolerances required to obtain good imaging and polarimetric performance under all temperature conditions as well as telescope pointing positions. We present the optical and optomechanical design for WALOP-South which overcomes these challenges.

We aim to gain insight into the effect of network and faculae on solar irradiance from their apparent intensity. Taking full-disc observations from the Solar Dynamics Observatory, we examined the intensity contrast of network and faculae in the continuum and core of the Fe I 6173 {\AA} line and 1700 {\AA}, including the variation with magnetic flux density, distance from disc centre, nearby magnetic fields, and time. The brightness of network and faculae is believed to be suppressed by nearby magnetic fields from its effect on convection. The difference in intensity contrast between the quiet-Sun network and active region faculae, noted by various studies, arises because active regions are more magnetically crowded and is not due to any fundamental physical differences between network and faculae. These results highlight that solar irradiance models need to include the effect of nearby magnetic fields on network and faculae brightness. We found evidence that suggests that departures from local thermal equilibrium (LTE) might have limited effect on intensity contrast. This could explain why solar irradiance models that are based on the intensity contrast of solar surface magnetic features calculated assuming LTE reproduce the observed spectral variability even where the LTE assumption breaks down. Certain models of solar irradiance employ chromospheric indices as direct indications of the effect of network and faculae on solar irradiance. Based on past studies of the Ca II K line and on the intensity contrast measurements derived here, we show that the fluctuations in chromospheric emission from network and faculae are a reasonable estimate of the emission fluctuations in the middle photosphere, but not of those in the lower photosphere. The data set, which extends from 2010 to 2018, indicates that intensity contrast was stable to about 3% in this period.

We investigate the effect of the Biermann battery during the Epoch of Reionization (EoR) using cosmological Adaptive Mesh Refinement simulations within the framework of the SPHINX project. We develop a novel numerical technique to solve for the Biermann battery term in the Constrained Transport method, preserving both the zero divergence of the magnetic field and the absence of Biermann battery for isothermal flows. The structure-preserving nature of our numerical method turns out to be very important to minimise numerical errors during validation tests of the propagation of a Str\"omgren sphere and of a Sedov blast wave. We then use this new method to model the evolution of a 2.5 and 5 co-moving Mpc cosmological box with a state-of-the-art galaxy formation model within the RAMSES code. Contrary to previous findings, we show that three different Biermann battery channels emerge: the first one is associated with linear perturbations before the EoR, the second one is the classical Biermann battery associated with reionization fronts during the EoR, and the third one is associated with strong, supernova-driven outflows. While the two former channels generate spontaneously volume-filling magnetic fields with a strength on the order or below $10^{-20}$ G, the latter, owing to the higher plasma temperature and a marginally-resolved turbulent dynamo, reaches a field strength as high as $10^{-18}$ G in the intergalactic medium around massive haloes.

Summary statistics of the likelihood, such as the Bayesian evidence, offer a principled way of comparing models and assessing tension between, or within, the results of physical experiments. Noisy realisations of the data induce scatter in these model comparison statistics. For a realistic case of cosmological inference from large-scale structure we show that the logarithm of the Bayes factor attains scatter of order unity, increasing significantly with stronger tension between the models under comparison. We develop an approximate procedure that quantifies the sampling distribution of the evidence at small additional computational cost and apply it to real data to demonstrate the impact of the scatter, which acts to reduce the significance of any model discrepancies. Data compression is highlighted as a potential avenue to suppressing noise in the evidence to negligible levels, with a proof of concept on Planck cosmic microwave background data.

The axion provides a solution for the strong CP problem and is one of the leading candidates for dark matter. This paper proposes an axion detection scheme based on quantum nondemolition detection of magnon, i.e., quanta of collective spin excitations in solid, which is expected to be excited by the axion-electron interaction predicted by the Dine-Fischer-Srednicki-Zhitnitsky (DFSZ) model. The prototype detector is composed of a ferromagnetic sphere as an electronic spin target and a superconducting qubit. Both of these are embedded inside a microwave cavity, which leads to a coherent effective interaction between the uniform magnetostatic mode in the ferromagnetic crystal and the qubit. An upper limit for the coupling constant between an axion and an electron is obtained as $g_{aee}<1.6\times10^{-7}$ at the 95\% confidence level for the axion mass of $33.117$$\mu$eV$<m_{a}<33.130$$\mu$eV.

We reanalyse the solar eclipse linked to the Biblical passage about the military leader Joshua who ordered the sun to halt in the midst of the day (Joshua 10:12). Although there is agreement that the basic story is rooted in a real event, the date is subject to different opinions. We review the historical emergence of the text and confirm that the total eclipse of the sun of 30 September 1131 BCE is the most likely candidate. The Besselian Elements for this eclipse are re-computed. The error for the deceleration parameter of Earth's rotation, $\Delta T$, is improved by a factor of 2.