Papers for Thursday, Apr 22 2021

We study the Fundamental Plane (FP) for a volume- and luminosity-limited sample of 560 early-type galaxies from the SAMI survey. Using r-band sizes and luminosities from new Multi-Gaussian Expansion (MGE) photometric measurements, and treating luminosity as the dependent variable, the FP has coefficients a=1.294$\pm$0.039, b= 0.912$\pm$0.025, and zero-point c= 7.067$\pm$0.078. We leverage the high signal-to-noise of SAMI integral field spectroscopy, to determine how structural and stellar-population observables affect the scatter about the FP. The FP residuals correlate most strongly (8$\sigma$ significance) with luminosity-weighted simple-stellar-population (SSP) age. In contrast, the structural observables surface mass density, rotation-to-dispersion ratio, S\'ersic index and projected shape all show little or no significant correlation. We connect the FP residuals to the empirical relation between age (or stellar mass-to-light ratio $\Upsilon_\star$) and surface mass density, the best predictor of SSP age amongst parameters based on FP observables. We show that the FP residuals (anti-)correlate with the residuals of the relation between surface density and $\Upsilon_\star$. This correlation implies that part of the FP scatter is due to the broad age and $\Upsilon_\star$ distribution at any given surface mass density. Using virial mass and $\Upsilon_\star$ we construct a simulated FP and compare it to the observed FP. We find that, while the empirical relations between observed stellar population relations and FP observables are responsible for most (75%) of the FP scatter, on their own they do not explain the observed tilt of the FP away from the virial plane.

Orbit superposition models are a non-parametric dynamical modelling technique to determine the mass of a galaxy's central supermassive black hole (SMBH), its stars, or its dark-matter halo. One of the main problems is how to decide which model out of a large pool of trial models based on different assumed mass distributions represents the true structure of an observed galaxy best. We show that the traditional approach to judge models solely by their goodness-of-fit can lead to substantial biases in estimated galaxy properties caused by varying model flexibilities. We demonstrate how the flexibility of the models can be estimated using bootstrap iterations and present a model selection framework that removes these biases by taking the variable flexibility into account in the model evaluation. We extend the model selection approach to optimize the degree of regularisation directly from the data. Altogether, this leads to a significant improvement of the constraining power of the modeling technique. We show with simulations that one can reconstruct the mass, anisotropy and viewing angle of an axisymmetric galaxy with a few percent accuracy from realistic observational data with fully resolved line-of-sight velocity distributions (LOSVDs). In a first application, we reproduce a photometric estimate of the inclination of the disk galaxy NGC 3368 to within 5 degree accuracy from kinematic data that cover only a few sphere-of-influence radii around the galaxy's SMBH. This demonstrates the constraining power that can be achieved with orbit models based on fully resolved LOSVDs and a model selection framework.

An error in the gravitational force that the source of gravity induces on itself (a self-force error) violates both the conservation of linear momentum and the conservation of energy. If such errors are present in a self-gravitating system and are not sufficiently random to average out, the obtained numerical solution will become progressively more unphysical with time: the system will acquire or lose momentum and energy due to numerical effects. In this paper, we demonstrate how self-force errors can arise in the case where self-gravity is solved on an adaptively refined mesh when the refinement is nonuniform. We provide the analytical expression for the self-force error and numerical examples that demonstrate such self-force errors in idealized settings. We also show how these errors can be corrected to an arbitrary order by straightforward addition of correction terms at the refinement boundaries.

We show that measuring the velocity components of hypervelocity stars (HVSs) can discriminate between Modified Newtonian Dynamics (MOND) and Newtonian gravity. HVSs are ejected from the Galactic center on radial trajectories with a null tangential velocity component in the reference frame of the Galaxy. They acquire tangential components due to the non-spherical components of the Galactic gravitational potential. Axisymmetric potentials only affect the latitudinal components $v_\theta$, and non-null azimuthal components $v_\phi$ originate from non-axisymmetric matter distributions. For HVSs with sufficiently large ejection speed, the azimuthal velocity components are proportionate to the deviation of the gravitational potential from axial symmetry. The ejection velocity threshold is $\sim$ 750 km s$^{-1}$ for 4 $M_{\odot}$ stars and increases with decreasing HVS mass. We determine the upper limit of $v_\phi$ as a function of the Galactocentric distance for these high-speed HVSs, if MOND, in its QUMOND formulation, is the correct theory of gravity and the triaxial Galactic bulge is the primary source of the azimuthal component $v_\phi$. In Newtonian gravity, values of $v_\phi$ larger than this limit can easily appear if the dark matter halo is triaxial or if the dark matter halo and the baryonic components are axisymmetric but their two axes of symmetry are misaligned. Therefore, even a limited sample of high-speed HVSs could in principle distinguish between the QUMOND scenario and the dark matter model if the azimuthal components $v_\phi$ of the observed HVSs do not meet the QUMOND upper limit. However, the uncertainties on the Galactocentric azimuthal velocity components of the currently known HVSs are still at least a factor of $\sim 10$ too large to make this test conclusive. Astrometric measurements with micro-arcsecond precision would make this test feasible.

Ambipolar diffusion likely plays a pivotal role in the formation and evolution of dense cores in weakly-ionized molecular clouds. Linear analyses show that the evolutionary times and fragmentation scales are significantly greater than the hydrodynamic (Jeans) values even for clouds with mildly supercritical mass-to-flux ratio. We utilize values of fragmentation scales and growth times that correspond to typical ionization fractions within a molecular cloud, and apply to the context of the observed estimated lifetime of prestellar cores as well as the observed number of such embedded cores forming in a parent clump. By varying a single parameter, the mass-to-flux ratio, over the range of observationally measured densities, we fit the range of estimated prestellar core lifetimes ($\sim 0.1$ to a few Myr) identified with Herschel as well as the number of embedded cores formed in a parent clump measured in Perseus with the Submillimeter Array (SMA). Our model suggests that the prestellar cores are formed with a transcritical mass-to-flux ratio and higher densities correspond to somewhat higher mass-to-flux ratio but the normalized mass-to-flux ratio $\mu$ remains in the range $1 \lesssim \mu \lesssim 2$. Our best-fit model exhibits $B \propto n^{0.43}$ for prestellar cores, due to partial flux-freezing as a consequence of ambipolar diffusion.

