Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “STELLAR RADIATION”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

On large-scale dynamos with stable stratification and the application to stellar radiative zones

Our understanding of large-scale magnetic fields in stellar radiative zones remains fragmented and incomplete. Such magnetic fields, which must be produced by some form of dynamo mechanism, are thought to dominate angular-momentum transport, making them crucial to stellar evolution. A major difficulty is the effect of stable stratification, which generally suppresses dynamo action. We explore the effects of stable stratification on mean-field dynamo theory with a particular focus on a non-helical large-scale dynamo (LSD) mechanism known as the magnetic shear-current effect. We find that the mechanism is robust to increasing stable stratification as long as the original requirements for its operation are met: a source of shear and non-helical magnetic fluctuations (e.g. from a small-scale dynamo). Both are plausibly sourced in the presence of differential rotation. Our idealized direct numerical simulations, supported by mean-field theory, demonstrate the generation of near equipartition large-scale toroidal fields. Additionally, a scan over magnetic Reynolds number shows no change in the growth or saturation of the LSD, providing good numerical evidence of a dynamo mechanism resilient to catastrophic quenching, which has been an issue for helical dynamos. These properties – the absence of catastrophic quenching and robustness to stable stratification – make the mechanism a plausible candidate for generating in situ large-scale magnetic fields in stellar radiative zones.

79 ASTRONOMY AND ASTROPHYSICS↗

OUTWARD MOTION OF POROUS DUST AGGREGATES BY STELLAR RADIATION PRESSURE IN PROTOPLANETARY DISKS

We study the dust motion at the surface layer of protoplanetary disks. Dust grains in the surface layer migrate outward owing to angular momentum transport via gas-drag force induced by the stellar radiation pressure. In this study we calculate the mass flux of the outward motion of compact grains and porous dust aggregates by the radiation pressure. The radiation pressure force for porous dust aggregates is calculated using the T-Matrix Method for the Clusters of Spheres. First, we confirm that porous dust aggregates are forced by strong radiation pressure even if they grow to be larger aggregates, in contrast to homogeneous and spherical compact grains, for which radiation pressure efficiency becomes lower when their sizes increase. In addition, we find that the outward mass flux of porous dust aggregates with monomer size of 0.1 μm is larger than that of compact grains by an order of magnitude at the disk radius of 1 AU, when their sizes are several microns. This implies that large compact grains like calcium-aluminum-rich inclusions are hardly transported to the outer region by stellar radiation pressure, whereas porous dust aggregates like chondritic-porous interplanetary dust particles are efficiently transported to the comet formation region. Crystalline silicates are possibly transported in porous dust aggregates by stellar radiation pressure from the inner hot region to the outer cold cometary region in the protosolar nebula.

79 ASTRONOMY AND ASTROPHYSICS↗

Baroclinic instability in stellar radiation zones

Surfaces of constant pressure and constant density do not coincide in differentially rotating stars. Stellar radiation zones with baroclinic stratification can be unstable. Instabilities in radiation zones are of crucial importance for angular momentum transport, mixing of chemical species, and, possibly, for magnetic field generation. This paper performs linear analysis of baroclinic instability in differentially rotating stars. Linear stability equations are formulated for differential rotation of arbitrary shape and then solved numerically for rotation nonuniform in radius. As the differential rotation increases, r- and g-modes of initially stable global oscillations transform smoothly into growing modes of baroclinic instability. The instability can therefore be interpreted as stability loss to r- and g-modes excitation. Regions of stellar parameters where r- or g-modes are preferentially excited are defined. Baroclinic instability onsets at a very small differential rotation of below 1%. The characteristic time of instability growth is about 1000 rotation periods. Growing disturbances possess kinetic helicity. Magnetic field generation by the turbulence resulting from baroclinic instability in differentially rotating radiation zones is therefore possible.

79 ASTRONOMY AND ASTROPHYSICS↗

OBSERVATIONAL CONSTRAINTS ON THE STELLAR RADIATION FIELD IMPINGING ON TRANSITIONAL DISK ATMOSPHERES

