Papers for Tuesday, Jun 29 2021

We present an automated non-parametric light profile extraction pipeline called AutoProf. All steps for extracting surface brightness profiles are included in AutoProf, allowing streamlined analyses of galaxy images. AutoProf improves upon previous non-parametric ellipse fitting implementations with fit-stabilization procedures adapted from machine learning techniques. Additional advanced analysis methods are included in the flexible pipeline for the extraction of alternative brightness profiles (along radial or axial slices), smooth axisymmetric models, and the implementation of decision trees for arbitrarily complex pipelines. Detailed comparisons with widely used photometry algorithms (photutils, XVISTA, and GALFIT) are presented. These comparisons rely on a large collection of late-type galaxy images from the PROBES survey. The direct comparison of surface brightness profiles shows that AutoProf can reliably extract fainter isophotes than other methods on the same images, typically by >2 mag/arcsec^2. Contrasting non-parametric elliptical isophote fitting with simple parametric models also shows that two-component fits (e.g., S\'ersic plus exponential) are insufficient to describe late-type galaxies with high fidelity. It is established that elliptical isophote fitting, and in particular AutoProf, is ideally suited for a broad range of automated isophotal analysis tasks. AutoProf is freely available to the community at: https://github.com/ConnorStoneAstro/AutoProf

We present spatially resolved Hubble Space Telescope grism spectroscopy of 15 galaxies at $z\sim0.8$ drawn from the DEEP2 survey. We analyze H$\alpha$+[N II], [S II] and [S III] emission on kpc scales to explore which mechanisms are powering emission lines at high redshifts, testing which processes may be responsible for the well-known offset of high redshift galaxies from the $z\sim0$ locus in the [O III]/H$\beta$ versus [N II]/H$\alpha$ BPT (Baldwin-Phillips-Terlevich) excitation diagram. We study spatially resolved emission line maps to examine evidence for active galactic nuclei (AGN), shocks, diffuse ionized gas (DIG), or escaping ionizing radiation, all of which may contribute to the BPT offsets observed in our sample. We do not find significant evidence of AGN in our sample and quantify that, on average, AGN would need to contribute $\sim$25% of the H$\alpha$ flux in the central resolution element in order to cause the observed BPT offsets. We find weak ($2\sigma$) evidence of DIG emission at low surface brightnesses, yielding an implied total DIG emission fraction of $\sim$20%, which is not significant enough to be the dominant emission line driver in our sample. In general we find that the observed emission is dominated by star forming H II regions. We discuss trends with demographic properties and the possible role of $\alpha$-enhanced abundance patterns in the emission spectra of high redshift galaxies. Our results indicate that photo-ionization modeling with stellar population synthesis inputs is a valid tool to explore the specific star formation properties which may cause BPT offsets, to be explored in future work.

The level of star formation in elliptical galaxies is poorly constrained, due to difficulties in quantifying the contamination of flux-based estimates of star formation from unrelated phenomena, such as AGN and old stellar populations. We here utilise core-collapse supernovae (CCSNe) as unambiguous tracers of recent star formation in ellipticals within a cosmic volume. We firstly isolate a sample of 421 z < 0.2, r < 21.8 mag CCSNe from the SDSS-II Supernova Survey. We then introduce a Bayesian method of identifying ellipticals via their colours and morphologies in a manner unbiased by redshift and yet consistent with manual classification from Galaxy Zoo 1. We find ~ 25 % of z < 0.2 r < 20 mag galaxies in the Stripe 82 region are ellipticals (~ 28000 galaxies). In total, 36 CCSNe are found to reside in ellipticals. We demonstrate that such early-types contribute a non-negligible fraction of star formation to the present-day cosmic budget, at 11.2 $\pm$ 3.1 (stat) $^{+3.0}_{-4.2}$ (sys) %. Coupling this result with the galaxy stellar mass function of ellipticals, the mean specific star formation rate (SSFR; $\overline{S}$) of these systems is derived. The best-fit slope is given by log ($\overline{S}(M)$/yr) = - (0.80 $\pm$ 0.59) log ($M/10^{10.5}\rm{M}_{\odot}$) - 10.83 $\pm$ 0.18. The mean SSFR for all log ($M/\rm{M}_{\odot}$) > 10.0 ellipticals is found to be $\overline{S} = 9.2 \pm 2.4$ (stat) $^{+2.7}_{-2.3}$ (sys) $\times 10^{-12}$ yr$^{-1}$, which is consistent with recent estimates via SED-fitting, and is 11.8 $\pm$ 3.7 (stat) $^{+3.5}_{-2.9}$ (sys) % of the mean SSFR level on the main sequence as also derived from CCSNe. We find the median optical spectrum of elliptical CCSN hosts is statistically consistent with that of a control sample of ellipticals that do not host CCSNe, implying that these SN-derived results are well-representative of the total low-z elliptical population.

We present the results from a new search for candidate galaxies at z ~ 8.5-11 discovered over the 850 arcmin^2 area probed by the Cosmic Assembly Near-Infrared Deep Extragalactic Legacy Survey (CANDELS). We use a photometric redshift selection including both Hubble and Spitzer Space Telescope photometry to robustly identify galaxies in this epoch at F160W < 26.6. We use a detailed vetting procedure, including screening for persistence, stellar contamination, inclusion of ground-based imaging, and followup space-based imaging to build a robust sample of 11 candidate galaxies, three presented here for the first time. The inclusion of Spitzer/IRAC photometry in the selection process reduces contamination, and yields more robust redshift estimates than Hubble alone. We constrain the evolution of the rest-frame ultraviolet luminosity function via a new method of calculating the observed number densities without choosing a prior magnitude bin size. We find that the abundance at our brightest probed luminosities (M_UV=-22.3) is consistent with predictions from simulations which assume that galaxies in this epoch have gas depletion times at least as short as those in nearby starburst galaxies. Due to large Poisson and cosmic variance uncertainties we cannot conclusively rule out either a smooth evolution of the luminosity function continued from z=4-8, or an accelerate decline at z > 8. We calculate that the presence of seven galaxies in a single field (EGS) is an outlier at the 2-sigma significance level, implying the discovery of a significant overdensity. These scenarios will be imminently testable to high confidence within the first year of observations of the James Webb Space Telescope.

We study and model the properties of galaxy clusters in the normal-branch Dvali-Gabadadze-Porrati (nDGP) model of gravity, which is representative of a wide class of theories which exhibit the Vainshtein screening mechanism. Using the first cosmological simulations which incorporate both full baryonic physics and nDGP, we find that, despite being efficiently screened within clusters, the fifth force can raise the temperature of the intra-cluster gas, affecting the scaling relations between the cluster mass and three observable mass proxies: the gas temperature, the Compton $Y$-parameter of the Sunyaev-Zel'dovich effect and the X-ray analogue of the $Y$-parameter. Therefore, unless properly accounted for, this could lead to biased measurements of the cluster mass in tests that make use of cluster observations, such as cluster number counts, to probe gravity. Using a suite of dark-matter-only simulations, which span a wide range of box sizes and resolutions, and which feature very different strengths of the fifth force, we also calibrate general fitting formulae which can reproduce the nDGP halo concentration at percent accuracy for $0\leq z\leq1$, and halo mass function with $\lesssim3\%$ accuracy at $0\leq z\leq1$ (increasing to $\lesssim5\%$ for $1\leq z\leq 2$), over a halo mass range spanning four orders of magnitude. Our model for the concentration can be used for converting between halo mass overdensities and predicting statistics such as the nonlinear matter power spectrum. The results of this work will form part of a framework for unbiased constraints of gravity using the data from ongoing and upcoming cluster surveys.

We examine possible environmental sources of the enhanced star formation and active galactic nucleus (AGN) activity in the $z = 3.09$ SSA22 protocluster using Hubble WFC3 F160W ($\sim1.6\ \rm \mu m$) observations of the SSA22 field, including new observations centered on eight X-ray selected protocluster AGN. To investigate the role of mergers in the observed AGN and star formation enhancement, we apply both quantitative (S\'ersic-fit and Gini-$M_{20}$) and visual morphological classifications to F160W images of protocluster Lyman break galaxies (LBGs) in the fields of the X-ray AGN and $z \sim 3$ field LBGs in SSA22 and GOODS-N. We find no statistically significant differences between the morphologies and merger fractions of protocluster and field LBGs, though we are limited by small number statistics in the protocluster. We also fit the UV-to-near-IR spectral energy distributions (SED) of F160W-detected protocluster and field LBGs to characterize their stellar masses and star formation histories (SFH). We find that the mean protocluster LBG is by a factor of $\sim2$ times more massive and more attenuated than the mean $z \sim 3$ field LBG. We take our results to suggest that ongoing mergers are not more common among protocluster LBGs than field LBGs, though protocluster LBGs appear to be more massive. We speculate that the larger mass of the protocluster LBGs contributes to the enhancement of SMBH mass and accretion rate in the protocluster, which in turn drives the observed protocluster AGN enhancement.

Context. Characterizing the evolution of protoplanetary disks is necessary to improve our understanding of planet formation. Constraints on both dust and gas are needed to determine the dominant disk dissipation mechanisms. Aims. We aim to compare the disk dust masses in the Chamaeleon II (Cha II) star-forming region with other regions with ages between 1 and 10Myr. Methods. We use ALMA band 6 observations (1.3 mm) to survey 29 protoplanetary disks in Cha II. Dust mass estimates are derived from the continuum data. Results. Out of our initial sample of 29 disks, we detect 22 sources in the continuum, 10 in 12CO, 3 in 13CO, and none in C18O (J=2-1). Additionally, we detect two companion candidates in the continuum and 12CO emission. Most disk dust masses are lower than 10Mearth, assuming thermal emission from optically thin dust. We compare consistent estimations of the distributions of the disk dust mass and the disk-to-stellar mass ratios in Cha II with six other low mass and isolated star-forming regions in the age range of 1-10Myr: Upper Sco, CrA, IC 348, Cha I, Lupus, and Taurus. When comparing the dust-to-stellar mass ratio, we find that the masses of disks in Cha II are statistically different from those in Upper Sco and Taurus, and we confirm that disks in Upper Sco, the oldest region of the sample, are statistically less massive than in all other regions. Performing a second statistical test of the dust mass distributions from similar mass bins, we find no statistical differences between these regions and Cha II. Conclusions. We interpret these trends, most simply, as a sign of decline in the disk dust masses with time or dust evolution. Different global initial conditions in star-forming regions may also play a role, but their impact on the properties of a disk population is difficult to isolate in star-forming regions lacking nearby massive stars.

Intermediate-mass black holes (IMBHs) by definition have masses of $M_{\rm IMBH} \sim 10^{2-5}~M_\odot$, a range with few observational constraints. Finding IMBHs in globular star clusters (GCs) would validate a formation channel for massive black-hole seeds in the early universe. Here, we simulate a 60-hour observation with the next-generation Very Large Array (ngVLA) of 728 GC candidates in the Virgo Cluster galaxy NGC\,4472. Interpreting the radio detection thresholds as signatures of accretion onto IMBHs, we benchmark IMBH mass thresholds in three scenarios and find the following: (1) Radio analogs of ESO\,243-49 HLX-1, a strong IMBH candidate with $M_{\rm IMBH}^{\rm HLX} \sim 10^{4-5}~M_\odot$ in a star cluster, are easy to access in all 728 GC candidates. (2) For the 30 GC candidates with extant X-ray detections, the empirical fundamental-plane relation involving black hole mass plus X-ray and radio luminosities suggests access to $M_{\rm IMBH}^{\rm FP} \sim 10^{1.7-3.6}~M_\odot$, with an uncertainty of 0.44 dex. (3) A fiducial Bondi accretion model was applied to all 728 GC candidates and to radio stacks of GC candidates. This model suggests access to IMBH masses, with a statistical uncertainty of 0.39 dex, of $M_{\rm IMBH}^{\rm B} \sim 10^{4.9-5.1}~M_\odot$ for individual GC candidates and $M_{\rm IMBH}^{\rm B,stack} \sim 10^{4.5}~M_\odot$ for radio stacks of about 100-200 GC candidates. The fiducial Bondi model offers initial guidance, but is subject to additional systematic uncertainties and should be superseded by hydrodynamical simulations of gas flows in GCs.

A high sensitivity, 7mm Very Long Baseline Array image of M\,87 is analyzed in order to estimate the jet velocity within 0.65 mas of the point of origin. The image captured a high signal to noise, double-ridged, counter-jet extending $\sim 1$ mas from the nucleus. After defining conditions and requirements that justify approximate time averaged bilateral symmetry, a continuous set of Lorentz transformations are found that map the double-ridged counter-jet intensity profile into the double-ridged jet intensity profile. The mapping is realized by a uniformly accelerating flow with intrinsic velocity of $\sim 0.27$c at 0.4 mas (a de-projected distance of 0.38 lt-yrs) to $0.38$c at 0.65 mas (a de-projected distance of 0.61 lt-yrs) from the nucleus. Since the velocity field is derived from the global surface brightness profile and does not depend on the motion of enhanced features, it is most likely a bulk flow velocity as opposed to a pattern velocity. This interpretation is corroborated by the fact that the distribution of the apparent velocities of previously identified individual features (from the literature) within 0.65 mas of the nucleus are consistent with local hydrodynamic shocks being advected with the local bulk flow velocity. The bulk flow velocity of the visible inner jet is a constraint that can potentially break degeneracies between numerical simulations that are designed to replicate both the annulus that was imaged by the Event Horizon Telescope as well as the base of the inner jet.

Recent near-IR Surveys have discovered a number of new bulge globular cluster (GC) candidates that need to be further investigated. Our main objective is to use public data from the Gaia Mission, VVV, 2MASS and WISE in order to measure the physical parameters of Minni48, a new candidate GC located in the inner bulge of the Galaxy at l=359.35 deg, b=2.79 deg. Even though there is a bright foreground star contaminating the field, the cluster appears quite bright in near- and mid-IR images. We obtain deep decontaminated optical and near-IR colour-magnitude diagrams (CMDs) for this cluster. The heliocentric cluster distance is determined from the red clump (RC) and the red giant branch (RGB) tip magnitudes in the near-IR CMD, while the cluster metallicity is estimated from the RGB slope and the fit to theoretical isochrones. The GC size is found to be r = 6' +/- 1', while reddening and extinction values are E(J-Ks)=0.60 +/- 0.05 mag, A_G=3.23 +/- 0.10 mag, A_Ks=0.45 +/- 0.05 mag. The resulting mean Gaia proper motions are PMRA=-3.5 +/- 0.5 mas/yr, PMDEC=-6.0 +/- 0.5 mas/yr. The IR magnitude of the RC yields an accurate distance modulus estimate of (m-M)_0=14.61 mag, equivalent to a distance D=8.4 +/- 1.0 kpc. This is consistent with the optical distance estimate: (m-M)_0=14.67 mag, D=8.6 +/- 1.0 kpc, and with the RGB tip distance: (m-M)_0=14.45 mag, D=7.8 +/- 1.0 kpc. The derived metallicity is [Fe/H]=-0.20 +/- 0.30 dex. A good fit to the PARSEC stellar isochrones is obtained in all CMDs using Age = 10 +/- 2 Gyr. The total absolute magnitude of this GC is estimated to be M_Ks= -9.04 +/- 0.66 mag. Based on its position, kinematics, metallicity and age, we conclude that Minni48 is a genuine GC, similar to other well known metal-rich bulge GCs. It is located at a projected Galactocentric angular distance of 2.9 deg, equivalent to 0.4 kpc, being one of the closest GCs to the Galactic centre.

I report new orbital periods (P) for 13 classical novae, based on light curves from TESS, AAVSO, and other public archives. These new nova periods now constitute nearly one-seventh of all known nova periods. Five of my systems have P>1 day, which doubles the number of such systems that must have evolved companion stars. (This is simply because ground-based time series have neither the coverage nor the stability required to discover these small-amplitude long periods.) V1016 Sgr has a rare P below the period gap, and suddenly becomes useful for current debates on evolution of novae. Five of the novae (FM Cir, V399 Del, V407 Lup, YZ Ret, and V549 Vel) have the orbital modulations in the tail of the eruption after the transition phase. Soon after the transition, YZ Ret shows a unique set of aperiodic diminishing oscillations, plus YZ Ret shows two highly-significant transient periods, 1.1% and 4.5% longer than the orbital period, much like for the superhump phenomenon. I also report an optical 591.27465 second periodicity for V407 Lup, which is coherent and must be tied to the white dwarf spin period. The new orbital periods in days are 0.1883907 +- 0.0000048 for V1405 Cas, 3.4898 +- 0.0072 for FM Cir, 0.162941 +- 0.000060 for V339 Del, 3.513 +- 0.020 for V407 Lup, 1.32379 +- 0.00048 for V2109 Oph, 3.21997 +- 0.00039 for V392 Per, 0.1628714 +- 0.0000110 for V598 Pup, 0.1324539 +- 0.0000098 for YZ Ret, 0.07579635 +- 0.00000017 for V1016 Sgr, 7.101 +- 0.016 for V5583 Sgr, 0.61075 +- 0.00071 for V1534 Sco, 0.40319 +- 0.00005 for V549 Vel, and 0.146501 +- 0.000058 for NQ Vul.

The narrow-line region (NLR) consists of gas clouds ionized by the strong radiation from the active galactic nucleus (AGN), distributed in the spatial scale of AGN host galaxies. The strong emission lines from the NLR are useful to diagnose physical and chemical properties of the interstellar medium in AGN host galaxies. However, the origin of the NLR is unclear; the gas clouds in NLRs may be originally in the host and photoionized by the AGN radiation, or they may be transferred from the nucleus with AGN-driven outflows. For studying the origin of the NLR, we systematically investigate the gas density and velocity dispersion of NLR gas clouds using a large spectroscopic data set taken from the Sloan Digital Sky Survey. The [S II] emission-line flux ratio and [O III] velocity width of 9,571 type-2 Seyfert galaxies and 110,041 star-forming galaxies suggest that the gas density and velocity dispersion of NLR clouds in Seyfert galaxies (ne ~ 194 cm-3 and sigma([O III]) ~ 147 km s-1) are systematically larger than those of clouds in H II regions of star-forming galaxies (ne ~ 29 cm-3 and sigma([O III]) ~ 58 km s-1). Interestingly, the electron density and velocity dispersion of NLR gas clouds are larger for Seyfert galaxies with a higher [O III]/Hbeta flux ratio, i.e., with a more active AGN. We also investigate the spatially-resolved kinematics of ionized gas clouds using the Mapping Nearby Galaxies at the Apache Point Observatory (MaNGA) survey data for 90 Seyfert galaxies and 801 star-forming galaxies. We find that the velocity dispersion of NLR gas clouds in Seyfert galaxies is larger than that in star-forming galaxies at a fixed stellar mass, at both central and off-central regions. These results suggest that gas clouds in NLRs come from the nucleus, probably through AGN outflows.

The success of primordial nucleosynthesis has been limited by the long-standing Lithium problem. We use a self-consistent perturbative analysis of the effects of relevant theoretical parameters on primordial nucleosynthesis, including variations of nature's fundamental constants, to explore the problem and its possible solutions, in the context of the latest observations and theoretical modeling. We quantify the amount of depletion needed to solve the Lithium problem, and show that transport processes of chemical elements in stars are able to account for it. Specifically, the combination of atomic diffusion, rotation and penetrative convection allows us to reproduce the lithium surface abundances of Population II stars, starting from the primordial Lithium abundance. We also show that even with this depletion factor there is a preference for a value of the fine-structure constant at this epoch that is larger than the current laboratory one by a few parts per million of relative variation, at the two to three standard deviations level of statistical significance. This preference is driven by the recently noticed discrepancy between the best-fit values for the baryon-to-photon ratio (or equivalently the Deuterium abundance) inferred from cosmic microwave background and primordial nucleosynthesis analyses, and is largely insensitive to the Helium-4 abundance. We thus conclude that the Lithium problem most likely has an astrophysical solution, while the Deuterium discrepancy provides a possible hint of new physics.

We investigate the inflationary consequences of the oscillating dark energy model proposed by Ti\'an [\href{https://doi.org/10.1103/PhysRevD.101.063531}{Phys. Rev. D {\bf 101}, 063531 (2020)}], which aims to solve the cosmological coincidence problem with multi-accelerating Universe (MAU). We point out that the inflationary dynamics belong to slow-roll inflation. The spectral index of scalar perturbations and the tensor-to-scalar ratio $r$ are shown to be consistent with current \textit{Planck} measurements. Especially, this model predicts $r\sim10^{-7}$, which is far below the observation limits. This result motivates us to explore the smallness of $r$ in the general MAU. We propose a quintessential generalization of the original model and prove $r<0.01$ in general. The null detection to date of primordial gravitational waves provides a circumstantial evidence for the MAU. After the end of inflation, the scalar field rolls toward infinity instead of a local minimum, and meanwhile its equation of state is oscillating with an average value larger than $1/3$. In this framework, we show that gravitational particle creation at the end of inflation is capable of reheating the Universe.

We have observed the young protostellar system NGC 2264 CMM3 in the 1.3 mm and 2.0 mm bands at a resolution of about 0.1$"$ (70 au) with ALMA. The structures of two distinct components, CMM3A and CMM3B, are resolved in the continuum images of both bands. CMM3A has an elliptical structure extending along the direction almost perpendicular to the known outflow, while CMM3B reveals a round shape. We have fitted two 2D-Gaussian components to the elliptical structure of CMM3A and CMM3B, and have separated the disk and envelope components for each source. The spectral index $\alpha$ between 2.0 mm and 0.8 mm is derived to be 2.4-2.7 and 2.4-2.6 for CMM3A and CMM3B, respectively, indicating the optically thick dust emission and/or the grain growth. A velocity gradient in the disk/envelope direction is detected for CMM3A in the CH$_3$CN, CH$_3$OH, and $^{13}$CH$_3$OH lines detected in the 1.3 mm band, which can be interpreted as the rotation of the disk/envelope system. From this result, the protostellar mass of CMM3A is roughly evaluated to be $0.1- 0.5$ $M_\odot$ by assuming Keplerian rotation. The mass accretion rate is thus estimated to be $5\times10^{-5}$ - 4 $\times$ $10^{-3}$ $M_\odot$ yr$^{-1}$, which is higher than typical mass accretion rate of low-mass protostars. The OCS emission line shows a velocity gradient in both outflow direction and disk/envelope direction. A hint of outflow rotation is found, and the specific angular momentum of the outflow is estimated to be comparable to that of the disk. These results provide us with novel information on the initial stage of a binary/multiple system.

In this study, we investigate an extreme ultraviolet (EUV) wave event on 2010 February 11, which occurred as a limb event from the Earth viewpoint and a disk event from the STEREO--B viewpoint. We use the data obtained by the Atmospheric Imaging Assembly (AIA) aboard the Solar Dynamics Observatory (SDO) in various EUV channels. The EUV wave event was launched by a partial prominence eruption. Similar to some EUV wave events in previous works, this EUV wave event contains a faster wave with a speed of $\sim$445$\pm$6 km s$^{-1}$, which we call coronal Moreton wave, and a slower wave with a speed of $\sim$298$\pm$5 km s$^{-1}$, which we call "EIT wave". The coronal Moreton wave is identified as a fast-mode wave and the "EIT wave" is identified as an apparent propagation due to successive field-line stretching. We also observe a stationary front associated with the fast mode EUV wave. This stationary front is explained as mode conversion from the coronal Moreton wave to a slow-mode wave near a streamer.

We describe the deployment and first tests on Sky of CONCERTO, a large field-of-view (18.6arc-min) spectral-imaging instrument. The instrument operates in the range 130-310GHz from the APEX 12-meters telescope located at 5100m a.s.l. on the Chajnantor plateau. Spectra with R=1-300 are obtained using a fast (2.5Hz mechanical frequency) Fourier Transform Spectrometer (FTS), coupled to a continuous dilution cryostat with a base temperature of 60mK. Two 2152-pixels arrays of Lumped Element Kinetic Inductance Detectors (LEKID) are installed in the cryostat that also contains the cold optics and the front-end electronics. CONCERTO, installed in April 2021, generates more than 20k spectra per second during observations. We describe the final development phases, the installation and the first results obtained on Sky.

The core accretion scenario of planet formation assumes that planetesimals and planetary embryos are formed during the primordial, gaseous phases of the protoplanetary disk. However, how the dust particles overcome the traditional growth barriers is not well understood. The recently proposed viscous ring-instability may explain the concentric rings observed in protoplanetary disks by assuming that the dust grains can reduce the gas conductivity, which can weaken the magneto-rotational instability. We present an analysis of this model with the help of GPU-based numerical hydrodynamic simulations of coupled gas and dust in the thin-disk limit. During the evolution of the disk the dusty rings become Rossby unstable and break up into a cascade of small-scale vortices. The vortices form secularly stable dusty structures, which could be sites of planetesimal formation by the streaming instability as well as direct gravitational collapse. The phenomenon of self-sustaining vortices is consistent with observational constraints of exoplanets and sets a favorable environment for planetary system formation.

We present a comprehensive radiative magnetohydrodynamic simulation of the quiet Sun and large solar active regions. The 197 Mm wide simulation domain spans from the uppermost convection zone to over 100 Mm in the solar corona. Sophisticated treatments of radiative transfer and conduction transport provide the necessary realism for synthesizing observables to compare with remote sensing observations of the Sun. This model self-consistently reproduces observed features of the quiet Sun, emerging and developed active regions, and solar flares up to M class. Here, we report an overview on the first results. The surface magnetoconvection yields an upward Poynting flux that is dissipated in the corona and heats the plasma to over one million K. The quiescent corona also presents ubiquitous propagating waves, jets, and bright points with sizes down to 2 Mm. Magnetic flux bundles generated in a solar convective dynamo emerge into the photosphere and gives rise to strong and complex active regions with Over $10^{23}$ Mx magnetic flux. The coronal free magnetic energy, which is about 18\% of the total magnetic energy, accumulates to about $10^{33}$ erg. The coronal magnetic field is not forcefree, as the Lorentz force needs to balance the pressure force and viscous stress as well as to drive magnetic field evolution. Emission measure from $\log_{10}T = 4.5$ to $\log_{10}T > 7$ provides a comprehensive view on structures and dynamics in the active region corona, such as coronal loops in various lengths and temperatures, mass circulation by evaporation and condensation, and eruptions from jets to large-scale mass ejections.

Gamma-ray bursts (GRBs) are the most luminous explosions and can be detectable out to the edge of Universe. It has long been thought they can extend the Hubble diagram to very high redshifts. Several correlations between temporal or spectral properties and GRB luminosities have been proposed to make GRBs cosmological tools. However, those correlations cannot be properly standardized. In this paper, we select a long GRB sample with X-ray plateau phases produced by electromagnetic dipole emissions from central new-born magnetars. A tight correlation is found between the plateau luminosity and the end time of the plateau in X-ray afterglows out to the redshift $z=5.91$. We standardize these long GRBs X-ray light curves to a universal behavior by this correlation for the first time, with a luminosity dispersion of 0.5 dex. The derived distance-redshift relation of GRBs is in agreement with the standard $\Lambda$CDM model both at low and high redshifts. The evidence of accelerating universe from this GRB sample is $3\sigma$, which is the highest statistical significance from GRBs to date.

In a recent study, we developed a method to model the impact of photometric redshift uncertainty on the two-point correlation function (2PCF). In this method, we can obtain both the intrinsic clustering strength and the photometric redshift errors simultaneously by fitting the projected 2PCF with two integration depths along the line-of-sight. Here we apply this method to the DESI Legacy Imaging Surveys Data Release 8 (LS DR8), the largest galaxy sample currently available. We separate galaxies into 20 samples in 8 redshift bins from $z=0.1$ to $z=1.0$, and a few $\rm z$-band absolute magnitude bins, with $M_{\rm z} \le -20$. These galaxies are further separated into red and blue sub-samples according to their $M^{0.5}_{\rm r}-M^{0.5}_{\rm z}$ colors. We measure the projected 2PCFs for all these galaxy (sub-)samples, and fit them using our photometric redshift 2PCF model. We find that the photometric redshift errors are smaller in red sub-samples than the overall population. On the other hand, there might be some systematic photometric redshift errors in the blue sub-samples, so that some of the sub-samples show significantly enhanced 2PCF at large scales. Therefore, focusing only on the red and all (sub-)samples, we find that the biases of galaxies in these (sub-)samples show clear color, redshift and luminosity dependencies, in that red brighter galaxies at higher redshift are more biased than their bluer and low redshift counterparts. Apart from the best fit set of parameters, $\sigma_{z}$ and $b$, from this state-of-the-art photometric redshift survey, we obtain high precision intrinsic clustering measurements for these 40 red and all galaxy (sub-)samples. These measurements on large and small scales hold important information regarding the cosmology and galaxy formation, which will be used in our subsequent probes in this series.

During April and May 2020, SGR J1935+2154 emitted hundreds of short bursts and became one of the most prolific transient magnetars. At the onset of the active bursting period, a 130-s burst ``forest", which included some bursts with peculiar time profiles, were observed with the $Fermi$/Gamma-ray Burst Monitor. In this paper, we present the results of time-resolved spectral analysis of this burst ``forest" episode, which occurred on April 27, 2020. We identify thermal spectral components prevalent during the entire 130-s episode; high-energy maxima appear during the photon flux peaks, which are modulated by the spin period of the source. Moreover, the evolution of the $\nu F_{\nu}$ spectral hardness (represented by $E_{\rm peak}$ or blackbody temperature) within the lightcurve peaks is anti-correlated with the pulse phases extrapolated from the pulsation observed within the persistent soft X-ray emission of the source six hours later. Throughout the episode, the emitting area of the high-energy (hotter) component is 1--2 orders of magnitude smaller than that for the low-energy component. We interpret this with a geometrical viewing angle scenario, inferring that the high-energy component likely originates from a low-altitude hotspot located within closed toroidal magnetic field lines.

For the next generation of neutrino telescopes at the South Pole new optical sensors with a segmented photosensitive area are being developed. These sensors will significantly increase the detectors' sensitivity not only to high-energy astrophysical neutrinos, but also to neutrinos in the MeV energy range, such as those produced in supernovae during core collapse. These low-energy neutrinos can provide a detailed picture of the events that follow the collapse of the stellar core, thus verifying and improving our understanding of these massive explosions. The new sensor design has the potential to enable event-based detection of MeV neutrinos with a single sensor while effectively suppressing background. This paper presents the results of studies on the sensitivity of such a segmented sensor to MeV neutrinos and, for the first time, the potential of a corresponding detector in the deep ice at the South Pole for the detection of extra-galactic core-collapse supernovae (CCSN). We find that using coincidence conditions between the photocathode segments within a sensor we can detect a CCSN with a progenitor mass of $27\ \mathrm{M}_{\odot}$ up to a distance of $370\,\mathrm{kpc}$ with a false detection rate of $0.4$ per year. If the arrival burst time is known from an independent observation with $\delta t = 1\,\mathrm{h}$, such CCSN can be detected with a probability of $50\%$ at [407, 341]$\,$kpc distance with a certainty of [3.2, 4.9]$\,\sigma$ that the signal was not produced by background fluctuations.

We present results of phase- and time-resolved study of iron spectral features in the emission of the Be/X-ray transient pulsar V0332+53, during its type II outburst in 2004 using archival RXTE/PCA data. Coherent pulsations of both fluorescent iron line at {\approx}6.4 keV and neutral iron K-edge at {\approx}7.1 keV have been detected throughout the entire outburst. The pulsating iron K-edge is reported for the first time for this object. Near the peak of the outburst, the 3-12 keV pulse profile shows two deep, F_max / F_min ~ 2, and narrow dips of nearly identical shape, separated by exactly {\Delta}{\phi}=0.5 in phase. The dip spectra are nearly identical to each other and very similar in shape to the spectra outside the dips. The iron K-edge peaks at the phase intervals corresponding to the dips, although its optical depth {\tau}_K ~ 0.05 is by far insufficient to explain the dips. The iron line shows pulsations with a complex pulse profile without any obvious correlation with the total flux or optical depth of the K-edge. Accounting for the component associated with reprocessing of the pulsar emission by the surface of the donor star and circumstellar material, we find a very high pulsation amplitude of the iron line flux, F_max / F_min ~ 10. We demonstrate that these properties of V0332+53 can not be explained by contemporary emission models for accreting X-ray pulsars and speculate about the origin of the observed iron spectral features.

We report an ultra-bright lensed submillimeter galaxy at $z_{spec}=1.4796$, identified as a result of a full-sky cross-correlation of the AllWISE and Planck compact source catalogs aimed to search for bright submillimeter galaxies at $z \sim 1.5-2.8$. APEX/LABOCA observations of the candidate galaxy reveal a source with flux (S$_{870 \mu m}= 54\pm 8$ mJy). The position of the APEX source coincides with the position of the AllWISE mid-IR source, and with the Einstein ring GAL-CLUS-022058s, observed with the HST. Archival VLT/FORS observations reveal the redshift of this Einstein ring, $z_{spec}=1.4796$, and detection of the CO(5-4) line at $z_{spec} = 1.4802$ with APEX/nFLASH230 confirms the redshift of the submillimeter emission. The lensed source appears to be gravitationally magnified by a massive foreground galaxy cluster lens at $z = 0.36$. We use Lenstool to model the gravitational lensing, which is near to a "fold arc" configuration for an elliptical mass distribution of the central halo, where four images of the lensed galaxy are seen; the mean magnification is $\mu_{\rm L} =18\pm 4$. We have determined an intrinsic rest-frame infrared luminosity of $L_{IR} \approx 10^{12} L_\odot $ and a likely star formation rate of $\sim 70-170$ $M_\odot\ yr^{-1}$. The molecular gas mass is $M_{mol} \sim 2.6 \times 10^{10} M_\odot$ and the gas fraction is $f = 0.34\pm 0.07$. We also obtain a stellar mass log$(M_\ast/M_\odot) = 10.7 \pm 0.1$ and a specific star formation rate log$(sSFR/Gyr^{-1})=0.15 \pm 0.03$. This galaxy lies on the so-called main sequence of star-forming galaxies at this redshift.

We use the data-driven method, The Cannon, to bring 21,000 stars from the ARGOS bulge survey, including 10,000 red clump stars, onto the parameter and abundance scales of the cross-Galactic survey, APOGEE, obtaining rms precisions of 0.10 dex, 0.07 dex, 74 K, and 0.18 dex for [Fe/H], [Mg/Fe], Teff, and log(g), respectively. The re-calibrated ARGOS survey - which we refer to as the A2A survey - is combined with the APOGEE survey to investigate the abundance structure of the Galactic bulge. We find X-shaped [Fe/H] and [Mg/Fe] distributions in the bulge that are more pinched than the bulge density, a signature of its disk origin. The mean abundance along the major axis of the bar varies such that the stars are more [Fe/H]-poor and [Mg/Fe]-rich near the Galactic center than in the long bar/outer bulge region. The vertical [Fe/H] and [Mg/Fe] gradients vary between the inner bulge and long bar with the inner bulge showing a flattening near the plane that is absent in the long bar. The [Fe/H]-[Mg/Fe] distribution shows two main maxima, an ``[Fe/H]-poor [Mg/Fe]- rich'' maximum and an ``[Fe/H]-rich [Mg/Fe]-poor'' maximum, that vary in strength with position in the bulge. In particular, the outer long bar close to the Galactic plane is dominated by super-solar [Fe/H], [Mg/Fe]-normal stars. Stars composing the [Fe/H]-rich maximum show little kinematic dependence on [Fe/H], but for lower [Fe/H] the rotation and dispersion of the bulge increase slowly. Stars with [Fe/H]<-1 dex have a very different kinematic structure than stars with higher [Fe/H]. Comparing with recent models for the Galactic boxy-peanut bulge, the abundance gradients and distribution, and the relation between [Fe/H] and kinematics suggest that the stars comprising each maximum have separate disk origins with the ``[Fe/H]-poor [Mg/Fe]-rich'' stars originating from a thicker disk than the ``[Fe/H]-rich [Mg/Fe]-poor'' stars.

Homogeneity and isotropy are geometrical properties of the metric of space-time, independent of the amplitude of perturbations, and can be characterized using the fractal (or Hausdorff) dimension. We present such measurement in the Early Universe using the Planck CMB temperature variations in the sky. This new measurement of the isotropy scale is model independent and purely geometrical. We find strong evidence of homogeneity for scales larger than $\pi/3 \simeq 60$ degrees. Surprisingly, this finding is at odds with the $\Lambda$CDM prediction, which assumes a scale invariant (i.e. fractal) universe to larger scales. We estimate the significance of our finding using a principal component analysis from the sampling variations in the sky. The Planck maps are in good agreement with the null hypothesis of isotropy at scales larger than $\pi/3$. While the best fit $\Lambda CDM$ prediction has a $\chi^2$ which is 20 times larger. Assuming translation invariance (and flat space $k=0$) we can convert the isotropy scale $\pi/3$ into an homogeneity scale of $3.3 c/H_0$, which is close to the distance to the last scattering surface.

A solar sail propelled small satellite mission concept to intercept and potentially rendezvous with newly discovered Interstellar Objects (ISOs) is described. The mission concept derives from the proposal for a technology demonstration mission (TDM) for exiting the solar system at high velocity, eventually to reach the focal region of the solar gravitational lens (SGL). The ISO mission concept is to fly a solar sail inward toward a holding orbit around the Sun and when the ISO orbit is confirmed, target the sailcraft to reach an escape velocity of over 6\, AU/year. This would permit rapid response to a new ISO discovery and an intercept within 10 AU from the Sun. Two new proven interplanetary technologies are utilized to enable such a mission: i) interplanetary Smallsats, such as those demonstrated by the MarCO mission, and ii) solar sails, such as demonstrated by LightSail and IKAROS missions and developed for NEA Scout and Solar Cruiser missions. Current technology work suggests such a mission could fly and reach an ISO moving through the solar system within this decade. Such a mission might enable the first encounter with an ISO to allow for imaging and spectroscopy, measurements of size and mass, potentially giving a unique information about the object's origin and composition.