We propose a method of searching for radio sources exhibiting the Sunyaev-Zeldovich effect in the multi-frequency emission maps from the Planck mission data using a convolutional neural network. A catalog for recognizing radio sources is compiled using the GLESP pixelation scheme at the frequencies of 100, 143, 217, 353, and 545 GHz. The quality of the proposed approach is evaluated and the quality of the dependence of model data on the S/N ratio is estimated. We show that the presented neural network approach allows the detection of sources with the Sunyaev-Zeldovich effect. The proposed method can be used to find the most likely galaxy cluster candidates at large redshifts.

Non-linear effects in the early Universe generate non-zero bispectra of the cosmic microwave background (CMB) temperature and polarization, even in the absence of primordial non-Gaussianity. In this paper, we compute the contributions from isocurvature modes to the CMB bispectra using a modified version of the second-order Boltzmann solver SONG. We investigate the ability of current and future CMB experiments to constrain these modes with observations of the bispectrum. Our results show that the enhancement due to single isocurvature modes mixed with the adiabatic mode is negligible for the parameter ranges currently allowed by the most recent Planck results. However, we find that a large compensated isocurvature mode can produce a detectable bispectrum when its correlation with the adiabatic mode is appreciable. The non-observation of this contribution in searches for the lensing bispectrum from Planck allows us to place a new constraint on the relative amplitude of the correlated part of the compensated isocurvature mode of $f_{\rm CIP}=1\pm100$. We compute forecasts for future observations by COrE, SO, CMB-S4 and an ideal experiment and conclude that a dedicated search for the bispectrum from compensated modes could rule out a number of scenarios realised in the curvaton model. In addition, the CMB-S4 experiment could detect the most extreme of those scenarios ($f_{\rm CIP}=16.5$) at 2 to 3-$\sigma$ significance.

Fast outflows and their interaction with slow shells (generally known as the fossil circumstellar envelope of asymptotic giant branch stars) play an important role in the structure and kinematics of protoplanetary and planetary nebulae (pPNe, PNe). To properly study their effects within these objects, we also need to observe the intermediate-temperature gas, which is only detectable in the far-infrared (FIR) and submillimetre (submm) transitions. We study the physical conditions of the outflows presented in a number of pPNe and PNe, with a focus on their temperature and excitation states. We carried out Herschel/HIFI observations in the submm lines of 12CO in nine pPNe and nine PNe and complemented them with low-J CO spectra obtained with the IRAM 30m telescope and taken from the literature. The spectral resolution of HIFI allows us to identify and measure the different nebular components in the line profiles. The comparison with large velocity gradient (LVG) model predictions was used to estimate the physical conditions of the warm gas in the nebulae, such as excitation conditions, temperature, and density. We found high kinetic temperatures for the fast winds of pPNe, typically reaching between 75 K and 200 K. In contrast, the high-velocity gas in the sampled PNe is colder, with characteristic temperatures between 25 K and 75 K, and it is found in a lower excitation state. We interpret this correlation of the kinetic temperature and excitation state of fast outflows with the amount of time elapsed since their acceleration (probably driven by shocks) as a consequence of the cooling that occurred during the pPN phase.

We present the first results of the Galaxy Activity, Torus and Outflows Survey (GATOS), a project aimed at understanding the properties of the dusty molecular tori and their connection to the host galaxy in nearby Seyfert galaxies. Our project expands the range of AGN luminosities and Eddington ratios covered by previous surveys of Seyferts conducted by ALMA and allows us to study the gas feeding and feedback cycle in a combined sample of 19 Seyferts. We used ALMA to obtain new images of the emission of molecular gas and dust using the CO(3-2) and HCO+(4-3) lines as well as their underlying continuum emission at 870 microns with high spatial resolutions (0.1'' ~ 7 - 13 pc) in the CND of 10 nearby (D < 28 Mpc) Seyfert galaxies. Our new ALMA observations detect 870 micron continuum and CO line emission from spatially resolved disks located around the AGN in all the sources. The bulk of the continuum flux can be accounted for by thermal emission from dust in the majority of the targets. For most of the sources the disks show a preponderant orientation perpendicular to the AGN wind axes, as expected for dusty molecular tori. The median diameters and molecular gas masses of the tori are ~ 42 pc, and ~ 6 x 10**5 Msun, respectively. We find a positive correlation between the line-of-sight gas column densities responsible for the absorption of X-rays and the molecular gas column densities derived from CO towards the AGN in our sources. The radial distributions of molecular gas in the CND of our combined sample show signs of nuclear-scale molecular gas deficits. We also detect molecular outflows in the sources that show the most extreme nuclear-scale gas deficits in our sample. These observations find for the first time supporting evidence that the imprint of AGN feedback is more extreme in higher luminosity and/or higher Eddington ratio Seyfert galaxies.

