Papers for Friday, Jan 22 2021

Common-envelope evolution is important in the formation of neutron star binaries within the isolated binary formation channel. As a neutron star inspirals within the envelope of a primary massive star, it accretes and spins up. Because neutron stars are in the strong-gravity regime, they have a substantial relativistic mass deficit, i.e., their gravitational mass is less than their baryonic mass. This effect causes some fraction of the accreted baryonic mass to convert into neutron star binding energy. The relativistic mass deficit also depends on the nuclear equation of state, since more compact neutron stars will have larger binding energies. We model the mass growth and spin-up of neutron stars inspiraling within common-envelope environments and quantify how different initial binary conditions and hadronic equations of state affect the post-common-envelope neutron star's mass and spin. From these models, we find that neutron star mass growth is suppressed by $\approx 15-30\%$. We also find that for a given amount of accreted baryonic mass, more compact neutron stars will spin-up faster while gaining less gravitational mass, and vice versa. This work demonstrates that a neutron star's strong gravity and nuclear microphysics plays a role in neutron-star-common-envelope evolution, in addition to the macroscopic astrophysics of the envelope. Strong gravity and the nuclear equation of state may thus affect both the population properties of neutron star binaries and the cosmic double neutron star merger rate.

The direct imaging of extrasolar giant planets demands the highest possible contrasts (dH ~10 magnitudes) at the smallest angular separations (~0.1'') from the star. We present an adaptive optics observing method, called star-hopping, recently offered as standard queue observing for the SPHERE instrument at the VLT. The method uses reference difference imaging (RDI) but unlike earlier works, obtains images of a reference star for PSF subtraction, within minutes of observing the target star. We aim to significantly gain in contrast over the conventional angular differencing imaging (ADI) method, to search for a fifth planet at separations less than 10 au, interior to the four giant planets of the HR 8799 system. We obtained a total of 4.5 hours of simultaneous integral field spectroscopy (R~30, Y-H band with IFS) and dual-band imaging (K1 and K2-band with IRDIS) of the HR 8799 system and a reference star. The reference star was observed for ~1/3 of the total time, and should have dR~1 mag and separated on sky by ~1-2 deg. The star hops were made every 6-10 minutes, with only 1 minute gaps in on-sky integration per hop. We did not detect the hypothetical fifth planet at the most plausible separations, 7.5 and 9.7 au, down to mass limits of 3.6 MJup high signal-to-noise ratios. As noted in previous works, the planet spectra are matched very closely by some red field dwarfs. We also demonstrated that with star-hopping RDI, the contrast improvement at 0.1'' separation can be up to 2 magnitudes. Since ADI, meridian transit and the concomitant sky rotation are not needed, the time of observation can be chosen from within a 2-3 times larger window. In general, star-hopping can be used for stars fainter than R=4 magnitudes, since for these a reference star of suitable brightness and separation is usually available. The reduction software used in this paper has been made available online.

Molecular outflows contributing to the matter cycle of star forming galaxies are now observed in small and large systems at low and high redshift. Their physical origin is still unclear. In most theoretical studies only warm ionised/neutral and hot gas outflowing from the interstellar medium is generated by star formation. We investigate an in-situ H$_2$ formation scenario in the outflow using high-resolution simulations, including non-equilibrium chemistry and self-gravity, of turbulent, warm, and atomic clouds with densities 0.1, 0.5 and $1\,\mathrm{cm}^{-3}$ exposed to a magnetised hot wind. For cloud densities $\gtrsim 0.5\,\mathrm{cm}^{-3}$ a magnetised wind triggers H$_2$ formation before cloud dispersal. Up to 3 per cent of the initial cloud mass can become molecular on $\sim 10\,\mathrm{Myr}$ time scales. The effect is stronger for winds with perpendicular $B$-fields and intermediate density clouds ($n_\mathrm{c}\sim 0.5\,\mathrm{cm}^{-3}$). Here H$_2$ formation can be boosted by up to one order of magnitude compared to isolated cooling clouds independent of self-gravity. Self-gravity preserves the densest clouds way past their $\sim 15\,\mathrm{Myr}$ cloud crushing time scales. This model could provides a plausible in-situ origin for the observed molecular gas. Warm ionised gas is also generated, almost independent of the cloud density. The amount solely depend on the magnetic field configuration in the wind. For low density clouds ($0.1\,\mathrm{cm}^{-3}$), the forming warm ionised gas can be as much as 60 per cent of the initially atomic cloud mass. This could contribute to observations of outflows with ionised gas sensitive tracers.

We present observations of the dwarf galaxies GALFA Dw3 and GALFA Dw4 with the Advanced Camera for Surveys (ACS) on the Hubble Space Telescope (HST). These galaxies were initially discovered as optical counterparts to compact HI clouds in the GALFA survey. Both objects resolve into stellar populations which display an old red giant branch, younger helium burning, and massive main sequence stars. We use the tip of the red giant branch method to determine the distance to each galaxy, finding distances of 7.61$_{-0.29}^{+0.28}$ Mpc and 3.10$_{-0.17}^{+0.16}$ Mpc, respectively. With these distances we show that both galaxies are extremely isolated, with no other confirmed objects within ~1.5 Mpc of either dwarf. GALFA Dw4 is also found to be unusually compact for a galaxy of its luminosity. GALFA Dw3 and Dw4 contain HII regions with young star clusters and an overall irregular morphology; they show evidence of ongoing star formation through both ultraviolet and H$\alpha$ observations and are therefore classified as dwarf irregulars (dIrrs). The star formation histories of these two dwarfs show distinct differences: Dw3 shows signs of a recently ceased episode of active star formation across the entire dwarf, while Dw4 shows some evidence for current star formation in spatially limited HII regions. Compact HI sources offer a promising method for identifying isolated field dwarfs in the Local Volume, including GALFA Dw3 & Dw4, with the potential to shed light on the driving mechanisms of dwarf galaxy formation and evolution.

Binary neutron star mergers are thought to be one of the dominant sites of production for rapid neutron capture elements, and the source of platinum and gold in the Universe. Since the discovery of the binary neutron star merger GW170817, and its associated kilonova AT2017gfo, numerous works have attempted to determine the composition of its outflowing material, but they have been hampered by the lack of complete atomic data. Here we demonstrate how inclusion of new atomic data in synthetic spectra calculations can provide insights and constraints on the production of the heaviest elements. We employ theoretical atomic data for neutral, singly- and doubly-ionised platinum and gold, to generate photospheric and simple nebular-phase model spectra for kilonova-like ejecta properties. We make predictions for the locations of strong transitions, which could feasibly appear in the spectra of kilonovae that are rich in these species. We use GRASP0 to generate the underlying atomic structure and TARDIS to model the diffusion phase showing that the strongest features lie in the ultra-violet region. We identify low-lying electric quadrupole and magnetic dipole transitions that may give rise to forbidden lines when the ejecta becomes optically thin. The strongest lines lie beyond 8000 Angstroms, motivating high quality near-infrared spectroscopic follow-up of kilonova candidates. We compare our model spectra to the observed spectra of AT2017gfo, and conclude that no platinum or gold signatures are prominent in the ejecta. This work demonstrates how new atomic data of heavy elements can be included in radiative transfer calculations, and motivates future searches for elemental signatures.

Cosmography is a powerful tool to investigate the Universe kinematic and then to reconstruct dynamics in a model-independent way. However, recent new measurements of supernovae Ia and quasars have populated the Hubble diagram up to high redshifts ($z \sim 7.5$) and the application of the traditional cosmographic approach has become less straightforward due to the large redshifts implied. Here we investigate this issue through an expansion of the luminosity distance-redshift relation in terms of "orthogonal" logarithmic polynomials. In particular we point out the advantages of a new procedure of "orthogonalization" and we show that such an expansion provides a very good fit in the whole $z=0\div 7.5$ range to both real and mock data obtained assuming various cosmological models. Moreover, despite of the fact that the cosmographic series is tested well beyond its convergence radius, the parameters obtained expanding the luminosity distance - redshift relation for the $\Lambda$CDM model are broadly consistent with the results from a fit of mock data obtained with the same cosmological model. This provides a method to test the reliability of a cosmographic function to study cosmological models at high redshifts and it demonstrates that the logarithmic polynomial series can be used to test the consistency of the $\Lambda$CDM model with the current Hubble diagram of quasars and supernovae Ia. We confirm a strong tension (at $>4\sigma$) between the concordance cosmological model and the Hubble diagram at $z>1.5$. Such a tension is dominated by the contribution of quasars at $z>2$ and starts to be present also in the few supernovae Ia observed at $z>1$.