Mid-infrared atomic and ionic line ratios measured in spectra of pre-main-sequence stars are sensitive indicators of the hardness of the radiation field impinging on the disk surface. We present a low-resolution Spitzer IRS search for [Ar II] at 6.98 {mu}m, [Ne II] at 12.81 {mu}m, and [Ne III] 15.55 {mu}m lines in 56 transitional disks. These objects, characterized by reduced near-infrared but strong far-infrared excess emission, are ideal targets to set constraints on the stellar radiation field onto the disk, because their spectra are not contaminated by shock emission from jets/outflows or by molecular emission lines. After demonstrating that we can detect [Ne II] lines and recover their fluxes from the low-resolution spectra, here we report the first detections of [Ar II] lines toward protoplanetary disks. We did not detect [Ne III] emission in any of our sources. Our [Ne II]/[Ne III] line flux ratios combined with literature data suggest that a soft-EUV or X-ray spectrum produces these gas lines. Furthermore, the [Ar II]/[Ne II] line flux ratios point to a soft X-ray and/or soft-EUV stellar spectrum as the ionization source of the [Ar II] and [Ne II] emitting layer of the disk. If the soft X-ray component dominates over the EUV, then we would expect larger photoevaporation rates and, hence, a reduction of the time available to form planets.

79 ASTRONOMY AND ASTROPHYSICS↗

An iterative phase-space explicit discontinuous Galerkin method for stellar radiative transfer in extended atmospheres

In this work, a phase-space discontinuous Galerkin (PSDG) method is presented for the solution of stellar radiative transfer problems. It allows for greater adaptivity than competing methods without sacrificing generality. The method is extensively tested on a spherically symmetric, static, inverse-power-law scattering atmosphere. Results for different sizes of atmospheres and intensities of scattering agreed with asymptotic values. The exponentially decaying behavior of the radiative field in the diffusive-transparent transition region, and the forward peaking behavior at the surface of extended atmospheres were accurately captured. The integrodifferential equation of radiation transfer is solved iteratively by alternating between the radiative pressure equation and the original equation with the integral term treated as an energy density source term. In each iteration, the equations are solved via an explicit, flux-conserving, discontinuous Galerkin method. Finite elements are ordered in wave fronts perpendicular to the characteristic curves so that elemental linear algebraic systems are solved quickly by sweeping the phase space element by element. Two implementations of a diffusive boundary condition at the origin are demonstrated wherein the finite discontinuity in the radiation intensity is accurately captured by the proposed method. This allows for a consistent mechanism to preserve photon luminosity. The method was proved to be robust and fast, and a case is made for the adequacy of parallel processing. In addition to classical two-dimensional plots, results of normalized radiation intensity were mapped onto a log-polar surface exhibiting all distinguishing features of the problem studied.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

ROTATION AND STABILITY OF THE TOROIDAL MAGNETIC FIELD IN STELLAR RADIATION ZONES

The stability of the magnetic field in radiation zones is of crucial importance for mixing and angular momentum transport in the stellar interior. We consider the stability properties of stars containing a predominant toroidal field in spherical geometry by means of a linear stability in the Boussinesq approximation taking into account the effect of thermal conductivity. We calculate the growth rate of instability and analyze in detail the effects of stable stratification and heat transport. We argue that the stabilizing influence of gravity can never entirely suppress the instability caused by electric currents in radiation zones. However, the stable stratification can essentially decrease the growth rate of instability.

79 ASTRONOMY AND ASTROPHYSICS↗

Stochastic microhertz gravitational radiation from stellar convection

High Reynolds-number turbulence driven by stellar convection in main-sequence stars generates stochastic gravitational radiation. We calculate the wave-strain power spectral density as a function of the zero-age main-sequence mass for an individual star and for an isotropic, universal stellar population described by the Salpeter initial mass function and redshift-dependent Hopkins-Beacom star formation rate. The spectrum is a broken power law, which peaks near the turnover frequency of the largest turbulent eddies. The signal from the Sun dominates the universal background. For the Sun, the far-zone power spectral density peaks at S(f {sub peak}) = 5.2 × 10{sup –52} Hz{sup –1} at frequency f {sub peak} = 2.3 × 10{sup –7} Hz. However, at low observing frequencies f < 3 × 10{sup –4} Hz, the Earth lies inside the Sun's near zone and the signal is amplified to S {sub near}(f {sub peak}) = 4.1 × 10{sup –27} Hz{sup –1} because the wave strain scales more steeply with distance (∝d {sup –5}) in the near zone than in the far zone (∝d {sup –1}). Hence the Solar signal may prove relevant for pulsar timing arrays. Other individual sources and the universal background fall well below the projected sensitivities of the Laser Interferometer Space Antenna and next-generation pulsar timing arrays. Stellar convection sets a fundamental noise floor for more sensitive stochastic gravitational-wave experiments in the more distant future.

