Search SciRate

When the solar wind interacts with the ionosphere of an unmagnetized planet, it induces currents that form an induced magnetosphere. These currents and their associated magnetic fields play a pivotal role in controlling the movement of charged particles, which is essential for understanding the escape of planetary ions. Unlike the well-documented magnetospheric current systems, the ionospheric current systems on unmagnetized planets remain less understood, which constrains the quantification of electrodynamic energy transfer from stars to these planets. Here, utilizing eight years of data from the Mars Atmosphere and Volatile EvolutioN (MAVEN) mission, we investigate the global distribution of ionospheric currents on Mars. We have identified two distinct current systems in the ionosphere: one aligns with the solar wind electric field yet exhibits hemispheric asymmetry perpendicular to the electric field direction; the other corresponds to the flow pattern of annually-averaged neutral winds. We propose that these two current systems are driven by the solar wind and atmospheric neutral winds, respectively. Our findings reveal that Martian ionospheric dynamics are influenced by the neutral winds from below and the solar wind from above, highlighting the complex and intriguing nature of current systems on unmagnetized planets.

The star formation history (SFH) is a key issue in the evolution of galaxies. In this work, we developed a model based on a Gaussian and gamma function mixture to fit SFHs with varying numbers of components. Our primary objective was to use this model to reveal the shape of SFHs and the corresponding physical driving factors. Specifically, we applied this model to fit SFHs from the TNG100-1 simulation. Our study led to the following findings: 1) Our model fits with TNG star formation histories well, especially for high-mass and red galaxies; 2) A clear relationship exists between the number and shape of fitted components and the mass and color of galaxies, with notable differences observed between central/isolated and satellite galaxies. 3) Our model allowed us to extract different episodes of star formation within star formation histories with ease and analyze the duration and timing of each star formation episode. Our findings indicated a strong relationship between the timing of each star formation episode and galaxy mass and color.

The intergalactic medium (IGM) in the vicinity of galaxy protoclusters are interesting testbeds to study complex baryonic effects such as gravitational shocks and feedback. Here, we utilize hydrodynamical simulations from the SIMBA and The Three Hundred suites to study the mechanisms influencing large-scale Lyman-$\alpha$ transmission in $2<z<2.5$ protoclusters. We focus on the matter overdensity-Lyman-$\alpha$ transmission relation $(\delta_m-\delta_F)$ on Megaparsec-scales in these protoclusters, which is hypothesized to be sensitive to the feedback implementations. The lower-density regions represented by the SIMBA-100 cosmological volume trace the power-law $\delta_m-\delta_F$ relationship often known as the fluctuating Gunn-Peterson approximation. This trend is continued into higher-density regions covered by simulations that implement stellar feedback only. Simulations with AGN thermal and AGN jet feedback , however, exhibit progressively more Lyman-$\alpha$ transmission at fixed matter overdensity. Compared with the 7 protoclusters observed in the COSMOS field, only 2 display the excess absorption expected from protoclusters. The others exhibit deviations: 4 show some increased transparency suggested by AGN X-ray thermal feedback models while the highly transparent COSTCO-I protocluster appears to reflect intense jet feedback. Discrepancies with the stellar-feedback-only model suggests processes at play beyond gravitational heating and/or stellar feedback as the cause of the protocluster transparencies. Some form of AGN feedback is likely at play in the observed protoclusters, and possibly long-ranged AGN jets in the case of COSTCO-I. While more detailed and resolved simulations are required to move forward, our findings open new avenues for probing AGN feedback at Cosmic Noon.

We report a rare regime of Earth's magnetosphere interaction with sub-Alfvénic solar wind in which the windsock-like magnetosphere transforms into one with Alfvén wings. In the magnetic cloud of a Coronal Mass Ejection (CME) on April 24, 2023, NASA's Magnetospheric Multiscale mission distinguishes the following features: (1) unshocked and accelerated cold CME plasma coming directly against Earth's dayside magnetosphere; (2) dynamical wing filaments representing new channels of magnetic connection between the magnetosphere and foot points of the Sun's erupted flux rope; (3) cold CME ions observed with energized counter-streaming electrons, evidence of CME plasma captured due to reconnection between magnetic-cloud and Alfvén-wing field lines. The reported measurements advance our knowledge of CME interaction with planetary magnetospheres, and open new opportunities to understand how sub-Alfvénic plasma flows impact astrophysical bodies such as Mercury, moons of Jupiter, and exoplanets close to their host stars.

Large-scale characterization of exoplanetary atmospheres is on the horizon, thereby making it possible in the future to extract their statistical properties. In this context, by using a well validated model in the solar system, we carry out three-dimensional magnetohydrodynamic simulations to compute nonthermal atmospheric ion escape rates of unmagnetized rocky exoplanets as a function of their radius based on fixed stellar radiation and wind conditions. We find that the atmospheric escape rate is, unexpectedly and strikingly, a nonmonotonic function of the planetary radius $R$ and that it evinces a maximum at $R \sim 0.7\,R_\oplus$. This novel nonmonotonic behavior may arise from an intricate tradeoff between the cross-sectional area of a planet (which increases with size, boosting escape rates) and its associated escape velocity (which also increases with size, but diminishes escape rates). Our results could guide forthcoming observations because worlds with certain values of $R$ (such as $R \sim 0.7\,R_\oplus$) might exhibit comparatively higher escape rates when all other factors are constant.

Mars lacks a global magnetic field, and instead possesses small-scale crustal magnetic fields, making its magnetic environment fundamentally different from intrinsic magnetospheres like those of Earth or Saturn. Here we report the discovery of magnetospheric ion drift patterns, typical of intrinsic magnetospheres, at Mars usingmeasurements fromMarsAtmosphere and Volatile EvolutioNmission. Specifically, we observewedge-like dispersion structures of hydrogen ions exhibiting butterfly-shaped distributions within the Martian crustal fields, a feature previously observed only in planetary-scale intrinsic magnetospheres. These dispersed structures are the results of driftmotions that fundamentally resemble those observed in intrinsic magnetospheres. Our findings indicate that the Martian magnetosphere embodies an intermediate case where both the unmagnetized and magnetized ion behaviors could be observed because of the wide range of strengths and spatial scales of the crustal magnetic fields around Mars.

Ultra-cool dwarf stars are abundant, long-lived, and uniquely suited to enable the atmospheric study of transiting terrestrial companions with JWST. Amongst them, the most prominent is the M8.5V star TRAPPIST-1 and its seven planets. While JWST Cycle 1 observations have started to yield preliminary insights into the planets, they have also revealed that their atmospheric exploration requires a better understanding of their host star. Here, we propose a roadmap to characterize the TRAPPIST-1 system -- and others like it -- in an efficient and robust manner. We notably recommend that -- although more challenging to schedule -- multi-transit windows be prioritized to mitigate the effects of stellar activity and gather up to twice more transits per JWST hour spent. We conclude that, for such systems, planets cannot be studied in isolation by small programs, but rather need large-scale, jointly space- and ground-based initiatives to fully exploit the capabilities of JWST for the exploration of terrestrial planets.

Heliophysics is the field that "studies the nature of the Sun, and how it influences the very nature of space - and, in turn, the atmospheres of planetary bodies and the technology that exists there." However, NASA's Heliophysics Division tends to limit study of planetary magnetospheres and atmospheres to only those of Earth. This leaves exploration and understanding of space plasma physics at other worlds to the purview of the Planetary Science and Astrophysics Divisions. This is detrimental to the study of space plasma physics in general since, although some cross-divisional funding opportunities do exist, vital elements of space plasma physics can be best addressed by extending the expertise of Heliophysics scientists to other stellar and planetary magnetospheres. However, the diverse worlds within the solar system provide crucial environmental conditions that are not replicated at Earth but can provide deep insight into fundamental space plasma physics processes. Studying planetary systems with Heliophysics objectives, comprehensive instrumentation, and new grant opportunities for analysis and modeling would enable a novel understanding of fundamental and universal processes of space plasma physics. As such, the Heliophysics community should be prepared to consider, prioritize, and fund dedicated Heliophysics efforts to planetary targets to specifically study space physics and aeronomy objectives.

