Papers for Thursday, Mar 25 2021

In this paper we investigate the effectiveness of direct statistical simulation (DSS) for two low-order models of dynamo action. The first model, which is a simple model of solar and stellar dynamo action, is third-order and has cubic nonlinearities whilst the second has only quadratic nonlinearities and describes the interaction of convection and an aperiodically reversing magnetic field. We show how DSS can be utilised to solve for the statistics of these systems of equations both in the presence and the absence of stochastic terms, by truncating the cumulant hierarchy at either second or third order. We compare two different techniques for solving for the statistics, timestepping -- which is able to locate only stable solutions of the equations for the statistics and direct detection of the fixed points. We develop a complete methodology and symbolic package in Python for deriving the statistical equations governing the Low-order dynamic systems in cumulant expansions. We demonstrate that although direct detection of the fixed points is efficient and accurate for DSS truncated at second order, the addition of higher order terms leads to the inclusion of many unstable fixed points that may be found by direct detection of the fixed point by iterative methods. In those cases timestepping is a more robust protocol for finding meaninful solutions to DSS.

Aims: We aim to investigate the consequences of a fast massive stellar remnant -- a black hole (BH) or a neutron star (NS) -- encountering a planetary system. Methods: We model a~close encounter between the actual Solar system (SS) and a $2\,M_\odot$ NS and a $10\,M_\odot$ BH, using a few-body symplectic integrator. We use a range of impact parameters, orbital phases at the start of the simulation derived from the current SS orbital parameters, encounter velocities, and incidence angles relative to the plane of the SS. Results: We give the distribution of possible outcomes, i.e., when the SS remains bound, when it suffers a partial or complete disruption, and in which cases the intruder is able to capture one or more planets, yielding planetary systems around a BH or a NS. We also show examples of the long-term stability of the captured planetary systems.

It is widely assumed that long-wavelength infrared (IR) emission from cold dust (T~20-40 K) is a reliable tracer of star formation even in the presence of a bright active galactic nucleus (AGN). Based on radiative transfer (RT) models of clumpy AGN tori, hot dust emission from the torus contributes negligibly to the galaxy spectral energy distribution (SED) at $\lambda>60\,\mu$m. However, these models do not include AGN heating of host-galaxy-scale diffuse dust, which may have far-IR (FIR) colors comparable to cold diffuse dust heated by stars. To quantify the contribution of AGN heating to host-galaxy-scale cold dust emission, we perform dust RT calculations on a simulated galaxy merger both including and excluding the bright AGN that it hosts. By differencing the SEDs yielded by RT calculations with and without AGN that are otherwise identical, we quantify the FIR cold dust emission arising solely from re-processed AGN photons. In extreme cases, AGN-heated host-galaxy-scale dust can increase galaxy-integrated FIR flux densities by factors of 2-4; star formation rates calculated from the FIR luminosity assuming no AGN contribution can overestimate the true value by comparable factors. Because the FIR colors of such systems are similar to those of purely star-forming galaxies and redder than torus models, broadband SED decomposition may be insufficient for disentangling the contributions of stars and heavily dust-enshrouded AGN in the most IR-luminous galaxies. We demonstrate how kpc-scale resolved observations can be used to identify deeply dust-enshrouded AGN with cool FIR colors when spectroscopic and/or X-ray detection methods are unavailable.

Long secondary periods (LSPs), observed in a third of pulsating red giant stars, are the only unexplained type of large-amplitude stellar variability known at this time. Here we show that this phenomenon is a manifestation of a substellar or stellar companion orbiting the red giant star. Our investigation is based on a sample of about 16,000 well-defined LSP variables detected in the long-term OGLE photometric database of the Milky Way and Magellanic Clouds, combined with the mid-infrared data extracted from the NEOWISE-R archive. From this collection, we selected about 700 objects with stable, large-amplitude, well-sampled infrared light curves and found that about half of them exhibit secondary eclipses, thus presenting an important piece of evidence that the physical mechanism responsible for LSPs is binarity. Namely, the LSP light changes are due to the presence of a dusty cloud orbiting the red giant together with the companion and obscuring the star once per orbit. The secondary eclipses, visible only in the infrared wavelength, occur when the cloud is hidden behind the giant. In this scenario, the low-mass companion is a former planet that has accreted a significant amount of mass from the envelope of its host star and grown into a brown dwarf.

To date, at least three comets -- 2I/Borisov, C/2016 R2 (PanSTARRS), and C/2009 P1 (Garradd) -- have been observed to have unusually high CO concentrations compared to water. We attempt to explain these observations by modeling the effect of drifting solid (ice and dust) material on the ice compositions in protoplanetary disks. We find that, independent of the exact disk model parameters, we always obtain a region of enhanced ice-phase CO/H2O that spreads out in radius over time. The inner edge of this feature coincides with the CO snowline. Almost every model achieves at least CO/H2O of unity, and one model reaches a CO/H2O ratio > 10. After running our simulations for 1 Myr, an average of 40% of the disk ice mass contains more CO than H2O ice. In light of this, a population of CO ice enhanced planetesimals are likely to generally form in the outer regions of disks, and we speculate that the aforementioned CO-rich comets may be more common, both in our own Solar System and in extrasolar systems, than previously expected.

We explore the connection between the kinematics, structures and stellar populations of massive galaxies at $0.6<z<1.0$ using the Fundamental Plane (FP). Combining stellar kinematic data from the Large Early Galaxy Astrophysics Census (LEGA-C) survey with structural parameters measured from deep Hubble Space Telescope imaging, we obtain a sample of 1419 massive ($\log(M_*/M_\odot) >10.5$) galaxies that span a wide range in morphology, star formation activity and environment, and therefore is representative of the massive galaxy population at $z\sim0.8$. We find that quiescent and star-forming galaxies occupy the parameter space of the $g$-band FP differently and thus have different distributions in the dynamical mass-to-light ratio ($M_{\rm dyn}/L_g$), largely owing to differences in the stellar age and recent star formation history, and, to a lesser extent, the effects of dust attenuation. In contrast, we show that both star-forming and quiescent galaxies lie on the same mass FP at $z\sim 0.8$, with a comparable level of intrinsic scatter about the plane. We examine the variation in $M_{\rm dyn}/M_*$ through the thickness of the mass FP, finding no significant residual correlations with stellar population properties, S\'ersic index, or galaxy overdensity. Our results suggest that, at fixed size and velocity dispersion, the variations in $M_{\rm dyn}/L_g$ of massive galaxies reflect an approximately equal contribution of variations in $M_*/L_g$, and variations in the dark matter fraction or initial mass function.

The successive discoveries of binary merger events by Advanced LIGO-Virgo have been revealing the statistical properties of binary black hole (BBH) populations. A stochastic gravitational wave background (GWB) is a useful tool to probe the cosmological evolution of those compact mergers. In this paper, we study the upper bound on a GWB produced by BBH mergers, whose stellar progenitors dominate the reionization process at the cosmic dawn. Since early reionization by those progenitors yields a high optical depth of the universe inconsistent with the {\it Planck} measurements, the cumulative mass density is limited to $\rho_\star \lesssim 10^7~M_\odot~{\rm Mpc}^{-3}$. Even with this upper bound, the amplitude of a GWB owing to the high-$z$ BBH mergers is expected to be as high as $\Omega_{\rm gw}\simeq 1.48_{-1.27}^{+1.80}\times 10^{-9}$ at $f\simeq 25$ Hz, while their merger rate at the present-day is consistent or lower than the observed GW event rate. This level of GWB is detectable at the design sensitivity of Advanced LIGO-Virgo and would indicate a major contribution of the high-$z$ BBH population to the local GW events. The spectral index is expected to be substantially flatter than the canonical value of $\simeq 2/3$ generically produced by lower-redshift and less massive BBHs. Moreover, if their mass function is more top-heavy than in the local universe, the GWB spectrum is even more skewed toward lower frequencies, which would allow us to extract information on the mass function of merging BBHs at high redshifts.