The origin of the slow solar wind is still an open issue. It has been suggested that upflows at the edge of active regions (AR) can contribute to the slow solar wind. Here, we compared the upflow region and the AR core and studied how the plasma properties change from the chromosphere via the transition region to the corona. We studied limb-to-limb observations NOAA 12687 (14th - 25th Nov 2017). We analysed spectroscopic data simultaneously obtained from IRIS and Hinode/EIS in six spectral lines. We studied the mutual relationships between the plasma properties for each emission line, as well as comparing the plasma properties between the neighbouring formation temperature lines. To find the most characteristic spectra, we classified the spectra in each wavelength using the machine learning technique k-means. We found that in the upflow region the Doppler velocities of the coronal lines are strongly correlated, but the transition region and coronal lines show no correlation. However, their fluxes are strongly correlated. The upflow region has lower density and lower temperature than the AR core. In the upflow region, the Doppler and non-thermal velocity show a strong correlation in the coronal lines, but the correlation is not seen in the AR core. At the boundary between the upflow region and the AR core, the upflow region shows an increase in the coronal non-thermal velocity, the emission obtained from the DEM, and the domination of the redshifted regions in the chromosphere. The obtained results suggest that at least three parallel mechanisms generate the plasma upflow: (1) the reconnection between closed loops and open magnetic field lines in the lower corona or upper chromosphere; (2) the reconnection between the chromospheric small-scale loops and open magnetic field; (3) the expansion of the magnetic field lines that allows the chromospheric plasma to escape to the solar corona.

Using a generalized function of the stellar spin-down law, we investigate the age dependence of the magnetic braking index ($q$). Our survey includes 9 open clusters aged lower than 1 Gyr and ranged in mass from 0.7 to 1.1$M_{\odot}$. Our aim is to verify the time behavior of the nonextensive braking index $q$ which brings the cumulative distribution of the rotational velocities of the stars of the youngest cluster ($\alpha{\rm Per}$) taken at the future age of an older cluster. As a result, the $q$-index is calculated over time $t-t_{\alpha{\rm Per}}$, where $t$ is the age of older open cluster used to estimate the future cumulative distribution of the rotational velocity of the $\alpha{\rm Per}$ cluster with present-day age $t_{\alpha{\rm Per}}$. We find that the values of $q$ are slightly constant around 1.36 and 1.38 according to the mass bin. In conclusion, the results seem to indicate that the mechanism that controls the rotational decay of stars in open clusters does not depend on the increment of time.

We revisit Wagoner et al. (1967), a classic contribution in the development of Big Bang Nucleosynthesis. We demonstrate that it presents an incorrect expression for the temperature of the early universe as a function of time in the high temperature limit, $T \gtrsim 10^{10}$K. As this incorrect expression has been reproduced elsewhere, we present a corrected form for the initial conditions required for calculating the formation of the primordial elements in the Big Bang.

Observations conducted over the last few decades show that the expansion of the Universe is accelerating. In the standard model of cosmology, this accelerated expansion is attributed to a dark energy in the form of a cosmological constant. It is conceivable, however, for the dark energy to exhibit mild dynamics (so that its energy density changes with time rather than having a constant value), or for the accelerated expansion of the Universe to be caused by some mechanism other than dark energy. In this work I will investigate both of these possibilities by using observational data to place constraints on the parameters of simple models of dynamical dark energy as well as cosmological models without dark energy. I find that these data favor the standard model while leaving some room for dynamical dark energy. The standard model also holds that the Universe is flat on large spatial scales. The same observational data used to test dark energy dynamics can be used to constrain the large-scale curvature of the Universe, and these data generally favor spatial flatness, with some mild preference for spatial curvature in some data combinations.

Carbon monoxide (CO) is the most abundant molecule after molecular hydrogen and is important for the chemistry in circumstellar envelopes around evolved stars. When modelling the strength and shape of molecular lines, the size of the CO envelope is an input parameter and influences the derived mass-loss rates. In particular the low-J transition CO lines are sensitive to the CO photodissociation radius. Recently, new CO photodissociation radii have been published using different formalisms that differ considerably. One set of calculations is based on an escape-probability formalisms that uses numerical approximations derived in the early-eighties. The accuracy of these approximations is investigated and it is shown that they are less accurate than claimed. Improved formalism are derived. Nevertheless, the changes in CO envelope size are small to moderate, less than 2\% for models with $10^{-7}< \dot{M}< 10^{-4}$ \msolyr\ and at most 7\% for model with $\dot{M} = 10^{-8}$ \msolyr.

The nuclear stellar disc (NSD) is, together with the nuclear star cluster (NSC) and the central massive black hole, one of the main components in the central parts of our Milky Way. However, until recently, only few studies of the stellar content of the NSD have been obtained due to extreme extinction and stellar crowding. With a dedicated KMOS (VLT, ESO) spectroscopic survey, we study the kinematics and global metallicities of the NSD based on the observations of K/M giant stars. We trace radial velocities and metallicities which were derived based on spectral indices (Na I and CO) along the NSD and compare those with a Galactic Bulge sample of APOGEE (DR16) and data from the NSC. We find that the metallicity distribution function and the fraction of metal-rich and metal-poor stars in the NSD are different from the corresponding distributions and ratios of the NSC and the Galactic Bulge. By tracing the velocity dispersion as a function of metallicity, we clearly see that the NSD is kinematically cool and that the velocity dispersion decreases with increasing metallicity contrary to the inner Bulge sample of APOGEE ($\rm |b| < 4^{o}$). Using molecular gas tracers ($\rm H_{2}CO$, CO(4-3)) of the Central Molecular Zone (CMZ) we find an astonishing agreement between the gas rotation and the rotation of the metal-rich population indicating that the metal-rich stars could have formed from gas in the CMZ. On the other hand, the metal-poor stars show a much slower rotation profile with signs of counter-rotation indicating a different origin of these stars. Coupling kinematics with global metallicities, our results demonstrate that the NSD is chemically and kinematically distinct with respect to the inner Bulge indicating a different formation scenario.

