Search

Educational research highlights that improved careers education is needed to increase participation in science, technology, engineering, and mathematics (STEM). Current UK careers resources concerning the space sector, however, are found to perhaps not best reflect the diversity of roles present and may in fact perpetuate misconceptions about the usefulness of science. We, therefore, compile a more diverse set of space-related jobs, which will be used in the development of a new space careers resource.

Many observations of protoplanetary discs studied with ALMA have revealed the complex substructure present in the discs. Rings and gaps in the dust continuum are now a common sight in many discs; however, their origins still remain unknown. We look at all protoplanetary disc images taken with ALMA from cycles 0 to 5 and find that 56 discs show clear substructure. We further study the 56 discs and classify the morphology seen according to four categories: Rim, Ring, Horseshoe, and Spiral. We calculate the ages of the host stars using stellar isochrones and investigate the relation between the morphology of the substructure seen in the protoplanetary discs and the age of the host stars. We find that there is no clear evolutionary sequence in the protoplanetary discs as the stars increase in age, although there is a slight tendency for spirals to appear in younger systems and horseshoes to be seen in more evolved systems. We also show that majority of the images of protoplanetary discs made by ALMA may not have had a sufficiently high resolution or sensitivity to resolve substructure in the disc. We show that angular resolution is important in detecting substructure within protoplanetary discs, with sensitivity distinguishing between the different types of substructure. We compare the substructure seen in protoplanetary discs at sub-mm to those seen in scattered light. We find that cavities are a common substructure seen in discs at both sub-mm wavelengths and in scattered light.

R-mode oscillations in a rotating star produce characteristic signatures in a Fourier amplitude spectrum at frequencies related with the rotation frequency, which can be, in turn, used to obtain the surface rotation rate of the star. Some binary stars observed by Kepler indicate the presence of r~modes that are probably excited by the tidal effect. In this paper, we have obtained stellar rotation periods in 20 eccentric (heartbeat) binaries with r-mode signatures. The majority of the rotation periods are found to be comparable to pseudo-synchronous periods, in which the angular velocity of rotation is similar to the angular orbital motion of the companion at periastron. In particular, for the heartbeat stars with orbital periods longer than about 8\,d, all but one agree with pseudo-synchronous rotation. In contrast to a previous investigation by Zimmerman et al., our result supports the pseudo-synchronisation theory developed by Hut., Comment: 13 pages, to appear in MNRAS

We use a sample of $z=0$ galaxies visually classified as slow rotators (SRs) in the EAGLE hydrodynamical simulations to explore the effect of galaxy mergers on their formation, characterise their intrinsic galaxy properties, and study the connection between quenching and kinematic transformation. SRs that have had major or minor mergers (mass ratios $\ge 0.3$ and $0.1-0.3$, respectively) tend to have a higher triaxiality parameter and ex-situ stellar fractions than those that had exclusively very minor mergers or formed in the absence of mergers ("no-merger" SRs). No-merger SRs are more compact, have lower black hole-to-stellar mass ratios and quenched later than other SRs, leaving imprints on their $z=0$ chemical composition. For the vast majority of SRs we find that quenching, driven by active galactic nuclei feedback, precedes kinematic transformation, except for satellite SRs, in which these processes happen in tandem. However, in $\approx 50$\% of these satellites, satellite-satellite mergers are responsible for their SR fate, while environment (i.e. tidal field and interactions with the central) can account for the transformation in the rest. By splitting SRs into kinematic sub-classes, we find that flat SRs prefer major mergers; round SRs prefer minor or very minor mergers; prolate SRs prefer gas-poor mergers. Flat and prolate SRs are more common among satellites hosted by massive halos ($>10^{13.6}\,\rm M_{\odot}$) and centrals of high masses ($M_{\star} > 10^{10.5}\, \rm M_{\odot}$). Although EAGLE galaxies display kinematic properties that broadly agree with observations, there are areas of disagreement, such as inverted stellar age and velocity dispersion profiles. We discuss these and how upcoming simulations can solve them., Accepted for publications in MNRAS. 22 pages and 20 figures. Main difference with previous version is additional panel in Fig. 7 showing Sersic index vs stellar mass for different merger histories

The response of the Earth’s magnetotail current sheet to the external solar wind driver is highly time-dependent and asymmetric. For example, the current sheet twists in response to variations in the By component of the interplanetary magnetic field (IMF), and is hinged by the dipole tilt. Understanding the timescales over which these asymmetries manifest is of particular importance during geomagnetic storms when the dynamics of the tail control substorm activity. To investigate this, we use the Gorgon MHD model to simulate a geomagnetic storm which commenced on 3 May 2014, and was host to multiple By and Bz reversals and a prolonged period of southward IMF driving. We find that the twisting of the current sheet is well-correlated to IMF By throughout the event, with the angle of rotation increasing linearly with downtail distance and being more pronounced when the tail contains less open flux. During periods of southward IMF the twisting of the central current sheet responds most strongly at a timelag of ∼ 100 min for distances beyond 20 RE, consistent with the 1–2 h convection timescale identified in the open flux content. Under predominantly northward IMF the response of the twisting is bimodal, with the strongest correlations between 15 and 40 RE downtail being at a shorter timescale of ∼ 30 min consistent with that estimated for induced By due to wave propagation, compared to a longer timescale of ∼ 3 h further downtail again attributed to convection. This indicates that asymmetries in the magnetotail communicated by IMF By are influenced mostly by global convection during strong solar wind driving, but that more prompt induced By effects can dominate in the near-Earth tail and during periods of weaker driving. These results provide new insight into the characteristic timescales of solar wind-magnetosphere-ionosphere coupling.

We present a statistical investigation (September 2012 - September 2017) of pitch angle distribution (PAD) of energetic electrons (∼30 keV - 1 MeV) in the outer radiation belt (L ≥ 3) during CME- and CIR-driven geomagnetic storms using Van Allen Probe measurements. We selected geomagnetic storms based on minimum of SYM-H being less than -50 nT and classified the storms according to their drivers. Thus, we obtained 23 CME- and 24 CIR-driven storms. During the storm intervals, pitch angle resolved electron flux measurements are obtained from the MagEIS instrument on-board Van Allen Probe-A spacecraft. We assume symmetric pitch angle distributions around 90° pitch angle and fit the observed PADs with Legendre polynomials after propagating them to the magnetic equator. Legendre coefficients c2 and c4, and the ratio R = |c2/c4| are used to categorize the different PAD types. To resolve the spatio-temporal distribution of PADs, these coefficients are binned in 5 L-shell bins, 12 MLT bins for seven energy channels and four storm phases. We found that several hundreds of keV electrons exhibit clear dependence on local time, storm phases and storm drivers, with increased anisotropy for CME-driven storms during main and early recovery phases. On the contrary, we found that tens of keV electrons do not exhibit significant dependence on these parameters. We have discussed the different physical mechanisms responsible for the observed MLT dependent PADs and found drift-shell splitting to be the major contributor.

Since the beginning of the 21st century, the fields of astronomy and astrophysics have experienced significant growth at observational and computational levels, leading to the acquisition of increasingly huge volumes of data. In order to process this vast quantity of information, artificial intelligence (AI) techniques are being combined with data mining to detect patterns with the aim of modeling, classifying or predicting the behavior of certain astronomical phenomena or objects. Parallel to the exponential development of the aforementioned techniques, the scientific output related to the application of AI and machine learning (ML) in astronomy and astrophysics has also experienced considerable growth in recent years. Therefore, the increasingly abundant articles make it difficult to monitor this field in terms of which research topics are the most prolific or novel, or which countries or authors are leading them. In this article, a text-mining-based scientometric analysis of scientific documents published over the last three decades on the application of AI and ML in the fields of astronomy and astrophysics is presented. The VOSviewer software and data from the Web of Science (WoS) are used to elucidate the evolution of publications in this research field, their distribution by country (including co-authorship), the most relevant topics addressed, and the most cited elements and most significant co-citations according to publication source and authorship. The obtained results demonstrate how application of AI/ML to the fields of astronomy/astrophysics represents an established and rapidly growing field of research that is crucial to obtaining scientific understanding of the universe.

Magnetic reconnection is a sporadic process responsible for energy release in space and laboratory plasmas. It is believed that the tearing mode instability may be responsible for the onset of reconnection in the magnetotail. However, due to its elusive nature, there is an absence of in-situ observations of the tearing instability prior to magnetic reconnection in our nearest natural plasma laboratory. Using neural network outlier detection methods in conjunction with Cluster spacecraft data, we find unique electron pitch angle distributions that are consistent with simulation predictions of the tearing instability and the subsequent evolution of plasma electrons and reconnection. We evaluate tearing stability criterion for the events identified via our neural network outlier method, and find signatures of magnetic reconnection minutes after the majority of tearing observations. Our analysis of the tearing instability provides new insights into the fundamental understanding of the mechanism responsible for reconnection, a process that is ubiquitous in different astrophysical plasma regimes across the Universe and in laboratory experiments on Earth.