We present observations of two bright M dwarfs (TOI-1634 and TOI-1685: $J=9.5-9.6$) hosting ultra-short period (USP) planet candidates, identified by the TESS mission. The two stars are similar in temperature, mass, and radius, and the planet candidates are both super-Earth-sized ($1.25\,R_\oplus<R_p<2.0\,R_\oplus$). For both systems, light curves from the ground-based photometry exhibit planetary transits, whose depths are consistent with those by the TESS photometry. We also refine the transit ephemerides based on the ground-based photometry, finding the orbital periods of $P=0.9893457 \pm 0.0000034$ day and $P=0.6691425 \pm 0.0000037$ day for TOI-1634b and TOI-1685b, respectively. Through intensive radial velocity (RV) observations using the IRD spectrograph on the Subaru 8.2m telescope, we confirm the planetary nature of the TOIs, and measure their masses: $9.94\pm0.94\,M_\oplus$ and $3.84\pm0.96\,M_\oplus$ for TOI-1634b and TOI-1685b, respectively, when the observed RVs are fitted with a single-planet circular-orbit model. Combining those with the planet radii of $R_p=1.773\pm 0.077\,R_\oplus$ (TOI-1634b) and $1.463\pm0.065\,R_\oplus$ (TOI-1685b), we find that both USP planets have mean densities consistent with an Earth-like internal composition, which is typical for small USP planets. TOI-1634b is currently the most massive USP planet in this category, and it resides near the radius valley, which makes it a benchmark planet in the context of testing the formation scenarios for USP planets. Excess scatter in the RV residuals for TOI-1685 suggests the presence of a possible secondary planet or unknown activity/instrumental noise in the RV data, but further observations are required to check those possibilities.

We report the detection of an Al II line at 2669.155 Angstroms in 11 metal-poor stars, using ultraviolet spectra obtained with the Space Telescope Imaging Spectrograph on board the Hubble Space Telescope. We derive Al abundances from this line using a standard abundance analysis, assuming local thermodynamic equilibrium (LTE). The mean [Al/Fe] ratio is -0.06 +/- 0.04 (sigma = 0.22) for these 11 stars spanning -3.9 < [Fe/H] < -1.3, or [Al/Fe] = -0.10 +/- 0.04 (sigma = 0.18) for 9 stars spanning -3.0 < [Fe/H] < -1.3 if two carbon-enhanced stars are excluded. We use these abundances to perform an empirical test of non-LTE (NLTE) abundance corrections predicted for resonance lines of Al I, including the commonly-used optical Al I line at 3961 Angstroms. The Al II line is formed in LTE, and the abundance derived from this line matches that derived from high-excitation Al I lines predicted to have minimal NLTE corrections. The differences between the abundance derived from the Al II line and the LTE abundance derived from Al I resonance lines are +0.4 to +0.9 dex, which match the predicted NLTE corrections for the Al I resonance lines. We conclude that the NLTE abundance calculations are approximately correct and should be applied to LTE abundances derived from Al I lines.

Spin-orbit coupling of planetary systems plays an important role in the dynamics and habitability of planets. However, symplectic integrators that can accurately simulate not only how orbit affects spin but also how spin affects orbit have not been constructed for general systems. Thus, we develop symplectic Lie-group integrators to simulate systems consisting gravitationally interacting rigid bodies. A user friendly package (GRIT) is provided and external forcings such as tidal interactions are also included. As a demonstration, this package is applied to Trappist-I. It shows that the differences in transit timing variations due to spin-orbit coupling could reach a few min in ten year measurements, and strong planetary perturbations can push Trappist-I f, g and h out of the synchronized states.

While dust is a major player in galaxy evolution, its relationship with gas and stellar radiation in the early universe is still not well understood. We combine 3D-HST emission line fluxes with far-UV through far-IR photometry in a sample of 669 emission-line galaxies (ELGs) between 1.2 < z < 1.9 and use the MCSED spectral energy distribution fitting code to constrain the galaxies' physical parameters, such as their star formation rates (SFRs), stellar masses, and dust masses. We find that the assumption of energy balance between dust attenuation and emission is likely unreasonable in many cases. We highlight a relationship between the mass-specific star formation rate (sSFR), stellar mass, and dust mass, although its exact form is still unclear. Finally, a stacking of H$\alpha$ and H$\beta$ fluxes shows that nebular attenuation increases with stellar mass and SFR for IR-bright ELGs.

We investigate the properties, composition, and dynamics of dust formation and growth for a diverse set of core-collapse supernovae (CCSNe), varying the progenitor mass, explosion energy, and engine type. These explosions are evolved with a 1-D Lagrangian hydrodynamics code out to several hundred days to model the ejecta as it expands and cools. A multigrain dust nucleation and growth model is applied to these results. We find that higher explosion energies lead to an earlier onset of dust formation, smaller grain sizes, and larger silicate abundances. Further, we see that nuclear burning during the explosion leads to enhanced formation of silicate dust. Finally, we build composite models from our suite to predict the efficiency of CCSNe dust production as a function of metallicity.

We present the discovery of eclipses in the X-ray light curves of the X-ray binary Swift J1858.6-0814. From these, we find an orbital period of $P=76841.3_{-1.4}^{+1.3}$ s ($\approx21.3$ hours) and an eclipse duration of $t_{\rm ec}=4098_{-18}^{+17}$ s ($\approx1.14$ hours). We also find several absorption dips during the pre-eclipse phase. From the eclipse duration to orbital period ratio, the inclination of the binary orbit is constrained to $i>70^\circ$. The most likely range for the companion mass suggests that the inclination is likely to be closer to this value than $90^\circ$. The eclipses are also consistent with earlier data, in which strong variability ('flares') and the long orbital period prevent clear detection of the period or eclipses. We also find that the bright flares occurred preferentially in the post-eclipse phase of the orbit, likely due to increased thickness at the disc-accretion stream interface preventing flares being visible during the pre-eclipse phase. This supports the notion that variable obscuration is responsible for the unusually strong variability in Swift J1858.6-0814.

We investigate the possibility that radio-bright active galactic nuclei (AGN) are responsible for the TeV--PeV neutrinos detected by IceCube. We use an unbinned maximum-likelihood-ratio method, 10 years of IceCube muon-track data, and 3388 radio-bright AGN selected from the Radio Fundamental Catalog. None of the AGN in the catalog have a large global significance. The two most significant sources have global significance of $\simeq$ 1.5$\sigma$ and 0.8$\sigma$, though 4.1$\sigma$ and 3.8$\sigma$ local significance. Our stacking analyses show no significant correlation between the whole catalog and IceCube neutrinos. We infer from the null search that this catalog can account for at most 30\% (95\% CL) of the diffuse astrophysical neutrino flux measured by IceCube. Moreover, our results disagree with recent work that claimed a 4.1$\sigma$ detection of neutrinos from the sources in this catalog, and we discuss the reasons of the difference.