We report a $z=2.30$ galaxy protocluster (COSTCO-I) in the COSMOS field, where the Lyman-$\alpha$ forest as seen in the CLAMATO IGM tomography survey does not show significant absorption. This departs from the transmission-density relationship (often dubbed the fluctuating Gunn-Peterson approximation; FGPA) usually expected to hold at this epoch, which would lead one to predict strong Ly$\alpha$ absorption at the overdensity. For comparison, we generate mock Lyman-$\alpha$ forest maps by applying FGPA to constrained simulations of the COSMOS density field, and create mocks that incorporate the effects of finite sightline sampling, pixel noise, and Wiener filtering. Averaged over $r=15\,h^{-1}\,\mathrm{Mpc}$ around the protocluster, the observed Lyman-$\alpha$ forest is consistently more transparent in the real data than in the mocks, indicating a rejection of the null hypothesis that the gas in COSTCO-I follows FGPA ($p=0.0026$, or $2.79 \sigma$ significance). It suggests that the large-scale gas associated with COSTCO-I is being heated above the expectations of FGPA, which might be due to either large-scale AGN jet feedback or early gravitational shock heating. COSTCO-I is the first known large-scale region of the IGM that is observed to be transitioning from the optically-thin photoionized regime at Cosmic Noon, to eventually coalesce into an intra-cluster medium (ICM) by $z=0$. Future observations of similar structures will shed light on the growth of the ICM and allow constraints on AGN feedback mechanisms.

Collisionless magnetic reconnection typically requires kinetic treatments that are, in general, computationally expensive compared to fluid-based models. In this study, we use the magnetohydrodynamics with adaptively embedded particle-in-cell (MHD-AEPIC) model to study the interaction of two magnetic flux ropes. This innovative model embeds one or more adaptive PIC regions into a global MHD simulation domain such that the kinetic treatment is only applied in regions where kinetic physics is prominent. We compare the simulation results among three cases: 1) MHD with adaptively embedded PIC regions, 2) MHD with statically (or fixed) embedded PIC regions, and 3) a full PIC simulation. The comparison yields good agreement when analyzing their reconnection rates and magnetic island separations, as well as the ion pressure tensor elements and ion agyrotropy. In order to reach a good agreement among the three cases, large adaptive PIC regions are needed within the MHD domain, which indicates that the magnetic island coalescence problem is highly kinetic in nature where the coupling between the macro-scale MHD and micro-scale kinetic physics is important.

Space weather observations and modeling at Mars have begun but they must be significantly increased to support the future of Human Exploration on the Red Planet. A comprehensive space weather understanding of a planet without a global magnetosphere and a thin atmosphere is very different from our situation at Earth so there is substantial fundamental research remaining. It is expected that the development of suitable models will lead to a comprehensive operational Mars space weather alert (MSWA) system that would provide rapid dissemination of information to Earth controllers, astronauts in transit, and those in the exploration zone (EZ) on the surface by producing alerts that are delivered rapidly and are actionable. To illustrate the importance of such a system, we use a magnetohydrodynamic code to model an extreme Carrington-type coronal mass ejection (CME) event at Mars. The results show a significant induced surface field of nearly 3000 nT on the dayside that could radically affect unprotected electrical systems that would dramatically impact human survival on Mars. Other associated problems include coronal mass ejection (CME) shock-driven acceleration of solar energetic particles producing large doses of ionizing radiation at the Martian surface. In summary, along with working more closely with international partners, the next Heliophysics Decadal Survey must include a new initiative to meet expected demands for space weather forecasting in support of humans living and working on the surface of Mars. It will require significant effort to coordinate NASA and the international community contributions.

Kinetic approaches are generally accurate in dealing with microscale plasma physics problems but are computationally expensive for large-scale or multiscale systems. One of the long-standing problems in plasma physics is the integration of kinetic physics into fluid models, which is often achieved through sophisticated analytical closure terms. In this paper, we successfully construct a multi-moment fluid model with an implicit fluid closure included in the neural network using machine learning. The multi-moment fluid model is trained with a small fraction of sparsely sampled data from kinetic simulations of Landau damping, using the physics-informed neural network (PINN) and the gradient-enhanced physics-informed neural network (gPINN). The multi-moment fluid model constructed using either PINN or gPINN reproduces the time evolution of the electric field energy, including its damping rate, and the plasma dynamics from the kinetic simulations. In addition, we introduce a variant of the gPINN architecture, namely, gPINN$p$ to capture the Landau damping process. Instead of including the gradients of all the equation residuals, gPINN$p$ only adds the gradient of the pressure equation residual as one additional constraint. Among the three approaches, the gPINN$p$-constructed multi-moment fluid model offers the most accurate results. This work sheds light on the accurate and efficient modeling of large-scale systems, which can be extended to complex multiscale laboratory, space, and astrophysical plasma physics problems.

Magnetohydrodynamic turbulence regulates the transfer of energy from large to small scales in many astrophysical systems, including the solar atmosphere. We perform three-dimensional magnetohydrodynamic simulations with unprecedentedly large magnetic Reynolds number to reveal how rapid reconnection of magnetic field lines changes the classical paradigm of the turbulent energy cascade. By breaking elongated current sheets into chains of small magnetic flux ropes (or plasmoids), magnetic reconnection leads to a new range of turbulent energy cascade, where the rate of energy transfer is controlled by the growth rate of the plasmoids. As a consequence, the turbulent energy spectra steepen and attain a spectral index of -2.2 that is accompanied by changes in the anisotropy of turbulence eddies. The omnipresence of plasmoids and their consequences on, e.g., solar coronal heating, can be further explored with current and future spacecraft and telescopes.

Deriving governing equations of complex physical systems based on first principles can be quite challenging when there are certain unknown terms and hidden physical mechanisms in the systems. In this work, we apply a deep learning architecture to learn fluid partial differential equations (PDEs) of a plasma system based on the data acquired from a fully kinetic model. The learned multi-moment fluid PDEs are demonstrated to incorporate kinetic effects such as Landau damping. Based on the learned fluid closure, the data-driven, multi-moment fluid modeling can well reproduce all the physical quantities derived from the fully kinetic model. The calculated damping rate of Landau damping is consistent with both the fully kinetic simulation and the linear theory. The data-driven fluid modeling of PDEs for complex physical systems may be applied to improve fluid closure and reduce the computational cost of multi-scale modeling of global systems.

Heliophysics theory and modeling build understanding from fundamental principles to motivate, interpret, and predict observations. Together with observational analysis, they constitute a comprehensive scientific program in heliophysics. As observations and data analysis become increasingly detailed, it is critical that theory and modeling develop more quantitative predictions and iterate with observations. Advanced theory and modeling can inspire and greatly improve the design of new instruments and increase their chance of success. In addition, in order to build physics-based space weather forecast models, it is important to keep developing and testing new theories, and maintaining constant communications with theory and modeling. Maintaining a sustainable effort in theory and modeling is critically important to heliophysics. We recommend that all funding agencies join forces and consider expanding current and creating new theory and modeling programs--especially, 1. NASA should restore the HTMS program to its original support level to meet the critical needs of heliophysics science; 2. a Strategic Research Model program needs to be created to support model development for next-generation basic research codes; 3. new programs must be created for addressing mission-critical theory and modeling needs; and 4. enhanced programs are urgently required for training the next generation of theorists and modelers.

Magnetic reconnection is prevalent in magnetized plasmas in space and laboratories. Despite significant investigations on reconnection in electron-ion plasmas, studies of reconnection in magnetized plasmas with negatively charged dust grains are quite sparse. Here we report the first fully kinetic simulations of collisionless reconnection in a three-species (i.e., electron, proton, and negatively charged dust grain) dusty plasma, through which the discovery of double Hall pattern is made. The double Hall pattern consists of a traditional Hall quadruple current in between the ion and electron diffusion region, and a reversed Hall current in between the boundary of the ion and dust diffusion region. The analysis of the reconnection rate is also given. This study may be applicable to explain observations of planetary magnetospheres and the astrophysical objects, and may be realized in the laboratory studies of dusty plasmas.