We present the results of a temporal and spectral analysis of the transient source ULX-4 in the galaxy M51. The data used were drawn from {\it Chandra}, {\it XMM-Newton} and {\it XRT} archives, spanning the years 2000-2019. The X-ray flux of the source is seen to vary by two orders of magnitudes within a month but a short-term variability was not observed. We find some evidence for the existence of a bi-modality feature in the flux distribution of ULX-4. The X-ray spectra from most of the data can be fitted by a simple {\it power-law} with a photon index in the range 1.44-2.0 suggesting that the source is in a hard state characterized by non-thermal emission. The presence of a soft component, in two {\it XMM-Newton} observations, is best described by the addition of a mekal component suggesting a contribution from a diffuse background source. We identified two optical sources as possible counterparts within an error radius of 0.18 arcsec at 95% confidence level for ULX-4 based on the archival {\it HST}/ACS and {\it HST}/WFC3 data. The absolute magnitude of these counterparts is -6.6 mag. Blackbody fits of the spectral energy distributions indicate the spectral type to be B-type stars. One of these counterparts exhibits a low-amplitude optical periodicity of 264 $\pm$ 37 days in the F606W filter; If we assume this apparent periodicity is associated with the orbital motion of the donor, then it is more likely that the donor is a red supergiant satisfying the long periodicity and accretion via Roche-lobe overflow. Consequently, the SED would then have to be interpreted as a superposition of emissions from a cold donor and a hot flow component, most likely from an accretion disk. If, on the other hand, the periodicity is super orbital in nature i.e., due to possible interactions of the compact object with a circumstellar disk, the donor could then be a Be/X star hosting a neutron star.

