Search

Flows driven by photons have been studied for almost a century, and a quantitative description of the radiative forces on atoms and ions is important for understanding a wide variety of systems with outflows and accretion disks, such as active galactic nuclei. Quantifying the associated forces is crucial to determining how these outflows enable interactive mechanisms within these environments, such as AGN feedback. The total number of spectral lines in any given ion of the outflow material must be tabulated in order to give a complete characterization of this force. Here we provide calculations of the dimensionless line force multiplier for AGN environments. For a wide array of representative AGN sources, we explicitly calculate the photoionization balance at the proposed wind-launching region above the accretion disk, compute the strength of the line-driving force on the gas, and revisit and formalize the role of the commonly-used ionization parameter $\xi$ in ultimately determining the line-driving force. We perform these computations and analyses for a variety of AGN central source properties, such as black hole mass, initial wind velocity, and number density. We find that, while useful, the ionization parameter provides an incomplete description of the overall ionization state of the outflow material. We use these findings to provide an updated method for calculating the strength of the radiative line-driving using both the X-ray spectral index $\Gamma_X$ and the ionization parameter.

Magnetic bright points on the solar photosphere mark the footpoints of kilogauss magnetic flux tubes extending toward the corona. Convective buffeting of these tubes is believed to excite magnetohydrodynamic waves, which can propagate to the corona and there deposit heat. Measuring wave excitation via bright-point motion can thus constrain coronal and heliospheric models, and this has been done extensively with centroid tracking, which can estimate kink-mode wave excitation. DKIST is the first telescope to provide well-resolved observations of bright points, allowing shape and size measurements to probe the excitation of other wave modes that have been difficult, if not impossible, to study to date. In this work, we demonstrate a method of automatic bright-point tracking that robustly identifies the shapes of bright points, and we develop a technique for interpreting measured bright-point shape changes as the driving of a range of thin-tube wave modes. We demonstrate these techniques on a MURaM simulation of DKIST-like resolution. These initial results suggest that modes other than the long-studied kink mode could increase the total available energy budget for wave-heating by 50%. Pending observational verification as well as modeling of the propagation and dissipation of these additional wave modes, this could represent a significant increase in the potency of wave-turbulence heating models., Comment: Accepted in the Astrophysical Journal. 24 pages, 11 figures. arXiv admin note: substantial text overlap with arXiv:2108.10987

The solar wind is the extension of the Sun's hot and ionized corona, and it exists in a state of continuous expansion into interplanetary space. The radial distance at which the wind's outflow speed exceeds the phase speed of Alfvenic and fast-mode magnetohydrodynamic (MHD) waves is called the Alfven radius. In one-dimensional models, this is a singular point beyond which most fluctuations in the plasma and magnetic field cannot propagate back down to the Sun. In the multi-dimensional solar wind, this point can occur at different distances along an irregularly shaped "Alfven surface." In this article, we review the properties of this surface and discuss its importance in models of solar-wind acceleration, angular-momentum transport, MHD waves and turbulence, and the geometry of magnetically closed coronal loops. We also review the results of simulations and data analysis techniques that aim to determine the location of the Alfven surface. Combined with recent perihelia of Parker Solar Probe, these studies seem to indicate that the Alfven surface spends most of its time at heliocentric distances between about 10 and 20 solar radii. It is becoming apparent that this region of the heliosphere is sufficiently turbulent that there often exist multiple (stochastic and time-dependent) crossings of the Alfven surface along any radial ray. Thus, in many contexts, it is more useful to make use of the concept of a topologically complex "Alfven zone" rather than one closed surface. This article also reviews how the Polarimeter to Unify the Corona and Heliosphere (PUNCH) mission will measure the properties of the Alfven surface and provide key constraints on theories of solar-wind acceleration., Comment: Accepted for publication in Solar Physics, part of the topical collection titled "The Polarimeter to Unify the Corona and Heliosphere (PUNCH) Mission: An Overview." 30 pages, 8 figures

We observe an enhanced stellar wind mass-loss rate from low-mass stars exhibiting higher X-ray flux. This trend, however, does not align with the Sun, where no evident correlation between X-ray flux and mass-loss rate is present. To reconcile these observations, we propose a hybrid model for the stellar wind from solar-type stars, incorporating both Alfv\'en wave dynamics and flux emergence-driven interchange reconnection, an increasingly studied concept guided by the latest heliospheric observations. For establishing a mass-loss rate scaling law, we perform a series of magnetohydrodynamic simulations across varied magnetic activities. Through a parameter survey concerning the surface (unsigned) magnetic flux ($\Phi^{\rm surf}$) and the open-to-surface magnetic flux ratio ($\xi^{\rm open} = \Phi^{\rm open}/\Phi^{\rm surf}$), we derive a scaling law of the mass-loss rate given by $\dot{M}_w/\dot{M}_{w,\odot} = \left( \Phi^{\rm surf} / \Phi^{\rm surf}_\odot \right)^{0.52}\left( \xi^{\rm open} / \xi^{\rm open}_\odot \right)^{0.86}$, where $\dot{M}_{w,\odot} = 2.0 \times 10^{-14} \ M_\odot {\rm \ yr}^{-1}$, $\Phi^{\rm surf}_\odot = 3.0 \times 10^{23} {\rm \ Mx}$, and $\xi^{\rm open}_\odot = 0.2$. By comparing cases with and without flux emergence, we find that the increase in the mass-loss rate with the surface magnetic flux can be attributed to the influence of flux emergence. Our scaling law demonstrates an agreement with solar wind observations spanning 40 years, exhibiting superior performance when compared to X-ray-based estimations. Our findings suggest that flux emergence may play a significant role in the stellar winds of low-mass stars, particularly those originating from magnetically active stars., Comment: Accepted for publication in The Astrophysical Journal

The solar atmosphere is known to contain many different types of wavelike oscillation. Waves and other fluctuations (e.g., turbulent eddies) are believed to be responsible for at least some of the energy transport and dissipation that heats the corona and accelerates the solar wind. Thus, it is important to understand the behavior of magnetohydrodynamic (MHD) waves as they propagate and evolve in different regions of the Sun's atmosphere. In this paper, we investigate how MHD waves can affect the overall plasma state when they reflect and refract at sharp, planar interfaces in density. First, we correct an error in a foundational paper (Stein 1971) that affects the calculation of wave energy-flux conservation. Second, we apply this model to reflection-driven MHD turbulence in the solar wind, where the presence of density fluctuations can enhance the generation of inward-propagating Alfven waves. This model reproduces the time-averaged Elsasser imbalance fraction (i.e., ratio of inward to outward Alfvenic power) from several published numerical simulations. Lastly, we model how the complex magnetic field threading the transition region between the chromosphere and corona helps convert a fraction of upward-propagating Alfven waves into fast-mode and slow-mode MHD waves. These magnetosonic waves dissipate in a narrow region around the transition region and produce a sharp peak in the heating rate. This newly found source of heating sometimes exceeds the expected heating rate from Alfvenic turbulence by an order of magnitude. It may explain why some earlier models seemed to require an additional ad-hoc heat source at this location., Comment: Accepted for publication in the Astrophysical Journal. 19 pages, 10 figures, 1 table