It has long been suspected that powerful radio sources may lower the efficiency with which stars form from the molecular gas in their host galaxy, but so far, alternative mechanisms, in particular related to the stellar mass distribution in the massive bulges of their host galaxies, are not ruled out. We present new ALMA CO(1-0) interferometry of cold molecular gas in the nearby (z=0.0755), massive (M_stellar=4x10^11 M_sun), isolated, late-type spiral galaxy 2MASSX J23453269-044925, which is outstanding for having two pairs of powerful, giant radio jets, and a bright X-ray halo. The molecular gas is in a massive (M_gas=2x10^10 M_sun), 24 kpc wide, rapidly rotating ring, which is associated with the inner stellar disk. Broad (FWHM=70-180 km s^-1) lines with complex profiles associated with the radio source are seen over large regions in the ring, indicating gas velocities that are high enough to keep the otherwise marginally Toomre-stable gas from fragmenting into gravitationally bound, star-forming clouds. About 1-2% of the jet kinetic energy are required to power these motions. Resolved star-formation rate surface densities fall factors 50-75 short of expectations from the Kennicutt-Schmidt law of star-forming galaxies, and near gas-rich early-type galaxies with signatures of star formation lowered by jet feedback. We argue that radio AGN feedback is the only plausible mechanism to explain the low star-formation rates in this galaxy. Previous authors have already noted that the X-ray halo of J2345-0449 implies a baryon fraction near the cosmic average, which is very high for a galaxy. We contrast this finding with other, equally massive baryon-rich spiral galaxies without prominent radio sources. Most of the baryons in these galaxies are in stars, not in the halos. We also discuss the implications of our results for our general understanding of AGN feedback in massive galaxies.

We report the discovery of a circumstellar debris disk viewed nearly edge-on and associated with the young, K1 star BD+45$^{\circ}$598 using high-contrast imaging at 2.2$\mu$m obtained at the W.M.~Keck Observatory. We detect the disk in scattered light with a peak significance of $\sim$5$\sigma$ over three epochs, and our best-fit model of the disk is an almost edge-on $\sim$70 AU ring, with inclination angle $\sim$87$^\circ$. Using the NOEMA interferometer at the Plateau de Bure Observatory operating at 1.3mm, we find resolved continuum emission aligned with the ring structure seen in the 2.2$\mu$m images. We estimate a fractional infrared luminosity of $L_{IR}/L_{tot}$ $\simeq6^{+2}_{-1}$$\times$$10^{-4}$, higher than that of the debris disk around AU Mic. Several characteristics of BD+45$^{\circ}$598, such as its galactic space motion, placement in a color-magnitude diagram, and strong presence of Lithium, are all consistent with its membership in the $\beta$ Pictoris Moving Group with an age of 23$\pm$3 Myr. However, the galactic position for BD+45$^{\circ}$598 is slightly discrepant from previously-known members of the $\beta$ Pictoris Moving Group, possibly indicating an extension of members of this moving group to distances of at least 70pc. BD+45$^{\circ}$598 appears to be an example from a population of young circumstellar debris systems associated with newly identified members of young moving groups that can be imaged in scattered light, key objects for mapping out the early evolution of planetary systems from $\sim$10-100 Myr. This target will also be ideal for northern-hemisphere, high-contrast imaging platforms to search for self-luminous, planetary mass companions residing in this system.

In this study we follow up our recent paper (Monteiro et al. 2020) and present a homogeneous sample of fundamental parameters of open clusters in our Galaxy, entirely based on Gaia DR2 data. We used published membership probability of the stars derived from Gaia DR2 data and applied our isochrone fitting code, updated as in Monteiro et al. (2020), to GB and GR Gaia DR2 data for member stars. In doing this we take into account the nominal errors in the data and derive distance, age, and extinction of each cluster. This work therefore provides parameters for 1743 open clusters and, as a byproduct, a list of likely not physical or dubious open clusters is provided as well for future investigations. Furthermore, it was possible to estimate the mean radial velocity of 831 clusters (198 of which are new and unpublished so far) using stellar radial velocities from Gaia DR2 catalog. By comparing the open cluster distances obtained from isochrone fitting with those obtained from a maximum likelihood estimate of individual member parallaxes, we found a systematic offset of $(-0.05\pm0.04)$mas.

High resolution precision spectroscopy provides a multitude of robust techniques for probing exoplanetary atmospheres. We present multiple VLT/ESPRESSO transit observations of the hot-Jupiter exoplanet WASP-19b with previously published but disputed atmospheric features from low resolution studies. Through spectral synthesis and modeling of the Rossiter-McLaughlin (RM) effect we calculate stellar, orbital and physical parameters for the system. From narrow-band spectroscopy we do not detect any of H\,I, Fe\,I, Mg\,I, Ca\,I, Na\,I and K\,I neutral species, placing upper limits on their line contrasts. Through cross correlation analyses with atmospheric models, we do not detect Fe\,I and place a 3$\sigma$ upper limit of $\log\,(X_{\textrm{Fe}}/X_\odot) \approx -1.83\,\pm\,0.11$ on its mass fraction, from injection and retrieval. We show the inability to detect the presence of H$_2$O for known abundances, owing to lack of strong absorption bands, as well as relatively low S/N ratio. We detect a barely significant peak (3.02\,$\pm$\,0.15\,$\sigma$) in the cross correlation map for TiO, consistent with the sub-solar abundance previously reported. This is merely a hint for the presence of TiO and does \textit{not} constitute a confirmation. However, we do confirm the presence of previously observed enhanced scattering towards blue wavelengths, through chromatic RM measurements, pointing to a hazy atmosphere. We finally present a reanalysis of low resolution transmission spectra of this exoplanet, concluding that unocculted starspots alone cannot explain previously detected features. Our reanalysis of the FORS2 spectra of WASP-19b finds a $\sim$\,100$\times$ sub-solar TiO abundance, precisely constrained to $\log\,X_{\textrm{TiO}} \approx -7.52 \pm 0.38$, consistent with the TiO hint from ESPRESSO. We present plausible paths to reconciliation with other seemingly contradicting results.

The eccentricity of a planet's orbit and the inclination of its orbital plane encode important information about its formation and history. However, exoplanets detected via direct-imaging are often only observed over a very small fraction of their period, making it challenging to perform reliable physical inferences given wide, unconstrained posteriors. The aim of this project is to investigate biases (deviation of the median and mode of the posterior from the true values of orbital parameters, and the width and coverage of their credible intervals) in the estimation of orbital parameters of directly-imaged exoplanets, particularly their eccentricities, and to define general guidelines to perform better estimations of uncertainty. For this, we constructed various orbits and generated mock data for each spanning $\sim 0.5 \%$ of the orbital period. We used the Orbits For The Impatient (OFTI) algorithm to compute orbit posteriors, and compared those to the true values of the orbital parameters. We found that the inclination of the orbital plane is the parameter that most affects our estimations of eccentricity, with orbits that appear near edge-on producing eccentricity distributions skewed away from the true values, and often bi-modal. We also identified a degeneracy between eccentricity and inclination that makes it difficult to distinguish posteriors of face-on, eccentric orbits and edge-on, circular orbits. For the exoplanet-imaging community, we propose practical recommendations, guidelines and warnings relevant to orbit-fitting.

We present results from high signal-to-noise optical spectropolarimetric observations of the Seyfert 1 galaxies NGC783 and Mrk509 in the 3500-7000 A range. We find complex structure in the polarized emission for both objects. In particular, Position Angle (PA) changes across the Balmer lines show a distinctive 'M'-shaped profile that had not been observed in this detail before, but could represent a common trait in Seyfert 1 galaxies. In fact, while this shape is observed in all Balmer lines in NGC3783, Mrk509 transitions into a 'M'-shaped PA profile for higher transitions lines. We have modeled the observed profiles using the STOKES radiative transfer code and assuming that the scattering region is co-spatial with the BLR and outflowing. The results give compelling new evidence for the presence of nuclear winds in these two Seyfert 1 galaxies.