Oscillatory reconnection can manifest through the interaction between the ubiquitous MHD waves and omnipresent null points in the solar atmosphere and is characterized by an inherent periodicity. In the current study, we focus on the relationship between the period of oscillatory reconnection and the strength of the wave pulse initially perturbing the null point, in a hot coronal plasma. We use the PLUTO code to solve the fully compressive, resistive MHD equations for a 2D magnetic X-point. Using wave pulses with a wide range of amplitudes, we perform a parameter study to obtain values for the period, considering the presence and absence of anisotropic thermal conduction separately. In both cases, we find that the resulting period is independent of the strength of the initial perturbation. The addition of anisotropic thermal conduction only leads to an increase in the mean value for the period, in agreement with our previous study. We also consider a different type of initial driver and we obtain an oscillation period matching the independent trend previously mentioned. Thus, we report for the first time on the independence between the type and strength of the initializing wave pulse and the resulting period of oscillatory reconnection in a hot coronal plasma. This makes oscillatory reconnection a promising mechanism to be used within the context of coronal seismology., 19 pages, 14 figures, accepted for publication in ApJ

Context. Accurate detections of frequent small-scale extreme ultraviolet (EUV) brightenings are essential to the investigation of the physical processes heating the corona. Aims. We detected small-scale brightenings, termed campfires, using their morphological and intensity structures as observed in coronal EUV imaging observations for statistical analysis. Methods. We applied a method based on Zernike moments and a support vector machine (SVM) classifier to automatically identify and track campfires observed by Solar Orbiter/Extreme Ultraviolet Imager (EUI) and Solar Dynamics Observatory (SDO)/Atmospheric Imaging Assembly (AIA). Results. This method detected 8678 campfires (with length scales between 400 km and 4000 km) from a sequence of 50 High Resolution EUV telescope (HRIEUV) 174 Å images. From 21 near co-temporal AIA images covering the same field of view as EUI, we found 1131 campfires, 58% of which were also detected in HRIEUV images. In contrast, about 16% of campfires recognized in HRIEUV were detected by AIA. We obtain a campfire birthrate of 2 × 10−16 m−2 s−1. About 40% of campfires show a duration longer than 5 s, having been observed in at least two HRIEUV images. We find that 27% of campfires were found in coronal bright points and the remaining 73% have occurred out of coronal bright points. We detected 23 EUI campfires with a duration greater than 245 s. We found that about 80% of campfires are formed at supergranular boundaries, and the features with the highest total intensities are generated at network junctions and intense H I Lyman-α emission regions observed by EUI/HRILya. The probability distribution functions for the total intensity, peak intensity, and projected area of campfires follow a power law behavior with absolute indices between 2 and 3. This self-similar behavior is a possible signature of self-organization, or even self-organized criticality, in the campfire formation process.

A detailed study of the morphology of lenticular galaxies is an important way to understand how this type of galaxy formed and evolves over time. Decomposing a galaxy into its components (disc, bulge, bar, ...) allows recovering the colour gradients present in each system, its star formation history, and its assembly history. We use GALFITM to perform a multi-wavelength structural decomposition of the closest lenticular galaxy, NGC 3115, resulting in the description of its stellar light into several main components: a bulge, a thin disc, a thick disc and also evidence of a bar. We report the finding of central bluer stellar populations in the bulge, as compared to the colour of the galaxy outskirts, indicating either the presence of an Active Galactic Nucleus (AGN) and/or recent star formation activity. From the spectral energy distribution results, we show that the galaxy has a low luminosity AGN component, but even excluding the effect of the nuclear activity, the bulge is still bluer than the outer-regions of the galaxy, revealing a recent episode of star formation. Based on all of the derived properties, we propose a scenario for the formation of NGC 3115 consisting of an initial gas-rich merger, followed by accretions and feedback that quench the galaxy, until a recent encounter with the companion KK084 that reignited the star formation in the bulge, provoked a core displacement in NGC 3115 and generated spiral-like features. This result is consistent with the two-phase formation scenario, proposed in previous studies of this galaxy., Comment: Accepted for publication in MNRAS: 22 pages, 15 figures, 7 tables, 1 appendix

We present a new library of semi-empirical stellar spectra that is based on the empirical MILES library. A new, high resolution library of theoretical stellar spectra is generated that is specifically designed for use in stellar population studies. We test these models across their full wavelength range against other model libraries and find reasonable agreement in their predictions of spectral changes due to atmospheric $\alpha$-element variations, known as differential corrections. We also test the models against the MILES and MaStar libraries of empirical stellar spectra and also find reasonable agreements, as expected from previous work. We then use the abundance pattern predictions of the new theoretical stellar spectra to differentially correct MILES spectra to create semi-empirical MILES (sMILES) star spectra with abundance patterns that differ from those present in the Milky Way. The final result is 5 families of 801 sMILES stars with [$\alpha$/Fe] abundances ranging from $-$0.20 to 0.60 dex at MILES resolution (FWHM=$2.5\,${\AA}) and wavelength coverage ($3540.5-7409.6\,${\AA}). We make the sMILES library publicly available., Comment: 26 pages, 18 figures plus Supplementary Material. Accepted for publication in MNRAS. Corrected a minor font size error and two figure caption placements from the previous version

The spectrum of eight-times ionized iron, Fe IX, was studied in the 110-200 {\AA} region. A low inductance vacuum spark and a 3-m grazing incidence spectrograph were used for the excitation and recording of the spectrum. Previous analyses of Fe IX have been greatly extended and partly revised. The numbers of known lines in the 3p^53d - 3p^54f and 3p^53d - 3p^43d^2 transition arrays are extended to 25 and 81, respectively. Most of the identifications of the Fe IX lines from the 3p^53d - 3p^43d^2 transition array in the solar spectrum have been confirmed and several new identifications are suggested., Comment: Accepted by ApJ, 24 pages, 7 figures, 5 tables

A leap forward in our understanding of particle energization in plasmas throughout the heliosphere is essential to answer longstanding questions in heliophysics, including the heating of the solar corona, acceleration of the solar wind, and energization of particles that lead to observable phenomena, such as the Earth’s aurora. The low densities and high temperatures of typical heliospheric environments lead to weakly collisional plasma conditions. Under these conditions, the energization of particles occurs primarily through collisionless interactions between the electromagnetic fields and the individual plasma particles with energies characteristic of a particular interaction. To understand how the plasma heating and particle acceleration impacts the macroscopic evolution of the heliosphere, impacting phenomena such as extreme space weather, it is critical to understand these collisionless wave-particle interactions on the characteristic ion and electron kinetic timescales. Such understanding requires high-cadence measurements of both the electromagnetic fields and the three-dimensional particle velocity distributions. Although existing instrument technology enables these measurements, a major challenge to maximize the scientific return from these measurements is the limited amount of data that can be transmitted to the ground due to telemetry constraints. A valuable, but underutilized, approach to overcome this limitation is to compute on-board correlations of the maximum-cadence field and particle measurements to improve the sampling time by several orders of magnitude. Here we review the fundamentals of the innovative field-particle correlation technique, present a formulation of the technique that can be implemented as an on-board wave-particle correlator, and estimate results that can be achieved with existing instrumental capabilities for particle velocity distribution measurements.

We present the first observation of a solar prominence at $84-116$ GHz using the high resolution interferometric imaging of ALMA. Simultaneous observations in H$\alpha$ from Bia{\l}kaw Observatory and with SDO/AIA reveal similar prominence morphology to the ALMA observation. The contribution functions of 3 mm and H$\alpha$ emission are shown to have significant overlap across a range of gas pressures. We estimate the maximum millimetre-continuum optical thickness to be $\tau_\mathrm{3mm}\approx 2$, and the brightness temperature from the observed H$\alpha$ intensity. The brightness temperature measured by ALMA is $\sim 6000-7000$ K in the prominence spine, which correlates well with the estimated brightness temperature for a gas temperature of 8000 K., Comment: 5 pages, 7 figures

Using a combination of Juno magnetometer and plasma data, we show evidence of Alfvénic turbulence within the mid‐to‐high latitude magnetosphere with sufficient conditions to trigger auroral particle acceleration. We analyze 12 events that, in agreement with theoretical results, are found to be dissipative at the electron inertial scale. Furthermore, these events contain significant Poynting flux in the range ∼0.8–20 mW/m2 at ionospheric altitudes. This is sufficient to generate auroral emissions. We confirm that such events are incompressible, confirming their Alfvénicity, occur at dissipative scales, have intermittent features present and are multifractal in nature. These results illustrate the importance of turbulence in the mid‐to‐high latitudes of Jupiter's magnetosphere as a driver of particle acceleration.

