Papers for Friday, Apr 02 2021

Ultra-compact dwarf galaxies (UCDs) were serendipitously discovered by spectroscopic surveys in the Fornax cluster twenty years ago. Nowadays, it is commonly accepted that many bright UCDs are the nuclei of galaxies that have been stripped. However, this conclusion might be driven by biased samples of UCDs in high-density environments, on which most searches are based. With the deep optical images of the Fornax Deep Survey, combined with public near-infrared data, we revisit the UCD population of the Fornax cluster and search for UCD candidates, for the first time, systematically out to the virial radius of the galaxy cluster. Our search is complete down to magnitude m$_g$ = 21 mag or M$_g$ $\sim$ -10.5 mag at the distance of the Fornax cluster. The UCD candidates are identified and separated from foreground stars and background galaxies by their optical and near-infrared colours. This primarily utilizes the $u-i$/$i-Ks$ diagram and a machine learning technique is employed to incorporate other colour combinations to reduce the number of contaminants. The newly identified candidates (44) in addition to the spectroscopically confirmed UCDs (61), increases the number of known Fornax UCD considerably (105). Almost all of the new UCD candidates are located outside the Fornax cluster core (360 kpc), where all of the known UCDs were found. The distribution of UCDs within the Fornax cluster shows that a population of UCDs may form in low-density environments. This most likely challenges the current models of UCD formation.

The dark matter content of the Fornax dwarf spheroidal galaxy inferred from its kinematics is substantially lower than expected from LCDM cosmological simulations. We use N-body simulations to examine whether this may be the result of Galactic tides. We find that, despite improved proper motions from the Gaia mission, the pericentric distance of Fornax remains poorly constrained, mainly because its largest velocity component is roughly anti-parallel to the solar motion. Translating Fornax's proper motion into a Galactocentric velocity is thus sensitively dependent on Fornax's assumed distance: the observed distance uncertainty, $\pm 8\%$, implies pericentric distances that vary between $r_{\rm peri}\sim 50$ and $r_{\rm peri}\sim 150$ kpc. Our simulations show that for $r_{\rm peri}$ in the lower range of that estimate, an LCDM subhalo with maximum circular velocity $V_{\rm max}=40$ km s$^{-1}$ (or virial mass $M_{200}\approx 10^{10} M_\odot$, as expected from LCDM) would be tidally stripped to $V_{\rm max} \sim 23$ km s$^{-1}$ over $10$ Gyr. This would reduce the dark mass within the Fornax stellar half-mass radius to about half its initial value, bringing it into agreement with observations. Tidal stripping affects mainly Fornax's dark matter halo; its stellar component is affected little, losing less than $5\%$ of its initial mass in the process. We also explore the effect of Galactic tides on the dynamical friction decay times of Fornax's population of globular clusters (GC) and find little evidence for substantial changes, compared with models run in isolation. A population of GCs with initial orbital radii between $1$ and $2$ kpc is consistent with the present-day spatial distribution of Fornax GCs, despite assuming a cuspy halo. Both the dark matter content and the GC population of Fornax seem consistent with current models of dwarf galaxy formation in LCDM.

Using data from the HAWC gamma-ray Telescope, we have studied a sample of 37 millisecond pulsars (MSPs), selected for their spindown power and proximity. From among these MSP, we have identified four which favor the presence of very high-energy gamma-ray emission at a level of $(2\Delta \ln \mathcal{L})^{1/2} \ge 2.5$. Adopting a correlation between the spindown power and gamma-ray luminosity of each pulsar, we performed a stacked likelihood analysis of these 37 MSPs, finding that the data supports the conclusion that these sources emit very high-energy gamma-rays at a level of $(2\Delta \ln \mathcal{L})^{1/2} = 4.24$. Among sets of randomly selected sky locations within HAWC's field-of-view, less than 1\% of such realizations yielded such high statistical significance. Our analysis suggests that MSPs produce very high-energy gamma-ray emission with a similar efficiency to that observed from the Geminga TeV-halo, $\eta_{\rm MSP} = (0.39-1.08) \times \eta_{\rm Geminga}$. This conclusion poses a significant challenge for pulsar interpretations of the Galactic Center gamma-ray excess, as it suggests that any population of MSPs potentially capable of producing the GeV excess would also produce TeV-scale emission in excess of that observed by HESS from this region. Future observations by CTA will be able to substantially clarify this situation.

We report photometric estimates of effective temperature, $T_{\rm eff}$, metallicity, [Fe/H], carbonicity, [C/Fe], and absolute carbon abundances, $A{\rm (C)}$, for over 700,000 stars from the Southern Photometric Local Universe Survey (S-PLUS) Data Release 2, covering a substantial fraction of the equatorial Sloan Digital Sky Survey Stripe 82. We present an analysis for two stellar populations: 1) halo main-sequence turnoff stars and 2) K-dwarf stars of mass $0.58 < M/M_{\odot} <0.75$ in the Solar Neighborhood. Application of the Stellar Photometric Index Network Explorer (SPHINX) to the mixed-bandwidth (narrow- plus wide-band) filter photometry from S-PLUS produces robust estimates of the metallicities and carbon abundances in stellar atmospheres over a wide range of temperature, $4250 < T_{\rm eff} \textrm{(K)} < 7000$. The use of multiple narrow-band S-PLUS filters enables SPHINX to achieve substantially lower levels of "catastrophic failures" (large offsets in metallicity estimates relative to spectroscopic determinations) than previous efforts using a single metallicity-sensitive narrow-band filter. We constrain the exponential slope of the Milky Way's K-dwarf halo metallicity distribution function (MDF), $\lambda_{10, \textrm{[Fe/H]}} = 0.85 \pm 0.21$, over the metallicity range $-2.5 < \textrm{[Fe/H]} < -1.0$; the MDF of our local-volume K-dwarf sample is well-represented by a gamma distribution with parameters $\alpha=2.8$ and $\beta=4.2$. S-PLUS photometry obtains absolute carbon abundances with a precision of $\sim 0.35$dex for stars with $T_{\rm eff} < 6500$K. We identify 364 candidate carbon-enhanced metal-poor stars, obtain assignments of these stars into the Yoon-Beers morphological groups in the $A$(C)-[Fe/H] space, and derive the CEMP frequencies.

Broadband tomography statistically extracts the redshift distribution of frequency dependent emission from the cross correlation of intensity maps with a reference catalog of galaxy tracers. We make forecasts for the performance of future all-sky UV experiments doing broadband tomography. We consider the Cosmological Advanced Survey Telescope for Optical-UV Research (CASTOR) and the Spectro-Photometer for the History of the Universe, Epoch of Reionization, and Ices Explorer (SPHEREx). We show that under reasonable error models, CASTOR measures the UV background SED 2-10 times better than existing data. It will also expand the applicable redshift range from the current z < 1 to z~0-3 with CASTOR and z=5-9 with SPHEREx. We show that CASTOR can provide competitive constraints on the EBL monopole to those available from galaxy number counts and direct measurement techniques. At high redshift especially, these results will help understand galaxy formation and reionization. Our modelling code and chains are publicly available.

The morphological classification of galaxies is a relevant probe for galaxy evolution and unveils its connection with cosmological structure formation. To this scope, it is fundamental to recover galaxy morphologies over large areas of the sky. In this paper, we present a morphological catalogue for galaxies in the Stripe-82 area, observed with S-PLUS, till a magnitude limit of $r\le17$, using the state-of-the-art of Convolutional Neural Networks (CNNs) for computer vision. This analysis will then be extended to the whole S-PLUS survey data, covering $\simeq 9300$ $deg^{2}$ of the celestial sphere in twelve optical bands. We find that the network's performance increases with 5 broad bands and additional 3 narrow bands compared to our baseline with 3 bands. However, it does lose performance when using the full $12$ band image information. Nevertheless, the best result is achieved with 3 bands, when using pre-trained network weights in an ImageNet dataset. These results enhance the importance of previous knowledge in the neural network weights based on training in non related extensive datasets. Thus, we release a model pre-trained in several bands that could be adapted to other surveys. We develop a catalogue of 3274 galaxies in Stripe-82 that are not presented in Galaxy Zoo 1 (GZ1). We also add classification to 4686 galaxies considered ambiguous in GZ1 dataset. Finally, we present a prospect of a novel way to take advantage of $12$ bands information for morphological classification using multiband morphometric features. The morphological catalogues are publicly available.

