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 · Search NASA
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.
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.
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.
A three-dimensional, multifluid simulation of a giant planet’s magnetospheric interaction with steady-state stellar wind from a Sun-like star was performed for four different orbital semimajor axes—10, 5, 1, and 0.2 au. We simulate the effect of the increasing, steady-state stellar wind pressure related to the planetary orbital semimajor axis on the global magnetospheric dynamics for a Saturn-like planet, including an Enceladus-like plasma torus. Mass-loss processes are shown to vary with orbital distance, with the centrifugal interchange instability displayed only in the 10 and 5 au cases, which reach a state of mass-loss equilibrium more slowly than the 1 or 0.2 au cases. The compression of the magnetosphere in the 1 and 0.2 au cases contributes to the quenching of the interchange process by increasing the ratio of total plasma thermal energy to corotational energy. The strength of field-aligned currents, associated with auroral radio emissions, is shown to increase in magnitude and latitudinal coverage with a corresponding shift equatorward from increased dynamic ram pressure experienced in the hotter orbits. Similar to observed hot Jovian planets, the warm exo-Saturn simulated in the current work shows enhanced ion density in the magnetosheath and magnetopause regions, as well as the plasma torus, which could contribute to altered transit signals, suggesting that for planets in warmer (>0.1 au) orbits, planetary magnetic field strengths and possibly exomoons—via the plasma torus—could be observable with future missions.
The measured spindown rates in quiescence of the transient accreting millisecond pulsars IGR J00291+5934, XTE J1751–305, SAX J1808.4–3658, and Swift J1756.9–2508 have been used to estimate the magnetic moments of these objects assuming standard magnetic dipole braking. It is shown that this approach leads to an overestimate if the amount of residual accretion is enough to distort the magnetosphere away from a force-free configuration through magnetospheric mass loading or crushing, so that the lever arm of the braking torque migrates inside the light cylinder. We derive an alternative spindown formula and calculate the residual accretion rates where the formula is applicable. As a demonstration we apply the alternative spindown formula to produce updated magnetic moment estimates for the four objects above. We note that based on current uncertain observations of quiescent accretion rates, magnetospheric mass loading and crushing are neither firmly indicated nor ruled out in these four objects. Because quiescent accretion rates are not measured directly (only upper limits are placed), without more data it is impossible to be confident about whether the thresholds for magnetospheric mass loading or crushing are reached or not.
We report on Nuclear Spectroscopic Telescope Array (NuSTAR) hard X-ray observations of the young rotation-powered radio pulsar PSR B1509−59 in the supernova remnant MSH 15−52. We confirm the previously reported curvature in the hard X-ray spectrum, showing that a log parabolic model provides a statistically superior fit to the spectrum compared with the standard power law. The log parabolic model describes the NuSTAR data, as well as previously published γ-ray data obtained with COMPTEL and AGILE, all together spanning 3 keV through 500 MeV. Our spectral modeling allows us to constrain the peak of the broadband high energy spectrum to be at 2.6 ± 0.8 MeV, an improvement of nearly an order of magnitude in precision over previous measurements. In addition, we calculate NuSTAR spectra in 26 pulse phase bins and confirm previously reported variations of photon indices with phase. Finally, we measure the pulsed fraction of PSR B1509−58 in the hard X-ray energy band for the first time. Using the energy resolved pulsed fraction results, we estimate that the pulsar’s off-pulse emission has a photon index value between 1.26 and 1.96. Our results support a model in which the pulsar’s lack of GeV emission is due to viewing geometry, with the X-rays originating from synchrotron emission from secondary pairs in the magnetosphere.
With an instantaneous view of 70% of the sky, the Fermi Gamma-ray Burst Monitor (GBM) is an excellent partner in the search for electromagnetic counterparts to gravitational-wave (GW) events. GBM observations at the time of the Laser Interferometer Gravitational-wave Observatory (LIGO) event GW150914 reveal the presence of a weak transient above 50 keV, 0.4 s after the GW event, with a false-alarm probability of 0.0022 (2.9σ). This weak transient lasting 1 s was not detected by any other instrument and does not appear to be connected with other previously known astrophysical, solar, terrestrial, or magnetospheric activity. Its localization is ill-constrained but consistent with the direction of GW150914. The duration and spectrum of the transient event are consistent with a weak short gamma-ray burst (GRB) arriving at a large angle to the direction in which Fermi was pointing where the GBM detector response is not optimal. If the GBM transient is associated with GW150914, then this electromagnetic signal from a stellar mass black hole binary merger is unexpected. We calculate a luminosity in hard X-ray emission between 1 keV and 10 MeV of 1.8{sub −1.0}{sup +1.5}×10{sup 49} erg s{sup −1}. Future joint observations of GW events by LIGO/Virgo and Fermi GBM could reveal whether the weak transient reported here is a plausible counterpart to GW150914 or a chance coincidence, and will further probe the connection between compact binary mergers and short GRBs.
We demonstrate that particles are regularly accelerated while experiencing curvature drift in flows driven by magnetic tension. Some examples of such flows include spontaneous turbulent reconnection and decaying magnetohydrodynamic turbulence, where a magnetic field relaxes to a lower-energy configuration and transfers part of its energy to kinetic motions of the fluid. We show that this energy transfer, which normally causes turbulent cascade and heating of the fluid, also results in a first-order acceleration of non-thermal particles. Since it is generic, this acceleration mechanism is likely to play a role in the production of non-thermal particle distribution in magnetically dominant environments such as the solar chromosphere, pulsar magnetospheres, jets from supermassive black holes, and γ-ray bursts.
Supermassive black holes (BHs) are believed to be the central powerhouse of active galactic nuclei. Applying the pulsar outer-magnetospheric particle accelerator theory to BH magnetospheres, we demonstrate that an electric field is exerted along the magnetic field lines near the event horizon of a rotating BH. In this particle accelerator (or a gap), electrons and positrons are created by photon–photon collisions and accelerated in the opposite directions by this electric field, efficiently emitting gamma-rays via curvature and inverse-Compton processes. It is shown that a gap arises around the null-charge surface formed by the frame-dragging effect, provided that there is no current injection across the gap boundaries. The gap is dissipating a part of the hole’s rotational energy, and the resultant gamma-ray luminosity increases with decreasing plasma accretion from the surroundings. Considering an extremely rotating supermassive BH, we show that such a gap reproduces the significant very-high-energy (VHE) gamma-ray flux observed from the radio galaxy IC 310, provided that the accretion rate becomes much less than the Eddington rate particularly during its flare phase. It is found that the curvature process dominates the inverse-Compton process in the magnetosphere of IC 310, and that the observed power-law-like spectrum in VHE gamma-rays can be explained to some extent by a superposition of the curvature emissions with varying curvature radius. It is predicted that the VHE spectrum extends into higher energies with increasing VHE photon flux.
The coronal magnetic field directly or indirectly affects a majority of the phenomena studied in the heliosphere. It provides energy for coronal heating, controls the release of coronal mass ejections, and drives heliospheric and magnetospheric activity, yet the coronal magnetic field itself has proven difficult to measure. This difficulty has prompted a decades-long effort to develop accurate, timely, models of the field—an effort that continues today. We have developed a method for improving global coronal magnetic field models by incorporating the type of morphological constraints that could be derived from coronal images. Here we report promising initial tests of this approach on two theoretical problems, and discuss opportunities for application.