We develop and apply a bespoke fitting routine to a large volume of solar wind electron distribution data measured by Parker Solar Probe (PSP) over its first five orbits, covering radial distances from 0.13 to 0.5 au. We characterise the radial evolution of the electron core, halo and strahl populations in the slow solar wind during these orbits. The fractional densities of these three electron populations provide evidence for the growth of the combined suprathermal halo and strahl populations from 0.13 to 0.17 au. Moreover, the growth in the halo population is not matched by a decrease of the strahl population at these distances, as has been reported for previous observations at distances greater than 0.3 au. We also find that the halo is negligible at small heliocentric distances. The fractional strahl density remains relatively constant ~1 % below 0.2 au, suggesting that the rise in the relative halo density is not solely due to the transfer of strahl electrons into the halo., Published in ApJ

We present an axi-symmetric model for the ultraviolet (UV)-to-submillimetre (submm) images of the nearly face-on spiral galaxy NGC 628. It was calculated using a radiative transfer (RT) code, accounting for the absorption and re-emission of starlight by dust in the interstellar medium of this galaxy. The code incorporates emission from Polycyclic Aromatic Hydrocarbons, anisotropic scattering and stochastic heating of the grains. This is the second successful modelling of a face-on spiral galaxy with RT methods, whereby the large-scale geometry of stars and dust is self-consistently determined. The solution was obtained by fitting azimuthally averaged profiles in the UV, optical and submm. The model predicts remarkably well all characteristics of the profiles, including the increase by a factor of 1.8 of the scale-length of the infrared emissivity between 70 and 500 $μ$m. We find that NGC 628 did not undergo an efficient inside-out disk growth, as predicted by semi-analytical hierarchical models for galaxy formation. We also find large amounts of dust grains at large radii, that could involve efficient transport mechanisms from the inner disk. Our results show that 71% of the dust emission in NGC 628 is powered by the young stellar populations, with the old stellar populations from the bulge contributing 65% to the heating of the dust in the central region ($R, 27 pages, 18 figures. Accepted for publication in MNRAS

Gaia has revealed clear evidence of bending waves in the vertical kinematics of stars in the Solar Neighbourhood. We study bending waves in two simulations, one warped, with the warp due to misaligned gas inflow, and the other unwarped. We find slow, retrograde bending waves in both models, with the ones in the warped model having larger amplitudes. We also find fast, prograde bending waves. Prograde bending waves in the unwarped model are very weak, in agreement with the expectation that these waves should decay on short, ~ crossing, timescales, due to strong winding. However, prograde bending waves are much stronger for the duration of the warped model, pointing to irregular gas inflow along the warp as a continuous source of excitation. We demonstrate that large amplitude bending waves that propagate through the Solar Neighbourhood give rise to a correlation between the mean vertical velocity and the angular momentum, with a slope consistent with that found by Gaia. The bending waves affect populations of all ages, but the sharpest features are found in the young populations, hinting that short wavelength waves are not supported by the older, kinematically hotter, populations. Our results demonstrate the importance of misaligned gas accretion as a recurrent source of vertical perturbations of disc galaxies, including in the Milky Way., Comment: 21 pages, 20 figures; Published on MNRAS

We present new imaging of a sub-kpc-scale binary active galactic nucleus (AGN) candidate from the Karl G. Jansky Very Large Array (VLA) and the Multi-Element Radio Linked Interferometer Network (e-MERLIN). Two unresolved radio sources of similar luminosity around 1022 WHz−1 are identified in ∼35 h of e-MERLIN 6 cm imaging. These radio sources have an angular separation of 0.19 ± 0.06 arcsec and position angle (PA) of 22° ± 10°, corresponding to a projected separation of 0.95 ± 0.29 kpc at the epoch of the source. Our results suggest the presence of a kpc-scale active black hole pair hosted by two galaxies in the late stage of a merger at z = 0.35. This work follows Woo et al., which presented two optical sources with a similar separation and PA, and a velocity separation of 200 km s−1. Our target adds to the currently limited sample of close-separation binary AGNs, which will aid in understanding the frequency of mergers and the stochastic gravitational wave background.

The near-relativistic electron population in the outer Van Allen radiation belt is highly dynamic and strongly coupled to geomagnetic activity such as storms and substorms, which are driven by the interaction of the magnetosphere with the solar wind. The energy, content and spatial extent of electrons in the outer radiation belt can vary on timescales of hours to days, dictated by the continuously evolving influence of acceleration and loss processes. While net changes in the electron population are directly observable, the relative influence of different processes is far from fully understood. Using a continuous 12-year dataset from the Proton Electron Telescope (PET) on board the Solar Anomalous Magnetospheric Particle Explorer (SAMPEX), we statistically compare the relative variations of trapped electrons to those in the bounce loss cone. Our results show that there is a proportional increase in flux entering the bounce loss cone outside the plasmapause during storm main phase and early recovery phase. Loss enhancement is sustained on the dawnside throughout the recovery phase while loss on the duskside is enhanced around minimum Sym-H and quickly diminishes. Spatial variations are also examined in relation to geomagnetic activity, making comparisons to possible causal wave modes such as whistler-mode chorus and plasmaspheric hiss.

The major-axis density profiles of bars are known to be either exponential or 'flat'. We develop an automated non-parametric algorithm to detect flat profiles and apply it to a suite of simulations (with and without gas). We demonstrate that flat profiles are a manifestation of a bar's secular growth, producing a 'shoulder' region (an overdensity above an exponential) in its outskirts. Shoulders are not present when bars form, but develop as the bar grows. If the bar does not grow, shoulders do not form. Shoulders are often accompanied by box/peanut bulges, but develop separately from them and are independent tracers of a bar's growth. They can be observed at a wide range of viewing orientations with only their slope varying significantly with inclination. We present evidence that shoulders are produced by looped x1 orbits. Since the growth rate of the bar moderately correlates with the growth rate of the shoulder strength, these orbits are probably recently trapped. Shoulders therefore are evidence of bar growth. The properties of the shoulders do not, however, establish the age of a bar, because secondary buckling or strong spirals may destroy shoulders, and also because shoulders do not form if the bar does not grow much. In particular, our results show that an exponential profile is not necessarily an indication of a young bar., Comment: 20 pages, 23 figures, 2 tables. Published in MNRAS

Correlated orientations of quasar optical and radio polarisation, and of radio jets, have been reported on Gpc scales, possibly arising from intrinsic alignment of spin axes. Optical quasar polarisation appears to be preferentially either aligned or orthogonal to the host large-scale structure, specifically large quasar groups (LQGs). Using a sample of 71 LQGs at redshifts $1.0 \leq z \leq 1.8$, we investigate whether LQGs themselves exhibit correlated orientation. We find that LQG position angles (PAs) are unlikely to be drawn from a uniform distribution ($p$-values $0.008 \lesssim p \lesssim 0.07$). The LQG PA distribution is bimodal, with median modes at $\bar{\theta}\sim45\pm2^{\circ}, 136\pm2^{\circ}$, remarkably close to the mean angles of quasar radio polarisation reported in two regions coincident with our LQG sample. We quantify the degree of alignment in the PA data, and find that LQGs are aligned and orthogonal across very large scales. The maximum significance is $\simeq 0.8\%$ ($2.4\sigma$) at typical angular (proper) separations of $\sim 30^{\circ}$ (1.6 Gpc). If the LQG orientation correlation is real, it represents large-scale structure alignment over scales larger than those predicted by cosmological simulations and at least an order of magnitude larger than any so far observed, with the exception of quasar-polarisation / radio-jet alignment. We conclude that LQG alignment helps explain quasar-polarisation / radio-jet alignment, but raises challenging questions about the origin of the LQG correlation and the assumptions of the concordance cosmological model., Comment: Accepted for MNRAS

Full list of authors: Baldi, R. D.; Williams, D. R. A.; McHardy, I. M.; Beswick, R. J.; Brinks, E.; Dullo, B. T.; Knapen, J. H.; Argo, M. K.; Aalto, S.; Alberdi, A.; Baan, W. A.; Bendo, G. J.; Corbel, S.; Fenech, D. M.; Gallagher, J. S.; Green, D. A.; Kennicutt, R. C.; Klöckner, H. -R.; Körding, E.; Maccarone, T. J.; Muxlow, T. W. B.; Mundell, C. G.; Panessa, F.; Peck, A. B.; Pérez-Torres, M. A.; Romero-Cañizales, C.; Saikia, P.; Shankar, F.; Spencer, R. E.; Stevens, I. R.; Varenius, E.; Ward, M. J.; Yates, J.; Uttley, P., We present the second data release of high-resolution (≤0.2 arcsec) 1.5-GHz radio images of 177 nearby galaxies from the Palomar sample, observed with the e-MERLIN array, as part of the Legacy e-MERLIN Multi-band Imaging of Nearby Galaxies Sample (LeMMINGs) survey. Together with the 103 targets of the first LeMMINGs data release, this represents a complete sample of 280 local active (LINER and Seyfert) and inactive galaxies (H ii galaxies and absorption line galaxies, ALG). This large program is the deepest radio survey of the local Universe, ≳1017.6 W Hz-1, regardless of the host and nuclear type: we detect radio emission ≳0.25 mJy beam-1 for 125/280 galaxies (44.6 per cent) with sizes of typically ≲100 pc. Of those 125, 106 targets show a core which coincides within 1.2 arcsec with the optical nucleus. Although we observed mostly cores, around one third of the detected galaxies features jetted morphologies. The detected radio core luminosities of the sample range between ∼1034 and 1040 erg s-1. LINERs and Seyferts are the most luminous sources, whereas H ii galaxies are the least. LINERs show FR I-like core-brightened radio structures while Seyferts reveal the highest fraction of symmetric morphologies. The majority of H ii galaxies have single radio core or complex extended structures, which probably conceal a nuclear starburst and/or a weak active nucleus (seven of them show clear jets). ALGs, which are typically found in evolved ellipticals, although the least numerous, exhibit on average the most luminous radio structures, similar to LINERs. © 2020 The Author(s). Published by Oxford University Press on behalf of Royal Astronomical Society., The authors thank the anonymous referee for his/her helpful comments to improve the manuscript. AA and MAPT acknowledge support from the Spanish MCIU through grant PGC2018-098915-B-C21 and from the State Agency for Research of the Spanish MCIU through the 'Center of Excellence Severo Ochoa' award for the Instituto de Astrof ' isica de Andalucia (SEV-2017-0709). BTD acknowledges support from a Spanish postdoctoral fellowship 'Ayudas 1265 para la atraccion del talento investigador. Modalidad 2: jovenes investigadores.' funded by Comunidad de Madrid under grant number 2016-T2/TIC-2039. BTD also acknowledges support from grant 'Ayudas para la realizaci on de proyectos de I + D para jovenes doctores 2019.' funded by Comunidad de Madrid and Universidad Complutense de Madrid under grant number PR65/19-22417. JHK acknowledges financial support from the European Union's Horizon 2020 research and innovation programme under Marie Sklodowska-Curie grant agreement No. 721463 to the SUNDIAL ITN network, from the State Research Agency (AEI-MCINN) of the Spanish Ministry of Science and Innovation under the grant 'The structure and evolution of galaxies and their central regions' with reference PID2019-105602GB-I00/10.13039/501100011033, and from IAC project P/300724, financed by the Ministry of Science and Innovation, through the State Budget and by the Canary Islands Department of Economy, Knowledge and Employment, through the Regional Budget of the Autonomous Community. JSG thanks the University of Wisconsin-Madison and its Foundation for support of this research through his Rupple Bascom Professorship. FS acknowledges partial support from a Leverhulme Trust Research fellowship. CGM acknowledges support from the University of Bath and Jim and Hiroko Sherwin. e-MERLIN is a National Facility operated by the University of Manchester at Jodrell Bank Observatory on behalf of STFC, part of UK Research and Innovation.

Solar coronal loops are the building blocks of the solar corona. These dynamic structures are shaped by the magnetic field that expands into the solar atmosphere. They can be observed in X-ray and extreme ultraviolet (EUV), revealing the high plasma temperature of the corona. However, the dissipation of magnetic energy to heat the plasma to millions of degrees and, more generally, the mechanisms setting the mass and energy circulation in the solar atmosphere are still a matter of debate. Furthermore, multi-dimensional modelling indicates that the very concept of a coronal loop as an individual entity and its identification in EUV images is ill-defined due to the expected stochasticity of the solar atmosphere with continuous magnetic connectivity changes combined with the optically thin nature of the solar corona. In this context, the recent discovery of ubiquitous long-period EUV pulsations, the observed coronal rain properties and their common link in between represent not only major observational constraints for coronal heating theories but also major theoretical puzzles. The mechanisms of thermal non-equilibrium (TNE) and thermal instability (TI) appear in concert to explain these multi-scale phenomena as evaporation-condensation cycles. Recent numerical efforts clearly illustrate the specific but large parameter space involved in the heating and cooling aspects, and the geometry of the loop affecting the onset and properties of such cycles. In this review we will present and discuss this new approach into inferring coronal heating properties and understanding the mass and energy cycle based on the multi-scale intensity variability and cooling properties set by the TNE-TI scenario. We further discuss the major numerical challenges posed by the existence of TNE cycles and coronal rain, and similar phenomena at much larger scales in the Universe.

Tidally tilted pulsators (TTPs) are an intriguing new class of oscillating stars in binary systems; in such stars, the pulsation axis coincides with the line of apsides, or semi-major axis, of the binary. All three TTPs discovered so far have been $\delta$~Scuti stars. In this Letter, we report the first conclusive discovery of tidally tilted pulsations in a subdwarf B (sdB) star. HD 265435 is an sdB--white dwarf binary with a 1.65-hr period that has been identified and characterized as the nearest potential Type Ia supernova progenitor. Using TESS 20-s cadence data from Sectors 44 and 45, we show that the pulsation axis of the sdB star has been tidally tilted into the orbital plane and aligned with the tidal axis of the binary. We identify 31 independent pulsation frequencies, 27 of which have between 1 and 7 sidebands separated by the orbital frequency ($\nu_{\rm orb}$), or multiples thereof. Using the observed amplitude and phase variability due to tidal tilting, we assign $\ell$ and $m$ values to most of the observed oscillation modes and use these mode identifications to generate preliminary asteroseismic constraints. Our work significantly expands our understanding of TTPs, as we now know that (i) they can be found in stars other than $\delta$~Scuti pulsators, especially highly-evolved stars that have lost their H-rich envelopes, and (ii) tidally tilted pulsations can be used to probe the interiors of stars in very tight binaries., Comment: 17 pages, 5 figures + 2 appendix figures, 1 table; accepted for publication in ApJL

We present a study based on the high-resolution spectroscopy and K2 space photometry of five chemically peculiar stars in the region of the open cluster M44. The analysis of the high-precision photometric K2 data reveals that the light variations in HD 73045 and HD 76310 are rotational in nature and caused by spots or cloud-like co-rotating structures, which are non-stationary and short-lived. The time-resolved radial velocity measurements, in combination with the K2 photometry, confirm that HD 73045 does not show any periodic variability on timescales shorter than 1.3 d, contrary to previous reports in the literature. In addition to these new rotational variables, we discovered a new heartbeat system, HD 73619, where no pulsational signatures are seen. The spectroscopic and spectropolarimetric analyses indicate that HD 73619 belongs to the peculiar Am class, with either a weak or no magnetic field considering the 200 G detection limit of our study. The Least-Squares Deconvolution (LSD) profiles for HD 76310 indicate a complex structure in its spectra suggesting that this star is either part of a binary system or surrounded by a cloud shell. When placed in the Hertzsprung-Russell diagram, all studied stars are evolved from main-sequence and situated in the $\delta$ Scuti instability strip. The present work is relevant for further detailed studies of CP stars, such as inhomogeneities (including spots) in the absence of magnetic fields and the origin of the pulsational variability in heartbeat systems., Comment: Accepted for publication in MNRAS

Gaia DR2 has revealed breathing motions in the Milky Way, with stars on both sides of the Galactic mid-plane moving coherently towards or away from it. The generating mechanism of these breathing motions is thought to be spiral density waves. Here we test this hypothesis. Using a self-consistent, high-resolution simulation with star formation, and which hosts prominent spirals, we first study the signatures of breathing motions excited by spirals. In the model, the breathing motions induced by the spiral structure have an increasing amplitude with distance from the mid-plane, pointing to an internal cause for them. We then show that, at fixed height, the breathing motion amplitude decreases with age. Next, we investigate the signature of the breathing motions in the Gaia DR2 dataset. We demonstrate that, at the location with a consistently large breathing motion, the corresponding amplitude increases monotonically with distance from the mid-plane, in agreement with the model. Furthermore, we show that at the same location, the breathing motion amplitude decreases with age, again similar to what we find in the model. This strengthens the case that the observed breathing motions are driven by spiral density waves., 17 pages, 19 figures, accepted for publication in MNRAS

We present 850 $��$m polarimetric observations toward the Serpens Main molecular cloud obtained using the POL-2 polarimeter on the James Clerk Maxwell Telescope (JCMT) as part of the B-fields In STar-forming Region Observations (BISTRO) survey. These observations probe the magnetic field morphology of the Serpens Main molecular cloud on about 6000 au scales, which consists of cores and six filaments with different physical properties such as density and star formation activity. Using the histogram of relative orientation (HRO) technique, we find that magnetic fields are parallel to filaments in less dense filamentary structures where $N_{H_2} < 0.93\times 10^{22}$ cm$^{-2}$ (magnetic fields perpendicular to density gradients), while being perpendicular to filaments (magnetic fields parallel to density gradients) in dense filamentary structures with star formation activity. Moreover, applying the HRO technique to denser core regions, we find that magnetic field orientations change to become perpendicular to density gradients again at $N_{H_2} \approx 4.6 \times 10^{22}$ cm$^{-2}$. This can be interpreted as a signature of core formation. At $N_{H_2} \approx 16 \times 10^{22}$ cm$^{-2}$ magnetic fields change back to being parallel to density gradients once again, which can be understood to be due to magnetic fields being dragged in by infalling material. In addition, we estimate the magnetic field strengths of the filaments ($B_{POS} = 60-300~��$G)) using the Davis-Chandrasekhar-Fermi method and discuss whether the filaments are gravitationally unstable based on magnetic field and turbulence energy densities., 18 pages, accepted for publication in ApJ

