Papers for Tuesday, Sep 28 2021

Here we show that precise Gaia EDR3 proper motions have provided robust estimates of 3D velocities, angular momentum and total energy for 40 Milky Way dwarfs. The results are statistically robust and are independent of the Milky Way mass profile. Dwarfs do not behave like long-lived satellites of the Milky Way because of their excessively large velocities, angular momenta, and total energies. Comparing them to other MW halo population, we find that many are at first passage, $\le$ 2 Gyr ago, i.e., more recently than the passage of Sagittarius, $\sim$ 4-5 Gyr ago. We suggest that this is in agreement with the stellar populations of all dwarfs, for which we find that a small fraction of young stars cannot be excluded. We also find that dwarf radial velocities contribute too little to their kinetic energy when compared to satellite systems with motions only regulated by gravity, and some other mechanism must be at work such as ram pressure. The latter may have preferentially reduced radial velocities when dwarf progenitors entered the halo until they lost their gas. It could also explain why most dwarfs lie near their pericenter. We also discover a novel large scale structure perpendicular to the Milky Way disk, which is made by 20% of dwarfs orbiting or counter orbiting with the Sagittarius dwarf.

Cosmic ray muons detected by deep underground or underwater detectors have served as an information source on the high-energy cosmic ray spectrum and hadronic interactions in air showers for almost a century. The theoretical interest in underground muons has nearly faded because space-borne experiments probe the cosmic ray spectrum more directly, and accelerators provide precise measurements of hadron yields. However, underground muons probe unique hadron interaction energies and phase space, which are still inaccessible to present accelerator experiments. The cosmic ray nucleon energies reach the hundred-TeV and PeV ranges, which are barely accessible with space-borne experiments. Our new calculation combines two modern computational tools, MCEq, for surface muon fluxes, and PROPOSAL, for underground transport. We demonstrate excellent agreement with measurements of cosmic ray muon intensities underground within estimated errors. Beyond that, the precision of historical data turns out to be significantly smaller than our error estimates. This result shows that the sources of high-energy atmospheric lepton flux uncertainties at the surface or underground can be significantly constrained without taking more data or building new detectors. The reduction of uncertainties can be expected to impact data analyses at large-volume neutrino telescopes and for the design of future ton-scale direct dark matter detectors.

Characterizing the sub-mm Galactic emission has become increasingly critical especially in identifying and removing its polarized contribution from the one emitted by the Cosmic Microwave Background (CMB). In this work, we present a parametric foreground removal performed onto sub-patches identified in the celestial sphere by means of spectral clustering. Our approach takes into account efficiently both the geometrical affinity and the similarity induced by the measurements and the accompanying errors. The optimal partition is then used to parametrically separate the Galactic emission encoding thermal dust and synchrotron from the CMB one applied on two nominal observations of forthcoming experiments from the ground and from the space. Performing the parametric fit singularly on each of the clustering derived regions results in an overall improvement: both controlling the bias and the uncertainties in the CMB $B-$mode recovered maps. We finally apply this technique using the map of the number of clouds along the line of sight, $\mathcal{N}_c$, as estimated from HI emission data and perform parametric fitting onto patches derived by clustering on this map. We show that adopting the $\mathcal{N}_c$ map as a tracer for the patches related to the thermal dust emission, results in reducing the $B-$mode residuals post-component separation. The code is made publicly available.

The distribution of stellar metallicities within and across galaxies is an excellent relic of the chemical evolution across cosmic time. We present a detailed analysis of spatially resolved stellar populations based on $>2.6$ million spatial bins from 7439 nearby galaxies in the SDSS-IV MaNGA survey. To account for accurate inclination corrections, we derive an equation for morphology dependent determination of galaxy inclinations. Our study goes beyond the well-known global mass-metallicity relation and radial metallicity gradients by providing a statistically sound exploration of local relations between stellar metallicity $[Z/H]$, stellar surface mass density $\Sigma_\star$ and galactocentric distance in the global mass-morphology plane. We find a significant resolved mass density-metallicity relation $\rm r\Sigma_\star ZR$ for galaxies of all types and masses above $10^{9.8}\,\mathrm{M_\odot}$. Different radial distances make an important contribution to the spread of the relation. Particularly, in low and intermediate mass galaxies, we find that at fixed $\Sigma_\star$ metallicity increases with radius independently of morphology. For high masses, this radial dependence is only observed in high $\Sigma_\star$ regions of spiral galaxies. This result calls for a driver of metallicity, in addition to $\Sigma_\star$ that promotes chemical enrichment in the outer parts of galaxies more strongly than in the inner parts. We discuss gas accretion, outflows, recycling and radial migration as possible scenarios.

Linking atmospheric characteristics of planets to their formation pathways is a central theme in the study of extrasolar planets. Although the 12C/13C isotope ratio shows little variation in the solar system, the atmosphere of a super-Jupiter was recently shown to be rich in 13CO, possibly as a result of dominant ice accretion beyond the CO snowline during its formation. In this paper we aim to measure the 12CO/13CO isotopologue ratio of a young, isolated brown dwarf. While the general atmospheric characteristics of young, low mass brown dwarfs are expected to be very similar to those of super-Jupiters, their formation pathways may be different, leading to distinct isotopologue ratios. We analyse archival K-band spectra of the L dwarf 2MASS J03552337+1133437 taken with NIRSPEC at the Keck telescope. A free retrieval analysis is applied to the data to determine the isotopologue ratio 12CO/13CO in its atmosphere. The isotopologue 13CO, is detected in the atmosphere through the cross-correlation method at a signal-to-noise of ~8.4. The detection significance is determined to be ~9.5 sigma using Bayesian model comparison between two retrieval models. We retrieve an isotopologue ratio 12CO/13CO of 97(+25/-18), marginally higher than the local interstellar standard. Its C/O ratio of ~0.56 is consistent with the solar value. Although only one super-Jupiter and one brown dwarf have now a measured 12CO/13CO ratio, it is intriguing that they are different, possibly hinting to distinct formation pathways. Albeit spectroscopic similarities, isolated brown dwarfs may experience a top-down formation via gravitational collapse, which resembles star formation, while giant exoplanets favorably form through core accretion, which potentially alters isotopologue ratios in their atmospheres depending on the material they accrete from protoplanetary disks. (abridged)

We measure the low-J CO line ratio R21=CO(2-1)/CO(1-0), R32=CO(3-2)/CO(2-1), and R31 = CO(3-2)/CO(1-0) using whole-disk CO maps of nearby galaxies. We draw CO(2-1) from PHANGS--ALMA, HERACLES, and follow-up IRAM surveys; CO(1-0) from COMING and the Nobeyama CO Atlas of Nearby Spiral Galaxies; and CO(3-2) from the JCMT NGLS and APEX LASMA mapping. Altogether this yields 76, 47, and 29 maps of R21, R32, and R31 at 20" \sim 1.3 kpc resolution, covering 43, 34, and 20 galaxies. Disk galaxies with high stellar mass, log10 M_* [Msun]=10.25-11 and star formation rate, SFR=1-5 Msun/yr, dominate the sample. We find galaxy-integrated mean values and 16%-84% range of R21 = 0.65 (0.50-0.83), R32=0.50 (0.23-0.59), and R31=0.31 (0.20-0.42). We identify weak trends relating galaxy-integrated line ratios to properties expected to correlate with excitation, including SFR/M_* and SFR/L_CO. Within galaxies, we measure central enhancements with respect to the galaxy-averaged value of \sim 0.18^{+0.09}_{-0.14} dex for R21, 0.27^{+0.13}_{-0.15} dex for R31, and 0.08^{+0.11}_{-0.09} dex for R32. All three line ratios anti-correlate with galactocentric radius and positively correlate with the local star formation rate surface density and specific star formation rate, and we provide approximate fits to these relations. The observed ratios can be reasonably reproduced by models with low temperature, moderate opacity, and moderate densities, in good agreement with expectations for the cold ISM. Because the line ratios are expected to anti-correlate with the CO(1-0)-to-H_2 conversion factor, alphaCO^(1-0), these results have general implications for the interpretation of CO emission from galaxies.

We explore the stellar content of the Javalambre Photometric Local Universe Survey (J-PLUS) Data Release 2 and show its potential to identify low-metallicity stars using the Stellar Parameters Estimation based on Ensemble Methods (SPEEM) pipeline. SPEEM is a tool to provide determinations of atmospheric parameters for stars and separate stellar sources from quasars, using the unique J-PLUS photometric system. The adoption of adequate selection criteria allows the identification of metal-poor star candidates suitable for spectroscopic follow-up. SPEEM consists of a series of machine learning models which uses a training sample observed by both J-PLUS and the SEGUE spectroscopic survey. The training sample has temperatures Teff between 4\,800 K and 9\,000 K; $\log g$ between 1.0 and 4.5, and $-3.1<[Fe/H]<+0.5$. The performance of the pipeline has been tested with a sample of stars observed by the LAMOST survey within the same parameter range. The average differences between the parameters of a sample of stars observed with SEGUE and J-PLUS, which were obtained with the SEGUE Stellar Parameter Pipeline and SPEEM, respectively, are $\Delta Teff\sim 41$ K, $\Delta \log g\sim 0.11$ dex, and $\Delta [Fe/H]\sim 0.09$ dex. A sample of 177 stars have been identified as new candidates with $[Fe/H]<-2.5$ and 11 of them have been observed with the ISIS spectrograph at the William Herschel Telescope. The spectroscopic analysis confirms that $64\%$ of stars have $[Fe/H]<-2.5$, including one new star with $[Fe/H]<-3.0$. SPEEM in combination with the J-PLUS filter system has shown the potential to estimate the stellar atmospheric parameters (Teff, $\log g$, and [Fe/H]). The spectroscopic validation of the candidates shows that SPEEM yields a success rate of $64\%$ on the identification of very metal-poor star candidates with $[Fe/H]<-2.5$.