Dynamical models allow us to connect the motion of a set of tracers to the underlying gravitational potential, and thus to the total (luminous and dark) matter distribution. They are particularly useful for understanding the mass and spatial distribution of dark matter (DM) in a galaxy. Globular clusters (GCs) are an ideal tracer population in dynamical models, since they are bright and can be found far out into the halo of galaxies. We aim to test how well Jeans-Anisotropic-MGE (JAM) models using GCs (positions and line-of-sight velocities) as tracers can constrain the mass and radial distribution of DM halos. For this, we use the E-MOSAICS suite of 25 zoom-in simulations of L* galaxies. We find that the DM halo properties are reasonably well recovered by the JAM models. There is, however, a strong correlation between how well we recover the mass and the radial distribution of the DM and the number of GCs in the galaxy: the constraints get exponentially worse with fewer GCs, and at least 150 GCs are needed in order to guarantee that the JAM model will perform well. We find that while the data quality (uncertainty on the radial velocities) can be important, the number of GCs is the dominant factor in terms of the accuracy and precision of the measurements. This work shows promising results for these models to be used in extragalactic systems with a sample of more than 150 GCs.

We observed the high-mass star-forming region G335.579-0.292 with the Atacama Large Millimeter/submillimeter Array (ALMA) at 226 GHz with an angular resolution of 0.3'' ($\sim 1000$ au resolution at the source distance). G335.579-0.292 hosts one of the most massive cores in the Galaxy (G335-MM1). The continuum emission shows that G335-MM1 fragments into at least five sources, while molecular line emission is detected in two of the continuum sources (ALMA1 and ALMA3). We found evidence of large and small scale infall in ALMA1 revealed by an inverse P-Cygni profile and the presence of a blue-shifted spot at the center of the first moment map of the CH$_3$CN emission. In addition, hot gas expansion in the innermost region is unveiled by a red-shifted spot in the first moment map of HDCO and (CH$_3$)$_2$CO (both with $E_u > 1100$ K). Our modeling reveals that this expansion motion originates close to the central source, likely due to reversal of the accretion flow induced by the expansion of the HII region, while infall and rotation motions originate in the outer regions. ALMA3 shows clear signs of rotation, with a rotation axis inclination with respect to the line of sight close to $90^\circ$, and a system mass (disk + star) in the range of 10-30 M$_\odot$.