The Southern Photometric Local Universe Survey (S-PLUS) is an ongoing survey of $\sim$9300 deg$^2$ in the southern sky in a 12-band photometric system. This paper presents the second data release (DR2) of S-PLUS, consisting of 514 tiles covering an area of 950 deg$^2$. The data has been fully calibrated using a new photometric calibration technique suitable for the new generation of wide-field multi-filter surveys. This technique consists of a $\chi^2$ minimisation to fit synthetic stellar templates to already calibrated data from other surveys, eliminating the need for standard stars and reducing the survey duration by $\sim$15\%. We compare the template-predicted and S-PLUS instrumental magnitudes to derive the photometric zero-points (ZPs). We show that these ZPs can be further refined by fitting the stellar templates to the 12 S-PLUS magnitudes, which better constrain the models by adding the narrow-band information. We use the STRIPE82 region to estimate ZP errors, which are $\lesssim10$ mmags for filters J0410, J0430, $g$, J0515, $r$, J0660, $i$, J0861 and $z$; $\lesssim 15$ mmags for filter J0378; and $\lesssim 25$ mmags for filters $u$ and J0395. We describe the complete data flow of the S-PLUS/DR2 from observations to the final catalogues and present a brief characterisation of the data. We show that, for a minimum signal-to-noise threshold of 3, the photometric depths of the DR2 range from 19.9 mag to 21.3 mag (measured in Petrosian apertures), depending on the filter. The S-PLUS DR2 can be accessed from the website: https://splus.cloud}{https://splus.cloud.

In this work we establish and test methods for implementing dynamical friction for massive black hole pairs that form in large volume cosmological hydrodynamical simulations which include galaxy formation and black hole growth. We verify our models and parameters both for individual black hole dynamics and for the black hole population in cosmological volumes. Using our model of dynamical friction (DF) from collisionless particles, black holes can effectively sink close to the galaxy center, provided that the black hole's dynamical mass is at least twice that of the lowest mass resolution particles in the simulation. Gas drag also plays a role in assisting the black holes' orbital decay, but it is typically less effective than that from collisionless particles, especially after the first billion years of the black hole's evolution. DF from gas becomes less than $1\%$ of DF from collisionless particles for BH masses $> 10^{7}$ M$_{\odot}$. Using our best DF model, we calculate the merger rate down to $z=1.1$ using an $L_{\rm box}=35$ Mpc$/h$ simulation box. We predict $\sim 2$ mergers per year for $z>1.1$ peaking at $z\sim 2$. These merger rates are within the range obtained in previous work using similar-resolution hydro-dynamical simulations. We show that the rate is enhanced by factor of $\sim 2$ when DF is taken into account in the simulations compared to the no-DF run. This is due to $>40\%$ more black holes reaching the center of their host halo when DF is added.

NASA's Great Observatories have opened up the electromagnetic spectrum from space, providing sustained access to wavelengths not accessible from the ground. Together, Hubble, Compton, Chandra, and Spitzer have provided the scientific community with an agile and powerful suite of telescopes with which to attack broad scientific questions, and react to a rapidly changing scientific landscape. As the existing Great Observatories age, or are decommissioned, community access to these wavelengths will diminish, with an accompanying loss of scientific capability. This report, commissioned by the NASA Cosmic Origins, Physics of the Cosmos and Exoplanet Exploration Program Analysis Groups (PAGs), analyzes the importance of multi-wavelength observations from space during the epoch of the Great Observatories, providing examples that span a broad range of astrophysical investigations.

The cold classical Kuiper Belt Objects (KBOs) possess a high, $\gtrsim30\%$ binary fraction. Widely separated and dynamically fragile, these binary systems have been useful in tracing the origins of KBOs. A new class of binaries was recently identified by their colours. The so-called blue binaries are unanimously members of the less red compositional class, and exhibit a 100% binary fraction. They appear to be push-out survivors, emplaced in the classical region during Neptune's phases of outward migration. The presence of these binary systems implies that the majority of objects that formed near the cold classical region formed as binaries. Here we present new optical colour measurements of cold classical KBOs from the Colours of the Outer Solar System Origins Survey, including colours of a blue binary discovered by the Solar System Origins Legacy Survey -- 2015 RJ277. The increased size of the colours sample has resulted in order-of-magnitude decrease in the probability that the binaries and singles sample share the same colour distribution. From the Anderson-Darling statistic, this probability is only a 0.3%, while it is only 0.002% when utilizing the difference of means statistic. We find a hint that the blue binaries have inflated free inclinations compared to their red counterparts, consistent with the push-out origin for these bodies.

In this work we propose a novel campaign for constraining relativistically compact MACHO dark matter, such as primordial black holes (PBHs), using the moon as a detector. PBHs of about $10^{19} \textrm{ g}$ to $10^{22} \textrm{ g}$ may be sufficiently abundant to have collided with the moon in the history of the solar system. We show that the crater profiles of a PBH collision differ from traditional impactors and may be detectable in high resolution lunar surface scans now available. Any candidates may serve as sites for in situ measurements to identify high pressure phases of matter which may have formed near the PBH during the encounter. While we primarily consider PBH dark matter, the discussion generalises to the entire family of MACHO candidates with relativistic compactness. Moreover, we focus on the Moon since it has been studied well, but the same principles can be applied to other rocky bodies in our solar system without an atmosphere.

X-ray observations provide a potentially powerful tool to study starburst feedback. The analysis and interpretation of such observations remain challenging, however, due to various complications, including the non-isothermality of the diffuse hot plasma and the inhomogeneity of the foreground absorption. We here illustrate such complications and a way to mitigate their effects by presenting an X-ray spectroscopy of the 30 Doradus nebula in the Large Magellanic Clouds, based on a 100 ks Suzaku observation. We measure the thermal and chemical properties of the hot plasma and quantitatively confront them with the feedback expected from embedded massive stars. We find that our spatially resolved measurements can be well reproduced by a global modeling of the nebula with a log-normal temperature distribution of the plasma emission measure and a log-normal foreground absorption distribution. The metal abundances and total mass of the plasma are consistent with the chemically enriched mass ejection expected from the central OB association and a ~55% mass-loading from the ambient medium. The total thermal energy of the plasma is smaller than what is expected from a simple superbubble model, demonstrating that important channels of energy loss are not accounted for. Our analysis indeed shows tentative evidence for a diffuse non-thermal X-ray component, indicating that cosmic-ray acceleration needs to be considered in such a young starburst region. Finally, we suggest that the log-normal modeling may be suitable for the X-ray spectral analysis of other giant HII regions, especially when spatially resolved spectroscopy is not practical.

From 22 to 26 June 2020, we hosted ESO's first live e-conference, #H02020, from within ESO headquarters in Garching, Germany. Every day, between 200 and 320 researchers around the globe tuned in to discuss the nature and implications of the discord between precise determinations of the Universe's expansion rate, H0. Originally planned as an in-person meeting, we moved to the virtual domain to maintain strong scientific discourse despite the SARS-CoV-2 (COVID-19) pandemic. Here, we describe our conference setup, participants feedback gathered before and after the meeting, and lessons learned from this unexpected exercise. As e-conferencing will become increasingly common in the future, we provide our perspective on how e-conferences can make scientific exchange more effective and inclusive, in addition to climate friendly.