We present results from end-to-end simulations of observations designed to constrain the rate of change in the expansion history of the Universe using the redshift drift of the Lyman-$\alpha$ forest absorption lines along the lines-of-sight toward bright quasars. For our simulations we take Lyman-$\alpha$ forest lines extracted from Keck/HIRES spectra of bright quasars at $z>3$, and compare the results from these real quasar spectra with mock spectra generated via Monte Carlo realizations. We use the results of these simulations to assess the potential for a dedicated observatory to detect redshift drift, and quantify the telescope and spectrograph requirements for these observations. Relative to Liske et al. (2008), two main refinements in the current work are inclusion of quasars from more recent catalogs and consideration of a realistic observing strategy for a dedicated redshift drift experiment that maximizes $\dot{v}/\sigma_{\dot{v}}$. We find that using a dedicated facility and our designed observing plan, the redshift drift can be detected at $3\sigma$ significance in 15 years with a 25m telescope, given a spectrograph with long term stability with $R=50,000$ and 25% total system efficiency. To achieve this significance, the optimal number of targets is four quasars, with observing time weighted based upon $\dot{v}/\sigma_{\dot{v}}$ and object visibility. This optimized strategy leads to a 9% decrease in the telescope diameter or a 6% decrease in the required time to achieve the same S/N as for the idealized case of uniformly distributing time to the same quasars.

The ESA-JAXA BepiColombo mission will provide simultaneous measurements from two spacecraft, offering an unprecedented opportunity to investigate magnetospheric and exospheric dynamics at Mercury as well as their interactions with the solar wind, radiation, and interplanetary dust. Many scientific instruments onboard the two spacecraft will be completely, or partially devoted to study the near-space environment of Mercury as well as the complex processes that govern it. Many issues remain unsolved even after the MESSENGER mission that ended in 2015. The specific orbits of the two spacecraft, MPO and Mio, and the comprehensive scientific payload allow a wider range of scientific questions to be addressed than those that could be achieved by the individual instruments acting alone, or by previous missions. These joint observations are of key importance because many phenomena in Mercury's environment are highly temporally and spatially variable. Examples of possible coordinated observations are described in this article, analysing the required geometrical conditions, pointing, resolutions and operation timing of different BepiColombo instruments sensors.

SPT0346-52 (z=5.7) is the most intensely star-forming galaxy discovered by the South Pole Telescope, with Sigma_SFR ~ 4200 Msol yr^-1 kpc^-2. In this paper, we expand on previous spatially-resolved studies, using ALMA observations of dust continuum, [NII]205 micron, [CII]158 micron, [OI]146 micron, and undetected [NII]122 micron and [OI]63 micron emission to study the multi-phase interstellar medium (ISM) in SPT0346-52. We use pixelated, visibility-based lens modeling to reconstruct the source-plane emission. We also model the source-plane emission using the photoionization code CLOUDY and find a supersolar metallicity system. We calculate T_dust = 48.3 K and lambda_peak = 80 micron, and see line deficits in all five lines. The ionized gas is less dense than comparable galaxies, with n_e < 32 cm^-3, while ~20% of the [CII]158 emission originates from the ionized phase of the ISM. We also calculate the masses of several phases of the ISM. We find that molecular gas dominates the mass of the ISM in SPT0346-52, with the molecular gas mass ~4x higher than the neutral atomic gas mass and ~100x higher than the ionized gas mass.

Study Analysis Group 21 (SAG21) of NASA's Exoplanet Exploration Program Analysis Group (ExoPAG) was organized to study the effect of stellar contamination on space-based transmission spectroscopy, a method for studying exoplanetary atmospheres by measuring the wavelength-dependent radius of a planet as it transits its star. Transmission spectroscopy relies on a precise understanding of the spectrum of the star being occulted. However, stars are not homogeneous, constant light sources but have temporally evolving photospheres and chromospheres with inhomogeneities like spots, faculae, plages, granules, and flares. This SAG brought together an interdisciplinary team of more than 100 scientists, with observers and theorists from the heliophysics, stellar astrophysics, planetary science, and exoplanetary atmosphere research communities, to study the current research needs that can be addressed in this context to make the most of transit studies from current NASA facilities like HST and JWST. The analysis produced 14 findings, which fall into three Science Themes encompassing (1) how the Sun is used as our best laboratory to calibrate our understanding of stellar heterogeneities ("The Sun as the Stellar Benchmark"), (2) how stars other than the Sun extend our knowledge of heterogeneities ("Surface Heterogeneities of Other Stars") and (3) how to incorporate information gathered for the Sun and other stars into transit studies ("Mapping Stellar Knowledge to Transit Studies"). In this invited review, we largely reproduce the final report of SAG21 as a contribution to the peer-reviewed literature.

Present and upcoming NASA missions will be intensively observing a selected, partially overlapping set of stars for exoplanet studies. Key physical and chemical information about these stars and their systems is needed for planning observations and interpreting the results. A target star archive of such data would benefit a wide cross-section of the exoplanet community by enhancing the chances of mission success and improving the efficiency of mission observatories. It would also provide a common, accessible resource for scientific analysis based on standardized assumptions, while revealing gaps or deficiencies in existing knowledge of stellar properties necessary for exoplanetary system characterization.

If humanity is ever to consider substantial, long-term colonization of Mars, the resources needed are going to be extensive. For a long-term human presence on Mars to be established, serious thought would need to be given to terraforming the planet. One major requirement for such terraforming is having the protection of a planetary magnetic field which Mars currently does not have. In this article we explore comprehensively for the first time, the practical and engineering challenges that affect the feasibility of creating an artificial magnetic field capable of encompassing Mars. This includes the concerns that define the design, where to locate the magnetic field generator and possible construction strategies. The rationale here is not to justify the need for a planetary magnetosphere but to put figures on the practicalities so as to be able to weigh the pros and cons of the different engineering approaches. The optimum solution proposed is completely novel, although inspired by natural situations and fusion plasma techniques. The solution with the lowest power, assembly and mass is to create an artificial charged particle ring (similar in form to a "radiation belt"), around the planet possibly formed by ejecting matter from one of the moons of Mars (in fashion similar to that that forms the Io-Jupiter plasma torus), but using electromagnetic and plasma waves to drive a net current in the ring(s) that results in an overall magnetic field. With a new era of space exploration underway, this is the time to start thinking about these new and bold future concepts and to begin filling strategic knowledge gaps. Furthermore, the principles explored here are also applicable to smaller scale objects like manned spacecraft, space stations or moon bases, which would benefit from the creation of protective mini-magnetospheres.

Exoplanets orbiting M-dwarfs within habitable zones are exposed to stellar environments more extreme than that terrestrial planets experience in our Solar System, which can significantly impact the atmospheres of the exoplanets and affect their habitability and sustainability. This study provides the first prediction of hot oxygen corona structure and the associated photochemical loss from a 1 bar CO2-dominated atmosphere of a Venus-like rocky exoplanet, where dissociative recombination of O2+ ions is assumed to be the major source reaction for the escape of neutral O atoms and formation of the hot O corona (or exospheres) as on Mars and Venus. We employ a 3D Monte Carlo code to simulate the exosphere of Proxima Centauri b (PCb) based on the ionosphere simulated by a 3D magnetohydrodynamic model. Our simulation results show that variability of the stellar wind dynamic pressure over one orbital period of PCb does not affect the overall spatial structure of the hot O corona but contributes to the change in the global hot O escape rate that varies by an order of magnitude. The escape increases dramatically when the planet possesses its intrinsic magnetic fields as the ionosphere becomes more extended with the presence of a global magnetic field. The extended hot O corona may lead to a more extended H exosphere through collisions between thermal H and hot O, which exemplifies the importance of considering nonthermal populations in exospheres to interpret future observations.

