Showing 1–36 of 36 results for author: Hahn, M

Observations of Non-thermal Velocity and Comparison with Alfvén Wave Turbulence Model in Solar Active Regions

Authors: M. Asgari-Targhi, D. H. Brooks, M. Hahn, S. Imada, E. Tajfirouze, D. W. Savin

Abstract: We present a study of spectral line width measurements from the Extreme Ultraviolet Imaging Spectrometer (EIS) on {\it Hinode}. We used spectral line profiles of Fe {\sc xvi} 262.984 Å, Fe {\sc xiv} 264.787 Å, Fe {\sc xiv} 270.519 Å, Fe {\sc xiv} 274.203 Å, and Fe {\sc xv} 284.160 Å, and studied 11 active regions. Previous studies of spectral line widths have shown that in hot loops in the cores o… ▽ More We present a study of spectral line width measurements from the Extreme Ultraviolet Imaging Spectrometer (EIS) on {\it Hinode}. We used spectral line profiles of Fe {\sc xvi} 262.984 Å, Fe {\sc xiv} 264.787 Å, Fe {\sc xiv} 270.519 Å, Fe {\sc xiv} 274.203 Å, and Fe {\sc xv} 284.160 Å, and studied 11 active regions. Previous studies of spectral line widths have shown that in hot loops in the cores of active regions, the observed non-thermal velocities are smaller than predicted from models of reconnection jets in the corona or shock heating associated with Alfvén waves. The observed line widths are also inconsistent with models of chromospheric evaporation due to coronal nanoflares. We show that recent advances in higher resolution Alfvén wave turbulence modeling enables us to obtain non-thermal velocities similar to those measured in active regions. The observed non-thermal velocities for the 11 active regions in our study are in the range of 17$-$30 $\rm km ~ s^{-1}$, consistent with the spectral line non-thermal widths predicted from our model of 16 interacting flux tubes, which are in the range of ~15$-$37 $\rm km ~ s^{-1}$. △ Less

Submitted 25 April, 2024; originally announced April 2024.

Experimental study of Alfvén wave reflection from an Alfvén-speed gradient relevant to the solar coronal holes

Authors: Sayak Bose, Jason M. TenBarge, Troy Carter, Michael Hahn, Hantao Ji, James Juno, Daniel Wolf Savin, Shreekrishna Tripathi, Stephen Vincena

Abstract: We report the first experimental detection of a reflected Alfvén wave from an Alfvén-speed gradient under conditions similar to those in coronal holes. The experiments were conducted in the Large Plasma Device at the University of California, Los Angeles. We present the experimentally measured dependence of the coefficient of reflection versus the wave inhomogeneity parameter, i.e., the ratio of t… ▽ More We report the first experimental detection of a reflected Alfvén wave from an Alfvén-speed gradient under conditions similar to those in coronal holes. The experiments were conducted in the Large Plasma Device at the University of California, Los Angeles. We present the experimentally measured dependence of the coefficient of reflection versus the wave inhomogeneity parameter, i.e., the ratio of the wave length of the incident wave to the length scale of the gradient. Two-fluid simulations using the Gkeyll code qualitatively agree with and support the experimental findings. Our experimental results support models of wave heating that rely on wave reflection at low heights from a smooth Alfvén-speed gradient to drive turbulence. △ Less

Submitted 9 February, 2024; originally announced February 2024.

Are Non-thermal Velocities in Active Region Coronal Loops Anisotropic?

Authors: Michael Hahn, Mahboubeh Asgari-Targhi, Daniel Wolf Savin

Abstract: We have measured line widths in active region coronal loops in order to determine whether the non-thermal broadening is anisotropic with respect to the magnetic field direction. These non-thermal velocities are caused by unresolved fluid motions. Our analysis method combines spectroscopic data and a magnetic field extrapolation. We analyzed spectra from the Extreme Ultraviolet Imaging Spectrometer… ▽ More We have measured line widths in active region coronal loops in order to determine whether the non-thermal broadening is anisotropic with respect to the magnetic field direction. These non-thermal velocities are caused by unresolved fluid motions. Our analysis method combines spectroscopic data and a magnetic field extrapolation. We analyzed spectra from the Extreme Ultraviolet Imaging Spectrometer on Hinode. A differential emission measure analysis showed that many spectral lines that are commonly considered to be formed in the active region have a substantial contribution from the background quiet Sun. From these spectra we identified lines whose emission was dominated by the active region loops rather than background sources. Using these lines, we constructed maps of the non-thermal velocity. With data from the Helioseismic Magnetic Imager on the Solar Dynamics Observatory and the Coronal Modeling System nonlinear force-free magnetic field reconstruction code, we traced several of the magnetic field lines through the active region. Comparing the spectroscopic and magnetic data, we looked for correlations of non-thermal velocity with the viewing angle between the line of sight and the magnetic field. We found that non-thermal velocities show a weak anti-correlation with the viewing angle. That is, the tendency is for the non-thermal velocity to be slightly larger in the parallel direction. This parallel broadening may be due to acoustic waves or unresolved parallel flows. △ Less

Submitted 5 May, 2023; originally announced May 2023.

Comments: Submitted to the Astrophysical Journal

Deriving instrumental point spread functions from partially occulted images

Authors: Stefan Johann Hofmeister, Michael Hahn, Daniel Wolf Savin

Abstract: The point-spread function (PSF) of an imaging system describes the response of the system to a point source. Accurately determining the PSF enables one to correct for the combined effects of focussing and scattering within the imaging system, and thereby enhance the spatial resolution and dynamic contrast of the resulting images. We present a semi-empirical semi-blind methodology to derive a PSF f… ▽ More The point-spread function (PSF) of an imaging system describes the response of the system to a point source. Accurately determining the PSF enables one to correct for the combined effects of focussing and scattering within the imaging system, and thereby enhance the spatial resolution and dynamic contrast of the resulting images. We present a semi-empirical semi-blind methodology to derive a PSF from partially occulted images. We partition the two-dimensional PSF into multiple segments, set up a multi-linear system of equations, and directly fit the system of equations to determine the PSF weight in each segment. The algorithm is guaranteed to converge towards the correct instrumental PSF for a large class of occultations, does not require a predefined functional form of the PSF, and can be applied to a large variety of partially occulted images, such as within laboratory settings, regular calibrations within a production line or in the field, astronomical images of distant clusters of stars, or partial solar eclipse images. We show that the central weight of the PSF, which gives the percentage of photons that are not scattered by the instrument, is accurate to bettern than 1.2%. The mean absolute percentage error between the reconstructed and true PSF is usually between 0.5% and 5% for the entire PSF, between 0.5% and 5% for the PSF core, and between 0.5% and 3% for the PSF tail. △ Less

Submitted 28 October, 2022; v1 submitted 6 July, 2022; originally announced July 2022.

Evidence for Parameteric Decay Instability in the Lower Solar Atmosphere

Authors: Michael Hahn, Xiangrong Fu, Daniel Wolf Savin

Abstract: We find evidence for the first observation of the parametric decay instability (PDI) in the lower solar atmosphere. Specifically, we find that the power spectrum of density fluctuations near the solar transition region resembles the power spectrum of the velocity fluctuations, but with the frequency axis scaled up by about a factor of two. These results are from an analysis of the Si IV lines obse… ▽ More We find evidence for the first observation of the parametric decay instability (PDI) in the lower solar atmosphere. Specifically, we find that the power spectrum of density fluctuations near the solar transition region resembles the power spectrum of the velocity fluctuations, but with the frequency axis scaled up by about a factor of two. These results are from an analysis of the Si IV lines observed by the Interface Region Imaging Spectrometer (IRIS) in the transition region of a polar coronal hole. We also find that the density fluctuations have radial velocity of about 75 km/s and that the velocity fluctuations are much faster with an estimated speed of 250 km/s, as is expected for sound waves and Alfvén waves, respectively, in the transition region. Theoretical calculations show that this frequency relationship is consistent with those expected from PDI for the plasma conditions of the observed region. These measurements suggest an interaction between sound waves and Alfvén waves in the transition region that is evidence for the parametric decay instability. △ Less

