Three-dimensional, Time-dependent MHD Simulation of Disk–Magnetosphere–Stellar Wind Interaction in a T Tauri, Protoplanetary System
Explore the source record for details and available documents.
SEARCH · Engineering Papers
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.
Explore the source record for details and available documents.
We analyze three epochs of ultraviolet, optical, and near-infrared (NIR) observations of the Taurus transitional disk GM Aur using the Hubble Space Telescope Imaging Spectrograph (STIS) and the Infrared Telescope Facility SpeX spectrograph. Observations were separated by one week and three months in order to study variability over multiple timescales. We calculate accretion rates for each epoch of observations using the STIS spectra and find that those separated by one week had similar accretion rates (∼1×10{sup −8} M{sub ⊙} yr{sup −1}) while the epoch obtained three months later had a substantially lower accretion rate (∼4×10{sup −9} M{sub ⊙} yr{sup −1}). We find that the decline in accretion rate is caused by lower densities of material in the accretion flows, as opposed to a lower surface coverage of the accretion columns. During the low accretion rate epoch, we also observe lower fluxes at both far-ultraviolet (FUV) and IR wavelengths, which trace molecular gas and dust in the disk, respectively. We find that this can be explained by a lower dust and gas mass in the inner disk. We attribute the observed variability to inhomogeneities in the inner disk, near the corotation radius, where gas and dust may co-exist near the footprints of the magnetospheric flows. These FUV–NIR data offer a new perspective on the structure of the inner disk, the stellar magnetosphere, and their interaction.
The theory of binary star formation predicts that close binaries ( a < 100 au) will experience periodic pulsed accretion events as streams of material form at the inner edge of a circumbinary disk (CBD), cross a dynamically cleared gap, and feed circumstellar disks or accrete directly onto the stars. The archetype for the pulsed accretion theory is the eccentric, short-period, classical T Tauri binary DQ Tau. Low-cadence (∼daily) broadband photometry has shown brightening events near most periastron passages, just as numerical simulations would predict for an eccentric binary. Magnetic reconnection events (flares) during the collision of stellar magnetospheres near periastron could, however, produce the same periodic, broadband behavior when observed at a one-day cadence. To reveal the dominant physical mechanism seen in DQ Tau’s low-cadence observations, we have obtained continuous, moderate-cadence, multiband photometry over 10 orbital periods, supplemented with 27 nights of minute-cadence photometry centered on four separate periastron passages. While both accretion and stellar flares are present, the dominant timescale and morphology of brightening events are characteristic of accretion. On average, the mass accretion rate increases by a factor of five near periastron, in good agreement with recent models. Large variability is observed in the morphology and amplitude of accretion events from orbit to orbit. We argue that this is due to the absence of stable circumstellar disks around each star, compounded by inhomogeneities at the inner edge of the CBD and within the accretion streams themselves. Quasiperiodic apastron accretion events are also observed, which are not predicted by binary accretion theory.
The transfer of matter between a circumbinary disk and a young binary system remains poorly understood, obscuring the interpretation of accretion indicators. To explore the behavior of these indicators in multiple systems, we have performed the first systematic time-domain study of young binaries in the ultraviolet. We obtained far- and near-ultraviolet HST/COS spectra of the young spectroscopic binaries DQ Tau and UZ Tau E. Here we focus on the continuum from 2800 to 3200 Å and on the C iv doublet (λλ1548.19, 1550.77 Å) as accretion diagnostics. Each system was observed over three or four consecutive binary orbits, at phases ∼0, 0.2, 0.5, and 0.7. Those observations are complemented by ground-based U-band measurements. Contrary to model predictions, we do not detect any clear correlation between accretion luminosity and phase. Further, we do not detect any correlation between C iv flux and phase. For both stars the appearance of the C iv line is similar to that of single Classical T Tauri Stars (CTTSs), despite the lack of stable long-lived circumstellar disks. However, unlike the case in single CTTSs, the narrow and broad components of the C iv lines are uncorrelated, and we argue that the narrow component is powered by processes other than accretion, such as flares in the stellar magnetospheres and/or enhanced activity in the upper atmosphere. We find that both stars contribute equally to the narrow component C iv flux in DQ Tau, but the primary dominates the narrow component C iv emission in UZ Tau E. The C iv broad component flux is correlated with other accretion indicators, suggesting an accretion origin. However, the line is blueshifted, which is inconsistent with its origin in an infall flow close to the star. It is possible that the complicated geometry of the region, as well as turbulence in the shock region, are responsible for the blueshifted line profiles.
We calculate the abundances of electrons and ions in the hot (≳500 K), dusty parts of protoplanetary disks, treating for the first time the effects of thermionic and ion emission from the dust grains. High-temperature ionization modeling has involved simply assuming that alkali elements such as potassium occur as gas-phase atoms and are collisionally ionized following the Saha equation. We show that the Saha equation often does not hold, because free charges are produced by thermionic and ion emission and destroyed when they stick to grain surfaces. This means the ionization state depends not on the first ionization potential of the alkali atoms, but rather on the grains’ work functions. The charged species’ abundances typically rise abruptly above about 800 K, with little qualitative dependence on the work function, gas density, or dust-to-gas mass ratio. Applying our results, we find that protoplanetary disks’ dead zone, where high diffusivities stifle magnetorotational turbulence, has its inner edge located where the temperature exceeds a threshold value ≈1000 K. The threshold is set by ambipolar diffusion except at the highest densities, where it is set by Ohmic resistivity. We find that the disk gas can be diffusively loaded onto the stellar magnetosphere at temperatures below a similar threshold. We investigate whether the “short-circuit” instability of current sheets can operate in disks and find that it cannot, or works only in a narrow range of conditions; it appears not to be the chondrule formation mechanism. We also suggest that thermionic emission is important for determining the rate of Ohmic heating in hot Jupiters.
We propose a generalized model for stellar spin-down, disk accretion, and truncation, and the origin of winds, jets, and bipolar outflows from young stellar objects. We consider the steady state dynamics of accretion of matter from a viscous and imperfectly conducting disk onto a young star with a strong magnetic field. For an aligned stellar magnetosphere, shielding currents in the surface layers of the disk prevent stellar field lines from penetrating the disk everywhere except for a range of radii about pi = R(sub x), where the Keplerian angular speed of rotation Omega(sub x) equals the angular speed of the star Omega(sub *). For the low disk accretion rates and high magnetic fields associated with typical T Tauri stars, R(sub x) exceeds the radius of the star R(sub *) by a factor of a few, and the inner disk is effectively truncated at a radius R(sub t) somewhat smaller than R(sub x). Where the closed field lines between R(sub t) and R(sub x) bow sufficiently inward, the accreting gas attaches itself to the field and is funneled dynamically down the effective potential (gravitational plus centrifugal) onto the star. Contrary to common belief, the accompanying magnetic torques associated with this accreting gas may transfer angular momentum mostly to the disk rather than to the star. Thus, the star can spin slowly as long as R(sub x) remains significantly greater than R(sub *). Exterior to R(sub x) field lines threading the disk bow outward, which makes the gas off the mid-plane rotate at super-Keplerian velocities. This combination drives a magnetocentrifugal wind with a mass-loss rate M(sub w) equal to a definite fraction f of the disk accretion rate M(sub D). For high disk accretion rates, R(sub x) is forced down to the stellar surface, the star is spun to breakup, and the wind is generated in a manner identical to that proposed by Shu, Lizano, Ruden, & Najita in a previous communication to this journal. In two companion papers (II and III), we develop a detailed but idealized theory of the magnetocentrifugal acceleration process.
Low-mass population III (PopIII) stars of ≲0.8 M{sub ⊙} could survive up until the present. The nondetection of low-mass PopIII stars in our Galaxy has already put a stringent constraint on the initial mass function (IMF) of PopIII stars, suggesting that PopIII stars have a top-heavy IMF. On the other hand, some claim that the lack of such stars stems from metal enrichment of their surfaces by the accretion of heavy elements from the interstellar medium (ISM). We investigate the effects of the stellar wind on metal accretion onto low-mass PopIII stars because accretion of the local ISM onto the Sun is prevented by the solar wind, even for neutrals. The stellar wind and radiation of low-mass PopIII stars are modeled based on knowledge of nearby low-mass stellar systems, including our Sun. We find that low-mass PopIII stars traveling across the Galaxy form a stellar magnetosphere in most of their life. Once the magnetosphere is formed, most of the neutral interstellar particles are photoionized before reaching the stellar surface and are blown away by the wind. Especially, the accretion abundance of iron will be reduced by a factor of <10{sup −12} compared with Bondi–Hoyle–Lyttleton accretion. The metal accretion can enhance iron abundance [Fe/H] only up to ∼−14. This demonstrates that low-mass PopIII stars remain pristine and will be found as metal-free stars and that further searches for them are valuable in constraining the IMF of PopIII stars.
We consider the structure and emission spectra of nonthermal accretion disks around both strongly and weakly magnetized neutron stars. Such disks may be dissipating their gravitational binding energy and transferring their angular momentum via semicontinuous magnetic reconnections. We consider specifically the structure of the disk-stellar magnetospheric boundary where magnetic pressure balances the disk pressure. We consider energy dissipation via reconnection of the stellar field and small-scale disk turbulent fields of opposite polarity. Constraints on the disk emission spectrum are discussed.
We argue that the strong emission lines of T Tauri stars are generally produced in infalling envelopes. Simple models of infall constrained to a dipolar magnetic field geometry explain many peculiarities of observed line profiles that are difficult, if not impossible, to reproduce with wind models. Radiative transfer effects explain why certain lines can appear quite symmetric while other lines simultaneously exhibit inverse P Cygni profiles, without recourse to complicated velocity fields. The success of the infall models in accounting for qualitative features of observed line profiles supports the proposal that stellar magnetospheres disrupt disk accretion in T Tauri stars, that true boundary layers are not usually present in T Tauri stars, and that the observed 'blue veiling' emission arises from the base of the magnetospheric accretion column.
We investigate how the radio emission of a pulsar interacts with plasma derived from a stellar companion. Various physical mechanisms that can cause radio pulse eclipse are discussed, and predictions are made for the polarization properties of the emergent radio wave. We consider eclipses by a wind from the stellar companion, by a stellar magnetosphere, or by material entrained in the pulsar wind. Eclipses due to refraction require either a relatively high plasma density or a sharp edge to the plasma distribution. The conditions that must prevail for free-free absorption to be effective in eclipsing a radio beam are also outlined. Pulse smearing may be important at higher frequencies; related eclipse mechanisms include pulse spreading due to a rapidly changing electron column, and scattering by Langmuir turbulence. The high brightness temperature radio beam can generate its own plasma turbulence via a number of nonlinear parametric instabilities, such as the instability associated with stimulated Raman scattering. When the plasma turbulence is heavily damped, the radio bean can still undergo induced Compton scattering. Stimulated scattering effects such as these are very sensitive to the presence of narrow-band substructure in the pulsar radio emission. Finally, we consider the possibility that plasma derived from a stellar companion may mix with the relativistic pulsar wind and cause cyclotron absorption at low radio frequencies. Even if the cyclotron optical depth is small, fluctuations in the emergent polarization of the radio beam on the timescale of a few seconds are a very sensitive probe of the spatial structure of the magnetic field in the pulsar wind. The current observational properties of two known eclipsing pulsar systems, PSR 1957+20 and PSR 1744-24A, are used to construct tentative eclipse models. The favored model for PSR 1957+20 is cyclotron or synchrotron absorption by plasma embedded in the pulsar wind combined with pulse smearing at high frequency, and the favored model for PSR 1744-24A is backscattering off plasma turbulence generated by the stimulated Raman scattering parametric instability. Pulsar eclipses promise to provide a good diagnostic of pulsar winds and possible of the pulse emission mechanism.
Magnetars are the most highly magnetized neutron stars in the cosmos (with magnetic field 10 13 –10 15 G). Giant flares from magnetars are rare, short-duration (about 0.1 s) bursts of hard X-rays and soft γ rays. Owing to the limited sensitivity and energy coverage of previous telescopes, no magnetar giant flare has been detected at gigaelectronvolt (GeV) energies. Here, we report the discovery of GeV emission from a magnetar giant flare on 15 April 2020 (refs. 3,4 and A. J. Castro-Tirado et al., manuscript in preparation). The Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope detected GeV γ rays from 19 s until 284 s after the initial detection of a signal in the megaelectronvolt (MeV) band. Our analysis shows that these γ rays are spatially associated with the nearby (3.5 megaparsecs) Sculptor galaxy and are unlikely to originate from a cosmological γ-ray burst. Thus, we infer that the γ rays originated with the magnetar giant flare in Sculptor. Furthermore, we suggest that the GeV signal is generated by an ultra-relativistic outflow that first radiates the prompt MeV-band photons, and then deposits its energy far from the stellar magnetosphere. After a propagation delay, the outflow interacts with environmental gas and produces shock waves that accelerate electrons to very high energies; these electrons then emit GeV γ rays as optically thin synchrotron radiation. This observation implies that a relativistic outflow is associated with the magnetar giant flare, and suggests the possibility that magnetars can power some short γ-ray bursts.
Thermal X ray sources associated with rotating collapsed stars with surrounding plasma shells, discussing plasma density profile and electron distribution in stellar magnetosphere
This paper reports the observation and interpretation of the ultraviolet spectrum variations of nine helium-strong stars: HD 36485, 37017, 37479, 37776, 58260, 60344, 64740, 96446, and 133518. A unified model is developed to account for the observed correlation among three stellar properties: the line profile characteristics of the C IV and Si IV resonance doublets, the variations in the strength of these lines, and the inferred magnetic field geometry. It is proposed that circumstellar plasma is trapped in the stellar magnetosphere near the magnetic equator or is channeled to form jetlike outflows from the magnetic polar regions. These results, together with those of a previous study of the helium-weak sn stars, show that both helium-weak and helium-strong stars can possess magnetospherically trapped plasma, notwithstanding their different photospheric properties. New results for radii and temperatures of the helium-strong stars in Ori OB 1 and for HD 64740 from low-dispersion IUE spectra are also presented.
Among young binary stars whose magnetospheres are expected to collide, only two systems have been observed near periastron in the X-ray band: the low-mass DQ Tau and the older and more massive HD 152404. Both exhibit elevated levels of X-ray emission at periastron. Our goal is to determine whether colliding magnetospheres in young high-eccentricity binaries commonly produce elevated average levels of X-ray activity. This work is based on Chandra snapshots of multiple periastron and non-periastron passages in four nearby young eccentric binaries (Parenago 523, RX J1622.7-2325 Nw, UZ Tau E, and HD 152404). We find that for the merged sample of all four binaries the current X-ray data show an increasing average X-ray flux near periastron (at a ∼2.5-sigma level). Further comparison of these data with the X-ray properties of hundreds of young stars in the Orion Nebula Cluster, produced by the Chandra Orion Ultradeep Project (COUP), indicates that the X-ray emission from the merged sample of our binaries cannot be explained within the framework of the COUP-like X-ray activity. However, due to the inhomogeneities of the merged binary sample and the relatively low statistical significance of the detected flux increase, these findings are regarded as tentative only. More data are needed to prove that the flux increase is real and is related to the processes of colliding magnetospheres.
We test analytic predictions from different models of magnetospheric accretion, which invoke disk-locking, using stellar and accretion parameters derived from models of low-resolution optical spectra of 36 T Tauri stars (TTSs) in NGC 2264 (age {approx}3 Myr). Little evidence is found for models that assume purely dipolar field geometries; however, strong support is found in the data for a modified version of the X-wind model which allows for non-dipolar field geometries. The trapped flux concept in the X-wind model is key to making the analytic predictions which appear supported in the data. By extension, our analysis provides support for the outflows predicted by the X-wind as these also originate in the trapped flux region. In addition, we find no support in the data for accretion-powered stellar winds from young stars. By comparing the analysis presented here of NGC 2264 with a similar analysis of stars in Taurus (age {approx}1-2 Myr), we find evidence that the equilibrium interaction between the magnetic field and accretion disk in TTS systems evolves as the stars grow older, perhaps as the result of evolution of the stellar magnetic field geometry. We compare the accretion rates we derive with accretion rates based on U-band excess, finding good agreement. In addition, we use our accretion parameters to determine the relationship between accretion and H{beta} luminosity, again finding good agreement with previously published results; however, we also find that care must be taken when applying this relationship due to strong chromospheric emission in young stars, which can lead to erroneous results in some cases.
Global-scale solutions for the magnetosphere of a pulsar consist of a region of low-lying, closed magnetic field near the star, bounded by opposite-polarity regions of open magnetic field along which the pulsar wind flows into space. Separating these open-field regions is a magnetic discontinuity—an electric current sheet—consisting of generally nonneutral plasma. We have developed a self-consistent model for the internal equilibrium structure of the sheet by generalizing the charge-neutral Vlasov/Maxwell equilibria of Harris and Hoh to allow for net electric charge. The resulting equations for the electromagnetic field are solved analytically and numerically. Our results show that the internal thermal pressure needed to establish equilibrium force balance, and the associated effective current-sheet thickness and magnetization, can differ by orders of magnitude from the Harris/Hoh charge-neutral limit. The new model provides a starting point for kinetic or fluid investigations of instabilities that can cause magnetic reconnection and flaring in pulsar magnetospheres.
We perform ''first-principles'' relativistic particle-in-cell simulations of aligned pulsar magnetosphere. We allow free escape of particles from the surface of a neutron star and continuously populate the magnetosphere with neutral pair plasma to imitate pair production. As pair plasma supply increases, we observe the transition from a charge-separated ''electrosphere'' solution with trapped plasma and no spin-down to a solution close to the ideal force-free magnetosphere with electromagnetically dominated pulsar wind. We calculate the magnetospheric structure, current distribution, and spin-down power of the neutron star. We also discuss particle acceleration in the equatorial current sheet.
We study the efficiency of pair production in polar caps of young pulsars under a variety of conditions to estimate the maximum possible multiplicity of pair plasma in pulsar magnetospheres. We develop a semi-analytic model for calculation of cascade multiplicity which allows efficient exploration of the parameter space and corroborate it with direct numerical simulations. Pair creation processes are considered separately from particle acceleration in order to assess different factors affecting cascade efficiency, with acceleration of primary particles described by recent self-consistent non-stationary model of pair cascades. We argue that the most efficient cascades operate in the curvature radiation/synchrotron regime, the maximum multiplicity of pair plasma in pulsar magnetospheres is ∼few × 10{sup 5}. The multiplicity of pair plasma in magnetospheres of young energetic pulsars weakly depends on the strength of the magnetic field and the radius of curvature of magnetic field lines and has a stronger dependence on pulsar inclination angle. This result questions assumptions about very high pair plasma multiplicity in theories of pulsar wind nebulae.