SPT0311-58 is the most massive infrared luminous system discovered so far during the Epoch of Reionization (EoR). In this paper, we present a detailed analysis of the molecular interstellar medium at z = 6.9, through high-resolution observations of the CO(6-5), CO(7-6), CO(10-9), [CI](2-1), and p-H2O(211-202) lines and dust continuum emission with the Atacama Large Millimeter/submillimeter Array (ALMA). The system consists of a pair of intensely star-forming gravitationally lensed galaxies (labelled West and East). The intrinsic far-infrared luminosity is (16 $\pm$ 4)$\times\rm 10^{12} \ \rm L_{\odot}$ in West and (27 $\pm$ 4)$\times\rm 10^{11} \ \rm L_{\odot}$ in East. We model the dust, CO, and [CI] using non-local thermodynamic equilibrium radiative transfer models and estimate the intrinsic gas mass to be (5.4 $\pm$ 3.4)$\times\rm 10^{11} \ \rm M_{\odot}$ in West and (3.1 $\pm$ 2.7)$\times\rm 10^{10} \ \rm M_{\odot}$ in East. We find that the CO spectral line energy distribution in West and East are typical of high-redshift sub-millimeter galaxies (SMGs). The CO-to-H2 conversion factor ($\alpha_{CO}$) and the gas depletion time scales estimated from the model are consistent with the high-redshift SMGs in the literature within the uncertainties. We find no evidence of evolution of depletion time with redshift in SMGs at z > 3. This is the most detailed study of molecular gas content of a galaxy in the EoR to-date, with the most distant detection of H2O in a galaxy without any evidence for active galactic nuclei in the literature.

Characterizing habitable exoplanets and/or their moons is of paramount importance. Here we show the results of our magnetic field topological modeling which demonstrate that terrestrial exoplanet-exomoon coupled magnetospheres work together to protect the early atmospheres of both the exoplanet and the exomoon. When exomoon magnetospheres are within the exoplanet's magnetospheric cavity, the exomoon magnetosphere acts like a protective magnetic bubble providing an additional magnetopause confronting the stellar winds when the moon is on the dayside. In addition, magnetic reconnection would create a critical pathway for the atmosphere exchange between the early exoplanet and exomoon. When the exomoon's magnetosphere is outside of the exoplanet's magnetosphere it then becomes the first line of defense against strong stellar winds, reducing the exoplanet's atmospheric loss to space. A brief discussion is given on how this type of exomoon would modify radio emissions from magnetized exoplanets.

Parker Solar Probe (PSP) aims at exploring the nascent solar wind close to the Sun. Meanwhile, PSP is also expected to encounter small objects like comets and asteroids. In this work, we survey the ephemerides to find a chance of recent encounter, and then model the interaction between released dusty plasmas and solar wind plasmas. On 2019 September 2, a comet-like object 322P/SOHO just passed its perihelion flying to a heliocentric distance of 0.12 au, and swept by PSP at a relative distance as close as 0.025 au. We present the dynamics of dust particles released from 322P, forming a curved dust tail. Along the PSP path in the simulated inner heliosphere, the states of plasma and magnetic field are sampled and illustrated, with the magnetic field sequences from simulation results being compared directly with the in-situ measurements from PSP. Through comparison, we suggest that 322P might be at a deficient activity level releasing limited dusty plasmas during its way to becoming a "rock comet". We also present images of solar wind streamers as recorded by WISPR, showing an indication of dust bombardment for the images superposed with messy trails. We observe from LASCO coronagraph that 322P was transiting from a dimming region to a relatively bright streamer during its perihelion passage, and simulate to confirm that 322P was flying from relatively faster to slower solar wind streams, modifying local plasma states of the streams.

Massive galaxy-scale outflows of gas are one of the most commonly-invoked mechanisms to regulate the growth and evolution of galaxies throughout the universe. While the gas in outflows spans a large range of temperatures and densities, the cold molecular phase is of particular interest because molecular outflows may be capable of suppressing star formation in galaxies by removing the star-forming gas. We have conducted the first survey of molecular outflows at z > 4, targeting 11 strongly-lensed dusty, star-forming galaxies (DSFGs) with high-resolution Atacama Large Millimeter Array (ALMA) observations of OH 119um absorption as an outflow tracer. In this first paper, we give an overview of the survey, focusing on the detection rate and structure of molecular outflows. We find unambiguous evidence for outflows in 8/11 (73%) galaxies, more than tripling the number known at z > 4. This implies that molecular winds in z > 4 DSFGs must have both a near-unity occurrence rate and large opening angles to be detectable in absorption. Lensing reconstructions reveal that 500pc-scale clumpy structures in the outflows are common. The individual clumps are not directly resolved, but from optical depth arguments we expect that future observations will require 50-200pc spatial resolution to do so. We do not detect high-velocity [CII] wings in any of the sources with clear OH outflows, indicating that [CII] is not a reliable tracer of molecular outflows. Our results represent a first step toward characterizing molecular outflows at z > 4 at the population level, demonstrating that large-scale outflows are ubiquitous among early massive, dusty galaxies.

Galactic outflows of molecular gas are a common occurrence in galaxies and may represent a mechanism by which galaxies self-regulate their growth, redistributing gas that could otherwise have formed stars. We previously presented the first survey of molecular outflows at z > 4 towards a sample of massive, dusty galaxies. Here we characterize the physical properties of the molecular outflows discovered in our survey. Using low-redshift outflows as a training set, we find agreement at the factor-of-two level between several outflow rate estimates. We find molecular outflow rates 150-800Msun/yr and infer mass loading factors just below unity. Among the high-redshift sources, the molecular mass loading factor shows no strong correlations with any other measured quantity. The outflow energetics are consistent with expectations for momentum-driven winds with star formation as the driving source, with no need for energy-conserving phases. There is no evidence for AGN activity in our sample, and while we cannot rule out deeply-buried AGN, their presence is not required to explain the outflow energetics, in contrast to nearby obscured galaxies with fast outflows. The fraction of the outflowing gas that will escape into the circumgalactic medium (CGM), though highly uncertain, may be as high as 50%. This nevertheless constitutes only a small fraction of the total cool CGM mass based on a comparison to z~2-3 quasar absorption line studies, but could represent >~10% of the CGM metal mass. Our survey offers the first statistical characterization of molecular outflow properties in the very early universe.

We present APEX-LABOCA 870 micron observations of the fields surrounding the nine brightest, high-redshift, unlensed objects discovered in the South Pole Telescope's (SPT) 2500 square degrees survey. Initially seen as point sources by SPT's 1-arcmin beam, the 19-arcsec resolution of our new data enables us to deblend these objects and search for submillimetre (submm) sources in the surrounding fields. We find a total of 98 sources above a threshold of 3.7 sigma in the observed area of 1300 square arcminutes, where the bright central cores resolve into multiple components. After applying a radial cut to our LABOCA sources to achieve uniform sensitivity and angular size across each of the nine fields, we compute the cumulative and differential number counts and compare them to estimates of the background, finding a significant overdensity of approximately 10 at 14 mJy. The large overdensities of bright submm sources surrounding these fields suggest that they could be candidate protoclusters undergoing massive star-formation events. Photometric and spectroscopic redshifts of the unlensed central objects range from 3 to 7, implying a volume density of star-forming protoclusters of approximately 0.1 per giga-parsec cube. If the surrounding submm sources in these fields are at the same redshifts as the central objects, then the total star-formation rates of these candidate protoclusters reach 10,000 solar masses per year, making them much more active at these redshifts than what has been seen so far in both simulations and observations.

Magnetic reconnection underlies many explosive phenomena in the heliosphere and in laboratory plasmas. The new research capabilities in theory/simulations, observations, and laboratory experiments provide the opportunity to solve the grand scientific challenges summarized in this whitepaper. Success will require enhanced and sustained investments from relevant funding agencies, increased interagency/international partnerships, and close collaborations of the solar, heliospheric, and laboratory plasma communities. These investments will deliver transformative progress in understanding magnetic reconnection and related explosive phenomena including space weather events.