Gamma-ray Bursts (GRBs) are the most powerful transients in the Universe, over-shining for a few seconds all other $\gamma$-ray sky sources. Their emission is produced within narrowly collimated relativistic jets launched after the core-collapse of massive stars or the merger of compact binaries. THESEUS will open a new window for the use of GRBs as cosmological tools by securing a statistically significant sample of high-$z$ GRBs, as well as by providing a large number of GRBs at low-intermediate redshifts extending the current samples to low luminosities. The wide energy band and unprecedented sensitivity of the Soft X-ray Imager (SXI) and X-Gamma rays Imaging Spectrometer (XGIS) instruments provide us a new route to unveil the nature of the prompt emission. For the first time, a full characterisation of the prompt emission spectrum from 0.3 keV to 10 MeV with unprecedented large count statistics will be possible revealing the signatures of synchrotron emission. SXI spectra, extending down to 0.3 keV, will constrain the local metal absorption and, for the brightest events, the progenitors' ejecta composition. Investigation of the nature of the internal energy dissipation mechanisms will be obtained through the systematic study with XGIS of the sub-second variability unexplored so far over such a wide energy range. THESEUS will follow the spectral evolution of the prompt emission down to the soft X-ray band during the early steep decay and through the plateau phase with the unique ability of extending above 10 keV the spectral study of these early afterglow emission phases.

Context. Hot subdwarfs experienced strong mass loss on the Red Giant Branch (RGB) and are now hot and small He-burning objects. Aims. In this project we aim to perform a transit survey in all available light curves of hot subdwarfs from space-based telescopes (Kepler, K2, TESS, and CHEOPS), with our custom-made pipeline SHERLOCK, in order to determine the occurrence rate of planets around these stars, as a function of orbital period and planetary radius. Methods. In this first paper, we perform injection-and-recovery tests of synthetic transits for a selection of representative Kepler, K2 and TESS light curves, to determine which transiting bodies, in terms of object radius and orbital period, we will be able to detect with our tools. We also provide such estimates for CHEOPS data, which we analyze with the pycheops package. Results. Transiting objects with a radius $\lesssim$ 1.0 $R_{\Earth}$ can be detected in most of Kepler, K2 and CHEOPS targets for the shortest orbital periods (1 d and below), reaching values as small as $\sim$0.3 $R_{\Earth}$ in the best cases. Reaching sub-Earth-sized bodies is achieved only for the brightest TESS targets, and the ones observed during a significant number of sectors. We also give a series of representative results for farther and bigger planets, for which the performances strongly depend on the target magnitude, the length and the quality of the data. Conclusions. The TESS sample will provide the most important statistics for the global aim of measuring the planet occurrence rate around hot subdwarfs. The Kepler, K2 and CHEOPS data will allow us to search for planetary remnants, i.e. very close and small (possibly disintegrating) objects, which would have partly survived the engulfment in their red giant host.

Context. The nitrogen reservoir in planetary systems is a long standing problem. Part of the N-bearing molecules is probably incorporated into the ice bulk during the cold phases of the stellar evolution, and may be gradually released into the gas phase when the ice is heated, such as in active comets. The chemical nature of the N-reservoir should greatly influence how, when and in what form N returns to the gas phase, or is incorporated into the refractory material forming planetary bodies. Aims. We present the study the thermal desorption of two ammonium salts: ammonium formate and ammonium acetate from a gold surface and from a water ice substrate. Methods. Temperature-programmed desorption experiments and Fourier transform infrared reflection spectroscopy were conducted to investigate the desorption behavior of ammonium salts. Results. Ammonium salts are semi-volatile species releasing neutral species as major components upon desorption, that is ammonia and the corresponding organic acid (HCOOH and CH3COOH), at temperatures higher than the temperature of thermal desorption of water ice. Their desorption follows a first-order Wigner-Polanyi law. We find the first order kinetic parameters A = 7.7 $\pm$ 0.6 $\times$ 10$^{15}$ s$^{-1}$ and E$_{bind}$ = 68.9 $\pm$ 0.1 kJ~mol$^{-1}$ for ammonium formate and A = 3.0 $\pm$ 0.4 $\times$ 10$^{20}$ s$^{-1}$ and E$_{bind}$ = 83.0 $\pm$ 0.2 kJ~mol$^{-1}$ for ammonium acetate. The presence of a water ice substrate does not influence the desorption kinetics. Ammonia molecules locked in salts desorb as neutral molecules at temperatures much higher than previously expected that are usually attributed to refractory materials. Conclusions. Ammonia snow-line has a smaller radius than the water snow-line. As a result, the NH3/H2O ratio content in solar system bodies can be a hint as to where they formed and subsequently migrated.

The hydrogen cyanide (HCN) molecule in the planetary atmosphere is key to the formation of building blocks of life. We present the spectroscopic detection of the rotational molecular line of nitrile species hydrogen cyanide (HCN) in the atmosphere of Saturn using the archival data of the Atacama Large Millimeter/Submillimeter Array (ALMA) in band 7 observation. The strong rotational emission line of HCN is detected at frequency $\nu$ = 354.505 GHz (>4$\sigma$ statistical significance). We also detect the rotational emission line of carbon monoxide (CO) at frequency $\nu$ = 345.795 GHz. The statistical column density of hydrogen cyanide and carbon monoxide emission line is N(HCN)$\sim$2.42$\times$10$^{16}$ cm$^{-2}$ and N(CO)$\sim$5.82$\times$10$^{17}$ cm$^{-2}$. The abundance of HCN and CO in the atmosphere of Saturn relative to the H$_{2}$ is estimated to be f(HCN)$\sim$1.02$\times$10$^{-9}$ and f(CO)$\sim$2.42$\times$10$^{-8}$. We discussed possible chemical pathways to the formation of the detected nitrile gas HCN in the atmosphere of Saturn.