This work provides the first full set of vibrational and rotational spectral data needed to aid in the detection of AlH$_3$OH$_2$, SiH$_3$OH, and SiH$_3$NH$_2$ in astrophysical or simulated laboratory environments through the use of quantum chemical computations at the CCSD(T)-F12b level of theory employing quartic force fields for the three molecules of interest. Previous work has shown that SiH$_3$OH and SiH$_3$NH$_2$ contain some of the strongest bonds of the most abundant elements in space. AlH$_3$OH$_2$ also contains highly abundant atoms and represents an intermediate along the reaction pathway from H$_2$O and AlH$_3$ to AlH$_2$OH. All three of these molecules are also polar with AlH$_3$OH$_2$ having the largest dipole of 4.58 D and the other two having dipole moments in the 1.10-1.30 D range, large enough to allow for the detection of these molecules in space through rotational spectroscopy. The molecules also have substantial infrared intensities with many of the frequencies being over 90 km mol$^{-1}$ and falling within the currently uncertain 12-17 $\mu$m region of the spectrum. The most intense frequency for AlH$_3$OH$_2$ is $\nu_9$ which has an intensity of 412 km mol$^{-1}$ at 777.0 cm$^{-1}$ (12.87 $\mu$m). SiH$-3$OH has an intensity of 183 km mol$^{-1}$ at 1007.8 cm$^{-1}$ (9.92 $\mu$m) for $\nu_5$, and SiH$_3$NH$_2$ has an intensity of 215 km mol$^{-1}$ at 1000.0 cm$^{-1}$ (10.00 $\mu$m) for $\nu_7$.

The population of satellite galaxies in a host galaxy is a combination of the cumulative accretion of subhaloes and their associated star formation efficiencies, therefore, the luminosity distribution of satellites provides valuable information of both dark matter properties and star formation physics. Recently, the luminosity function of satellites in nearby Milky Way-mass galaxies has been well measured to satellites as faint as Leo I with $M_{V} \sim -8$. In addition to the finding of the diversity in the satellite luminosity functions, it has been noticed that there is a big gap among the magnitude of satellites in some host galaxies, such as M101, where the gap is around 5 in magnitude, noticeably larger than the prediction from the halo abundance matching method. The reason of this gap is still unknown. In this paper, we use a semi-analytical model of galaxy formation, combined with high-resolution N-body simulation, to investigate the probability and origin of such big gap in M101-alike galaxies. We found that, although M101 analogues are very rare with probability of \sim 0.1%-0.2% in the local universe, their formation is a natural outcome of the CDM model. The gap in magnitude is mainly due to the mass of the accreted subhaloes, not from the stochastic star formation in them. We also found that the gap is correlated with the total satellite mass and host halo mass. By tracing the formation history of M101 type galaxies, we find that they likely formed after $z \sim 1$ due to the newly accreted bright satellites. The gap is not in a stable state, and it will disappear in ~7 Gyr due to mergers of bright satellites with the central galaxy.

In the context of Cosmic Microwave Background data analysis, we study the solution to the equation that transforms scanning data into a map. As originally suggested in "messenger" methods for solving linear systems, we split the noise covariance into uniform and non-uniform parts and adjust their relative weights during the iterative solution. With simulations, we study mock instrumental data with different noise properties, and find that this "cooling" or perturbative approach is particularly effective when there is significant low-frequency noise in the timestream. In such cases, a conjugate gradient algorithm applied to this modified system converges faster and to a higher fidelity solution than the standard conjugate gradient approach. We give an analytic estimate for the parameter that controls how gradually the linear system should change during the course of the solution.

In this work, we investigate upper limits on the global escape fraction of ionizing photons ($f_{\rm esc/global}^{\rm abs}$) from a sample of galaxies probed for Lyman-continuum (LyC) emission characterized as non-LyC and LyC leakers. We present a sample of 9 clean non-contaminated (by low redshift interlopers, CCD problems and internal reflections of the instrument) galaxies which do not show significant ($>$ $3\sigma$) LyC flux between 880\AA\ $<\lambda_{rest}<$ 910\AA. The 9 galaxy stacked spectrum reveals no significant LyC flux with an upper limit of $f_{\rm esc}^{\rm abs} \leq 0.06$. In the next step of our analysis, we join all estimates of $f_{\rm esc}^{\rm abs}$ upper limits derived from different samples of $2\lesssim z < 6$ galaxies from the literature reported in last $\sim$20 years and include the sample presented in this work. We find the $f_{\rm esc}^{\rm abs}$ upper limit $\leq$ 0.084 for the galaxies recognized as non-LyC leakers. After including all known detections from literature $f_{\rm esc/global}^{\rm abs}$ upper limit $\leq$ 0.088 for all galaxies examined for LyC flux. Furthermore, $f_{\rm esc}^{\rm abs}$ upper limits for different groups of galaxies indicate that the strongest LyC emitters could be galaxies classified as Lyman alpha emitters. We also discuss the possible existence of a correlation among the observed flux density ratio $(F_{\nu}^{LyC}/F_{\nu}^{UV})_{\rm obs}$ and Lyman alpha equivalent width EW(Ly$\alpha)$, where we confirm the existence of moderately significant correlation among galaxies classified as non-LyC leakers.

The Galileo Near Infrared Mapping Spectrometer (NIMS) collected spectra of Europa in the 0.7-5.2 $\mu$m wavelength region, which have been critical to improving our understanding of the surface composition of this moon. However, most of the work done to get constraints on abundances of species like water ice, hydrated sulfuric acid, hydrated salts and oxides have used proxy methods, such as absorption strength of spectral features or fitting a linear mixture of laboratory generated spectra. Such techniques neglect the effect of parameters degenerate with the abundances, such as the average grain-size of particles, or the porosity of the regolith. In this work we revisit three Galileo NIMS spectra, collected from observations of the trailing hemisphere of Europa, and use a Bayesian inference framework, with the Hapke reflectance model, to reassess Europa's surface composition. Our framework has several quantitative improvements relative to prior analyses: (1) simultaneous inclusion of amorphous and crystalline water ice, sulfuric-acid-octahydrate (SAO), CO$_2$, and SO$_2$; (2) physical parameters like regolith porosity and radiation-induced band-center shift; and (3) tools to quantify confidence in the presence of each species included in the model, constrain their parameters, and explore solution degeneracies. We find that SAO strongly dominates the composition in the spectra considered in this study, while both forms of water ice are detected at varying confidence levels. We find no evidence of either CO$_2$ or SO$_2$ in any of the spectra; we further show through a theoretical analysis that it is highly unlikely that these species are detectable in any 1-2.5 $\mu$m Galileo NIMS data.

All circumbinary planets currently detected are in orbits that are almost coplanar to the binary orbit. While misaligned circumbinary planets are more difficult to detect, observations of polar aligned circumbinary gas and debris disks around eccentric binaries suggest that polar planet formation may be possible. A polar aligned planet has a stable orbit that is inclined by 90 degrees to the orbital plane of the binary with an angular momentum vector that is aligned to the binary eccentricity vector. With n- body simulations we model polar terrestrial planet formation using hydrodynamic gas disk simulations to motivate the initial particle distribution. Terrestrial planet formation around an eccentric binary is more likely in a polar alignment than in a coplanar alignment. Similar planetary systems form in a polar alignment around an eccentric binary and a coplanar alignment around a circular binary. The polar planetary systems are stable even with the effects of general relativity. Planetary orbits around an eccentric binary exhibit tilt and eccentricity oscillations at all inclinations, however, the oscillations are larger in the coplanar case than the polar case. We suggest that polar aligned terrestrial planets will be found in the future.

Direct imaging of widely separated exoplanets from space will obtain their reflected light spectra and measure their atmospheric properties, and small and temperate planets will be the focus of the next generation telescopes. In this work, we used our Bayesian retrieval algorithm ExoReL$^{\Re}$ to determine the constraints on the atmospheric properties of sub-Neptune planets from observations taken with a HabEx-like telescope. Small and temperate planets may have a non-H$_2$-dominated atmosphere, and therefore, we introduced the compositional analysis technique in our framework to explore the bulk atmospheric chemistry composition without any prior knowledge about it. We have developed a novel set of prior functions for the compositional analysis free parameters. We compared the performances of the framework with the flat prior and the novel prior and we reported a better performance when using the novel priors set. We found that the retrieval algorithm can not only identify the dominant gas of the atmosphere but also to constraint other less abundant gases with high statistical confidence without any prior information on the composition. The results presented here demonstrates that reflected light spectroscopy can characterize small exoplanets with diverse atmospheric composition. The Bayesian framework should be applied to design the instrument and the observation plan of exoplanet direct imaging experiments in the future.