We examine the equilibrium of a magnetized neutron-star-crust. We calculate axially symmetric models in which an elastic force balances solenoidal motion driven by a Lorentz force. A large variety of equilibrium models are allowed by incorporating the elastic shear deformation; in addition, toroidal-magnetic-field dominated models are available. These results remarkably differ from those in barotropic fluid stars. We demonstrate some models wherein the magnetic energy exceeds the elastic energy. The excess comes from the fact that a large amount of magnetic energy is associated with the irrotational part of the magnetic force, which is balanced with gravity and pressure. It is sufficient for equilibrium models that the minor solenoidal part is balanced by a weak elastic force. We find that the elasticity in the crust plays an important role on the magnetic-field confinement. Further, we present the spatial distribution of the shear-stress at the elastic limit, by which the crust-fracture location can be identified. The result has useful implications for realistic crust-quake models.

Jets, collimated outflows of particles and fields, are observed in a wide variety of astrophysical systems, including Active Galactic Nuclei of various types, microquasars, gamma-ray bursts, and young stellar objects. Despite intensive efforts along several decades of observational and theoretical research, there are still many uncertainties and open questions related to how jets are produced and what is their composition. In this review I offer an outline of some current views on the content and basic properties of astrophysical jets.

Wide field small aperture telescopes (WFSATs) are mainly used to obtain scientific information of point--like and streak--like celestial objects. However, qualities of images obtained by WFSATs are seriously affected by the background noise and variable point spread functions. Developing high speed and high efficiency data processing method is of great importance for further scientific research. In recent years, deep neural networks have been proposed for detection and classification of celestial objects and have shown better performance than classical methods. In this paper, we further extend abilities of the deep neural network based astronomical target detection framework to make it suitable for photometry and astrometry. We add new branches into the deep neural network to obtain types, magnitudes and positions of different celestial objects at the same time. Tested with simulated data, we find that our neural network has better performance in photometry than classical methods. Because photometry and astrometry are regression algorithms, which would obtain high accuracy measurements instead of rough classification results, the accuracy of photometry and astrometry results would be affected by different observation conditions. To solve this problem, we further propose to use reference stars to train our deep neural network with transfer learning strategy when observation conditions change. The photometry framework proposed in this paper could be used as an end--to--end quick data processing framework for WFSATs, which can further increase response speed and scientific outputs of WFSATs.

We perform a deep wide-field imaging survey of nearby galaxies using H$\alpha$ and broadband filters to investigate the characteristics of star formation in galaxies. Motivated by the finding that star formation rates (SFRs) derived from H$\alpha$ fluxes in dwarf galaxies are systematically lower than those inferred from far-ultraviolet (FUV) fluxes, we attempt to determine whether the same trend exists in the extended disks of two star-forming galaxies. We perform spatially resolved photometry using grid-shaped apertures to measure the FUV and H$\alpha$ fluxes of star-forming regions. We also perform spectral energy distribution (SED) fittings using 11 photometric data (FUV-to-MIR) including data from the literature to estimate the local properties such as internal attenuation of individual star-forming clumps. Comparing SFR$_\mathrm{FUV}$ and SFR$_\mathrm{H\alpha}$, which are converted from the H$\alpha$ and FUV fluxes corrected for the local properties, we find that SFR$_\mathrm{H\alpha}$/SFR$_\mathrm{FUV}$ tends to decrease as the SFR decreases. We evaluate possible causes of this discrepancy between the two SFRs by restricting parameters in the SED fitting and conclude that deficient H$\alpha$ fluxes in the extended disks of galaxies are tightly correlated with recent starbursts. The strong and short starburst which is being rapidly suppressed over the last 10 Myr seems to induce a significant discrepancy between the SFR$_\mathrm{H\alpha}$ and SFR$_\mathrm{FUV}$. In addition, the recent bursts in the extended disks of galaxies appear to have occurred azimuth-symmetrically, implying that these were likely triggered by gas accretion or internal processes rather than external perturbation.