We present Gemini-S and \it Spitzer-IRAC optical-through-near-IR observations in the field of the SPT2349-56 proto-cluster at $z=4.3$. We detect optical/IR counterparts for only nine of the 14 submillimetre galaxies (SMGs) previously identified by ALMA in the core of SPT2349-56. In addition, we detect four $z\sim4$ Lyman-break galaxies (LBGs) in the 30 arcsec diameter region surrounding this proto-cluster core. Three of the four LBGs are new systems, while one appears to be a counterpart of one of the nine observed SMGs. We identify a candidate brightest cluster galaxy (BCG) with a stellar mass of $(3.2^{+2.5}_{-1.4})\times10^{11}\,{\rm M}_{\odot}$. The stellar masses of the eight other SMGs place them on, above, and below the main sequence of star formation at $z\approx4.5$. The cumulative stellar mass for the SPT2349-56 core is at least $(11.5\pm2.9)\times10^{11}\,{\rm M}_{\odot}$, a sizeable fraction of the stellar mass in local BCGs, and close to the universal baryon fraction (0.16) relative to the virial mass of the core ($10^{13}\,{\rm M}_{\odot}$). As all 14 of these SMGs are destined to quickly merge, we conclude that the proto-cluster core has already developed a significant stellar mass at this early stage, comparable to $z=1$ BCGs. Importantly, we also find that the SPT2349-56 core structure would be difficult to uncover in optical surveys, with none of the ALMA sources being easily identifiable or constrained through $g,r,$ and $i$ colour-selection in deep optical surveys and only a modest overdensity of LBGs over the extended core structure. SPT2349-56 therefore represents a truly dust-obscured phase of a massive cluster core under formation.

The South Pole Telescope (SPT) has systematically identified 81 high-redshift, strongly gravitationally lensed, dusty star-forming galaxies (DSFGs) in a 2500 square degree cosmological mm-wave survey. We present the final spectroscopic redshift survey of this flux-limited ($S_{870\, \mathrm{\mu m}} > 25\, \mathrm{mJy}$) sample, initially selected at $1.4$ mm. The redshift survey was conducted with the Atacama Large Millimeter/submillimeter Array across the $3$ mm spectral window, targeting carbon monoxide line emission. By combining these measurements with ancillary data, the SPT sample is now spectroscopically complete, with redshifts spanning $1.9$$<$$z$$<$$6.9$ and a median of $z=3.9 \pm 0.2$. We present the mm through far-infrared photometry and spectral energy density fits for all sources, along with their inferred intrinsic properties. Comparing the properties of the SPT sources to the unlensed DSFG population, we demonstrate that the SPT-selected DSFGs represent the most extreme infrared-luminous galaxies, even after accounting for strong gravitational lensing. The SPT sources have a median star formation rate of $2.3(2)\times 10^3\, \mathrm{M_\odot yr^{-1}}$ and a median dust mass of $1.4(1)\times10^9\, \mathrm{M_\odot}$. However, the inferred gas depletion timescales of the SPT sources are comparable to those of unlensed DSFGs, once redshift is taken into account. This SPT sample contains roughly half of the known spectroscopically confirmed DSFGs at $z$$>$$5$, making this the largest sample of high-redshift DSFGs to-date, and enabling the "high-redshift tail" of extremely luminous DSFGs to be measured. Though galaxy formation models struggle to account for the SPT redshift distribution, the larger sample statistics from this complete and well-defined survey will help inform future theoretical efforts.

The recent discovery of an Earth-sized planet (TOI-700 d) in the habitable zone of an early-type M-dwarf by the Transiting Exoplanet Survey Satellite constitutes an important advance. In this Letter, we assess the feasibility of this planet to retain an atmosphere -- one of the chief ingredients for surface habitability -- over long timescales by employing state-of-the-art magnetohydrodynamic models to simulate the stellar wind and the associated rates of atmospheric escape. We take two major factors into consideration, namely, the planetary atmospheric composition and magnetic field. In all cases, we determine that the atmospheric ion escape rates are potentially a few orders of magnitude higher than the inner Solar system planets, but TOI-700 d is nevertheless capable of retaining a $1$ bar atmosphere over gigayear timescales for certain regions of the parameter space. The simulations show that the unmagnetized TOI-700 d with a 1 bar Earth-like atmosphere could be stripped away rather quickly ($<$ 1 gigayear), while the unmagnetized TOI-700 d with a 1 bar CO$_2$-dominated atmosphere could persist for many billions of years; we find that the magnetized Earth-like case falls in between these two scenarios. We also discuss the prospects for detecting radio emission of the planet (thereby constraining its magnetic field) and discerning the presence of an atmosphere.

This white paper summarizes major scientific challenges and opportunities in understanding magnetic reconnection and related explosive phenomena as a fundamental plasma process.

We present an extensive ALMA spectroscopic follow-up programme of the $z\,{=}\,4.3$ structure SPT2349$-$56, one of the most actively star-forming proto-cluster cores known, to identify additional members using their [C\sc ii] 158\,$\mu$m and \mboxCO(4--3) lines. In addition to robustly detecting the 14 previously published galaxies in this structure, we identify a further 15 associated galaxies at $z\,{=}\,4.3$, resolving 55$\,{\pm}\,$5\u2009per cent of the 870-$\mu$m flux density at 0.5\u2009arcsec resolution compared to 21\u2009arcsec single-dish data. These galaxies are distributed into a central core containing 23 galaxies extending out to 300\u2009kpc in diameter, and a northern extension, offset from the core by 400\u2009kpc, containing three galaxies. We discovered three additional galaxies in a red \it Herschel\/-SPIRE source 1.5\u2009Mpc from the main structure, suggesting the existence of many other sources at the same redshift as SPT2349$-$56 that are not yet detected in the limited coverage of our data. An analysis of the velocity distribution of the central galaxies indicates that this region may be virialized with a mass of (9$\pm$5)$\,{\times}\,$10$^{12}$\u2009M$_{\odot}$, while the two offset galaxy groups are about 30 and 60\u2009per cent less massive and show significant velocity offsets from the central group. We calculate the [C\sc ii] and far-infrared number counts, and find evidence for a break in the [C\sc ii] luminosity function. We estimate the average SFR density within the region of SPT2349$-$56 containing single-dish emission (a proper diametre of 720\u2009kpc), assuming spherical symmetry, to be roughly 4$\,{\times}\,10^4$\u2009M$_{\odot}$\u2009yr$^{-1}$\u2009Mpc$^{-3}$; this may be an order of magnitude greater than the most extreme examples seen in simulations.

Fundamental coding and software development skills are increasingly necessary for success in nearly every aspect of astronomical and astrophysical research as large surveys and high resolution simulations become the norm. However, professional training in these skills is inaccessible or impractical for many members of our community. Students and professionals alike have been expected to acquire these skills on their own, apart from formal classroom curriculum or on-the-job training. Despite the recognized importance of these skills, there is little opportunity to develop them - even for interested researchers. To ensure a workforce capable of taking advantage of the computational resources and the large volumes of data coming in the next decade, we must identify and support ways to make software development training widely accessible to community members, regardless of affiliation or career level. To develop and sustain a technology capable astronomical and astrophysical workforce, we recommend that agencies make funding and other resources available in order to encourage, support and, in some cases, require progress on necessary training, infrastructure and policies. In this white paper, we focus on recommendations for how funding agencies can lead in the promotion of activities to support the astronomy and astrophysical workforce in the 2020s.

We present a three-species (H$^+$, O$^+$ and e$^-$) multi-fluid magnetohydrodynamic (MHD) model, endowed with the requisite upper atmospheric chemistry, that is capable of accurately quantifying the magnitude of oxygen ion losses from "Earth-like" exoplanets in habitable zones, whose magnetic and rotational axes are roughly coincidental with one another. We apply this model to investigate the role of planetary obliquity in regulating atmospheric losses from a magnetic perspective. For Earth-like exoplanets orbiting solar-type stars, we demonstrate that the dependence of the total atmospheric ion loss rate on the planetary (magnetic) obliquity is relatively weak; the escape rates are found to vary between $2.19 \times 10^{26}$ s$^{-1}$ to $2.37 \times 10^{26}$ s$^{-1}$. In contrast, the obliquity can influence the atmospheric escape rate ($\sim$ $10^{28}$ s$^{-1}$) by more than a factor of $2$ (or $200\%$) in the case of Earth-like exoplanets orbiting late-type M-dwarfs. Thus, our simulations indicate that planetary obliquity may play a weak-to-moderate role insofar as the retention of an atmosphere (necessary for surface habitability) is concerned.