The recent Parker Solar Probe (PSP) observations of type III radio bursts show that the effects of finite background magnetic field can be an important factor in the interpretation of data. In the present paper, the effects of background magnetic field on the plasma emission process, which is believed to be the main emission mechanism for solar coronal and interplanetary type III radio bursts, are investigated by means of the particle-in-cell simulation method. The effects of ambient magnetic field are systematically surveyed by varying the ratio of plasma frequency to electron gyro-frequency. The present study shows that for a sufficiently strong ambient magnetic field, the wave-particle interaction processes lead to a highly field-aligned longitudinal mode excitation and anisotropic electron velocity distribution function, accompanied by a significantly enhanced plasma emission at the second harmonic plasma frequency. For such a case, the polarization of the harmonic emission is almost entirely in the sense of extraordinary mode. On the other hand, for moderate strengths of the ambient magnetic field, the interpretation of the simulation result is less than clear. The underlying nonlinear mode coupling processes indicate that to properly understand and interpret the simulation results require sophisticated analyses involving interactions among magnetized plasma normal modes including the two transverse modes of the magneto-active plasma, namely, extraordinary and ordinary modes, as well as electron-cyclotron-whistler, plasma oscillation, and upper-hybrid modes. At present, a nonlinear theory suitable for quantitatively analyzing such complex mode-coupling processes in magnetized plasmas is incomplete, which calls for further theoretical research, but the present simulation results could provide a guide for future theoretical efforts.

A small variation of the circular shape of the hodograph theorem states that for every elliptical solution of the two-body problem, it is possible to find an appropriate inertial frame such that the speed of the bodies is constant. We use this result and data from the NASA JPL Horizon Web Interface to find the best fitting ellipse for the trajectory of Mercury, Venus, Earth, Mars, and Jupiter. The process requires us to find procedures to obtain the plane and ellipse that best fit a collection of points in space. We show that if we aim for the plane that minimizes the sum of the square distances from the given points to the unknown plane, we obtain three planes that appear to divide the set of points equally into octants, one of these being our desired plane of best fit. We provide a detailed proof of the hodograph theorem.

The Ice Chamber for Astrophysics-Astrochemistry (ICA) is a new laboratory end-station located at the Institute for Nuclear Research (Atomki) in Debrecen, Hungary. The ICA has been specifically designed for the study of the physico-chemical properties of astrophysical ice analogues and their chemical evolution when subjected to ionising radiation and thermal processing. The ICA is an ultra-high vacuum compatible chamber containing a series of IR-transparent substrates mounted in a copper holder connected to a closed-cycle cryostat capable of being cooled down to 20 K, itself mounted on a 360{\deg} rotation stage and a z-linear manipulator. Ices are deposited onto the substrates via background deposition of dosed gases. Ice structure and chemical composition are monitored by means of FTIR absorbance spectroscopy in transmission mode, although use of reflectance mode is possible by using metallic substrates. Pre-prepared ices may be processed in a variety of ways. A 2 MV Tandetron accelerator is capable of delivering a wide variety of high-energy ions into the ICA, which simulates ice processing by cosmic rays, the solar wind, or magnetospheric ions. The ICA is also equipped with an electron gun which may be used for electron impact radiolysis of ices. Thermal processing of both deposited and processed ices may be monitored by means of both FTIR spectroscopy and quadrupole mass spectrometry. In this paper, we provide a detailed description of the ICA set-up, as well as an overview of preliminary results obtained and future plans.

The modelling of molecular excitation and dissociation processes relevant to astrochemistry requires the validation of theories by comparison with data generated from laboratory experimentation. The newly commissioned Ice Chamber for Astrophysics-Astrochemistry (ICA) allows for the study of astrophysical ice analogues and their evolution when subjected to energetic processing, thus simulating the processes and alterations interstellar icy grain mantles and icy outer Solar System bodies undergo. ICA is an ultra-high vacuum compatible chamber containing a series of IR-transparent substrates upon which the ice analogues may be deposited at temperatures of down to 20 K. Processing of the ices may be performed in one of three ways: (i) ion impacts with projectiles delivered by a 2 MV Tandetron-type accelerator, (ii) electron irradiation from a gun fitted directly to the chamber, and (iii) thermal processing across a temperature range of 20-300 K. The physico-chemical evolution of the ices is studied in situ using FTIR absorbance spectroscopy and quadrupole mass spectrometry. In this paper, we present an overview of the ICA facility with a focus on characterising the electron beams used for electron impact studies, as well as reporting the preliminary results obtained during electron irradiation and thermal processing of selected ices.

On 4 March 2021 at 9 UTC a 30-m in diameter near-Earth asteroid 2021 DW1 passed the Earth at a distance of 570000 km, reaching the maximum brightness of V=14.6 mag. We observed it photometrically from 2 March, when it was visible at V=16.5 mag, until 7 March (V=18.2 mag). During that time 2021 DW1 swept a 170 degrees long arc in the northern sky, spanning solar phase angles in the range from 36 to 86 degrees. This made it an excellent target for physical characterisation, including spin axis and shape derivation. Convex inversion of the asteroid lightcurves gives a sidereal period of rotation P=0.013760 +/- 0.000001 h, and two solutions for the spin axis ecliptic coordinates: (A) lambda_1=57 +/- 10, beta_1=29 +/- 10, and (B) lambda_2=67 +/- 10, beta_2=-40 +/- 10. The magnitude-phase curve can be fitted with a standard H, G function with H=24.8 +/- 0.5 mag and an assumed G=0.24. The asteroid colour indices are g-i=0.79 +/- 0.01 mag, and i-z=0.01 +/- 0.02 mag which indicates an S taxonomic class, with an average geometric albedo p_V=0.23 +/- 0.02. The asteroid effective diameter, derived from H and p_V, is D=30 +/- 10 m. It was found that the inclination of the spin axis of 2021 DW1 is not perpendicular to the orbital plane (obliquity epsilon=54 +/- 10 or epsilon=123 +/- 10). More spin axes of VSAs should be determined to check, if 2021 DW1 is an exception or a typical case.

Magnetic fields in several astrophysical objects are amplified and maintained by a dynamo mechanism, which is the conversion of the turbulent kinetic energy to magnetic energy. A dynamo that amplifies magnetic fields at scales $<$ the driving scale of turbulence is known as the fluctuation dynamo. We study the properties of the fluctuation dynamo in supersonic turbulent plasmas, which is of relevance to the ISM, structure formation, and lab experiments of laser-plasma turbulence. Using simulations, we explore the properties of the exponentially growing and saturated state of the fluctuation dynamo for subsonic and supersonic turbulence. We confirm that the fluctuation dynamo efficiency decreases with compressibility. We show that the fluctuation dynamo generated magnetic fields are spatially intermittent and the level of intermittency decreases as the field saturates. We find a stronger back reaction of the magnetic field on the velocity for the subsonic case as compared to the supersonic case. Locally, we find that the level of alignment between vorticity and velocity, velocity and magnetic field, and current density and magnetic field in the saturated stage is enhanced in comparison to the exponentially growing phase for the subsonic case, but only the current density and magnetic field alignment is enhanced for the supersonic case. We show that both the magnetic field amplification (due to weaker stretching of field lines) and diffusion decreases when the field saturates, but the diffusion is enhanced relative to amplification. This occurs throughout the volume in the subsonic turbulence, but primarily in the strong-field regions for the supersonic case. This leads to the saturation of the fluctuation dynamo. Overall, both the amplification and diffusion of magnetic fields are affected and thus a drastic change in either of them is not required for the saturation. [Abstract abridged]

We numerically investigate non-Gaussianities in the late-time cosmological density field in Fourier space. We explore various statistics, including the two-point and three-point probability distribution function (PDF) of phase and modulus, and two \& three-point correlation function of of phase and modulus. We detect significant non-Gaussianity for certain configurations. We compare the simulation results with the theoretical expansion series of \citet{2007ApJS..170....1M}. We find that the $\mathcal{O}(V^{-1/2})$ order term alone is sufficiently accurate to describe all the measured non-Gaussianities in not only the PDFs, but also the correlations. We also numerically find that the phase-modulus cross-correlation contributes $\sim 50\%$ to the bispectrum, further verifying the accuracy of the $\mathcal{O}(V^{-1/2})$ order prediction. This work demonstrates that non-Gaussianity of the cosmic density field is simpler in Fourier space, and may facilitate the data analysis in the era of precision cosmology.

Axion is one of the most popular candidates of the cosmological dark matter. Recent studies considering the misalignment production of axions suggest some benchmark axion mass ranges near $m_a \sim 20$ $\mu$eV. For such axion mass, the spontaneous decay of axions can give photons in radio band frequency $\nu \sim 1-3$ GHz, which can be detected by radio telescopes. In this article, we show that using radio data of galaxy clusters would be excellent to constrain axion dark matter. Specifically, by using radio data of the Bullet cluster (1E 0657-55.8), we find that the upper limit of the axion-photon coupling constant can be constrained to $g_{a \gamma \gamma} \sim 10^{-12}-10^{-11}$ GeV$^{-1}$ for $m_a \sim 20$ $\mu$eV, which is tighter than the limit obtained by the CERN Axion Solar Telescope (CAST).

We provide a method for estimating the projected density distribution $\bar{n}_2w_p(r_p)$ of photometric objects around spectroscopic objects in a redshift survey. This quantity describes the distribution of Photometric sources with certain physical properties (e.g. luminosity, mass, color etc) Around Cosmic webs (PAC) traced by the spectroscopic objects. The method can make full use of current and future deep and wide photometric surveys to explore the formation of galaxies up to medium redshift ($z_s < 2$), with the aid of cosmological redshift surveys that sample only a fairly limited species of objects (e.g. Emission Line Galaxies). As an example, we apply the PAC method to the CMASS spectroscopic and HSC-SSP PDR2 photometric samples to explore the distribution of galaxies for a wide range of stellar mass from $10^{9.0}{\rm M_\odot}$ to $10^{12.0}{\rm M_\odot}$ around massive ones at $z_s\approx 0.6$. Using the abundance matching method, we model $\bar{n}_2w_p(r_p)$ in N-body simulation using MCMC sampling, and accurately measure the stellar-halo mass relation (SHMR) and stellar mass function (SMF) for the whole mass range. We can also measure the conditional stellar mass function (CSMF) of satellites for central galaxies of different mass. The PAC method has many potential applications for studying the evolution of galaxies.