FU Orionis-type objects (FUors) are low-mass pre-main sequence stars undergoing a temporary, but significant increase of mass accretion rate from the circumstellar disk onto the protostar. It is not yet clear what triggers the accretion bursts and whether the disks of FUors are in any way different from disks of non-bursting young stellar objects. Motivated by this, we conducted a 1.3 mm continuum survey of ten FUors and FUor-like objects with ALMA, using both the 7 m array and the 12 m array in two different configurations to recover emission at the widest possible range of spatial scales. We detected all targeted sources and several nearby objects as well. To constrain the disk structure, we fit the data with models of increasing complexity from 2D Gaussian to radiative transfer, enabling comparison with other samples modeled in a similar way. The radiative transfer modeling gives disk masses that are significantly larger than what is obtained from the measured millimeter fluxes assuming optically thin emission, suggesting that the FUor disks are optically thick at this wavelength. In comparison with samples of regular Class II and Class I objects, the disks of FUors are typically a factor of 2.9-4.4 more massive and a factor of 1.5-4.7 smaller in size. A significant fraction of them (65-70%) may be gravitationally unstable.

The Galactic black hole candidate (BHC) 4U~1630-472 has gone through several outbursts (13 to be particular) in the last two and a half decades starting from the RXTE era till date. Like the outbursts of other transient BHCs, the outbursts of this source show variations in duration, peak numbers, highest peak flux, etc. However, unlike any other soft X-ray transients, this source showed outbursts of two types, such as normal and super. The normal outbursts of duration $\sim 100-200$~days are observed quasi periodically at an average recurrence/quiescence period of $\sim 500$~days. The super outbursts of duration $\sim 1.5-2.5$~years contain one or more normal outbursts other than one mega outburst. We make an effort to separate flux contribution of the normal and the mega outbursts from the super outbursts, and tried to understand the nature of evolution of both types (normal and mega) of outbursts, based on the quiescent period prior to the outbursts. Archival data of RXTE/ASM from January 1996 to June 2011, and MAXI/GSC from August 2009 to July 2020 are used for our study. A possible linear relation between the quiescent and outburstsing periods for both types of outbursts are observed. This makes the BHC a special source, and it may contain two companion binaries. Two companions might be responsible for two types of outbursts.

In agreement with the gravitational-wave events which are constantly increasing, new aspects of the internal structure of compact stars have come to light. A scenario in which a first order transition takes place inside these stars is of particular interest as it can lead, under conditions, to a third gravitationally stable branch (besides white dwarfs and neutron stars). This is known as the twin star scenario. The new branch yields stars with the same mass as normal compact stars but quite different radii. In the current work, we focus on hybrid stars undergone a hadron to quark phase transition near their core and how this new stable configuration arises. Emphasis is to be given especially in the aspects of the phase transition and its parameterization in two different ways, namely with Maxwell construction and with Gibbs construction. Qualitative findings of mass-radius relations of these stars will also be presented.

Exoplanets transiting bright nearby stars are key objects for advancing our knowledge of planetary formation and evolution. The wealth of photons from the host star gives detailed access to the atmospheric, interior, and orbital properties of the planetary companions. $\nu^2$ Lupi (HD 136352) is a naked-eye ($V = 5.78$) Sun-like star that was discovered to host three low-mass planets with orbital periods of 11.6, 27.6, and 107.6 days via radial velocity monitoring (Udry et al. 2019). The two inner planets (b and c) were recently found to transit (Kane et al. 2020), prompting a photometric follow-up by the brand-new $CHaracterising\:ExOPlanets\:Satellite\:(CHEOPS)$. Here, we report that the outer planet d is also transiting, and measure its radius and mass to be $2.56\pm0.09$ $R_{\oplus}$ and $8.82\pm0.94$ $M_{\oplus}$, respectively. With its bright Sun-like star, long period, and mild irradiation ($\sim$5.7 times the irradiation of Earth), $\nu^2$ Lupi d unlocks a completely new region in the parameter space of exoplanets amenable to detailed characterization. We refine the properties of all three planets: planet b likely has a rocky mostly dry composition, while planets c and d seem to have retained small hydrogen-helium envelopes and a possibly large water fraction. This diversity of planetary compositions makes the $\nu^2$ Lupi system an excellent laboratory for testing formation and evolution models of low-mass planets.

The impact of cosmic reionization on the Ly$\alpha$ forest power spectrum has recently been shown to be significant even at low redshifts ($z \sim 2$). This memory of reionization survives cosmological time scales because high-entropy mean-density gas is heated to $\sim 3\times10^4$ K by reionization, which is inhomogeneous, and subsequent shocks from denser regions. In the near future, the first measurements of the Ly$\alpha$ forest 3D power spectrum will be very likely achieved by upcoming observational efforts such as the Dark Energy Spectroscopic Instrument (DESI). In addition to abundant cosmological information, these observations have the potential to extract the astrophysics of reionization from the Ly$\alpha$ forest. We forecast, for the first time, the accuracy with which the measurements of Ly$\alpha$ forest 3D power spectrum can place constraints on the reionization parameters with DESI. Specifically, we demonstrate that the constraints on the ionization efficiency, $\zeta$, and the threshold mass for haloes that host ionizing sources, $m_{\rm turn}$, will have the $1\sigma$ error at the level of $\zeta = 25.0 \pm 11.9$ and $\log_{10} (m_{\rm turn}/{\rm M}_\odot) = 8.7^{+0.35}_{-0.70}$, respectively. The Ly$\alpha$ forest 3D power spectrum will thus provide an independent, probably even earlier in detection with DESI, probe of reionization with a sensitivity only slightly worse than the upcoming 21~cm power spectrum measurement with the Hydrogen Epoch of Reionization Array (HERA).

Context. The eROSITA X-ray telescope onboard the Spectrum-Roentgen-Gamma (SRG) observatory combines a large field of view and collecting area in the energy range $\sim$0.2 to $\sim$8.0 keV with the capability to perform uniform scanning observations of large sky areas. Aims. SRG/eROSITA performed scanning observations of the $\sim$140 square degrees eROSITA Final Equatorial Depth Survey (eFEDS) field as part of its performance verification phase. The observing time was chosen to slightly exceed the depth of equatorial fields after the completion of the eROSITA all-sky survey. We present a catalog of detected X-ray sources in the eFEDS field providing source positions and extent information, as well as fluxes in multiple energy bands and document a suite of tools and procedures developed for eROSITA data processing and analysis, validated and optimized by the eFEDS work. Methods. A multi-stage source detection procedure was optimized and calibrated by performing realistic simulations of the eROSITA eFEDS observations. We cross-matched the eROSITA eFEDS source catalog with previous XMM-ATLAS observations, confirming excellent agreement of the eROSITA and XMM-ATLAS source fluxes. Result. We present a primary catalog of 27910 X-ray sources, including 542 with significant spatial extent, detected in the 0.2-2.3 keV energy range with detection likelihoods $\ge 6$, corresponding to a point source flux limit of $\approx 7 \times 10^{-15}$ erg/cm$^2$/s in the 0.5-2.0 keV energy band. A supplementary catalog contains 4774 low-significance source candidates with detection likelihoods between 5 and 6. In addition, a hard band sample of 246 sources detected in the energy range 2.3-5.0 keV above a detection likelihood of 10 is provided. A description of the dedicated data analysis software, calibration database and standard calibrated data products is provided in appendix.

The eROSITA Final Equatorial-Depth Survey has been carried out during the PV phase of the SRG/eROSITA telescope and completed in November 2019. This survey is designed to provide the first eROSITA-selected sample of galaxy clusters and groups and to test the predictions for the all-sky survey in the context of cosmological studies with clusters. In the 140 deg$^2$ area covered by eFEDS, 542 candidate clusters and groups are detected as extended X-ray sources, down to a flux of $\sim10^{-14} $erg/s/cm$^2$ in the soft band (0.5-2 keV) within 1'. In this work, we provide the catalog of candidate galaxy clusters and groups in eFEDS. We perform imaging and spectral analysis on the eFEDS clusters with eROSITA X-ray data, and study the properties of the sample. The clusters are distributed in the redshift range [0.01, 1.3], with the median redshift at 0.35. We obtain the ICM temperature measurement with $>2\sigma$ c.l. for $\sim$1/5 (102/542) of the sample. The average temperature of these clusters is $\sim$2 keV. Radial profiles of flux, luminosity, electron density, and gas mass are measured from the precise modeling of the imaging data. The selection function, the purity and completeness of the catalog are examined and discussed in detail. The contamination fraction is $\sim1/5$ in this sample, dominated by misidentified point sources. The X-ray Luminosity Function of the clusters agrees well with the results obtained from other recent X-ray surveys. We also find 19 supercluster candidates in eFEDS, most of which are located at redshifts between 0.1 and 0.5. The eFEDS cluster and group catalog provides a benchmark proof-of-concept for the eROSITA All-Sky Survey extended source detection and characterization. We confirm the excellent performance of eROSITA for cluster science and expect no significant deviations from our pre-launch expectations for the final All-Sky Survey.

The eROSITA Final Equatorial-Depth Survey (eFEDS), covering ~140 square degrees, was performed as part of the performance verification phase of the eROSITA telescope on board of the Russian-German satellite Spectrum-Roentgen-Gamma (SRG). In this paper we present the optical follow-up of 542 X-ray extent selected galaxy group and cluster candidates providing redshifts and cluster confirmation for the full sample. We use optical imaging data from the Hyper Suprime-Cam Subaru Strategic Program and from the Legacy Survey to run the cluster redshift and confirmation tool MCMF as well as the optical cluster finder CAMIRA at the location of the X-ray candidates. While providing redshift estimates for all 542 candidates, we construct an optically confirmed sample of 477 clusters and groups with a residual contamination of 6%. Of these, 470 (98.5%) are confirmed using MCMF and 7 systems are added through cross matching with spectroscopic group catalogs. Using observable to observable scaling and the applied confirmation threshold, we predict 8 +/- 2 real systems have been excluded with the MCMF cut required to build this low contamination sample. This number is in good agreement with the 7 systems recovered through cross matching. Thus, we expect those 477 systems to include >99% of all true clusters in the candidate list. Using an MCMF independent method, we confirm the catalog contamination of the confirmed subsample to be 6 +/- 3% and find 17 +/- 3% contamination for the full X-ray sample. The estimated contamination of the fulls sample is in agreement with MCMF dependent estimate of ~17% and the expectation from dedicated X-ray simulations of ~20%. We further present a sample of optically selected merging cluster candidates.

In Fall 2019, eROSITA on board of SRG observatory started to map the entire sky in X-rays. After the 4-year survey program, it will reach flux limits about 25 times deeper than ROSAT. During the SRG Performance Verification phase, eROSITA observed a contiguous 140 deg^2 area of the sky down to the final depth of the eROSITA all-sky survey ("eROSITA Final Equatorial-Depth Survey": eFEDS), with the goal of getting a census of the X-ray emitting populations that will be discovered over the entire sky. This paper presents the identification of the counterparts to the point-sources detected in eFEDS in the Main and Hard samples described in Brunner et al., and their multi-wavelength properties, including redshift. For the identification of the counterparts we combined the results from two independent methods NWAY and ASTROMATCH, trained on the multi-wavelength properties of a sample of 23k XMM-Newton sources detected in the DESI Legacy Imaging Survey DR8. Then spectroscopic redshifts and photometry from ancillary surveys are collated for the computation of photometric redshifts. The eFEDS sources with a reliable counterparts are 24774/27369 (90.5\%) in the Main sample and 231/246 (93.9\%) in the Hard sample, including 2514 (3) sources for which a second counterpart is equally likely. By means of reliable spectra, Gaia parallaxes, and/or multiwavelength properties we have classified the counterparts in both samples as 'Galactic' (2822) and 'extragalactic' (21952). For about 340 of the extragalactic sources we cannot rule out the possibility that they are unresolved clusters or belong to clusters. Inspection of the distributions of the X-ray sources in various optical/IR color-magnitude spaces reveal a rich variety of diverse classes of objects. The photometric redshifts are most reliable within the KiDS/VIKING area, where also deep near-infrared data is available.(abridged)

Stars are ubiquitous X-ray emitters and will be a substantial fraction of the X-ray sources detected in the on-going all-sky survey performed by the eROSITA instrument aboard the Spectrum Roentgen Gamma (SRG) observatory. We use the X-ray sources in the eROSITA Final Equatorial-Depth Survey (eFEDS) field observed during the SRG performance verification phase to investigate different strategies to identify the stars among other source categories. We focus here on Support Vector Machine (SVM) and Bayesian approaches, and our approaches are based on a cross-match with the Gaia catalog, which will eventually contain counterparts to virtually all stellar eROSITA sources. We estimate that 2060 stars are among the eFEDS sources based on the geometric match distance distribution, and we identify the 2060 most likely stellar sources with the SVM and Bayesian methods, the latter being named HamStars in the eROSITA context. Both methods reach completeness and reliability percentages of almost 90%, and the agreement between both methods is, incidentally, also about 90%. Knowing the true number of stellar sources allowed us to derive association probabilities $p_{ij}$ for the SVM method similar to the Bayesian method so that one can construct samples with defined completeness and reliability properties using appropriate cuts in $p_{ij}$. The thus identified stellar sources show the typical characteristics known for magnetically active stars, specifically, they are generally compatible with the saturation level, show a large spread in activity for stars of spectral F to G, and have comparatively high fractional X-ray luminosities for later spectral types.

Context: After the successful launch of the Spectrum-Roentgen-Gamma (SRG) mission in July 2019, eROSITA, the soft X-ray instrument aboard SRG, performed scanning observations of a large contiguous field, namely the eROSITA Final Equatorial Depth Survey (eFEDS), ahead of the planned four-year all-sky survey. eFEDS yielded a large sample of X-ray sources with very rich multi-band photometric and spectroscopic coverage. Aims: We present here the eFEDS Active Galactic Nuclei (AGN) catalog and the eROSITA X-ray spectral properties of the eFEDS sources. Methods: Using a Bayesian method, we perform a systematic X-ray spectral analysis for all eFEDS sources. The appropriate model is chosen based on the source classification and the spectral quality, and, in the case of AGN, including the possibility of intrinsic (rest-frame) absorption and/or soft excess emission. Hierarchical Bayesian modeling (HBM) is used to estimate the spectral parameter distribution of the sample. Results: X-ray spectral properties are presented for all eFEDS X-ray sources. There are 21952 candidate AGN, which comprise 79% of the eFEDS sample. Despite a large number of faint sources with low photon counts, our spectral fitting provides meaningful measurements of fluxes, luminosities, and spectral shapes for a majority of the sources. This AGN catalog is dominated by X-ray unobscured sources, with an obscured (logNH>21.5) fraction of 10% derived by HBM. The power-law slope of the catalog can be described by a Gaussian distribution of 1.94+-0.22. Above a photon counts threshold of 500, nine out of 50 AGN have soft excess detected. For the sources with blue UV to optical color (type-I AGN), the X-ray emission is well correlated with the UV emission with the usual anti-correlation between the X-ray to UV spectral slope {\alpha}_{OX} and the UV luminosity.