Galactic short orbital period black hole candidate (BHC) XTE~J1752-223 was discovered on 2009 Oct 21 by the Rossi X-ray Timing Explorer (RXTE). We study the spectral properties of this outburst using transonic flow solution based two component advective flow (TCAF) model. TCAF model fitted spectrum gives an estimation of the physical flow parameters, such as the Keplerian disk rate, sub-Keplerian halo rate, properties of the so-called {\it{Compton cloud}}, other than the mass of the source and normalization ($N$). $N$ is a standardized ratio of emitted to observed photon flux in TCAF which does not include X-ray emission from jets. In the presence of jets, this ratio changes and this deviation is used to obtain the estimation of X-ray contribution from the jets. Nature of the jet is found to be compact during low luminous hard state and discrete or blobby during high luminous intermediate states. We find a correlation between the radio (5.5 GHz) and X-ray ($2.5-25$ keV) fluxes from different components. The radio ($F_R$) and jet X-ray ($F_{ouf}$) fluxes are found to be correlated within the acceptable range of the standard correlation ($0.6$ to $0.7$). A similar correlation indices were reported by our group for three other short orbital period transient BHCs (Swift~J1753.5-0127, MAXI~J1836-194 \& XTE~J1118+480).

We have studied stellar candidates for close (within 1 pc) encounters with the Solar system. For all of the stars under consideration the kinematic characteristics have been taken from the Gaia EDR3 catalogue. The parameters of the encounters of these stars with the Solar system have been calculated using three methods: (1) the linear one, (2) by integrating the orbits in an axisymmetric potential, and (3) by integrating the orbits in a potential with a spiral density wave. All three methods are shown to yield similar results. We have selected five stars that are good candidates for reaching the boundaries of the Oort cloud and passing through it. Based on the second method, in good agreement with the other two methods, we have obtained the following estimates of the encounter parameters for the star GJ 710: $t_{min}=1.320\pm0.028$ Myr and $d_{min}=0.020\pm0.007$ pc. It is also interesting to note the star Gaia EDR3 510911618569239040 with parameters $t_{min}=-2.863\pm0.046$ Myr and $d_{min}=0.057\pm0.079$ pc.

This paper uses linear magnetohydrodynamics to model resonant absorption in coronal plasma with a Cartesian coordinate system. We impose line-tied boundary conditions and tilt the background magnetic field to be oblique to the transition region. Halberstadt & Goedbloed (1993, 1995); Goedbloed& Halberstadt (1994); Arregui et al. (2003) show that line-tied boundary conditions cause their resonant absorption models to produce steep boundary layers/evanescent fast waves. We aim to study the importance of the boundary layers and assess their significance in a solar context. We calculate the solutions in a model where we impose line-tied boundary conditions and compare this with a model where we include the chromosphere instead. Results are calculated analytically and then verified numerically. We show that line-tied boundary conditions can cause the model to overestimate the boundary layers' amplitude significantly. If the fast waves can propagate in the chromosphere, then the line-tied model accurately predicts the boundary layers' amplitude. However, if the fast waves are evanescent, then the boundary layers' size is reduced significantly, and the line-tied model overestimates their amplitude. This leads to the counter-intuitive result that the length scales tangential to the transition region can play an essential role in determining line-tied boundary conditions' validity. The results suggest that line-tied boundary conditions can cause the model to generate unphysically large boundary layers. However, researchers may wish to continue to use them in their models for their simplicity and ability to significantly reduce computation time if they understand and are aware of their flaws.

We report on Chandra, NuSTAR, and MDM observations of two INTEGRAL sources, namely IGR J17528-2022 and IGR J20063+3641. IGR J17528-2022 is an unidentified INTEGRAL source, while IGR J20063+3641 was recently identified as a magnetic cataclysmic variable (mCV) by Halpern et al. (2018). The Chandra observation of IGR J17528-2022 has allowed us to locate the optical counterpart to the source and to obtain its optical spectrum, which shows a strong H$\alpha$ emission line. The optical spectrum and flickering observed in the optical time-series photometry in combination with the X-ray spectrum, which is well fit by an absorbed partially covered thermal bremsstrahlung model, suggests that this source is a strong mCV candidate. The X-ray observations of IGR J20063+3641 reveal a clear modulation with a period of 172.46$\pm0.01$ s, which we attribute to the white dwarf spin period. Additional MDM spectroscopy of the source has also allowed for a clear determination of the orbital period at 0.731$\pm0.015$ d. The X-ray spectrum of this source is also well fit by an absorbed partially covered thermal bremsstrahlung model. The X-ray spectrum, spin periodicity, and orbital periodicity allow this source to be further classified as an intermediate polar.

In a failed supernova, partial ejection of the progenitor's outer envelope can occur due to weakening of the core's gravity by neutrino emission in the protoneutron star phase. We consider emission when this ejecta sweeps up the circumstellar material, analogous to supernova remnants (SNRs). We focus on failed explosions of blue supergiants, and find that the emission can be bright in soft X-rays. Due to its soft emission, we find that sources in the Large Magellanic Cloud (LMC) are more promising to detect than those in the Galactic disk. These remnants are characteristic in smallness ($\lesssim 10$ pc) and slowness (100s of ${\rm km\ s^{-1}}$) compared to typical SNRs. Although the expected number of detectable sources is small (up to a few by eROSITA 4-year all-sky survey), prospects are better for deeper surveys targeting the LMC. Detection of these failed SNRs will realize observational studies of mass ejection upon black hole formation.