We report on the multi-frequency multi-epoch radio observations of the magnetar, XTE J1810-197, which exhibited a radio outburst from December 2018 after its 10-year quiescent period. We performed quasi-simultaneous observations with VERA (22 GHz), Hitachi (6.9 GHz and 8.4 GHz), Kashima (2.3 GHz), and Iitate (0.3 GHz) radio telescopes located in Japan to trace the variability of the magnetar radio pulsations during the observing period from 13 December 2018 to 12 June 2019. The pulse width goes narrower as the observing frequency goes higher, analogous to the general profile narrowing behavior of ordinary pulsars. When assuming a simple power law in the range of 2.3 GHz and 8.7 GHz, the radio spectrum of the magnetar goes steeper with the average spectral index $ \langle \alpha \rangle \approx -0.85$ for the first four months. The wide-band radio spectra inferred from our observations and the literature suggest that XTE J1810-197 would have a double-peaked spectrum with a valley point in 22 - 150 GHz, where the first spectral peak infers a gigahertz-peaked spectrum (GPS) feature with a peak at a few GHz. The GPS and the high-frequency peak have been identified in the spectra of other radio-loud magnetars, thus they may be intrinsic features that can give a new insight to understand various emission mechanisms and surrounding environments of radio magnetars. Our study emphasizes the importance of simultaneous long-term broad-band observations toward radio-loud magnetars to capture the puzzling spectral features and establish a link to other types of neutron stars.

We report results of a spectropolarimetric and photometric monitoring of the weak-line T Tauri star V410 Tau based on data collected mostly with SPIRou, the near-infrared (NIR) spectropolarimeter recently installed at the Canada-France-Hawaii Telescope, as part of the SPIRou Legacy Survey large programme, and with TESS between October and December 2019. Using Zeeman-Doppler Imaging (ZDI), we obtained the first maps of photospheric brightness and large-scale magnetic field at the surface of this young star derived from NIR spectropolarimetric data. For the first time, ZDI is also simultaneously applied to high-resolution spectropolarimetric data and very-high-precision photometry. V410 Tau hosts both dark and bright surface features and magnetic regions similar to those previously imaged with ZDI from optical data, except for the absence of a prominent dark polar spot. The brightness distribution is significantly less contrasted than its optical equivalent, as expected from the difference in wavelength. The large-scale magnetic field (~410 G), found to be mainly poloidal, features a dipole of ~390 G, again compatible with previous studies at optical wavelengths. NIR data yield a surface differential rotation slightly weaker than that estimated in the optical at previous epochs. Finally, we measured the radial velocity of the star and filtered out the stellar activity jitter using both ZDI and Gaussian Process Regression down to a precision of ~0.15 and 0.08 $\mathrm{km\,s^{-1}}$ RMS, respectively, confirming the previously published upper limit on the mass of a potential close-in massive planet around V410 Tau.

We have calculated the stellar $\beta$-decay rate of the important s-process branching point ${}^{134}$Cs based on the state of the art shell model calculations. At typical $s$-process temperatures ($T\sim$ 0.2-0.3 GK), our new rate is one order of magnitude lower than the widely-used rate from Takahashi and Yokoi (hereafter TY87). The impact on the nucleosynthesis in AGB stars is investigated with various masses and metallicities. Our new decay rate leads to an overall decrease in the ${}^{134}$Ba/${}^{136}$Ba ratio, and well explains the measured ratio in meteorities without introducing the $i$ process. We also derive the elapsed time from the last AGB nucleosynthetic event that polluted the early Solar System to be $>$28 Myr based on the ${}^{135}$Cs/${}^{133}$Cs ratio, which is consistent with the elapsed times derived from ${}^{107}$Pd and ${}^{182}$Hf. The $s$-process abundance sum of ${}^{135}$Ba and ${}^{135}$Cs is found to increase, resulting in a smaller $r$-process contribution of ${}^{135}$Ba in the Solar System.

We present the first comparison between properties of clusters of galaxies detected by the eROSITA Final Equatorial-Depth Survey (eFEDS) and the XMM-Newton X-ray telescope. We have compared, in an ensemble fashion, properties from the eFEDS X-ray catalogue (542 candidates) with those from the Ultimate XMM eXtragaLactic (or XXL) survey project (100 clusters). We find the distributions of redshift and $T_{\rm X}$ to be broadly similar between the two surveys, with a larger proportion of clusters above 4 keV in the XXL-100-GC sample. However, the fractional $\Delta T_{\rm X}$ values are significantly larger in eFEDS compared to XXL. We construct a sample of 37 clusters that have eFEDS and XMM data, and we calculate minimum contamination fractions of ${\sim}$18% and ${\sim}$1% in the eFEDS X-ray and optically confirmed samples respectively, in agreement with eFEDS findings. We then make direct comparisons of properties measured by eFEDS and by XMM for individual clusters using tools developed for the XMM Cluster Survey (XCS). We compare 8 clusters that have $T_{\rm X}$ measured by both telescopes, and find that XMM temperatures are (on average) 25% larger than the corresponding eROSITA temperatures. We compare 29 X-ray luminosities ($L_{\rm X}$) measured by eFEDS and XCS and find them to be in excellent agreement. We also construct $L_{\rm X}$ - $T_{\rm X}$ scaling relations based on eFEDS and XCS measurements, which are in tension; the tension is decreased when we measure a third scaling relation with calibrated XCS temperatures.

These protoplanetary disks in T Tauri stars play a central role in star and planet formation. We spatially resolve at sub-au scales the innermost regions of a sample of T Tauri's disks to better understand their morphology and composition. We extended our homogeneous data set of 27 Herbig stars and collected near-IR K-band observations of 17 T Tauri stars, spanning effective temperatures and luminosities in the ranges of ~4000-6000 K and ~0.4-10 Lsun. We focus on the continuum emission and develop semi-physical geometrical models to fit the interferometric data and search for trends between the properties of the disk and the central star. The best-fit models of the disk's inner rim correspond to wide rings. We extend the Radius-luminosity relation toward the smallest luminosities (0.4-10 Lsun) and find the R~L^(1/2) trend is no longer valid, since the K-band sizes measured with GRAVITY are larger than the predicted sizes from sublimation radius computation. No clear correlation between the K-band half-flux radius and the mass accretion rate is seen. Having magnetic truncation radii in agreement with the K-band GRAVITY sizes would require magnetic fields as strong as a few kG, which should have been detected, suggesting that accretion is not the main process governing the location of the half-flux radius of the inner dusty disk. Our measurements agree with models that take into account the scattered light. The N-to-K band size ratio may be a proxy for disentangling disks with silicate features in emission from disks with weak and/or in absorption silicate features. When comparing inclinations and PA of the inner disks to those of the outer disks (ALMA) in nine objects of our sample, we detect misalignments for four objects.

We study the intermediate inflation in the mimetic Dirac-Born-Infeld model. By considering the scale factor as $a=a_{0}\exp(bt^{\beta})$, we show that in some ranges of the intermediate parameters $b$ and $\beta$, the model is free of the ghost and gradient instabilities. We study the scalar spectral index, tensor spectral index, and the tensor-to-scalar ratio in this model and compare the results with Planck2018 TT, TE, EE+lowE+lensing +BAO +BK14 data at $68\%$ and $95\%$ CL. In this regard, we find some constraints on the intermediate parameters that lead to the observationally viable values of the perturbation parameters. We also seek the non-gaussian features of the primordial perturbations in the equilateral configuration. By performing the numerical analysis on the nonlinearity parameter in this configuration, we show that the amplitude of the non-gaussianity in the intermediate mimetic DBI model is predicted to be in the range $-16.7<f^{equil}<-12.5$. We show that, with $0<b\leq 10$ and $0.345<\beta<0.387$, we have an instabilities-free intermediate mimetic DBI model that gives the observationally viable perturbation and non-gaussianity parameters.

The asymmetry in the cosmic microwave background (CMB) towards several nearby galaxies detected by Planck data is probably due to the rotation of "cold gas" clouds present in the galactic halos. In 1995 it had been proposed that galactic halos are populated by pure molecular hydrogen clouds which are in equillibrium with the CMB. More recently, it was shown that the equillibrium could be stable. Nevertheless, the cloud chemical composition is still a matter to be studied. To investigate this issue we need to trace the evolution of these virial cloud from the time of their formation to the present, and to confront the model with the observational data. The present paper is a short summary of a paper [1]. Here we only concentrate on the evolution of these clouds from the last scattering surface (LSS) up to the formation of first generation of stars (population-III stars).

We revisit the differential energy distribution of steady-state dynamical models. It has been shown that the differential energy distribution of steady-state spherical models does not vary strongly with the anisotropy profile, and that it is hence mainly determined by the density distribution of the model. We explore this similarity in more detail. Through a worked example and a simple proof, we show that the mean binding energy per unit mass $\langle{\cal{E}}\rangle$, or equivalently the total integrated binding energy $B_{\text{tot}} = M\langle{\cal{E}}\rangle$, is independent of the orbital structure, not only for spherical models but for any steady-state dynamical model. Only the higher-order moments of the differential energy distribution depend on the details of the orbital structure. We show that the standard deviation of the differential energy distribution of spherical dynamical models varies systematically with the anisotropy profile: radially anisotropic models tend to prefer more average binding energies, whereas models with a more tangential orbital distribution slightly favour more extreme binding energies. Finally, we find that the total integrated binding energy supplements the well-known trio consisting of total kinetic energy, total potential energy, and total energy on an equal footing. Knowledge of any one out of these four energies suffices to calculate the other three.