Massive Early-Type Galaxies (ETG) in the local Universe are believed to be the most mature stage of galaxy evolution. Their stellar population content reveals the evolutionary history of these galaxies. However, while state-of-the-art Stellar Population Synthesis (SPS) models provide an accurate description of observed galaxy spectra in the optical range, the modelling in the Near-Infrared (NIR) is still in its infancy. Here we focus on NIR CO absorption features to show, in a systematic and comprehensive manner, that for massive ETGs, all CO indices, from H through to K band, are significantly stronger than currently predicted by SPS models. We explore and discuss several possible explanations of this "CO mismatch", including the effect of intermediate-age, AGB-dominated, stellar populations, high metallicity populations, non-solar abundance ratios and the initial mass function. While none of these effects is able to reconcile models and observations, we show that ad-hoc "empirical" corrections, taking into account the effect of CO-strong giant stars in the low-temperature regime, provide model predictions that are closer to the observations. Our analysis points to the effect of carbon abundance as the most likely explanation of NIR CO line-strengths, indicating possible routes for improving the SPS models in the NIR., 24 pages, 12 figures, accepted for publication in MNRAS

New spectrograms of multiply ionized iron have been recorded and analyzed, targeting the Fe vii spectrum. As a result, several previously unknown spectral lines and energy levels have been identified in this spectrum. These new data have been analyzed together with all previously published laboratory and astrophysical data on this spectrum. The energy levels have been interpreted using parametric calculations with Cowan codes. Radiative transition rates calculated in this work supplemented other previously published calculations in constructing a complete set of recommended transition probabilities. The ionization energy of Fe vii has been redetermined with a fivefold improvement in accuracy. Its new value is 1,007,928(20) cm−1, corresponding to 124.9671(25) eV.