This research was stimulated by the recent studies of damping solutions in dynamically stable spherical stellar systems. Using the simplest model of the homogeneous stellar medium, we discuss nontrivial features of stellar systems. Taking them into account will make it possible to correctly interpret the results obtained earlier and will help to set up decisive numerical experiments in the future. In particular, we compare the initial value problem versus the eigenvalue problem. It turns out that in the unstable regime, the Landau-damped waves can be represented as a superposition of van Kampen modes {\it plus} a discrete damped mode, usually ignored in the stability study. This mode is a solution complex conjugate to the unstable Jeans mode. In contrast, the Landau-damped waves are not genuine modes: in modes, eigenfunctions depend on time as $\exp (-{\rm i} \omega t)$, while the waves do not have eigenfunctions on the real $v$-axis at all. However, `eigenfunctions' on the complex $v$-contours do exist. Deviations from the Landau damping are common and can be due to singularities or cut-off of the initial perturbation above some fixed value in the velocity space.

The study of low surface brightness light in large, deep imaging surveys is still uncharted territory as automated data reduction pipelines over-subtract or eliminate this light. Using archival data of the Abell 85 cluster of galaxies taken with Hyper Suprime-Cam on the Subaru Telescope, we show that using careful data processing can unveil the diffuse light within the cluster, the intracluster light. We reach surface brightness limits of $\mu_{g}^{limit}$(3$\sigma$, 10"x10") = 30.9 mag/arcsec$^2$, and $\mu_{i}^{limit}$(3$\sigma$, 10"x10") = 29.7 mag/arcsec$^2$. We measured the radial surface brightness profiles of the brightest cluster galaxy out to the intracluster light (radius $\sim215$ kpc), for the g and i bands. We found that both the surface brightness and the color profiles become shallower beyond $\sim75$ kpc suggesting that a distinct component, the intracluster light, starts to dominate at that radius. The color of the profile at $\sim100$ kpc suggests that the buildup of the intracluster light of Abell 85 occurs by the stripping of massive ($\sim10^{10}M_{\odot}$) satellites. The measured fraction of this light ranges from 8% to 30% in g, depending on the definition of intracluster light chosen.

We present detailed observations of photoionization conditions and galaxy kinematics in eleven z$=1.39-2.59$ radio-loud quasar host galaxies. Data was taken with OSIRIS integral field spectrograph (IFS) and the adaptive optics system at the W.M. Keck Observatory that targeted nebular emission lines (H$\beta$,[OIII],H$\alpha$,[NII]) redshifted into the near-infrared (1-2.4 \micron). We detect extended ionized emission on scales ranging from 1-30 kpc photoionized by stars, shocks, and active galactic nuclei (AGN). Spatially resolved emission-line ratios indicate that our systems reside off the star formation and AGN-mixing sequence on the Baldwin, Phillips $\&$ Terlevich (BPT) diagram at low redshift. The dominant cause of the difference between line ratios of low redshift galaxies and our sample is due to lower gas-phase metallicities, which are 2-5$\times$ less compared to galaxies with AGN in the nearby Universe. Using gas velocity dispersion as a proxy to stellar velocity dispersion and dynamical mass measurement through inclined disk modeling we find that the quasar host galaxies are under-massive relative to their central supermassive black hole (SMBH) mass, with all systems residing off the local scaling ($M_{\bullet}-\sigma~$,$M_{\bullet}-M_{*}~$) relationship. These quasar host galaxies require substantial growth, up to an order of magnitude in stellar mass, to grow into present-day massive elliptical galaxies. Combining these results with part I of our sample paper (Vayner et al. 2021) we find evidence for winds capable of causing feedback before the AGN host galaxies land on the local scaling relation between black hole and galaxy stellar mass, and before the enrichment of the ISM to a level observed in local galaxies with AGN.

Central compact objects are young neutron stars emitting thermal X-rays with bolometric luminosities $L_X$ in the range $10^{32}$-$10^{34}$ erg/s. Gourgouliatos, Hollerbach and Igoshev recently suggested that peculiar emission properties of central compact objects can be explained by tangled magnetic field configurations formed in a stochastic dynamo during the proto-neutron star stage. In this case the magnetic field consists of multiple small-scale components with negligible contribution of global dipolar field. We study numerically three-dimensional magneto-thermal evolution of tangled crustal magnetic fields in neutron stars. We find that all configurations produce complicated surface thermal patterns which consist of multiple small hot regions located at significant separations from each other. The configurations with initial magnetic energy of $2.5-10\times 10^{47}$ erg have temperatures of hot regions that reach $\approx 0.2$ keV, to be compared with the bulk temperature of $\approx 0.1$ keV in our simulations with no cooling. A factor of two in temperature is also seen in observations of central compact objects. The hot spots produce periodic modulations in light curve with typical amplitudes of $\leq 9-11$ %. Therefore, the tangled magnetic field configuration can explain thermal emission properties of some central compact objects.

A small fraction of gamma-ray bursts (GRBs) exhibit blackbody emission in the X-ray spectra, the origin of which is debated. In order to gain a more complete understanding of this phenomenon, we present a search for blackbody components in 116 GRBs with known redshifts observed by {\it Swift}~XRT. A time-resolved spectral analysis is carried out and the significance of the blackbody is assessed with respect to an absorbed power-law model. We report nine new detections and confirm the previously reported blackbody in GRB~171205A. Together with our previous results, there are a total of 19 GRBs with significant blackbody emission in a sample of 199 GRBs observed by {\it Swift} over 13 years. The detections include one short GRB and two low-luminosity GRBs. We estimate fireball parameters from the blackbody components and note that the blackbody luminosity is correlated with the temperature and inferred Lorentz factor. There is a large spread in the properties of the blackbody components and the light curves, which points to different origins for the emission. In about a third of the GRBs, the blackbody is clearly associated with late prompt emission from the jet. The rest of the sample includes cases that are fully consistent with the expectations from a cocoon, as well cases that may be explained by high-latitude emission or more energetic cocoons. These results indicate that thermal emission is associated with all parts of the jet.

We have carried out a search for massive white dwarfs (WDs) in the direction of young open star clusters using the Gaia DR2 database. The aim of this survey was to provide robust data for new and previously known high-mass WDs regarding cluster membership, to highlight WDs previously included in the Initial Final Mass Relation (IFMR) that are unlikely members of their respective clusters according to Gaia astrometry and to select an unequivocal WD sample that could then be compared with the host clusters' turnoff masses. All promising WD candidates in each cluster CMD were followed up with spectroscopy from Gemini in order to determine whether they were indeed WDs and derive their masses, temperatures and ages. In order to be considered cluster members, white dwarfs were required to have proper motions and parallaxes within 2, 3, or 4-$\sigma$ of that of their potential parent cluster based on how contaminated the field was in their region of the sky, have a cooling age that was less than the cluster age and a mass that was broadly consistent with the IFMR. A number of WDs included in current versions of the IFMR turned out to be non-members and a number of apparent members, based on Gaia's astrometric data alone, were rejected as their mass and/or cooling times were incompatible with cluster membership. In this way, we developed a highly selected IFMR sample for high mass WDs that, surprisingly, contained no precursor masses significantly in excess of ${\sim}$6 $M_{\odot}$.

The multi analytical study of terrestrial analogues is a useful strategy to deepen the knowledge about the geological and environmental evolution of Mars and other extraterrestrial bodies. In spite of the increasing importance that LIBS, NIR and Raman techniques are acquiring in the field of space exploration, there is a lack web-based platform providing free access to a wide multi-spectral database of terrestrial analogue materials. The Planetary Terrestrial Analogue Library (PTAL) project aims at responding to this critical need by developing and providing free web accessibility to LIBS, NIR and Raman data from more than 94 terrestrial analogues selected according to their congruence with Martian geological contexts. In this framework, the present manuscript provides the scientific community with a complete overview of the over 4500 Raman spectra collected to feed the PTAL database. Raman data, obtained through the complementary use of laboratory and spacecraft-simulator systems, confirmed the effectiveness of this spectroscopic technique for the detection of major and minor mineralogical phases of the samples, the latter being of critical importance for the recognition of geological processes that could have occurred on Mars and other planets. In light of the forthcoming missions to Mars, the results obtained through the RLS ExoMars Simulator offer a valuable insight on the scientific outcome that could derive from the RLS spectrometer that will soon land on Mars as part of the ExoMars rover payload.

We present and validate a new network of atmospheric thermo-chemical and photo-chemical sulfur reactions. We use a 1-D chemical kinetics model to investigate these reactions as part of a broader HCNO chemical network in a series of hot and warm Jupiters. We find that temperatures approaching $1400\,\mathrm{K}$ are favourable for the production of H2S and HS around $\mathrm{10^{-3}\,bar}$, the atmospheric level where detection by transit spectroscopy may be possible, leading to mixing ratios of around $10^{-6}$. At lower temperatures, down to $1000\,\mathrm{K}$, the abundance of S2 can be up to a mixing ratio of $10^{-5}$ at the same pressure, at the expense of H2S and HS, which are depleted down to a mixing ratio of $10^{-6}$. We also investigate how the inclusion of sulfur can manifest in an atmosphere indirectly, by its effect on the abundance of non-sulfur-bearing species. We find that in a model of the atmosphere of HD 209458 b, the inclusion of sulfur can lower the abundance of NH3, CH4 and HCN by up to two orders of magnitude around $\mathrm{10^{-3}\,bar}$. In the atmosphere of the warm Jupiter 51 Eri b, we additionally find the inclusion of sulphur depletes the peak abundance of CO2 by a factor of five, qualitatively consistent with prior models. We note that many of the reactions used in the network have poorly determined rates, especially at higher temperatures. To obtain a truly accurate idea of the impact of sulfur chemistry in hot and warm Jupiter atmospheres, new measurements of these reaction rates must take place.

We present the fourth catalog of serendipitously discovered compact extragalactic emission-line sources -- H$\alpha$ Dots. A total of 454 newly discovered objects are included in the current survey list. These objects have been detected in searches of moderately deep narrow-band images acquired for the ALFALFA H$\alpha$ project (Van Sistine et al. 2016). The catalog of H-alpha Dots presented in the current paper was derived from searches carried out using ALFALFA H$\alpha$ images obtained with the KPNO 2.1 m telescope. This results in a substantially deeper sample of Dots compared to our previous lists, which were all discovered in images taken with the WIYN 0.9 m telescope. The median R-band magnitude of the current catalog is 21.59, more than 1.6 magnitudes fainter than the median for the 0.9~m sample (factor of 4.4x fainter). Likewise, the median emission-line flux of the detected sources is a factor of 4.3x fainter. The line-flux completeness limit of the current sample is approximately 3 x 10$^{-16}$ erg/s/cm$^2$. We present accurate coordinates, apparent magnitudes and narrow-band line fluxes for each object in the sample. Unlike our previous lists of H$\alpha$ Dots, the current sample does not include follow-up spectroscopy.

We explore changes in the adiabatic low-order g-mode pulsation periods of 0.526, 0.560, and 0.729 M$_\odot$ carbon-oxygen white dwarf models with helium-dominated envelopes due to the presence, absence, and enhancement of $^{22}$Ne in the interior. The observed g-mode pulsation periods of such white dwarfs are typically given to 6$-$7 significant figures of precision. Usually white dwarf models without $^{22}$Ne are fit to the observed periods and other properties. The root-mean-square residuals to the $\simeq$ 150$-$400 s low-order g-mode periods are typically in the range of $\sigma_{rms}$ $\lesssim$ 0.3 s, for a fit precision of $\sigma_{rms}/ P$ $\lesssim$ 0.3 %. We find average relative period shifts of $\Delta P/P$ $\simeq$ $\pm$ 0.5 % for the low-order dipole and quadrupole g-mode pulsations within the observed effective temperature window, with the range of $\Delta P/P$ depending on the specific g-mode, abundance of $^{22}$Ne, effective temperature, and mass of the white dwarf model. This finding suggests a systematic offset may be present in the fitting process of specific white dwarfs when $^{22}$Ne is absent. As part of the fitting processes involves adjusting the composition profiles of a white dwarf model, our study on the impact of $^{22}$Ne can provide new inferences on the derived interior mass fraction profiles. We encourage routinely including $^{22}$Ne mass fraction profiles, informed by stellar evolution models, to future generations of white dwarf model fitting processes.

As host to two accreting planets, PDS 70 provides a unique opportunity to probe the chemical complexity of atmosphere-forming material. We present ALMA Band 6 observations of the PDS~70 disk and report the first chemical inventory of the system. With a spatial resolution of $0.4''-0.5''$ ( $\sim50\,$au), 12 species are detected, including CO isotopologues and formaldehyde, small hydrocarbons, HCN and HCO$^+$ isotopologues, and S-bearing molecules. SO and CH$_3$OH are not detected. All lines show a large cavity at the center of the disk, indicative of the deep gap carved by the massive planets. The radial profiles of the line emission are compared to the (sub-)mm continuum and infrared scattered light intensity profiles. Different molecular transitions peak at different radii, revealing the complex interplay between density, temperature and chemistry in setting molecular abundances. Column densities and optical depth profiles are derived for all detected molecules, and upper limits obtained for the non detections. Excitation temperature is obtained for H$_2$CO. Deuteration and nitrogen fractionation profiles from the hydro-cyanide lines show radially increasing fractionation levels. Comparison of the disk chemical inventory to grids of chemical models from the literature strongly suggests a disk molecular layer hosting a carbon to oxygen ratio C/O$>$1, thus providing for the first time compelling evidence of planets actively accreting high C/O ratio gas at present time.

We argue that a natural explanation for a variety of robust galaxy scaling relations comes from the perspective of pattern formation and self-organization as a result of symmetry breaking. We propose a simple Lagrangian model that combines a conventional model for normal matter in a galaxy with a conventional model for stripe pattern formation in systems that break continuous translation invariance. We show that the energy stored in the pattern field acts as an effective dark matter. Our theory reproduces the gross features of elliptic galaxies as well as disk galaxies (HSB and LSB) including their detailed rotation curves, the radial acceleration relation (RAR), and the Freeman law. We investigate the stability of disk galaxies in the context of our model and obtain scaling relations for the central dispersion for elliptical galaxies. A natural interpretation of our results is that (1) `dark matter' is potentially a collective, emergent phenomenon and not necessarily an as yet undiscovered particle, and (2) MOND is an effective theory for the description of a self-organized complex system rather than a fundamental description of nature that modifies Newton's second law.

The baryon acoustic oscillation feature can be used as a standard cosmological ruler. In practice, for sub-percent level accuracy on the distance scale, it must be standardized. The physical reason why is understood, so we use this to develop an algorithm which improves the estimated scale. The algorithm exploits the fact that, over the range of scales where the initial correlation function is well-fit by a polynomial, the leading order effects which distort the length of the ruler can be accounted for analytically. Tests of the method in numerical simulations show that it provides simple and fast reconstruction of the full shape of the BAO feature, as well as subpercent determination of the linear point in the correlation function of biased tracers with minimal assumptions about the underlying cosmological model or the nature of the observed tracers.

Starting in 2008, NASA has provided the exoplanet community an observational program aimed at obtaining the highest resolution imaging available as part of its mission to validate and characterize exoplanets, as well as their stellar environments, in search of life in the universe. Our current program uses speckle interferometry in the optical (320-1000 nm) with new instruments on the 3.5-m WIYN and both 8-m Gemini telescopes. Starting with Kepler and K2 follow-up, we now support TESS and other space- and ground-based exoplanet related discovery and characterization projects. The importance of high-resolution imaging for exoplanet research comes via identification of nearby stellar companions that can dilute the transit signal and confound derived exoplanet and stellar parameters. Our observations therefore provide crucial information allowing accurate planet and stellar properties to be determined. Our community program obtains high-resolution imagery, reduces the data, and provides all final data products, without any exclusive use period, to the community via the Exoplanet Follow-Up Observation Program (ExoFOP) website maintained by the NASA Exoplanet Science Institute. This paper describes the need for high-resolution imaging and gives details of the speckle imaging program, highlighting some of the major scientific discoveries made along the way.