To fully characterize the atmospheres, or lack thereof, of terrestrial exoplanets we must include the high-energy environments provided by their host stars. The nearby mid-M dwarf LHS 3844 hosts a terrestrial world which lacks a substantial atmosphere. We present a time series UV spectrum of LHS 3844 from 1131-3215A captured by HST/COS. We detect one flare in the FUV, which has an absolute energy of 8.96+/-0.79e28 erg and an equivalent duration of 355+/-31 s. We extract the flare and quiescent UV spectra separately. For each spectrum we estimate the Ly-alpha flux using correlations between UV line strengths. We use Swift-XRT to place an upper limit on the soft X-ray flux and construct a differential emission model (DEM) to estimate flux that is obscured by the interstellar medium. We compare the DEM flux estimates in the XUV to other methods that rely on scaling from the Ly-alpha, Si IV, and N V lines in the UV. The XUV, FUV, and NUV flux of LHS 3844 relative to its bolometric luminosity is log10(Lband/LBol) = -3.65, -4.16, and -4.56, respectively, for the quiescent state. These values agree with trends in high-energy flux as a function of stellar effective temperature found by the MUSCLES survey for a sample of early-M dwarfs. Many of the most spectroscopically accessible terrestrial exoplanets orbit inactive mid- to late-M dwarfs like LHS 3844. Measurements of M dwarf high-energy spectra are preferable for exoplanet characterization, but are not always possible. The spectrum of LHS 3844 is a useful proxy for the current radiation environment for these worlds.

Stars of spectral types O and B produce neutron stars (NSs) after supernova explosions. Most of NSs are strongly magnetised including normal radio pulsars with $B \propto 10^{12}$ G and magnetars with $B\propto 10^{14}$ G. A fraction of 7-12 per cent of massive stars are also magnetised with $B\propto 10^3$ G and some are weakly magnetised with $B\propto 1$ G. It was suggested that magnetic fields of NSs could be the fossil remnants of magnetic fields of their progenitors. This work is dedicated to study this hypothesis. First, we gather all modern precise measurements of surface magnetic fields in O, B and A stars. Second, we estimate parameters for log-normal distribution of magnetic fields in B stars and found $\mu_B = 2.83\pm 0.1$ $\log_{10}$ (G), $\sigma_B=0.65\pm 0.09$ for strongly magnetised and $\mu_B = 0.14\pm 0.5$ $\log_{10}$ (G), $\sigma=0.7_{-0.27}^{+0.57}$ for weakly magnetised. Third, we assume that the magnetic field of pulsars and magnetars have $2.7$ DEX difference in magnetic fields and magnetars represent 10 per cent of all young NSs and run population synthesis. We found that it is impossible to simultaneously reproduce pulsars and magnetars populations if the difference in their magnetic fields is 2.7 DEX. Therefore, we conclude that the simple fossil origin of the magnetic field is not viable for NSs.

We report that Ly$\alpha$-emitting galaxies (LAEs) may not faithfully trace the cosmic web of neutral hydrogen (HI), but their distribution is likely biased depending on the viewing direction. We calculate the cross-correlation (CCF) between galaxies and Ly$\alpha$ forest transmission fluctuations on the near and far sides of the galaxies separately, for three galaxy samples at $z\sim2$: LAEs, [OIII] emitters (O3Es), and continuum-selected galaxies. We find that only LAEs have anisotropic CCFs, with the near side one showing lower signals up to $r=3-4~h^{-1}$ comoving Mpc. This means that the average HI density on the near side of LAEs is lower than that on the far-side by a factor of $2.1$ under the Fluctuating Gunn-Peterson Approximation. Mock LAEs created by assigning Ly$\alpha$ equivalent width ($EW_\text{Ly$\alpha$}^\text{obs}$) values to O3Es with an empirical relation also show similar, anisotropic CCFs if we use only objects with higher $EW_\text{Ly$\alpha$}^\text{obs}$ than a certain threshold. These results indicate that galaxies on the far side of a dense region are more difficult to be detected ("hidden") in Ly$\alpha$ because Ly$\alpha$ emission toward us is absorbed by dense neutral hydrogen. If the same region is viewed from a different direction, a different set of LAEs will be selected as if galaxies are playing hide-and-seek using HI gas. Care is needed when using LAEs to search for overdensities.

We extend the linear analysis of the drag instability in a 1D perpendicular isothermal C-shock by Gu \& Chen to 2D perpendicular and oblique C-shocks in the typical environment of star-forming clouds. Simplified dispersion relations are derived for the unstable modes. We find that the mode property of the drag instability depends generally on the ratio of the transverse (normal to the shock flow) to longitudinal (along the shock flow) wavenumber. For the transversely large-scale mode, the growth rate and wave frequency of the drag instability in a 2D shock resemble those in a 1D shock. For the transversely small-scale mode, the drag instability is characterized by an unstable mode coupled with an acoustic mode primarily along the transverse direction. When the shock is perpendicular or less oblique, there exists a slowly propagating mode, which can potentially grow into a nonlinear regime and contribute to the maximum growth of the instability. In contrast, when the shock is more oblique, this slowly propagating unstable mode disappears and the maximum growth of the drag instability is likely contributed from the transversely large-scale mode (i.e., almost 1D mode). In all cases that we consider, the magnitude of the density perturbations is significantly larger than that of the velocity and magnetic field perturbations, implying that the density enhancement governs the dynamics in the linear regime of the instability. A few issues in the linear analysis as well as the possible astrophysical implications are also briefly discussed.