Oscillatory reconnection (a relaxation mechanism with periodic changes in connectivity) has been proposed as a potential physical mechanism underpinning several periodic phenomena in the solar atmosphere including, but not limited to, quasi-periodic pulsations (QPPs). Despite its importance, however, the mechanism has never been studied within a hot, coronal plasma. We investigate oscillatory reconnection in a one million Kelvin plasma by solving the fully-compressive, resistive MHD equations for a 2D magnetic X-point under coronal conditions using the PLUTO code. We report on the resulting oscillatory reconnection including its periodicity and decay rate. We observe a more complicated oscillating profile for the current density compared to that found for a cold plasma, due to mode-conversion at the equipartition layer. We also consider, for the first time, the effect of adding anisotropic thermal conduction to the oscillatory reconnection mechanism, and we find this simplifies the spectrum of the oscillation profile and increases the decay rate. Crucially, the addition of thermal conduction does not prevent the oscillatory reconnection mechanism from manifesting. Finally, we reveal a relationship between the equilibrium magnetic field strength, decay rate, and period of oscillatory reconnection, which opens the tantalising possibility of utilizing oscillatory reconnection as a seismological tool., 17 pages, 13 figures, accepted for publication in ApJ

We analyze probability distributions of Birkeland current densities measured by the Active Magnetosphere and Planetary Electrodynamics Response Experiment (AMPERE). We find that the distributions are leptokurtic rather than normal and they are sometimes heavy-tailed. We fit q-exponential functions to the distributions and use these to estimate where the largest currents are likely to occur. The shape and scale parameters of the fitted q-exponential distribution vary with location: The scale parameter maximises for current densities with the same polarity and in the same location as the average Region 1 current, whereas the shape parameter maximises for current densities with the same polarity and in the same location as the average Region 2 current. We find that current densities |J|≥ 0.2 μA m−2 are most likely to occur in the average Region 1 current region, and second most likely to occur in the average Region 2 current region. However, for extreme currents (|J|≥ 4.0 μA m−2), we find that the most likely location is colocated with the average Region 2 current region on the dayside, at a colatitude of 18° − 22°.