We present results of sub-arcsec ALMA observations of CO(2-1) and CO(5-4) toward a massive main sequence galaxy at z = 1.45 in the SXDS/UDS field, aiming at examining the internal distribution and properties of molecular gas in the galaxy. Our target galaxy consists of the bulge and disk, and has a UV clump in the HST images. The CO emission lines are clearly detected and the CO(5-4)/CO(2-1) flux ratio (R_52) is ~1, similar to that of the Milky Way. Assuming a metallicity dependent CO-toH_2 conversion factor and a CO(2-1)/CO(1-0) flux ratio of 2 (the Milky Way value), the molecular gas mass and the gas mass fraction (f_gas = molecular gas mass / (molecular gas mass + stellar mass)) are estimated to be ~1.5x10^11 M_Sun and ~0.55, respectively. We find that R_52 peak coincides with the position of the UV clump and its value is approximately two times higher than the galactic average. This result implies high gas density and/or high temperature in the UV clump, which qualitatively agrees with a numerical simulation of a clumpy galaxy. The CO(2-1) distribution is well represented by a rotating disk model and its half-light radius is ~2.3 kpc. Compared to the stellar distribution, the molecular gas is more concentrated in the central region of the galaxy. We also find that f_gas decreases from ~0.6 at the galactic center to ~0.2 at 3xhalf-light radius, indicating that the molecular gas is distributed in more central region of the galaxy than stars and seems to associate with the bulge rather than the stellar disk.

The modern Very Long Baseline Interferometry (VLBI) relativistic delay model, as documented in the IERS Conventions refers to the time epoch when the signal passes one of two stations of an interferometer baseline (selected arbitrarily from the pair of stations and called the 'reference station', or 'station 1'). This model consists of the previous correlation procedure used before the year 2002. However, since 2002 a new correlation procedure that produces the VLBI group delays referring to the time epoch of signal passage at the geocenter has been used. A corresponding correction to the conventional VLBI model delay has to be introduced. However, this correction has not been thoroughly presented in peer reviewed journals, and different approaches are used at the correlators to calculate the final group delays officially published in the IVS database. This may cause an inconsistency up to 6 ps for ground-based VLBI experiments between the group delay obtained by the correlator and the geometrical model delay from the IERS Conventions used in data analysis software. Moreover, a miscalculation of the signal arrival moment to the 'reference station' could result a larger modelling error (up to 50 ps). The paper presents the justification of the correction due to transition between two epochs elaborated from the Lorentz transformation, and the approach to model the uncertainty of the calculation of the signal arrival moment. The both changes are particularly essential for upcoming broadband technology geodetic VLBI observations.

The LiteBIRD satellite is planned to be launched by JAXA in the late 2020s. Its main purpose is to observe the large-scale B-mode polarization in the Cosmic Microwave Background (CMB) anticipated from the Inflation theory. LiteBIRD will observe the sky for three years at the second Lagrangian point (L2) of the Sun-Earth system. Planck was the predecessor for observing the CMB at L2, and the onboard High Frequency Instrument (HFI) suffered contamination by glitches caused by the cosmic-ray (CR) hits. We consider the CR hits can also be a serious source of the systematic uncertainty for LiteBIRD. Thus, we have started a comprehensive end-to-end simulation study to assess impact of the CR hits for the LiteBIRD detectors. Here, we describe procedures to make maps and power spectra from the simulated time-ordered data, and present initial results. Our initial estimate is that $C_l^{BB}$ by CR is $\sim 2 \times 10^{-6}~\mu$K$_{\mathrm{CMB}}^{2}$ in a one-year observation with 12 detectors assuming that the noise is 1~aW/$\sqrt{\mathrm{Hz}}$ for the differential mode of two detectors constituting a polarization pair.

We present the new results of our long-term observational project to detect the small variations in the orbital periods of low-mass and short-period eclipsing binaries. About 120 new precise mid-eclipse times were obtained for three relatively well-known dwarf eclipsing binaries: SDSS J143547.87+373338.5 (P = 0.126 d), NSVS 07826147 (0.162 d), and NSVS 14256825 (0.110 d). Observed-minus-calculated (O-C) diagrams of these systems were analyzed using all accurate timings, and, where possible, new parameters of the light-time effect were calculated. For the first time, we derive (or improve upon previous findings with regard to) the short orbital periods of 13 and 10 years of possible third bodies for SDSS J143547.87+373338.5 and NSVS 07826147, respectively. In these binaries, our data show that period variations can be modeled simply on the basis of a single circumbinary object. For the first two objects, we calculated the minimum mass of the third components to be 17 MJ, and 1.4 MJ respectively, which corresponds to the mass of a brown dwarf or a giant planet. For NSVS 14256825, the cyclical period changes caused by a single additional body cannot be confirmed by our recent eclipse time measurements. More complex behavior connected with two orbiting bodies, or yet unknown effects, should be taken into account.

Theoretical investigations predicted that high spatio-temporal resolution observations in the Sr I 4607 A line must show a conspicuous scattering polarization pattern at the solar disk-center, which encodes information on the unresolved magnetism of the inter-granular photospheric plasma. Here we present a study of the impact of limited time resolution on the observability of such forward scattering (disk-center) polarization signals. Our investigation is based on three-dimensional radiative transfer calculations in a time-dependent magneto-convection model of the quiet solar photosphere, taking into account anisotropic radiation pumping and the Hanle effect. This type of radiative transfer simulation is computationally costly, reason why the time variation had not been investigated before for this spectral line. We compare our theoretical results with recent disk-center filter polarimetric observations in the Sr I 4607 A line, showing that there is good agreement in the polarization patterns. We also show what we can expect to observe with the Visible Spectro-Polarimeter at the upcoming Daniel K. Inouye Solar Telescope.

We report the ubiquitous occurrence of nighttime temperature inversions in the tropical martian atmosphere during the dusty season, as observed by the Mars Climate Sounder. The inversions are linked to the occurrence of large-scale regional dust storms, with their strengths largely correlated to the strengths of the dust storms. Inversions strengthen between 2 am and 4 am, with the bases of the inversions getting cooler, and the tops of the inversions getting warmer. The inversions are strongest around Tharsis and Terra Sabaea, which are higher-elevation regions, suggesting they are forming due to a combination of topographically-excited tides and cloud radiative cooling. However, inversions are also observed over the flat plains, and are likely associated with stronger tides resulting from the increased dust abundance. These results highlight an important interplay between the dust distribution, water ice clouds and thermal tides.

The building of a stellar structure requires knowing the Rosseland mean opacity at each layer of the model. This mean opacity is very often interpolated in pre-computed tables due to the overwhelming time to compute it from monochromatic cross sections. The main drawback to using tables is that the opacities can be inconsistent with the actual local chemical composition, for instance in the regions of the star where nucleosynthesis occurs. We developed a strategy that allows very fast calculations of Rosseland opacities from monochromatic cross sections. This method has been implemented in the Toulouse-Geneva evolution code, which we used to compute evolutionary tracks with models whose Rosseland opacities are fully consistent with the chemical mix everywhere in the star. Our self-consistent models show very small structural differences compared to models where the Rosseland opacity is computed with a fixed metal mixture. The main-sequence evolutionary tracks are almost the same for models of mass ranging from 2 to 8 solar masses. At a given surface gravity the relative difference in age is lower than 2% and generally below 1% between the two kinds of calculations, the self-consistent model being younger most of the time. Unless such a precision in age is sought out, the use of tabulated Rosseland opacities with a metal content defined globally is still acceptable, at least in main-sequence stars where the chemical mix changes only through nucleosynthesis.

The planets HR8799bc display nearly identical colours and spectra as variable young exoplanet analogues such as VHS 1256-1257ABb and PSO J318.5-22, and are likely to be similarly variable. Here we present results from a 5-epoch SPHERE IRDIS broadband-$H$ search for variability in these two planets. HR 8799b aperture photometry and HR 8799bc negative simulated planet photometry share similar trends within uncertainties. Satellite spot lightcurves share the same trends as the planet lightcurves in the August 2018 epochs, but diverge in the October 2017 epochs. We consider $\Delta(mag)_{b} - \Delta(mag)_{c}$ to trace non-shared variations between the two planets, and rule out non-shared variability in $\Delta(mag)_{b} - \Delta(mag)_{c}$ to the 10-20$\%$ level over 4-5 hours. To quantify our sensitivity to variability, we simulate variable lightcurves by inserting and retrieving a suite of simulated planets at similar radii from the star as HR 8799bc, but offset in position angle. For HR 8799b, for periods $<$10 hours, we are sensitive to variability with amplitude $>5\%$. For HR 8799c, our sensitivity is limited to variability $>25\%$ for similar periods.