The halo model formalism is widely adopted in cosmological studies for predicting the growth of large-scale structure in the Universe. However, to date there have been relatively few direct comparisons of the halo model with more accurate (but much more computationally expensive) cosmological simulations. We test the accuracy of the halo model in reproducing the non-linear matter power spectrum, P(k), when the main inputs of the halo model (specifically the matter density profiles, halo mass function, and linear bias) are taken directly from the BAHAMAS simulations and we assess how well the halo model reproduces P(k) from the same simulations. We show that the halo model generally reproduces P(k) in the deep non-linear regime (1-halo) to typically a few percent accuracy, but struggles to reproduce (approx. 15% error) P(k) at intermediate scales of 0.1 < k [h/ Mpc] < 3 at z=0, marking the transition between the 1-halo and 2-halo terms. We show that the magnitude of this error is a strong function of the halo mass definition (through its effects on radial extent of haloes) and of redshift. Furthermore, we test the accuracy of the halo model in recovering the relative impact of baryons on P(k). We show that the systematic errors in recovering the absolute P(k) largely cancel when considering the relative impact of baryons. This suggests that the halo model can make precise predictions for the baryonic suppression, offering a fast and accurate way to adjust collisionless matter power spectra for the presence of baryons and associated processes.

Most solar flares demonstrate a prolonged, hourlong post-flare (or gradual) phase, characterized by arcade-like, post-flare loops (PFLs) visible in many extreme ultraviolet (EUV) passbands. These coronal loops are filled with hot -- $\sim 30 \,\mathrm{MK}$ -- and dense plasma, evaporated from the chromosphere during the impulsive phase of the flare, and they very gradually recover to normal coronal density and temperature conditions. During this gradual cooling down to $\sim 1 \,\mathrm{MK}$ regimes, much cooler -- $\sim 0.01 \,\mathrm{MK}$ -- and denser coronal rain is frequently observed inside PFLs. Understanding PFL dynamics in this long-duration, gradual phase is crucial to the entire corona-chromosphere mass and energy cycle. Here we report a simulation in which a solar flare evolves from pre-flare, over impulsive phase all the way into its gradual phase, which successfully reproduces post-flare coronal rain. This rain results from catastrophic cooling caused by thermal instability, and we analyse the entire mass and energy budget evolution driving this sudden condensation phenomenon. We find that the runaway cooling and rain formation also induces the appearance of dark post-flare loop systems, as observed in EUV channels. We confirm and augment earlier observational findings, suggesting that thermal conduction and radiative losses alternately dominate the cooling of PFLs.

A comparative analysis of the observational characteristics of fast radio bursts at the frequencies 111 and 1400 MHz is carried out. The distributions of radio bursts by the dispersion measure are constructed. At both frequencies, they are described by a lognormal distribution with the parameters $\mu =6.2$ $\sigma = 0.7$. The dependence $\tau_{sc}(DM)$ of the scattering value on the dispersion measure at 111 MHz and 1400 MHz is also constructed. This dependence is fundamentally different from the dependence for pulsars. A comparative analysis of the relationship between the scattering of pulses and the dispersion measure at 1400 MHz and 111 MHz showed that for both frequencies it has the form $\tau_{sc}(DM)\sim DM^k$, where $k = 0.49 \pm 0.18$ and $k = 0.43 \pm 0.15$ for the frequencies 111 and 1400 MHz, respectively. The obtained dependence is explained within the framework of the assumption of the extragalactic occurrence of fast radio bursts and an almost uniform distribution of matter in intergalactic space. From the dependence $\tau_{sc}(DM)$ a total estimate of the contribution to the matter of the halo of our and the host galaxy to $DM$ is obtained $DM_{halo} + \frac{DM_{host}}{1+z}\approx 60\;{\rm pc/cm}^3$. Based on the LogN - LogS dependence, the average spectral index of radio bursts is derived $\alpha = - 0.63 \pm 0.20$ provided that the statistical properties of these samples at 111 and 1400 MHz are the same.

Combining numerical simulations and analytical modeling, we investigate whether close binary systems form by the effect of magnetic braking. Using magnetohydrodynamics simulations, we calculate the cloud evolution with a sink, for which we do not resolve the binary system or binary orbital motion to realize long-term time integration. Then, we analytically estimate the binary separation using the accreted mass and angular momentum obtained from the simulation. In unmagnetized clouds, wide binary systems with separations of >100 au form, in which the binary separation continues to increase during the main accretion phase. In contrast, close binary systems with separations of <100 au can form in magnetized clouds. Since the efficiency of magnetic braking strongly depends on both the strength and configuration of the magnetic field, they also affect the formation conditions of a close binary. In addition, the protostellar outflow has a negative impact on close binary formation, especially when the rotation axis of the prestellar cloud is aligned with the global magnetic field. The outflow interrupts the accretion of gas with small angular momentum, which is expelled from the cloud, while gas with large angular momentum preferentially falls from the side of the outflow onto the binary system and widens the binary separation. This study shows that a cloud with a magnetic field that is not parallel to the rotation axis is a favorable environment for the formation of close binary systems.

The results of gamma-ray observations of the binary system HESS J0632+057 collected during 450 hours over 15 years, between 2004 and 2019, are presented. Data taken with the atmospheric Cherenkov telescopes H.E.S.S., MAGIC, and VERITAS at energies above 350 GeV were used together with observations at X-ray energies obtained with Swift-XRT, Chandra, XMM-Newton, NuSTAR, and Suzaku. Some of these observations were accompanied by measurements of the H{\alpha} emission line. A significant detection of the modulation of the VHE gamma-ray fluxes with a period of 316.7+-4.4 days is reported, consistent with the period of 317.3+-0.7 days obtained with a refined analysis of X-ray data. The analysis of data of four orbital cycles with dense observational coverage reveals short timescale variability, with flux-decay timescales of less than 20 days at very high energies. Flux variations observed over the time scale of several years indicate orbit-to-orbit variability. The analysis confirms the previously reported correlation of X-ray and gamma-ray emission from the system at very high significance, but can not find any correlation of optical H{\alpha} parameters with X-ray or gamma-ray energy fluxes in simultaneous observations. The key finding is that the emission of HESS J0632+057 in the X-ray and gamma-ray energy bands is highly variable on different time scales. The ratio of gamma-ray to X-ray flux shows the equality or even dominance of the gamma-ray energy range. This wealth of new data is interpreted taking into account the insufficient knowledge of the ephemeris of the system, and discussed in the context of results reported on other gamma-ray binary systems.

The surface-density profiles of dense filaments, in particular those traced by dust emission, appear to be well fit with Plummer profiles, i.e. Sigma(b)=Sigma_B+Sigma_O{1+[b/w_O]^2}^{[1-p]/2}. Here Sigma_B is the background surface-density; Sigma_B+Sigma_O is the surface-density on the filament spine; b is the impact parameter of the line-of-sight relative to the filament spine; w_O is the Plummer scale-length (which for fixed p is exactly proportional to the full-width at half-maximum, w_O=FWHM/2{2^{2/[p-1]}-1}^{1/2}); and p is the Plummer exponent (which reflects the slope of the surface-density profile away from the spine).} In order to improve signal-to-noise it is standard practice to average the observed surface-densities along a section of the filament, or even along its whole length, before fitting the profile. We show that, if filaments do indeed have intrinsic Plummer profiles with exponent p_INTRINSIC, but there is a range of w_O values along the length of the filament (and secondarily a range of Sigma_B values), the value of the Plummer exponent, p_FIT, estimated by fitting the averaged profile, may be significantly less than p_INTRINSIC. The decrease, Delta p = p_INTRINSIC - p_FIT, increases monotonically with increasing p_INTRINSIC; with increasing range of w_O values; and -- if, but only if, there is a finite range of w_O values -- with increasing range of Sigma_B values. For typical filament parameters the decrease is insignificant if p_INTRINSIC = 2 (0.05 <~ Delta p <~ 0.10), but for p_INTRINSIC =3 it is larger (0.18 <~ Delta p <~ 0.50), and for p_INTRINSIC =4 it is substantial (0.50 <~ Delta p <~ 1.15). On its own this effect is probably insufficient to support a value of p_INTRINSIC much greater than p_FIT ~ 2, but it could be important in combination with other effects.