We analyze the micro-kinetic stability of the electron strahl in the solar wind depending on heliocentric distance. The oblique fast-magnetosonic/whistler (FM/W) instability has emerged in the literature as a key candidate mechanism for the effective scattering of the electron strahl into the electron halo population. Using data from Parker Solar Probe (PSP) and Helios, we compare the measured strahl properties with the analytical thresholds for the oblique FM/W instability in the low- and high-$\beta_{\parallel c}$ regimes, where $\beta_{\parallel c}$ is the ratio of the core parallel thermal pressure to the magnetic pressure. Our PSP and Helios data show that the electron strahl is on average stable against the oblique FM/W instability in the inner heliosphere. Our analysis suggests that the instability, if at all, can only be excited sporadically and on short timescales. We discuss the caveats of our analysis and potential alternative explanations for the observed scattering of the electron strahl in the solar wind. Furthermore, we recommend the numerical evaluation of the stability of individual distributions in the future to account for any uncertainties in the validity of the analytical expressions for the instability thresholds., Comment: Accepted in ApJL

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted reuse, distribution, and reproduction in any medium, provided the original work is properly cited.--Full list of authors: Satapathy, Kaushik; Psaltis, Dimitrios; Özel, Feryal; Medeiros, Lia; Dougall, Sean T.; Chan, Chi-Kwan; Wielgus, Maciek; Prather, Ben S.; Wong, George N.; Gammie, Charles F.; Akiyama, Kazunori; Alberdi, Antxon; Alef, Walter; Algaba, Juan Carlos; Anantua, Richard; Asada, Keiichi; Azulay, Rebecca; Baczko, Anne-Kathrin; Ball, David; Baloković, Mislav; Barrett, John; Benson, Bradford A.; Bintley, Dan; Blackburn, Lindy; Blundell, Raymond; Boland, Wilfred; Bouman, Katherine L.; Bower, Geoffrey C.; Boyce, Hope; Bremer, Michael; Brinkerink, Christiaan D.; Brissenden, Roger; Britzen, Silke; Broderick, Avery E.; Broguiere, Dominique; Bronzwaer, Thomas; Bustamente, Sandra; Byun, Do-Young; Carlstrom, John E.; Chael, Andrew; Chatterjee, Koushik; Chatterjee, Shami; Chen, Ming-Tang; Chen, Yongjun; Cho, Ilje; Christian, Pierre; Conway, John E.; Cordes, James M.; Crawford, Thomas M.; Crew, Geoffrey B.; Cruz-Osorio, Alejandro; Cui, Yuzhu; Davelaar, Jordy; De Laurentis, Mariafelicia; Deane, Roger; Dempsey, Jessica; Desvignes, Gregory; Dexter, Jason; Doeleman, Sheperd S.; Eatough, Ralph P.; Falcke, Heino; Farah, Joseph; Fish, Vincent L.; Fomalont, Ed; Ford, H. Alyson; Fraga-Encinas, Raquel; Friberg, Per; Fromm, Christian M.; Fuentes, Antonio; Galison, Peter; García, Roberto; Gentaz, Olivier; Georgiev, Boris; Goddi, Ciriaco; Gold, Roman; Gómez-Ruiz, Arturo I.; Gómez, José L.; Gu, Minfeng; Gurwell, Mark; Hada, Kazuhiro; Haggard, Daryl; Hecht, Michael H.; Hesper, Ronald; Ho, Luis C.; Ho, Paul; Honma, Mareki; Huang, Chih-Wei L.; Huang, Lei; Hughes, David H.; Ikeda, Shiro; Inoue, Makoto; Issaoun, Sara; James, David J.; Jannuzi, Buell T.; Janssen, Michael; Jeter, Britton; Jiang, Wu; Jimenez-Rosales, Alejandra; Johnson, Michael D.; Jorstad, Svetlana; Jung, Taehyun; Karami, Mansour; Karuppusamy, Ramesh; Kawashima, Tomohisa; Keating, Garrett K.; Kettenis, Mark; Kim, Dong-Jin; Kim, Jae-Young; Kim, Jongsoo; Kim, Junhan; Kino, Motoki; Koay, Jun Yi; Kofuji, Yutaro; Koch, Patrick M.; Koyama, Shoko; Kramer, Carsten; Kramer, Michael; Krichbaum, Thomas P.; Kuo, Cheng-Yu; Lauer, Tod R.; Lee, Sang-Sung; Levis, Aviad; Li, Yan-Rong; Li, Zhiyuan; Lindqvist, Michael; Lico, Rocco; Lindahl, Greg; Liu, Jun; Liu, Kuo; Liuzzo, Elisabetta; Lo, Wen-Ping; Lobanov, Andrei P.; Loinard, Laurent; Lonsdale, Colin; Lu, Ru-Sen; MacDonald, Nicholas R.; Mao, Jirong; Marchili, Nicola; Markoff, Sera; Marrone, Daniel P.; Marscher, Alan P.; Martí-Vidal, Iván; Matsushita, Satoki; Matthews, Lynn D.; Menten, Karl M.; Mizuno, Izumi; Mizuno, Yosuke; Moran, James M.; Moriyama, Kotaro; Moscibrodzka, Monika; Müller, Cornelia; Mejías, Alejandro Mus; Musoke, Gibwa; Nagai, Hiroshi; Nagar, Neil M.; Nakamura, Masanori; Narayan, Ramesh; Narayanan, Gopal; Natarajan, Iniyan; Nathanail, Antonios; Neilsen, Joey; Neri, Roberto; Ni, Chunchong; Noutsos, Aristeidis; Nowak, Michael A.; Okino, Hiroki; Olivares, Héctor; Ortiz-León, Gisela N.; Oyama, Tomoaki; Palumbo, Daniel C. M.; Park, Jongho; Patel, Nimesh; Pen, Ue-Li; Pesce, Dominic W.; Piétu, Vincent; Plambeck, Richard; PopStefanija, Aleksandar; Porth, Oliver; Pötzl, Felix M.; Preciado-López, Jorge A.; Pu, Hung-Yi; Ramakrishnan, Venkatessh; Rao, Ramprasad; Rawlings, Mark G.; Raymond, Alexander W.; Rezzolla, Luciano; Ripperda, Bart; Roelofs, Freek; Rogers, Alan; Ros, Eduardo; Rose, Mel; Roshanineshat, Arash; Rottmann, Helge; Roy, Alan L.; Ruszczyk, Chet; Rygl, Kazi L. J.; Sánchez, Salvador; Sánchez-Arguelles, David; Sasada, Mahito; Savolainen, Tuomas; Schloerb, F. Peter; Schuster, Karl-Friedrich; Shao, Lijing; Shen, Zhiqiang; Small, Des; Sohn, Bong Won; SooHoo, Jason; Sun, He; Tazaki, Fumie; Tetarenko, Alexandra J.; Tiede, Paul; Tilanus, Remo P. J.; Titus, Michael; Toma, Kenji; Torne, Pablo; Traianou, Efthalia; Trent, Tyler; Trippe, Sascha; van Bemmel, Ilse; van Langevelde, Huib Jan; van Rossum, Daniel R.; Wagner, Jan; Ward-Thompson, Derek; Wardle, John; Weintroub, Jonathan; Wex, Norbert; Wharton, Robert; Wiik, Kaj; Wu, Qingwen; Yoon, Doosoo; Young, André; Young, Ken; Younsi, Ziri; Yuan, Feng; Yuan, Ye-Fei; Zensus, J. Anton; Zhao, Guang-Yao; Zhao, Shan-Shan., The black hole images obtained with the Event Horizon Telescope (EHT) are expected to be variable at the dynamical timescale near their horizons. For the black hole at the center of the M87 galaxy, this timescale (5–61 days) is comparable to the 6 day extent of the 2017 EHT observations. Closure phases along baseline triangles are robust interferometric observables that are sensitive to the expected structural changes of the images but are free of station-based atmospheric and instrumental errors. We explored the day-to-day variability in closure-phase measurements on all six linearly independent nontrivial baseline triangles that can be formed from the 2017 observations. We showed that three triangles exhibit very low day-to-day variability, with a dispersion of ∼3°–5°. The only triangles that exhibit substantially higher variability (∼90°–180°) are the ones with baselines that cross the visibility amplitude minima on the u–v plane, as expected from theoretical modeling. We used two sets of general relativistic magnetohydrodynamic simulations to explore the dependence of the predicted variability on various black hole and accretion-flow parameters. We found that changing the magnetic field configuration, electron temperature model, or black hole spin has a marginal effect on the model consistency with the observed level of variability. On the other hand, the most discriminating image characteristic of models is the fractional width of the bright ring of emission. Models that best reproduce the observed small level of variability are characterized by thin ring-like images with structures dominated by gravitational lensing effects and thus least affected by turbulence in the accreting plasmas. © 2022. The Author(s). Published by the American Astronomical Society., This work was supported, in part, by the NSF PIRE award 1743747 and NASA ATP award 80NSSC20K0521. L.M. acknowledges support from an NSF Astronomy and Astrophysics Postdoctoral Fellowship under award No. AST-1903847. M.W. acknowledges the support of the Black Hole Initiative at Harvard University, which is funded by grants from the John Templeton Foundation and the Gordon and Betty Moore Foundation to Harvard University. All ray tracing calculations for Set A were performed with the El Gato GPU cluster at the University of Arizona that is funded by NSF award 1228509. All analyses for Set B were performed on CyVerse, supported by NSF grants DBI-0735191, DBI-1265383, and DBI-1743442. The EHT Collaboration thanks the following organizations and programs: the Academy of Finland (projects 274477, 284495, 312496, and 315721); the Agencia Nacional de Investigación y Desarrollo, Chile via NCN19058 (TITANs) and Fondecyt 3190878; the Alexander von Humboldt Stiftung; an Alfred P. Sloan Research Fellowship; Allegro, the European ALMA Regional Centre node in the Netherlands, the NL astronomy research network NOVA, and the astronomy institutes of the University of Amsterdam, Leiden University, and Radboud University; the black hole Initiative at Harvard University, through a grant (60477) from the John Templeton Foundation; the China Scholarship Council; Consejo Nacional de Ciencia y Tecnología (Mexico, projects U0004-246083, U0004-259839, F0003-272050, M0037-279006, F0003-281692, 104497, 275201, and 263356); the Delaney Family via the Delaney Family John A. Wheeler Chair at Perimeter Institute; Dirección General de Asuntos del Personal Académico-Universidad Nacional Autónoma de México (projects IN112417 and IN112820); the European Research Council Synergy Grant "BlackHoleCam: Imaging the Event Horizon of Black Holes" (grant 610058); the Generalitat Valenciana postdoctoral grant APOSTD/2018/177 and GenT Program (project CIDEGENT/2018/021); MICINN Research Project PID2019-108995GB-C22; the Gordon and Betty Moore Foundation (grant GBMF-3561); the Istituto Nazionale di Fisica Nucleare sezione di Napoli, iniziative specifiche TEONGRAV; the International Max Planck Research School for Astronomy and Astrophysics at the Universities of Bonn and Cologne; Joint Princeton/Flatiron and Joint Columbia/Flatiron Postdoctoral Fellowships; research at the Flatiron Institute is supported by the Simons Foundation; the Japanese Government (Monbukagakusho: MEXT) Scholarship; the Japan Society for the Promotion of Science (JSPS) Grant-in-Aid for JSPS Research Fellowship (JP17J08829); the Key Research Program of Frontier Sciences, Chinese Academy of Sciences (CAS; grants QYZDJ-SSW-SLH057, QYZDJSSW-SYS008, and ZDBS-LY-SLH011); the Leverhulme Trust Early Career Research Fellowship; the Max-Planck-Gesellschaft (MPG); the Max-Planck Partner Group of the MPG and the CAS; the MEXT/JSPS KAKENHI (grants 18KK0090, JP18K13594, JP18K03656, JP18H03721, 18K03709, 18H01245, and 25120007); the Malaysian Fundamental Research Grant Scheme (FRGS). FRGS/1/2019/STG02/UM/02/6; the MIT International Science and Technology Initiatives Funds; the Ministry of Science and Technology (MOST) of Taiwan (105-2112-M-001-025-MY3, 106-2112-M-001-011, 106-2119- M-001-027, 107-2119-M-001-017, 107-2119-M-001-020, 107-2119-M-110-005, 108-2112-M-001-048, and 109-2124-M-001-005); the National Aeronautics and Space Administration (NASA; Fermi Guest Investigator grant 80NSSC20K1567, NASA Astrophysics Theory Program grant 80NSSC20K0527, and NASA NuSTAR award 80NSSC20K0645); the National Institute of Natural Sciences of Japan; the National Key Research and Development Program of China (grants 2016YFA0400704 and 2016YFA0400702); the National Science Foundation (NSF; grants AST-0096454, AST-0352953, AST-0521233, AST-0705062, AST-0905844, AST-0922984, AST-1126433, AST-1140030, DGE-1144085, AST-1207704, AST-1207730, AST-1207752, MRI-1228509, OPP-1248097, AST-1310896, AST-1555365, AST-1615796, AST-1715061, AST-1716327, AST-1903847, and AST-2034306); the Natural Science Foundation of China (grants 11573051, 11633006, 11650110427, 10625314, 11721303, 11725312, 11933007, 11991052, and 11991053); a fellowship of China Postdoctoral Science Foundation (2020M671266); the Natural Sciences and Engineering Research Council of Canada (NSERC; including a Discovery Grant and the NSERC Alexander Graham Bell Canada Graduate Scholarships-Doctoral Program); the National Youth Thousand Talents Program of China; the National Research Foundation of Korea (the Global PhD Fellowship Grant: grants NRF-2015H1A2A1033752 and 2015- R1D1A1A01056807; the Korea Research Fellowship Program: NRF-2015H1D3A1066561 and Basic Research Support Grant 2019R1F1A1059721); the Netherlands Organization for Scientific Research VICI award (grant 639.043.513) and Spinoza Prize SPI 78-409; the New Scientific Frontiers with Precision Radio Interferometry Fellowship awarded by the South African Radio Astronomy Observatory, which is a facility of the National Research Foundation, an agency of the Department of Science and Technology of South Africa; the Onsala Space Observatory (OSO) national infrastructure, for the provisioning of its facilities/observational support (OSO receives funding through the Swedish Research Council under grant 2017–00648) the Perimeter Institute for Theoretical Physics (research at Perimeter Institute is supported by the Government of Canada through the Department of Innovation, Science, and Economic Development and by the Province of Ontario through the Ministry of Research, Innovation, and Science); the Spanish Ministerio de Economía y Competitividad (grants PGC2018-098915-B-C21, AYA2016-80889-P, and PID2019-108995GB-C21); the State Agency for Research of the Spanish MCIU through the "Center of Excellence Severo Ochoa" award for the Instituto de Astrofísica de Andalucía (SEV-2017–0709); the Toray Science Foundation; the Consejería de Economía, Conocimiento, Empresas y Universidad of the Junta de Andalucía (grant P18-FR-1769), the Consejo Superior de Investigaciones Científicas (grant 2019AEP112); the US Department of Energy (USDOE) through the Los Alamos National Laboratory (operated by Triad National Security, LLC), for the National Nuclear Security Administration of the USDOE (Contract 89233218CNA000001); the European Union's Horizon 2020 research and innovation program under grant agreement No. 730562 RadioNet; ALMA North America Development Fund; the Academia Sinica; Chandra DD7-18089X and TM6-17006X; the GenT Program (Generalitat Valenciana) Project CIDEGENT/2018/021. This work used the Extreme Science and Engineering Discovery Environment (XSEDE), supported by NSF grant ACI-1548562, and CyVerse, supported by NSF grants DBI-0735191, DBI-1265383, and DBI-1743442. XSEDE Stampede2 resource at TACC was allocated through TG-AST170024 and TG-AST080026N. XSEDE JetStream resource at PTI and TACC was allocated through AST170028. The simulations were performed in part on the SuperMUC cluster at the LRZ in Garching, on the LOEWE cluster in CSC in Frankfurt, and on the HazelHen cluster at the HLRS in Stuttgart. This research was enabled, in part, by support provided by Compute Ontario (http://computeontario.ca), Calcul Quebec (http://www.calculquebec.ca), and Compute Canada (http://www.computecanada.ca). We thank the staff at the participating observatories, correlation centers, and institutions for their enthusiastic support. This paper makes use of the following ALMA data: ADS/JAO.ALMA#2016.1.01154.V. ALMA is a partnership of the European Southern Observatory (ESO; Europe, representing its member states), NSF, and the National Institutes of Natural Sciences of Japan, together with National Research Council (Canada), MOST (Taiwan), the Academia Sinica Institute of Astronomy and Astrophysics (ASIAA; Taiwan), and the Korea Astronomy and Space Science Institute (KASI; Republic of Korea), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by ESO, Associated Universities, Inc. (AUI)/NRAO, and the National Astronomical Observatory of Japan (NAOJ). The NRAO is a facility of the NSF operated under cooperative agreement by AUI. APEX is a collaboration between the Max-Planck-Institut für Radioastronomie (Germany), ESO, and the OSO (Sweden). The SMA is a joint project between the SAO and ASIAA and is funded by the Smithsonian Institution and the Academia Sinica. The JCMT is operated by the East Asian Observatory on behalf of the NAOJ, ASIAA, and KASI, as well as the Ministry of Finance of China, CAS, and the National Key R&D Program (No. 2017YFA0402700) of China. Additional funding support for the JCMT is provided by the Science and Technologies Facility Council (UK) and participating universities in the UK and Canada. The LMT is a project operated by the Instituto Nacional de Astrófisica, Óptica, y Electrónica (Mexico) and the University of Massachusetts at Amherst (USA). The IRAM 30 m telescope in Pico Veleta, Spain is operated by IRAM and supported by CNRS (Centre National de la Recherche Scientifique, France), MPG (Max-Planck Gesellschaft, Germany), and IGN (Instituto Geográfico Nacional, Spain). The SMT is operated by the Arizona Radio Observatory, a part of the Steward Observatory of the University of Arizona, with financial support of operations from the State of Arizona and financial support for instrumentation development from the NSF. The SPT is supported by the NSF through grant PLR-1248097. Partial support is also provided by the NSF Physics Frontier Center grant PHY-1125897 to the Kavli Institute of Cosmological Physics at the University of Chicago, the Kavli Foundation, and the Gordon and Betty Moore Foundation grant GBMF 947. The SPT hydrogen maser was provided on loan from the GLT, courtesy of ASIAA. The EHTC has received generous donations of FPGA chips from Xilinx Inc., under the Xilinx University Program. The EHTC has benefited from technology shared under an open-source license by the Collaboration for Astronomy Signal Processing and Electronics Research. The EHT project is grateful to T4Science and Microsemi for their assistance with hydrogen masers. This research has made use of NASA's Astrophysics Data System. We gratefully acknowledge the support provided by the extended staff of the ALMA, both from the inception of the ALMA Phasing Project through the observational campaigns of 2017 and 2018.