79 ASTRONOMY AND ASTROPHYSICS↗

Formation of primordial supermassive stars by rapid mass accretion

Supermassive stars (SMSs) forming via very rapid mass accretion ( M-dot {sub ∗}≳0.1 M{sub ⊙} yr{sup −1}) could be precursors of supermassive black holes observed beyond a redshift of about six. Extending our previous work, here we study the evolution of primordial stars growing under such rapid mass accretion until the stellar mass reaches 10{sup 4–5} M {sub ☉}. Our stellar evolution calculations show that a star becomes supermassive while passing through the 'supergiant protostar' stage, whereby the star has a very bloated envelope and a contracting inner core. The stellar radius increases monotonically with the stellar mass until ≅ 100 AU for M {sub *} ≳ 10{sup 4} M {sub ☉}, after which the star begins to slowly contract. Because of the large radius, the effective temperature is always less than 10{sup 4} K during rapid accretion. The accreting material is thus almost completely transparent to the stellar radiation. Only for M {sub *} ≳ 10{sup 5} M {sub ☉} can stellar UV feedback operate and disturb the mass accretion flow. We also examine the pulsation stability of accreting SMSs, showing that the pulsation-driven mass loss does not prevent stellar mass growth. Observational signatures of bloated SMSs should be detectable with future observational facilities such as the James Webb Space Telescope. Our results predict that an inner core of the accreting SMS should suffer from the general relativistic instability soon after the stellar mass exceeds 10{sup 5} M {sub ☉}. An extremely massive black hole should form after the collapse of the inner core.

79 ASTRONOMY AND ASTROPHYSICS↗

Reionization with galaxies and active galactic nuclei

In this work we investigate the properties of the sources that reionized the intergalactic medium (IGM) in the high-redshift Universe. Using a semi-analytical model aimed at reproducing galaxies and black holes in the first ~1.5 Gyr of the Universe, we revisit the relative role of star formation and black hole accretion in producing ionizing photons that can escape into the IGM. Both star formation and black hole accretion are regulated by supernova feedback, resulting in black hole accretion being stunted in low-mass haloes. We explore a wide range of combinations for the escape fraction of ionizing photons (redshift-dependent, constant, and scaling with stellar mass) from both star formation ( |$\langle f_{\rm esc}^{\rm sf} \rangle$| ) and AGN ( |$f_{\rm esc}^{\rm bh}$| ) to find: (i) the ionizing budget is dominated by stellar radiation from low stellar mass ( |$M_*\lt 10^9 \, {\rm \rm M_\odot }$| ) galaxies at z > 6 with the AGN contribution (driven by |$M_{bh}\gt 10^6 \, {\rm \rm M_\odot }$| black holes in |$M_* \gtrsim 10^9\, {\rm \rm M_\odot}$| galaxies) dominating at lower redshifts; (ii) AGN only contribute |$10-25{{\ \rm per\ cent}}$| to the cumulative ionizing emissivity by z = 4 for the models that match the observed reionization constraints; (iii) if the stellar mass dependence of |$\langle f_{\rm esc}^{\rm sf} \rangle$| is shallower than |$f_{\rm esc}^{\rm bh}$| , at z < 7 a transition stellar mass exists above which AGN dominate the escaping ionizing photon production rate; (iv) the transition stellar mass decreases with decreasing redshift. While AGN dominate the escaping emissivity above the knee of the stellar mass function at z ~ 6.8, they take-over at stellar masses that are a tenth of the knee mass by z = 4.

79 ASTRONOMY AND ASTROPHYSICS↗

The Survey of Lines in M31 (SLIM): The Drivers of the [C ii]/TIR Variation

The ratio of the [C ii] 158 μm emission line over the total infrared emission (TIR) is often used as a proxy for the photoelectric (PE) heating efficiency (ϵ{sub PE}) of the far-ultraviolet (FUV) photons absorbed by dust in the interstellar medium. In the nearby galaxy M31, we measure a strong radial variation of [C ii]/TIR that we rule out as being due to an intrinsic variation in ϵ{sub PE}. [C ii]/TIR fails as a proxy for ϵ{sub PE}, because the TIR measures all dust heating, not just the contribution from FUV photons capable of ejecting electrons from dust grains. Using extensive multi-wavelength coverage from the FUV to far-infrared, we infer the attenuated FUV emission (UV{sub att}), and the total attenuated flux (TOT{sub att}). We find [C ii]/TIR to be strongly correlated with UV{sub att}/TOT{sub att}, indicating that, in M31 at least, one of the dominant drivers for [C ii]/TIR variation is the relative hardness of the absorbed stellar radiation field. We define ϵ{sub PE}{sup UV}, [C ii]/UV{sub att} which should be more closely related to the actual PE efficiency, which we find to be essentially constant (1.85±0.8%) in all explored fields in M31. This suggests that part of the observed variation of [C ii]/TIR in other galaxies is likely due to a change in the relative hardness of the absorbed stellar radiation field, caused by a combination of variations in the stellar population, dust opacity, and galaxy metallicity, though PE efficiency may also vary across a wider range of environments.