Software is a critical part of modern research, and yet there are insufficient mechanisms in the scholarly ecosystem to acknowledge, cite, and measure the impact of research software. The majority of academic fields rely on a one-dimensional credit model whereby academic articles (and their associated citations) are the dominant factor in the success of a researcher's career. In the petabyte era of astronomical science, citing software and measuring its impact enables academia to retain and reward researchers that make significant software contributions. These highly skilled researchers must be retained to maximize the scientific return from petabyte-scale datasets. Evolving beyond the one-dimensional credit model requires overcoming several key challenges, including the current scholarly ecosystem and scientific culture issues. This white paper will present these challenges and suggest practical solutions for elevating the role of software as a product of the research enterprise.

Water ($\rm H_{2}O$), one of the most ubiquitous molecules in the universe, has bright millimeter-wave emission lines easily observed at high-redshift with the current generation of instruments. The low excitation transition of $\rm H_{2}O$, p$-$$\rm H_{2}O$(202 $-$ 111) ($\nu_{rest}$ = 987.927 GHz) is known to trace the far-infrared (FIR) radiation field independent of the presence of active galactic nuclei (AGN) over many orders-of-magnitude in FIR luminosity (L$_{\rm FIR}$). This indicates that this transition arises mainly due to star formation. In this paper, we present spatially ($\sim$0.5 arcsec corresponding to $\sim$1 kiloparsec) and spectrally resolved ($\sim$100 kms$^{-1}$) observations of p$-$$\rm H_{2}O$(202 $-$ 111) in a sample of four strong gravitationally lensed high-redshift galaxies with the Atacama Large Millimeter/submillimeter Array (ALMA). In addition to increasing the sample of luminous ($ > $ $10^{12}$L$_{\odot}$) galaxies observed with $\rm H_{2}O$, this paper examines the L$_{\rm H_{2}O}$/L$_{\rm FIR}$ relation on resolved scales for the first time at high-redshift. We find that L$_{\rm H_{2}O}$ is correlated with L$_{\rm FIR}$ on both global and resolved kiloparsec scales within the galaxy in starbursts and AGN with average L$_{\rm H_{2}O}$/L$_{\rm FIR}$ =$2.76^{+2.15}_{-1.21}\times10^{-5}$. We find that the scatter in the observed L$_{\rm H_{2}O}$/L$_{\rm FIR}$ relation does not obviously correlate with the effective temperature of the dust spectral energy distribution (SED) or the molecular gas surface density. This is a first step in developing p$-$$\rm H_{2}O$(202 $-$ 111) as a resolved star formation rate (SFR) calibrator.

We present Atacama Compact Array and Atacama Pathfinder Experiment observations of the [N II] 205 $\mu$m fine-structure line in 40 sub-millimetre galaxies lying at redshifts z = 3 to 6, drawn from the 2500 deg$^2$ South Pole Telescope survey. This represents the largest uniformly selected sample of high-redshift [N II] 205 $\mu$m measurements to date. 29 sources also have [C II] 158 $\mu$m line observations allowing a characterization of the distribution of the [C II] to [N II] luminosity ratio for the first time at high-redshift. The sample exhibits a median L$_{[C II]}$ /L$_{[N II]}$ $\approx$ 11 and interquartile range of 5.0 to 24.7. These ratios are similar to those observed in local (U)LIRGs, possibly indicating similarities in their interstellar medium. At the extremes, we find individual sub-millimetre galaxies with L$_{[C II]}$ /L$_{[N II]}$ low enough to suggest a smaller contribution from neutral gas than ionized gas to the [C II] flux and high enough to suggest strongly photon or X-ray region dominated flux. These results highlight a large range in this line luminosity ratio for sub-millimetre galaxies, which may be caused by variations in gas density, the relative abundances of carbon and nitrogen, ionization parameter, metallicity, and a variation in the fractional abundance of ionized and neutral interstellar medium.

Aims: We present and study spatially resolved imaging obtained with the Atacama Large Millimeter/submillimeter Array (ALMA) of multiple $^{12}$CO($J=$6$-$5, 8$-$7 and 9$-$8) and two H$_2$O(2$_{02}-$1$_{11}$ and 2$_{11}-$2$_{02}$) emission lines and cold dust continuum toward the gravitationally lensed dusty star forming galaxy SPT0346-52 at z=$5.656$. Methods: Using a visibility-domain source-plane reconstruction we probe the structure and dynamics of the different components of the interstellar medium (ISM) in this galaxy down to scales of 1 kpc in the source plane. Results: Measurements of the intrinsic sizes of the different CO emission lines indicate that the higher J transitions trace more compact regions in the galaxy. Similarly, we find smaller dust continuum intrinsic sizes with decreasing wavelength, based on observations at rest-frame 130, 300 and 450$\mu$m. The source shows significant velocity structure, and clear asymmetry where an elongated structure is observed in the source plane with significant variations in their reconstructed sizes. This could be attributed to a compact merger or turbulent disk rotation. The differences in velocity structure through the different line tracers, however, hint at the former scenario in agreement with previous [CII] line imaging results. Measurements of the CO line ratios and magnifications yield significant variations as a function of velocity, suggesting that modeling of the ISM using integrated values could be misinterpreted. Modeling of the ISM in SPT0346-52 based on delensed fluxes indicate a highly dense and warm medium, qualitatively similar to that observed in high redshift quasar hosts.

We are now on a clear trajectory for improvements in exoplanet observations that will revolutionize our ability to characterize their atmospheric structure, composition, and circulation, from gas giants to rocky planets. However, exoplanet atmospheric models capable of interpreting the upcoming observations are often limited by insufficiencies in the laboratory and theoretical data that serve as critical inputs to atmospheric physical and chemical tools. Here we provide an up-to-date and condensed description of areas where laboratory and/or ab initio investigations could fill critical gaps in our ability to model exoplanet atmospheric opacities, clouds, and chemistry, building off a larger 2016 white paper, and endorsed by the NAS Exoplanet Science Strategy report. Now is the ideal time for progress in these areas, but this progress requires better access to, understanding of, and training in the production of spectroscopic data as well as a better insight into chemical reaction kinetics both thermal and radiation-induced at a broad range of temperatures. Given that most published efforts have emphasized relatively Earth-like conditions, we can expect significant and enlightening discoveries as emphasis moves to the exotic atmospheres of exoplanets.

Over the last two decades, the discovery of exoplanets has fundamentally changed our perception of the universe and humanity's place within it. Recent work indicates that a solar system's X-ray and high energy particle environment is of fundamental importance to the formation and development of the atmospheres of close-in planets such as hot Jupiters, and Earth-like planets around M stars. X-ray imaging and spectroscopy provide powerful and unique windows into the high energy flux that an exoplanet experiences, and X-ray photons also serve as proxies for potentially transfigurative coronal mass ejections. Finally, if the host star is a bright enough X-ray source, transit measurements akin to those in the optical and infrared are possible and allow for direct characterization of the upper atmospheres of exoplanets. In this brief white paper, we discuss contributions to the study of exoplanets and their environs which can be made by X-ray data of increasingly high quality that are achievable in the next 10--15 years.

For the first time, we explore the tightly coupled interior-magnetosphere system of Mercury by employing a three-dimensional ten-moment multifluid model. This novel fluid model incorporates the non-ideal effects including the Hall effect, inertia, and tensorial pressures that are critical for collisionless magnetic reconnection; therefore, it is particularly well suited for investigating $collisionless$ magnetic reconnection in Mercury's magnetotail and at the planet's magnetopause. The model is able to reproduce the observed magnetic field vectors, field-aligned currents, and cross-tail current sheet asymmetry (beyond the MHD approach) and the simulation results are in good agreement with spacecraft observations. We also study the magnetospheric response of Mercury to a hypothetical extreme event with an enhanced solar wind dynamic pressure, which demonstrates the significance of induction effects resulting from the electromagnetically-coupled interior. More interestingly, plasmoids (or flux ropes) are formed in Mercury's magnetotail during the event, indicating the highly dynamic nature of Mercury's magnetosphere.