The 140 square degree Final Equatorial-Depth Survey (eFEDS) field, observed with the extended ROentgen Survey with an Imaging Telescope Array (eROSITA) aboard the Spectrum-Roentgen-Gamma (SRG) mission, provides a first look at the variable eROSITA sky. We analyze the intrinsic X-ray variability of the eFEDS sources, provide X-ray light curves and tables with variability test results in the 0.2-2.3 keV (soft) and 2.3-5.0 keV (hard) bands. respectively. We perform variability tests using the normalized excess variance and maximum amplitude variability methods as performed for the 2RXS catalogue and add results from the Bayesian excess variance and the Bayesian block methods. In total 65 sources have been identified as being significantly variable in the soft band. In the hard band only one source is found to vary significantly. For the most variable sources fits to stellar flare events reveal extreme flare properties. A few highly variable AGN have also been detected. About half of the variable eFEDS sources have been detected at X-rays with eROSITA for the first time. Comparison with 2RXS and XMM observations provide variability information on timescales of years to decades.

During the performance verification phase of the SRG/eROSITA telescope, the eROSITA Final Equatorial-Depth Survey (eFEDS) has been carried out. It covers a 140 deg$^2$ field located at 126$^\circ <$ R.A. $< 146^\circ$ and -3$^\circ <$ Dec. $< +6^\circ$ with a nominal exposure over the field of 2.2 ks. 542 candidate clusters were detected in this field, down to a flux limit $F_X \sim 10^{-14}$ erg s$^{-1}$ cm$^{-2}$ in the 0.5-2 keV band. In order to understand radio-mode feedback in galaxy clusters, we study the radio emission of brightest cluster galaxies of eFEDS clusters, and we relate it to the X-ray properties of the host cluster. Using LOFAR we identify 227 radio galaxies hosted in the BCGs of the 542 galaxy clusters and groups detected in eFEDS. We treat non-detections as radio upper limits. We analyse the properties of radio galaxies, such as redshift and luminosity distribution, offset from the cluster centre, largest linear size and radio power. We study their relation to the intracluster medium of the host cluster. We perform statistical tests to deal with upper limits on the radio luminosities. BCGs with radio-loud AGN are more likely to lie close to the cluster centre than radio-quiet BCGs. There is a clear relation between the cluster's X-ray luminosity and the radio power of the BCG. Statistical tests indicate that this correlation is not produced by selection effects in the radio band. We see no apparent link between largest linear size of the radio galaxy and central density of the host cluster. Converting the radio luminosity to kinetic luminosity, we find that radiative losses of the intracluster medium are in an overall balance with the heating provided by the central AGN. Finally, we tentatively classify our objects into disturbed and relaxed, and we show that the link between the AGN and the ICM apparently holds regardless of the dynamical state of the cluster.

The $\it{eROSITA}$ Final Equatorial Depth Survey (eFEDS), completed during the calibration and performance verification phase of the $\it{eROSITA}$ instrument on $\it{Spectrum\, Roentgen\, Gamma}$, delivers data at and beyond the final depth of the four-year $\it{eROSITA}$ all-sky survey (eRASS:8), $f_{0.5-2\,\text{ keV}}$ = $1.1\times10^{-14}$ erg s$^{-1}$ cm$^{2}$, over 140 deg$^{2}$. It provides the first view of normal galaxy X-ray emission from X-ray binaries (XRBs) and the hot interstellar medium at the full depth of eRASS:8. We use the Heraklion Extragalactic Catalogue (HECATE) of galaxies to correlate with eFEDS X-ray sources and identify 94 X-ray detected normal galaxies. We classify galaxies as star-forming, early-type, composite, and AGN using SDSS and 6dF optical spectroscopy. The eFEDS field harbours 37 normal galaxies: 36 late-type (star-forming) galaxies and 1 early-type galaxy. There are 1.9 times as many normal galaxies as predicted by scaling relations via SIXTE simulations, with an overabundance of late-type galaxies and a dearth of early-type galaxies. Dwarf galaxies with high specific star formation rate (SFR) have elevated L$_{\text{X}}$/SFR when compared with specific SFR and metallicity, indicating an increase in XRB emission due to low-metallicity. We expect that eRASS:8 will detect 12,500 normal galaxies, the majority of which will be star-forming, with the caveat that there are unclassified sources in eFEDS and galaxy catalogue incompleteness issues that could increase the actual number of detected galaxies over these current estimates. eFEDS observations detected a rare population of galaxies -- the metal-poor dwarf starbursts -- that do not follow known scaling relations. eRASS is expected to discover significant numbers of these high-redshift analogues, which are important for studying the heating of the intergalactic medium at high-redshift.

Context: The eROSITA X-ray telescope onboard the Spectrum-Roentgen-Gamma (SRG) satellite has started to observe new X-ray sources over the full sky at an unprecedented rate. Understanding the selection function of the source detection is important to the subsequent scientific analysis of the eROSITA catalogs. Aims: Through simulations, we test and optimize the eROSITA source detection procedures, and characterize the detected catalog quantitatively. Methods: Taking the eROSITA Final Equatorial-Depth Survey (eFEDS) as an example, we run extensive photon event simulations using our best knowledge of the instrument characteristics, the background spectrum, and the population of astronomical X-ray sources. We analyze the source detection results based on the origin of each photon. Results. The source detection procedure is optimized according to the source detection efficiency. We choose a two-pronged strategy to build the eFEDS X-ray catalogs, creating a main catalog using only the most sensitive band (0.2-2.3 keV) and an independent hard-band selected catalog using multi-band detection in a range up to 5 keV. From the mock catalogs (available with this paper), we measure the catalog completeness and purity, which can be used in both choosing the sample selection thresholds and in further studies of AGN and cluster demography.

The reliability of detecting source variability in sparsely and irregularly sampled X-ray light curves is investigated. This is motivated by the unprecedented survey capabilities of eROSITA onboard SRG, providing light curves for many thousand sources in its final-depth equatorial deep field survey. Four methods for detecting variability are evaluated: excess variance, amplitude maximum deviations, Bayesian blocks and a new Bayesian formulation of the excess variance. We judge the false detection rate of variability based on simulated Poisson light curves of constant sources, and calibrate significance thresholds. Simulations with flares injected favour the amplitude maximum deviation as most sensitive at low false detections. Simulations with white and red stochastic source variability favour Bayesian methods. The results are applicable also for the million sources expected in eROSITA's all-sky survey.

In the first months after the launch in July 2019, eROSITA onboard Spektr-RG (SRG) performed long-exposure observations in the regions around SN 1987A and SNR N132D in the Large Magellanic Cloud (LMC). We analyse the distribution and the spectrum of the diffuse X-ray emission in the observed fields to determine the physical properties of the hot phase of the interstellar medium (ISM). The eROSITA data are complemented by newly derived column density maps for the Milky Way and the LMC, 888 MHz radio continuum map from the Australian Square Kilometer Array Pathfinder (ASKAP), and optical images of the Magellanic Cloud Emission Line Survey (MCELS). We detect significant emission from thermal plasma with kT=0.2 keV in all the regions. There is also an additional higher-temperature emission component from a plasma with kT = 0.7 keV. In addition, non-thermal X-ray emission is significantly detected in the superbubble 30 Dor C. The absorbing column density NH in the LMC derived from the analysis of the X-ray spectra taken with eROSITA is consistent with the NH obtained from the emission of the cold medium over the entire area. Neon abundance is enhanced in the regions in and around 30 Dor and SN 1987A, indicating that the ISM has been chemically enriched by the young stellar population. Emission from the stellar cluster RMC 136 and the Wolf-Rayet stars RMC 139 and RMC 140 is best modelled with a high-temperature (kT>1 keV) non-equilibrium ionisation plasma emission and a non-thermal component with a photon index of {\Gamma} =1.3. In addition, the optical SNR candidate J0529-7004 is also detected with eROSITA and we thus confirm the source as an SNR.

The nearby young open cluster eta Chamaeleontis has been observed by eROSITA/SRG during its CalPV phase for 150 ks. The eROSITA data were taken in the field-scan mode, an observing mode of Spectrum-Roentgen-Gamma (SRG) that follows a rectangular grid-like pattern, here covering a 5x5 deg field with an exposure depth of about 5 ks. We study the known members in X-rays and search for potential new members of the anticipated dispersed low-mass cluster population. Detected sources were identified by cross-matching X-ray sources with Gaia and 2MASS, and young stars were identified by their X-ray activity, the position in the color-magnitude diagram, and by their astrometric and kinematic properties. X-ray-luminosities, light curves, and spectra of cluster members were obtained and compared with previous X-ray data. Literature results of other member searches were used to verify our new member candidates in the observed field. We determine X-ray properties of virtually all known eta Cha members and identify five additional stellar systems that show basically identical characteristics, but are more dispersed. Four of them were previously proposed as potential members; this status is supported by our X-ray study. Based on their spatial distribution, further members are expected beyond the sky region we surveyed. The identified stellar systems very likely belong to the ejected halo population, which brings the total number of eta Cha cluster members to at least 23.