Submitted 20 April, 2022; originally announced April 2022.

Comments: Submitted to the Astrophysical Journal

Critical Science Plan for the Daniel K. Inouye Solar Telescope (DKIST)

Authors: Mark P. Rast, Nazaret Bello González, Luis Bellot Rubio, Wenda Cao, Gianna Cauzzi, Edward DeLuca, Bart De Pontieu, Lyndsay Fletcher, Sarah E. Gibson, Philip G. Judge, Yukio Katsukawa, Maria D. Kazachenko, Elena Khomenko, Enrico Landi, Valentin Martínez Pillet, Gordon J. D. Petrie, Jiong Qiu, Laurel A. Rachmeler, Matthias Rempel, Wolfgang Schmidt, Eamon Scullion, Xudong Sun, Brian T. Welsch, Vincenzo Andretta, Patrick Antolin , et al. (62 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 20 August, 2020; v1 submitted 18 August, 2020; originally announced August 2020.

Initial results from the New Horizons exploration of 2014 MU69, a small Kuiper Belt Object

Authors: S. A. Stern, H. A. Weaver, J. R. Spencer, C. B. Olkin, G. R. Gladstone, W. M. Grundy, J. M. Moore, D. P. Cruikshank, H. A. Elliott, W. B. McKinnon, J. Wm. Parker, A. J. Verbiscer, L. A. Young, D. A. Aguilar, J. M. Albers, T. Andert, J. P. Andrews, F. Bagenal, M. E. Banks, B. A. Bauer, J. A. Bauman, K. E. Bechtold, C. B. Beddingfield, N. Behrooz, K. B. Beisser , et al. (180 additional authors not shown)

Abstract: The Kuiper Belt is a distant region of the Solar System. On 1 January 2019, the New Horizons spacecraft flew close to (486958) 2014 MU69, a Cold Classical Kuiper Belt Object, a class of objects that have never been heated by the Sun and are therefore well preserved since their formation. Here we describe initial results from these encounter observations. MU69 is a bi-lobed contact binary with a fl… ▽ More The Kuiper Belt is a distant region of the Solar System. On 1 January 2019, the New Horizons spacecraft flew close to (486958) 2014 MU69, a Cold Classical Kuiper Belt Object, a class of objects that have never been heated by the Sun and are therefore well preserved since their formation. Here we describe initial results from these encounter observations. MU69 is a bi-lobed contact binary with a flattened shape, discrete geological units, and noticeable albedo heterogeneity. However, there is little surface color and compositional heterogeneity. No evidence for satellites, ring or dust structures, gas coma, or solar wind interactions was detected. By origin MU69 appears consistent with pebble cloud collapse followed by a low velocity merger of its two lobes. △ Less

Submitted 2 April, 2020; originally announced April 2020.

Comments: 43 pages, 8 figure

Journal ref: Science 364, eaaw9771 (2019)

Laboratory Calibrations of Fe XII-XIV Line-Intensity Ratios for Electron Density Diagnostics

Authors: Thusitha P. Arthanayaka, Peter Beiersdorfer, Gregory V. Brown, Ming Feng Gu, Michael Hahn, Natalie Hell, Tom E. Lockard, Daniel Wolf Savin

Abstract: We have used an electron beam ion trap to measure electron-density-diagnostic line-intensity ratios for extreme ultraviolet lines from F XII, XIII, and XIV at wavelengths of 185-205 255-276 Angstroms. These ratios can be used as density diagnostics for astrophysical spectra and are especially relevant to solar physics. We found that density diagnostics using the Fe XIII 196.53/202.04 and the Fe XI… ▽ More We have used an electron beam ion trap to measure electron-density-diagnostic line-intensity ratios for extreme ultraviolet lines from F XII, XIII, and XIV at wavelengths of 185-205 255-276 Angstroms. These ratios can be used as density diagnostics for astrophysical spectra and are especially relevant to solar physics. We found that density diagnostics using the Fe XIII 196.53/202.04 and the Fe XIV 264.79/274.21 and 270.52A/274.21 line ratios are reliable using the atomic data calculated with the Flexible Atomic Code. On the other hand, we found a large discrepancy between the FAC theory and experiment for the commonly used Fe XII (186.85 + 186.88)/195.12 line ratio. These FAC theory calculations give similar results to the data tabulated in CHIANTI, which are commonly used to analyze solar observations. Our results suggest that the discrepancies seen between solar coronal density measurements using the Fe XII (186.85 + 186.88)/195.12 and Fe XIII 196.54/202.04 line ratios are likely due to issues with the atomic calculations for Fe XII. △ Less

Submitted 21 August, 2019; originally announced August 2019.

Measured reduction in Alfvén wave energy propagating through longitudinal gradients scaled to match solar coronal holes

Authors: Sayak Bose, Troy Carter, Michael Hahn, Shreekrishna Tripathi, Stephen Vincena, Daniel Wolf Savin

Abstract: We have explored the effectiveness of a longitudinal gradient in Alfvén speed in reducing the energy of propagating Alfvén waves under conditions scaled to match solar coronal holes. The experiments were conducted in the Large Plasma Device at the University of California, Los Angeles. Our results show that the energy of the transmitted Alfvén wave decreases as the inhomogeneity parameter,… ▽ More We have explored the effectiveness of a longitudinal gradient in Alfvén speed in reducing the energy of propagating Alfvén waves under conditions scaled to match solar coronal holes. The experiments were conducted in the Large Plasma Device at the University of California, Los Angeles. Our results show that the energy of the transmitted Alfvén wave decreases as the inhomogeneity parameter, $λ/L_{\rm A}$, increases. Here, $λ$ is the wavelength of the Alfvén wave and $L_{\rm A}$ is the scale length of Alfvén speed gradient. For gradients similar to those in coronal holes, the waves are observed to lose a factor of $\approx 5$ more energy than they do when propagating through a uniform plasma without a gradient. We have carried out further experiments and analyses to constrain the cause of wave energy reduction in the gradient. The loss of Alfvén wave energy from mode coupling is unlikely, as we have not detected any other modes. Contrary to theoretical expectations, the reduction in the energy of the transmitted wave is not accompanied by a detectable reflected wave. Nonlinear effects are ruled out as the amplitude of the initial wave is too small and the wave frequency well below the ion cyclotron frequency. Since the total energy must be conserved, it is possible that the lost wave energy is being deposited in the plasma. Further studies are needed to explore where the energy is going. △ Less

Submitted 4 July, 2019; v1 submitted 26 April, 2019; originally announced April 2019.

Journal ref: ApJ 882 183 (2019)

Perspectives on Astrophysics Based on Atomic, Molecular, and Optical (AMO) Techniques

Authors: Daniel Wolf Savin, James F. Babb, Paul M. Bellan, Crystal Brogan, Jan Cami, Paola Caselli, Lia Corrales, Gerardo Dominguez, Steven R. Federman, Chris J. Fontes, Richard Freedman, Brad Gibson, Leon Golub, Thomas W. Gorczyca, Michael Hahn, Sarah M. Hörst, Reggie L. Hudson, Jeffrey Kuhn, James E. Lawler, Maurice A. Leutenegger, Joan P. Marler, Michael C. McCarthy, Brett A. McGuire, Stefanie N. Milam, Nicholas A. Murphy , et al. (13 additional authors not shown)

Abstract: About two generations ago, a large part of AMO science was dominated by experimental high energy collision studies and perturbative theoretical methods. Since then, AMO science has undergone a transition and is now dominated by quantum, ultracold, and ultrafast studies. But in the process, the field has passed over the complexity that lies between these two extremes. Most of the Universe resides i… ▽ More About two generations ago, a large part of AMO science was dominated by experimental high energy collision studies and perturbative theoretical methods. Since then, AMO science has undergone a transition and is now dominated by quantum, ultracold, and ultrafast studies. But in the process, the field has passed over the complexity that lies between these two extremes. Most of the Universe resides in this intermediate region. We put forward that the next frontier for AMO science is to explore the AMO complexity that describes most of the Cosmos. △ Less

Submitted 14 November, 2018; originally announced November 2018.

Comments: White paper submission to the Decadal Assessment and Outlook Report on Atomic, Molecular, and Optical (AMO) Science (AMO 2020)

Density Fluctuations in a Polar Coronal Hole

Authors: Michael Hahn, Elke D'Huys, Daniel Wolf Savin

Abstract: We have measured the root-mean-square (RMS) amplitude of intensity fluctuations, $ΔI$, in plume and interplume regions of a polar coronal hole. These intensity fluctuations correspond to density fluctuations. Using data from the \rev{Sun Watcher using Active Pixel System detector and Image Processing} (SWAP) on \textit{Project for Onboard Autonomy (Proba2)}, our results extend up to a height of ab… ▽ More We have measured the root-mean-square (RMS) amplitude of intensity fluctuations, $ΔI$, in plume and interplume regions of a polar coronal hole. These intensity fluctuations correspond to density fluctuations. Using data from the \rev{Sun Watcher using Active Pixel System detector and Image Processing} (SWAP) on \textit{Project for Onboard Autonomy (Proba2)}, our results extend up to a height of about 1.35~$R_{\sun}$. One advantage of the RMS analysis is that it does not rely on a detailed evaluation of the power spectrum, which is limited by noise levels to low heights in the corona. The RMS approach can be performed up to larger heights where the noise level is greater, provided that the noise itself can be quantified. At low heights, both the absolute $ΔI$, and the amplitude relative to the mean intensity, $ΔI/I$, decrease with height. However, starting at about 1.2~$R_{\sun}$ $ΔI/I$ increases, reaching 20--40\% by 1.35~$R_{\sun}$. This corresponds to density fluctuations of $Δn_{\mathrm{e}}/n_{\mathrm{e}} \approx$~10--20\%. The increasing relative amplitude implies that the density fluctuations are generated in the corona itself. One possibility is that the density fluctuations are generated by an instability of Alfvén waves. This generation mechanism is consistent with some theoretical models and with observations of Alfvén wave amplitudes in coronal holes. Although we find that the energy of the observed density fluctuations is small, these fluctuations are likely to play an important indirect role in coronal heating by promoting the reflection of Alfvén waves and driving turbulence. △ Less

Submitted 26 April, 2018; originally announced April 2018.

Comments: Accepted for Astrophysical Journal

Spectroscopic Measurements of the Ion Velocity Distribution at the Base of the Fast Solar Wind

Authors: Natasha L. S. Jeffrey, Michael Hahn, Daniel W. Savin, Lyndsay Fletcher

Abstract: In situ measurements of the fast solar wind reveal non-thermal distributions of electrons, protons and, minor ions extending from $0.3$ AU to the heliopause. The physical mechanisms responsible for these non-thermal properties and the location where these properties originate remain open questions. Here we present spectroscopic evidence, from extreme ultraviolet spectroscopy, that the velocity dis… ▽ More In situ measurements of the fast solar wind reveal non-thermal distributions of electrons, protons and, minor ions extending from $0.3$ AU to the heliopause. The physical mechanisms responsible for these non-thermal properties and the location where these properties originate remain open questions. Here we present spectroscopic evidence, from extreme ultraviolet spectroscopy, that the velocity distribution functions (VDFs) of minor ions are already non-Gaussian at the base of the fast solar wind in a coronal hole, at altitudes of $< 1.1 R_{\odot}$. Analysis of Fe, Si, and Mg spectral lines reveal a peaked line-shape core and broad wings that can be characteristed by a kappa VDF. A kappa distribution fit gives very small kappa indices off-limb of $κ\approx1.9-2.5$, indicating either (a) ion populations far from thermal equilibrium, (b) fluid motions such as non-Gaussian turbulent fluctuations or non-uniform wave motions, or (c) some combination of both. These observations provide important empirical constraints for the source region of the fast solar wind and for the theoretical models of the different acceleration, heating, and energy deposition processes therein. To the best of our knowledge, this is the first time that the ion VDF in the fast solar wind has been probed so close to its source region. The findings are also a timely precursor to the upcoming 2018 launch of the Parker Solar Probe, which will provide the closest in situ measurements of the solar wind at approximately $0.04$ AU ($8.5$ solar radii). △ Less

Submitted 20 February, 2018; originally announced February 2018.

Comments: 9 pages, 4 figures, accepted for publication in the Astrophysical Journal Letters

Electron-Impact Multiple Ionization Cross Sections for Atoms and Ions of Helium through Zinc

Authors: Michael Hahn, Alfred Mueller, Daniel Wolf Savin

Abstract: We have compiled a set of electron-impact multiple ionization (EIMI) cross sections for astrophysically relevant ions. EIMI can have a significant effect on the ionization balance of non-equilibrium plasmas. For example, it can be important if there is a rapid change in the electron temperature or if there is a non-thermal electron energy distribution, such as a kappa distribution. Cross sections… ▽ More We have compiled a set of electron-impact multiple ionization (EIMI) cross sections for astrophysically relevant ions. EIMI can have a significant effect on the ionization balance of non-equilibrium plasmas. For example, it can be important if there is a rapid change in the electron temperature or if there is a non-thermal electron energy distribution, such as a kappa distribution. Cross sections for EIMI are needed in order to account for these processes in plasma modeling and for spectroscopic interpretation. Here, we describe our comparison of proposed semiempirical formulae to the available experimental EIMI cross section data. Based on this comparison, we have interpolated and extrapolated fitting parameters to systems that have not yet been measured. A tabulation of the fit parameters is provided for 3466 EIMI cross sections. We also highlight some outstanding issues that remain to be resolved. △ Less

Submitted 7 August, 2017; originally announced August 2017.

Comments: Submitted to Astrophysical Journal Supplement. The cross section database (Table 2 in the manuscript) is available upon request

Inferring the Coronal Density Irregularity from EUV Spectra

Authors: Michael Hahn, Daniel Wolf Savin

Abstract: Understanding the density structure of the solar corona is important for modeling both coronal heating and the solar wind. Direct measurements are difficult because of line-of-sight integration and possible unresolved structures. We present a new method for quantifying such structure using density-sensitive EUV line intensities to derive a density irregularity parameter, a relative measure of the… ▽ More Understanding the density structure of the solar corona is important for modeling both coronal heating and the solar wind. Direct measurements are difficult because of line-of-sight integration and possible unresolved structures. We present a new method for quantifying such structure using density-sensitive EUV line intensities to derive a density irregularity parameter, a relative measure of the amount of structure along the line of sight. We also present a simple model to relate the inferred irregularities to physical quantities, such as the filling factor and density contrast. For quiet Sun regions and interplume regions of coronal holes, we find a density contrast of at least a factor of three to ten and corresponding filling factors of about 10-20%. Our results are in rough agreement with other estimates of the density structures in these regions. The irregularity diagnostic provides a useful relative measure of unresolved structure in various regions of the corona. △ Less

Submitted 20 June, 2016; originally announced June 2016.

Comments: Submitted to the Astrophysical Journal

The Pluto system: Initial results from its exploration by New Horizons

Authors: S. A. Stern, F. Bagenal, K. Ennico, G. R. Gladstone, W. M. Grundy, W. B. McKinnon, J. M. Moore, C. B. Olkin, J. R. Spencer, H. A. Weaver, L. A. Young, T. Andert, J. Andrews, M. Banks, B. Bauer, J. Bauman, O. S. Barnouin, P. Bedini, K. Beisser, R. A. Beyer, S. Bhaskaran, R. P. Binzel, E. Birath, M. Bird, D. J. Bogan , et al. (126 additional authors not shown)

Abstract: The Pluto system was recently explored by NASA's New Horizons spacecraft, making closest approach on 14 July 2015. Pluto's surface displays diverse landforms, terrain ages, albedos, colors, and composition gradients. Evidence is found for a water-ice crust, geologically young surface units, surface ice convection, wind streaks, volatile transport, and glacial flow. Pluto's atmosphere is highly ext… ▽ More The Pluto system was recently explored by NASA's New Horizons spacecraft, making closest approach on 14 July 2015. Pluto's surface displays diverse landforms, terrain ages, albedos, colors, and composition gradients. Evidence is found for a water-ice crust, geologically young surface units, surface ice convection, wind streaks, volatile transport, and glacial flow. Pluto's atmosphere is highly extended, with trace hydrocarbons, a global haze layer, and a surface pressure near 10 microbars. Pluto's diverse surface geology and long-term activity raise fundamental questions about how small planets remain active many billions of years after formation. Pluto's large moon Charon displays tectonics and evidence for a heterogeneous crustal composition, its north pole displays puzzling dark terrain. Small satellites Hydra and Nix have higher albedos than expected. △ Less

Submitted 26 October, 2015; originally announced October 2015.

Comments: 8 pages - Initial Science paper from NASA's New Horizons Pluto Encounter

Journal ref: Science, vol 350, 292, 2015

The Shape of Saturn's Huygens Ringlet Viewed by Cassini ISS

Authors: Joseph N. Spitale, Joseph M. Hahn

Abstract: A new model for the shape of the prominent eccentric ringlet in the gap exterior to Saturn's B-ring is developed based on Cassini imaging observations taken over about 8 years. Unlike previous treatments, the new model treats each edge of the ringlet separately. The Keplerian component of the model is consistent with results derived from Voyager observations, and $m=2$ modes forced by the nearby M… ▽ More A new model for the shape of the prominent eccentric ringlet in the gap exterior to Saturn's B-ring is developed based on Cassini imaging observations taken over about 8 years. Unlike previous treatments, the new model treats each edge of the ringlet separately. The Keplerian component of the model is consistent with results derived from Voyager observations, and $m=2$ modes forced by the nearby Mimas 2:1 Lindblad resonance are seen. Additionally, a free $m=2$ mode is seen on the outer edge of the ringlet. Significant irregular structure that cannot be described using normal-mode analysis is seen on the ringlet edges as well. Particularly on the inner edge, that structure remains coherent over multi-year intervals, moving at the local Keplerian rate. We interpret the irregular structure as the signature of embedded massive bodies. The long coherence time suggests the responsible bodies are concentrated near the edge of the ringlet. Long wake-like structures originate from two locations on the inner edge of the ringlet, revealing the locations of the two most massive embedded bodies in that region. As with the Voyager observations, the Cassini data sets showed no correlation between the width and the radius of the ringlet as would be expected for a self-gravitating configuration, except for a brief interval during late 2006, when the width-radius relation was similar to those seen in most other narrow eccentric ringlets in the Solar System. △ Less

Submitted 7 October, 2015; v1 submitted 23 September, 2015; originally announced September 2015.

Comments: 36 pages, 13 figures

Storage Ring Cross Section Measurements for Electron Impact Ionization of Fe 7+

Authors: M. Hahn, A. Becker, D. Bernhardt, M. Grieser, C. Krantz, M. Lestinsky, A. Müller, O. Novotný, R. Repnow, S. Schippers, K. Spruck, A. Wolf, D. W. Savin

Abstract: We have measured electron impact ionization (EII) for Fe 7+ from the ionization threshold up to 1200 eV. The measurements were performed using the TSR heavy ion storage ring. The ions were stored long enough prior to measurement to remove most metastables, resulting in a beam of 94% ground state ions. Comparing with the previously recommended atomic data, we find that the Arnaud & Raymond (1992) c… ▽ More We have measured electron impact ionization (EII) for Fe 7+ from the ionization threshold up to 1200 eV. The measurements were performed using the TSR heavy ion storage ring. The ions were stored long enough prior to measurement to remove most metastables, resulting in a beam of 94% ground state ions. Comparing with the previously recommended atomic data, we find that the Arnaud & Raymond (1992) cross section is up to about 40\% larger than our measurement, with the largest discrepancies below about 400~eV. The cross section of Dere (2007) agrees to within 10%, which is about the magnitude of the experimental uncertainties. The remaining discrepancies between measurement and the most recent theory are likely due to shortcomings in the theoretical treatment of the excitation-autoionization contribution. △ Less

Submitted 15 July, 2015; originally announced July 2015.

Comments: Submitted to the Astrophysical Journal

A Simple Method for Modeling Collision Processes in Plasmas with a Kappa Energy Distribution

Authors: Michael Hahn, Daniel Wolf Savin

Abstract: We demonstrate that a nonthermal distribution of particles described by a kappa distribution can be accurately approximated by a weighted sum of Maxwell-Boltzmann distributions. We apply this method to modeling collision processes in kappa-distribution plasmas, with a particular focus on atomic processes important for solar physics. The relevant collision process rate coefficients are generated by… ▽ More We demonstrate that a nonthermal distribution of particles described by a kappa distribution can be accurately approximated by a weighted sum of Maxwell-Boltzmann distributions. We apply this method to modeling collision processes in kappa-distribution plasmas, with a particular focus on atomic processes important for solar physics. The relevant collision process rate coefficients are generated by summing appropriately weighted Maxwellian rate coefficients. This method reproduces the rate coefficients for a kappa distribution to an estimated accuracy of better than 5%. This is equal to or better than the accuracy of rate coefficients generated using "reverse engineering" methods, which attempt to extract the needed cross sections from the published Maxwellian rate coefficient data and then reconvolve the extracted cross sections with the desired kappa distribution. Our approach of summing Maxwellian rate coefficients is easy to implement using existing spectral analysis software. Moreover, the weights in the sum of the Maxwell-Boltzmann distribution rate coefficients can be found for any value of the parameter kappa, thereby enabling one to model plasmas with a time-varying kappa. Tabulated Maxwellian fitting parameters are given for specific values of kappa from 1.7 to 100. We also provide polynomial fits to these parameters over this entire range. Several applications of our technique are presented, including the plasma equilibrium charge state distribution (CSD), predicting line ratios, modeling the influence of electron impact multiple ionization on the equilibrium CSD of kappa-distribution plasmas, and calculating the time-varying CSD of plasmas during a solar flare. △ Less

Submitted 21 July, 2015; v1 submitted 23 June, 2015; originally announced June 2015.

Comments: Accepted in the Astrophysical Journal

Relative Abundance Measurements in Plumes and Interplumes

Authors: Chloé Guennou, Michael Hahn, Daniel Wolf Savin

Abstract: We present measurements of relative elemental abundances in plumes and interplumes. Plumes are bright, narrow structures in coronal holes that extend along open magnetic field lines far out into the corona. Previous work has found that in some coronal structures the abundances of elements with a low first ionization potential (FIP) < 10 eV are enhanced relative to their photospheric abundances. Th… ▽ More We present measurements of relative elemental abundances in plumes and interplumes. Plumes are bright, narrow structures in coronal holes that extend along open magnetic field lines far out into the corona. Previous work has found that in some coronal structures the abundances of elements with a low first ionization potential (FIP) < 10 eV are enhanced relative to their photospheric abundances. This coronal-to-photospheric abundance ratio, commonly called the FIP bias, is typically 1 for element with a high-FIP (> 10 eV). We have used EIS spectroscopic observations made on 2007 March 13 and 14 over an ~24 hour period to characterize abundance variations in plumes and interplumes. To assess their elemental composition, we have used a differential emission measure (DEM) analysis, which accounts for the thermal structure of the observed plasma. We have used lines from ions of iron, silicon, and sulfur. From these we have estimated the ratio of the iron and silicon FIP bias relative to that for sulfur. From the results, we have created FIP-bias-ratio maps. We find that the FIP-bias ratio is sometimes higher in plumes than in interplumes and that this enhancement can be time dependent. These results may help to identify whether plumes or interplumes contribute to the fast solar wind observed in situ and may also provides constraints on the formation and heating mechanisms of plumes. △ Less

Submitted 20 March, 2015; originally announced March 2015.

Comments: 21 pages; 3 tables; 12 figures

Influence of Electron-Impact Multiple Ionization on Equilibrium and Dynamic Charge State Distributions: A Case Study Using Iron

Authors: Michael Hahn, Daniel Wolf Savin

Abstract: We describe the influence of electron-impact multiple ionization (EIMI) on the ionization balance of collisionally ionized plasmas. We are unaware of any previous ionization balance calculations that have included EIMI, which is usually assumed to be unimportant. Here, we incorporate EIMI cross-section data into calculations of both equilibrium and non-equilibrium charge-state distributions (CSDs)… ▽ More We describe the influence of electron-impact multiple ionization (EIMI) on the ionization balance of collisionally ionized plasmas. We are unaware of any previous ionization balance calculations that have included EIMI, which is usually assumed to be unimportant. Here, we incorporate EIMI cross-section data into calculations of both equilibrium and non-equilibrium charge-state distributions (CSDs). For equilibrium CSDs, we find that EIMI has only a small effect and can usually be ignored. However, for non-equilibrium plasmas the influence of EIMI can be important. In particular, we find that for plasmas in which the temperature oscillates there are significant differences in the CSD when including versus neglecting EIMI. These results have implications for modeling and spectroscopy of impulsively heated plasmas, such as nanoflare heating of the solar corona. △ Less

Submitted 17 October, 2014; originally announced October 2014.

Comments: Submitted to the Astrophysical Journal

Evidence for Wave Heating of the Quiet Sun Corona

Authors: Michael Hahn, Daniel Wolf Savin

Abstract: We have measured the energy and dissipation of Alfvenic waves in the quiet Sun. A magnetic field was used to infer the location and orientation of the magnetic field lines along which the waves are expected to travel. The waves were measured using spectral lines to infer the wave amplitude. The waves cause a non-thermal broadening of the spectral lines, which can be expressed as a non-thermal velo… ▽ More We have measured the energy and dissipation of Alfvenic waves in the quiet Sun. A magnetic field was used to infer the location and orientation of the magnetic field lines along which the waves are expected to travel. The waves were measured using spectral lines to infer the wave amplitude. The waves cause a non-thermal broadening of the spectral lines, which can be expressed as a non-thermal velocity v_nt. By combining the spectroscopic measurements with this magnetic field model we were able to trace the variation of v_nt along the magnetic field. At the footpoints of the quiet Sun loops we find that waves inject an energy flux in the range of 1.2-5.2 x 10^5 erg cm^-2 s^-1. At the minimum of this range, this amounts to more than 80% of the energy needed to heat the quiet Sun. We also find that these waves are dissipated over a region centered on the top of the loops. The position along the loop where the damping begins is strongly correlated with the length of the loop, implying that the damping mechanism depends on the global loop properties rather than on local collisional dissipation. △ Less

Submitted 11 July, 2014; originally announced July 2014.

Comments: Submitted to the Astrophysical Journal

An N-body Integrator for Gravitating Planetary Rings, and the Outer Edge of Saturn's B Ring

Authors: Joseph M. Hahn, Joseph N. Spitale

Abstract: A new symplectic N-body integrator is introduced, one designed to calculate the global 360 degree evolution of a self-gravitating planetary ring that is in orbit about an oblate planet. This freely-available code is called epi_int, and it is distinct from other such codes in its use of streamlines to calculate the effects of ring self-gravity. The great advantage of this approach is that the pertu… ▽ More A new symplectic N-body integrator is introduced, one designed to calculate the global 360 degree evolution of a self-gravitating planetary ring that is in orbit about an oblate planet. This freely-available code is called epi_int, and it is distinct from other such codes in its use of streamlines to calculate the effects of ring self-gravity. The great advantage of this approach is that the perturbing forces arise from smooth wires of ring matter rather than discreet particles, so there is very little gravitational scattering and so only a modest number of particles are needed to simulate, say, the scalloped edge of a resonantly confined ring or the propagation of spiral density waves. The code is applied to the outer edge of Saturn's B ring, and a comparison of Cassini measurements of the ring's forced response to simulations of Mimas' resonant perturbations reveals that the B ring's surface density at its outer edge is 195+-60 gm/cm^2 which, if the same everywhere across the ring would mean that the B ring's mass is about 90% of Mimas' mass. Cassini observations show that the B ring-edge has several free normal modes, which are long-lived disturbances of the ring-edge that are not driven by any known satellite resonances. Although the mechanism that excites or sustains these normal modes is unknown, we can plant such a disturbance at a simulated ring's edge, and find that these modes persist without any damping for more than ~10^5 orbits or ~100 yrs despite the simulated ring's viscosity of 100 cm^2/sec. These simulations also indicate that impulsive disturbances at a ring can excite long-lived normal modes, which suggests that an impact in the recent past by perhaps a cloud of cometary debris might have excited these disturbances which are quite common to many of Saturn's sharp-edged rings. △ Less

Submitted 5 June, 2013; v1 submitted 5 June, 2013; originally announced June 2013.

Comments: 55 pages, 13 figures, accepted for publication in the Astrophysical Journal

Observational Quantification of the Energy Dissipated by Alfvén Waves in a Polar Coronal Hole: Evidence that Waves Drive the Fast Solar Wind

Authors: Michael Hahn, Daniel W. Savin

Abstract: We present a measurement of the energy carried and dissipated by Alfvén waves in a polar coronal hole. Alfvén waves have been proposed as the energy source that heats the corona and drives the solar wind. Previous work has shown that line widths decrease with height in coronal holes, which is a signature of wave damping, but have been unable to quantify the energy lost by the waves. This is becaus… ▽ More We present a measurement of the energy carried and dissipated by Alfvén waves in a polar coronal hole. Alfvén waves have been proposed as the energy source that heats the corona and drives the solar wind. Previous work has shown that line widths decrease with height in coronal holes, which is a signature of wave damping, but have been unable to quantify the energy lost by the waves. This is because line widths depend on both the non-thermal velocity v_nt and the ion temperature T_i. We have implemented a means to separate the T_i and v_nt contributions using the observation that at low heights the waves are undamped and the ion temperatures do not change with height. This enables us to determine the amount of energy carried by the waves at low heights, which is proportional to v_nt. We find the initial energy flux density present was 6.7 +/- 0.7 x 10^5 erg cm^-2 s^-1, which is sufficient to heat the coronal hole and acccelerate the solar wind during the 2007 - 2009 solar minimum. Additionally, we find that about 85% of this energy is dissipated below 1.5 R_sun, sufficiently low that thermal conduction can transport the energy throughout the coronal hole, heating it and driving the fast solar wind. The remaining energy is roughly consistent with what models show is needed to provide the extended heating above the sonic point for the fast solar wind. We have also studied T_i, which we found to be in the range of 1 - 2 MK, depending on the ion species. △ Less

Submitted 29 August, 2013; v1 submitted 21 February, 2013; originally announced February 2013.

Comments: Accepted for the Astrophysical Journal

Measurements of Anisotropic Ion Temperatures, Non-Thermal Velocities, and Doppler Shifts in a Coronal Hole

Authors: M. Hahn, D. W. Savin

Abstract: We present a new diagnostic allowing one to measure the anisotropy of ion temperatures and non-thermal velocities as well as Doppler shifts with respect to the ambient magnetic field. This method provides new results, as well as independent test for previous measurements obtained with other techniques. Our spectral data come from observations of a low latitude, on-disk coronal hole. A potential fi… ▽ More We present a new diagnostic allowing one to measure the anisotropy of ion temperatures and non-thermal velocities as well as Doppler shifts with respect to the ambient magnetic field. This method provides new results, as well as independent test for previous measurements obtained with other techniques. Our spectral data come from observations of a low latitude, on-disk coronal hole. A potential field source surface model was used to calculate the angle between the magnetic field lines and the line of sight for each spatial bin of the observation. A fit was performed to determine the line widths and Doppler shifts parallel and perpendicular to the magnetic field. For each line width component we derived parallel and perpendicular ion temperatures and non-thermal velocities. The perpendicular ion temperature was cooler than off-limb polar coronal hole measurements. The parallel ion temperature was consistent with a uniform temperature of 1.8 +/- 0.2 x 10^{6} K for each ion. Since parallel ion heating is expected to be weak, this ion temperature should reflect the proton temperature. A comparison between our results and others implies a large proton temperature gradient around 1.02 R_sun. The non-thermal velocities are thought to be proportional to the amplitudes of various waves. Our results for the perpendicular non-thermal velocity agree with Alfvén wave amplitudes inferred from off-limb polar coronal hole line width measurements. Our parallel non-thermal velocity results are consistent with slow magnetosonic wave amplitudes inferred from Fourier analysis of time varying intensity fluctuations. Doppler shift measurements yield outflows of ~5 km s^-1 for ions formed over a broad temperature range. This differs from other studies which found a strong Doppler shift dependence on formation temperature. △ Less

Submitted 14 December, 2012; v1 submitted 5 November, 2012; originally announced November 2012.

Comments: Accepted for The Astrophysical Journal

Electron-ion Recombination of Fe XII forming Fe XI: Laboratory Measurements and Theoretical Calculations

Authors: O. Novotný, N. R. Badnell, D. Bernhardt, M. Grieser, M. Hahn, C. Krantz, M. Lestinsky, A. Müller, R. Repnow, S. Schippers, A. Wolf, D. W. Savin

Abstract: We have measured electron-ion recombination for Fe XII forming Fe XI using a merged beams configuration at the heavy-ion storage ring TSR located at the Max Planck Institute for Nuclear Physics in Heidelberg, Germany. The measured merged beams recombination rate coefficient (MBRRC) for collision energies from 0 to 1500 eV is presented. This work uses a new method for determining the absolute MBRRC… ▽ More We have measured electron-ion recombination for Fe XII forming Fe XI using a merged beams configuration at the heavy-ion storage ring TSR located at the Max Planck Institute for Nuclear Physics in Heidelberg, Germany. The measured merged beams recombination rate coefficient (MBRRC) for collision energies from 0 to 1500 eV is presented. This work uses a new method for determining the absolute MBRRC based on a comparison of the ion beam decay rate with and without the electron beam on. For energies below 75 eV, the spectrum is dominated by dielectronic recombination (DR) resonances associated with 3s-3p and 3p-3d core excitations. At higher energies we observe contributions from 3-N' and 2-N' core excitations DR. We compare our experimental results to state-of-the-art multi-configuration Breit-Pauli (MCBP) calculations and find significant differences, both in resonance energies and strengths. We have extracted the DR contributions from the measured MBRRC data and transformed them into a plasma recombination rate coefficient (PRRC) for temperatures in the range of 10^3 to 10^7 K. We show that the previously recommended DR data for Fe XII significantly underestimate the PRRC at temperatures relevant for both photoionized plasmas (PPs) and collisionaly ionized plasmas (CPs). This is to be contrasted with our MCBP PRRC results which agree with the experiment to within 30% at PP temperatures and even better at CP temperatures. We find this agreement despite the disagreement shown by the detailed comparison between our MCBP and experimental MBRRC results. Lastly, we present a simple parameterized form of the experimentally derived PRRC for easy use in astrophysical modelling codes. △ Less

Submitted 27 April, 2012; originally announced April 2012.

Comments: April 24, 2012: accepted for publication in The Astrophysical Journal

Journal ref: ApJ 753 (2012) 57

Evidence of Wave Damping at Low Heights in a Polar Coronal Hole

Authors: Michael Hahn, Enrico Landi, Daniel Wolf Savin

Abstract: We have measured the widths of spectral lines from a polar coronal hole using the Extreme Ultraviolet Imaging Spectrometer onboard Hinode. Polar coronal holes are regions of open magnetic field and the source of the fast solar wind. We find that the line widths decrease at relatively low heights. Previous observations have attributed such decreases to systematic effects, but we find that such effe… ▽ More We have measured the widths of spectral lines from a polar coronal hole using the Extreme Ultraviolet Imaging Spectrometer onboard Hinode. Polar coronal holes are regions of open magnetic field and the source of the fast solar wind. We find that the line widths decrease at relatively low heights. Previous observations have attributed such decreases to systematic effects, but we find that such effects are too small to explain our results. We conclude that the line narrowing is real. The non-thermal line widths are believed to be proportional to the amplitude of Alfven waves propagating along these open field lines. Our results suggest that Alfven waves are damped at unexpectedly low heights in a polar coronal hole. We derive an estimate on the upper limit for the energy dissipated between 1.1 and 1.3 solar radii and find that it is enough to account for up to 70% of that required to heat the polar coronal hole and accelerate the solar wind. △ Less

Submitted 30 April, 2012; v1 submitted 8 February, 2012; originally announced February 2012.

Comments: Accepted for publication in the Astrophysical Journal, April 2012

Storage Ring Cross Section Measurements for Electron Impact Ionization of Fe^11+ Forming Fe^12+ and Fe^13+

Authors: M. Hahn, D. Bernhardt, M. Grieser, C. Krantz, M. Lestinsky, A. Mueller, O. Novotny, R. Repnow, S. Schippers, A. Wolf, D. W. Savin

Abstract: We report ionization cross section measurements for electron impact single ionization (EISI) of Fe^11+$ forming Fe^12+ and electron impact double ionization (EIDI) of Fe^11+ forming Fe^13+. The measurements cover the center-of-mass energy range from approximately 230 eV to 2300 eV. The experiment was performed using the heavy ion storage ring TSR located at the Max-Planck-Institut fur Kernphysik i… ▽ More We report ionization cross section measurements for electron impact single ionization (EISI) of Fe^11+$ forming Fe^12+ and electron impact double ionization (EIDI) of Fe^11+ forming Fe^13+. The measurements cover the center-of-mass energy range from approximately 230 eV to 2300 eV. The experiment was performed using the heavy ion storage ring TSR located at the Max-Planck-Institut fur Kernphysik in Heidelberg, Germany. The storage ring approach allows nearly all metastable levels to relax to the ground state before data collection begins. We find that the cross section for single ionization is 30% smaller than was previously measured in a single pass experiment using an ion beam with an unknown metastable fraction. We also find some significant differences between our experimental cross section for single ionization and recent distorted wave (DW) calculations. The DW Maxwellian EISI rate coefficient for Fe^11+ forming Fe^12+ may be underestimated by as much as 25% at temperatures for which Fe^11+ is abundant in collisional ionization equilibrium. This is likely due to the absence of 3s excitation-autoionization (EA) in the calculations. However, a precise measurement of the cross section due to this EA channel was not possible because this process is not distinguishable experimentally from electron impact excitation of an n=3 electron to levels of n > 44 followed by field ionization in the charge state analyzer after the interaction region. Our experimental results also indicate that the double ionization cross section is dominated by the indirect process in which direct single ionization of an inner shell 2l electron is followed by autoionization resulting in a net double ionization. △ Less

Submitted 11 January, 2011; originally announced January 2011.

Comments: Accepted for publication in Astrophysical Journal

Journal ref: Astrophys. J. 729 (2011) 76

Properties of a Polar Coronal Hole During the Solar Minimum in 2007

Authors: M. Hahn, P. Bryans, E. Landi, M. P. Miralles, D. W. Savin

Abstract: We report measurements of a polar coronal hole during the recent solar minimum using the Extreme Ultraviolet Imaging Spectrometer on Hinode. Five observations are analyzed that span the polar coronal hole from the central meridian to the boundary with the quiet Sun corona. We study the observations above the solar limb in the height range of 1.03 - 1.20 solar radii. The electron temperature Te and… ▽ More We report measurements of a polar coronal hole during the recent solar minimum using the Extreme Ultraviolet Imaging Spectrometer on Hinode. Five observations are analyzed that span the polar coronal hole from the central meridian to the boundary with the quiet Sun corona. We study the observations above the solar limb in the height range of 1.03 - 1.20 solar radii. The electron temperature Te and emission measure EM are found using the Geometric mean Emission Measure (GEM) method. The EM derived from the elements Fe, Si, S, and Al are compared in order to measure relative coronal-to-photospheric abundance enhancement factors. We also studied the ion temperature Ti and the non-thermal velocity Vnt using the line profiles. All these measurements are compared to polar coronal hole observations from the previous (1996-1997) solar minimum and to model predictions for relative abundances. There are many similarities in the physical properties of the polar coronal holes between the two minima at these low heights. We find that Te, electron density Ne, and Ti are comparable in both minima. Te shows a comparable gradient with height. Both minima show a decreasing Ti with increasing charge-to-mass ratio q/M. A previously observed upturn of Ti for ions above q/M > 0.25 was not found here. We also compared relative coronal-to-photospheric elemental abundance enhancement factors for a number of elements. These ratios were about 1 for both the low first ionization potential (FIP) elements Si and Al and the marginally high FIP element S relative to the low FIP element Fe, as is expected based on earlier observations and models for a polar coronal hole. These results are consistent with no FIP effect in a polar coronal hole. △ Less

Submitted 20 October, 2010; originally announced October 2010.

Comments: Accepted for publication in the Astrophysical Journal

Diagnosing Circumstellar Debris Disks

Authors: Joseph M. Hahn

Abstract: (Abridged) A numerical model of a circumstellar debris disk is developed and applied to observations of the circumstellar dust orbiting beta Pictoris. The model accounts for the rates at which dust is produced by collisions among unseen planetesimals, and the rate at which dust grains are destroyed due to collisions. The model also accounts for the effects of radiation pressure, which is the domin… ▽ More (Abridged) A numerical model of a circumstellar debris disk is developed and applied to observations of the circumstellar dust orbiting beta Pictoris. The model accounts for the rates at which dust is produced by collisions among unseen planetesimals, and the rate at which dust grains are destroyed due to collisions. The model also accounts for the effects of radiation pressure, which is the dominant perturbation on the disk's smaller but abundant dust grains. Solving the resulting system of rate equations then provides the dust abundances versus grain size and over time. Those solutions also provide the dust grains' collisional lifetime versus grain size, and the debris disk's optical depth and surface brightness versus distance from the star. Comparison to observations then yields estimates of the unseen planetesimal disk's radius, and the rate at which the disk sheds mass due to planetesimal grinding. The model is then applied to optical observations of the edge-on dust disk orbiting beta Pictoris, and good agreement is achieved when the unseen planetesimal disk is broad, with 75<r<150 AU. If it is assumed that the dust grains are bright like Saturn's icy rings, then the cross section of dust in the disk is A_d~2x10^20 km^2 and its mass is M_d~11 lunar masses. In this case the planetesimal disk's dust production rate is quite heavy, dM_d/dt~9 earth-masses per Myr, implying that there is or was a substantial amount of planetesimal mass there, at least 110 earth-masses. But if the dust grains are darker than assumed, then the planetesimal disk's mass-loss rate and its total mass are heavier. In fact, the apparent dearth of any major planets in this region, plus the planetesimal disk's heavy mass-loss rate, suggests that the 75<r<150 AU zone at beta Pic might be a region of planetesimal destruction, rather than a site of ongoing planet formation. △ Less

Submitted 22 June, 2010; originally announced June 2010.

Comments: Accepted for publication in the Astrophysical Journal

Dynamics of the Sharp Edges of Broad Planetary Rings

Authors: Joseph M. Hahn, Joseph N. Spitale, Carolyn C. Porco

Abstract: (Abridged) The following describes a model of a broad planetary ring whose sharp edge is confined by a satellite's m^th Lindblad resonance (LR). This model uses a streamline formalism to calculate the ring's internal forces, namely, ring gravity, pressure, viscosity, as well as a hypothetical drag force. The model calculates the streamlines' forced orbit elements and surface density throughout t… ▽ More (Abridged) The following describes a model of a broad planetary ring whose sharp edge is confined by a satellite's m^th Lindblad resonance (LR). This model uses a streamline formalism to calculate the ring's internal forces, namely, ring gravity, pressure, viscosity, as well as a hypothetical drag force. The model calculates the streamlines' forced orbit elements and surface density throughout the perturbed ring. The model is then applied to the outer edge of Saturn's B ring, which is maintained by an m=2 inner LR with the satellite Mimas. Ring models are used to illustrate how a ring's perturbed state depends on the ring's physical properties: surface density, viscosity, dispersion velocity, and the hypothetical drag force. A comparison of models to the observed outer B ring suggests that the ring's surface density there is between 10 and 280 gm/cm^2. The ring's edge also indicates where the viscous torque counterbalances the perturbing satellite's gravitational torque on the ring. But an examination of seemingly conventional viscous B ring models shows that they all fail to balance these torques at the ring's edge. This is due ring self-gravity and the fact that a viscous ring tends to be nearly peri-aligned with the satellite, which reduces the satellite's torque on the ring and makes the ring's edge more difficult to maintain. Nonetheless, the following shows that a torque balance can still be achieved in a viscous B ring, but only in an extreme case where the ratio of the ring's bulk/shear viscosities satisfy ~10^4. However, if the dissipation of the ring's forced motions is instead dominated by a weak drag force, then the satellite can exert a much stronger torque that can counterbalance the ring's viscous torque. △ Less

Submitted 4 April, 2009; originally announced April 2009.

Comments: Accepted for publication in the Astrophysical Journal on April 3, 2009

Journal ref: Astrophys.J.699:686-710,2009

The Secular Evolution of a Close Ring-Satellite System: The Excitation of Spiral Density Waves at a Nearby Gap Edge

Authors: Joseph M. Hahn

Abstract: The Lagrange planetary equations are used to study to secular evolution of a small, eccentric satellite that orbits within a narrow gap in a broad, self-gravitating planetary ring. These equations show that the satellite's secular perturbations of the ring will excite a very long-wavelength spiral density wave that propagates away from the gap's outer edge. The amplitude of these waves, as well… ▽ More The Lagrange planetary equations are used to study to secular evolution of a small, eccentric satellite that orbits within a narrow gap in a broad, self-gravitating planetary ring. These equations show that the satellite's secular perturbations of the ring will excite a very long-wavelength spiral density wave that propagates away from the gap's outer edge. The amplitude of these waves, as well as their dispersion relation, are derived here. That dispersion relation reveals that a planetary ring can sustain two types of density waves: long waves that, in Saturn's A ring, would have wavelengths of order 100 km, and short waves that tend to be very nonlinear and are expected to quickly damp. The excitation of these waves also transports angular momentum from the ring to the satellite in a way that damps the satellite's eccentricity e, which also tends to reduce the amplitude of subsequent waves. The rate of eccentricity damping due to this wave action is then compared to the rates at which the satellite's Lindblad and corotation resonances alter the satellite's e. These results are then applied to the gap-embedded Saturnian satellites Pan and Daphnis, and the long-term stability of their eccentricities is assessed. △ Less

Submitted 4 March, 2008; originally announced March 2008.

Comments: Accepted for publication in the Astrophysical Journal

The Secular Evolution of a Close Ring-Satellite System: The Excitation of Spiral Bending Waves at a Nearby Gap Edge

Authors: J. M. Hahn

Abstract: The secular perturbations exerted by an inclined satellite orbiting in a gap in a broad planetary ring tends to excite the inclinations of the nearby ring particles, and the ring's self-gravity can allow that disturbance to propagate away in the form of a spiral bending wave. The amplitude of this spiral bending wave is determined, as well as the wavelength, which shrinks as the waves propagate… ▽ More The secular perturbations exerted by an inclined satellite orbiting in a gap in a broad planetary ring tends to excite the inclinations of the nearby ring particles, and the ring's self-gravity can allow that disturbance to propagate away in the form of a spiral bending wave. The amplitude of this spiral bending wave is determined, as well as the wavelength, which shrinks as the waves propagate outwards due to the effects of the central planet's oblateness. The excitation of these bending waves also damps the satellite's inclination I. This secular I damping is also compared to the inclination excitation that is due to the satellite's many other vertical resonances in the ring, and the condition for inclination damping is determined. The secular I damping is likely responsible for confining the orbits of Saturn's two known gap-embedded moons, Pan and Daphnis, to the ring plane. △ Less

Submitted 25 May, 2007; originally announced May 2007.

Comments: 31 pages, 8 figures

Journal ref: Astrophys.J.665:856-865,2007

Neptune's Migration into a Stirred-Up Kuiper Belt: A Detailed Comparison of Simulations to Observations

Authors: Joseph M. Hahn, Renu Malhotra

Abstract: Nbody simulations are used to examine the consequences of Neptune's outward migration into the Kuiper Belt, with the simulated endstates being compared rigorously and quantitatively to the observations. These simulations confirm the findings of Chiang et al. (2003), who showed that Neptune's migration into a previously stirred-up Kuiper Belt can account for the Kuiper Belt Objects (KBOs) known t… ▽ More Nbody simulations are used to examine the consequences of Neptune's outward migration into the Kuiper Belt, with the simulated endstates being compared rigorously and quantitatively to the observations. These simulations confirm the findings of Chiang et al. (2003), who showed that Neptune's migration into a previously stirred-up Kuiper Belt can account for the Kuiper Belt Objects (KBOs) known to librate at Neptune's 5:2 resonance. We also find that capture is possible at many other weak, high-order mean motion resonances, such as the 11:6, 13:7, 13:6, 9:4, 7:3, 12:5, 8:3, 3:1, 7:2, and the 4:1. The more distant of these resonances, such as the 9:4, 7:3, 5:2, and the 3:1, can also capture particles in stable, eccentric orbits beyond 50 AU, in the region of phase space conventionally known as the Scattered Disk. Indeed, 90% of the simulated particles that persist over the age of the Solar System in the so-called Scattered Disk zone never had a close encounter with Neptune, but instead were promoted into these eccentric orbits by Neptune's resonances during the migration epoch. This indicates that the observed Scattered Disk might not be so scattered. This model also produced only a handful of Centaurs, all of which originated at Neptune's mean motion resonances in the Kuiper Belt. We also report estimates of the abundances and masses of the Belt's various subpopulations (e.g., the resonant KBOs, the Main Belt, and the so-called Scattered Disk), and also provide upper limits on the abundance of Centaurs and Neptune's Trojans, as well as upper limits on the sizes and abundances of hypothetical KBOs that might inhabit the a>50 AU zone. △ Less

Submitted 13 July, 2005; originally announced July 2005.

Comments: 60 pages, 16 figures. Accepted for publication in the Astronomical Journal

Journal ref: Astron.J.130:2392-2414,2005

The Secular Evolution of the Primordial Kuiper Belt

Authors: Joseph M. Hahn

Abstract: A model that computes the secular evolution of a gravitating disk-planet system is developed. The disk is treated as a set of gravitating rings, with the rings'/planets' time-evolution governed by the classical Laplace-Lagrange solution for secular evolution but modified to account for the disk's finite thickness h. This system's Lagrange planetary equations yield a particular class of spiral wa… ▽ More A model that computes the secular evolution of a gravitating disk-planet system is developed. The disk is treated as a set of gravitating rings, with the rings'/planets' time-evolution governed by the classical Laplace-Lagrange solution for secular evolution but modified to account for the disk's finite thickness h. This system's Lagrange planetary equations yield a particular class of spiral wave solutions, usually denoted as apsidal density waves and nodal bending waves. There are two varieties of apsidal waves:long waves and short waves. Planets typically launch long density waves at the disk's nearer edge or else at a secular resonance in the disk, and these waves ultimately reflect downstream at a more distant disk edge or else at a Q-barrier in the disk, whereupon they return as short density waves. Planets also launch nodal bending waves, and these have the property that they can stall in the disk, that is, their group velocity plummets to zero upon approaching a disk region too thick to support the continued propagation of bending waves. The rings model is used to compute the secular evolution of a Kuiper Belt having a variety of masses, and it is shown that the early massive Belt was very susceptible to the propagation of low-amplitude apsidal and nodal waves launched by the giant planets. △ Less

Submitted 30 May, 2003; originally announced May 2003.

Comments: 45 pages, 6 figures

Report number: LPI Contribution #1165

Journal ref: Astrophys.J. 595 (2003) 531-549

Clementine Observations of the Zodiacal Light and the Dust Content of the Inner Solar System

Authors: J. M. Hahn, H. A. Zook, B. Cooper, B. Sunkara

Abstract: Using the Moon to occult the Sun, the Clementine spacecraft used its navigation cameras to map the inner zodiacal light at optical wavelengths over elongations of 3-30 degrees from the Sun. This surface brightness map is then used to infer the spatial distribution of interplanetary dust over heliocentric distances of about 10 solar radii to the orbit of Venus. We also apply a simple model that a… ▽ More Using the Moon to occult the Sun, the Clementine spacecraft used its navigation cameras to map the inner zodiacal light at optical wavelengths over elongations of 3-30 degrees from the Sun. This surface brightness map is then used to infer the spatial distribution of interplanetary dust over heliocentric distances of about 10 solar radii to the orbit of Venus. We also apply a simple model that attributes the zodiacal light as being due to three dust populations having distinct inclination distributions, namely, dust from asteroids and Jupiter-family comets (JFCs), dust from Halley-type comets, and an isotropic cloud of dust from Oort Cloud comets. The best-fitting scenario indicates that asteroids + JFCs are the source of about 45% of the optical dust cross-section seen in the ecliptic at 1 AU, but that at least 89% of the dust cross-section enclosed by a 1 AU radius sphere is of a cometary origin. When these results are extrapolated out to the asteroid belt, we find an upper limit on the mass of the light-reflecting asteroidal dust that is equivalent to a 12 km asteroid, and a similar extrapolation of the isotropic dust cloud out to Oort Cloud distances yields a mass equivalent to a 30 km comet, although the latter mass is uncertain by orders of magnitude. △ Less

Submitted 7 April, 2002; originally announced April 2002.

Comments: To be published in Icarus

Orbital Evolution of Planets Embedded in a Planetesimal Disk

Authors: Joseph M. Hahn, Renu Malhotra

Abstract: The existence of the Oort Comet Cloud, the Kuiper Belt, and plausible inefficiencies in planetary core formation, all suggest that there was once a residual planetesimal disk of mass 10-100 Earth-masses in the vicinity of the giant planets following their formation. Since removal of this disk requires an exchange of orbital energy and angular momentum with the planets, significant planetary migr… ▽ More The existence of the Oort Comet Cloud, the Kuiper Belt, and plausible inefficiencies in planetary core formation, all suggest that there was once a residual planetesimal disk of mass 10-100 Earth-masses in the vicinity of the giant planets following their formation. Since removal of this disk requires an exchange of orbital energy and angular momentum with the planets, significant planetary migration can ensue. The planet migration phenomenon is examined numerically by evolving the orbits of the giant planets while they are embedded in a planetesimal disk having a mass of M_d=10 to 200 Earth-masses. We find that Saturn, Uranus, and Neptune evolve radially outwards as they scatter the planetesimals, while Jupiter's orbit shrinks as it ejects mass. Higher-mass disks result in more rapid and extensive planet migration. If orbit expansion and resonance trapping by Neptune is invoked to explain the eccentricities of Pluto and its cohort of Kuiper Belt Objects at Neptune's 3:2 mean-motion resonance, then our simulations suggest that a disk mass of order M_d~50 Earth-masses is required to expand Neptune's orbit by ~7 AU in order to pump up Plutino eccentricities to e~0.3. Such planet migration implies that the initial Solar System was more compact in the past, with the Jupiter-Neptune separation having been smaller by about 30%. The planetesimal disk is also the source of the Oort Cloud of comets. Using the results of our simulations together with a simple treatment of Oort Cloud dynamics, we estimate that ~12 Earth-masses of disk material was initially deposited in the Oort Cloud, of which ~4 Earth-masses will persist over the age of the Solar System. The majority of these comets originated from the Saturn-Neptune region of the solar nebula. △ Less

Submitted 25 February, 1999; originally announced February 1999.

Comments: to be published in Astronomical Journal

Report number: 966