The long term dreams to approach Mars requires numerous spacecraft attempts for exploration as well as to understand the perception of the red planet. Before launching a mission, the space probe undergoes critical ground testing and effective preparation. Though probes were carefully tested and validated, many experiences temporary or permanent setbacks prior to their final state of mission accomplishment, resulting in the failure of the mission. In order to figure out the problems concerning probe malfunction or failure, we conducted a study on failed Mars probes that are launched between 1960 to 2020. The probes were characterized to determine various modes of failure and their impact on the missions. The results of our study from past probes showed effective integration and testing, sterling fabrication and validation of space probes, adequate software design, feasible recovery options, and novel guidance to probe computers and communication systems.

We present the first 850 $\mu$m polarization observations in the most active star-forming site of the Rosette Molecular Cloud (RMC, $d\sim$1.6 kpc) in the wall of the Rosette Nebula, imaged with the SCUBA-2/POL-2 instruments of the JCMT, as part of the B-Fields In Star-Forming Region Observations 2 (BISTRO-2) survey. From the POL-2 data we find that the polarization fraction decreases with the 850 $\mu$m continuum intensity with $\alpha$ = 0.49 $\pm$ 0.08 in the $p \propto I^{\rm -\alpha}$ relation, which suggests that some fraction of the dust grains remain aligned at high densities. The north of our 850 $\mu$m image reveals a "gemstone ring" morphology, which is a $\sim$1 pc-diameter ring-like structure with extended emission in the "head" to the south-west. We hypothesize that it might have been blown by feedback in its interior, while the B-field is parallel to its circumference in most places. In the south of our SCUBA-2 field the clumps are apparently connected with filaments which follow Infrared Dark Clouds (IRDCs). Here, the POL-2 magnetic field orientations appear bimodal with respect to the large-scale Planck field. The mass of our effective mapped area is $\sim$ 174 $M_\odot$ that we calculate from 850 $\mu$m flux densities. We compare our results with masses from large-scale emission-subtracted Herschel 250 $\mu$m data, and find agreement within 30%. We estimate the POS B-field strength in one typical subregion using the Davis-Chandrasekhar-Fermi (DCF) technique and find 80 $\pm$ 30 $\mu$G toward a clump and its outskirts. The estimated mass-to-flux ratio of $\lambda$ = 2.3 $\pm$ 1.0 suggests that the B-field is not sufficiently strong to prevent gravitational collapse in this subregion.

In this work we explore the new catalog of galactic open clusters that became available recently, containing 1750 clusters that have been re-analysed using the Gaia DR2 catalog to determine the stellar memberships. We used the young open clusters as tracers of spiral arms and determined the spiral pattern rotation speed of the Galaxy and the corotation radius, the strongest Galactic resonance. The sample of open clusters used here increases the last one from Dias et al. (2019) used in the previous determination of the pattern speed by dozens objects. In addition, the distances and ages values are better determined, using improvements to isochrone fitting and including an updated extinction polynomial for the Gaia DR2 photometric band-passes, and the Galactic abundance gradient as a prior for metallicity. In addition to the better age determinations, the catalog contains better positions in the Galactic plane and better proper motions. This allow us to discuss not only the present space distribution of the clusters, but also the space distribution of the clusters's birthplaces, obtained by integration of the orbits for a time equal to their age. The value of the rotation velocity of the arms ($28.5 \pm 1.0$ km s$^{-1}$ kpc$^{-1}$) implies that the corotation radius ($R_c$) is close to the solar Galactic orbit ($R_c/R_0 = 1.01\pm0.08$), which is supported by other observational evidence discussed in this text. A simulation is presented, illustrating the motion of the clusters in the reference frame of corotation. We also present general statistics of the catalog of clusters, like spatial distribution, distribution relative to height from the Galactic plane, and distribution of ages and metallicity. An important feature of the space distribution, the corotation gap in the gas distribution and its consequences for the young clusters, is discussed.