The precise astrometric measurements of the Gaia Data Release 2 have opened the door to detailed tests of the predictions of white dwarf cooling models. Significant discrepancies between theory and observations have been identified, the most striking affecting ultramassive white dwarfs. Cheng et al. (2019) found that a small fraction of white dwarfs on the so-called Q branch must experience an extra cooling delay of $\sim 8\,$Gyr not predicted by current models. $^{22}$Ne phase separation in a crystallizing C/O white dwarf can lead to a distillation process that efficiently transports $^{22}$Ne toward its center, thereby releasing a considerable amount of gravitational energy. Using state-of-the-art Monte Carlo simulations, we show that this mechanism can largely resolve the ultramassive cooling anomaly if the delayed population consists of white dwarfs with moderately above-average $^{22}$Ne abundances. We also argue that $^{22}$Ne phase separation can account for the smaller cooling delay currently missing for models of white dwarfs with more standard compositions.

We performed a radial velocity (RV) monitoring of the 40 Myr old star DS Tuc A with HARPS at the ESO-3.6m to determine the planetary mass of its 8.14-days planet, first revealed by TESS. We also observed two planetary transits with HARPS and ESPRESSO at ESO-VLT, to measure the Rossiter-McLaughlin (RM) effect and characterise the planetary atmosphere. We measured the high-energy emission of the host with XMM observations to investigate models for atmospheric evaporation. We employed Gaussian Processes (GP) regression to model the high level of the stellar activity, which is more than 40 times larger than the expected RV planetary signal. We extracted the transmission spectrum of DS Tuc A b from the ESPRESSO data and searched for atmospheric elements/molecules either by single-line retrieval and by performing cross-correlation with a set of theoretical templates. Through a set of simulations, we evaluated different scenarios for the atmospheric photo-evaporation of the planet induced by the strong XUV stellar irradiation. While the stellar activity prevented us from obtaining a clear detection of the planetary signal from the RVs, we set a robust mass upper limit of 14.4 M_e for DS Tuc A b. We also confirm that the planetary system is almost (but not perfectly) aligned. The strong level of stellar activity hampers the detection of any atmospheric compounds, in line with other studies presented in the literature. The expected evolution of DS Tuc A b from our grid of models indicates that the planetary radius after the photo-evaporation phase will fall within the Fulton gap. The comparison of the available parameters of known young transiting planets with the distribution of their mature counterpart confirms that the former are characterised by a low density, with DS Tuc A b being one of the less dense.

PolarLight is a space-borne X-ray polarimeter that measures the X-ray polarization via electron tracking in an ionization chamber. It is a collimated instrument and thus suffers from the background on the whole detector plane. The majority of background events are induced by high energy charged particles and show ionization morphologies distinct from those produced by X-rays of interest. Comparing on-source and off-source observations, we find that the two datasets display different distributions on image properties. The boundaries between the source and background distributions are obtained and can be used for background discrimination. Such a means can remove over 70% of the background events measured with PolarLight. This approaches the theoretical upper limit of the background fraction that is removable and justifies its effectiveness. For observations with the Crab nebula, the background contamination decreases from 25% to 8% after discrimination, indicative of a polarimetric sensitivity of around 0.2 Crab for PolarLight. This work also provides insights into future X-ray polarimetric telescopes.

In order to gain possible hints for planet formation from the current data of known extra-solar planets, the period-ratios and mass-ratios of adjacent planet pairs in multi-planet systems are determined. A moderate period-ratio-mass-ratio correlation is found to have a correlation coefficient r=0.5779 with 99% confidence interval (0.464, 0.672). In contrast, for non-adjacent planet pairs, the correlation coefficient is r=0.2820 with 99% confidence interval (0.133, 0.419). Our results reveal the imprint of planet-planet interactions of the adjacent planet pairs in a certain fraction of the multi-planet systems during the stage of planet formation.

We presented a multi-wavelength study of AFGL 333-Ridge. The molecular line data reveals that the AFGL 333-Ridge has two independent velocity components at -50.5 km/s and -48.0 km/s. In the Position-Velocity diagram, the bridge feature connects with two parts that are spatially correlated but separated in velocity. These observational evidences support the scenario that the two velocity components have collided and merged into one molecular cloud. The majority of Class I YSOs are distributed within the collision region, suggesting that the cloud-cloud collision has induced the YSOs formation in the ridge. Using the radio recombination line (RRL) data obtained by the Five-hundred-meter Aperture Spherical radio Telescope (FAST), the RRL velocities of three HII regions are consistent with that of the AFGL 333-Ridge. By comparing the three HII regions' dynamical ages with the collision timescale of the two components, we conclude that the influence of the three HII regions may not drive the two clouds to merge. The formation of the AFGL 333-Ridge is probably due to the expansion of the giant HII region W4.

We present one of the first Shanghai Tian Ma Radio Telescope (TMRT) K Band observations towards a sample of 26 infrared dark clouds (IRDCs). We observed the (1,1), (2,2), (3,3), and (4,4) transitions of NH$_{3}$ together with CCS (2$_{1}$-1$_{0}$) and HC$_{3}$N $J\,$=2-1, simultaneously. The survey dramatically increases the existing CCS-detected IRDC sample from 8 to 23, enabling a better statistical study of the ratios of carbon-chain molecules (CCM) to N-bearing molecules in IRDCs. With the newly developed hyperfine group ratio (HFGR) method of fitting NH$_{3}$ inversion lines, we found the gas temperature to be between 10 and 18 K. The column density ratios of CCS to NH$_{3}$ for most of the IRDCs are less than 10$^{-2}$, distinguishing IRDCs from low-mass star-forming regions. We carried out chemical evolution simulations based on a three-phase chemical model NAUTILUS. Our measurements of the column density ratios between CCM and NH$_{3}$ are consistent with chemical evolutionary ages of $\lesssim$10$^{5}$ yr in the models. Comparisons of the data and chemical models suggest that CCS, HC$_{3}$N, and NH$_{3}$ are sensitive to the chemical evolutionary stages of the sources.

Mars today has no active volcanism and its atmosphere is oxidizing, dominated by the photochemistry of CO2 and H2O. Using a one-dimensional photochemical model, we consider whether plausible volcanic gas fluxes could have switched the redox-state of the past martian atmosphere to reducing conditions. In our model, the total quantity and proportions of volcanic gases depend on the water content, outgassing pressure, and oxygen fugacity of the source melt. We find that with reasonable melt parameters the past martian atmosphere (~3.5 Gyr to present) could have easily reached reducing and anoxic conditions with modest levels of volcanism, >0.14 km^3/yr, well within the range of prior estimates. Counter-intuitively we also find that more reducing melts with lower oxygen fugacity require greater amounts of volcanism to switch a paleo-atmosphere from oxidizing to reducing. The reason is that sulfur is more stable in such melts and lower absolute fluxes of sulfur-bearing gases more than compensate for increases in the proportions of H2 and CO. These results imply that ancient Mars should have experienced periods with anoxic and reducing atmospheres even through the mid-Amazonian whenever volcanic outgassing was sustained at sufficient levels. Reducing anoxic conditions are potentially conducive to the synthesis of prebiotic organic compounds, such as amino acids, and are therefore relevant to the possibility of life on Mars. Also, anoxic reducing conditions should have influenced the type of minerals that were formed on the surface or deposited from the atmosphere such as elemental polysulfur (S8) as a signature of past reducing atmospheres. Finally, our models allow us to estimate the amount of volcanically sourced atmospheric sulfate deposited over Mars' history, approximately 10^6 to 10^9 Tmol, with a spread depending on assumed outgassing rate history and magmatic source conditions.