79 ASTRONOMY AND ASTROPHYSICS↗

A DISK-WIND MODEL FOR THE NEAR-INFRARED EXCESS EMISSION IN PROTOSTARS

Protostellar systems, ranging from low-luminosity classical T Tauri and Herbig Ae stars to high-luminosity Herbig Be stars, exhibit a near-infrared (NIR) excess in their spectra that is dominated by a bump in the monochromatic luminosity with a peak near 3 {mu}m. The bump can be approximated by a thermal emission component of temperature {approx}1500 K that is of the order of the sublimation temperature of interstellar dust grains. In the currently popular 'puffed-up rim' scenario, the bump represents stellar radiation that propagates through the optically thin inner region of the surrounding accretion disk and is absorbed and reemitted by the dust that resides just beyond the dust sublimation radius r {sub sub}. However, this model cannot account for the strongest bumps measured in these sources, and it predicts a pronounced secondary bounce in the interferometric visibility curve that is not observed. In this paper we present an alternative interpretation, which attributes the bump to reemission of stellar radiation by dust that is uplifted from the disk by a centrifugally driven wind. Winds of this type are a leading candidate for the origin of the strong outflows associated with protostars, and there is observational evidence for disk winds originating on scales {approx}r {sub sub}. Using a newly constructed Monte Carlo radiative transfer code and focusing on low-luminosity sources, we show that this model can account for the NIR excess emission even in bright Herbig Ae stars such as AB Auriga and MWC 275, and that it successfully reproduces the basic features of the visibilities measured in these protostars. We argue that a robust dusty outflow in these sources could be self-limiting-through shielding of the stellar FUV photons-to a relatively narrow launching region between r {sub sub} and {approx}2 r {sub sub}. We also suggest that the NIR and scattered-light variability exhibited by a source like MWC 275 can be attributed in this picture to the uplifting of dust clouds from the disk.

79 ASTRONOMY AND ASTROPHYSICS↗

CO and [C ii] line emission of molecular clouds: the impact of stellar feedback and non-equilibrium chemistry

We analyse synthetic 12 CO, 13 CO, and [C ii] emission maps of molecular cloud (MC) simulations from the SILCC-Zoom project. We present radiation, magnetohydrodynamic zoom-in simulations of individual clouds, both with and without radiative stellar feedback, forming in a turbulent multiphase interstellar medium following on-the-fly the evolution of e.g. H 2 , CO, and C + . We introduce a novel post-processing routine based on cloudy which accounts for higher ionization states of carbon due to stellar radiation in H ii regions. Synthetic emission maps of [C ii] in and around feedback bubbles show that the bubbles are largely devoid of [C ii], as recently found in observations, which we attribute to the further ionization of C+ into C 2+ . For both 12 CO and 13 CO, the cloud-averaged luminosity ratio, $L_\rm {CO}/L_\rm {[C\, \small {II}]}$, can neither be used as a reliable measure of the H 2 mass fraction nor of the evolutionary stage of the clouds. We note a relation between the $I_\rm {CO}/I_\rm {[C\, \small {II}]}$ intensity ratio and the H 2 mass fraction for individual pixels of our synthetic maps. The scatter, however, is too large to reliably infer the H 2 mass fraction. Finally, the assumption of chemical equilibrium overestimates H2 and CO masses by up to 150 and 50 per cent, respectively, and $L_\rm {CO}$ by up to 60 per cent. The masses of H and C + would be underestimated by 65 and 30 per cent, respectively, and $L_\rm {[C\, \small {II}]}$ by up to 35 per cent. Hence, the assumption of chemical equilibrium in MC simulations introduces intrinsic errors of a factor of 2 in chemical abundances, luminosities, and luminosity ratios.

79 ASTRONOMY AND ASTROPHYSICS↗