We present a revised SALT2 surface (`SALT2-2021') for fitting the light curves of Type Ia supernovae (SNe Ia), which incorporates new measurements of zero-point calibration offsets and Milky Way extinction. The most notable change in the new surface occurs in the UV region. This new surface alters the distance measurements of SNe~Ia, which can be used to investigate the nature of dark energy by probing the expansion history of the Universe. Using the revised SALT2 surface on public data from the first three years of the Dark Energy Survey Supernova Program (combined with an external low-$z$ SNe Ia sample) and combining with cosmic microwave background constraints, we find a change in the dark energy equation of state parameter, $\Delta w = 0.015 \pm 0.004$. This result highlights the continued importance of controlling and reducing systematic uncertainties, particularly with the next generation of supernova analyses aiming to improve constraints on dark energy properties.

Modeling collisionless magnetic reconnection rate is an outstanding challenge in basic plasma physics research. While the seemingly universal rate of an order $\mathcal{O}(0.1)$ is often reported in the low-$\beta$ regime, it is not clear how reconnection rate scales with a higher plasma $\beta$. Due to the complexity of the pressure tensor, the available reconnection rate model is limited to the low plasma-$\beta$ regime, where the thermal pressure is arguably negligible. However, the thermal pressure effect becomes important when $\beta \gtrsim \mathcal{O}(1)$. Using first-principle kinetic simulations, we show that both the reconnection rate and outflow speed drop as $\beta$ gets larger. A simple analytical framework is derived to take account of the self-generated pressure anisotropy and pressure gradient in the force-balance around the diffusion region, explaining the varying trend of key quantities and reconnection rates in these simulations with different $\beta$. The predicted scaling of the normalized reconnection rate is $\simeq \mathcal{O}(0.1/\sqrt{\beta_{i0}})$ in the high $\beta$ limit, where $\beta_{i0}$ is the ion $\beta$ of the inflow plasma.

Expanding on our prior efforts to search for Lorentz invariance violation (LIV) using linear optical polarimetry of extragalactic objects, we propose a new method that combines linear and circular polarization measurements. While existing work focused on the tendency of LIV to reduce the linear polarization degree, this new method additionally takes into account the coupling between photon helicities induced by some models. This coupling can generate circular polarization as light propagates, even if there is no circular polarization at the source. Combining significant detections of linear polarization of light from extragalactic objects with the absence of a detection of circular polarization in most measurements results in significantly tighter constraints on LIV. The analysis is carried out in the framework of the Standard-Model Extension (SME), an effective field theory framework to describe low-energy effects of an underlying more fundamental quantum gravity theory. We evaluate the performance of our method by deriving constraints on mass dimension $d=4$ CPT-even SME coefficients from a small set of archival circular and linear optical polarimetry and compare them to similar constraints derived in previous work with far larger sample sizes and based on linear polarimetry only. The new method yields constraints that are an order of magnitude tighter even for our modest sample size of 21 objects. Based on the demonstrated gain in constraining power from scarce circular data, we advocate for the need in future extragalactic circular polarization surveys.

In this study, a novel type of Fourier transform radio spectrometer (termed as all-digital radio spectrometer; ADRS) has been developed in which all functionalities comprising a radio spectrometer including a sampler and Fourier computing unit were implemented as a soft-core on a field-programmable gate array (FPGA). A delay-line-based ramp-compare analog-to-digital converter (ADC), one of completely digital ADC, was used, and two primary elements of the ADC, an analog-to-time converter (ATC) and a time-to-digital converter (TDC), were implemented on the FPGA. The sampling rate of the ADRS $f$ and the quantization bit rate $n$ are limited by the relation, $\tau = \frac{1}{2^{n}f}$, where $\tau$ is the latency of the delay element of the delay-line. Given that the typical latency of the delay element implemented on FPGAs is $\sim10$ ps, adoption of a low quantization bit rate, which satisfies the requirements for radio astronomy, facilitates the realization of a high sampling rate up to $\sim$100 GSa/s. In addition, as the proposed \ADRS does not require a discrete ADC and can be implemented on mass-produced evaluation boards, its fabrication cost is much lower than that of conventional spectrometers. The ADRS prototype was fabricated with values of $f$ = 600 MSa/s and $n$ = 6.6 using a PYNQ-Z1 evaluation board, with a $\tau$ of 16.7 ps. The performance of the prototype, including its linearity and stability, was measured, and a test observation was conducted using the Osaka Prefecture University 1.85-m mm-submm telescope; this confirmed the potential application of the prototype in authentic radio observations. With 10 times better cost performance ($\sim$800 USD GHz$^{-1}$) than conventional radio spectrometers, the prototype facilitates cost-effective coverage of intermediate frequency (IF) bandwidths of $\sim100$ GHz in modern receiver systems.

Methods. We use both a Monte Carlo (MC) sampling and a Bayesian analysis to find the effects of the photospheric composition and the touchdown radius. We apply these two methods to six LMXBs that show PRE XRBs. In both methods, we solve the Eddington flux equation and the apparent angular area equation both of which include the correction terms. For the MC sampling, we have developed an iterative method in order to solve these two equations more efficiently. Results. We confirm that the effects of the photospheric composition and the touchdown radius are similar in the statistical and analytical estimation of mass and radius even when the correction terms are considered. Furthermore, in all of the six sources, we find that a H-poor photosphere and a large touchdown radius are favored statistically regardless of the statistical method. Our Bayesian analysis also hints that touchdown can occur farther from the neutron star surface when the photosphere is more H-poor. This correlation could be qualitatively understood with the Eddington flux equation. We propose a physical explanation for this correlation between the photospheric composition and the touchdown radius. Our results show that when accounting for the uncertainties of the photospheric composition and the touchdown radius, most likely radii of the neutron stars in these six LMXBs are less than 12.5 km, which is similar to the bounds for the neutron star radius placed with the tidal deformability measured from the gravitational wave signal.

Stellar, substellar, and planetary atmosphere models are all highly sensitive to the input opacities. Generational differences between various state-of-the-art stellar/planetary models are primarily because of incomplete and outdated atomic/molecular line-lists. Here we present a database of pre-computed absorption cross-sections for all isotopologues of key atmospheric molecules relevant to late-type stellar, brown dwarf, and planetary atmospheres: MgH, AlH, CaH, TiH, CrH, FeH, SiO, TiO, VO, and H2O. The pressure and temperature ranges of the computed opacities are between 10$^{-6}$--3000~bar and 75--4000~K, and their spectral ranges are 0.25--330~$\mu$m for many cases where possible. For cases with no pressure-broadening data, we use collision theory to bridge the gap. We also probe the effect of absorption cross-sections calculated from different line lists in the context of Ultra-Hot Jupiter and M-dwarf atmospheres. Using 1-D self-consistent radiative-convective thermochemical equilibrium models, we report significant variations in the theoretical spectra and thermal profiles of substellar atmospheres. With a 2000 K representative Ultra-Hot Jupiter, we report variations of up to 320 and 80 ppm in transmission and thermal emission spectra, respectively. For a 3000 K M-dwarf, we find differences of up to 125$\%$ in the spectra. We find that the most significant differences arise due to the choice of TiO line-lists, primarily below 1$\mu$m. In sum, we present (1) a database of pre-computed molecular absorption cross-sections, and (2) quantify biases that arise when characterizing substellar/exoplanet atmospheres due to line list differences, therefore highlighting the importance of correct and complete opacities for eventual applications to high precision spectroscopy and photometry.

We present evidence that multiple accretion events are required to explain the origin of the $Gaia$-Sausage and Enceladus (GSE) structures, based on an analysis of dynamical properties of main-sequence stars from the Sloan Digital Sky Survey Data Release 12 and $Gaia$ Data Release 2. GSE members are selected to have eccentricity ($e$) $>$ 0.7 and [Fe/H] $<$ -1.0, and separated into low and high orbital-inclination (LOI/HOI) groups. We find that the LOI stars mainly have $e < 0.9$ and are clearly separable into two groups with prograde and retrograde motions. The LOI stars exhibit prograde motions in the inner-halo region and strong retrograde motions in the outer-halo region. We interpret the LOI stars in these regions to be stars accreted from two massive dwarf galaxies with low-inclination prograde and retrograde orbits, affected to different extents by dynamical friction due to their different orbital directions. In contrast, the majority of the HOI stars have $e > 0.9$, and exhibit a globally symmetric distribution of rotational velocities ($V_{\rm \phi}$) near zero, although there is evidence for a small retrograde motion for these stars ($V_{\rm \phi}$ $\sim$ -15 $\rm{km~s^{-1}}$) in the outer-halo region. We consider these stars to be stripped from a massive dwarf galaxy on a high-inclination orbit. We also find that the LOI and HOI stars on highly eccentric and tangential orbits with clear retrograde motions exhibit different metallicity peaks at [Fe/H] = -1.7 and -1.9, respectively, and argue that they are associated with two low-mass dwarf galaxies accreted in the outer-halo region of the Galaxy.

We investigate the gravitational microlensing event KMT-2019-BLG-1715, of which light curve shows two short-term anomalies from a caustic-crossing binary-lensing light curve: one with a large deviation and the other with a small deviation. We identify five pairs of solutions, in which the anomalies are explained by adding an extra lens or source component in addition to the base binary-lens model. We resolve the degeneracies by applying a method, in which the measured flux ratio between the first and second source stars is compared with the flux ratio deduced from the ratio of the source radii. Applying this method leaves a single pair of viable solutions, in both of which the major anomaly is generated by a planetary-mass third body of the lens, and the minor anomaly is generated by a faint second source. A Bayesian analysis indicates that the lens comprises three masses: a planet-mass object with $\sim 2.6~M_{\rm J}$ and binary stars of K and M dwarfs lying in the galactic disk. We point out the possibility that the lens is the blend, and this can be verified by conducting high-resolution followup imaging for the resolution of the lens from the source.

Mt. Abu Faint Object Spectrograph and Camera - Pathfinder (MFOSC-P) is an imager-spectrograph developed for the Physical Research Laboratory (PRL) 1.2m telescope at Gurushikhar, Mt. Abu, India. MFOSC-P is based on a focal reducer concept and provides seeing limited imaging (with a sampling of 3.3 pixels per arc-second) in Bessell's B, V, R, I and narrow-band H-$\alpha$ filters. The instrument uses three plane reflection gratings, covering the spectral range of 4500-8500$\AA$, with three different resolutions of 500, 1000, and 2000 around their central wavelengths. MFOSC-P was conceived as a pathfinder instrument for a next-generation instrument on the PRL's 2.5m telescope which is coming up at Mt. Abu. The instrument was developed during 2015-2019 and successfully commissioned on the PRL 1.2m telescope in February 2019. The designed performance has been verified with laboratory characterization tests and on-sky commissioning observations. Different science programs covering a range of objects are being executed with MFOSC-P since then, e.g., spectroscopy of M-dwarfs, novae $\&$ symbiotic systems, and detection of H-$\alpha$ emission in star-forming regions. MFOSC-P presents a novel design and cost-effective way to develop a FOSC (Faint Object Spectrograph and Camera) type of instrument on a shorter time-scale of development. The design and development methodology presented here is most suitable in helping the small aperture telescope community develop such a versatile instrument, thereby diversifying the science programs of such observatories.

Here we reprise an initial orbit determination method introduced by O. F. Mossotti \cite{mossotti} employing four geocentric sky-plane observations and linear equations to compute the angular momentum of the observed body. We then extend the method to topocentric observations, yielding quadratic equations for the angular momentum. The performance of the two versions are compared through numerical tests with synthetic asteroid data using different time intervals between consecutive observations and different astrometric errors.

Classical novae are among the most frequent transient events in the Milky Way, and key agents of ongoing nucleosynthesis. Despite their large numbers, they have never been observed in soft $\gamma$-ray emission. Measurements of their $\gamma$-ray signatures would provide both, insights on explosion mechanism as well as nucleosynthesis products. Our goal is to constrain the ejecta masses of $\mathrm{^7Be}$ and $\mathrm{^{22}Na}$ from classical novae through their $\gamma$-ray line emissions at 478 and 1275 keV. We extract posterior distributions on the line fluxes from archival data of the INTEGRAL/SPI spectrometer telescope. We then use a Bayesian hierarchical model to link individual objects and diffuse emission and infer ejecta masses from the whole population of classical novae in the Galaxy. Individual novae are too dim to be detectable in soft $\gamma$-rays, and the upper bounds on their flux and ejecta mass uncertainties cover several orders of magnitude. Within the framework of our hierarchical model, we can, nevertheless, infer tight upper bounds on the $\mathrm{^{22}Na}$ ejecta masses, given all uncertainties from individual objects as well as diffuse emission, of $<2.0 \times 10^{-7}\,\mathrm{M_{\odot}}$ (99.85th percentile). In the context of ONe nucleosynthesis, the $\mathrm{^{22}Na}$ bounds are consistent with theoretical expectations, and exclude that most ONe novae happen on white dwarfs with masses around $1.35\,\mathrm{M_{\odot}}$. The upper bounds from $\mathrm{^{7}Be}$ are uninformative. From the combined ejecta mass estimate of $\mathrm{^{22}Na}$ and its $\beta^+$-decay, we infer a positron production rate of $<5.5 \times 10^{42}\,\mathrm{e^+\,s^{-1}}$, which would make at most 10% of the total annihilation rate in the Milky Way.

In this paper we present recent Low Frequency Array (LOFAR) observations of the pulsar J1740+1000. We confirm that its spectrum has a turnover at 260 MHz, which is unusual for a typical pulsar. We argue that in this case interferometric imaging provides more accurate pulsar flux estimates than other, more traditional, means such as beamformed observations. We conclude that existing calibration and imaging techniques can be used for a more comprehensive study of the influence of the interstellar medium on the point-like sources at very low frequencies in the near future.

We report here results on the analysis of correlated flux variations between the optical and GeV $\gamma$-ray bands in three bright BL Lac objects, namely AO\, 0235+164, OJ 287 and PKS 2155$-$304. This was based on the analysis of about 10 years of data from the {\it Fermi} Gamma-ray Space Telescope covering the period between 08 August 2008 to 08 August 2018 along with optical data covering the same period. For all the sources, during the flares analysed in this work, the optical and $\gamma$-ray flux variations are found to be closely correlated. From broad band spectral energy distribution modelling of different epochs in these sources using the one zone leptonic emission model, we found that the optical-UV emission is dominated by synchrotron emission from the jet. The $\gamma$-ray emission in the low synchrotron peaked sources AO\, 0235+164 and OJ 287 are found to be well fit with external Compton (EC) component, while, the $\gamma$-ray emission in the high synchrotron peaked source PKS 2155$-$304 is well fit with synchrotron self Compton component. Further we note that the $\gamma$-ray emission during the high flux state of AO 0235+164 (epochs A and B) requires seed photons from both the dusty torus and broad line region, while the $\gamma$-ray emission in OJ 287 and during epochs C and D of AO\,0235+164 can be modelled by EC scattering of infra-red photons from the torus.

Winged radio sources are a small sub-class of extragalactic radio sources which display a pair of low surface brightness radio lobes, known as `wings' aligned at a certain angle with the primary jets. Depending on the location of wings, these galaxies look like X or Z and are known as X-shaped Radio Galaxy (XRG) or Z-shaped Radio Galaxy (ZRG). We report the identification of 33 winged radio sources from the LOFAR Two-metre Sky Survey First Data Release (LoTSS DR1) out of which 21 sources are identified as X-shaped radio galaxies and 12 as Z-shaped radio galaxies. Optical counterparts are identified for 14 XRGs (67 per cent) and 12 ZRGs (100 per cent). We studied various physical parameters of these sources like spectral index, radio luminosity, and power. The radio spectrum of the majority of XRGs and ZRGs is steep ($\alpha_{1400}^{144}>0.5$), which is typical of lobe dominated radio galaxies. The statistical studies are done on the relative size of the major and minor axes and the angle between the major axis and minor axis for XRGs.

The small but measurable effect of weak gravitational lensing on the cosmic microwave background radiation provide information about the large-scale distribution of matter in the universe. We use the all sky distribution of matter, as represented by the {\em convergence map} that is inferred from CMB lensing measurement by Planck survey, to test the fundamental assumption of Statistical Isotropy (SI) of the universe. For the analysis we use the $\alpha$ statistic that is devised from the contour Minkowski tensor, a tensorial generalization of the scalar Minkowski functional, the contour length. In essence, the $\alpha$ statistic captures the ellipticity of isofield contours at any chosen threshold value of a smooth random field and provides a measure of anisotropy. The SI of the observed convergence map is tested against the suite of realistic simulations of the convergence map provided by the Planck collaboration. We first carry out a global analysis using the full sky data after applying the galactic and point sources mask. We find that the observed data is consistent with SI. Further we carry out a local search for departure from SI in small patches of the sky using $\alpha$. This analysis reveals several sky patches which exhibit deviations from simulations with statistical significance higher than 95\% confidence level (CL). Our analysis indicates that the source of the anomalous behaviour of most of the outlier patches is inaccurate estimation of noise. We identify two outlier patches which exhibit anomalous behaviour originating from departure from SI at higher than 95\% CL. Most of the anomalous patches are found to be located roughly along the ecliptic plane or in proximity to the ecliptic poles.

Accurate component separation of full-sky maps in the radio and microwave frequencies, such as the cosmic microwave background (CMB), relies on a thorough understanding of the statistical properties of the Galactic foreground emissions. Using scalar Minkowski functionals and their tensorial generalization known as Minkowski tensors, we analyze the statistical properties of one of the major foreground components, namely the Galactic synchrotron given by the full sky 408 MHz Haslam map. We focus on understanding the nature of non-Gaussianity and statistical isotropy of the cooler regions of the map as a function of angular scale. We find that the overall level of the non-Gaussian deviations does decrease as more high emission regions are masked and as we go down to smaller scales, in agreement with the results obtained in earlier works. However, they remain significantly high, of order 3.3$\sigma$, at the smallest angular scales relevant for the Haslam map. We carry out a detailed examination of the non-Gaussian nature using the generalized skewness and kurtosis cumulants that arise in the perturbative expansion of Minkowski functionals for weakly non-Gaussian fields. We find that the leading sources of non-Gaussianity are the kurtosis terms which are considerably larger than the skewness terms at all angular scales. Further, for the cooler regions of the Haslam map, we find that the non-Gaussian deviations of the Minkowski functionals can be well explained by the perturbative expansion up to second-order (up to kurtosis terms), with first-order terms being sub-dominant. Lastly, we test the statistical isotropy of the Haslam map and find that it becomes increasingly more isotropic at smaller scales.

Aim: To study the multimessenger nature of the signal that can result from the phase transition of a neutron star to a quark star and their corresponding astrophysical observations. Methods: The phase transition process is initiated by the abrupt pressure and density changes at the star center, giving rise to a shock which deconfines matter followed by a weak front converting excess down to strange quarks to attain absolute stability. This process's effects are investigated by understanding how the energy escapes from the star in the form of neutrino-antineutrino annihilation. For such annihilation process, the corresponding energy deposition rate is calculated. Structural changes due to the energy loss have been investigated in the likes of misalignment angle evolution of the star and its astrophysical observation through gravitational waves. Results: The energy and time signature for the neutrino-antineutrino annihilation is compared with the observed isotropic energy for a short gamma-ray burst. The misalignment angle evolves to align the star's tilt axis, which can lead to the sudden increase or decrease of radio intensity from the pulsar. The corresponding gravitational wave emission, both continuous and burst, all lead towards multimessenger signals coming from the phase transition.

Planets are thought to form at the early stage of stellar evolution when the mass accretion is still ongoing. RY Tau is a T Tauri type star at the age of a few Myr, with accretion disc seen at high inclination, so that line of sight crosses both the wind and the accretion gas flows. In a long series of spectroscopic monitoring of the star in 2013-2020, we detected variations in H-alpha and NaI D absorptions at radial velocities of infall (accretion) and outflow (wind) with a period of about 22 days. The absorptions in the infalling and the outflowing gas streams vary in anti-phase: an increase of infall is accompanied by a decrease of outflow, and vice versa. These flip-flop oscillations retain phase over several years of observations. We suggest that this may result from the MHD processes at the disk-magnetosphere boundary in the propeller mode. Another possibility is that a massive planet modulates some processes in the disc and provides the observed effects. The period, if Keplerian, corresponds to a distance of 0.2 AU, which is close to the dust sublimation radius in this star. The presence of the putative planet may be confirmed by radial velocity measurements: expected amplitude is > 90 m/s if a planet mass is > 2 Mj.

We investigate of the systems of proper motions of stars in the ground-based catalogues HSOY, UCAC5, GPS1 and PMA derived by combining with the Gaia DR1 space data. Assuming the systematic differences of stellar proper motions of two catalogues to be caused by the mutual solid-body rotation and glide of the coordinate systems produced by the data of the catalogues under comparison, we analyse the components of the mutual rotation vector and displacement of the origins of these systems. The equatorial components of the vector of mutual rotation velocity of the compared coordinate systems, as well as velocities of the mutual displacement of their origins, varying within the range from 0.2 to 2.9 mas/yr, were derived from a comparison of proper motions of the sources that are common for Gaia EDR3 and the TGAS, UCAC5, HSOY, GPS1 and PMA catalogues, respectively. The systematic errors of proper motions of stars in the HSOY, GPS1, PMA and Gaia EDR3 catalogues in the range of faint stellar magnitudes were estimated by analysing the formal proper motions of extragalactic objects contained in these catalogues. The coordinate system realised by the Gaia EDR3 data at the level of < 0.1 mas/yr is shown to have no rotation and glide relative to the LQAC-5, ALLWISEAGN, Milliquas extragalactis sources within the range from 15 to 21 stellar G magnitude. Among the ground-based catalogues, the system of proper motions of the PMA stars, which is independent of the Gaia EDR3 data, is the closest to the Gaia EDR3 system of proper motions in G magnitude range from 15 to 21.

The Radar Echo Telescope for Cosmic Rays (RET-CR) is a recently initiated experiment designed to detect the englacial cascade of a cosmic-ray initiated air shower via in-ice radar, toward the goal of a full-scale, next-generation experiment to detect ultra high energy neutrinos in polar ice. For cosmic rays with a primary energy greater than 10 PeV, roughly 10% of an air-shower's energy reaches the surface of a high elevation ice-sheet ($\gtrsim$2 km) concentrated into a radius of roughly 10 cm. This penetrating shower core creates an in-ice cascade many orders of magnitude more dense than the preceding in-air cascade. This dense cascade can be detected via the radar echo technique, where transmitted radio is reflected from the ionization deposit left in the wake of the cascade. RET-CR will test the radar echo method in nature, with the in-ice cascade of a cosmic-ray initiated air-shower serving as a test beam. We present the projected event rate and sensitivity based upon a three part simulation using CORSIKA, GEANT4, and RadioScatter. RET-CR expects $\sim$1 radar echo event per day.

We use N-body/hydrodynamic simulations to study the evolution of the spin of a Milky Way-like galaxy through interactions. We perform a controlled experiment of co-planner galaxy-galaxy encounters and study the evolution of disk spins of interacting galaxies. Specifically, we consider the cases where the late-type target galaxy encounters an equally massive companion galaxy, which has either a late or an early-type morphology, with the closest approach distance of about 50 kpc, in prograde or retrograde sense. By examining the time change of the circular velocity of the disk material of the target galaxy from each case, we find that the target galaxy tends to lose the spin through prograde collisions but hardly through retrograde collisions, regardless of the companion galaxy type. The decrease of the spin results mainly from the deflection of the orbit of the disk material by tidal disruption. Although there is some disk material which gains the circular velocity through hydrodynamic as well as gravitational interactions or by transferring material from the companion galaxy, it turns out that the amount of the material is generally insufficient to increase the overall galactic spin under the conditions we set. It is found that the spin angular momentum of the target galaxy disk decreases by 15 - 20% after a prograde collision. We conclude that the accumulated effects of galaxy-galaxy interactions will play an important role in determining the total angular momentum of late-type galaxies at current stage.

The products of magnetic reconnection in Saturn's magnetotail are identified in magnetometer observations primarily through characteristic deviations in the north-south component of the magnetic field. These magnetic deflections are caused by travelling plasma structures created during reconnection rapidly passing over the observing spacecraft. Identification of these signatures have long been performed by eye, and more recently through semi-automated methods, however these methods are often limited through a required human verification step. Here, we present a fully automated, supervised learning, feed forward neural network model to identify evidence of reconnection in the Kronian magnetosphere with the three magnetic field components observed by the Cassini spacecraft in Kronocentric radial-theta-phi (KRTP) coordinates as input. This model is constructed from a catalogue of reconnection events which covers three years of observations with a total of 2093 classified events, categorized into plasmoids, travelling compression regions and dipolarizations. This neural network model is capable of rapidly identifying reconnection events in large time-span Cassini datasets, tested against the full year 2010 with a high level of accuracy (87%), true skill score (0.76), and Heidke skill score (0.73). From this model, a full cataloguing and examination of magnetic reconnection events in the Kronian magnetosphere across Cassini's near Saturn lifetime is now possible.

In this work, retrograde mean motion resonances (MMRs) are investigated by means of analytical and numerical approaches. Initially, we define a new resonant angle to describe the retrograde MMRs and then perform a series of canonical transformations to formulate the resonant model, in which the phase portrait, resonant centre and resonant width can be analytically determined. To validate the analytical developments, the non-perturbative analysis is made by taking advantage of Poincar\'e surfaces of section. Some modifications are introduced in the production of Poincar\'e sections and, in particular, it becomes possible to make direct comparisons between the analytical and numerical results. It is found that there exists an excellent correspondence between the phase portraits and the associated Poincar\'e sections, and the analytical results agree well with the numerical results in terms of the resonant width and the location of resonant centre. Finally, the numerical approach is utilized to determine the resonant widths and resonant centres over the full range of eccentricity. In particular, seven known examples of retrograde asteroids including 2015 BZ509, 2008 SO218, 1999 LE31, 2000 DG8, 2014 AT28, 2016 LS and 2016 JK24 are found inside the libration zones of retrograde MMRs with Jupiter. The results obtained in this work may be helpful for understanding the dynamical evolution for asteroids inside retrograde MMRs.

The winter polar vortices on Mars are annular in terms of their potential vorticity (PV) structure, a phenomenon identified in observations, reanalysis and some numerical simulations. Some recent modeling studies have proposed that condensation of atmospheric carbon dioxide at the winter pole is a contributing factor to maintaining the annulus through the release of latent heat. Dust and topographic forcing are also known to be causes of internal and interannual variability in the polar vortices. However, coupling between these factors remains uncertain, and previous studies of their impact on vortex structure and variability have been largely limited to a single Martian global climate model (MGCM). Here, by further developing a novel MGCM, we decompose the relative roles of latent heat and dust as drivers for the variability and structure of the northern Martian polar vortex. We also consider how Martian topography modifies the driving response. By also analyzing a reanalysis dataset we show that there is significant dependence in the polar vortex structure and variability on the observations assimilated. In both model and reanalysis, high atmospheric dust loading (such as that seen during a global dust storm) can disrupt the vortex, cause the destruction of PV in the low-mid altitudes (> 0.1 hPa), and significantly reduce spatial and temporal vortex variability. Through our simulations, we find that the combination of dust and topography primarily drives the eddy activity throughout the Martian year, and that although latent heat release can produce an annular vortex, it has a relatively minor effect on vortex variability.

In the dynamical models of gamma-ray burst (GRB) afterglows, the uniform assumption of the shocked region is known as provoking total energy conservation problem. In this work we consider shocks originating from magnetized ejecta, extend the energy-conserving hydrodynamical model of Yan et al. (2007) to the MHD limit by applying the magnetized jump conditions from Zhang & Kobayashi (2005). Compared with the non-conservative models, our Lorentz factor of the whole shocked region is larger by a factor $\lesssim\sqrt{2}$. The total pressure of the forward shocked region is higher than the reversed shocked region, in the relativistic regime with a factor of about 3 in our interstellar medium (ISM) cases while ejecta magnetization degree $\sigma<1$, and a factor of about 2.4 in the wind cases. For $\sigma\le 1$, the non-conservative model loses $32-42$% of its total energy for ISM cases, and for wind cases $25-38$%, which happens specifically in the forward shocked region, making the shock synchrotron emission from the forward shock less luminous than expected. Once the energy conservation problem is fixed, the late time light curves from the forward shock become nearly independent of the ejecta magnetization. The reverse shocked region doesn't suffer from the energy conservation problem since the changes of the Lorentz factor are recompensed by the changes of the shocked particle number density. The early light curves from the reverse shock are sensitive to the magnetization of the ejecta, thus are an important probe of the magnetization degree.

In the currently accepted model for cosmic baryon evolution, Cosmic Dawn and the Epoch of Reionization are significant times when first light from the first luminous objects emerged, transformed and subsequently ionized the primordial gas. The 21 cm hyperfine transition of neutral hydrogen, redshifted from these cosmic times to a frequency range of 40 to 200 MHz, has been recognized as an important probe of the physics of CD/EoR. The global 21-cm signal is predicted to be a spectral distortion of a few 10's to a few 100's of mK, which is expected to be present in the cosmic radio background as a trace additive component. SARAS, Shaped Antenna measurement of the background RAdio Spectrum, is a spectral radiometer purpose designed to detect the weak 21-cm signal from CD/EoR. An important subsystem of the radiometer, the digital correlation spectrometer, is developed around a high speed digital signal processing platform called pSPEC. pSPEC is built around two quad 10 bit analog-to-digital converters and a Virtex 6 field programmable gate array, with provision for multiple Gigabit Ethernet and 4.5 Gbps fibre optic interfaces. Here we describe the system design of the digital spectrometer, the pSPEC board, and the adaptation of pSPEC to implement a high spectral resolution of about 61 kHz, high dynamic range correlation spectrometer covering the entire CD/EoR band. As the SARAS radiometer is required to be deployed in remote locations where terrestrial radio frequency interference is a minimum, the spectrometer is designed to be compact, portable and operating off internal batteries. The paper includes an evaluation of the spectrometer's susceptibility to radio frequency interference and capability to detect signals from CD/EoR.

The relevance of non-parametric reconstructions of cosmological functions lies in the possibility of analyzing the observational data independently of any theoretical model. Several techniques exist and, recently, Artificial Neural Networks have been incorporated to this type of analysis. By using Artificial Neural Networks we present a new strategy to perform non-parametric data reconstructions without any preliminary statistical or theoretical assumptions and even for small observational datasets. In particular, we reconstruct cosmological observables from cosmic chronometers, $f\sigma_8$ measurements and the distance modulus of the Type Ia supernovae. In addition, we introduce a first approach to generate synthetic covariance matrices through a variational autoencoder, for which we employ the covariance matrix of the Type Ia supernovae compilation. To test the usefulness of our developed methods, with the neural network models we generated random data points mostly absent in the original datasets and performed a Bayesian analysis on some simple dark energy models. Some of our findings point out to slight deviations from the $\Lambda$CDM standard model, contrary to the expected values coming from the original datasets.

We present an analysis of the kinematics of fourteen satellites of the Milky Way (MW). We use proper motions (PMs) from the $Gaia$ Early Data Release 3 (EDR3) and line-of-sight velocities ($v_{\mathrm{los}}$) available in the literature to derive the systemic 3D motion of these systems. For six of them, namely the Carina, Draco, Fornax, Sculptor, Sextans, and Ursa Minor dwarf spheroidal galaxies (dSph), we study the internal kinematics projecting the stellar PMs into radial, $V_R$ (expansion/contraction), and tangential, $V_T$ (rotation), velocity components with respect to the centre-of-mass. We find significant rotation in the Carina ($|V_T| = 9.6 \pm 4.5 \ {\rm{km \ s^{-1}}}\>$), Fornax ($|V_T| = 2.8 \pm 1.3 \ {\rm{km \ s^{-1}}}\>$), and Sculptor ($|V_T| = 3.0 \pm 1.0 \ {\rm{km \ s^{-1}}}\>$) dSphs. Besides the Sagittarius dSph, these are the first measurements of internal rotation in the plane of the sky in the MW's classical dSphs. All galaxies except Carina show $|V_T| / \sigma_v < 1$. We find that slower rotators tend to show, on average, larger sky-projected ellipticity (as expected for a sample with random viewing angles), and are located at smaller Galactocentric distances (as expected for tidal stirring scenarios in which rotation is transformed into random motions as satellites sink into the parent halo). However, these trends are small and not statistically significant, indicating that rotation has not played a dominant role in shaping the 3D structure of these galaxies. Either tidal stirring had a weak impact on the evolution of these systems, or it perturbed them with similar efficiency regardless of their current Galactocentric distance.

Nearby massive compact elliptical galaxies (MCEGs) are strong candidates for relic galaxies (i.e. local analogs of red nuggets at high redshifts). It is expected that the globular cluster (GC) systems of relic galaxies are dominated by red (metal-rich) GCs. NGC 1277 is known as a unique example of such a galaxy in the previous study. In this study, we search for GCs in 12 nearby MCEGs at distances of $\lesssim 100$ Mpc from the Hubble Space Telescope/Wide Field Camera 3 F814W($I_{814}$)/F160W($H_{160}$) archival images. We find that most of these MCEGs host a rich population of GCs with a color range of $0.0<(I_{814}-H_{160})_0<1.1$. The fractions of their red GCs range from $f_{RGC} =0.2$ to 0.7 with a mean of $f_{RGC} =0.48\pm0.14$. We divide the MCEG sample into two groups: one in clusters and the other in groups/fields. The mean red GC fraction of the cluster MCEGs is $0.60\pm0.06$, which is 0.2 larger than the value of the group/field MCEGs, $0.40\pm0.10$. The value for the cluster MCEGs is $\sim$0.3 larger than the mean value of giant early-type galaxies with similar stellar mass in the Virgo cluster ($f_{RGC} =0.33\pm0.13$). Our results show that most of the MCEGs in our sample are indeed relic galaxies. This further implies that a majority of the red GCs in MCEGs are formed early in massive galaxies and that most MCEGs in the local universe have rarely undergone mergers after they became red nuggets about 10 Gyr ago.

The coupling between the angular momentum of a compact object and an external tidal field gives rise to the "rotational" tidal Love numbers, which affect the tidal deformability of a spinning self-gravitating body and enter the gravitational waveform of a binary inspiral at high post-Newtonian order. We provide numerical evidence for a surprising "hidden" symmetry among the rotational tidal Love numbers with opposite parities, which are associated to perturbations belonging to separate sectors. This symmetry, whose existence had been suggested on the basis of a Lagrangian description of the tidal interaction in a binary system, holds independently of the equation of state of the star.

The manuscript New Theory of Sunspots was identified and attributed to the Jesuit astronomer and hydrographer from Marseille Esprit Pezenas (1692-1776) during the year 2002, during the systematic search for papers concerning him in the various French funds. The manuscript, in large part unpublished, New Theory of Sunspots was written by astronomer and professor of Marseilles Jesuit hydrography, Father Esprit Pezenas (1692-1776). The text was composed and revised between the years 1766 and 1772. In this manuscript, of which we are going to study the conditions of its composition, P. Pezenas gives one of the last geometric methods, related graphical methods, allowing the inclination of the Sun's equator to be deduced on ecliptic, using three observations of a sunspot, at given times, and in taking into account the movement of the Earth with respect to the Sun during the duration of observations. This text contains all the digital elements allowing to follow and understand the application of geometric processing methods of observations and observation techniques mentioned above in the introduction: passages of edges and Sun spots at the meridian, use of objective micrometers and wire micrometers. At detour of his calculations, Father Pezenas announces a duration of about 26 days and 9 hours for the rotation of the Sun around its axis, comments on and corrects some observation data published in Le Monnier's Heavenly History.

We study the dynamical Chern-Simons gravity as an effective quantum field theory, and discuss a broad range of its parameter space where the theory is valid. Within that validity range, we show that slowly rotating black holes acquire novel geometric structures due to the gravitational dynamical Chern-Simons term. In particular, the rotating black hole solutions get endowed with two, cap-like domains, emanating from the north and south poles in the standard Boyer-Lindquist coordinates. The domains extend out to a distance that is approximately a few percent of the black hole's size. The cap-like domains have an unusual equation of state, pointing to non-standard dynamics within the caps. In particular, the focusing condition for geodesics is violated in those domains. This in turn implies that the Hawking-Penrose singularity theorem cannot be straightforwardly applied to hypothetical probe matter placed within the Chern-Simons caps.

The XENON collaboration recently reported an excess of electron recoil events in the low energy region with a significance of around $3.3\sigma$. An explanation of this excess in terms of thermal dark matter seems challenging. We propose a scenario where dark matter in the Milky Way halo gets boosted as a result of scattering with the diffuse supernova neutrino background. This interaction can accelerate the dark-matter to semi-relativistic velocities, and this flux, in turn, can scatter with the electrons in the detector, thereby providing a much better fit to the data. We identify regions in the parameter space of dark-matter mass and interaction cross-section which satisfy the excess. Furthermore, considering the data only hypothesis, we also impose bounds on the dark-matter scattering cross-section, which are competitive with bounds from other experiments.

Just a few decades after the discovery of the Charon Relay, and the ensuing First Contact War, relatively little is known about the population of planets linked by the Prothean mass relays. Understanding the nature of these systems and how they may differ from the broader population of planetary systems in our galaxy is key to both continued human habitation across the broader Galaxy, as well as to our understanding of the Prothean civilization. What factors motivated their choices of planetary systems? Characterizing these systems allows us to peer into Prothean society and culture, and make inferences about the preferences that drove their expansion throughout the Galaxy. In this study, we undertake a broad analysis of the systems recorded in the Systems Alliance Planetary Survey, examining their dynamical stability, orbital properties, and the climates of the inhabited worlds. We find that the Alliance data is inconsistent with both a modern understanding of planetary system dynamics, as well as with our understanding of Earth-like climate dynamics. We suggest this is due in part to security-related data obfuscation by the Alliance, and in part due to the real preferences of the Protheans.

It is shown that ions can be accelerated to MeV energy range in the direction perpendicular to the magnetic field by the ExB mechanism of electrostatic waves. The acceleration occurs in discrete steps of duration being a small fraction the gyroperiod and can explain observations of ion energization to 10 keV at quasi-perpendicular shocks and to 100-1000 keV at quasi-parallel shocks. A general expression is provided for the maximum energy of ions accelerated in shocks of arbitrary configuration. The waves involved in the acceleration are related to three cross-field current-driven instabilities: the lower hybrid drift (LHD) instability induced by the density gradients in shocks and shocklets, followed by the modified two-stream (MTS) and electron cyclotron drift (ECD) instabilities, induced by the ExB drift of electrons in the strong LHD wave electric field. The ExB wave mechanism accelerates heavy ions to energies proportional to the atomic mass number, which is consistent with satellite observations upstream of the bow shock and also with observations of post-shocks in supernovae remnants.

Considering the finite IR behavior of quantum chromodynamics (QCD) running coupling constant in some experiments, we intend to investigate different models presenting running coupling with finite values in the IR region. Using analytic and background perturbation theories, we obtain some equation of states (EoSs) of strange quark matter which satisfy necessary conditions of suitable EoSs. Then we evaluate the properties of strange quark stars in general relativity. Our results for maximum gravitational mass is comparable with the recent LIGO data for the compact binary merger, GW190425.

Gravitational counterpart of the chiral magnetic effect, which is referred as the chiral gravitational effect, can also be of interest in a cosmological setup. In this study, we investigate this effect in the time-dependent chiral asymmetric fermion background and in the expanding spacetime by formulating the effective action of gravitational waves. We also analyze the anomaly equation to see how the backreaction from gravitational waves to thermal chiral plasma occurs. We find that the non-trivial time dependence of chiral chemical potential, which can be induced in some scenarios of baryogenesis, is the key ingredient of the chiral gravitational effect. It turns out that the "memory" of the effect is imprinted on the high frequency gravitational waves propagating in the plasma. Cosmological implications and potential effects on the gravitational wave observation are briefly discussed.

Accurate phase diagrams of multicomponent plasmas are required for the modeling of dense stellar plasmas, such as those found in the cores of white dwarf stars and the crusts of neutron stars. Those phase diagrams have been computed using a variety of standard techniques, which suffer from physical and computational limitations. Here, we present an efficient and accurate method that overcomes the drawbacks of previously used approaches. In particular, finite-size effects are avoided as each phase is calculated separately; the plasma electrons and volume changes are explicitly taken into account; and arbitrary analytic fits to simulation data are avoided. Furthermore, no simulations at uninteresting state conditions, i.e., away from the phase coexistence curves, are required, which improves the efficiency of the technique. The method consists of an adaptation of the so-called Gibbs-Duhem integration approach to electron-ion plasmas, where the coexistence curve is determined by direct numerical integration of its underlying Clapeyron equation. The thermodynamics properties of the coexisting phases are evaluated separately using Monte Carlo simulations in the isobaric semi-grand canonical ensemble. We describe this Monte Carlo-based Clapeyron integration method, including its basic principles, our extension to electron-ion plasmas, and our numerical implementation. We illustrate its applicability and benefits with the calculation of the melting curve of dense C/O plasmas under conditions relevant for white dwarf cores and provide analytic fits to implement this new melting curve in white dwarf models. While this work focuses on the liquid-solid phase boundary of dense two-component plasmas, a wider range of physical systems and phase boundaries are within the scope of the Clapeyron integration method, which had until now only been applied to simple model systems of neutral particles.

We consider the cosmology obtained using scalar fields with a negative potential energy, such as employed to obtain an Ekpyrotic phase of contraction. Applying the covariant entropy bound to the tower of states dictated by the distance conjecture, we find that the relative slope of the potential $|V^{\prime}| / |V|$ is bounded from below by a constant of the order one in Planck units. This is consistent with the requirement to obtain slow Ekpyrotic contraction. We also derive a refined condition on the potential which holds near local minima of a negative potential.