The LIGO-Virgo gravitational-wave (GW) observation unveiled the new population of black holes (BHs) that appears to have an extended mass spectrum up to around $70M_\odot$, much heavier than the previously-believed mass range ($\sim 8M_\odot$). In this paper, we study the capability of a microlensing observation of stars in the Milky Way (MW) bulge region to identify BHs of GW mass scales, taking into account the microlensing parallax characterized by the parameter $\pi_{\rm E}\propto M^{-1/2}$ ($M$ is the mass of a lens), which is a dimension-less quantity defined by the ratio of the astronomical unit to the projected Einstein radius. First, assuming that BHs follow the same spatial and velocity distributions of stars as predicted by the standard MW model, we show that microlensing events with long light curve timescales, $t_{\rm E}\gtrsim 100~{\rm days}$, and small parallax effects, $\pi_{\rm E}\sim 10^{-2}$, are dominated by BH lenses compared to stellar-mass lenses. Second, using a Markov chain Monte Carlo analysis of the simulated light curve, we show that BH lens candidates are securely identified on individual basis, if the parallax effect is detected or well constrained to the precision of a percent level in $\pi_{\rm E}$. We also discuss that a microlensing event of an intermediate-mass BH of $\sim 1000M_\odot$, if it occurs, can be identified in a distinguishable way from stellar-mass BHs.

The Central Molecular Zone (CMZ) is a ring-like structure sitting at the center of the Milky Way. Using the 870 $\mu$m continuum map from the APEX Telescope Large Area Survey of the Galaxy (ATLASGAL), we study the anisotropy of the density structure of the gas in the Central Molecular Zone (CMZ) using the 2D correlation function. To quantify the spatial anisotropy, we define the critical angle $\theta_{\rm half}$, as well as the anisotropy parameter $A\equiv \frac{\pi}{4\theta_{\rm half}}-1$. We find that the density structure is strongly anisotropic at the large scale (~ 100 pc), and the degree of spatial anisotropy decreases with the decreasing scale. At the scale of ~ 10 pc, the structure is still mildly anisotropic. In our analyses, we provide a quantitative description of the anisotropic density structure of gas in the CMZ, and the formalism can be applied to different regions to study their differences.

Heisenberg's uncertainty principle tells us that it is impossible to determine simultaneously the position of a photon crossing a telescope's aperture and its momentum. Super-resolution imaging techniques rely on modification of the observed sample, or on entangling photons. In astronomy we have no access to the object, but resolution may be improved by optical amplification. Unfortunately, spontaneous emission contributes noise and negates the possible gain from stimulated emissions. We show that it is possible to increase the weight of the stimulated photons by considering photon statistics, and observe an improvement in resolution. Most importantly, we demonstrate a method which can apply for all imaging purposes.

New surveys such as ESA's Euclid mission and NASA's Nancy Grace Roman Space Telescope are planned to map with unprecedented precision the large-scale structure of the Universe by measuring the 3D positions of tens of millions of galaxies. It is necessary to develop theoretically modelled galaxy catalogues to estimate the expected performance and to optimise the analysis strategy of these surveys. We populate two pairs of 1 (Gpc/h)$^3$ volume dark-matter-only simulations from the UNIT project with galaxies using the SAGE semi-analytic model of galaxy formation, coupled to the photoionisation model GET_EMLINES to estimate their Halpha emission. These catalogues represent a unique suite that includes galaxy formation physics and - thanks to the fixed-pair technique used - an effective volume of ~ (5 Gpc/h)$^3$, which is several times larger than the Euclid survey. We study the abundance and clustering of those model star-forming Halpha emission-line galaxies (ELGs). For scales greater than ~5 Mpc/h, we find for the ELGs a scale-independent bias with values in the range b $\in$ [1,4.5], increasing with redshift over the interval z $\in$ [0.5,2]. Model galaxy properties, including their emission-line fluxes are publicly available.

The evolution of the atmospheres of low and intermediate-mass planets is strongly connected to the physical properties of their host stars. The types and the past activities of planet-hosting stars can, therefore, affect the overall planetary population. In this paper, we perform a comparative study of sub-Neptune-like planets orbiting stars of different masses and different evolutionary histories. We discuss the general patterns of the evolved population as a function of parameters and environments of planets. As a model of the atmospheric evolution, we employ the own framework combining planetary evolution in MESA with the realistic prescription of the escape of hydrogen-dominated atmospheres. {We find that the final populations look qualitatively similar in terms of the atmospheres survival around different stars, but qualitatively different, with this difference accentuated for planets orbiting more massive stars. We show that a planet has larger chances of keeping its primordial atmosphere in the habitable zone of a solar mass star compared to M or K dwarfs and if it starts the evolution having a relatively compact envelope.} We also address the problem of the uncertain initial temperatures (luminosities) of planets and show that this issue is only of particular importance for planets exposed to extreme atmospheric mass-losses.

Supermassive black holes in the centre of galaxies dominate the gravitational potential of their surrounding stellar clusters. In these dense environments, stars follow nearly Keplerian orbits, which get slowly distorted as a result of the potential fluctuations generated by the stellar cluster itself as a whole. In particular, stars undergo a rapid relaxation of their eccentricities through both resonant and non-resonant processes. An efficient implementation of the resonant diffusion coefficients allows for detailed and systematic explorations of the parameter space describing the properties of the stellar cluster. In conjunction with recent observations of the S-cluster orbiting SgrA*, this framework can be used to jointly constrain the distribution of the unresolved, old, background stellar cluster and the characteristics of a putative dark cluster. Specifically, we show how this can be used to estimate the typical mass and cuspide exponent of intermediate-mass black holes consistent with the relaxed state of the distribution of eccentricities in the observed S-cluster. This should prove useful in constraining super massive black hole formation scenarios.

The transition region between the North Equatorial Band (NEBn) and North Tropical Zone (NTrZ) in Jupiter is home to convective storms, systems of cyclones and anticyclones and atmospheric waves. A large anticyclone formed in the year 2006 at planetographic latitude 19N and persists since then after a complex dynamic history, being possibly the third longest-lived oval in the planet after Jupiter's Great Red Spot and oval BA. This anticyclone has experienced close interactions with other ovals, merging with another oval in February 2013; it has also experienced color changes, from white to red (September 2013). The oval survived the effects of the closely located North Temperate Belt Disturbance, which occurred in October 2016 and fully covered the oval, rendering it unobservable for a short time. When it became visible again at its expected longitude from its previous longitudinal track, it reappeared as a white large oval keeping this color and the same morphology since 2017 at least until the onset of the new convective disturbance in Jupiter's North Temperate Belt in August 2020. Here we describe the historic evolution of the properties of this oval. We use JunoCam and Hubble Space Telescope (HST) images to measure its size and its internal rotation. We also used HST and PlanetCam-UPV/EHU multi-wavelength observations to characterize its color changes and Junocam images to unveil its detailed structure. The color and the altitude-opacity indices show that the oval is higher and has redder clouds than its environment but has lower cloud tops than other large ovals like the GRS, and it is less red than the GRS and oval BA. We show that in spite of the dramatic environmental changes suffered by the oval during all these years, its main characteristics are stable in time and therefore must be related with the atmospheric dynamics below the observable cloud decks.