Long-baseline interferometry of microlensing events can resolve the individual images of the source produced by the lens, which combined with the modeling of the microlensing light curve, lead to the exact lens mass and distance. Interferometric observations thus offer a unique opportunity to constrain the mass of exoplanets detected by microlensing, and to precisely measure the mass of distant isolated objects such as stars, brown dwarfs and stellar remnants like white dwarfs, neutron stars or stellar black holes. Having accurate models and reliable numerical methods is of particular importance, as the number of targets is expected to increase significantly in the near future. In this work, we discuss the different approaches to calculating the fringe complex visibility for the important case of a single lens. We propose a robust integration scheme to calculate the exact visibility, and introduce a novel approximation, that we call `thin-arcs approximation', which applies over a wide range of lens-source separations. We find that this approximation runs $\times6$ to $\times10$ times faster than the exact calculation, depending of the characteristics of the event and the required accuracy. This approximation provides accurate results for microlensing events of medium to higher magnification observed around the peak, i.e. a large fraction of potential observational targets.

Using the absolute astrometric positions and proper motions for common stars in the Hipparcos and Gaia catalogs separated by 24.75 years in the mean epoch, we compute mass ratios for long-period, resolved binary systems without any astrophysical assumptions or dependencies except the presence of inner binary subsystems that may perturb the observed mean proper motions. The mean epoch positions of binary companions from the Hipparcos Double and Multiple System Annex (DMSA) are used as the first epoch. The mean positions and proper motions of carefully cross-matched counterparts in Gaia EDR3 comprise the second epoch data. Selecting only results with sufficiently high signal-to-noise ratio and discarding numerous optical pairs, we construct a catalog of 248 binary systems, which is published online. Several cases with unusual properties or results are also discussed.

We investigate the constraining power of forthcoming Stage-IV weak lensing surveys (Euclid, LSST, and NGRST) for extensions to the LCDM model on small scales, via their impact on the cosmic shear power spectrum. We use high-resolution cosmological simulations to calculate how warm dark matter (WDM), self-interacting dark matter (SIDM) and a running of the spectral index affect the non-linear matter power spectrum, P(k), as a function of scale and redshift. We evaluate the cosmological constraining power using synthetic weak lensing observations derived from these power spectra and that take into account the anticipated source densities, shape noise and cosmic variance errors of upcoming surveys. We show that upcoming Stage-IV surveys will be able to place useful, independent constraints on both WDM models (ruling out models with a particle mass of < 0.5 keV) and SIDM models (ruling out models with a velocity-independent cross-section of > 10 cm^2 g^-1) through their effects on the small-scale cosmic shear power spectrum. Similarly, they will be able to strongly constrain cosmologies with a running spectral index. Finally, we explore the error associated with the cosmic shear cross-spectrum between tomographic bins, finding that it can be significantly affected by Poisson noise (the standard assumption is that the Poisson noise cancels between tomographic bins). We provide a new analytic form for the error on the cross-spectrum which accurately captures this effect.

BayesicFitting is a comprehensive, general-purpose toolbox for simple and standardized model fitting. Its fitting options range from simple least-squares methods, via maximum likelihood to fully Bayesian inference, working on a multitude of available models. BayesicFitting is open source and has been in development and use since the 1990s. It has been applied to a variety of science applications, chiefly in astronomy. BayesicFitting consists of a collection of PYTHON classes that can be combined to solve quite complicated inference problems. Amongst the classes are models, fitters, priors, error distributions, engines, samples, and of course NestedSampler, our general-purpose implementation of the nested sampling algorithm. Nested sampling is a novel way to perform Bayesian calculations. It can be applied to inference problems, that consist of a parameterized model to fit measured data to. NestedSampler calculates the Bayesian evidence as the numeric integral over the posterior probability of (hyper)parameters of the problem. The solution in terms of the parameters is obtained as a set of weighted samples drawn from the posterior. In this paper, we emphasize nested sampling and all classes that are directly connected to it. Additionally, we present the fitters, which fit the data by the least-squares method or the maximum likelihood method. They can also calculate the Bayesian evidence as a Gaussian approximation. We will discuss the architecture of the toolbox. Which classes are present, what is their function, how they are related and implementational details where it gets complicated.

The low frequency component of the Square Kilometre Array (SKA1-Low) will be an aperture phased array located at the Murchison Radio-astronomy Observatory (MRO) site in Western Australia. It will be composed of 512 stations, each of them consisting of 256 log-periodic dual polarized antennas, and will operate in the low frequency range (50 MHz - 350 MHz) of the SKA bandwidth. The Aperture Array Verification System 2 (AAVS2), operational since late 2019, is the last full-size engineering prototype station deployed at the MRO site before the start of the SKA1-Low construction phase. The aim of this paper is to characterize the station performance through commissioning observations at six different frequencies (55, 70, 110, 160, 230 and 320 MHz) collected during its first year of activities. We describe the calibration procedure, present the resulting all-sky images and their analysis, and discuss the station calibratability and system stability. Using the difference imaging method, we also derive estimates of the SKA1-Low sensitivity for the same frequencies, and compare them to those obtained through electromagnetic simulations across the entire telescope bandwidth, finding good agreement (within $\leq 13%$). Moreover, our estimates exceed the SKA1-Low requirements at all the considered frequencies, by up to a factor of $\sim$2.3. Our results are very promising and allow an initial validation of the AAVS2 prototype station performance, which is an important step towards the upcoming SKA-Low telescope construction and science.

In this paper, we investigate the early-time optical$-$near-infrared (NIR) spectral energy distributions (SEDs) of four Type Ibn supernovae (SNe). We find that the SEDs of SN~2010al, LSQ13ddu, and SN~2015G can be well explained by the single-component blackbody model, while the SEDs of OGLE-2012-SN-006 cannot. We invoke the double-component model assuming that the SEDs were produced by the SN photosphere and the heated dust to fit the optical$-$NIR SEDs of the four SNe Ibn, finding that the derived temperatures of the dust associated with OGLE-2012-SN-006 favor the scenario that the dust consists of the graphite grains, and the mass and temperature of dust are $\sim$$0.5-2.0\times10^{-3}~M_\odot$ and $\sim$ $1200-1300$ K, respectively. Moreover, our fits for SN~2010al, LSQ13ddu, and SN~2015G show that the upper limits of the masses of the dust associated with the three SNe Ibn are respectively $1.45\times 10^{-5}~M_\odot$, $5.9\times 10^{-7}~M_\odot$, and $2.4\times 10^{-7}~M_\odot$. A further analysis demonstrates that the inferred radius of the dust shell surrounding OGLE-2012-SN-006 is significantly larger than that of the SN ejecta at early epochs, indicating that the NIR excesses of the SEDs of OGLE-2012-SN-006 were produced by a preexisting dust shell. Our study for the early-time SEDs of four SNe Ibn, together with the previous studies and the fact that some SNe showed the evidence of dust formation at the late-time SEDs, indicates that at least $\sim$1/3 of SNe Ibn show evidence for dust formation.

We present pulse width measurements for a sample of radio pulsars observed with the MeerKAT telescope as part of the Thousand-Pulsar-Array (TPA) programme in the MeerTime project. For a centre frequency of 1284 MHz, we obtain 762 $W_{10}$ measurements across the total bandwidth of 775 MHz, where $W_{10}$ is the width at the 10% level of the pulse peak. We also measure about 400 $W_{10}$ values in each of the four or eight frequency sub-bands. Assuming, the width is a function of the rotation period P, this relationship can be described with a power law with power law index $\mu=-0.29\pm 0.03$. However, using orthogonal distance regression, we determine a steeper power law with $\mu=-0.63\pm 0.06$. A density plot of the period-width data reveals such a fit to align well with the contours of highest density. Building on a previous population synthesis model, we obtain population-based estimates of the obliquity of the magnetic axis with respect to the rotation axis for our pulsars. Investigating the width changes over frequency, we unambiguously identify a group of pulsars that have width broadening at higher frequencies. The measured width changes show a monotonic behaviour with frequency for the whole TPA pulsar population, whether the pulses are becoming narrower or broader with increasing frequency. We exclude a sensitivity bias, scattering and noticeable differences in the pulse component numbers as explanations for these width changes, and attempt an explanation using a qualitative model of five contributing Gaussian pulse components with flux density spectra that depend on their rotational phase.

As modern-day precision cosmology aims for statistical uncertainties of the percent level or lower, it becomes increasingly important to reconsider estimator assumptions at each step of the process, and their consequences on the statistical variability of the scientific results. We compare $L^1$ regression statistics to the weighted mean, the canonical $L^2$ method based on Gaussian assumptions, for inference of the weak gravitational shear signal from a catalog of background ellipticity measurements around a sample of clusters, in many recent analyses a standard step in the process. We use the shape measurements of background sources around 6925 AMICO clusters detected in the KiDS 3rd data release. We investigate the robustness of our results and the dependence of uncertainties on the signal-to-noise ratios of the background source detections. Using a halo model approach, we derive lensing masses from the estimated excess surface density profiles. The highly significant shear signal allows us to study the scaling relation between the $r$-band cluster luminosity $L_{200}$, and the derived lensing mass $M_{200}$. We show the results of the scaling relations derived in 13 bins in $L_{200}$, with a tightly constrained power law slope of $\sim 1.24\pm 0.08$. We observe a small, but significant relative bias of a few percent in the recovered excess surface density profiles between the two regression methods, which translates to a $1\sigma$ difference in $M_{200}$. The efficiency of $L^1$ is at least that of the weighted mean, relatively increasing with higher signal-to-noise shape measurements. Our results indicate the relevance of optimizing the estimator for infering the gravitational shear from a distribution of background ellipticities. The interpretation of measured relative biases can be gauged by deeper observations, while increased computation times remain feasible.

The 1-D evolution equations for warped discs come in two flavors: For very viscous discs the internal torque vector G is uniquely determined by the local conditions in the disc, and warps tend to damp out rapidly if they are not continuously driven. For very inviscid discs, on the other hand, G becomes a dynamic quantity, and a warp will propagate through the disc as a wave. The equations governing both regimes are usually treated separately. A unified set of equations was postulated recently by Martin et al. (2019), but not yet derived from the underlying physics. The standard method for deriving these equations is based on a perturbation series expansion, which is a powerful, but somewhat abstract technique. A more straightforward method is to employ the warped shearing box framework of Ogilvie and Latter (2013), which so far has not yet been used to derive the equations for the wavelike regime. The goal of this paper is to analyze the warped disc equations in both regimes using the warped shearing box framework, to derive a unified set of equations, valid for small warps, and to discuss how our results can be interpreted in terms of the affine tilted-slab approach of Ogilvie (2018).

Ultra-short period (USP) planets are exoplanets which have orbital periods of less than one day and are unique because they orbit inside the nominal magnetic truncation gap of their host stars. In some cases, USP planets have also been observed to exhibit unique dynamical parameters such as significant misalignments in inclination angle with respect to nearby planets. In this paper, we explore how the geometry of a multi-planet system hosting a USP planet can be expected to evolve as a star ages. In particular, we explore the relationship between the mutual inclination of the USP planet and the quadrupole moment ($J_2$) of the host star. We use secular perturbation theory to predict the past evolution of the example TOI-125 system, and then confirm the validity of our results using long-term N-body simulations. Through investigating how the misalignment between the candidate USP planet and the three other short-period planets in the TOI-125 system arose, we intend to derive a better understanding of the population of systems with misaligned USP planets and how their observed parameters can be explained in the context of their dynamical histories.

Satisfactory description of gravitational and gravity potentials is needed for a proper modelling of a wide spectrum of physical problems on various size scales, ranging from atmosphere dynamics up to the movements of stars in a galaxy. In certain cases, Similar Oblate Spheroidal (SOS) coordinate system can be of advantage for such modelling tasks, mainly inside or in the vicinity of oblate spheroidal objects (planets, stars, galaxies). Although the solution of the relevant expressions for the SOS system cannot be written in a closed form, it can be derived as analytical expressions -- convergent infinite power series. Explicit analytical expressions for the Cartesian coordinates in terms of the curvilinear Similar Oblate Spheroidal coordinates are derived in the form of infinite power series with generalized binomial coefficients. The corresponding partial derivatives are found in a suitable form, further enabling derivation of the metric scale factors necessary for differential operations. The terms containing derivatives of the metric scale factors in the velocity advection term of the momentum equation in SOS coordinate system are expressed. The Jacobian determinant is derived as well. The presented analytical solution of SOS coordinate system solution is a tool applicable for a broad variety of objects exhibiting density, gravity or gravitation levels resembling similar oblate spheroids. Such objects range from the bodies with small oblateness (the Earth itself on the first place), through elliptical galaxies up to significantly flattened objects like disk galaxies.

The Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) will soon carry out an unprecedented wide, fast, and deep survey of the sky in multiple optical bands. The data from LSST will open up a new discovery space in astronomy and cosmology, simultaneously providing clues toward addressing burning issues of the day, such as the origin of dark energy and and the nature of dark matter, while at the same time yielding data that will, in turn, pose fresh new questions. To prepare for the imminent arrival of this remarkable data set, it is crucial that the associated scientific communities be able to develop the software needed to analyze it. Computational power now available allows us to generate synthetic data sets that can be used as a realistic training ground for such an effort. This effort raises its own challenges -- the need to generate very large simulations of the night sky, scaling up simulation campaigns to large numbers of compute nodes across multiple computing centers with different architectures, and optimizing the complex workload around memory requirements and widely varying wall clock times. We describe here a large-scale workflow that melds together Python code to steer the workflow, Parsl to manage the large-scale distributed execution of workflow components, and containers to carry out the image simulation campaign across multiple sites. Taking advantage of these tools, we developed an extreme-scale computational framework and used it to simulate five years of observations for 300 square degrees of sky area. We describe our experiences and lessons learned in developing this workflow capability, and highlight how the scalability and portability of our approach enabled us to efficiently execute it on up to 4000 compute nodes on two supercomputers.

We present a study of the relationships and environmental dependencies between stellar mass, star formation rate, and gas metallicity for more than 700 galaxies in groups up to redshift 0.35 from the Galaxy And Mass Assembly (GAMA) survey. To identify the main drivers, our sample was analyzed as a function of group-centric distance, projected galaxy number density, and stellar mass. By using control samples of more than 16000 star-forming field galaxies and volume limited samples, we find that the highest enhancement in SFR (0.3 dex) occurs in galaxies with the lowest local density. In contrast to previous work, our data show small enhancements of $\sim$0.1 dex in SFR for galaxies at the highest local densities or group-centric distances. Our data indicates quenching in SFR only for massive galaxies, suggesting that stellar mass might be the main driver of quenching processes for star forming galaxies. We can discard a morphological driven quenching, since the S\'ersic index distribution for group and control galaxies are similar. The gas metallicity does not vary drastically. It increases $\sim$0.08 dex for galaxies at the highest local densities, and decreases for galaxies at the highest group-centric distances, in agreement with previous work. Altogether, the local density, rather than group-centric distance, shows the stronger impact in enhancing both, the SFR and gas metallicity. We applied the same methodology to galaxies from the IllustrisTNG simulations, and although we were able to reproduce the general observational trends, the differences between group and control samples only partially agree with the observations

Axions and other light particles appear ubiquitously in physics beyond the Standard Model, with a variety of possible couplings to ordinary matter. Cosmology offers a unique probe of these particles as they can thermalize in the hot environment of the early universe for any such coupling. For sub-MeV particles, their entropy must leave a measurable cosmological signal, usually via the effective number of relativistic particles, $N_\mathrm{eff}$. In this paper, we will revisit the cosmological constraints on the couplings of axions and other pseudo-Nambu-Goldstone bosons to Standard Model fermions from thermalization below the electroweak scale, where these couplings are marginal and give contributions to the radiation density of $\Delta N_\mathrm{eff} > 0.027$. We update the calculation of the production rates to eliminate unnecessary approximations and find that the cosmological bounds on these interactions are complementary to astrophysical constraints, e.g. from supernova SN 1987A. We additionally provide quantitative explanations for these bounds and their relationship.

We study the gravitational-wave signal stemming from strongly coupled models featuring both, dark chiral and confinement phase transitions. We therefore identify strongly coupled theories that can feature a first-order phase transition. Employing the Polyakov-Nambu-Jona-Lasinio model, we focus our attention on SU(3) Yang-Mills theories featuring fermions in fundamental, adjoint, and two-index symmetric representations. We discover that for the gravitational-wave signals analysis, there are significant differences between the various representations. Interestingly we also observe that the two-index symmetric representation leads to the strongest first-order phase transition and therefore to a higher chance of being detected by the Big Bang Observer experiment. Our study of the confinement and chiral phase transitions is further applicable to extensions of the Standard Model featuring composite dynamics.

In Chain Inflation the universe tunnels along a series of false vacua of ever-decreasing energy. The main goal of this paper is to embed Chain Inflation in high energy fundamental physics. We begin by illustrating a simple effective formalism for calculating Cosmic Microwave Background (CMB) observables in Chain Inflation. Density perturbations seeding the anisotropies emerge from the probabilistic nature of tunneling (rather than from quantum fluctuations of the inflation). To obtain the correct normalization of the scalar power spectrum and the scalar spectral index, we find an upper limit on the scale of inflation at horizon crossing of CMB scales, $V_*^{1/4}< 10^{12}$~GeV. We then provide an explicit realization of chain inflation, in which the inflaton is identified with an axion in supergravity. The axion enjoys a perturbative shift symmetry which is broken to a discrete remnant by instantons. The model, which we dub `natural chain inflation' satisfies all cosmological constraints and can be embedded into a standard $\Lambda$CDM cosmology. Our work provides a major step towards the ultraviolet completion of chain inflation in string theory.