Isocyanic acid (HNCO) is a simple molecule with a potential to form prebiotic and complex organic species. Using a spectral survey collected with the Atacama Pathfinder EXperiment (APEX), in this work we report the detection of 42 transitions of HNCO in the hot molecular core/outflow G331.512-0.103 (hereafter G331). The spectral lines were observed in the frequency interval $\sim$ 160 - 355 GHz. By means of Local Thermodynamic Equilibrium (LTE) analyses, applying the rotational diagram method, we studied the excitation conditions of HNCO. The excitation temperature and column density are estimated to be $T_{ex}$ = 58.8 $\pm$ 2.7 K and $N$ = (3.7 $\pm$ 0.5) $\times$ 10$^{15}$ cm$^{-2}$, considering beam dilution effects. The derived relative abundance is between (3.8 $\pm$ 0.5) $\times $10$^{-9}$ and (1.4 $\pm$ 0.2) $\times $10$^{-8}$. In comparison with other hot molecular cores, our column densities and abundances are in agreement. An update of the internal partition functions of the four CHNO isomers: HNCO; cyanic acid, HOCN; fulminic acid, HCNO; and isofulminic acid, HONC is provided. We also used the astrochemical code Nautilus to model and discuss HNCO abundances. The simulations could reproduce the abundances with a simple zero-dimensional model at a temperature of 60 K and for a chemical age of $\sim$ 10$^5$ years, which is larger than the estimated dynamical age for G331. This result could suggest the need for a more robust model and even the revision of chemical reactions associated with HNCO.

The melt productivity of a differentiated planet's mantle is primarily controlled by its iron content, which is itself approximated by the planet's core mass fraction (CMF). Here we show that estimates of an exo-planet's CMF allows robust predictions of the thickness, composition and mineralogy of the derivative crust. These predicted crustal compositions allow constraints to be placed on volatile cycling between surface and the deep planetary interior, with implications for the evolution of habitable planetary surfaces. Planets with large, terrestrial-like, CMFs ($\geq$0.32) will exhibit thin crusts that are inefficient at transporting surface water and other volatiles into the underlying mantle. By contrast, rocky planets with smaller CMFs ($\leq$0.24) and higher, Mars-like, mantle iron contents will develop thick crusts capable of stabilizing hydrous minerals, which can effectively sequester volatiles into planetary interiors and act to remove surface water over timescales relevant to evolution. The extent of core formation has profound consequences for the subsequent planetary surface environment and may provide additional constraints in the hunt for habitable, Earth-like exo-planets.

The prevailing regulatory framework for light pollution control is based on establishing conditions on individual light sources or single installations (regarding features like ULOR, spectrum, illuminance levels, glare, ...), in the hope that an ensemble of individually correct lighting installations will be effective to somehow solve this problem. This "local sources" approach is indeed necessary, and shall no doubt be enforced; however, it seems to be clearly insufficient for curbing the actual process of degradation of the night, and for effectively attaining the necessary remediation goals. In this paper we describe a complementary (not substitutive) 'red-lines' strategy that should in our opinion be adopted as early as possible in the policies for light pollution control. This top-down approach seeks to set definite limits on the allowable degradation of the night, providing the methodological tools required for making science-informed public policy decisions and for managing the transition processes. Light pollution abatement should routinely be included as an integral part of any territorial management plan. A practical application case-study is described to illustrate these concepts.

Technological advances in instrumentation have led to an exponential increase in exoplanet detection and scrutiny of stellar features such as spots and faculae. While the spots and faculae enable us to understand the stellar dynamics, exoplanets provide us with a glimpse into stellar evolution. While the ubiquity of noise (e.g., telluric, instrumental, or photonic) is unavoidable, combining this with increased spectrographic resolution compounds technological challenges. To account for these noise sources and resolution issues, we use a temporal multifractal framework to study data from the SOAP 2.0 tool, which simulates a stellar spectrum in the presence of a spot, a facula or a planet. Given these controlled simulations, we vary the resolution as well as the signal-to-noise (S/N) ratio to obtain a lower limit on the resolution and S/N required to robustly detect features. We show that a spot and a facula with a 1% coverage of the stellar disk can be robustly detected for a S/N (per pixel) of 35 and 60 respectively, for any spectral resolution above 20,000, while a planet with a radial velocity (RV) of 10 m/s can be detected for a S/N (per pixel) of 600. Rather than viewing noise as an impediment, our approach uses noise as a source of information.

We update Halo Zeldovich Perturbation Theory (HZPT), an analytic model for the two-point statistics of dark matter, to describe halo and galaxy clustering, and galaxy-matter cross-correlation on nonlinear scales. The model correcting Zeldovich has an analytic Fourier transform, and therefore is valid in both configuration space and Fourier space. The model is accurate at the $2\%$-level or less for $P_{mm}$ (k < 1 h/Mpc), $P_{hm}$ (k < 1 h/Mpc), $P_{hh}$ (k < 2 h/Mpc), $P_{gm}$ (k < 1 h/Mpc), $P_{gg}$ (k < 1 h/Mpc), $\xi_{mm}$ (r > 1 Mpc/h), $\xi_{hm}$ (r > 2 Mpc/h), $\xi_{hh}$ (r > 2 Mpc/h), $\xi_{gm}$ (r > 1 Mpc/h), $\xi_{gg}$ (r > 2 Mpc/h), for LRG-like mock galaxies. We show that the HZPT model for matter correlators can account for the effects of a wide range of baryonic feedback models and provide extended dark matter models which are of $1\% ~(3\%)$ accuracy for k < 10 (8) h/Mpc. We explicitly model the non-perturbative features of halo exclusion for the halo-halo and galaxy-galaxy correlators, as well as the presence of satellites for galaxy-matter and galaxy-galaxy correlation functions. We perform density estimation using N-body simulations and a wide range of HOD galaxy mocks to obtain correlations of model parameters with the cosmological parameters $\Omega_{m}$ and $\sigma_{8}$. HZPT can provide a fast, interpretable, and analytic model for combined-probe analyses of redshift surveys using scales well into the non-linear regime.