We present updated measurements of the Crab pulsar glitch of 2019 July 23 using a dataset of pulse arrival times spanning $\sim$5 months. On MJD 58687, the pulsar underwent its seventh largest glitch observed to date, characterised by an instantaneous spin-up of $\sim$1 $\mu$Hz. Following the glitch the pulsar's rotation frequency relaxed exponentially towards pre-glitch values over a timescale of approximately one week, resulting in a permanent frequency increment of $\sim$0.5 $\mu$Hz. Due to our semi-continuous monitoring of the Crab pulsar, we were able to partially resolve a fraction of the total spin-up. This delayed spin-up occurred exponentially over a timescale of $\sim$18 hours. This is the fourth Crab pulsar glitch for which part of the initial rise was resolved in time and this phenomenon has not been observed in any other glitching pulsars, offering a unique opportunity to study the microphysical processes governing interactions between the neutron star interior and the crust.

For astrophysical reacting flows, operator splitting is commonly used to couple hydrodynamics and reactions. Each process operates independent of one another, but by staggering the updates in a symmetric fashion (via Strang splitting) second order accuracy in time can be achieved. However, approximations are often made to the reacting system, including the choice of whether or not to integrate temperature with the species. Here we demonstrate through a simple convergence test that integrating an energy equation together with reactions achieves the best convergence when modeling reactive flows with Strang splitting. Additionally, second order convergence cannot be achieved without integrating an energy or temperature equation.

Results from optical CCD photometric observations of 13 pre-main-sequence stars collected during the period from February 2007 to November 2020 are presented. These stars are located in the association Cepheus OB3, in the field of the young star V733 Cephei. Photometric observations, especially concerning the long-term variability of the stars, are missing in the literature. We present the first long-term $V(RI)_{c}$ monitoring for them, that cover 13 years. Results from our study indicate that all of the investigated stars manifest strong photometric variability. The presented paper is a part of our program for the photometric study of PMS stars located in active star-forming regions.

Dark matter might be made of "warm" particles, such as sterile neutrinos in the keV mass range, which can decay into photons through mixing and are consequently detectable by X-ray telescopes. Axionlike particles (ALPs) are detectable by X-ray telescopes too when coupled to standard model particles and decay into photons in the keV range. Both particles could explain the unidentified 3.5 keV line and, interestingly, XENON1T observed an excess of electron recoil events most prominent at 2-3 keV. One explanation could be an ALPs origin, which is not yet excluded by X-ray constraints in an anomaly-free symmetry model in which the photon production is suppressed. We study the diffuse emission coming from the Galactic halo, and calculate the sensitivity of all-sky X-ray survey performed by eROSITA to identify a sterile neutrino or ALP dark matter. We estimate bounds on the mixing angle of the sterile neutrinos and coupling strength of the ALPs. After four years of data-taking by eROSITA, we expect to set stringent constraints, and in particular, we expect to firmly probe mixing angle $\sin^2(2\theta)$ up to nearly two orders magnitude below the best-fit value for explaining the unidentified 3.5 keV line. Moreover, with eROSITA, we will be able to probe the ALP parameter space of couplings to photons and electrons, and potentially confirm an ALP origin of the XENON1T excess.

Dwarf spheroidal galaxies are dark matter dominated systems, and as such, ideal for indirect dark matter searches. If dark matter decays into high-energy photons in the dwarf galaxies, they will be a good target for current and future generations of X-ray and gamma-ray telescopes. By adopting the latest estimates of density profiles of dwarf galaxies in the Milky Way, we revise the estimates dark matter decay rates in dwarf galaxies; our results are more robust, but weaker than previous estimates. Applying these results, we study the detectability of dark matter decays with X-ray and very-high-energy gamma-ray telescopes, such as eROSITA, XRISM, Athena, HAWC, and CTA. Our projection shows that all of these X-ray telescopes will be able to critically assess the claim of the 7 keV sterile neutrino decays from stacked galaxy clusters and nearby galaxies. For TeV decaying dark matter, we can constrain its lifetime to be longer than $\sim$10$^{27}$-10$^{28}$ s. We also make projections for future dwarf galaxies that would be newly discovered with the Vera Rubin Observatory Legacy Survey of Space and Time, which will further improve the expected sensitivity to dark matter decays both in the keV and PeV mass ranges.

In this paper we examine the factors that shape the distribution of molecular gas surface densities on the 150 pc scale across 67 morphologically diverse star-forming galaxies in the PHANGS-ALMA CO (2-1) survey. Dividing each galaxy into radial bins, we measure molecular gas surface density contrasts, defined here as the ratio between a fixed high percentile of the CO distribution and a fixed reference level in each bin. This reference level captures the level of the faint CO floor that extends between bright filamentary features, while the intensity level of the higher percentile probes the structures visually associated with bright, dense ISM features like spiral arms, bars, and filaments. We compare these contrasts to matched percentile-based measurements of the 3.6 $\mu$m emission measured using Spitzer/IRAC imaging, which trace the underlying stellar mass density. We find that the logarithms of CO contrasts on 150 pc scales are 3-4 times larger than, and positively correlated with, the logarithms of 3.6 $\mu$m contrasts probing smooth non-axisymmetric stellar bar and spiral structures. The correlation appears steeper than linear, consistent with the compression of gas as it flows supersonically in response to large-scale stellar structures, even in the presence of weak or flocculent spiral arms. Stellar dynamical features appear to play an important role in setting the cloud-scale gas density in our galaxies, with gas self-gravity perhaps playing a weaker role in setting the 150 pc-scale distribution of gas densities.

The See Change survey was designed to make $z>1$ cosmological measurements by efficiently discovering high-redshift Type Ia supernovae (SNe Ia) and improving cluster mass measurements through weak lensing. This survey observed twelve galaxy clusters with the Hubble Space Telescope spanning the redshift range $z=1.13$ to $1.75$, discovering 57 likely transients and 27 likely SNe Ia at $z\sim 0.8-2.3$. As in similar previous surveys (Dawson et al. 2009), this proved to be a highly efficient use of HST for SN observations; the See Change survey additionally tested the feasibility of maintaining, or further increasing, the efficiency at yet higher redshifts, where we have less detailed information on the expected cluster masses and star-formation rates. We find that the resulting number of SNe Ia per orbit is a factor of $\sim 8$ higher than for a field search, and 45% of our orbits contained an active SN Ia within 22 rest-frame days of peak, with one of the clusters by itself yielding 6 of the SNe Ia. We present the survey design, pipeline, and SN discoveries. Novel features include fully blinded SN searches, the first random forest candidate classifier for undersampled IR data (with a 50% detection threshold within 0.05 magnitudes of human searchers), real-time forward-modeling photometry of candidates, and semi-automated photometric classifications and follow-up forecasts. We also describe the spectroscopic follow-up, instrumental in measuring host-galaxy redshifts. The cosmology analysis of our sample will be presented in a companion paper.

We investigate possible reasons for the absence of historical records of the supernova of 1054 in Europe. At the same time, we search for the new evidences as well. We establish that the previously acclaimed 'Arabic' records from ibn Butlan originate from Europe. As one of the most prominent scientists of the era, he was in Constantinople at the time of the supernova and actively participated in the medieval Church feud known as the Great Schism. Next, we reconstruct the European sky at the time of the event and find that the 'new star' (SN 1054) was in the west while the planet Venus was on the opposite side of the sky (in the east) with the Sun sited directly between these two equally bright objects, as documented in East-Asian records.