Funding: OA acknowledges financial support from the University of Exeter and from the United Kingdom Natural Environment Research Council (NERC) Independent Research Fellowship NE/V013963/1. OA and CW acknowledge financial support from the NERC Highlight Topic Grant NE/P017274/1 (Rad-Sat), and from United Kingdom Science and Technology Facilities Council (STFC) via Consolidated Grant ST/W000369/1. TE acknowledges financial support from an Early Career Fellowship, split jointly by the Leverhulme Trust (ECF2019-155) and the University of Leicester in the first instance (2019-21), but presently the University of Glasgow (2021-). TN acknowledges financial support from the STFC via Consolidated Grant ST/S000402/1. The University of Exeter cover the Open Access Publication Fee via a UKRI block grant. We derive weak turbulence and quasilinear models for relativistic charged particle dynamics in pitch-angle and energy space, due to interactions with electromagnetic waves propagating (anti-)parallel to a uniform background magnetic field. We use a Markovian approach that starts from the consideration of single particle motion in a prescribed electromagnetic field. This Markovian approach has a number of benefits, including: (i) the evident self-consistent relationship between a more general weak turbulence theory and the standard resonant diffusion quasilinear theory (as is commonly used in e.g. radiation belt and solar wind modelling); (ii) the general nature of the Fokker-Planck equation that can be derived without any prior assumptions regarding its form; (iii) the clear dependence of the form of the Fokker-Planck equation and the transport coefficients on given specific timescales. The quasilinear diffusion coefficients that we derive are not new in and of themselves, but this concise derivation and discussion of the weak turbulence and quasilinear theories using the Markovian framework is physically very instructive. The results presented herein form fundamental groundwork for future studies that consider phenomena for which some of the assumptions made in this manuscript may be relaxed. Publisher PDF

Funding: T. Elsden was supported by a Leverhulme Trust Early Career Fellowship (ECF-2019-155), the University of Leicester and the University of Glasgow. I. J. Rae was supported by NERC Grants NE/P017185/1 and NE/V002554/2 and STFC Grant ST/V006320/1. T. K. Yeoman and M. K. James were supported by STFC Grant ST/S000429/1. J. K. Sandhu was supported by NERC Grants NE/P017185/2 and NE/V002554/2. M-T Walach was supported by NERC Grant NE/T000937/1. Previous observational studies have shown that the natural Alfvén frequencies of geomagnetic field lines vary significantly over the course of a geomagnetic storm, decreasing by up to 50% from their quiet time values outside the plasmasphere. This was recently demonstrated statistically using ground magnetometer observations across 132 geomagnetic storm events (Wharton et al., 2020). This then brings into question where field line resonances (FLRs) will form in storm-time conditions relative to quiet times. With storm-time radiation belt dynamics depending heavily upon wave-particle interactions, understanding how FLR locations change over the course of a storm will have important implications for this area. Using 3D magnetohydrodynamic (MHD) simulations, we investigate how changes in the Alfvén frequency continuum of the Earth's dayside magnetosphere over the course of a geomagnetic storm affect the fast-Alfvén wave coupling. By setting the model Alfvén frequencies consistent with the observations, and permitting a modest change in the plasmapause/magnetopause locations consistent with storm-time behavior, we show that FLR locations can change substantially during storms. The combined effects of higher fast waveguide frequencies and lower Alfvén frequencies during storm main phases, act together to move the FLR locations radially inwards compared to quiet times. FLRs outside of the plasmasphere are moved radially inward by 1.7 Earth radii for the cases considered. Publisher PDF