We present Yebes 40m telescope observations of the three most stable C4H3N isomers towards the cyanopolyyne peak of TMC-1. We have detected 13 transitions from CH3C3N (A and E species), 16 lines from CH2CCHCN, and 27 lines (a-type and b-type) from HCCCH2CN. We thus provide a robust confirmation of the detection of HCCCH2CN and CH2CCHCN in space. We have constructed rotational diagrams for the three species, and obtained rotational temperatures between 4-8 K and similar column densities for the three isomers, in the range (1.5-3)e12 cm-2. Our chemical model provides abundances of the order of the observed ones, although it overestimates the abundance of CH3CCCN and underestimates that of HCCCH2CN. The similarity of the observed abundances of the three isomers suggests a common origin, most probably involving reactions of the radical CN with the unsaturated hydrocarbons methyl acetylene and allene. Studies of reaction kinetics at low temperature and further observations of these molecules in different astronomical sources are needed to draw a clear picture of the chemistry of C4H3N isomers in space.

Using the state-of-the-art SKA precursor, the MeerKAT radio telescope, we explore the limits to precision pulsar timing of millisecond pulsars achievable due to pulse stochasticity (jitter). We report new jitter measurements in 15 of the 29 pulsars in our sample and find that the levels of jitter can vary dramatically between them. For some, like the 2.2~ms pulsar PSR J2241--5236, we measure an implied jitter of just $\sim$ 4~ns/hr, while others like the 3.9~ms PSR J0636--3044 are limited to $\sim$ 100 ns/hr. While it is well known that jitter plays a central role to limiting the precision measurements of arrival times for high signal-to-noise ratio observations, its role in the measurement of dispersion measure (DM) has not been reported, particularly in broad-band observations. Using the exceptional sensitivity of MeerKAT, we explored this on the bright millisecond pulsar PSR J0437--4715 by exploring the DM of literally every pulse. We found that the derived single pulse DMs vary by typically 0.0085 cm$^{-3}$ pc from the mean, and that the best DM estimate is limited by the differential pulse jitter across the band. We postulate that all millisecond pulsars will have their own limit on DM precision which can only be overcome with longer integrations. Using high-time resolution filterbank data of 9 $\mu$s, we also present a statistical analysis of single pulse phenomenology. Finally, we discuss optimization strategies for the MeerKAT pulsar timing program and its role in the context of the International Pulsar Timing Array (IPTA).

We present the compositional analysis of three terrestrial analogues of Martian olivine-bearing rocks derived from both laboratory and flight-derived analytical instruments. In the first step, state-of-the-art spectroscopic (XRF, NIR and Raman) and diffractometric (XRD) laboratory systems were complementary used. Besides providing a detailed mineralogical and geochemical characterization of the samples, results comparison shed light on the advantages ensured by the combined use of Raman and NIR techniques, being these the spectroscopic instruments that will soon deploy (2021) on Mars as part of the ExoMars/ESA rover payload. In order to extrapolate valuable indicators of the mineralogical data that could derive from the ExoMars/Raman Laser Spectrometer (RLS), laboratory results were then compared with the molecular data gathered through the RLS ExoMars Simulator. Beside correctly identifying all major phases (feldspar, pyroxene and olivine), the RLS ExoMars Simulator confirmed the presence of additional minor compounds (i.e. hematite and apatite) that were not detected by complementary techniques. Furthermore, concerning the in-depth study of olivine grains, the RLS ExoMars simulator was able to effectively detect the shifting of the characteristic double peak around 820 and 850 cm-1, from which the Fe-Mg content of the analysed crystals can be extrapolated. Considering that olivine is one of the main mineral phases of the ExoMars landing site (Oxia Planum), this study suggests that the ExoMars/RLS system has the potential to provide detailed information about the elemental composition of olivine on Mars.

In the present work, NIR, LIBS, Raman and XRD techniques have been complementarily used to carry out a comprehensive characterization of a terrestrial analogue selected from the Chesapeake Bay Impact Structure (CBIS). The obtained data clearly highlight the key role of Raman spectroscopy in the detection of minor and trace compounds, through which inferences about geological processes occurred in the CBIS can be extrapolated. Beside the use of commercial systems, further Raman analyses were performed by the Raman Laser Spectrometer (RLS) ExoMars Simulator. This instrument represents the most reliable tool to effectively predict the scientific capabilities of the ExoMars/Raman system that will be deployed on Mars in 2021. By emulating the analytical procedures and operational restrictions established by the ExoMars mission rover design, it was proved that the RLS ExoMars Simulator is able to detect the amorphization of quartz, which constitutes an analytical clue of the impact origin of craters. On the other hand, the detection of barite and siderite, compounds crystallizing under hydrothermal conditions, helps to indirectly confirm the presence of water in impact targets. Furthermore, the RLS ExoMars Simulator capability of performing smart molecular mappings was also evaluated. According to the obtained results, the algorithms developed for its operation provide a great analytical advantage over most of the automatic analysis systems employed by commercial Raman instruments, encouraging its application for many additional scientific and commercial purposes.

Galaxies can be classified as passive ellipticals or star-forming discs. Ellipticals dominate at the high end of the mass range, and therefore there must be a mechanism responsible for the quenching of star-forming galaxies. This could either be due to the secular processes linked to the mass and star formation of galaxies or to external processes linked to the surrounding environment. In this paper, we analytically model the processes that govern galaxy evolution and quantify their contribution. We have specifically studied the effects of mass quenching, gas stripping, and mergers on galaxy quenching. To achieve this, we first assumed a set of differential equations that describe the processes that shape galaxy evolution. We then modelled the parameters of these equations by maximising likelihood. These equations describe the evolution of galaxies individually, but the parameters of the equations are constrained by matching the extrapolated intermediate-redshift galaxies with the low-redshift galaxy population. In this study, we modelled the processes that change star formation and stellar mass in massive galaxies from the GAMA survey between z~0.4 and the present. We identified and quantified the contributions from mass quenching, gas stripping, and mergers to galaxy quenching. The quenching timescale is on average 1.2 Gyr and a closer look reveals support for the slow-then-rapid quenching scenario. The major merging rate of galaxies is about once per 10~Gyr, while the rate of ram pressure stripping is significantly higher. In galaxies with decreasing star formation, we show that star formation is lost to fast quenching mechanisms such as ram pressure stripping and is countered by mergers, at a rate of about 41% Gyr$^{-1}$ and to mass quenching 49% Gyr$^{-1}$. (abridged)

Classical Cepheids (DCEPs) are the most important primary indicators for the extragalactic distance scale, but they are also important objects per se, allowing us to put constraints on the physics of intermediate-mass stars and the pulsation theories. We have investigated the peculiar DCEP HD 344787, which is known to exhibit the fastest positive period change among DCEPs along with a quenching amplitude of the light variation. We have used high-resolution spectra obtained with HARPS-N@TNG for HD 344787 and the more famous Polaris DCEP, to infer their detailed chemical abundances. Results from the analysis of new time-series photometry of HD 344787 obtained by the TESS satellite are also reported. The double mode nature of HD344787 pulsation is confirmed by analysis of the TESS light curve, although with rather tiny amplitudes of a few tens of millimag. This is an indication that HD344787 is on the verge of quenching the pulsation. Analysis of the HARPS-N@TNG spectra reveals an almost solar abundance and no depletion of carbon and oxygen. Hence, the star appears to have not gone through the first dredge-up. Similar results are obtained for Polaris. Polaris and HD344787 are confirmed to be both most likely at their first crossing of the instability strip (IS). The two stars are likely at the opposite borders of the IS for first overtone DCEPs with metal abundance Z=0.008. A comparison with other DCEPs which are also thought to be at their first crossing allows us to speculate that the differences we see in the Hertzsprung-Russell diagram might be due to differences in the properties of the DCEP progenitors during the main sequence phase.

We evaluate what will be the effectiveness of the ExoMars Raman Laser Spectrometer (RLS) to determine the degree of serpentinization of olivine-rich units on Mars. We selected terrestrial analogues of martian ultramafic rocks from the Leka Ophiolite Complex (LOC) and analyzed them with both laboratory and flight-like analytical instruments. We first studied the mineralogical composition of the samples (mostly olivine and serpentine) with state-of-the-art diffractometric and spectroscopic laboratory systems. We compared these results with those obtained using our RLS ExoMars Simulator. Our work shows that the RLS ExoMars Simulator successfully identified all major phases. Moreover, when emulating the automatic operating mode of the flight instrument, the RLS ExoMars simulator also detected several minor compounds, some of which were not observed by NIR and XRD. Thereafter, we produced RLS dedicated calibration curves (R2 between 0.9993 and 0.9995 with an uncertainty between 3.0% and 5.2% with a confidence interval of 95%) to estimate the relative content of olivine and serpentine in the samples. Our results show that RLS can be very effective to identify serpentine, a scientific target of primary importance for the potential detection of biosignatures on Mars the main objective of the ExoMars rover mission.