ENERGY CONSERVATION AND GRAVITY WAVES IN SOUND-PROOF TREATMENTS OF STELLAR INTERIORS. II. LAGRANGIAN CONSTRAINED ANALYSIS

The speed of sound greatly exceeds typical flow velocities in many stellar and planetary interiors. To follow the slow evolution of subsonic motions, various sound-proof models attempt to remove fast acoustic waves while retaining stratified convection and buoyancy dynamics. In astrophysics, anelastic models typically receive the most attention in the class of sound-filtered stratified models. Generally, anelastic models remain valid in nearly adiabatically stratified regions like stellar convection zones, but may break down in strongly sub-adiabatic, stably stratified layers common in stellar radiative zones. However, studying stellar rotation, circulation, and dynamos requires understanding the complex coupling between convection and radiative zones, and this requires robust equations valid in both regimes. Here we extend the analysis of equation sets begun in Brown et al., which studied anelastic models, to two types of pseudo-incompressible models. This class of models has received attention in atmospheric applications, and more recently in studies of white-dwarf supernova progenitors. We demonstrate that one model conserves energy but the other does not. We use Lagrangian variational methods to extend the energy conserving model to a general equation of state, and dub the resulting equation set the generalized pseudo-incompressible (GPI) model. We show that the GPI equations suitably capture low-frequency phenomena in both convection and radiative zones in stars and other stratified systems, and we provide recommendations for converting low-Mach number codes to this equation set.

79 ASTRONOMY AND ASTROPHYSICS↗

Detailed opacity comparison for an improved stellar modeling of the envelopes of massive stars

Seismic observations have led to doubts or ambiguities concerning the opacity calculations used in stellar physics. Here, we concentrate on the iron-group opacity peak, due to iron, nickel, and chromium, located around T = 200,000 K for densities from 10{sup −8} to 10{sup −4} g cm{sup −3}, which creates some convective layers in stellar radiative envelopes for masses between 3 and 18 M{sub ⊙}. These conditions were extensively studied in the 1980s. More recently, inconsistencies between OP and OPAL opacity calculations have complicated the interpretation of seismic observations as the iron-group opacity peak excites acoustic and gravity modes in SPB, β Cephei, and sdB stars. We investigate the reliability of the theoretical opacity calculations using the modern opacity codes ATOMIC and SCO-RCG. We show their temperature and density dependence for conditions that are achievable in the laboratory and equivalent to astrophysical conditions. We also compare new theoretical opacity spectra with OP spectra and quantify how different approximations impact the Rosseland mean calculations.This detailed study estimates new ATOMIC and SCO-RCG Rosseland mean values for astrophysical conditions which we compare to OP values. Some puzzling questions are still under investigation for iron, but we find a strong increase in the Rosseland mean nickel opacity of a factor between 2 and 6 compared to OP. This appears to be due to the use of extrapolated atomic data for the Ni opacity within the OP calculations. A study on chromium is also shown.

79 ASTRONOMY AND ASTROPHYSICS↗

THE ORION FINGERS: NEAR-IR SPECTRAL IMAGING OF AN EXPLOSIVE OUTFLOW

We present near-IR (1.1–2.4 μ m) position–position–velocity cubes of the 500 year old Orion BN/KL explosive outflow with spatial resolution 1″ and spectral resolution 86 km s{sup −1}. We construct integrated intensity maps free of continuum sources of 15 H{sub 2} and [Fe ii] lines while preserving kinematic information of individual outflow features. Included in the detected H{sub 2} lines are the 1-0 S(1) and 1-0 Q(3) transitions, allowing extinction measurements across the outflow. Additionally, we present dereddened flux ratios for over two dozen outflow features to allow for the characterization of the true excitation conditions of the BN/KL outflow. All of the ratios show the dominance of the shock excitation of the H{sub 2} emission, although some features exhibit signs of fluorescent excitation from stellar radiation or J-type shocks. We also detect tracers of the PDR/ionization front north of the Trapezium stars in [O i] and [Fe ii] and analyze other observed outflows not associated with the BN/KL outflow.

79 ASTRONOMY AND ASTROPHYSICS↗

X-Ray-induced Deuterium Enrichment of N-rich Organics in Protoplanetary Disks: An Experimental Investigation Using Synchrotron Light