Nearly a century after the discovery that we live in an expanding Universe, and two decades after the discovery of accelerating cosmic expansion, there remains no direct detection of this acceleration via redshift drift - a change in the cosmological expansion velocity versus time. Because cosmological redshift drift directly determines the Hubble parameter H(z), it is arguably the cleanest possible measurement of the expansion history, and has the potential to constrain dark energy models (e.g. Kim et al. 2015). The challenge is that the signal is small - the best observational constraint presently has an uncertainty several orders of magnitude larger than the expected signal (Darling 2012). Nonetheless, direct detection of redshift drift is becoming feasible, with upcoming facilities such as the ESO-ELT and SKA projecting possible detection within two to three decades. This timescale is uncomfortably long given the potential of this cosmological test. With dedicated experiments it should be possible to rapidly accelerate progress and detect redshift drift with only a five-year observational baseline. Such a facility would also be ideal for precision radial velocity measurements of exoplanets, which could be obtained as a byproduct of the ongoing calibration measurements for the experiment.

It is now recognized that energetic stellar photon and particle radiation evaporates and erodes planetary atmospheres and controls upper atmospheric chemistry. Key exoplanet host stars will be too faint at X-ray wavelengths for accurate characterization using existing generation and future slated X-ray telescopes. Observation of stellar coronal mass ejections and winds are also beyond current instrumentation. In line with theCommittee on an Exoplanet Science Strategy recognition that holistic observational approaches are needed, we point out here that a full understanding of exoplanet atmospheres, their evolution and determination of habitability requires a powerful high-resolution X-ray imaging and spectroscopic observatory. This is the only capability that can: (1) characterize by proxy the crucial, difficult to observe, EUV stellar flux, its history and its variations for planet hosting stars; (2) observe the stellar wind; (3) detect the subtle Doppler signatures of coronal mass ejections.

The field of exoplanetary science is making rapid progress both in statistical studies of exoplanet properties as well as in individual characterization. As space missions provide an emerging picture of formation and evolution of exoplanetary systems, the search for habitable worlds becomes one of the fundamental issues to address. To tackle such a complex challenge, we need to specify the conditions favorable for the origin, development and sustainment of life as we know it. This requires the understanding of global (astrospheric) and local (atmospheric, surface and internal) environments of exoplanets in the framework of the physical processes of the interaction between evolving planet-hosting stars along with exoplanetary evolution over geological timescales, and the resulting impact on climate and habitability of exoplanets. Feedbacks between astrophysical, physico-chemical atmospheric and geological processes can only be understood through interdisciplinary studies with the incorporation of progress in heliophysics, astrophysics, planetary, Earth sciences, astrobiology, and the origin of life communities. The assessment of the impacts of host stars on the climate and habitability of terrestrial (exo)planets and potential exomoons around them may significantly modify the extent and the location of the habitable zone and provide new directions for searching for signatures of life. Thus, characterization of stellar ionizing outputs becomes an important task for further understanding the extent of habitability in the universe. The goal of this white paper is to identify and describe promising key research goals to aid the theoretical characterization and observational detection of ionizing radiation from quiescent and flaring upper atmospheres of planet hosts as well as properties of stellar coronal mass ejections and stellar energetic particle events.

Transmission spectra probe the atmospheres of transiting exoplanets, but these observations are also subject to signals introduced by magnetic active regions on host stars. Here we outline scientific opportunities in the next decade for providing useful constraints on stellar photospheres and inform interpretations of transmission spectra of the smallest ($R<4\,R_{\odot}$) exoplanets. We identify and discuss four primary opportunities: (1) refining stellar magnetic active region properties through exoplanet crossing events; (2) spectral decomposition of active exoplanet host stars; (3) joint retrievals of stellar photospheric and planetary atmospheric properties with studies of transmission spectra; and (4) continued visual transmission spectroscopy studies to complement longer-wavelength studies from $\textit{JWST}$. We make five recommendations to the Astro2020 Decadal Survey Committee: (1) identify the transit light source (TLS) effect as a challenge to precise exoplanet transmission spectroscopy and an opportunity ripe for scientific advancement in the coming decade; (2) include characterization of host star photospheric heterogeneity as part of a comprehensive research strategy for studying transiting exoplanets; (3) support the construction of ground-based extremely large telescopes (ELTs); (4) support multi-disciplinary research teams that bring together the heliophysics, stellar physics, and exoplanet communities to further exploit transiting exoplanets as spatial probes of stellar photospheres; and (5) support visual transmission spectroscopy efforts as complements to longer-wavelength observational campaigns with $\textit{JWST}$.

The processes that transform gas and dust in circumstellar disks into diverse exoplanets remain poorly understood. One key pathway is to study exoplanets as they form in their young ($\sim$few~Myr) natal disks. Extremely Large Telescopes (ELTs) such as GMT, TMT, or ELT, can be used to establish the initial chemical conditions, locations, and timescales of planet formation, via (1)~measuring the physical and chemical conditions in protoplanetary disks using infrared spectroscopy and (2)~studying planet-disk interactions using imaging and spectro-astrometry. Our current knowledge is based on a limited sample of targets, representing the brightest, most extreme cases, and thus almost certainly represents an incomplete understanding. ELTs will play a transformational role in this arena, thanks to the high spatial and spectral resolution data they will deliver. We recommend a key science program to conduct a volume-limited survey of high-resolution spectroscopy and high-contrast imaging of the nearest protoplanetary disks that would result in an unbiased, holistic picture of planet formation as it occurs.

We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of high-J CO lines ($J_\mathrm{up}=6$, 7, 8) and associated dust continuum towards five strongly lensed, dusty, star-forming galaxies (DSFGs) at redshift $z = 2.7$-5.7. These galaxies, discovered in the South Pole Telescope survey, are observed at $0.2''$-$0.4''$ resolution with ALMA. Our high-resolution imaging coupled with the lensing magnification provides a measurement of the structure and kinematics of molecular gas in the background galaxies with spatial resolutions down to kiloparsec scales. We derive visibility-based lens models for each galaxy, accurately reproducing observations of four of the galaxies. Of these four targets, three show clear velocity gradients, of which two are likely rotating disks. We find that the reconstructed region of CO emission is less concentrated than the region emitting dust continuum even for the moderate-excitation CO lines, similar to what has been seen in the literature for lower-excitation transitions. We find that the lensing magnification of a given source can vary by 20-50% across the line profile, between the continuum and line, and between different CO transitions. We apply Large Velocity Gradient (LVG) modeling using apparent and intrinsic line ratios between lower-J and high-J CO lines. Ignoring these magnification variations can bias the estimate of physical properties of interstellar medium of the galaxies. The magnitude of the bias varies from galaxy to galaxy and is not necessarily predictable without high resolution observations.

SPT0346-52 is one of the most most luminous and intensely star-forming galaxies in the universe, with L_FIR > 10^13 L_sol and Sigma_SFR ~ 4200 M_sol yr^-1 kpc^-2. In this paper, we present ~0.15'' ALMA observations of the [CII]158micron emission line in this z=5.7 dusty star-forming galaxy. We use a pixellated lensing reconstruction code to spatially and kinematically resolve the source-plane [CII] and rest-frame 158 micron dust continuum structure at ~700 pc (~0.12'') resolution. We discuss the [CII] deficit with a pixellated study of the L_[CII]/L_FIR ratio in the source plane. We find that individual pixels within the galaxy follow the same trend found using unresolved observations of other galaxies, indicating that the deficit arises on scales <700 pc. The lensing reconstruction reveals two spatially and kinematically separated components (~1 kpc and ~500 km s^-1 apart) connected by a bridge of gas. Both components are found to be globally unstable, with Toomre Q instability parameters << 1 everywhere. We argue that SPT0346-52 is undergoing a major merger, which is likely driving the intense and compact star formation.