Motivated by Hubble tension, we have recently witnessed a number of "model independent" $H_0$ determinations at cosmological scales. Here we compare two "model independent" techniques, Taylor expansion and Gaussian Processes (GP). While Taylor expansion is truly model independent in a limited range, we show that one can reduce the $H_0$ errors by increasing the range of the expansion, but the approximation suffers. For GP, we confirm for the Mat\'ern class kernels that the errors on $H_0$ decrease as the parameter $\nu \rightarrow \infty$, where we recover the Gaussian kernel. The errors from GP are typically smaller than Taylor and by mapping the GP analysis back into the Taylor expansion, we show that GP explores a smaller portion of the parameter space. In a direct comparison of GP with the CPL model, the simplest model of dynamical dark energy, we see that correlations are suppressed by GP relative to CPL. Therefore, GP cannot be model independent. We emphasise that if a truly model independent statement of Hubble tension exists, then it will have serious consequences for the FLRW framework.

Geochemical and astronomical evidence demonstrate that planet formation occurred in two spatially and temporally separated reservoirs. The origin of this dichotomy is unknown. We use numerical models to investigate how the evolution of the solar protoplanetary disk influenced the timing of protoplanet formation and their internal evolution. Migration of the water snow line can generate two distinct bursts of planetesimal formation that sample different source regions. These reservoirs evolve in divergent geophysical modes and develop distinct volatile contents, consistent with constraints from accretion chronology, thermo-chemistry, and the mass divergence of inner and outer Solar System. Our simulations suggest that the compositional fractionation and isotopic dichotomy of the Solar System was initiated by the interplay between disk dynamics, heterogeneous accretion, and internal evolution of forming protoplanets.

We present a detailed three-dimensional (3D) view of a prominence eruption, coronal loop expansion, and coronal mass ejections (CMEs) associated with an M4.4 flare that occurred on 2011 March 8 in the active region NOAA 11165. Full-disk H$\alpha$ images of the flare and filament ejection were successfully obtained by the Flare Monitoring Telescope (FMT) following its relocation to Ica University, Peru. Multiwavelength observation around the H$\alpha$ line enabled us to derive the 3D velocity field of the H$\alpha$ prominence eruption. Features in extreme ultraviolet were also obtained by the Atmospheric Imager Assembly onboard the {\it Solar Dynamic Observatory} and the Extreme Ultraviolet Imager on board the {\it Solar Terrestrial Relations Observatory - Ahead} satellite. We found that, following collision of the erupted filament with the coronal magnetic field, some coronal loops began to expand, leading to the growth of a clear CME. We also discuss the succeeding activities of CME driven by multiple interactions between the expanding loops and the surrounding coronal magnetic field.

The depletion of SO$_2$ and H$_2$O in and above the clouds of Venus (45 -- 65 km) cannot be explained by known gas-phase chemistry and the observed composition of the atmosphere. We apply a full-atmosphere model of Venus to investigate three potential explanations for the SO$_2$ and H$_2$O depletion: (1) varying the below-cloud water vapor (H$_2$O), (2) varying the below-cloud sulfur dioxide (SO$_2$), and (3) the incorporation of chemical reactions inside the sulfuric acid cloud droplets. We find that increasing the below-cloud H$_2$O to explain the SO$_2$ depletion results in a cloud top that is 20 km too high, above-cloud O$_2$ three orders of magnitude greater than observational upper limits and no SO above 80 km. The SO$_2$ depletion can be explained by decreasing the below-cloud SO$_2$ to 20 ppm. The depletion of SO$_2$ in the clouds can also be explained by the SO$_2$ dissolving into the clouds, if the droplets contain hydroxide salts. These salts buffer the cloud pH. The amount of salts sufficient to explain the SO$_2$ depletion entail a droplet pH of $\sim 1$ at 50 km. Since sulfuric acid is constantly condensing out into the cloud droplets, there must be a continuous and pervasive flux of salts of ~1e-13 mol cm$^{-2}$ s$^{-1}$ driving the cloud droplet chemistry. An atmospheric probe can test both of these explanations by measuring the pH of the cloud droplets and the concentrations of gas-phase SO$_2$ below the clouds.

We study the cosmological consequences of a class of Dirac-Born-Infeld models, and assess their viability as a candidate for the recent acceleration of the Universe. The model includes both the rolling tachyon field and the generalized Chaplygin gas models as particular limits, and phenomenologically each of these provides a possible mechanism for a deviation of the value of the dark energy equation of state from its canonical (cosmological constant) value. The field-dependent potential that is characteristic of the rolling tachyon also leads to variations of the fine-structure constant $\alpha$, implying that the model can be constrained both by standard cosmological probes and by astrophysical measurements of $\alpha$. Our analysis, using the latest available low-redshfit data and local constraints from atomic clock and weak equivalence principle experiments, shows that the two possible deviations of the dark energy equation of state are constrained to be $\log_{10}{(1+w_0)_V}<-7.85$ and $\log_{10}{(1+w_0)_C}<-0.85$, respectively for the rolling tachyon and Chaplygin components, both being at the $95.4\%$ confidence level (although the latter depends on the choice of priors, in a way that we quantify). Alternatively, the $95.4\%$ confidence level bound on the dimensionless slope of the potential is $\log_{10}{\lambda}<-5.36$. This confirms previous analyses indicating that in these models the potential needs to be extremely flat.

Using our HST/ACS observations of the recently found isolated dwarf spheroidal galaxies, we homogeneously measured their star formation histories. We determined star formation rate as a function of time, as well as age and metallicity of the stellar populations. All objects demonstrate complex star formation history, with a significant portion of stars formed 10-13 Gyr ago. Nevertheless, stars of middle ages (1-8 Gyr) are presented. In order to understand how the star formation parameters influence the evolution of dSphs, we also studied a sample of nearest dSphs in different environment: isolated (d < 2 Mpc); beyond the Local Group virial radius (but within the Local Group zero velocity sphere); and the satellites of M 31 located within the virial zone (300 kpc). Using archival HST/ACS observations, we measured their star formation histories. A comparative analysis of the parameters obtained allow us to distinguish a possible effect of the spatial segregation on the dSphs evolution scenario.

The Event Horizon Telescope (EHT) will soon provide the first high-resolution images of the Galactic Centre supermassive black hole (SMBH) candidate Sagittarius A* (Sgr A*), enabling us to probe gravity in the strong-field regime. Besides studying the accretion process in extreme environments, the obtained data and reconstructed images could be used to investigate the underlying spacetime structure. In its current configuration, the EHT is able to distinguish between a rotating Kerr black hole and a horizon-less object like a boson star. Future developments can increase the ability of the EHT to tell different spacetimes apart. We investigate the capability of an advanced EHT concept, including an orbiting space antenna, to image and distinguish different spacetimes around Sgr A*. We use GRMHD simulations of accreting compact objects (Kerr and dilaton black holes, as well as boson stars) and compute their radiative signatures via general relativistic radiative transfer calculations. To facilitate comparison with upcoming and future EHT observations we produce realistic synthetic data including the source variability, diffractive and refractive scattering while incorporating the observing array, including a space antenna. From the generated synthetic observations we dynamically reconstructed black hole shadow images using regularised Maximum Entropy methods. We employ a genetic algorithm to optimise the orbit of the space antenna with respect to improved imaging capabilities and u-v-plane coverage of the combined array (ground array and space antenna and developed a new method to probe the source variability in Fourier space. The inclusion of an orbiting space antenna improves the capability of the EHT to distinguish the spin of Kerr black holes and dilaton black holes based on reconstructed radio images and complex visibilities.

We present multispectral analysis (radio, H$\alpha$, UV/EUV, and hard X-ray) of a confined flare from 2015 March 12. This flare started within the active region NOAA 12 297 and then it expanded into a large preexisting magnetic rope embedded with a cold filament. The expansion started with several brightenings located along the rope. This process was accompanied by a group of slowly positively drifting bursts in the 0.8--2 GHz range. The frequency drift of these bursts was 45 -- 100 MHz s$^{-1}$. One of the bursts had an S-like form. During the brightening of the rope we observed an unique bright EUV structure transverse to the rope axis. The structure was observed in a broad range of temperatures and it moved along the rope with the velocity of about 240 km s$^{-1}$. When the structure dissipated, we saw a plasma further following twisted threads in the rope. The observed slowly positively drifting bursts were interpreted considering particle beams and we show that one with the S-like form could be explained by the beam propagating through the helical structure of the magnetic rope. The bright structure transverse to the rope axis was interpreted considering line-of-sight effects and the dissipation-spreading process, which we found to be more likely.

A dark-energy which behaves as the cosmological constant until a sudden phantom transition at very-low redshift ($z<0.1$) seems to solve the >4$\sigma$ disagreement between the local and high-redshift determinations of the Hubble constant, while maintaining the phenomenological success of the $\Lambda$CDM model with respect to the other observables. Here, we show that such a hockey-stick dark energy cannot solve the $H_0$ crisis. The basic reason is that the supernova absolute magnitude $M_B$ that is used to derive the local $H_0$ constraint is not compatible with the $M_B$ that is necessary to fit supernova, BAO and CMB data, and this disagreement is not solved by a sudden phantom transition at very-low redshift. Finally, we encourage the community to adopt in the statistical analyses the prior on the supernova absolute magnitude $M_B$ as an alternative to the prior on $H_0$. The three reasons are: i) one avoids double counting of low-redshift supernovae, ii) one avoids fixing the deceleration parameter to the standard model value $q_0=-0.55$, iii) one includes in the analysis the fact that $M_B$ is constrained by local calibration, an information which would otherwise be neglected in the analysis, biasing both model selection and parameter constraints.