We develop the mathematical theory for an automatic, space-based system to deflect NEOs by virtue of missiles shot from the Earth-Moon L1 and L3 Lagrangian Points. A patent application has been filed for the relevant code dubbed AsterOFF (=Asteroids OFF !). This code was already implemented, and a copyright for it was registered. In a paper published in Acta Astronautica, this author proved mathematically the following theorem: "Within the sphere of influence of the Earth, any NEO could be hit by a missile at just an angle of 90 degrees, was the missile shot from the Lagrangian Points L1 or L3 of the Earth-Moon system, rather than from the surface of the Earth". As a consequence, the hitting missile would have to move along a "confocal ellipse" (centered at the Earth) uniquely determined by the NEO's incoming hyperbola. The author further shows in this paper that: 1) The proposed defense system would be ideal to deflect NEOs that are small, i.e. less than one kilometer in diameter. 2) The traditional theory of Keplerian orbits can successfully be applied to get an excellent first-order approximation of the mathematical formulae of the energy-momentum requested to achieve the NEO deflection. Many engineering details about the missiles shot from L1 and L3, however, still have to be implemented into our simulations. 3) Was one missile not enough to deflect the NEO completely, it is a great advantage of the "confocal conics" used here that the new, slightly deflected NEO's hyperbola would certainly be hit at nearly 90 degrees by another and slightly more eccentric elliptical missile trajectory. A sufficient number of missiles could thus be launched in a sequence from the Earth-Moon Lagrangian points L1 and L3 with the result that the SUM of all these small and repeated deflections will finally throw the NEO off its collision hyperbola with the Earth.

We present the first linear polarization measurements from the 2015 long-duration balloon flight of SPIDER, an experiment designed to map the polarization of the cosmic microwave background (CMB) on degree angular scales. Results from these measurements include maps and angular power spectra from observations of 4.8% of the sky at 95 and 150 GHz, along with the results of internal consistency tests on these data. While the polarized CMB anisotropy from primordial density perturbations is the dominant signal in this region of sky, Galactic dust emission is also detected with high significance; Galactic synchrotron emission is found to be negligible in the SPIDER bands. We employ two independent foreground-removal techniques in order to explore the sensitivity of the cosmological result to the assumptions made by each. The primary method uses a dust template derived from Planck data to subtract the Galactic dust signal. A second approach, employing a joint analysis of SPIDER and Planck data in the harmonic domain, assumes a modified-blackbody model for the spectral energy distribution of the dust with no constraint on its spatial morphology. Using a likelihood that jointly samples the template amplitude and $r$ parameter space, we derive 95% upper limits on the primordial tensor-to-scalar ratio from Feldman-Cousins and Bayesian constructions, finding $r<0.11$ and $r<0.19$, respectively. Roughly half the uncertainty in $r$ derives from noise associated with the template subtraction. New data at 280 GHz from SPIDER's second flight will complement the Planck polarization maps, providing powerful measurements of the polarized Galactic dust emission.

The decade-long Vista Variables in the Via Lactea (VVV) survey has detected numerous highly variable young stellar objects (YSOs). We present a study of 61 highly variable VVV YSOs ($\Delta K_s = 1$-5~mag), combining near infrared spectra from Magellan and VLT with VVV and NEOWISE light curves to investigate physical mechanisms behind eruptive events. Most sources are spectroscopically confirmed as eruptive variables (typically Class I YSOs) but variable extinction is also seen. Among them, magnetically controlled accretion, identified by H{\sc i} recombination emission (usually accompanied by CO emission), is observed in 46 YSOs. Boundary layer accretion, associated with FU~Ori-like outbursts identified by CO overtone and H$_2$O absorption, is observed only in longer duration events ($\ge$5~yr total duration). However, even in long duration events, the magnetically controlled accretion mode predominates, with amplitudes similar to the boundary layer mode. Shorter (100-700~days) eruptive events usually have lower amplitudes and these events are generally either periodic accretors or multiple timescale events, wherein large photometric changes occur on timescales of weeks and years. We find that the ratio of amplitudes in $K_s$ and $W2$ can distinguish between variable accretion and variable extinction. Several YSOs are periodic or quasi-periodic variables. We identify examples of periodic accretors and extinction-driven periodicity among them (with periods up to 5~yr) though more data are needed to classify some cases. The data suggest that dynamic interactions with a companion may control the accretion rate in a substantial proportion of eruptive systems, although star-disc interactions should also be considered.

According to planetary interior models, some giant planets contain large metal masses with large metal-mass fractions. HD 149026b and TOI-849b are characteristic examples of these giant planets. It has been suggested that the envelope mass loss during giant impacts plays a key role in the formation of such giant planets. The aim of the present letter is to propose a mechanism that can explain the origin of such giant planets. We investigate the formation of giant planets in a rapidly dissipating disk using N-body simulations that consider pebble accretion. The results show that although the pebble isolation mass is smaller than the metal mass (> 30 Earth masses) in some giant planets, the interior metal mass can be increased by giant impacts between planets with the isolation mass. Regarding the metal fraction, the cores accrete massive envelopes by runaway gas accretion during the disk-dissipation phase of 1-10 Myr in a disk that evolves without photoevaporation. Although a large fraction of the envelope can be lost during giant impacts, the planets can reaccrete the envelope after impacts in a slowly dissipating disk. Here, we demonstrate that, by photoevaporation in a rapidly dissipating disk, the runaway gas accretion is quenched in the middle, resulting in the formation of giant planets with large metal-mass fractions. The origins of HD 149026b and TOI-849b, which are characterized by their large metal-mass fractions, can be naturally explained by a model that considers a disk evolving with photoevaporation.

We investigate the photoionised X-ray emission line regions (ELRs) within the Seyfert 2 galaxy NGC 1068, to determine if there are any characteristic changes between observations taken fourteen years apart. We compare XMM-Newton observations collected in 2000 and 2014, simultaneously fitting the reflection grating spectrometer (RGS) and EPIC-pn spectra of each epoch, for the first time, with the photoionisation model, PION, in SPEX. We find that four PION components are required to fit the majority of the emission lines in the spectra of NGC 1068, with $\log \xi=1-4$, $\log N_H>26 m^{-2}$, and $v_{out}=-100$ to $-600 kms^{-1}$ for both epochs. Comparing the ionisation state of the components shows almost no difference between the two epochs, while there is an increase in the total column density. To estimate the locations of these plasma regions from the central black hole we compare distance methods, excluding the variability arguments as there is no spectral change between observations. Although the methods are unable to constrain the distances, the locations are consistent with the narrow line region, with the possibility of the higher ionised component being part of the broad line region, but we cannot conclude this for certain. In addition, we find evidence for emission from collisionally ionised plasma, while previous analysis had suggested that collisional plasma emission was unlikely. However, although PION is unable to account for the FeXVII emission lines at 15 and 17 \AA, we do not rule out that photoexcitation is a valid processes to produce these lines too. NGC 1068 has not changed, both in terms of the observed spectra or from our modelling, within the 14 year time period between observations. This suggests that the ELRs are fairly static relative to the 14 year time frame between observations, or there is no dramatic change in the black hole variability.