We study a dark $SU(2)_D$ gauge extension of the standard model (SM) with the possibility of a strong first order phase transition (FOPT) taking place below the electroweak scale in the light of NANOGrav 12.5 yr data. As pointed out recently by the NANOGrav collaboration, gravitational waves (GW) from such a FOPT with appropriate strength and nucleation temperature can explain their 12.5 yr data. We impose a classical conformal invariance on the scalar potential of $SU(2)_D$ sector involving only a complex scalar doublet with negligible couplings with the SM Higgs. While a FOPT at sub-GeV temperatures can give rise to stochastic GW around nano-Hz frequencies being in agreement with NANOGrav findings, the $SU(2)_D$ vector bosons which acquire masses as a result of the FOPT in dark sector, can also serve as dark matter (DM) in the universe. The relic abundance of such vector DM can be generated in a non-thermal manner from the SM bath via scalar portal mixing. We also discuss future sensitivity of gravitational wave experiments to the model parameter space.

Universal relations (i.e., insensitive to the equation of state) between macroscopic properties of neutron stars (NSs) have proven useful for a variety of applications -- from providing a direct means to extract observables from data to breaking degeneracies that hinder tests of general relativity. Similarly, equation-of-state insensitive relations directly connecting macroscopic and microscopic properties of NSs can potentially provide a clean window into the behavior of nuclear matter. In this Letter, we uncover a tight correlation between certain macroscopic properties of NSs -- its compactness, moment of inertia or tidal deformability -- and the ratio of central pressure to central energy density, which can be interpreted as a mean notion of the stiffness of nuclear matter inside a NS. We describe interesting properties of this stiffness measure, quantify its universality, and explore its consequences in the face of recent NS observations.

We develop a general formalism for treating radiative degrees of freedom near $\mathscr{I}^{+}$ in theories with an arbitrary Ricci-flat internal space. These radiative modes are encoded in a generalized news tensor which decomposes into gravitational, electromagnetic, and scalar components. We find a preferred gauge which simplifies the asymptotic analysis of the full nonlinear Einstein equations and makes the asymptotic symmetry group transparent. This asymptotic symmetry group extends the BMS group to include angle-dependent isometries of the internal space. We apply this formalism to study memory effects, which are expected to be observed in future experiments, that arise from bursts of higher-dimensional gravitational radiation. We outline how measurements made by gravitational wave observatories might probe properties of the compact extra dimensions.