The Cadmium Zinc Telluride Imager (CZTI) is an imaging instrument onboard AstroSat. This instrument operates as a nearly open all-sky detector above ~60 keV, making possible long integrations irrespective of the spacecraft pointing. We present a technique based on the AstroSat-CZTI data to explore the hard X-ray characteristics of the $\gamma$-ray pulsar population. We report highly significant ($\sim 30\sigma$) detection of hard X-ray (60--380 keV) pulse profile of the Crab pulsar using $\sim$5000 ks of CZTI observations within 5 to 70 degrees of Crab position in the sky, using a custom algorithm developed by us. Using Crab as our test source, we estimate the off-axis sensitivity of the instrument and establish AstroSat-CZTI as a prospective tool in investigating hard X-ray characteristics of $\gamma$-ray pulsars as faint as 10 mCrab.

We study the effects of additional cooling due to the emission of a dark matter candidate particle, the dark photon, on the final phases of the evolution of a $15\,M_\odot$ star and resulting modifications of the pre-supernova neutrino signal. For a substantial portion of the dark photon parameter space the extra cooling speeds up Si burning, which results in a reduced number of neutrinos emitted during the last day before core collapse. This reduction can be described by a systematic acceleration of the relevant timescales and the results can be estimated semi-analytically in good agreement with the numerical simulations. Outside the semi-analytic regime we find more complicated effects. In a narrow parameter range, low-mass dark photons lead to an increase of the number of emitted neutrinos because of additional shell burning episodes that delay core collapse. Furthermore, relatively strong couplings produce a thermonuclear runaway during O burning, which could result in a complete disruption of the star but requires more detailed simulations to determine the outcome. Our results show that pre-supernova neutrino signals are a potential probe of the dark photon parameter space.

OPUS (Observatoire de Paris UWS System) is a job control system that aims at facilitating the access to analysis and simulation codes through an interoperable interface. The Universal Worker System pattern v1.1 (UWS) as defined by the International Virtual Observatory Alliance (IVOA) is implemented as a REST service to control the asynchronous execution of a job on a work cluster. OPUS also follows the recent IVOA Provenance Data Model recommendation to capture and expose the provenance information of jobs and results. By following well defined standards, the tool is interoperable and jobs can be run either through a web interface, or a script, and can be integrated to existing web platforms. Current instances are used in production by several projects at the Observatoire de Paris (CTA/H.E.S.S, MASER, CompOSE).

We present a case study of using a novel method to identify red supergiant (RSG) candidates in NGC 6822, based on their 1.6 $\mu$m "H-bump". We collected 32 bands of photometric data for NGC 6822 ranging from optical to MIR. By using the theoretical spectra from MARCS, we demonstrate that there is a prominent difference around 1.6 $\mu$m ("H-bump") between low-surface-gravity (LSG) and high-surface-gravity (HSG) targets. Taking advantage of this feature, we identify efficient color-color diagrams (CCDs) of rzH and rzK to separate HSG and LSG targets from crossmatching of optical and NIR data. Moreover, synthetic photometry from ATLAS9 also give similar results. Further separating RSG candidates from the rest of the LSG candidates is done by using semi-empirical criteria on NIR CMDs and resulted in 323 RSG candidates. Meanwhile, the simulation of foreground stars from Besan\c{c}on models also indicates that our selection criteria is largely free from the contamination of Galactic giants. In addition to the "H-bump" method, we also use the traditional BVR method as a comparison and/or supplement, by applying a slightly aggressive cut to select as much as possible RSG candidates (358 targets). Furthermore, the Gaia astrometric solution is used to constrain the sample, where 181 and 193 targets were selected from the "H-bump" and BVR method, respectively. The percentages of selected targets in both methods are similar as $\sim$60\%, indicating the comparable accuracy of the two methods. In total, there are 234 RSG candidates after combining targets from both methods with 140 ($\sim$60\%) of them in common. The final RSG candidates are in the expected locations on the MIR CMDs, while the spatial distribution is also coincident with the FUV-selected star formation regions, suggesting the selection is reasonable and reliable.

Recently the International Virtual Observatory Alliance (IVOA) released a standard to structure provenance metadata, and several implementations are in development in order to capture, store, access and visualize the provenance of astronomy data products. This BoF will be focused on practical needs for provenance in astronomy. A growing number of projects express the requirement to propose FAIR data (Findable, Accessible, Interoperable and Reusable) and thus manage provenance information to ensure the quality, reliability and trustworthiness of this data. The concepts are in place, but now, applied specifications and practical tools are needed to answer concrete use cases. During this session we discussed which strategies are considered by projects (observatories or data providers) to capture provenance in their context and how a end-user might query the provenance information to enhance her/his data selection and retrieval. The objective was to identify the development of tools and formats now needed to make provenance more practical needed to increase provenance take-up in the astronomical domain.

The generation of energy in a power grid which uses Photovoltaic (PV) systems depends on the projection of shadows from moving clouds in the Troposphere. This investigation proposes an efficient method of data processing for the statistical quantification of cloud features using long-wave infrared (IR) images and Global Solar Irradiance (GSI) measurements. The IR images are obtained using a data acquisition system (DAQ) mounted on a solar tracker. We explain how to remove cyclostationary biases in GSI measurements. Seasonal trends are removed from the GSI time series, using the theoretical GSI to obtain the Clear-Sky Index (CSI) time series. We introduce an atmospheric model to remove from IR images both the effect of atmosphere scatter irradiance and the effect of the Sun's direct irradiance. Scattering is produced by water spots and dust particles on the germanium lens of the enclosure. We explain how to remove the scattering effect produced by the germanium lens attached to the DAQ enclosure window of the IR camera. An atmospheric condition model classifies the sky-conditions in four different categories: clear-sky, cumulus, stratus and nimbus. When an IR image is classified in the category of clear-sky, it is used to model the scattering effect of the germanium lens.

We present forecasted cosmological constraints from combined measurements of galaxy cluster abundances from the Simons Observatory and galaxy clustering from a DESI-like experiment on two well-studied modified gravity models, the chameleon-screened $f(R)$ Hu-Sawicki model and the nDGP braneworld Vainshtein model. A Fisher analysis is conducted using $\sigma_8$ constraints derived from thermal Sunyaev-Zel'dovich (tSZ) selected galaxy clusters, as well as linear and mildly non-linear redshift-space 2-point galaxy correlation functions. We find that the cluster abundances drive the constraints on the nDGP model while $f(R)$ constraints are led by galaxy clustering. The two tracers of the cosmological gravitational field are found to be complementary, and their combination significantly improves constraints on the $f(R)$ in particular in comparison to each individual tracer alone. For a fiducial model of $f(R)$ with $\text{log}_{10}(f_{R0})=-6$ and $n=1$ we find combined constraints of $\sigma(\text{log}_{10}(f_{R0}))=0.48$ and $\sigma(n)=2.3$, while for the nDGP model with $n_{\text{nDGP}}=1$ we find $\sigma(n_{\text{nDGP}})=0.087$. Around a fiducial General Relativity (GR) model, we find a $95\%$ confidence upper limit on $f(R)$ of $f_{R0}\leq5.68\times 10^{-7}$. Our results present the exciting potential to utilize upcoming galaxy and CMB survey data available in the near future to discern and/or constrain cosmic deviations from GR.

Two-phase xenon detectors, such as that at the core of the forthcoming LZ dark matter experiment, use photomultiplier tubes to sense the primary (S1) and secondary (S2) scintillation signals resulting from particle interactions in their liquid xenon target. This paper describes a simulation study exploring two techniques to lower the energy threshold of LZ to gain sensitivity to low-mass dark matter and astrophysical neutrinos, which will be applicable to other liquid xenon detectors. The energy threshold is determined by the number of detected S1 photons; typically, these must be recorded in three or more photomultiplier channels to avoid dark count coincidences that mimic real signals. To lower this threshold: a) we take advantage of the double photoelectron emission effect, whereby a single vacuum ultraviolet photon has a $\sim20\%$ probability of ejecting two photoelectrons from a photomultiplier tube photocathode; and b) we drop the requirement of an S1 signal altogether, and use only the ionization signal, which can be detected more efficiently. For both techniques we develop signal and background models for the nominal exposure, and explore accompanying systematic effects, including the dependence on the free electron lifetime in the liquid xenon. When incorporating double photoelectron signals, we predict a factor of $\sim 4$ sensitivity improvement to the dark matter-nucleon scattering cross-section at $2.5$ GeV/c$^2$, and a factor of $\sim1.6$ increase in the solar $^8$B neutrino detection rate. Dropping the S1 requirement may allow sensitivity gains of two orders of magnitude in both cases. Finally, we apply these techniques to even lower masses by taking into account the atomic Migdal effect; this could lower the dark matter particle mass threshold to $80$ MeV/c$^2$.