We present new spectroscopic observations of Ly$\alpha$ Blob 2 (LAB2) in the SSA22 protocluster region ($z \sim$ 3.1). By creating a narrow-band Ly$\alpha$ image, we observed extended Ly$\alpha$ emission in three distinct regions, in which the highest Ly$\alpha$ surface brightness (SB) center is far away from the known continuum sources. We have searched through the MOSFIRE slits that cover the high Ly$\alpha$ SB regions, but are unable to detect any significant nebular emission near the highest SB center. We further map the blue-to-red flux ratio and find that it is anti-correlated with Ly$\alpha$ SB with a power-law index of $\sim$ -0.4. To decode the spatially-resolved Ly$\alpha$ profiles using Monte-Carlo radiative transfer (MCRT) modelling, we use both multiphase, clumpy models and shell models and successfully reproduced the diverse Ly$\alpha$ morphologies with reasonable physical parameters. Significant correlations exist between parameters of two different models, and our multiphase, clumpy model parameters naturally alleviated the previously reported discrepancies between the shell model parameters and data. In addition, we have modeled Ly$\alpha$ spectra at different spatial positions simultaneously, and we find that the variation of the inferred clump outflow velocity can be approximately explained by line-of-sight projection effects. Our results support the `central powering + scattering' scenario, i.e. the Ly$\alpha$ photons are generated by a central powering source and then scatter with outflowing, multiphase HI gas while propagating outwards. The infalling of cool gas near the blob outskirts shapes the observed blue-dominated Ly$\alpha$ profiles, but its energy contribution is likely to be minor compared to the photo-ionization by star-forming galaxies and/or AGNs.

The Gross-Pitaevskii-Poisson equations that govern the evolution of self-gravitating Bose-Einstein condensates, possibly representing dark matter halos, experience a process of gravitational cooling and violent relaxation. We propose a heuristic parametrization of this complicated process in the spirit of Lynden-Bell's theory of violent relaxation for collisionless stellar systems. We derive a generalized wave equation (that was introduced phenomenologically in [P.H. Chavanis, Eur. Phys. J. Plus {\bf 132}, 248 (2017)]) involving a logarithmic nonlinearity associated with an effective temperature $T_{\rm eff}$ and a damping term associated with a friction $\xi$. These terms can be obtained from a maximum entropy production principle and are linked by a form of Einstein relation expressing the fluctuation-dissipation theorem. The wave equation satisfies an $H$-theorem for the Lynden-Bell entropy and relaxes towards a stable equilibrium state which is a maximum of entropy at fixed mass and energy. This equilibrium state represents the most probable state of a Bose-Einstein condensate dark matter halo. It generically has a core-halo structure. The quantum core prevents gravitational collapse and may solve the core-cusp problem. The isothermal halo leads to flat rotation curves in agreement with the observations. These results are consistent with the phenomenology of dark matter halos. Furthermore, as shown in a previous paper [P.H. Chavanis, Phys. Rev. D {\bf 100}, 123506 (2019)], the maximization of entropy with respect to the core mass at fixed total mass and total energy determines a core mass--halo mass relation which agrees with the relation obtained in direct numerical simulations. We stress the importance of using a microcanonical description instead of a canonical one. We also explain how our formalism can be applied to the case of fermionic dark matter halos.

We reconsider the dynamical systems approach to analyze inflationary universe in the Jordan frame models of scalar field nonminimally coupled to curvature. The adopted set of variables allows us to clearly distinguish between different asymptotic states in the phase space, including the kinetic and inflationary regimes. Inflation is realized as a heteroclinic trajectory originating either at infinity from a nonhyperbolic asymptotic de Sitter point or from a regular saddle de Sitter point. We also present a comprehensive picture of possible initial conditions leading to sufficient inflationary expansion and show their extent on the phase diagrams. In addition we determine the correct slow roll conditions applicable in the Jordan frame and show how they approximate the leading inflationary "attractor solution". As particular examples we portrait quadratic and quartic potential models and note that increasing the nonminimal coupling diminishes the range of good initial conditions in the quadratic case, but enlarges is in the quartic case.

A method involving intensity correlation measurements is described, which allows for the complete removal of Doppler broadening in the emission of electromagnetic radiation from far-away sources that are inaccessible to conventional Doppler-free measurements. The technique, relying on a correction to g(2) of order N-1, probes the separation between neighboring spectral lines and is also applicable to the elimination of broadening due to collisions (N is the number of emitting particles and g(2) is the second-order field correlation function). Possible applications include a determination of cosmological parameters from red shifts of gravitationally-lensed quasars.

In this paper we perform systematic investigation of all possible exponential solutions in Einstein-Gauss-Bonnet gravity with the spatial section being a product of two subspaces. We describe a scheme which always allow to find solution for a given $\{p, q\} > 2$ (number of dimensions of two subspaces) and $\zeta$ (ratio of the expansion rates of these two subspaces). Depending on the parameters, for given $\{\alpha, \Lambda\}$ (Gauss-Bonnet coupling and cosmological constant) there could be up to four distinct solutions (with different $\zeta$'s). Stability requirement introduces relation between $\zeta$, $\{p, q\}$ and sign of the expansion rate. Nevertheless, for any $\{p, q\} > 2$ we can always choose sign for expansion rates so that the resulting solution would be stable. The scheme for finding solutions is described and the bounds on the parameters are drawn. Specific cases with $\{p, q\} = \{1, 2\}$ are also considered. Finally, we separately described physically sensible case with one of the subspaces being three-dimensional and expanding (resembling our Universe) while another to be contracting (resembling extra dimensions), describing successful compactification; for this case we also drawn bounds on the parameters where such regime occurs.

In this paper, we study the preheating duration modeled by an e-folding number $N_{\rm pre}$ and an effective equation-of-state parameter $\omega_{\rm pre}$. We show that a preheating equation-of-state parameter $\omega_{\rm pre} = 0$ implies more post-inflationary e-folds of expansion. Considering the Higgs and polynomial potentials, we conclude that the chaotic and the Higgs models show good compatibility with the observational data according to the Planck 2018 results.