The deuterium enrichment of organics in the interstellar medium, protoplanetary disks, and meteorites has been proposed to be the result of ionizing radiation. The goal of this study is to simulate and quantify the effects of soft X-rays (0.1–2 keV), an important component of stellar radiation fields illuminating protoplanetary disks, on the refractory organics present in the disks. We prepared tholins, nitrogen-rich organic analogs to solids found in several astrophysical environments, e.g., Titan’s atmosphere, cometary surfaces, and protoplanetary disks, via plasma deposition. Controlled irradiation experiments with soft X-rays at 0.5 and 1.3 keV were performed at the SEXTANTS beamline of the SOLEIL synchrotron, and were immediately followed by ex-situ infrared, Raman, and isotopic diagnostics. Infrared spectroscopy revealed the preferential loss of singly bonded groups (N–H, C–H, and R–N≡C) and the formation of sp{sup 3} carbon defects with signatures at ∼1250–1300 cm{sup −1}. Raman analysis revealed that, while the length of polyaromatic units is only slightly modified, the introduction of defects leads to structural amorphization. Finally, tholins were measured via secondary ion mass spectrometry to quantify the D, H, and C elemental abundances in the irradiated versus non-irradiated areas. Isotopic analysis revealed that significant D-enrichment is induced by X-ray irradiation. Our results are compared to previous experimental studies involving the thermal degradation and electron irradiation of organics. The penetration depth of soft X-rays in μ m-sized tholins leads to volume rather than surface modifications: lower-energy X-rays (0.5 keV) induce a larger D-enrichment than 1.3 keV X-rays, reaching a plateau for doses larger than 5 × 10{sup 27} eV cm{sup −3}. Synchrotron fluences fall within the expected soft X-ray fluences in protoplanetary disks, and thus provide evidence of a new non-thermal pathway to deuterium fractionation of organic matter.

79 ASTRONOMY AND ASTROPHYSICS↗

The role of stellar feedback in the dynamics of H II regions

Stellar feedback is often cited as the biggest uncertainty in galaxy formation models today. This uncertainty stems from a dearth of observational constraints as well as the great dynamic range between the small scales (≲1 pc) where the feedback originates and the large scales of galaxies (≳1 kpc) that are shaped by this feedback. To bridge this divide, in this paper we aim to assess observationally the role of stellar feedback at the intermediate scales of H II regions (∼10-100 pc). In particular, we employ multiwavelength data to examine several stellar feedback mechanisms in a sample of 32 H II regions (with ages ∼3-10 Myr) in the Large and Small Magellanic Clouds, respectively. Using optical, infrared, radio, and X-ray images, we measure the pressures exerted on the shells from the direct stellar radiation, the dust-processed radiation, the warm ionized gas, and the hot X-ray-emitting gas. We find that the warm ionized gas dominates over the other terms in all of the sources, although two have comparable dust-processed radiation pressures to their warm gas pressures. The hot gas pressures are comparatively weak, while the direct radiation pressures are one to two orders of magnitude below the other terms. We discuss the implications of these results, particularly highlighting evidence for hot gas leakage from the H II shells and regarding the momentum deposition from the dust-processed radiation to the warm gas. Furthermore, we emphasize that similar observational work should be done on very young H II regions to test whether direct radiation pressure and hot gas can drive the dynamics at early times.

79 ASTRONOMY AND ASTROPHYSICS↗

A PRIMER ON UNIFYING DEBRIS DISK MORPHOLOGIES

A “minimum model” for debris disks consists of a narrow ring of parent bodies, secularly forced by a single planet on a possibly eccentric orbit, colliding to produce dust grains that are perturbed by stellar radiation pressure. We demonstrate how this minimum model can reproduce a wide variety of disk morphologies imaged in scattered starlight. Five broad categories of disk shape can be captured: “rings,” “needles,” “ships-and-wakes,” “bars,” and “moths (a.k.a. fans),” depending on the viewing geometry. Moths can also sport “double wings.” We explain the origin of morphological features from first principles, exploring the dependence on planet eccentricity, disk inclination dispersion, and the parent body orbital phases at which dust grains are born. A key determinant in disk appearance is the degree to which dust grain orbits are apsidally aligned. Our study of a simple steady-state (secularly relaxed) disk should serve as a reference for more detailed models tailored to individual systems. We use the intuition gained from our guidebook of disk morphologies to interpret, informally, the images of a number of real-world debris disks. These interpretations suggest that the farthest reaches of planetary systems are perturbed by eccentric planets, possibly just a few Earth masses each.

79 ASTRONOMY AND ASTROPHYSICS↗