In Galaxy And Mass Assembly Data Release 4 (GAMA DR4), we make available our full spectroscopic redshift sample. This includes 248682 galaxy spectra, and, in combination with earlier surveys, results in 330542 redshifts across five sky regions covering ~250deg^2. The redshift density, is the highest available over such a sustained area, has exceptionally high completeness (95 per cent to r_KIDS=19.65mag), and is well suited for the study of galaxy mergers, galaxy groups, and the low redshift (z, Comment: Accepted for publication in MNRAS. GAMA Data Release 4 is available at: http://www.gama-survey.org/dr4/

The outer atmosphere of the Sun is composed of plasma heated to temperatures well in excess of the visible surface. We investigate short cool and warm (, Accepted for publication in the Astrophysical Journal, 24 pages, 18 figures

Thermal instability is a fundamental process of astrophysical plasmas. It is expected to occur whenever the cooling is dominated by radiation and cannot be compensated for by heating. In this work, we conduct 2.5D radiation MHD simulations with the Bifrost code of an enhanced activity network in the solar atmosphere. Coronal loops are produced self-consistently, mainly through Joule heating, which is sufficiently stratified and symmetric to produce thermal nonequilibrium. During the cooling and driven by thermal instability, coronal rain is produced along the loops. Due to flux freezing, the catastrophic cooling process leading to a rain clump produces a local enhancement of the magnetic field, thereby generating a distinct magnetic strand within the loop up to a few Gauss stronger than the surrounding coronal field. These strands, which can be considered fundamental, are a few hundred kilometers in width, span most of the loop leg, and emit strongly in the UV and extreme UV (EUV), thereby establishing a link between the commonly seen rain strands in the visible spectrum with the observed EUV coronal strands at high resolution. The compression downstream leads to an increase in temperature that generates a plume-like structure, a strongly emitting spicule-like feature, and short-lived brightening in the UV during the rain impact, providing an explanation for similar phenomena seen with IRIS. Thermal instability and nonequilibrium can therefore be associated with localized and intermittent UV brightening in the transition region and chromosphere, as well as contribute to the characteristic filamentary morphology of the solar corona in the EUV.

Attenuation by dust severely impacts our ability to obtain unbiased observations of galaxies, especially as the amount and wavelength dependence of the attenuation varies with the stellar mass M-*, inclination i, and other galaxy properties. In this study, we used the attenuation - inclination models in ultraviolet, optical, and near-infrared bands designed by Tuffs and collaborators to investigate the average global dust properties in galaxies as a function of M-*, the stellar mass surface density mu(*), the star-formation rate SFR, the specific star-formation rate sSFR, the star-formation main-sequence offset dMS, and the star-formation rate surface density Sigma(SFR) at redshifts z similar to 0 and z similar to 0.7. We used star-forming galaxies from the Sloan Digital Sky Survey (similar to 20 000) and Galaxy And Mass Assembly (similar to 2000) to form our low-z sample at 0.04 < z < 0.1 and star-forming galaxies from Cosmological Evolution Survey (similar to 2000) for the sample at 0.6 < z < 0.8. We found that galaxies at z similar to 0.7 have a higher optical depth tau(f)(B) and clumpiness F than galaxies at z similar to 0. The increase in F hints that the stars of z similar to 0.7 galaxies are less likely to escape their birth cloud, which might indicate that the birth clouds are larger. We also found that tau(f)(B) increases with M-* and mu(*), independent of the sample and therefore redshift. We found no clear trends in tau(f)(B) or F with the SFR, which could imply that the dust mass distribution is independent of the SFR. In turn, this would imply that the balance of dust formation and destruction is independent of the SFR. Based on an analysis of the inclination dependence of the Balmer decrement, we found that reproducing the Balmer line emission requires not only a completely optically thick dust component associated with star-forming regions, as in the standard model, but an extra component of an optically thin dust within the birth clouds. This new component implies the existence of dust inside H II regions that attenuates the Balmer emission before it escapes through gaps in the birth cloud and we found it is more important in high-mass galaxies. These results will inform our understanding of dust formation and dust geometry in star-forming galaxies across redshift., Australian Research Council FT140101202 639.042.611, Netherlands Organization for Scientific Research (NWO), Scientific Exchanges visitor fellowship IZSEZO_202357, Swiss National Science Foundation (SNSF), European Commission

Recent observations have shown that besides the characteristic multimillion degree component, the corona also contains a large amount of cool material called coronal rain, whose clumps are 10–100 times cooler and denser than the surroundings and are often organized in larger events, termed showers. Thermal instability (TI) within a coronal loop in a state of thermal nonequilibrium (TNE) is the leading mechanism behind the formation of coronal rain but no investigation on showers exists to date. In this study, we conduct a morphological and thermodynamic multiwavelength study of coronal rain showers observed in an active region (AR) off-limb with IRIS and the Solar Dynamics Observatory, spanning chromospheric to transition region and coronal temperatures. Rain showers were found to be widespread across the AR over the 5.45 hr observing time, with an average length, width, and duration of 27.37 ± 11.95 Mm, 2.14 ± 0.74 Mm, and 35.22 ± 20.35 minutes, respectively. We find a good correspondence between showers and the cooling coronal structures consistent with the TNE–TI scenario, thereby properly identifying coronal loops in the “coronal veil”, including the strong expansion at low heights and an almost zero expansion in the corona. This agrees with previous work suggesting that the observed zero expansion in the EUV is due to specific cross-field temperature distribution. We estimate the total number of showers to be 155 ± 40, leading to a TNE volume of 4.56 ± [3.71] × 1028 cm3, i.e., on the same order of the AR volume. This suggests a prevalence of TNE over the AR indicating strongly stratified and high-frequency heating on average.

Polarized dust emission is a key tracer in the study of interstellar medium and of star formation. The observed polarization, however, is a product of magnetic field structure, dust grain properties and grain alignment efficiency, as well as their variations in the line of sight, making it difficult to interpret polarization unambiguously. The comparison of polarimetry at multiple wavelengths is a possible way of mitigating this problem. We use data from HAWC+/SOFIA and from SCUBA-2/POL-2 (from the BISTRO survey) to analyse the NGC 2071 molecular cloud at 154, 214 and 850 $\mu$m. The polarization angle changes significantly with wavelength over part of NGC 2071, suggesting a change in magnetic field morphology on the line of sight as each wavelength best traces different dust populations. Other possible explanations are the existence of more than one polarization mechanism in the cloud or scattering from very large grains. The observed change of polarization fraction with wavelength, and the 214-to-154 $\mu$m polarization ratio in particular, are difficult to reproduce with current dust models under the assumption of uniform alignment efficiency. We also show that the standard procedure of using monochromatic intensity as a proxy for column density may produce spurious results at HAWC+ wavelengths. Using both long-wavelength (POL-2, 850 $\mu$m) and short-wavelength (HAWC+, $\lesssim 200\, \mu$m) polarimetry is key in obtaining these results. This study clearly shows the importance of multi-wavelength polarimetry at submillimeter bands to understand the dust properties of molecular clouds and the relationship between magnetic field and star formation., Comment: Main article: 18 pages, 11 figures. Online supplemental material: 2 pages, 3 figures

A grand-challenge problem at the forefront of physics is to understand how energy is transported and transformed in plasmas. This fundamental research priority encapsulates the conversion of plasma-flow and electromagnetic energies into particle energy, either as heat or some other form of energisation. The smallest characteristic scales, at which electron dynamics determines the plasma behaviour, are the next frontier in space and astrophysical plasma research. The analysis of astrophysical processes at these scales lies at the heart of the field of electron-astrophysics. Electron scales are the ultimate bottleneck for dissipation of plasma turbulence, which is a fundamental process not understood in the electron-kinetic regime. Since electrons are the most numerous and most mobile plasma species in fully ionised plasmas and are strongly guided by the magnetic field, their thermal properties couple very efficiently to global plasma dynamics and thermodynamics., Comment: White Paper for the Voyage 2050 Long-Term Plan in the ESA Science Programme; 27 pages

In this editorial I will demonstrate with newly discovered solar activity proxy-magnetic field that the Sun has entered into the modern Grand Solar Minimum (2020–2053) that will lead to a significant reduction of solar magnetic field and activity like during Maunder minimum leading to noticeable reduction of terrestrial temperature.

In this paper, a V-shape photonic crystal fiber (PCF) polarization filter based on surface plasmon resonance (SPR) effect is proposed. With the full vector finite element method, the V-shape photonic crystal fiber polarization filter is designed, and the coupling characteristics between the core mode and surface plasmon polariton mode are analyzed. The simulation results show that the fiber parameters have significant effects on the shift and strength of the SPR wavelength. Moreover, the losses of the core modes along the Y and X polarization directions are 68904 dB/m and 858 dB/m at the SPR wavelength of 1550 nm, respectively. Finally, the extinction ratio and error-tolerant rate of extinction ratio are also analyzed. The proposed V-shape PCF polarization filter can achieve good filtering effect in the communication band and is easy to integrate with existing optical fiber communication and sensing systems.