We use X-ray and infrared observations to study the properties of three classes of young stars in the Carina Nebula: intermediate-mass (2--8M$_\odot$) pre-main sequence stars (IMPS; i.e. intermediate-mass T Tauri stars), late-B and A stars on the zero-age main sequence (AB), and lower-mass T Tauri stars (TTS). We divide our sources among these three sub-classifications and further identify disk-bearing young stellar objects versus diskless sources with no detectable infrared (IR) excess emission using IR (1--8 $\mu$m) spectral energy distribution modeling. We then perform X-ray spectral fitting to determine the hydrogen absorbing column density ($N_{\rm H}$), absorption-corrected X-ray luminosity ($L_{\rm X}$), and coronal plasma temperature ($kT$) for each source. We find that the X-ray spectra of both IMPS and TTS are characterized by similar $kT$ and $N_{\rm H}$, and on average $L_{\rm X}$/$L_{\rm bol} \sim4\times10^{-4}$. IMPS are systematically more luminous in X-rays (by $\sim$0.3 dex) than all other sub-classifications, with median $L_{\rm X} = 2.5\times10^{31}$ erg s$^{-1}$, while AB stars of similar masses have X-ray emission consistent with TTS companions. These lines of evidence converge on a magneto-coronal flaring source for IMPS X-ray emission, a scaled-up version of the TTS emission mechanism. IMPS therefore provide powerful probes of isochronal ages for the first $\sim$10 Myr in the evolution of a massive stellar population, because their intrinsic, coronal X-ray emission decays rapidly after they commence evolving along radiative tracks. We suggest that the most luminous (in both X-rays and IR) IMPS could be used to place empirical constraints on the location of the intermediate-mass stellar birth line.

The modeling of fast flavor evolution of neutrinos in dense environments has been traditionally carried out by relying on a two flavor approximation for simplicity. In addition, vacuum mixing has been deemed negligible. For the first time, we highlight that the fast flavor evolution in three flavors is intrinsically different from the one obtained in the two flavor approximation. This is due to the exponential growth of flavor mixing in the $e$--$\mu$ and $e$--$\tau$ sectors generated by the vacuum term in the Hamiltonian. As a result, substantially larger flavor mixing is found in three flavors. Our findings highlight that the two flavor approximation is not justified for fast pairwise conversion, even if the angular distributions of non-electron type neutrinos are initially identical.

Due to the lower ionospheric thermal pressure and existence of the crustal magnetism at Mars, the Martian ionopause is expected to behave differently from the ionopause at Venus. We study the solar wind interaction and pressure balance at the ionopause of Mars using both in situ and remote sounding measurements from the MARSIS (Mars Advanced Radar for Subsurface and Ionosphere Sounding) instrument. We show that the magnetic pressure usually dominates the thermal pressure to hold off the solar wind in the ionopause at Mars, with only 13\% unmagnetized ionopauses observed over a 11-year period. We also find that the ionopause altitude decreases as the normal component of the solar wind dynamic pressure increases. Moreover, our results show that the ionopause thickness at Mars is mainly determined by the ion gyromotion and equivalent to about 5.7 ion gyroradii.

We motivate a minimal realization of slow-roll $k$-inflation by incorporating the local conformal symmetry and the broken global $\mathrm{SO}(1,1)$ symmetry in the metric-affine geometry. With use of the metric-affine geometry where both the metric and the affine connection are treated as independent variables, the local conformal symmetry can be preserved in each term of the Lagrangian and thus higher derivatives of scalar fields can be easily added in a conformally invariant way. Predictions of this minimal slow-roll $k$-inflation, $n_\mathrm{s}\sim 0.96$, $r\sim 0.005$, and $c_\mathrm{s}\sim 0.03$, are not only consistent with current observational data but also have a prospect to be tested by forthcoming observations.

Based on the recently formulated chiral radiation transport theory for left-handed neutrinos, we study the chiral transport of neutrinos near thermal equilibrium in core-collapse supernovae. We first compute the near-equilibrium solution of the chiral radiation transport equation under the relaxation time approximation, where the relaxation time is directly derived from the effective field theory of the weak interaction. By using such a solution, we systematically derive analytic expressions for the nonequilibrium corrections of the neutrino energy-momentum tensor and neutrino number current induced by magnetic fields via the neutrino absorption on nucleons. In particular, we find the nonequilibrium neutrino energy current proportional to the magnetic field. We also discuss its phenomenological consequences such as the possible relation to pulsar kicks.

In the presence of the fluid helicity $\boldsymbol{v} \cdot \boldsymbol{\omega}$, the magnetic field induces an electric current of the form $\boldsymbol{j} = C_{\rm HME} (\boldsymbol{v} \cdot \boldsymbol{\omega}) \boldsymbol{B}$. This is the helical magnetic effect (HME). We show that for massless Dirac fermions with charge $e=1$, the transport coefficient $C_{\rm HME}$ is fixed by the chiral anomaly coefficient $C=1/(2\pi^2)$ as $C_{\rm HME} = C/2$ independently of interactions. We show the conjecture that the coefficient of the magneto-vorticity coupling for the local vector charge, $n = C_{B \omega} \boldsymbol{B} \cdot \boldsymbol{\omega}$, is related to the chiral anomaly coefficient as $C_{B \omega} = C/2$. We also discuss the condition for the emergence of the helical plasma instability that originates from the HME.

We consider an axion cloud around a black hole with background magnetic fields. We calculate the decay rate of the axion cloud due to the axion-photon conversion associated with the axion-photon coupling. For simplicity, we consider the situation where the axion configuration is dominated by a solution for the eigenvalue equation equivalent to that for the Hydrogen atom, and the coupling term can be evaluated by a successive perturbation method. For the monopole background, we find the decay rate of the axion cloud is given by $\sim q^2\kappa^2(GM)^5\mu^8$, where $\mu$, $M$, $G$, $\kappa$ and $q$ are the axion mass, black hole mass, gravitational constant, coupling constant of the axion-photon coupling and monopole charge, respectively. For the uniform background magnetic field, we obtain the decay rate of the axion cloud $\sim B_0^2\kappa^2 (GM)^7\mu^6$, where $B_0$ is the magnetic field strength. Applying our formula to the central black hole in our galaxy, we find that the value of the decay rate for the case of the uniform magnetic field is comparable to the growth rate of the superradiant instability with $\kappa\sim 10^{-11}{\rm GeV^{-1}}$, $B_0\sim 10^2{\rm G}$ and $\mu\sim 10^{-18}{\rm eV}$. The ratio is $10^5$ times larger for the monopole magnetic field with the same values of the parameters.

We use 3D fully kinetic particle-in-cell simulations to study the occurrence of magnetic reconnection in a simulation of decaying turbulence created by anisotropic counter-propagating low-frequency Alfv\'en waves consistent with critical-balance theory. We observe the formation of small-scale current-density structures such as current filaments and current sheets as well as the formation of magnetic flux ropes as part of the turbulent cascade. The large magnetic structures present in the simulation domain retain the initial anisotropy while the small-scale structures produced by the turbulent cascade are less anisotropic. To quantify the occurrence of reconnection in our simulation domain, we develop a new set of indicators based on intensity thresholds to identify reconnection events in which both ions and electrons are heated and accelerated in 3D particle-in-cell simulations. According to the application of these indicators, we identify the occurrence of reconnection events in the simulation domain and analyse one of these events in detail. The event is related to the reconnection of two flux ropes, and the associated ion and electron exhausts exhibit a complex three-dimensional structure. We study the profiles of plasma and magnetic-field fluctuations recorded along artificial-spacecraft trajectories passing near and through the reconnection region. Our results suggest the presence of particle heating and acceleration related to small-scale reconnection events within magnetic flux ropes produced by the anisotropic Alfv\'enic turbulent cascade in the solar wind. These events are related to current structures of order a few ion inertial lengths in size.