We report on SRG/eROSITA, ZTF, ASAS-SN, Las Cumbres, NEOWISE-R, and Swift XRT/UVOT observations of the unique ongoing event AT 2019avd, located in the nucleus of a previously inactive galaxy at $z=0.029$. eROSITA first observed AT 2019avd on 2020-04-28 during its first all sky survey, when it was detected as an ultra-soft X-ray source ($kT\sim 85$ eV) that was $\gtrsim 90$ times brighter in the $0.2-2$ keV band than a previous 3$\sigma$ upper flux detection limit (with no archival X-ray detection at this position). The ZTF optical light curve in the $\sim 450$ days preceding the eROSITA detection is double peaked, and the eROSITA detection coincides with the rise of the second peak. Follow-up optical spectroscopy shows the emergence of a Bowen fluorescence feature and high-ionisation coronal lines ([\ion{Fe}{X}] 6375 {\AA}, [\ion{Fe}{XIV}] 5303 {\AA}), along with persistent broad Balmer emission lines (FWHM$\sim 1400$ km s$^{-1}$). Whilst the X-ray properties make AT 2019avd a promising tidal disruption event (TDE) candidate, the optical properties are atypical for optically selected TDEs. We discuss potential alternative origins that could explain the observed properties of AT 2019avd, such as a stellar binary TDE candidate, or a TDE involving a super massive black hole binary.

We have discovered a new candidate redback millisecond pulsar binary near the center of the error ellipse of the bright unassociated Fermi-LAT $\gamma$-ray source 4FGL J0940.3-7610. The candidate counterpart is a variable optical source that also shows faint X-ray emission. Optical photometric and spectroscopic monitoring with the SOAR telescope indicates the companion is a low-mass star in a 6.5-hr orbit around an invisible primary, showing both ellipsoidal variations and irradiation and consistent with the properties of known redback millisecond pulsar binaries. Given the orbital parameters, preliminary modeling of the optical light curves suggests an edge-on inclination and a low-mass ($\sim 1.2$ - $1.4\,M_{\odot}$) neutron star, along with a secondary mass somewhat more massive than typical $\gtrsim 0.4\,M_{\odot}$. This combination of inclination and secondary properties could make radio eclipses more likely for this system, explaining its previous non-discovery in radio pulsation searches. Hence 4FGL J0940.3-7610 may be a strong candidate for a focused search for $\gamma$-ray pulsations to enable the future detection of a millisecond pulsar.

We study the black hole's shadow for Schwarzschild - de Sitter and Kerr - de Sitter metrics with the contribution of the cosmological constant \Lambda. Based on the reported parameters of the M87* black hole shadow we obtain constraints for the $\Lambda$ and show the agreement with the cosmological data. It is shown that, the coupling of the \Lambda-term with the spin parameter reveals peculiarities for the photon spheres and hence for the shadows. Within the parametrized post-Newtonian formalism the constraint for the corresponding \Lambda-determined parameter is obtained.

Numerical relativity is central to the investigation of astrophysical sources in the dynamical and strong-field gravity regime, such as binary black hole and neutron star coalescences. Current challenges set by gravitational-wave and multi-messenger astronomy call for highly performant and scalable codes on modern massively-parallel architectures. We present GR-Athena++, a general-relativistic, high-order, vertex-centered solver that extends the oct-tree, adaptive mesh refinement capabilities of the astrophysical (radiation) magnetohydrodynamics code Athena++. To simulate dynamical space-times GR-Athena++ uses the Z4c evolution scheme of numerical relativity coupled to the moving puncture gauge. We demonstrate stable and accurate binary black hole merger evolutions via extensive convergence testing, cross-code validation, and verification against state-of-the-art effective-one-body waveforms. GR-Athena++ leverages the task-based parallelism paradigm of Athena++ to achieve excellent scalability. We measure strong scaling efficiencies above $95\%$ for up to $\sim 1.2\times10^4$ CPUs and excellent weak scaling is shown up to $\sim 10^5$ CPUs in a production binary black hole setup with adaptive mesh refinement. GR-Athena++ thus allows for the robust simulation of compact binary coalescences and offers a viable path towards numerical relativity at exascale.

The axion, a hypothetical elementary pseudoscalar, is expected to solve the strong CP problem of QCD and is also a promising candidate for dark matter. The most sensitive axion search experiments operate at millikelvin temperatures and hence rely on instrumentation that carries signals from a system at cryogenic temperatures to room temperature instrumentation. One of the biggest limiting factors affecting the parameter scanning speed of these detectors is the noise added by the components in the signal detection chain. Since the first amplifier in the chain limits the minimum noise, low-noise amplification is of paramount importance. This paper reports on the operation of a flux-driven Josephson parametric amplifier (JPA) operating at around 2.3 GHz with added noise approaching the quantum limit. The JPA was employed as a first stage amplifier in an experimental setting similar to the ones used in haloscope axion detectors. By operating the JPA at a gain of 19 dB and cascading it with two cryogenic amplifiers operating at 4 K, noise temperatures as low as 120 mK were achieved for the whole signal detection chain.

General relativistic effects in the spacetime around the massive astrophysical objects can be captured using a spinning test gyro orbiting around the object in a circular geodesic. This article discusses how the tidal disruption due to a companion object affects the precession frequency of a spinning gyro orbiting around a compact astrophysical object. The precession frequency is studied in a region of space around the central object using a perturbative approach. In this study, the central object is either a neutron star or a white dwarf. The gyro is any planetary or asteroid-like object orbiting the neutron star or a white dwarf. Moreover, the companion object that causes the tidal field can be a neutron star, white dwarf, a black hole, or a main-sequence star. The tidal effect significantly affects the spacetime around the host star, which affects the gyro precession frequency. The gyro's precession frequency increases with the mass of the companion object and decreases as the separation between the host star and the companion star increases. The tidal effect also varies with the stiffness of the equation of state of matter describing the neutron star. We also find that the tidal field affects the spacetime around a white dwarf more than that of the neutron star.

In this paper, we propose an extension of the Ricci-inverse gravity, which has been proposed recently as a very novel type of fourth-order gravity, by introducing a second order term of the so-called anticurvature scalar as a correction. The main purpose of this paper is that we would like to see whether the extended Ricci-inverse gravity model admits the homogeneous and isotropic Friedmann-Lemaitre-Robertson-Walker metric as its stable inflationary solution. However, a no-go theorem for inflation in this extended Ricci-inverse gravity is shown to appear through a stability analysis based on the dynamical system method. As a result, this no-go theorem implies that it is impossible to have such stable inflation in this extended Ricci-inverse gravity model.

The composition of dark matter is one of the puzzling topics in astrophysics. Since, the existence of axions would fill this gap of knowledge, several experiments for the search of axions have been designed in the last twenty years. Among all the others, light shining through walls experiments promise to push the exclusion limits to lower energies. To this end, effort is put for the development of single-photon detectors operating at frequencies $<100$ Ghz. Here, we review recent advancements in superconducting single-photon detection. In particular, we present two sensors based on one-dimensional Josephson junctions with the capability to be in situ tuned by simple current bias: the nanoscale transition edge sensor (nano-TES) and the Josephson escape sensor (JES). These two sensors seem to be the perfect candidates for the realization of microwave light shining through walls (LSW) experiments, since they show unprecedented frequency resolutions of about 100 GHz and 2 GHz for the nano-TES and JES, respectively.

Using computer simulations we study the geometry of a typical quantum universe, i.e. the geometry one might expect before a possible period of inflation. We display it using coordinates defined by means of four classical scalar fields satisfying the Laplace equation with non-trivial boundary conditions. The field configurations reveal cosmic web structures surprisingly similar to the ones observed in the present-day Universe. Inflation might make these structures relevant for our Universe.

The Aria project consists of a plant, hosting a 350 m cryogenic isotopic distillation column, the tallest ever built, which is currently in the installation phase in a mine shaft at Carbosulcis S.p.A., Nuraxi-Figus (SU), Italy. Aria is one of the pillars of the argon dark-matter search experimental program, lead by the Global Argon Dark Matter Collaboration. Aria was designed to reduce the isotopic abundance of $^{39}$Ar, a $\beta$-emitter of cosmogenic origin, whose activity poses background and pile-up concerns in the detectors, in the argon used for the dark-matter searches, the so-called Underground Argon (UAr). In this paper, we discuss the requirements, design, construction, tests, and projected performance of the plant for the isotopic cryogenic distillation of argon. We also present the successful results of isotopic cryogenic distillation of nitrogen with a prototype plant, operating the column at total reflux.