We present a discontinuous Galerkin-finite-difference hybrid scheme that allows high-order shock capturing with the discontinuous Galerkin method for general relativistic magnetohydrodynamics. The hybrid method is conceptually quite simple. An unlimited discontinuous Galerkin candidate solution is computed for the next time step. If the candidate solution is inadmissible, the time step is retaken using robust finite-difference methods. Because of its a posteriori nature, the hybrid scheme inherits the best properties of both methods. It is high-order with exponential convergence in smooth regions, while robustly handling discontinuities. We give a detailed description of how we transfer the solution between the discontinuous Galerkin and finite-difference solvers, and the troubled-cell indicators necessary to robustly handle slow-moving discontinuities and simulate magnetized neutron stars. We demonstrate the efficacy of the proposed method using a suite of standard and very challenging 1d, 2d, and 3d relativistic magnetohydrodynamics test problems. The hybrid scheme is designed from the ground up to efficiently simulate astrophysical problems such as the inspiral, coalescence, and merger of two neutron stars.

Dark matter (DM) particles captured by the Sun can produce high energy electrons outside the Sun through annihilating into meta-stable mediators. The corresponding cosmic-ray electron signals observed by the space-based experiments will be time dependent due to the orbital motion of the space-based detectors. The shape of this time dependence is predictable given the orbital information of the detectors. Since the high-energy CR electron (with energy E>100 GeV) fluxes are expected to be constant in time, non-observation of such time variation can be used to place upper limits on the DM annihilation cross section. We analyze the time dependence of dark matter cosmic-ray signals in three space-based experiments: AMS-02, DAMPE and CALET. Under the assumption that no time dependent signal is observed, we derive the 95% C.L. exclusion limits on the signal strength from the current data. We map our limits onto the parameter space of the dark photon model and find that the constraints are comparable with that derived from the supernova SN1987A.

We perform a general analysis of the cosmological viability of Geometric Inflation. We show that the evolution of the universe, from inflation to the present day, can be seen from the addition of an infinite tower of curvature invariants into the Hilbert-Einstein action. The main epochs of the Universe can be reproduced: Inflation, Big Bang Nucleosynthesis, and Late-time acceleration driven by the cosmological constant. The slow-roll condition is a robust prediction of the theory. Inflation possesses a graceful exit with enough number of $e-$folds between the limit imposed by the Planck density and the exit of the exponential expansion to solve the horizon problem and the absence of topological defects. We also provide some scenarios where the energy scale of the theory can be calibrated.

It has been known that a non-perfect fluid that accounts for dissipative viscous effects can evade a highly anisotropic chaotic mixmaster approach to a singularity. Viscosity is often simply parameterised in this context, so it remains unclear whether isotropisation can really occur in physically motivated contexts. We present a few examples of microphysical manifestations of viscosity in fluids that interact either gravitationally or, for a scalar field for instance, through a self-coupling term in the potential. In each case, we derive the viscosity coefficient and comment on the applicability of the approximations involved when dealing with dissipative non-perfect fluids. Upon embedding the fluids in a cosmological context, we then show the extent to which these models allow for isotropisation of the universe in the approach to a singularity. We first do this in the context of expansion anisotropy only, i.e., in the case of a Bianchi type-I universe. We then include anisotropic 3-curvature modelled by the Bianchi type-IX metric. It is found that a self-interacting scalar field at finite temperature allows for efficient isotropisation, whether in a Bianchi type-I or type-IX spacetime, although the model is not tractable all the way to a singularity. Mixmaster chaotic behaviour, which is well known to arise in anisotropic models including anisotropic 3-curvature, is found to be suppressed in the latter case as well. We find that the only model permitting an isotropic singularity is that of a dense gas of black holes.

On the 17th of August 2017, the thriving discovery of gravitational waves event GW170817 and its optical counterpart GRB170817A, owing to coalescing of two neutron starts divulged a very small amount of difference around ${\cal O}(10^{-16})$ between traveling speed of light and the velocity of gravitational waves ($C_{T}$). This small deviation can be used as a strong constraint on modified gravity models. We concentrate on the Higher-Order expansion of Mimetic Gravity (HOMimG) model, which exposes in the low-energy limit of projectable Ho\v{r}ava-Lifshitz gravity. Thence, we specify the parametric space of three parameters of our model ($a$, $b$, and $c$ ) utilizing the observational constraint from GW170817-GRB170817A on $C_{T}$, beside two theoretical constraints on $C_{T}^{2}$ and $C_{s}^{2}$ due to assurance of the stability of the model and subluminal promulgation of the scalar and tensor perturbations. Thereafter, we increase the accuracy of the parametric space with the aid of imposing further limitation of $\gamma$ parameter (related to the latest observational value for the age of the Universe) in addition to the previous ones. In pursuance of determining the other parameter of the model ($\lambda$), the potential of the model is specified, and another observational bound related to the Equation of State (EoS) parameter of dark energy is taken into account in addition to the erstwhile ones. In consequence, we attain a viable HOMimG model with four parameter ($a$, $b$, $c$, and $\lambda$) confined to numbers of observational and theoretical constraints. At the end, regarding the concluded numerical ranges for the model parameters, and cogitating two different potential (quadratic and quartic potentials) to specify $\lambda$ parameter, we illustrate that the values of the model parameters are independent of the form of potential.

A compressible liquid-drop approach adjusted to uniform matter many-body calculations based on chiral EFT interactions and to the experimental nuclear masses is used to investigate the neutron star crust properties. Eight chiral EFT hamiltonians and a representative phenomenological force (SLy4) are confronted. We show that some properties of the crust, e.g. clusters mass, charge, and asymmetry, are mostly determined by symmetric matter properties close to saturation density and are therefore mainly constrained by experimental nuclear masses, while other properties, e.g., energy per particle, pressure, sound speed, are mostly influenced by low-density predictions in neutron matter, where chiral EFT and phenomenological forces substantially differ.

We try to constraints some of the nuclear matter parameters such as symmetry energy ($J$) and its slope ($L$) from the recent inferred data of the PREX-2. Other nuclear matter parameters are adopted from {\bf [Phys. Rev. C 85 035201 (2012), Phys. Rev. C 90 055203 (2014)]} papers and the linear correlation among them are checked by using the Pearson's formula. We find the correlation between $J-L$, $K_\tau-J$ and $K_\tau-L$ with coefficients 0.85, 0.81 and 0.76 respectively. The neutron star properties such as mass and radius are calculated with 50 unified equation of states. The results are consistent with recently observed pulsars and NICER data except few exceptions. From the radii constraints, we find that the new NICER data allows a narrow radius range contrary to a large range of PREX-2 and the old NICER data leaving us an inconclusive determination of the neutron star radius.

Electron number densities in stars and the Earth are inhomogeneous because of atomic electrons. The large inhomogeneities on atomic-scale tend to form at tops of respective layers of stars, and 1s electrons of O locally produce weak potentials higher than that of the high MSW resonance. Then, supernova neutrinos experience vast numbers of non-adiabatic transitions. This inhomogeneous electron potential generates finite amplitudes of all three propagation eigenstates, and wave packets effectively separate. Then, spectral differences between three flavors significantly diminishes after propagation.

Discontinuous Galerkin methods are popular because they can achieve high order where the solution is smooth, because they can capture shocks while needing only nearest-neighbor communication, and because they are relatively easy to formulate on complex meshes. We perform a detailed comparison of various limiting strategies presented in the literature applied to the equations of general relativistic magnetohydrodynamics. We compare the standard minmod/$\Lambda\Pi^N$ limiter, the hierarchical limiter of Krivodonova, the simple WENO limiter, the HWENO limiter, and a discontinous Galerkin-finite-difference hybrid method. The ultimate goal is to understand what limiting strategies are able to robustly simulate magnetized TOV stars without any fine-tuning of parameters. Among the limiters explored here, the only limiting strategy we can endorse is a discontinous Galerkin-finite-difference hybrid method.

The cosmological term, $\Lambda$, in Einstein's equations is an essential ingredient of the `concordance' $\Lambda$CDM model of cosmology. In this mini-review presentation, we assess the possibility that $\Lambda$ can be a dynamical quantity, more specifically a `running quantity' in quantum field theory in curved spacetime. A great deal of phenomenological works have shown in the last few years that this option (sometimes accompanied with a running gravitational coupling) may cure some of the tensions afflicting the $\Lambda$CDM. The `running vacuum models' (RVM's) are characterized by the vacuum energy density, $\rho_{\rm vac}$, being a series of (even) powers of the Hubble rate and its time derivatives. Here we describe the technical quantum field theoretical origin of the RVM structure in FLRW spacetime, which goes well-beyond the original semi-qualitative renormalization group arguments. In particular, we compute the renormalized energy-momentum tensor using the adiabatic regularization procedure and show that it leads to the RVM form. In other words, we find that the renormalized vacuum energy density, $\rho_{vac}(H)$ evolves as a (constant) additive term plus a leading dynamical components ${\cal O}(H^2)$. There are also ${\cal O}(H^4)$ contributions, which can be relevant for the early universe. Remarkably enough, the renormalized $\rho_{\rm vac}(H)$ does not exhibit dangerous terms proportional to the quartic power of the masses ($\sim m^4$) of the fields. It is well-known that these terms have been the main source of trouble since they are responsible for the extreme fine tuning and ultimately for the cosmological constant problem. In this context, the current $\rho_{vac}(H)$ is dominated by a constant term, as it should be, but it acquires a mild dynamical component $\sim \nu H^2$ ($|\nu|\ll1$) which makes the RVM to mimic quintessence.