Acoustic wave heating is believed to contribute significantly to the missing energy input required to maintain the solar chromosphere in its observed state. We studied the propagation of waves above the acoustic cutoff in the upper photosphere into the chromosphere with ultraviolet and optical spectral observations interpreted through comparison with three dimensional radiative magnetohydrodynamic (rMHD) \emph{Bifrost} models to constrain the heating contribution from acoustic waves in the solar atmosphere. Sit-and-stare observations taken with the Interface Region Imaging Spectrograph (IRIS) and data from the Interferometric BIdimensional Spectrograph (IBIS) were used to provide the observational basis of this work. We compared the observations with synthetic observables derived from the Bifrost solar atmospheric model. Our analysis of the \emph{Bifrost} simulations show that internetwork and enhanced network regions exhibit significantly different wave propagation properties, which are important for the accurate wave flux estimates. The inferred wave energy fluxes based on our observations are not sufficient to maintain the solar chromosphere. We point out that the systematics of the modeling approaches in the literature lead to differences which could determine the conclusions of this type of studies, based on the same observations., Comment: Accepted for publication in ApJ

We present millimeter, optical, and soft X-ray observations of a stellar flare with an energy squarely in the regime of typical X1 solar flares. The flare was observed from Proxima Cen on 2019 May 6 as part of a larger multi-wavelength flare monitoring campaign and was captured by Chandra, LCOGT, du Pont, and ALMA. Millimeter emission appears to be a common occurrence in small stellar flares that had gone undetected until recently, making it difficult to interpret these events within the current multi-wavelength picture of the flaring process. The May 6 event is the smallest stellar millimeter flare detected to date. We compare the relationship between the soft X-ray and millimeter emission to that observed in solar flares. The X-ray and optical flare energies of 10$^{30.3\pm0.2}$ and 10$^{28.9\pm0.1}$ erg, respectively, the coronal temperature of T=11.0$\pm$2.1 MK, and the emission measure of 9.5$\pm$2.2 X 10$^{49}$ cm$^{-3}$ are consistent with M-X class solar flares. We find the soft X-ray and millimeter emission during quiescence are consistent with the Gudel-Benz Relation, but not during the flare. The millimeter luminosity is >100X higher than that of an equivalent X1 solar flare and lasts only seconds instead of minutes as seen for solar flares., Comment: 12 pages, 5 figures. Accepted to ApJ

In the solar corona and solar wind, electron heat conduction is an important process that transports energy over large distances and helps determine the spatial variation of temperature. High-density regions undergoing rapid particle-particle collisions exhibit a heat flux described well by classical Spitzer-Harm theory. However, much of the heliosphere is closer to a more collisionless state, and there is no standard description of heat conduction for fluid-based (e.g., magnetohydrodynamic) models that applies generally. Some proposed models rely on electron velocity distributions that exhibit negative values of the phase-space density. In this paper, we explore how positive-definite velocity distributions can be used in fluid-based conservation equations for the electron heat flux along magnetic-field lines in the corona and solar wind. We study both analytic forms of skewed distributions (e.g., skew-normal distributions, two-sided bi-Maxwellians, and constant-collision-time electrostatic solutions) and empirical fits to measurements of core, halo, and strahl electrons in interplanetary space. We also present example solutions to a generalized conservation equation for the heat flux in the solar wind, with some limiting cases found to resemble known free-streaming approximations. The resulting values of the electron heat flux vary as a function of radial distance and Knudsen number in ways that resemble observed data. We note that this model does not include the effects of kinetic instabilities (which may impose saturation limits when active), so for now its regime of applicability is limited to collisionless heat-flux evolution away from the known instability boundaries in parameter space., Comment: Accepted for publication in the Journal of Geophysical Research: Space Physics. 30 pages, 5 figures. Associated data are in a third-party repository at https://doi.org/10.5281/zenodo.4927859

We present observational constraints on the solar chromospheric heating contribution from acoustic waves with frequencies between 5 and 50 mHz. We utilize observations from the Dunn Solar Telescope in New Mexico complemented with observations from the Atacama Large Millimeter Array collected on 2017 April 23. The properties of the power spectra of the various quantities are derived from the spectral lines of Ca II 854.2 nm, H I 656.3 nm, and the millimeter continuum at 1.25 mm and 3 mm. At the observed frequencies the diagnostics almost all show a power law behavior, whose particulars (slope, peak and white noise floors) are correlated with the type of solar feature (internetwork, network, plage). In order to disentangle the vertical versus transverse plasma motions we examine two different fields of view; one near disk center and the other close to the limb. To infer the acoustic flux in the middle chromosphere, we compare our observations with synthetic observables from the time-dependent radiative hydrodynamic RADYN code. Our findings show that acoustic waves carry up to about 1 kW m$^{-2}$ of energy flux in the middle chromosphere, which is not enough to maintain the quiet chromosphere, contrary to previous publications., Comment: Accepted for publication in ApJ

We present the discovery of an extreme flaring event from Proxima Cen by ASKAP, ALMA, HST, TESS, and the du Pont Telescope that occurred on 2019 May 1. In the millimeter and FUV, this flare is the brightest ever detected, brightening by a factor of >1000 and >14000 as seen by ALMA and HST, respectively. The millimeter and FUV continuum emission trace each other closely during the flare, suggesting that millimeter emission could serve as a proxy for FUV emission from stellar flares and become a powerful new tool to constrain the high-energy radiation environment of exoplanets. Surprisingly, optical emission associated with the event peaks at a much lower level with a time delay. The initial burst has an extremely short duration, lasting for <10 sec. Taken together with the growing sample of millimeter M dwarf flares, this event suggests that millimeter emission is actually common during stellar flares and often originates from short burst-like events., Comment: 11 pages, 4 figures, 1 appendix, published in ApJ Letters

Light curves produced by the Kepler mission demonstrate stochastic brightness fluctuations (or "flicker") of stellar origin which contribute to the noise floor, limiting the sensitivity of exoplanet detection and characterization methods. In stars with surface convection, the primary driver of these variations on short (sub-eight-hour) timescales is believed to be convective granulation. In this work, we improve existing models of this granular flicker amplitude, or $F_8$, by including the effect of the Kepler bandpass on measured flicker, by incorporating metallicity in determining convective Mach numbers, and by using scaling relations from a wider set of numerical simulations. To motivate and validate these changes, we use a recent database of convective flicker measurements in Kepler stars, which allows us to more fully detail the remaining model--prediction error. Our model improvements reduce the typical misprediction of flicker amplitude from a factor of 2.5 to 2. We rule out rotation period and strong magnetic activity as possible explanations for the remaining model error, and we show that binary companions may affect convective flicker. We also introduce an "envelope" model which predicts a range of flicker amplitudes for any one star to account for some of the spread in numerical simulations, and we find that this range covers 78% of observed stars. We note that the solar granular flicker amplitude is lower than most Sun-like stars. This improved model of convective flicker amplitude can better characterize this source of noise in exoplanet studies as well as better inform models and simulations of stellar granulation., Comment: Accepted for publication in the Astrophysical Journal. 23 pages, 12 figures. Code and data available at https://doi.org/10.5281/zenodo.4444282

Although the solar wind flows primarily outward from the Sun to interplanetary space, there are times when small-scale plasma inflows are observed. Inward-propagating density fluctuations in polar coronal holes were detected by the COR2 coronagraph on board the STEREO-A spacecraft at heliocentric distances of 7 to 12 solar radii, and these fluctuations appear to undergo substantial deceleration as they move closer to the Sun. Models of linear magnetohydrodynamic waves have not been able to explain these deceleration patterns, so they have been interpreted more recently as jets from coronal sites of magnetic reconnection. In this paper, we develop a range of dynamical models of discrete plasma parcels with the goal of better understanding the observed deceleration trend. We found that parcels with a constant mass do not behave like the observed flows, and neither do parcels undergoing ablative mass loss. However, parcels that accrete mass in a snowplow-like fashion can become decelerated as observed. We also extrapolated OMNI in situ data down to the so-called Alfven surface and found that the initial launch-point for the observed parcels may often be above this critical radius. In other words, in order for the parcels to flow back down to the Sun, their initial speeds are probably somewhat nonlinear (i.e., supra-Alfvenic) and thus the parcels may be associated with structures such as shocks, jets, or shear instabilities., Comment: Accepted for publication in the Astrophysical Journal. 16 pages, 8 figures