The origin of the high SFR observed in high-z dusty star-forming galaxies is still unknown. Large fractions of dense molecular gas might provide part of the explanation, but there are few observational constraints on the amount of dense gas in high-z systems dominated by star formation. We present the results of our ALMA program targeting dense-gas tracers (HCN(5-4), HCO+(5-4), and HNC(5-4)) in 5 strongly lensed galaxies from the SPT SMG sample. We detected two of these lines (SNR>5) in SPT-125-47 at z=2.51 and tentatively detected all three (SNR~3) in SPT0551-50 at z=3.16. Since a significant fraction of our target lines is not detected, we developed a statistical method to derive unbiased mean properties taking into account both detections and non-detections. On average, the HCN(5-4) and HCO+(5-4) luminosities of our sources are a factor of ~1.7 fainter than expected, based on the local L'HCN(5-4)-LIR relation, but this offset corresponds to only ~2 sigma. We find that both the HCO+/HCN and HNC/HCN flux ratios are compatible with unity. The first ratio is expected for PDRs while the second is consistent with PDRs or XDRs and/or mid-IR pumping of HNC. Our sources are at the high end of the local relation between the star formation efficiency, determined using the LIR/[CI] and LIR/CO ratios, and the dense gas fraction, estimated using the HCN/[CI] and HCN/CO ratios. In SPT0125-47, we found that the velocity profiles of the lines tracing dense (HCN, HCO+) and lower-density (CO, [CI]) gas are similar. In addition to these lines, we obtained one robust and one tentative detection of 13CO(4-3) and found an average I12CO(4-3)/I13CO(4-3) flux ratio of 26.1$_{-3.5}^{+4.5}$, indicating a young but not pristine interstellar medium. We argue that the combination of large and slightly enriched gas reservoirs and high dense-gas fractions could explain the prodigious star formation in these systems.

Galaxies grow inefficiently, with only a few percent of the available gas converted into stars each free-fall time. Feedback processes, such as outflowing winds driven by radiation pressure, supernovae or supermassive black hole accretion, can act to halt star formation if they heat or expel the gas supply. We report a molecular outflow launched from a dust-rich star-forming galaxy at redshift 5.3, one billion years after the Big Bang. The outflow reaches velocities up to 800 km/s relative to the galaxy, is resolved into multiple clumps, and carries mass at a rate within a factor of two of the star formation rate. Our results show that molecular outflows can remove a large fraction of the gas available for star formation from galaxies at high redshift.

Close encounters (CEs) between celestial objects may exert significant influence on their orbits. The influence will be even enhanced when two groups of celestial objects are confined in stable orbital configurations, e.g. in adjacent mean motion resonances (MMRs). Plutinos and Neptune Trojans, trapped in the 2:3 and 1:1 MMRs with Neptune respectively, are such examples. %Meanwhile, many objects among these two groups have peculiar orbits, seemingly as the vestige of CEs. As the first part of our investigation, this paper provides a detailed description of CEs between Plutinos and Trojans and their potential influences on the Trojans' orbits. Statistical analyses of CE data from numerical simulations reveal the randomness lying in the CEs between the two planetesimals. The closest positions of CEs distribute symmetrically inside the given CE region and no particular bias is found between the positive and negative effects on the orbital elements of Trojans. Based on the Gaussian approximation on the distribution of the velocity orientation of Plutino, and the integral derivatives of Gaussian perturbation equations, a theoretical method is built to estimate the CE effects. To further verify the randomness of CEs, a Monte Carlo approach is applied, and it generates distribution features consistent with the numerical results. In summary, CEs brought by realistic Plutinos exert impartial effects and tiny total influence on the orbital elements of Trojans. However, driven by the random walk mechanism, tiny effects may accumulate to a prominent variation given sufficient CEs, which will be discussed in the accompanying paper.

The plasmoid instability in evolving current sheets has been widely studied due to its effects on the disruption of current sheets, the formation of plasmoids, and the resultant fast magnetic reconnection. In this Letter, we study the role of the plasmoid instability in two-dimensional magnetohydrodynamic (MHD) turbulence by means of high-resolution direct numerical simulations. At sufficiently large magnetic Reynolds number ($R_m=10^6$), the combined effects of dynamic alignment and turbulent intermittency lead to a copious formation of plasmoids in a multitude of intense current sheets. The disruption of current sheet structures facilitates the energy cascade towards small scales, leading to the breaking and steepening of the energy spectrum. In the plasmoid-mediated regime, the energy spectrum displays a scaling that is close to the spectral index $-2.2$ as proposed by recent analytic theories. We also demonstrate that the scale-dependent dynamic alignment exists in 2D MHD turbulence and the corresponding slope of the alignment angle is close to 0.25.

The geophysics of extrasolar planets is a scientific topic often regarded as standing largely beyond the reach of near-term observations. This reality in no way diminishes the central role of geophysical phenomena in shaping planetary outcomes, from formation, to thermal and chemical evolution, to numerous issues of surface and near-surface habitability. We emphasize that for a balanced understanding of extrasolar planets, it is important to look beyond the natural biases of current observing tools, and actively seek unique pathways to understand exoplanet interiors as best as possible during the long interim prior to a time when internal components are more directly accessible. Such pathways include but are not limited to: (a) enhanced theoretical and numerical modeling, (b) laboratory research on critical material properties, (c) measurement of geophysical properties by indirect inference from imprints left on atmospheric and orbital properties, and (d) the purpose-driven use of Solar System object exploration expressly for its value in comparative planetology toward exoplanet-analogs. Breaking down barriers that envision local Solar System exploration, including the study of Earth's own deep interior, as separate from and in financial competition with extrasolar planet research, may greatly improve the rate of needed scientific progress for exoplanet geophysics. As the number of known rocky and icy exoplanets grows in the years ahead, we expect demand for expertise in 'exogeoscience' will expand at a commensurately intense pace. We highlight key topics, including: how water oceans below ice shells may dominate the total habitability of our galaxy by volume, how free-floating nomad planets may often attain habitable subsurface oceans supported by radionuclide decay, and how deep interiors may critically interact with atmospheric mass loss via dynamo-driven magnetic fields.

We study roles of the thermosphere and exosphere on the Martian ionospheric structure and ion escape rates in the process of the solar wind-Mars interaction. We employ a four-species multifluid MHD (MF-MHD) model to simulate the Martian ionosphere and magnetosphere. The $cold$ thermosphere background is taken from the Mars Global Ionosphere Thermosphere Model (M-GITM) and the $hot$ oxygen exosphere is adopted from the Mars exosphere Monte Carlo model - Adaptive Mesh Particle Simulator (AMPS). A total of four cases with the combination of 1D (globally averaged) and 3D thermospheres and exospheres are studied. The ion escape rates calculated by adopting 1D and 3D atmospheres are similar; however, the latter are required to adequately reproduce MAVEN ionospheric observations. In addition, our simulations show that the 3D hot oxygen corona plays an important role in preventing planetary molecular ions (O$_2^+$ and CO$_2^+$) escaping from Mars, mainly resulting from the mass loading of the high-altitude exospheric O$^+$ ions. The $cold$ thermospheric oxygen atom, however, is demonstrated to be the primary neutral source for O$^+$ ion escape during the relatively weak solar cycle 24.

Jupiter's radio emission has been linked to its planetary-scale magnetic field, and spacecraft investigations have revealed that most planets, and some moons, have or had a global magnetic field. Generated by internal dynamos, magnetic fields are one of the few remote sensing means of constraining the properties of planetary interiors. For the Earth, its magnetic field has been speculated to be partially responsible for its habitability, and knowledge of an extrasolar planet's magnetic field may be necessary to assess its habitability. The radio emission from Jupiter and other solar system planets is produced by an electron cyclotron maser, and detections of extrasolar planetary electron cyclotron masers will enable measurements of extrasolar planetary magnetic fields. This white paper draws heavily on the W. M. Keck Institute for Space Studies report Planetary Magnetic Fields: Planetary Interiors and Habitability (Lazio, Shkolnik, Hallinan, et al.), it incorporates topics discussed at the American Astronomical Society Topical Conference "Radio Exploration of Planetary Habitability," it complements the Astrobiology Science Strategy white paper "Life Beyond the Solar System: Space Weather and Its Impact on Habitable Worlds" (Airapetian et al.), and it addresses aspects of planetary magnetic fields discussed in the NASA Astrobiology Strategy.