SN 1987A is the closest observed supernova in the last four centuries and provides a unique opportunity to witness the birth and evolution of a supernova remnant. The source has been monitored by XMM Newton EPIC-pn from 2007--2020. SRG/eROSITA also observed the source during its commissioning phase and the first light in Sept. and Oct. 2019. We investigated the spectral and flux evolution of SN 1987A in X-rays over the last fourteen years up to Nov. 2020 using XMM-Newton and eROSITA observations. We performed a detailed spectral analysis using a three-component plane-parallel shock model and analysed and modelled the EPIC-pn monitoring and eROSITA observations in a consistent manner. This paper reports a complete and most up to date flux evolution of SN 1987A in the soft (0.5-2 keV) and hard (3-10 keV) X-ray band. The flux in the soft band flattened around 9424 d and then displayed a turnover between 10142-10493 d after which it showed a continued decline. Around the same time, a break in the hard-band flux time evolution slope was detected. This implies that the blast wave has now passed beyond the dense structures of the equatorial ring and is expanding further into more tenuous circumstellar medium. The temporal evolution of the normalizations of the three shock components match well the results of hydrodynamical simulations predicting a blue supergiant progenitor scenario. The trend at recent epochs indicate that the emission caused by the forward shock after leaving the equatorial ring and by the reverse shock in the ejecta is becoming more dominant now. The elemental abundances in the hot plasma component are significantly higher than those in the `cooler' one, indicating its origin from the reverse shock propagating into the ejecta.

(abridged version) We present a detailed spectroscopic and timing analysis of X-ray observations of the bright radio-to-gamma-ray emitting pulsar PSR B0656+14, which were obtained simultaneously with eROSITA and XMM-Newton during the Calibration and Performance Verification phase of the Spektrum-Roentgen-Gamma mission (SRG) for 100 ks. Using XMM-Newton and NICER we firstly established an X-ray ephemeris for the time interval 2015 to 2020, which connects all X-ray observations in this period without cycle count alias and phase shifts. The mean eROSITA spectrum clearly reveals an absorption feature originating from the star at 570 eV with a Gaussian sigma of about 70 eV, tentatively identified earlier in a long XMM-Newton observation (Arumugasamy et al. 2018). A second absorption feature, described here as an absorption edge, occurs at 260-265 eV. It could be of atmospheric or of instrumental origin. These absorption features are superposed on various emission components, phenomenologically described as the sum of hot (120 eV) and cold (65 eV) blackbody components, both of photospheric origin, and a power-law with photon index Gamma=2. The phase-resolved spectroscopy reveals that the Gaussian absorption line at 570 eV is clearly present throughout ~60% of the spin cycle. The visibility of the line strength coincides in phase with the maximum flux of the hot blackbody. We also present three families of model atmospheres: a magnetised atmosphere, a condensed surface, and a mixed model, which were applied to the mean observed spectrum and whose continuum fit the observed data well. The atmosphere model, however, predicts too short distances. For the mixed model, the Gaussian absorption may be interpreted as proton cyclotron absorption in a field as high as 10^14 G, which is significantly higher than that derived from the moderate observed spin-down.

The galaxy cluster Abell 3266 is one of the X-ray brightest in the sky and is a well-known merging system. Using the ability of the eROSITA telescope onboard SRG (Spectrum R\"ontgen Gamma) to observe a wide field with a single pointing, we analyse a new observation of the cluster out to a radius of R_200. The X-ray images highlight substructures present in the cluster, including the northeast-southwest merger seen in previous ASCA, Chandra and XMM-Newton data, a merging group towards the northwest and filamentary structures between the core and one or more groups towards the west. We compute spatially-resolved spectroscopic maps of the thermodynamic properties of the cluster, including the metallicity. The merging subclusters are seen as low entropy material within the cluster. The filamentary structures could be the rims of a powerful AGN outburst, or most likely material stripped from the western group(s) as they passed through the cluster core. Seen in two directions is a pressure jump at a radius of 1.1 Mpc consistent with a shock with a Mach number of ~1.5-1.7. The eROSITA data confirm that the cluster is not a simple merging system, but is made up of several subclusters which are merging or will shortly merge. For the first time we find a radio halo associated with the system detected in GLEAM data. We compute a hydrostatic mass from the eROSITA data, finding good agreement with a previous XMM-Newton result. With this pointing we detect several extended sources, where we find for seven of them secure associations between z=0.36-1.0; i.e., background galaxy groups and clusters, highlighting the power of eROSITA to find such systems.

We present the results of the analysis of five observations of the globular clutser 47 Tucanae (47 Tuc) with eROSITA (extended Roentgen Survey with an Imaging Telescope Array) on board Spektrum-Roentgen-Gamma (Spektr-RG, SRG). The aim of the work is the study of the X-ray population in the field of one of the most massive globular clusters in our Milky Way. We focused on the classification of point-like sources in the field of 47 Tuc. The unresolved dense core of 47~Tuc (1'.7 radius) and also the sources, which show extended emission are excluded in this study. We applied different methods of X-ray spectral and timing analysis together with multi wavelength studies for the classification of the X-rays sources in the field of 47 Tuc. We detected 888 point-like sources in the energy range of 0.2-5.0 keV. We identified 92 background AGNs and 26 foreground stars. One of the foreground stars is classified as a variable M~dwarf. We also classified 23 X-ray sources as members of 47 Tuc, including 13 symbiotic stars, 3 quiescent low mass X-ray binaries, one millisecond pulsar candidate, and one cataclysmic variable. There are also 4 X-ray sources, which can be either a cataclysmic variable or a contact binary. Moreover, we calculated the X-ray luminosity function of 47 Tuc X-ray sources within a radius of 18'.8. It shows that the main population of X-ray sources in 47 Tuc has a luminosity <10$^{32}$erg s$^{-1}$ in the energy range of 0.5-2.0 keV. These sources can mainly be candidates for quiescent low mass X-ray binaries and different types of accreting white dwarfs, especially symbiotic stars.

Context. The Small Magellanic Cloud (SMC) is hosting many known high-mass X-ray binaries, all but one (SMC X-1) having Be companion stars. Through the calibration and verification phase of eROSITA on board the SRG spacecraft, the Be/X-ray binary XMMU J010429.4-723136 was in the field of view during observations of the supernova remnant 1E0102.2-7219, used as calibration standard. Aims. We report here a time and spectral analysis of XMMU J010429.4-723136, based on two eROSITA observations of the field, performed on 2019 November 7-9. We also reanalyse the OGLE light curve for that source, in order to determine the orbital period. Methods. The search for pulsations (from the X-ray data) and for the orbital period (from the OGLE data) is done via Lomb-Scargle periodogram analysis. X-ray spectral parameters and fluxes are retrieved from the best-fit model. Results. We detected for the first time, the pulsations of XMMU J010429.4-723136 at a period of 164 s, and therefore designate the source as SXP 164. From the spectral fitting, we derive a source flux of 1x10e-12 erg s-1 cm-2 for both observations, corresponding to a luminosity of 4x10e35 erg s-1 at the distance of the SMC. Furthermore, reanalysing the OGLE light curve including the latest observations, we found a significant periodic signal at 22.3d likely being the orbital period, which is shorter than the previously reported values.

We present simultaneous multiwavelength observations of the ultra-active star AB~Doradus obtained in the X-ray range with the eROSITA instrument on board the Russian--German Spectrum-Roentgen-Gamma mission (SRG), and in the optical range obtained with the Transiting Exoplanet Survey Satellite (TESS). Thanks to its fortuitous location in the vicinity of the southern ecliptic pole, AB~Dor was observed by these missions simultaneously for almost 20 days. With the hitherto obtained data we study the long-term evolution of the X-ray flux from AB~Dor and the relation between this observable and the photospheric activity of its spots. Over the 1.5 years of eROSITA survey observations, the "quiescent" X-ray flux of AB~Dor has not changed, and furthermore it appears unrelated to the photospheric modulations observed by TESS. During the simultaneous eROSITA and TESS coverage, an extremely large flare event with a total energy release of at least 4 $\times$ 10$^{36}$ erg in the optical was observed, the largest ever seen on AB~Dor. We show that the total X-ray output of this flare was far smaller than this, and discuss whether this maybe a general feature of flares on late-type stars.

We report the identification and follow-up of the transient SRG 062340.2-265715 detected with both instruments on board the Spektrum-Roentgen-Gamma mission. Optical spectroscopy of the G=12.5 counterpart firmly classifies the object as a novalike cataclysmic variable (CV) at a distance of 495 pc. A highly significant TESS period of 3.941 hours, tentatively identified with the orbital period of the binary, could not be found when the object was reobserved with TESS two years later. The newer high-cadence TESS data revealed quasi-periodic oscillations around 25 min, while ground-based photometry indicated periodic variability at 32 min. Located in very sparsely populated regions of color-magnitude diagrams involving X-ray and optical magnitudes and colors, the new object could be an X-ray underluminous magnetic CV, an intermediate polar, or an overluminous nonmagnetic CV. The lack of uniquely identified spin and orbital periods prevents a final classification. The site of X-ray production in the system, L(X, bol) = 4.8 x 10^{32} erg/s, remains to be understood given its high variability on long and short timescales.

Ongoing all-sky surveys by the the eROSITA and the Mikhail Pavlinsky ART-XC telescopes on-board the Spectrum Roentgen Gamma (SRG) mission have already revealed over a million of X-ray sources. One of them, SRGA J124404.1-632232/SRGU J124403.8-632231, was detected as a new source in the third (of the planned eight) consecutive X-ray surveys by ART-XC. Based on the properties of the identified optical counterpart it was classified as a candidate X-ray binary (XRB). We report on the follow-up observations of this source with Nuclear Spectroscopic Telescope Array (NuSTAR), Neil Gehrels Swift Observatory (Swift), and the Southern African Large Telescope (SALT), which allowed us to unambiguously confirm the initial identification and establish SRGU J124403.8-632231 as a new X-ray pulsar with a spin period of ~538 s and a Be-star companion, making it one of the first Galactic X-ray pulsars discovered by SRG.

We report the discovery of a bright (V ~ 15), eclipsing, two-pole accreting magnetic cataclysmic variable (CV), a polar, as counterpart of the eROSITA and Gaia transients eRASSt 192932.9-560346 and Gaia21bxo. Frequent large amplitude changes of its brightness at X-ray and optical wavelengths by more than 4 magnitudes was indicative of a CV nature of the source. Identification spectra obtained with the 10m SALT telescope revealed the typical features of a magnetic CV, strong, broad HeI, HeII and hydrogen Balmer emission lines superposed on a blue continuum. Time-resolved photoelectric polarimetry revealed circular polarization to vary from -20% to +20%, and linear polarization from 0% to 10% confirming the system to be magnetic CV of the polar subclass. High cadence photometry revealed deep, structured eclipses, indicating that the system is a two pole accretor. The orbital period determined from the eclipse times is 92.5094 +- 0.0002 min. The X-ray spectrum is thermal only and the implied luminosity is L_X=2.2 x 10^(31) erg/s at the Gaia-determined distance of 376 pc.

In 2018, XMM-Newton observed the awakening in X-rays of the Be/X-ray binary (Be/XRB) A0538-66. It showed bright and fast flares close to periastron with properties that had never been observed in other Be/XRBs before. We report the results from the observations of A0538-66 collected during the first all-sky survey of eROSITA, an X-ray telescope (0.2-10 keV) on board the Spektrum-Roentgen-Gamma (SRG) satellite. eROSITA caught two flares within one orbital cycle at orbital phases $\phi = 0.29$ and $\phi = 0.93$ (where $\phi=0$ corresponds to periastron), with peak luminosities of $\sim 2-4 \times 10^{36}$ erg/s (0.2-10 keV) and durations of $42 \leq \Delta t_{\rm fl} \leq 5.7\times 10^4$ s. The flare observed at $\phi \approx 0.29$ shows that the neutron star can accrete considerably far from periastron, although it is expected to be outside of the circumstellar disk, thus providing important new information about the plasma environment surrounding the binary system. We also report the results from the photometric monitoring of A0538-66 carried out with the REM, OGLE, and MACHO telescopes from January 1993 until March 2020. We found that the two sharp peaks that characterize the orbital modulation in the optical occur asymmetrically in the orbit, relative to the position of the donor star.

We investigated the origin and gas properties of a simulated galaxy cluster pair, to connect simulation predictions to the SRG/eROSITA X-ray observations of the Abell 3391/95 field. The simulated system has been extracted from the (352 cMpc/h)^3 volume of the Magneticum Pathfinder cosmological simulations at z=0.07. We tracked back in time the main progenitors of the pair clusters and surrounding groups to study the assembly history of the system and its evolution. Similarly to the observed A3391/95 system, the simulated pair is embedded in a complex network of gas filaments, with structures aligned over more than 20 projected Mpc and the whole region collapsing towards the central overdense node. The spheres of influence (3*R200) of the two main clusters already overlap at z=0.07, but their virial boundaries are still physically separated. The diffuse gas located in the interconnecting bridge closely reflects the WHIM, with typical temperature of ~1 keV and overdensity $\sim 100$, with respect to the mean baryon density of the Universe, and lower enrichment level compared to the ICM in clusters. We find that most of the bridge gas collapsed from directions roughly orthogonal to the intra-cluster gas accretion directions, and its origin is mostly unrelated to the two cluster progenitors. We find clear signatures in the surrounding groups of infall motion towards the pair, such as significant radial velocities and slowdown of gas compared to dark matter. These findings further support the picture of the Northern Clump (MCXC J0621.7-5242) cluster infalling along a cosmic gas filament towards Abell 3391, possibly merging with it. We conclude that, in such a configuration, the pair clusters of the A3391/95-like system are in a pre-merger phase, and did not interact yet. The diffuse gas in the interconnecting bridge is mostly warm filament gas, rather than tidally-stripped cluster gas.

SRG/eROSITA PV observations revealed the A3391/95 cluster system and the Northern Clump (MCXC J0621.7-5242 galaxy cluster) are aligning along a cosmic filament in soft X-rays, similarly to what has been seen in simulations before. We aim to understand the dynamical state of the Northern Clump as it enters the atmosphere ($3\times R_{200}$) of A3391. We analyzed joint eROSITA, XMM-Newton, and Chandra observations to probe the morphological, thermal, and chemical properties of the Northern Clump from its center out to a radius of 1077 kpc ($1.1\times R_{200}$). We utilized the ASKAP/EMU radio data, DECam optical image, and Planck y-map to study the influence of the WAT radio source on the Northern Clump central ICM. From the Magneticum simulation, we identified an analog of the A3391/95 system along with an infalling group resembling the Northern Clump. The Northern Clump is a WCC cluster centered on a WAT radio galaxy. The gas temperature over $0.2-0.5R_{500}$ is $k_BT_{500}=1.99\pm0.04$ keV. We employed the $M-T$ scaling relation and obtained a mass estimate of $M_{500}=(7.68\pm0.43)\times10^{13}M_{\odot}$ and $R_{500}=(636\pm12)$ kpc. Its atmosphere has a boxy shape and deviates from spherical symmetry. We identify a southern surface brightness edge, likely caused by subsonic motion relative to the filament gas. At $\sim\! R_{500}$, the southern atmosphere appears to be 42% hotter than its northern atmosphere. We detect a downstream tail pointing towards north with a projected length of $\sim318$ kpc, plausibly the result of ram pressure stripping. The analog group in the Magneticum simulation is experiencing changes in its gas properties and a shift between the position of the halo center and that of the bound gas, while approaching the main cluster pair.

The X-ray telescope eROSITA onboard the newly launched SRG mission serendipitously observed the galaxy cluster A3408 ($z=0.0420$) during the PV observation of the AGN 1H0707-495. Despite its brightness and large extent, it has not been observed by any modern X-ray observatory. A neighbouring cluster in NW direction, A3407 ($z=0.0428$), appears to be close at least in projection ($\sim 1.7$ Mpc). This cluster pair could be in a pre- or post-merger state. We perform a detailed X-ray analysis of A3408. We construct particle background subtracted and exposure corrected images and surface brightness profiles in different sectors. The spectral analysis is performed out to $1.4r_{500}$. Additionally, a temperature map is presented depicting the distribution of the ICM temperature. Furthermore, we make use of data from the RASS to estimate some bulk properties of A3408 and A3407, using the growth curve analysis method and scaling relations. The imaging analysis shows a complex morphology of A3408 with a strong elongation in SE-NW direction. This is quantified by comparing the surface brightness profiles of the NW, SW, SE and NE directions, where the NW and SE directions show a significantly higher surface brightness compared to the other directions. We determine a gas temperature ${\rm k_B}T_{500}=(2.23\pm0.09)$ keV. The T-profile reveals a hot core within $2'$ of the emission peak, ${\rm k_B}T=3.04^{+0.29}_{-0.25}$ keV. Employing a M-T relation, we obtain $M_{500}=(9.27\pm0.75)\times 10^{13}M_{\odot}$ iteratively. The $r_{200}$ of A3407 and A3408 are found to overlap in projection which makes ongoing interactions plausible. The 2d T-map reveals higher temperatures in W than in E direction. A3407 and A3408 are likely in a pre-merger state, affecting the ICM properties, i.e., increased temperatures in the direction of A3407 indicate adiabatic compression or shocks due to the interaction.

The most nearby clusters are the best places to study physical and enrichment effects in the faint cluster outskirts. A3158 located at z=0.059 is quite extended with a characteristic radius r$_{200}$=23.95 arcmin. In 2019, A3158 was observed as a calibration target in a pointed observation with the eROSITA telescope onboard the SRG mission. We determined 1d temperature, abundance and normalisation profiles from eROSITA and XMM-Newton and Chandra data as well as 2d maps of temperature distribution from eROSITA data. The velocity dispersion was determined and the cluster mass was calculated. The overall temperature was measured to be 4.725$\pm$ 0.035 keV. The profiles of eROSITA all agree on a ~10% level with those determined with XMM-Newton and Chandra data. From the temperature map we see that the cluster lacks a cool core, as noted before. The presence of a previously detected off-centre cool clump West of the central cluster region is observed. Furthermore there is a bow shaped edge near the location of the cool gas clump West of the cluster centre. An extension of gas is detected for the first time further out in the West. The velocity dispersion of the cluster was measured to be 1058$\pm$41 km s$^{-1}$. The total mass was determined as $M_{200}$=1.38$\pm$ 0.25x10$^{15}$ $M_{\odot}$. The mass estimate from the M-T relation is significantly lower at M$_{200}$=5.09$\pm$ 0.59x10$^{14}M_{\odot}$. An extended X-ray source located South of the cluster also coincides with a galaxy overdensity with redshifts in the range 0.05<z<0.07. These are indications that the cluster may be undergoing merger activity. Another extended source located North of the cluster is detected in X-rays and coincides with an overdensity of galaxies with redshifts in the range of 0.070<z<0.077. This is likely a background cluster not related to A3158. Additionally a known SPT cluster at z=0.53 was detected.

Stellar activity is inherently time variable, therefore simultaneous measurements are necessary to study the correlation between different activity indicators. In this study we compare X-ray fluxes measured within the first all-sky survey conducted by the extended ROentgen Survey with an Imaging Telescope Array (eROSITA) instrument on board the Spectrum-Roentgen-Gamma (SRG) observatory to Ca II H & K, excess flux measurements R+, using observations made with the robotic TIGRE telescope. We created the largest sample of simultaneous X-ray and spectroscopic Ca II H & K observations of late-type stars obtained so far, and in addition, previous measurements of Ca II H & K for all sample stars were obtained. We find the expected correlation between our log(L_X/L_bol) to log(R+) measurements, but when the whole stellar ensemble is considered, the correlation between coronal and chromospheric activity indicators does not improve when the simultaneously measured data are used. A more detailed analysis shows that the correlation of log(L_X/L_bol) to log(R+) measurements of the pseudo-simultaneous data still has a high probability of being better than that of a random set of non-simultaneous measurements with a long time baseline between the observations. Cyclic variations on longer timescales are therefore far more important for the activity flux-flux relations than short-term variations in the form of rotational modulation or flares, regarding the addition of "noise" to the activity flux-flux correlations. Finally, regarding the question of predictability of necessarily space-based log(L_X/L_bol) measurements by using ground-based chromospheric indices, we present a relation for estimating log(L_X/L_bol) from R+ values and show that the expected error in the calculated minus observed (C-O) log(L_X/L_bol) values is 0.35 dex.

We present the first X-ray detections of ultracool dwarfs (UCDs) from the first all-sky survey of the extended ROentgen Survey with an Imaging Telescope Array (eROSITA) onboard the Russian Spektrum-Roentgen-Gamma (SRG) mission. We use three publicly available input catalogs of spectroscopically confirmed UCDs and Gaia-selected UCD candidates that together comprise nearly 20000 objects. We first extracted all X-ray sources from the catalog of the first eROSITA survey, eRASS1, that have a UCD or candidate within three times their positional uncertainty. Then we examined all Gaia objects in the vicinity of these 96 X-ray sources and we associated them to the most plausible counterpart on the basis of their spatial separation to the X-ray position and their multiwavelength properties. This way we find 40 UCDs that have a secure identification with an X-ray source and 18 plausible UCD X-ray emitters. Twenty-one of these X-ray emitting UCDs have a spectroscopic confirmation, while the others have been selected based on Gaia photometry and we computed spectral types from the G-J color. The spectral types of the X-ray emitting UCDs and candidates range between M5 and M9, and the distances range from 3.5 to 190 pc. The majority of the UCDs from the eRASS1 sample show a ratio of X-ray to bolometric luminosity well above the canonical saturation limit of log (Lx/Lbol) ~ -3. For the two most extreme outliers, we showed through an analysis of the eRASS1 light curve that these high values are due to flaring activity. The X-ray spectra of the two brightest objects both reveal an emission-measure weighted plasma temperature of kT ~ 0.75 keV. These observations demonstrate the potential of eROSITA for advancing our knowledge on the faint coronal X-ray emission from UCDs by building statistical samples for which the average X-ray brightness, flares, and coronal temperatures can be derived.

We present the first study of nearby M dwarfs with the ROentgen Survey with an Imaging Telescope Array (eROSITA) on board the Russian Spektrum-Roentgen-Gamma mission (SRG). To this end we extracted the Gaia DR2 data for the ~9000 nearby M dwarfs in the superblink proper motion catalog and calculated their stellar parameters from empirical relations with optical-IR colors. We cross-matched this catalog with the eROSITA Final Equatorial Depth Survey (eFEDS) and the first eROSITA all-sky survey (eRASS1). Our sample consists of 704 stars (SpT = K5-M7). This unprecedented data base for X-ray emitting M dwarfs allowed to quantitatively constrain the mass dependence of the X-ray luminosity, and to determine the change in the activity level with respect to pre-main-sequence stars. We also combined these data with the Transiting Exoplanet Survey Satellite (TESS) observations that are available for 501 of 704 X-ray detected M dwarfs and determined the rotation period for 180 of them. With the joint eROSITA-TESS sample, and combining it with our historical X-ray and rotation data for M dwarfs, we examined the mass dependence in the saturated regime of the rotation-activity relation. A first comparison of eROSITA hardness ratios and spectra shows that 65% of our X-ray detected M dwarfs have coronal temperatures of $\sim 0.5$ keV. We investigated their long-term X-ray variability by comparing the eRASS1 and ROSAT all-sky survey (RASS) measurements. Evidence for X-ray flares is found in various parts of our analysis: directly from inspection of the eFEDS light curves, in the relation between RASS and eRASS1 X-ray luminosities, and in stars displaying X-ray emission hotter than the bulk of the sample according to the hardness ratios. Finally, we point out the need of X-ray spectroscopy for more M dwarfs to study the coronal temperature-luminosity relation, not well constrained by our eFEDS results.

We present the results of the first X-ray all-sky survey (eRASS1) performed by the eROSITA instrument on board the Spectrum-Roentgen-Gamma (SRG) observatory of the Sco-Cen OB association. Bona fide Sco-Cen member stars are young and are therefore expected to emit X-rays at the saturation level. The sensitivity limit of eRASS1 makes these stars detectable down to about a tenth of a solar mass. By cross-correlating the eRASS1 source catalog with the Gaia EDR3 catalog, we arrive at a complete identification of the stellar (i.e., coronal) source content of eROSITA\ in the Sco-Cen association, and in particular obtain for the first time a 3D view of the detected stellar X-ray sources. Focusing on the low-mass population and placing the optical counterparts identified in this way in a color-magnitude diagram, we can isolate the young stars out of the detected X-ray sources and obtain age estimates of the various Sco-Cen populations. A joint analysis of the 2D and 3D space motions, the latter being available only for a smaller subset of the detected stellar X-ray sources, reveals that the space motions of the selected population show a high degree of parallelism, but there is also an additional population of young, X-ray emitting and essentially cospatial stars that appears to be more diffuse in velocity space. Its nature is currently unclear. We argue that with our procedures, an identification of almost the whole stellar content of the Sco-Cen association will become possible once the final Gaia and eROSITA catalogs are available by the end of this decade. We furthermore call into question any source population classification scheme that relies on purely kinematic selection criteria.

High-energy irradiation is a driver for atmospheric evaporation and mass loss in exoplanets. This work is based on data from eROSITA, the soft X-ray instrument aboard SRG (Spectrum Roentgen Gamma) mission, as well as archival data from other missions, we aim to characterise the high-energy environment of known exoplanets and estimate their mass loss rates. We use X-ray source catalogues from eROSITA, XMM-Newton, Chandra and ROSAT to derive X-ray luminosities of exoplanet host stars in the 0.2-2 keV energy band with an underlying coronal, i.e. optically thin thermal spectrum. We present a catalogue of stellar X-ray and EUV luminosities, exoplanetary X-ray and EUV irradiation fluxes and estimated mass loss rates for a total of 287 exoplanets, 96 among them being characterised for the first time from new eROSITA detections. We identify 14 first time X-ray detections of transiting exoplanets that are subject to irradiation levels known to cause observable evaporation signatures in other exoplanets, which makes them suitable targets for follow-up observations.

In a previous paper we presented the results of applying machine learning to classify whether an HI 21-cm absorption spectrum arises in a source intervening the sight-line to a more distant radio source or within the host of the radio source itself. This is usually determined from an optical spectrum giving the source redshift. However, not only will this be impractical for the large number of sources expected to be detected with the Square Kilometre Array, but bright optical sources are the most ultra-violet luminous at high redshift and so bias against the detection of cool, neutral gas. Adding another 44, mostly newly detected absorbers, to the previous sample of 92, we test four different machine learning algorithms, again using the line properties (width, depth and number of Gaussian fits) as features. Of these algorithms, three gave a some improvement over the previous sample, with a logistic regression model giving the best results. This suggests that the inclusion of further training data, as new absorbers are detected, will further increase the prediction accuracy above the current 80%. We use the logistic regression model to classify the z = 0.42 absorption towards PKS 1657-298 and find this to be associated, which is consistent with a previous study which determined a similar redshift from the K-band magnitude-redshift relation.

Coronal loop observations have existed for many decades yet the precise shape of these fundamental coronal structures is still widely debated since the discovery that they appear to undergo negligible expansion between their footpoints and apex. In this work a selection of eight EUV loops and their twenty-two sub-element strands are studied from the second successful flight of NASA's High resolution Coronal Imager (Hi-C 2.1). Four of the loops correspond to open fan structures with the other four considered to be magnetically closed loops. Width analysis is performed on the loops and their sub-resolution strands using our method of fitting multiple Gaussian profiles to cross-sectional intensity slices. It is found that whilst the magnetically closed loops and their sub-element strands do not expand along their observable length, open fan structures may expand an additional 150% of their initial width. Following recent work, the Pearson correlation coefficient between peak intensity and loop/strand width are found to be predominantly positively correlated for the loops (~88%) and their sub-element strands (~80%). These results align with the hypothesis of Klimchuk & DeForest that loops and - for the first time - their sub-element strands have approximately circular cross-sectional profiles.

We report the results of combined analyses of X-ray and optical data of two galaxy clusters, CL 0024$+$1654 and RX J0152.7$-$1357 at redshift $z = 0.395$ and $z = 0.830$, respectively, which offer a holistic physical description of the two clusters. Our X-ray analysis yields temperature and density profile of the gas in the intra-cluster medium (ICM). Using optical photometric and spectroscopic data, complemented with mass distribution from gravitational lensing study, we investigate any possible correlation between the physical properties of the galaxy members, i.e., their color, morphology, and star formation rate (SFR) with their environments. We quantify the properties of the environment of each galaxy by galaxy number density, ICM temperature, and mass density. Although our result shows that the two clusters exhibit a weaker correlation compared to relaxed clusters, it still confirms the significant effect of the ICM on the SFRs in the galaxies. Various physical mechanisms are suggested to explain the relation between the properties of galaxies and their environment.

We present a study of the influence of magnetic field strength and morphology in Type Ia Supernovae and their late-time light curves and spectra. In order to both capture self-consistent magnetic field topologies as well evolve our models to late times, a two stage approach is taken. We study the early deflagration phase (1s) using a variety of magnetic field strengths, and find that the topology of the field is set by the burning, independent of the initial strength. We study late time (~1000 days) light curves and spectra with a variety of magnetic field topologies, and infer magnetic field strengths from observed supernovae. Lower limits are found to be 106G. This is determined by the escape, or lack thereof, of positrons that are tied to the magnetic field. The first stage employs 3d MHD and a local burning approximation, and uses the code Enzo. The second stage employs a hybrid approach, with 3D radiation and positron transport, and spherical hydrodynamics. The second stage uses the code HYDRA. In our models, magnetic field amplification remains small during the early deflagration phase. Late-time spectra bear the imprint of both magnetic field strength and morphology. Implications for alternative explosion scenarios are discussed.

Emission line observations together with photoionization models provide important information about the ionization mechanisms, densities, temperatures, and metallicities in AGN-ionized gas. Photoionization models usually assume Maxwell-Boltzmann (M-B) electron energy distributions (EED), but it has been suggested that using kappa distributions may be more appropriate and could potentially solve the discrepancies in temperatures and abundances found in HII regions and Planetary Nebulae (PNe). We consider the impact of the presence of kappa distributions in photoionized nebulae associated with AGN and study how this might affect spectral modelling and abundance analyses for such regions. Using the photoionization code MAPPINGS 1e we compute models adopting M-B and kappa distributions of electron energies, and compare the behaviour of emission line ratios for different values of kappa, gas metallicity, density, ionization parameter and SED slope. We find that the choice of EED can have a large impact on some UV and optical emission lines emitted by photoionized nebulae associated with AGN, and that the impact of adopting a kappa distribution is strongly dependent on gas metallicity and ionization parameter. We compile a sample of line ratios for 143 type 2 AGN and compare our models against the observed line ratios. We find that for 98 objects kappa distributions provide a better fit to the observed line ratios than M-B distributions. In addition, we find that adopting kappa-distributed electron energies results in significant changes in the inferred gas metallicity and ionization parameter in a significant fraction of objects.

We use the forward modeling approach to galaxy clustering combined with the likelihood from the effective-field theory of large-scale structure to measure assembly bias, i.e. the dependence of halo bias on properties beyond the total mass, in the linear ($b_1$) and second order bias parameters ($b_2$ and $b_{K^2}$) of dark matter halos in $N$-body simulations. This is the first time that assembly bias in the tidal bias parameter $b_{K^2}$ is measured. We focus on three standard halo properties: the concentration $c$, spin $\lambda$, and sphericity $s$, for which we find an assembly bias signal in $b_{K^2}$ that is opposite to that in $b_1$. Specifically, at fixed mass, halos that get more (less) positively biased in $b_1$, get less (more) negatively biased in $b_{K^2}$. We also investigate the impact of assembly bias on the $b_2(b_1)$ and $b_{K^2}(b_1)$ relations, and find that while the $b_2(b_1)$ relation stays roughly unchanged, assembly bias strongly impacts the $b_{K^2}(b_1)$ relation. This impact likely extends also to the corresponding relation for galaxies, which motivates future studies to design better priors on $b_{K^2}(b_1)$ for use in cosmological constraints from galaxy clustering data.

Coronal mass ejections (CMEs) that exhibit weak or no eruption signatures in the low corona, known as stealth CMEs, are problematic as upon arrival at Earth they can lead to geomagnetic disturbances that were not predicted by space weather forecasters. We investigate the origin and eruption of a stealth event that occurred on 2015 January 3 that was responsible for a strong geomagnetic storm upon its arrival at Earth. To simulate the coronal magnetic field and plasma parameters of the eruption we use a coupled approach. This approach combines an evolutionary nonlinear force-free field model of the global corona with a MHD simulation. The combined simulation approach accurately reproduces the stealth event and suggests that sympathetic eruptions occur. In the combined simulation we found that three flux ropes form and then erupt. The first two flux ropes, which are connected to a large AR complex behind the east limb, erupt first producing two near-simultaneous CMEs. These CMEs are closely followed by a third, weaker flux rope eruption in the simulation that originated between the periphery of AR 12252 and the southern polar coronal hole. The third eruption coincides with a faint coronal dimming, which appears in the SDO/AIA 211 A observations, that is attributed as the source responsible for the stealth event and later the geomagnetic disturbance at 1 AU. The incorrect interpretation of the stealth event being linked to the occurrence of a single partial halo CME observed by LASCO/C2 is mainly due to the lack of STEREO observations being available at the time of the CMEs. The simulation also shows that the LASCO CME is not a single event but rather two near-simultaneous CMEs. These results show the significance of the coupled data-driven simulation approach in interpreting the eruption and that an operational L5 mission is crucial for space weather forecasting.

We present ground-based optical transmission spectroscopy of the low-density hot Jupiter WASP-88b covering the wavelength range 4413-8333 {\AA} with the FORS2 spectrograph on the Very Large Telescope. The FORS2 white light curves exhibit a significant time-correlated noise which we model using a Gaussian Process and remove as a wavelength-independent component from the spectroscopic light curves. We analyse complementary photometric observations from the Transiting Exoplanet Survey Satellite and refine the system properties and ephemeris. We find a featureless transmission spectrum with increased absorption towards shorter wavelengths. We perform an atmospheric retrieval analysis with the AURA code, finding tentative evidence for haze in the upper atmospheric layers and a lower likelihood for a dense cloud deck. Whilst our retrieval analysis results point toward clouds and hazes, further evidence is needed to definitively reject a clear-sky scenario.

A classical approach to the restricted three-body problem is to analyze the dynamics of the massless body in the synodic reference frame. A different approach is represented by the perturbative treatment: in particular the averaged problem of a mean-motion resonance allows to investigate the long-term behavior of the solutions through a suitable approximation that focuses on a particular region of the phase space. In this paper, we intend to bridge a gap between the two approaches in the specific case of mean-motion resonant dynamics, establish the limit of validity of the averaged problem, and take advantage of its results in order to compute trajectories in the synodic reference frame. After the description of each approach, we develop a rigorous treatment of the averaging process, estimate the size of the transformation and prove that the averaged problem is a suitable approximation of the restricted three-body problem as long as the solutions are located outside the Hill's sphere of the secondary. In such a case, a rigorous theorem of stability over finite but large timescales can be proven. We establish that a solution of the averaged problem provides an accurate approximation of the trajectories on the synodic reference frame within a finite time that depend on the minimal distance to the Hill's sphere of the secondary. The last part of this work is devoted to the co-orbital motion (i.e., the dynamics in 1:1 mean-motion resonance) in the circular-planar case. In this case, an interpretation of the solutions of the averaged problem in the synodic reference frame is detailed and a method that allows to compute co-orbital trajectories is displayed.

The ability of a planet to maintain surface water, key to life as we know it, depends on solar and planetary energy. As a star ages, it delivers more energy to a planet. As a planet ages it produces less internal heat, which leads to cooling. For the Earth, interior cooling connects to plate tectonics - the surface manifestation of convection within the Earth's interior. This process cycles volatiles (CO2 and water) between surface and interior reservoirs, which affects planetary climate. Cycling rates depend on the efficiency of plate tectonic cooling. That efficiency remains debated and multiple hypotheses have been put forth. Geological proxy data allow us to validate these hypotheses accounting for model and data uncertainty. Multiple models pass the validation test. Those models define a distribution for terrestrial exoplanets akin to Earth, accounting for variations in tectonic efficiency. Feeding this distribution into climate models indicates that the time at which habitable conditions are established can vary by billions of years. Planets of the same absolute age and orbital distance can reside and not reside within the classic habitable zone due to differences in plate tectonic cooling efficiencies. The full model population allows a probability distribution to be constructed for the the time at which habitable conditions are established. The distribution indicates that Earth-like exoplanets, of the same age, can be at different evolutionary stages. It also indicates that planets around stars whose early evolution is unfavorable for life can become habitable later in their energetic histories.

Coronal loops, seen in solar coronal images, are believed to represent emission from magnetic flux tubes with compact cross-sections. We examine the 3D structure of plasma above an active region in a radiative magnetohydrodynamic simulation to locate volume counterparts for coronal loops. In many cases, a loop cannot be linked to an individual thin strand in the volume. While many thin loops are present in the synthetic images, the bright structures in the volume are fewer, and of complex shape. We demonstrate that this complexity can form impressions of thin bright loops, even in the absence of thin bright plasma strands. We demonstrate the difficulty of discerning from observations whether a particular loop corresponds to a strand in the volume, or a projection artifact. We demonstrate how apparently isolated loops could deceive observers, even when observations from multiple viewing angles are available. While we base our analysis on a simulation, the main findings are independent from a particular simulation setup and illustrate the intrinsic complexity involved in interpreting observations resulting from line-of-sight integration in an optically thin plasma. We propose alternative interpretation for strands seen in EUV images of the corona. The "coronal veil" hypothesis is mathematically more generic, and naturally explains properties of loops that are difficult to address otherwise -- such as their constant cross section and anomalously high density scale height. We challenge the paradigm of coronal loops as thin magnetic flux tubes, offering new understanding of solar corona and, by extension, of other magnetically confined bright, hot plasmas.

The recent star formation histories (SFHs) of post-starburst galaxies have been determined almost exclusively from detailed modeling of their composite star light. This has provided important but limited information on the number, strength, and duration of bursts of star formation. In this work, we present a direct and independent measure of the recent SFH of S12 (plate-mjd-fiber for SDSS 623-52051-207; designated EAS12 in Smercina et al. 2018) from its star cluster population. We detect clusters from high resolution, $UBR$ optical observations from HST, and compare their luminosities and colors with stellar population models to estimate the ages and masses of the clusters. No clusters younger than $\sim$70 Myr are found, indicating star formation shut off at this time. Clusters formed $\sim$120 Myr ago reach masses up to $\sim \mbox{few}\times10^7~M_{\odot}$, several times higher than similar age counterparts formed in actively merging galaxies like the Antennae and NGC 3256. We develop a new calibration based on known properties for 8 nearby galaxies to estimate the star formation rate (SFR) of a galaxy from the mass of the most massive cluster, $M_{\rm max}$. The cluster population indicates that S12 experienced an extremely intense but short-lived burst $\sim$120 Myr ago, with an estimated peak of $500^{+500}_{-250}~M_{\odot}~\mbox{yr}^{-1}$ and duration of $50\pm25$ Myr, one of the highest SFRs estimated for any galaxy in the nearby universe. Prior to the recent, intense burst, S12 was forming stars at a moderate rate of $\sim 3{-}5~M_{\odot}~\mbox{yr}^{-1}$, typical of spiral galaxies. However, the system also experienced an earlier burst approximately $1{-}3$ Gyr ago. While fairly uncertain, we estimate that the SFR during this earlier burst was $\sim20{-}30~M_{\odot}~\mbox{yr}^{-1}$, similar to the current SFR in the Antennae and NGC 3256.

Electromagnetic observations have provided strong evidence for the existence of massive black holes in the center of galaxies, but their origin is still poorly known. Different scenarios for the formation and evolution of massive black holes lead to different predictions for their properties and merger rates. LISA observations of coalescing massive black hole binaries could be used to reverse engineer the problem and shed light on these mechanisms. In this paper, we introduce a pipeline based on hierarchical Bayesian inference to infer the mixing fraction between different theoretical models by comparing them to LISA observations of massive black hole mergers. By testing this pipeline against simulated LISA data, we show that it allows us to accurately infer the properties of the massive black hole population as long as our theoretical models provide a reliable description of the Universe. We also show that measurement errors, including both instrumental noise and weak lensing errors, have little impact on the inference.

Gravastars are configurations of compact singularity-free gravitational objects which are interesting alternatives to classical solutions in the strong gravitational field regime. Although there are no static star-like solutions of the Einstein-Klein-Gordon equations for real scalar fields, we show that dynamical gravastars solutions arise through the direct interaction of a scalar field with matter. Two configurations presented here show that, within the internal zone, the scalar field plays a role similar to a cosmological constant, while the scalar field decays at large distances as the Yukawa potential. Like classical gravastars, these solutions exhibit small values of the temporal metric component near a transitional radial value, although this behaviour is not determined by the de Sitter nature of the internal space-time, but rather by a slowly-varying scalar field. The scalar field-matter interaction is able to define trapping forces that rigorously confine the polytropic gases to the interior of a sphere. At the surface of these spheres, pressures generated by the field-matter interaction play the role of "walls" preventing the matter from flowing out. These solutions predict a stronger scattering of the accreting matter with respect to Schwarzschild black holes.

In this work, we consider a recent novel regular black hole solution, which we denote the Simpson-Visser black-bounce model and investigate circular null geodesics to find the connection between the photon sphere, the event horizon and the black hole shadow radii. We also study the energy emission rate for this geometry and discuss how the parameters of the model affect the emission of particles around the black hole. Furthermore, we compare the resulting shadow of this regular black hole with observational data of the Event Horizon Telescope and find the allowed regions of the model parameters for which the obtained shadow is consistent with the data. Finally, we employ the correspondence between the quasinormal modes in the eikonal limit and shadow radius to study the scalar field perturbations in this background.

We study the effect of the cubic torus topology of the Universe on scalar cosmological perturbations which define the gravitational potential. We obtain three alternative forms of the solution for both the gravitational potential produced by point-like masses, and the corresponding force. The first solution includes the expansion of delta-functions into Fourier series, exploiting periodic boundary conditions. The second one is composed of summed solutions of the Helmholtz equation for the original mass and its images. Each of these summed solutions is the Yukawa potential. In the third formula, we express the Yukawa potentials via Ewald sums. We show that for the present Universe, both the bare summation of Yukawa potentials and the Yukawa-Ewald sums require smaller numbers of terms to yield the numerical values of the potential and the force up to desired accuracy. Nevertheless, the Yukawa formula is yet preferable owing to its much simpler structure.

In the solar wind plasma an excess of kinetic temperature along the background magnetic field stimulates proton firehose modes to grow if the parallel plasma beta parameter is sufficiently high, i.e., $\beta_{p \parallel}\gtrsim 1$. This instability can prevent the expansion-driven anisotropy from increasing indefinitely, and explain the observations. Moreover, such kinetic instabilities are expected to be even more effective in the presence of suprathermal Kappa-distributed populations, which are ubiquitous in the solar wind, are less affected by collisions than the core population, but contribute with an additional free energy. In this work we use both linear and extended quasi-linear (QL) frameworks to characterize the unstable periodic proton firehose modes (propagating parallel to the magnetic field) under the influence of suprathermal protons. Linear theory predicts a systematic stimulation of the instability, suprathermals amplifying the growth rates and decreasing the instability thresholds to lower anisotropies and lower plasma betas ($\beta_{p \parallel}<1$). In perfect agreement with these results, the QL approach reveals a significant enhancement of the resulting electromagnetic fluctuations up to the saturation with a stronger back reaction on protons, leading also to a faster and more efficient relaxation of the temperature anisotropy.

Motivated by the growing evidence for lepton flavour universality violation after the first results from Fermilab's muon $(g-2)$ measurement, we revisit one of the most widely studied anomaly free extensions of the standard model namely, gauged $L_{\mu}-L_{\tau}$ model, to find a common explanation for muon $(g-2)$ as well as baryon asymmetry of the universe via leptogenesis. The minimal setup allows TeV scale resonant leptogenesis satisfying light neutrino data while the existence of light $L_{\mu}-L_{\tau}$ gauge boson affects the scale of leptogenesis as the right handed neutrinos are charged under it. For $L_{\mu}-L_{\tau}$ gauge boson mass at GeV scale or above, the muon $(g-2)$ favoured parameter space is already ruled out by other experimental data while bringing down its mass to sub-GeV regime leads to vanishing lepton asymmetry due to highly restrictive structures of lepton mass matrices at the scale of leptogenesis. Extending the minimal model with two additional Higgs doublets can lead to a scenario consistent with successful resonant leptogenesis and muon $(g-2)$ while satisfying all relevant experimental data.

The deformability of a compact object under the presence of a tidal perturbation is encoded in the tidal Love numbers (TLNs), which vanish for isolated black holes in vacuum. We show that the TLNs of black holes surrounded by matter fields do not vanish and can be used to probe the environment around binary black holes. In particular, we compute the TLNs for the case of a black hole surrounded by a scalar condensate under the presence of scalar and vector tidal perturbations, finding a strong power-law behavior of the TLN in terms of the mass of the scalar field. Using this result as a proxy for gravitational tidal perturbations, we show that future gravitational-wave detectors like the Einstein Telescope and LISA can impose stringent constraints on the mass of ultralight bosons that condensate around black holes due to accretion or superradiance. Interestingly, LISA could measure the tidal deformability of dressed black holes across the range from stellar-mass ($\approx 10^2 M_\odot$) to supermassive ($\approx 10^7 M_\odot$) objects, providing a measurement of the mass of ultralight bosons in the range $(10^{-17} - 10^{-13}) \, {\rm eV}$ with less than $10\%$ accuracy, thus filling the gap between other superradiance-driven constraints coming from terrestrial and space interferometers. Altogether, LISA and Einstein Telescope can probe tidal effects from dressed black holes in the combined mass range $(10^{-17} - 10^{-11}) \, {\rm eV}$.

In the present investigation compact stellar models are dealt with in the framework of the modified gravity theory, specifically of $f(\mathbb{T},\mathcal{T})$ type. We have considered that the compact objects are following a spherically symmetric static metric and obtained the Einstein field equations in the spacetime of $f(\mathbb{T},\mathcal{T})$. To make the Einstein equations solvable we employ the methodology of conformal Killing vectors. Thereafter by using the MIT bag equation of state to the compact stars, considering that the stars are formed by strange quark, we find the solutions set. The solutions are examined via several physical tastings which exhibit viability of the model.

We study the optical appearance of a thin accretion disk around compact objects within the Einstein-Gauss-Bonnet gravity. Considering static spherically symmetric black holes and naked singularities we search for characteristic signatures which can arise in the observable images due to the modification of general relativity. While the images of the Gauss-Bonnet black holes closely resemble the Schwarzschild black hole, naked singularities possess a distinctive feature. A series of bright rings are formed in the central part of the images with observable radiation $10^3$ times larger than the rest of the flux making them observationally significant. We elucidate the physical mechanism, which causes the appearance of the central rings, showing that the image is determined by the light ring structure of the spacetime. In a certain region of the parametric space the Gauss-Bonnet naked singularities possess a stable and an unstable light ring. In addition the gravitational field becomes repulsive in a certain neighbourhood of the singularity. This combination of features leads to the formation of the central rings implying that the effect is not specific for the Einstein-Gauss-Bonnet gravity but would also appear for any other compact object with the same characteristics of the photon dynamics.

Chemical modelling serves two purposes in dynamical models: accounting for the effect of microphysics on the dynamics and providing observable signatures. Ideally, the former must be done as part of the hydrodynamic simulation but this comes with a prohibitive computational cost which leads to many simplifications being used in practice. To produce a statistical emulator that replicates a full chemical model capable of solving the temperature and abundances of a gas through time. This emulator should suffer only a minor loss of accuracy over including a full chemical solver in a dynamical model but would have a fraction of the computational cost. The gas-grain chemical code UCLCHEM was updated to include heating and cooling processes and a large dataset of model outputs from possible starting conditions was produced. A neural network was then trained to map directly from inputs to outputs. Chemulator replicates the outputs of UCLCHEM with an overall mean squared error (MSE) of 0.0002 for a single time step of 1000 yr and is shown to be stable over 1000 iterations with an MSE of 0.003 the log scaled temperature after one time step and 0.006 after 1000 time steps. Chemulator was found to be approximately 50,000 times faster than the time dependent model it emulates but can introduce a significant error to some models.

In this work we explore the effects that a possible primordial magnetic field can have on the inflaton effective potential, taking as the underlying model a warm inflation scenario, based on global supersymmetry with a new-inflation-type potential. The decay scheme for the inflaton field is a two-step process of radiation production, where the inflaton couples to heavy intermediate superfields, which in turn interact with light particles. In this context, we consider that both sectors, heavy and light, are charged and work in the strong magnetic field approximation for the light fields. We find an analytical expression for the one-loop effective potential, for an arbitrary magnetic field strength, and show that the trend of the magnetic contribution is to make the potential flatter, preserving the conditions for a successful inflationary process.