Radiation contributes to the acceleration of large-scale flows in various astrophysical environments because of the strong opacity in spectral lines. Quantification of the associated force is crucial to understanding these line-driven flows, and a large number of lines (due to the full set of elements and ionization stages) must be taken into account. Here we provide new calculations of the dimensionless line strengths and associated opacity-dependent force multipliers for an updated list of approximately 4.5 million spectral lines compiled from the NIST, CHIANTI, CMFGEN, and TOPbase databases. To maintain generality of application to different environments, we assume local thermodynamic equilibrium, illumination by a Planck function, and the Sobolev approximation. We compute the line forces in a two-dimensional grid of temperatures (i.e., values between 5,200 and 70,000 K) and densities (varying over 11 orders of magnitude). Historically, the force multiplier function has been described by a power-law function of optical depth. We revisit this assumption by fitting alternative functions that include a saturation to a constant value (Gayley's $\bar{Q}$ parameter) in the optically-thin limit. This alternate form is a better fit than the power-law form, and we use it to calculate example mass-loss rates for massive main-sequence stars. Because the power-law force multiplier does not continue to arbitrarily small optical depths, we find a sharp decrease, or quenching, of line-driven winds for stars with effective temperatures less than about 15,000 K., Comment: 20 pages, 14 figures, 1 table

One of the most important early results from the Parker Solar Probe (PSP) is the ubiquitous presence of magnetic switchbacks, whose origin is under debate. Using a three-dimensional direct numerical simulation of the equations of compressible magnetohydrodynamics from the corona to 40 solar radii, we investigate whether magnetic switchbacks emerge from granulation-driven Alfv\'en waves and turbulence in the solar wind. The simulated solar wind is an Alfv\'enic slow-solar-wind stream with a radial profile consistent with various observations, including observations from PSP. As a natural consequence of Alfv\'en-wave turbulence, the simulation reproduced magnetic switchbacks with many of the same properties as observed switchbacks, including Alfv\'enic v-b correlation, spherical polarization (low magnetic compressibility), and a volume filling fraction that increases with radial distance. The analysis of propagation speed and scale length shows that the magnetic switchbacks are large-amplitude (nonlinear) Alfv\'en waves with discontinuities in the magnetic field direction. We directly compare our simulation with observations using a virtual flyby of PSP in our simulation domain. We conclude that at least some of the switchbacks observed by PSP are a natural consequence of the growth in amplitude of spherically polarized Alfv\'en waves as they propagate away from the Sun., Comment: accepted for publication in The Astrophysical Journal

The high-speed solar wind is typically the simplest and least stochastic type of large-scale plasma flow in the heliosphere. For much of the solar cycle, it is connected magnetically to large polar coronal holes on the Sun's surface. Because these features are relatively well-known (and less complex than the multiple source-regions of the slow wind), the fast wind is often a useful testing-ground for theoretical models of coronal heating. In order to provide global empirical constraints to these models, here we collect together some older and more recent measurements of the temperatures of protons, electrons, and oxygen ions as a function of radial distance., Comment: Submitted to AAS journals. 3 pages, 1 figure

Strongly magnetic B-type stars with moderately rapid rotation form `centrifugal magnetospheres' (CMs), from the magnetic trapping of stellar wind material in a region above the Kepler co-rotation radius. A longstanding question is whether the eventual loss of such trapped material occurs from gradual drift and/or diffusive leakage, or through sporadic `{\em centrifugal break out}' (CBO) events, wherein magnetic tension can no longer contain the built-up mass. We argue here that recent empirical results for Balmer-$\alpha$ emission from such B-star CMs strongly favor the CBO mechanism. Most notably, the fact that the onset of such emission depends mainly on the field strength at the Kepler radius, and is largely {\em independent} of the stellar luminosity, strongly disfavors any drift/diffusion process, for which the net mass balance would depend on the luminosity-dependent wind feeding rate. In contrast, we show that in a CBO model the {\em maximum confined mass} in the magnetosphere is independent of this wind feeding rate, and has a dependence on field strength and Kepler radius that naturally explains the empirical scalings for the onset of H$\alpha$ emission, its associated equivalent width, and even its line profile shapes. However, the general lack of observed Balmer emission in late-B and A-type stars could still be attributed to a residual level of diffusive or drift leakage that does not allow their much weaker winds to fill their CMs to the breakout level needed for such emission; alternatively this might result from a transition to a metal-ion wind that lacks the requisite Hydrogen., Comment: 13 pages, 11 figures, accepted for MNRAS

Extreme-ultraviolet and X-ray emissions from stellar coronae drive mass loss from exoplanet atmospheres, and ultraviolet emission from stellar chromospheres drives photo-chemistry in exoplanet atmospheres. Comparisons of the spectral energy distributions of host stars are, therefore, essential for understanding the evolution and habitability of exoplanets. The large number of stars observed with the MUSCLES, Mega-MUSCLES, and other recent HST observing programs has provided for the first time a large sample (79 stars) of reconstructed Lyman-alpha fluxes that we compare with X-ray fluxes to identify significant patterns in the relative emission from these two atmospheric regions as a function of stellar age and effective temperature. We find that as stars age on the main sequence, the emissions from their chromospheres and coronae follow a pattern in response to the amount of magnetic heating in these atmospheric layers. A single trendline slope describes the pattern of X-ray vs. Lyman-alpha emission for G and K dwarfs, but the different trendlines for M dwarf stars show that the Lyman-alpha fluxes of M stars are significantly smaller than warmer stars with the same X-ray flux. The X-ray and Lyman-alpha luminosities divided by the stellar bolometric luminosities show different patterns depending on stellar age. The L(Lyman-alpha)/L(bol) ratios increase smoothly to cooler stars of all ages, but the L(X)/L(bol) ratios show different trends. For older stars, the increase in coronal emission with decreasing T(eff) is much steeper than chromospheric emission. We suggest a fundamental link between atmospheric properties and trendlines relating coronal and chromospheric heating, Comment: 42 pages, 9 figures. To appear in the Astrophysical Journal

When observing spectral lines in the optically-thin corona, line-of-sight (LOS) effects can strongly affect the interpretation of the data, especially in regions just above the limb. We present a semi-empirical forward model, called GHOSTS, to characterize these effects. GHOSTS uses inputs from several other models to compute non-equilibrium ionization states (which include the solar-wind freezing-in effect) for many ions. These are used to generate ensembles of simulated spectral lines that are examined in detail, with emphasis on: (1) relationships between quantities derived from observables and the radial variation of the observed quantities, (2) the behavior of thermal and non-thermal components of the line width, and (3) relative contributions of collisionally excited and radiatively scattered photons. We find that rapidly changing temperatures in the low corona can cause ion populations to vary dramatically with height. This can lead to line-width measurements that are constant with height (a "plateau" effect) even when the temperature is increasing rapidly, as the plane-of-sky becomes evacuated and the foreground/background plasma dominates the observation. We find that LOS effects often drive the velocity width to be close to the plane-of-sky value of the wind speed, despite it flowing perpendicularly to the LOS there. The plateau effect can also cause the non-thermal component of the line width to greatly exceed the solar wind velocity at the observation height. Lastly, we study how much of the LOS is significant to the observation, and the importance of including continuum in the solar spectrum when computing the radiatively scattered emission., Comment: 28 pages, 17 figures. Accepted for publication by the ApJ on 8/20/2020

The Daniel K. Inouye Solar Telescope (DKIST) will revolutionize our ability to measure, understand and model the basic physical processes that control the structure and dynamics of the Sun and its atmosphere. The first-light DKIST images, released publicly on 29 January 2020, only hint at the extraordinary capabilities which will accompany full commissioning of the five facility instruments. With this Critical Science Plan (CSP) we attempt to anticipate some of what those capabilities will enable, providing a snapshot of some of the scientific pursuits that the Daniel K. Inouye Solar Telescope hopes to engage as start-of-operations nears. The work builds on the combined contributions of the DKIST Science Working Group (SWG) and CSP Community members, who generously shared their experiences, plans, knowledge and dreams. Discussion is primarily focused on those issues to which DKIST will uniquely contribute.

Ultraviolet spectroscopy of the extended solar corona is a powerful tool for measuring the properties of protons, electrons, and heavy ions in the accelerating solar wind. The large coronal holes that expand up from the north and south poles at solar minimum are low-density collisionless regions in which it is possible to detect departures from one-fluid thermal equilibrium. An accurate characterization of these departures is helpful in identifying the kinetic processes ultimately responsible for coronal heating. In this paper, Ultraviolet Coronagraph Spectrometer (UVCS) measurements of the H I Lyman alpha line are analyzed to constrain values for the solar wind speed, electron density, electron temperature, proton temperature (parallel and perpendicular to the magnetic field) and Alfven-wave amplitude. The analysis procedure involves creating a large randomized ensemble of empirical models, simulating their Lyman alpha profiles, and building posterior probability distributions for only the models that agree with the UVCS data. The resulting temperatures do not exhibit a great deal of radial variation between heliocentric distances of 1.4 and 4 solar radii. Typical values for the electron, parallel proton, and perpendicular proton temperatures are 1.2, 1.8, and 1.9 MK, respectively. Resulting values for the "nonthermal" Alfven wave amplitude show evidence for weak dissipation, with a total energy-loss rate that agrees well with an independently derived total heating rate for the protons and electrons. The moderate Alfven-wave amplitudes appear to resolve some tension in the literature between competing claims of both higher (undamped) and lower (heavily damped) values., Comment: Accepted for publication in the Astrophysical Journal. 23 pages, 13 figures. Associated data are in a third-party repository at https://doi.org/10.5281/zenodo.3908519

Observations of stellar rotation show that low-mass stars lose angular momentum during the main sequence. We simulate the winds of Sun-like stars with a range of rotation rates, covering the fast and slow magneto-rotator regimes, including the transition between the two. We generalize an Alfv\'en-wave driven solar wind model that builds on previous works by including the magneto-centrifugal force explicitly. In this model, the surface-averaged open magnetic flux is assumed to scale as $B_\ast f^{\rm open}_\ast \propto {\rm Ro}^{-1.2}$, where $f^{\rm open}_\ast$ and ${\rm Ro}$ are the surface open-flux filling factor and Rossby number, respectively. We find that, 1. the angular momentum loss rate (torque) of the wind is described as $\tau_w \approx 2.59 \times 10^{30} {\rm \ erg} \ \left( \Omega_\ast / \Omega_\odot \right)^{2.82}$, yielding a spin-down law $\Omega_\ast \propto t^{-0.55}$. 2. the mass-loss rate saturates at $\dot{M}_w \sim 3.4 \times 10^{-14} M_\odot {\rm \ yr^{-1}}$, due to the strong reflection and dissipation of Alfv\'en waves in the chromosphere. This indicates that the chromosphere has a strong impact in connecting the stellar surface and stellar wind. Meanwhile, the wind ram pressure scales as $P_w \propto \Omega_\ast^{0.57}$, which is able to explain the lower-envelope of the observed stellar winds by Wood et al. 3. the location of the Alfv\'en radius is shown to scale in a way that is consistent with 1D analytic theory. Additionally, the precise scaling of the Alfv\'en radius matches previous works which used thermally-driven winds. Our results suggest that the Alfv\'en-wave driven magnetic rotator wind plays a dominant role in the stellar spin-down during the main-sequence., Comment: accepted for publication in The Astrophysical Journal

We summarize some of the compelling new scientific opportunities for understanding stars and stellar systems that can be enabled by sub-milliarcsec (sub-mas) angular resolution, UV-Optical spectral imaging observations, which can reveal the details of the many dynamic processes (e.g., evolving magnetic fields, accretion, convection, shocks, pulsations, winds, and jets) that affect stellar formation, structure, and evolution. These observations can only be provided by long-baseline interferometers or sparse aperture telescopes in space, since the aperture diameters required are in excess of 500 m (a regime in which monolithic or segmented designs are not and will not be feasible) and since they require observations at wavelengths (UV) not accessible from the ground. Such observational capabilities would enable tremendous gains in our understanding of the individual stars and stellar systems that are the building blocks of our Universe and which serve as the hosts for life throughout the Cosmos., Comment: Astro2020 Decadal Survey White Paper. arXiv admin note: substantial text overlap with arXiv:0903.2433

We present the first high-resolution, simultaneous observations of the solar chromosphere in the optical and millimeter wavelength ranges, obtained with ALMA and the IBIS instrument at the Dunn Solar Telescope. In this paper we concentrate on the comparison between the brightness temperature observed in ALMA Band 3 (3 mm; 100 GHz) and the core width of the H$\alpha$ 656.3 nm line, previously identified as a possible diagnostic of the chromospheric temperature. We find that in the area of plage, network and fibrils covered by our FOV the two diagnostics are well correlated, with similar spatial structures observed in both. The strength of the correlation is remarkable, given that the source function of the mm-radiation obeys local thermodynamic equilibrium, while the H$\alpha$ line has a source function that deviates significantly from the local Planck function. The observed range of ALMA brightness temperatures is sensibly smaller than the temperature range that was previously invoked to explain the observed width variations in H$\alpha$. We employ analysis from forward modeling with the RH code to argue that the strong correlation between H$\alpha$ width and ALMA brightness temperature is caused by their shared dependence on the population number $n_2$ of the first excited level of hydrogen. This population number drives millimeter opacity through hydrogen ionization via the Balmer continuum, and H$\alpha$ width through a curve-of-growth-like opacity effect. Ultimately, the $n_2$ population is regulated by the enhancement or lack of downward Ly$\alpha$ flux, which coherently shifts the formation height of both diagnostics to regions with different temperature, respectively., Comment: Accepted for publication in ApJ

We examine the different element abundances exhibited by the closed loop solar corona and the slow speed solar wind. Both are subject to the First Ionization Potential (FIP) Effect, the enhancement in coronal abundance of elements with FIP below 10 eV (e.g. Mg, Si, Fe) with respect to high FIP elements (e.g. O, Ne, Ar), but with subtle differences. Intermediate elements, S, P, and C, with FIP just above 10 eV, behave as high FIP elements in closed loops, but are fractionated more like low FIP elements in the solar wind. On the basis of FIP fractionation by the ponderomotive force in the chromosphere, we discuss fractionation scenarios where this difference might originate. Fractionation low in the chromosphere where hydrogen is neutral enhances the S, P and C abundances. This arises with nonresonant waves, which are ubiquitous in open field regions, and is also stronger with torsional Alfven waves, as opposed to shear (i.e. planar) waves. We discuss the bearing these findings have on models of interchange reconnection as the source of the slow speed solar wind. The outflowing solar wind must ultimately be a mixture of the plasma in the originally open and closed fields, and the proportions and degree of mixing should depend on details of the reconnection process. We also describe novel diagnostics in ultraviolet and extreme ultraviolet spectroscopy now available with these new insights, with the prospect of investigating slow speed solar wind origins and the contribution of interchange reconnection by remote sensing., Comment: 21 pages, accepted by ApJ

The corona is a layer of hot plasma that surrounds the Sun, traces out its complex magnetic field, and ultimately expands into interplanetary space as the supersonic solar wind. Although much has been learned in recent decades from advances in observations, theory, and computer simulations, we still have not identified definitively the physical processes that heat the corona and accelerate the solar wind. In this review, we summarize these recent advances and speculate about what else is required to finally understand the fundamental physics of this complex system. Specifically: (1) We discuss recent sub-arcsecond observations of the corona, some of which appear to provide evidence for tangled and braided magnetic fields, and some of which do not. (2) We review results from three-dimensional numerical simulations that, despite limitations in dynamic range, reliably contain sufficient heating to produce and maintain the corona. (3) We provide a new tabulation of scaling relations for a number of proposed coronal heating theories that involve waves, turbulence, braiding, nanoflares, and helicity conservation. An understanding of these processes is important not only for improving our ability to forecast hazardous space-weather events, but also for establishing a baseline of knowledge about a well-resolved star that is relevant to other astrophysical systems., Comment: To be published in the Annual Review of Astronomy and Astrophysics, volume 57 (2019). 30 pages, 6 figures. This is the authors' post-review preprint version; the final version will be available at https://www.annualreviews.org/journal/astro

From the the solar photosphere to the outer heliosphere, the Sun's plasma properties are fluctuating with a broad range of temporal and spatial scales. In fact, a turbulent cascade of energy from large to small scales is a frequently invoked explanation for heating the corona and accelerating the solar wind. NASA's Parker Solar Probe (PSP) is expected to revolutionize our understanding of coronal heating and magnetohydrodynamic (MHD) turbulence by performing in situ sampling closer to the Sun than any other prior space mission. This research note presents theoretical predictions for some properties of MHD turbulence (e.g., spacecraft-frame power spectra and variance anisotropies) in the regions to be explored by PSP. These results are derived from a previously published semi-empirical model of coupled Alfvenic and fast-mode turbulence in the fast solar wind. The primary reason for this research note is to show how straightforward it can be to extract useful predictions from existing theoretical models about measurable quantities that were not even considered when creating the models initially. The variance anisotropy, for example, may be an important quantity for distinguishing between different theoretical models for coronal heating and solar wind acceleration., Comment: Research Notes of the AAS, accepted for publication. "Data behind the figure" are included in ASCII machine-readable table form, both here and on the RNAAS site

The solar corona is a hot, dynamic, and highly magnetized plasma environment whose source of energy is not yet well understood. One leading contender for that energy source is the dissipation of magnetohydrodynamic (MHD) waves or turbulent fluctuations. Many wave-heating models for the corona and the solar wind presume that these fluctuations originate at or below the Sun's photosphere. However, this paper investigates the idea that magnetic reconnection may generate an additional source of MHD waves over a gradual range of heights in the low corona. A time-dependent Monte Carlo simulation of the mixed-polarity magnetic field is used to predict the properties of reconnection-driven coronal MHD waves. The total power in these waves is typically small in comparison to that of photosphere-driven waves, but their frequencies are much lower. Reconnection-driven waves begin to dominate the total power spectrum at periods longer than about 30 minutes. Thus, they may need to be taken into account in order to understand the low-frequency power-law spectra observed by both coronal spectropolarimetry and in-situ particle/field instruments. These low-frequency Alfven waves should carry more magnetic energy than kinetic energy, and thus they may produce less nonthermal Doppler broadening (in comparison to photosphere-driven high-frequency waves) in emission lines observed above the solar limb., Comment: Accepted for publication in the Astrophysical Journal. 13 pages, 6 figures

We present new analyses of ALMA 12-m and ACA observations at 233 GHz (1.3 mm) of the Proxima Centauri system with sensitivities of 9.5 and 47 $\mu$Jy beam$^{-1}$, respectively, taken from 2017 January 21 through 2017 April 25. These analyses reveal that the star underwent a significant flaring event during one of the ACA observations on 2017 March 24. The complete event lasted for approximately 1 minute and reached a peak flux density of $100\pm4$ mJy, nearly a factor of $1000\times$ brighter than the star's quiescent emission. At the flare peak, the continuum emission is characterized by a steeply falling spectral index with frequency, $F_\nu \propto \nu^\alpha$ with $\alpha = -1.77\pm0.45$, and a lower limit on the fractional linear polarization of $|Q/I| = 0.19\pm0.02$. Since the ACA observations do not show any quiescent excess emission, we conclude that there is no need to invoke the presence of a dust belt at $1-4$ AU. We also posit that the slight excess flux density of $101\pm9$ $\mu$Jy observed in the 12-m observations compared to the photospheric flux density of $74\pm4$ $\mu$Jy extrapolated from infrared wavelengths may be due to coronal heating from continual smaller flares, as is seen for AU Mic, another nearby, well-studied, M dwarf flare star. If this is true, then the need for warm dust at $\sim0.4$ AU is also removed., Comment: 7 pages, 3 figures; Submitted to ApJL 2017 December 22; Accepted to ApJL 2018 February 5

In the subset of luminous, early-type stars with strong, large-scale magnetic fields and moderate to rapid rotation, material from the star's radiatively driven stellar wind outflow becomes trapped by closed magnetic loops, forming a centrifugally supported, co-rotating magnetosphere. We present here a semi-analytic analysis of how this quasi-steady accumulation of wind mass can be balanced by losses associated with a combination of an outward, centrifugally driven drift in the region beyond the Kepler co-rotation radius, and an inward/outward diffusion near this radius. We thereby derive scaling relations for the equilibrium spatial distribution of mass, and the associated emission measure for observational diagnostics like Balmer line emission. We discuss the potential application of these relations for interpreting surveys of the emission line diagnostics for OB stars with centrifugally supported magnetospheres. For a specific model of turbulent field-line-wandering rooted in surface motions associated with the iron opacity bump, we estimate values for the associated diffusion and drift coefficients., Comment: 11 pages, 5 figures, accepted by MNRAS

Magnetic bright points in the solar photosphere, visible in both continuum and G-band images, indicate footpoints of kilogauss magnetic flux tubes extending to the corona. The power spectrum of bright-point motion is thus also the power spectrum of Alfven wave excitation, transporting energy up flux tubes into the corona. This spectrum is a key input in coronal and heliospheric models. We produce a power spectrum of bright-point motion using radiative magnetohydrodynamic simulations, exploiting spatial resolution higher than can be obtained in present-day observations, while using automated tracking to produce large data quantities. We find slightly higher amounts of power at all frequencies compared to observation-based spectra, while confirming the spectrum shape of recent observations. This also provides a prediction for observations of bright points with DKIST, which will achieve similar resolution and high sensitivity. We also find a granule size distribution in support of an observed two-population distribution, and we present results from tracking passive tracers which show a similar power spectrum to that of bright points. Finally, we introduce a simplified, laminar model of granulation, with which we explore the roles of turbulence and of the properties of the granulation pattern in determining bright-point motion., Comment: 13 pages, 8 figures. Accepted for publication in ApJ

The Sun's outer atmosphere is heated to temperatures of millions of degrees, and solar plasma flows out into interplanetary space at supersonic speeds. This paper reviews our current understanding of these interrelated problems: coronal heating and the acceleration of the ambient solar wind. We also discuss where the community stands in its ability to forecast how variations in the solar wind (i.e., fast and slow wind streams) impact the Earth. Although the last few decades have seen significant progress in observations and modeling, we still do not have a complete understanding of the relevant physical processes, nor do we have a quantitatively precise census of which coronal structures contribute to specific types of solar wind. Fast streams are known to be connected to the central regions of large coronal holes. Slow streams, however, appear to come from a wide range of sources, including streamers, pseudostreamers, coronal loops, active regions, and coronal hole boundaries. Complicating our understanding even more is the fact that processes such as turbulence, stream-stream interactions, and Coulomb collisions can make it difficult to unambiguously map a parcel measured at 1 AU back down to its coronal source. We also review recent progress -- in theoretical modeling, observational data analysis, and forecasting techniques that sit at the interface between data and theory -- that gives us hope that the above problems are indeed solvable., Comment: Accepted for publication in Space Science Reviews. Special issue connected with a 2016 ISSI workshop on "The Scientific Foundations of Space Weather." 44 pages, 9 figures

Coronal mass ejections (CMEs) are eruptive events that cause a solar-type star to shed mass and magnetic flux. CMEs tend to occur together with flares, radio storms, and bursts of energetic particles. On the Sun, CME-related mass loss is roughly an order of magnitude less intense than that of the background solar wind. However, on other types of stars, CMEs have been proposed to carry away much more mass and energy than the time-steady wind. Earlier papers have used observed correlations between solar CMEs and flare energies, in combination with stellar flare observations, to estimate stellar CME rates. This paper sidesteps flares and attempts to calibrate a more fundamental correlation between surface-averaged magnetic fluxes and CME properties. For the Sun, there exists a power-law relationship between the magnetic filling factor and the CME kinetic energy flux, and it is generalized for use on other stars. An example prediction of the time evolution of wind/CME mass-loss rates for a solar-mass star is given. A key result is that for ages younger than about 1 Gyr (i.e., activity levels only slightly higher than the present-day Sun), the CME mass loss exceeds that of the time-steady wind. At younger ages, CMEs carry 10 to 100 times more mass than the wind, and such high rates may be powerful enough to dispel circumstellar disks and affect the habitability of nearby planets. The cumulative CME mass lost by the young Sun may have been as much as 1% of a solar mass., Comment: Accepted for publication in the Astrophysical Journal. 10 pages, 5 figures

Coronal loops trace out bipolar, arch-like magnetic fields above the Sun's surface. Recent measurements that combine rotational tomography, extreme ultraviolet imaging, and potential-field extrapolation have shown the existence of large loops with inverted temperature profiles; i.e., loops for which the apex temperature is a local minimum, not a maximum. These "down loops" appear to exist primarily in equatorial quiet regions near solar minimum. We simulate both these and the more prevalent large-scale "up loops" by modeling coronal heating as a time-steady superposition of: (1) dissipation of incompressible Alfven-wave turbulence, and (2) dissipation of compressive waves formed by mode conversion from the initial population of Alfven waves. We found that when a large percentage (> 99%) of the Alfven waves undergo this conversion, heating is greatly concentrated at the footpoints and stable "down loops" are created. In some cases we found loops with three maxima that are also gravitationally stable. Models that agree with the tomographic temperature data exhibit higher gas pressures for "down loops" than for "up loops," which is consistent with observations. These models also show a narrow range of Alfven wave amplitudes: 3 to 6 km/s at the coronal base. This is low in comparison to typical observed amplitudes of 20 to 30 km/s in bright X-ray loops. However, the large-scale loops we model are believed to comprise a weaker diffuse background that fills much of the volume of the corona. By constraining the physics of loops that underlie quiescent streamers, we hope to better understand the formation of the slow solar wind., Comment: Accepted for publication in the Astrophysical Journal. 15 pages (emulateapj style), 11 figures

The Solar Probe Plus (SPP) mission will explore the Sun's corona and innermost solar wind starting in 2018. The spacecraft will also come close to a number of Mercury-crossing asteroids with perihelia less than 0.3 AU. At small heliocentric distances, these objects may begin to lose mass, thus becoming "active asteroids" with comet-like comae or tails. This paper assembles a database of 97 known Mercury-crossing asteroids that may be encountered by SPP, and it presents estimates of their time-dependent visible-light fluxes and mass loss rates. Assuming a similar efficiency of sky background subtraction as was achieved by STEREO, we find that approximately 80% of these asteroids are bright enough to be observed by the Wide-field Imager for SPP (WISPR). A model of gas/dust mass loss from these asteroids is developed and calibrated against existing observations. This model is used to estimate the visible-light fluxes and spatial extents of spherical comae. Observable dust clouds occur only when the asteroids approach the Sun closer than 0.2 AU. The model predicts that during the primary SPP mission between 2018 and 2025, there should be 113 discrete events (for 24 unique asteroids) during which the modeled comae have angular sizes resolvable by WISPR. The largest of these correspond to asteroids 3200 Phaethon, 137924, 155140, and 289227, all with angular sizes of roughly 15 to 30 arcminutes. We note that the SPP trajectory may still change, but no matter the details there should still be multiple opportunities for fruitful asteroid observations., Comment: Accepted for publication in Earth, Moon, and Planets. 28 pages, 12 figures

Recent IRIS observations have revealed a prevalence of intermittent small-scale jets with apparent speeds of 80 - 250 km s$^{-1}$, emanating from small-scale bright regions inside network boundaries of coronal holes. We find that these network jets appear not only in coronal holes but also in quiet-sun regions. Using IRIS 1330A (C II) slit-jaw images, we extract several parameters of these network jets, e.g. apparent speed, length, lifetime and increase in foot-point brightness. Using several observations, we find that some properties of the jets are very similar but others are obviously different between the quiet sun and coronal holes. For example, our study shows that the coronal-hole jets appear to be faster and longer than those in the quiet sun. This can be directly attributed to a difference in the magnetic configuration of the two regions with open magnetic field lines rooted in coronal holes and magnetic loops often present in quiet sun. We have also detected compact bright loops, likely transition region loops, mostly in quiet sun. These small loop-like regions are generally devoid of network jets. In spite of different magnetic structures in the coronal hole and quiet sun in the transition region, there appears to be no substantial difference for the increase in foot-point brightness of the jets, which suggests that the generation mechanism of these network jets is likely the same in both regions., Comment: 16 pages, 6 Figures, Solar Physics (in press)

Observations of debris disks offer a window into the physical and dynamical properties of planetesimals in extrasolar systems through the size distribution of dust grains. In particular, the millimeter spectral index of thermal dust emission encodes information on the grain size distribution. We have made new VLA observations of a sample of seven nearby debris disks at 9 mm, with 3" resolution and $\sim5$ $\mu$Jy/beam rms. We combine these with archival ATCA observations of eight additional debris disks observed at 7 mm, together with up-to-date observations of all disks at (sub)millimeter wavelengths from the literature to place tight constraints on the millimeter spectral indices and thus grain size distributions. The analysis gives a weighted mean for the slope of the power law grain size distribution, $n(a)\propto a^{-q}$, of $\langle q \rangle = 3.36\pm0.02$, with a possible trend of decreasing $q$ for later spectral type stars. We compare our results to a range of theoretical models of collisional cascades, from the standard self-similar, steady-state size distribution ($q=3.5$) to solutions that incorporate more realistic physics such as alternative velocity distributions and material strengths, the possibility of a cutoff at small dust sizes from radiation pressure, as well as results from detailed dynamical calculations of specific disks. Such effects can lead to size distributions consistent with the data, and plausibly the observed scatter in spectral indices. For the AU Mic system, the VLA observations show clear evidence of a highly variable stellar emission component; this stellar activity obviates the need to invoke the presence of an asteroid belt to explain the previously reported compact millimeter source in this system., Comment: 34 pages, 5 figures

The solar chromosphere contains thin, highly dynamic strands of plasma known as spicules. Recently, it has been suggested that the smallest and fastest (Type II) spicules are identical to intermittent jets observed by the Interface Region Imaging Spectrograph. These jets appear to expand out along open magnetic field lines rooted in unipolar network regions of coronal holes. In this paper we revisit a thirty-year-old idea that spicules may be caused by upward forces associated with Alfven waves. These forces involve the conversion of transverse Alfven waves into compressive acoustic-like waves that steepen into shocks. The repeated buffeting due to upward shock propagation causes nonthermal expansion of the chromosphere and a transient levitation of the transition region. Some older models of wave-driven spicules assumed sinusoidal wave inputs, but the solar atmosphere is highly turbulent and stochastic. Thus, we model this process using the output of a time-dependent simulation of reduced magnetohydrodynamic turbulence. The resulting mode-converted compressive waves are strongly variable in time, with a higher transition region occurring when the amplitudes are large and a lower transition region when the amplitudes are small. In this picture, the transition region bobs up and down by several Mm on timescales less than a minute. These motions produce narrow, intermittent extensions of the chromosphere that have similar properties as the observed jets and Type II spicules., Comment: Accepted for publication in the Astrophysical Journal. 8 pages (emulateapj style), 5 figures

We investigate several key questions of plasma heating in open-field regions of the corona that connect to the solar wind. We present results for a model of Alfven-wave-driven turbulence for three typical open magnetic field structures: a polar coronal hole, an open flux tube neighboring an equatorial streamer, and an open flux tube near a strong-field active region. We compare time-steady, one-dimensional turbulent heating models (Cranmer et al., 2007) against fully time-dependent three-dimensional reduced-magnetohydrodynamics modeling of BRAID (van Ballegooijen et al., 2011). We find that the time-steady results agree well with time-averaged results from BRAID. The time-dependence allows us to investigate the variability of the magnetic fluctuations and of the heating in the corona. The high-frequency tail of the power spectrum of fluctuations forms a power law whose exponent varies with height, and we discuss the possible physical explanation for this behavior. The variability in the heating rate is bursty and nanoflare-like in nature, and we analyze the amount of energy lost via dissipative heating in transient events throughout the simulation. The average energy in these events is 10^21.91 erg, within the "picoflare" range, and many events reach classical "nanoflare" energies. We also estimated the multithermal distribution of temperatures that would result from the heating-rate variability, and found good agreement with observed widths of coronal differential emission measure (DEM) distributions. The results of the modeling presented in this paper provide compelling evidence that turbulent heating in the solar atmosphere by Alfven waves accelerates the solar wind in open flux tubes., Comment: 13 pages, 15 figures; accepted to The Astrophysical Journal

There have been several ideas proposed to explain how the Sun's corona is heated and how the solar wind is accelerated. Some models assume that open magnetic field lines are heated by Alfven waves driven by photospheric motions and dissipated after undergoing a turbulent cascade. Other models posit that much of the solar wind's mass and energy is injected via magnetic reconnection from closed coronal loops. The latter idea is motivated by observations of reconnecting jets and also by similarities of ion composition between closed loops and the slow wind. Wave/turbulence models have also succeeded in reproducing observed trends in ion composition signatures versus wind speed. However, the absolute values of the charge-state ratios predicted by those models tended to be too low in comparison with observations. This letter refines these predictions by taking better account of weak Coulomb collisions for coronal electrons, whose thermodynamic properties determine the ion charge states in the low corona. A perturbative description of nonlocal electron transport is applied to an existing set of wave/turbulence models. The resulting electron velocity distributions in the low corona exhibit mild suprathermal tails characterized by "kappa" exponents between 10 and 25. These suprathermal electrons are found to be sufficiently energetic to enhance the charge states of oxygen ions, while maintaining the same relative trend with wind speed that was found when the distribution was assumed to be Maxwellian. The updated wave/turbulence models are in excellent agreement with solar wind ion composition measurements., Comment: Accepted for publication in Astrophysical Journal Letters. 5 pages (emulateapj style), 3 figures

Protons in the solar corona and heliosphere exhibit anisotropic velocity distributions, violation of magnetic moment conservation, and a general lack of thermal equilibrium with the other particle species. There is no agreement about the identity of the physical processes that energize non-Maxwellian protons in the solar wind, but a traditional favorite has been the dissipation of ion cyclotron resonant Alfven waves. This paper presents kinetic models of how ion cyclotron waves heat protons on their journey from the corona to interplanetary space. It also derives a wide range of new solutions for the relevant dispersion relations, marginal stability boundaries, and nonresonant velocity-space diffusion rates. A phenomenological model containing both cyclotron damping and turbulent cascade is constructed to explain the suppression of proton heating at low alpha-proton differential flow speeds. These effects are implemented in a large-scale model of proton thermal evolution from the corona to 1 AU. A Monte Carlo ensemble of realistic wind speeds, densities, magnetic field strengths, and heating rates produces a filled region of parameter space (in a plane described by the parallel plasma beta and the proton temperature anisotropy ratio) similar to what is measured. The high-beta edges of this filled region are governed by plasma instabilities and strong heating rates. The low-beta edges correspond to weaker proton heating and a range of relative contributions from cyclotron resonance. On balance, the models are consistent with other studies that find only a small fraction of the turbulent power spectrum needs to consist of ion cyclotron waves., Comment: Accepted for publication in the Astrophysical Journal Supplement Series. 26 pages (emulateapj style), 18 figures

Forecasting models of the solar wind often rely on simple parameterizations of the magnetic field that ignore the effects of the full magnetic field geometry. In this paper, we present the results of two solar wind prediction models that consider the full magnetic field profile and include the effects of Alfv\'en waves on coronal heating and wind acceleration. The one-dimensional MHD code ZEPHYR self-consistently finds solar wind solutions without the need for empirical heating functions. Another 1D code, introduced in this paper (The Efficient Modified-Parker-Equation-Solving Tool, TEMPEST), can act as a smaller, stand-alone code for use in forecasting pipelines. TEMPEST is written in Python and will become a publicly available library of functions that is easy to adapt and expand. We discuss important relations between the magnetic field profile and properties of the solar wind that can be used to independently validate prediction models. ZEPHYR provides the foundation and calibration for TEMPEST, and ultimately we will use these models to predict observations and explain space weather created by the bulk solar wind. We are able to reproduce with both models the general anticorrelation seen in comparisons of observed wind speed at 1 AU and the flux tube expansion factor. There is significantly less spread than comparing the results of the two models than between ZEPHYR and a traditional flux tube expansion relation. We suggest that the new code, TEMPEST, will become a valuable tool in the forecasting of space weather., Comment: 12 pages, 11 figures; accepted to ApJ

A large fraction of cool, low-mass stars exhibit brightness fluctuations that arise from a combination of convective granulation, acoustic oscillations, magnetic activity, and stellar rotation. Much of the short-timescale variability takes the form of stochastic noise, whose presence may limit the progress of extrasolar planet detection and characterization. In order to lay the groundwork for extracting useful information from these quasi-random signals, we focus on the origin of the granulation-driven component of the variability. We apply existing theoretical scaling relations to predict the star-integrated variability amplitudes for 508 stars with photometric light curves measured by the Kepler mission. We also derive an empirical correction factor that aims to account for the suppression of convection in F-dwarf stars with magnetic activity and shallow convection zones. So that we can make predictions of specific observational quantities, we performed Monte Carlo simulations of granulation light curves using a Lorentzian power spectrum. These simulations allowed us to reproduce the so-called "flicker floor" (i.e., a lower bound in the relationship between the full light-curve range and power in short-timescale fluctuations) that was found in the Kepler data. The Monte Carlo model also enabled us to convert the modeled fluctuation variance into a flicker amplitude directly comparable with observations. When the magnetic suppression factor described above is applied, the model reproduces the observed correlation between stellar surface gravity and flicker amplitude. Observationally validated models like these provide new and complementary evidence for a possible impact of magnetic activity on the properties of near-surface convection., Comment: Accepted for publication in the Astrophysical Journal. 7 pages (emulateapj style), 5 figures

Many low-mass pre-main-sequence stars exhibit strong magnetic activity and coronal X-ray emission. Even after the primordial accretion disk has been cleared out, the star's high-energy radiation continues to affect the formation and evolution of dust, planetesimals, and large planets. Young stars with debris disks are thus ideal environments for studying the earliest stages of non-accretion-driven coronae. In this paper we simulate the corona of AU Mic, a nearby active M dwarf with an edge-on debris disk. We apply a self-consistent model of coronal loop heating that was derived from numerical simulations of solar field-line tangling and magnetohydrodynamic turbulence. We also synthesize the modeled star's X-ray luminosity and thermal radio/millimeter continuum emission. A realistic set of parameter choices for AU Mic produces simulated observations that agree with all existing measurements and upper limits. This coronal model thus represents an alternative explanation for a recently discovered ALMA central emission peak that was suggested to be the result of an inner "asteroid belt" within 3 AU of the star. However, it is also possible that the central 1.3 mm peak is caused by a combination of active coronal emission and a bright inner source of dusty debris. Additional observations of this source's spatial extent and spectral energy distribution at millimeter and radio wavelengths will better constrain the relative contributions of the proposed mechanisms., Comment: Accepted for publication in the Astrophysical Journal. 8 pages (emulateapj style), 4 figures

The solar wind is connected to the Sun's atmosphere by flux tubes that are rooted in an ever-changing pattern of positive and negative magnetic polarities on the surface. Observations indicate that the magnetic field is filamentary and intermittent across a wide range of spatial scales. However, we do not know to what extent the complex flux tube topology seen near the Sun survives as the wind expands into interplanetary space. In order to study the possible long-distance connections between the corona and the heliosphere, we developed new models of turbulence-driven solar wind acceleration along empirically constrained field lines. We used a potential-field model of the Quiet Sun to trace field lines into the ecliptic plane with unprecedented spatial resolution at their footpoints. For each flux tube, a one-dimensional model was created with an existing wave/turbulence code that solves equations of mass, momentum, and energy conservation from the photosphere to 4 AU. To take account of stream-stream interactions between flux tubes, we used those models as inner boundary conditions for a time-steady MHD description of radial and longitudinal structure in the ecliptic. Corotating stream interactions smear out much of the smallest-scale variability, making it difficult to see how individual flux tubes on granular or supergranular scales can survive out to 1 AU. However, our models help clarify the level of "background" variability with which waves and turbulent eddies should be expected to interact. Also, the modeled fluctuations in magnetic field magnitude were seen to match measured power spectra quite well., Comment: Accepted for publication in the Astrophysical Journal. 18 pages (emulateapj style), 17 figures

In the parts of the solar corona and solar wind that experience the fewest Coulomb collisions, the component proton, electron, and heavy ion populations are not in thermal equilibrium with one another. Observed differences in temperatures, outflow speeds, and velocity distribution anisotropies are useful constraints on proposed explanations for how the plasma is heated and accelerated. This paper presents new predictions of the rates of collisionless heating for each particle species, in which the energy input is assumed to come from magnetohydrodynamic (MHD) turbulence. We first created an empirical description of the radial evolution of Alfven, fast-mode, and slow-mode MHD waves. This model provides the total wave power in each mode as a function of distance along an expanding flux tube in the high-speed solar wind. Next we solved a set of cascade advection-diffusion equations that give the time-steady wavenumber spectra at each distance. An approximate term for nonlinear coupling between the Alfven and fast-mode fluctuations is included. For reasonable choices of the parameters, our model contains enough energy transfer from the fast mode to the Alfven mode to excite the high-frequency ion cyclotron resonance. This resonance is efficient at heating protons and other ions in the direction perpendicular to the background magnetic field, and our model predicts heating rates for these species that agree well with both spectroscopic and in situ measurements. Nonetheless, the high-frequency waves comprise only a small part of the total Alfvenic fluctuation spectrum, which remains highly two-dimensional as is observed in interplanetary space., Comment: Accepted for publication in the Astrophysical Journal. 30 pages (emulateapj style), 18 figures

The basic mechanisms responsible for producing winds from cool, late-type stars are still largely unknown. We take inspiration from recent progress in understanding solar wind acceleration to develop a physically motivated model of the time-steady mass loss rates of cool main-sequence stars and evolved giants. This model follows the energy flux of magnetohydrodynamic turbulence from a subsurface convection zone to its eventual dissipation and escape through open magnetic flux tubes. We show how Alfven waves and turbulence can produce winds in either a hot corona or a cool extended chromosphere, and we specify the conditions that determine whether or not coronal heating occurs. These models do not utilize arbitrary normalization factors, but instead predict the mass loss rate directly from a star's fundamental properties. We take account of stellar magnetic activity by extending standard age-activity-rotation indicators to include the evolution of the filling factor of strong photospheric magnetic fields. We compared the predicted mass loss rates with observed values for 47 stars and found significantly better agreement than was obtained from the popular scaling laws of Reimers, Schroeder, and Cuntz. The algorithm used to compute cool-star mass loss rates is provided as a self-contained and efficient computer code. We anticipate that the results from this kind of model can be incorporated straightforwardly into stellar evolution calculations and population synthesis techniques., Comment: 23 pages (emulateapj style), 14 figures, ApJ, in press. A brief IDL subroutine that implements the model described in this paper will be distributed as "online-only material," and this code is also available at http://www.cfa.harvard.edu/~scranmer/cranmer_data.html

The physical processes that heat the solar corona and accelerate the solar wind remain unknown after many years of study. Some have suggested that the wind is driven by waves and turbulence in open magnetic flux tubes, and others have suggested that plasma is injected into the open tubes by magnetic reconnection with closed loops. In order to test the latter idea, we developed Monte Carlo simulations of the photospheric "magnetic carpet" and extrapolated the time-varying coronal field. These models were constructed for a range of different magnetic flux imbalance ratios. Completely balanced models represent quiet regions on the Sun and source regions of slow solar wind streams. Highly imbalanced models represent coronal holes and source regions of fast wind streams. The models agree with observed emergence rates, surface flux densities, and number distributions of magnetic elements. Despite having no imposed supergranular motions, a realistic network of magnetic "funnels" appeared spontaneously. We computed the rate at which closed field lines open up (i.e., recycling times for open flux), and we estimated the energy flux released in reconnection events involving the opening up of closed flux tubes. For quiet regions and mixed-polarity coronal holes, these energy fluxes were found to be much lower than required to accelerate the solar wind. For the most imbalanced coronal holes, the energy fluxes may be large enough to power the solar wind, but the recycling times are far longer than the time it takes the solar wind to accelerate into the low corona. Thus, it is unlikely that either the slow or fast solar wind is driven by reconnection and loop-opening processes in the magnetic carpet., Comment: 25 pages (emulateapj style), 13 figures, ApJ, in press