Experiments for Photoionized Plasma Emission from Accretion-Powered X-Ray Sources.
Abstract not provided.
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.
Abstract not provided.
We present global 3D radiation magnetohydrodynamic simulations of accretion onto a 6.62 solar-mass black hole, with quasi-steady-state accretion rates reaching 0.016–0.9 times the critical accretion rate, which is defined as the accretion rate for powering the Eddington luminosity, assuming a 10% radiative efficiency, in three different runs. The simulations show no sign of thermal instability over hundreds of thermal timescales at 10 r g . The energy dissipation occurs close to the mid-plane in the near-critical runs and near the disk surface in the low–accretion rate run. The total radiative luminosity inside ~20 r g is about 1%–30% of the Eddington limit, with radiative efficiencies of about 6% and 3%, respectively, in the sub- and near-critical accretion regimes. In both cases, self-consistent turbulence generated by the magnetorotational instability leads to angular momentum transfer, and the disk is supported by magnetic pressure. Outflows from the central low-density funnel, with a terminal velocity of ~0.1c, are seen only in the near-critical runs. We conclude that these magnetic pressure–dominated disks are thermally stable and thicker than the α disk, and that the effective temperature profiles are much flatter than those in the α disks. The magnetic pressures of these disks are comparable within an order of magnitude to the previous analytical magnetic pressure–dominated disk model.
Fueling and feedback couple supermassive black holes (SMBHs) to their host galaxies across many orders of magnitude in spatial and temporal scales, making this problem notoriously challenging to simulate. We use a multi-zone computational method based on the general relativistic magnetohydrodynamic (GRMHD) code KHARMA that allows us to span 7 orders of magnitude in spatial scale, to simulate accretion onto a non-spinning SMBH from an external medium with a Bondi radius of R B ≈ 2 × 10 5 GM•/c 2 , where M• is the SMBH mass. For the classic idealized Bondi problem, spherical gas accretion without magnetic fields, our simulation results agree very well with the general relativistic analytic solution. Meanwhile, when the accreting gas is magnetized, the SMBH magnetosphere becomes saturated with a strong magnetic field. The density profile varies as ~r -1 rather than r -3/2 and the accretion rate $\dot{M}$ is consequently suppressed by over 2 orders of magnitude below the Bondi rate $\dot{M}$ B . Here, we find continuous energy feedback from the accretion flow to the external medium at a level of $\thicksim$ 10 -2 $\dot{M}$c 2 $\thicksim$ 5 X 10 -5 $\dot{M}$ B c 2 . Energy transport across these widely disparate scales occurs via turbulent convection triggered by magnetic field reconnection near the SMBH. Thus, strong magnetic fields that accumulate on horizon scales transform the flow dynamics far from the SMBH and naturally explain observed extremely low accretion rates compared to the Bondi rate, as well as at least part of the energy feedback.
We present a comprehensive analysis of super-Eddington black hole accretion simulations that solve the GRMHD equations coupled with angle-discretized radiation transport. The simulations span a range of accretion rates, two black hole spins, and two magnetic field topologies, and include resolution studies as well as comparisons with nonradiative models. Super-Eddington accretion flows consistently develop geometrically thick disks supported by radiation pressure, regardless of magnetic field configuration. Radiation generated in the inner disk drives substantial outflows, forming conical funnel regions that limit photon escape and result in very low radiation efficiency. The accretion flows are highly turbulent, with thermal energy transport dominated by radiation advection rather than diffusion. Angular momentum is primarily carried outward by Maxwell stress, with turbulent Reynolds stress playing a subdominant role. Both strong and weak jets are produced. Strong jets arise from sufficient net vertical magnetic flux and rapid black hole spin, and they can effectively evacuate the funnel, enabling radiation to escape through strong geometric beaming. In contrast, weak jets fail to clear the funnel, which becomes obscured by radiation-driven outflows and leads to distinct observational signatures. Spiral structures are observed in the plunging region, behaving like density waves. These super-Eddington models are applicable to a variety of astronomical systems, including ultraluminous X-ray sources, little red dots, and black hole transients.
Recent observations of the Milky Way and galaxies at high redshifts suggest that galaxy discs were already in place soon after the Big Bang. While the gas infall history of the Milky Way in the inner disc has long been assumed to be characterised by a short accretion timescale, this has not been directly constrained using observations. Using data for the inner regions of the Milky Way recently produced by APOGEE and Gaia and of unprecedented quantity and quality, we aim to derive strong constraints on the infall history of the inner (< 6 kpc) Galaxy (with a focus on stars between 4 and 6 kpc, which we show is an appropriate proxy for the entire inner disc). We implemented gas infall into a chemical evolution model of the Galaxy disc, and used a Schmidt–Kennicutt law to connect the infall to the star formation. We explore a number of models, and two different formulations of the infall law. In one formulation, the infall is non-parametric, and in the other the infall has an explicitly exponential form. We fit the model parameters to the time–[Si/Fe] distribution of solar vicinity stars, and the metallicity and [Si/Fe] distribution function of stars with a galactocentric radius of between 4 and 6 kpc from APOGEE. Our results point to a fast, early gas accretion, and an upper limit on the accretion timescale of around 2 Gyr in the inner disc of the Milky Way. This suggests that at least half the baryons were in place within 2–3 Gyr of the Big Bang, and that half the stars of the inner disc formed within the first 5 Gyr, during the thick disc formation phase. This implies that the stellar mass of the inner disc is dominated by the thick disc, supporting our previous work, and that the gas accretion onto the inner disc was rapid and early.
We find that the proton separation energy, S(p), of 73 Rb is –640(40) keV, deduced from the observation of β-delayed ground-state protons following the decay of 73 Sr. This lower-limit determination of the proton separation energy of 73 Rb coupled with previous upper limits from nonobservation, provides a full constraint on the mass excess with ΔM ( 73 Rb) = –46.01 ± 0.04 MeV. With this new mass excess and the excitation energy of the J π = 5/2 – isobaric-analog state (T = 3/2) in 73 Rb, an improved constraint can be put on the mass excess of 73Sr using the isobaric-multiplet mass equation (IMME), and we find ΔM( 73 Sr) = –31.98 ± 0.37 MeV. Furthermore, these new data were then used to study the composition of ashes on accreting neutron stars following Type I x-ray bursts. Counterintuitively, we find that there should be an enhanced fraction of A > 102 nuclei with more negative proton separation energies at the 72 Kr rp-process waiting point. Larger impurities of heavier nuclei in the ashes of accreting neutron stars will impact the cooling models for such astrophysical scenarios.
Here, we present a new mechanism for deep crustal heating in accreting neutron stars. Charged pions (π + ) are produced in nuclear collisions on the neutron star surface during active accretion and upon decay they provide a flux of neutrinos into the neutron star crust. For massive and/or compact neutron stars, neutrinos deposit ≈ 1 – 2 MeV of heat per accreted nucleon into the inner crust. The strength of neutrino heating is comparable to the previously known sources of deep crustal heating, such as from pycnonuclear fusion reactions, and is relevant for studies of cooling neutron stars. We model the thermal evolution of a transient neutron star in a low-mass x-ray binary, and in the particular case of the neutron star MXB 1659-29 we show that additional deep crustal heating requires a higher thermal conductivity for the neutron star inner crust. A better knowledge of pion-production cross sections near threshold would improve the accuracy of our predictions.
The grain growth process in the Caltech water–ice dusty plasma experiment has been studied using a high-speed camera and a long-distance microscope lens. It is observed that (i) the ice grain number density decreases fourfold as the average grain major axis increases from 20 to 80 μ m, (ii) the major axis length has a log-normal distribution rather than a power-law dependence, and (iii) no collisions between ice grains are apparent. The grains have a large negative charge resulting in strong mutual repulsion and this, combined with the fractal character of the ice grains, prevents them from agglomerating. In order for the grain kinetic energy to be sufficiently small to prevent collisions between ice grains, the volumetric packing factor (i.e., ratio of the actual volume to the volume of a circumscribing ellipsoid) of the ice grains must be less than ∼0.1 depending on the exact relative velocity of the grains in question. Thus, it is concluded that direct accretion of water molecules is very likely to dominate the observed ice grain growth.
X-ray reverberation, which exploits the time delays between variability in different energy bands as a function of Fourier frequency, probes the structure of the inner accretion disks and X-ray coronae of active galactic nuclei. We present a systematic X-ray spectroscopic and reverberation study of the high-Eddington-ratio narrow-line Seyfert 1 galaxy Ark 564, using over 900 ks of XMM-Newton and NuSTAR observations spanning 13 yr. The time-averaged spectra can be well described by a model consisting of a coronal continuum, relativistic disk reflection, warm Comptonization, and three warm absorbers. Leveraging the high X-ray brightness of Ark 564, we are able to resolve the time evolution of the spectra and contemporaneous reverberation lags. The soft-band lag relative to the continuum increases with the X-ray flux, while Fe Kα lags are detected in only a subset of epochs and do not correlate with soft lags. Models based on a lamppost corona and reflection from a standard thin disk can broadly reproduce the observed lag-energy spectra of low-flux epochs; however, additional reverberation from the warm Comptonized atmosphere is required to explain the soft lags observed in high-flux epochs. A vertically puffed-up inner disk and a variable, vertically extended corona can better explain the observed evolution of the lags and covariance spectra. Our study underscores the importance of multiepoch, multiband analyses for a comprehensive understanding the inner accretion disk and corona.
One of the fundamental properties of an intermediate polar is the dynamical nature of the accretion flow as it encounters the white dwarf’s (WD’s) magnetosphere. Many works have presumed a dichotomy between disk-fed accretion, in which the WD accretes from a Keplerian disk, and stream-fed accretion, in which the matter stream from the donor star directly impacts the WD’s magnetosphere without forming a disk. However, there is also a third, poorly understood regime in which the accretion flow consists of a torus of diamagnetic blobs that encircles the WD. This mode of accretion is expected to exist at mass-transfer rates below those observed during disk-fed accretion, but above those observed during pure stream-fed accretion. We invoke the diamagnetic-blob regime to explain the exceptional Transiting Exoplanet Survey Satellite light curve of the intermediate polar TX Col, which transitioned into and out of states of enhanced accretion during Cycles 1 and 3. Power-spectral analysis reveals that the accretion was principally stream fed. However, when the mass-transfer rate spiked, large-amplitude quasi-periodic oscillations (QPOs) abruptly appeared and dominated the light curve for weeks. The QPOs have two striking properties: they appear in a stream-fed geometry at elevated accretion rates, and they occur preferentially within a well-defined range of frequencies (∼10–25 cycle day{sup −1}). We propose that during episodes of enhanced accretion, a torus of diamagnetic blobs forms near the binary’s circularization radius and that the QPOs are beats between the white dwarf’s spin frequency and unstable blob orbits within the WD’s magnetosphere. We discuss how such a torus could be a critical step in producing an accretion disk in a formerly diskless system.
TWA 3A is the most recent addition to a small group of young binary systems that both actively accrete from a circumbinary disk and have spectroscopic orbital solutions. As such, it provides a unique opportunity to test binary accretion theory in a well-constrained setting. To examine TWA 3A’s time-variable accretion behavior, we have conducted a two-year, optical photometric monitoring campaign, obtaining dense orbital phase coverage (∼20 observations per orbit) for ∼15 orbital periods. From U -band measurements we derive the time-dependent binary mass accretion rate, finding bursts of accretion near each periastron passage. On average, these enhanced accretion events evolve over orbital phases 0.85 to 1.05, reaching their peak at periastron. The specific accretion rate increases above the quiescent value by a factor of ∼4 on average but the peak can be as high as an order of magnitude in a given orbit. The phase dependence and amplitude of TWA 3A accretion is in good agreement with numerical simulations of binary accretion with similar orbital parameters. In these simulations, periastron accretion bursts are fueled by periodic streams of material from the circumbinary disk that are driven by the binary orbit. We find that TWA 3A’s average accretion behavior is remarkably similar to DQ Tau, another T Tauri binary with similar orbital parameters, but with significantly less variability from orbit to orbit. This is only the second clear case of orbital-phase-dependent accretion in a T Tauri binary.
We calculate the evolution of gas giant planets during the runaway gas accretion phase of formation, to understand how the luminosity of young giant planets depends on the accretion conditions. We construct steady-state envelope models, and run time-dependent simulations of accreting planets with the code Modules for Experiments in Stellar Astrophysics. We show that the evolution of the internal entropy depends on the contrast between the internal adiabat and the entropy of the accreted material, parametrized by the shock temperature T {sub 0} and pressure P {sub 0}. At low temperatures (T{sub 0}≲300–1000 K, depending on model parameters), the accreted material has a lower entropy than the interior. The convection zone extends to the surface and can drive a high luminosity, leading to rapid cooling and cold starts. For higher temperatures, the accreted material has a higher entropy than the interior, giving a radiative zone that stalls cooling. For T{sub 0}≳2000 K, the surface–interior entropy contrast cannot be accommodated by the radiative envelope, and the accreted matter accumulates with high entropy, forming a hot start. The final state of the planet depends on the shock temperature, accretion rate, and starting entropy at the onset of runaway accretion. Cold starts with L≲5×10{sup −6} L{sub ⊙} require low accretion rates and starting entropy, and the temperature of the accreting material needs to be maintained close to the nebula temperature. If instead the temperature is near the value required to radiate the accretion luminosity, 4πR{sup 2}σT{sub 0}{sup 4}∼(GM M-dot /R), as suggested by previous work on radiative shocks in the context of star formation, gas giant planets form in a hot start with L∼10{sup −4} L{sub ⊙}.
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.
A millisecond pulsar is a neutron star that has been substantially spun up by accretion from a binary companion. A previously unrecognized factor governing the spin evolution of such pulsars is the crucial effect of nonsteady or transient accretion. We numerically compute the evolution of accreting neutron stars through a series of outburst and quiescent phases, considering the drastic variation of the accretion rate and the standard disk–magnetosphere interaction. We find that, for the same long-term average accretion rate, X-ray transients can spin up pulsars to rates several times higher than can persistent accretors, even when the spin-down due to electromagnetic radiation during quiescence is included. We also compute an analytical expression for the equilibrium spin frequency in transients, by taking spin equilibrium to mean that no net angular momentum is transferred to the neutron star in each outburst cycle. We find that the equilibrium spin rate for transients, which depends on the peak accretion rate during outbursts, can be much higher than that for persistent sources. This explains our numerical finding. This finding implies that any meaningful study of neutron star spin and magnetic field distributions requires the inclusion of the transient accretion effect, since most accreting neutron star sources are transients. Our finding also implies the existence of a submillisecond pulsar population, which is not observed. This may point to the need for a competing spin-down mechanism for the fastest-rotating accreting pulsars, such as gravitational radiation.
It has been realized in recent years that the accretion of pebble-sized dust particles onto planetary cores is an important mode of core growth, which enables the formation of giant planets at large distances and assists planet formation in general. The pebble accretion theory is built upon the orbit theory of dust particles in a laminar protoplanetary disk (PPD). For sufficiently large core mass (in the “Hill regime”), essentially all particles of appropriate sizes entering the Hill sphere can be captured. However, the outer regions of PPDs are expected to be weakly turbulent due to the magnetorotational instability (MRI), where turbulent stirring of particle orbits may affect the efficiency of pebble accretion. We conduct shearing-box simulations of pebble accretion with different levels of MRI turbulence (strongly turbulent assuming ideal magnetohydrodynamics, weakly turbulent in the presence of ambipolar diffusion, and laminar) and different core masses to test the efficiency of pebble accretion at a microphysical level. We find that accretion remains efficient for marginally coupled particles (dimensionless stopping time τ{sub s}∼0.1--1) even in the presence of strong MRI turbulence. Though more dust particles are brought toward the core by the turbulence, this effect is largely canceled by a reduction in accretion probability. As a result, the overall effect of turbulence on the accretion rate is mainly reflected in the changes in the thickness of the dust layer. On the other hand, we find that the efficiency of pebble accretion for strongly coupled particles (down to τ{sub s}∼0.01) can be modestly reduced by strong turbulence for low-mass cores.
We introduce simulations of accretion on to point masses embedded in an isothermal gas that is magnetized and supersonically turbulent, as occurs for protostars in molecular clouds. We use the ORION2 adaptive mesh refinement (AMR) code to carry out ideal magnetohydrodynamic simulations for an rms Mach number $\mathcal{M}$ rms =5 and a wide range of Alfvén Mach numbers. We find that the probability density functions for accretion rates in all models are very wide (±0.5 dex) and asymmetric, tilted to low accretion rates; the mean accretion rate is about twice the median. We find that the results of Lee et al. for magnetized Bondi-Hoyle accretion with the relative velocity parallel to the field describe our results to within a factor of 2, and we suggest that this should be valid at least for $\mathcal{M}$ rms ≲ 10. Our results demonstrate that turbulent magnetic fields of the strength observed in molecular clouds reduce the accretion rate relative to the classical Bondi-Hoyle rate by a factor of a few for Alfvén Mach numbers of order unity, but this is comparable to the reduction due to supersonic hydrodynamic turbulence alone. This reduction in accretion rates should be taken into account in analytic models of competitive accretion and analytic estimates of the accretion luminosities of young stellar objects in molecular clouds.
We report on our continuing studies of Classical Nova explosions by following the evolution of thermonuclear runaways (TNRs) on carbon-oxygen (CO) white dwarfs (WDs). We have varied both the mass of the WD and the composition of the accreted material. Rather than assuming that the material has mixed from the beginning, we now rely on the results of the multidimensional (multi-D) studies of mixing as a consequence of the TNRs in WDs that accreted only Solar matter. The multi-D studies find that mixing with the core occurs after the TNR is well underway and reach enrichment levels in agreement with observations of the ejecta abundances. We report on 3 studies in this paper. First, simulations in which we accrete only Solar matter with NOVA (our 1-D, fully implicit, hydro code). Second, we use MESA for similar studies in which we accrete only Solar material and compare the results. Third, we accrete Solar matter until the TNR is initiated and then switch the composition in the accreted layers to a mixed composition: either 25% core and 75% Solar or 50% core and 50% Solar. The amount of accreted material is inversely proportional to the initial 12 C abundance. Thus, accreting Solar material results in more material to fuel the outburst - much larger than in the earlier studies where mixed materials were used from the beginning. We tabulate the amount of ejected gases, their velocities, and abundances. We predict the amount of 7 Li and 7 Be produced and ejected by the explosion and compare our predictions to our Large Binocular Telescope (LBT) high dispersion spectra which determined the abundance of 7 Li in nova V5668 Sgr. Finally, many of these simulations eject significantly less mass than accreted and, therefore, the WD is growing in mass toward the Chandrasekhar Limit.
Recent discoveries of young exoplanets within their natal disks offer exciting opportunities to study ongoing planet formation. In particular, a planet’s mass accretion rate can be constrained by observing the accretion-induced excess emission. So far, planetary accretion is only probed by the Hα line, which is then converted to a total accretion luminosity using correlations derived for stars. However, the majority of the accretion luminosity is expected to emerge from hydrogen continuum emission, and is best measured in the ultraviolet (UV). In this paper, we present HST/WFC3/UVIS F336W (UV) and F656N (Hα) high-contrast imaging observations of PDS 70. Applying a suite of novel observational techniques, we detect the planet PDS 70 b with signal-to-noise ratios of 5.3 and 7.8 in the F336W and F656N bands, respectively. This is the first time that an exoplanet has been directly imaged in the UV. Our observed Hα flux of PDS 70 b is higher by 3.5σ than the most recent published result. However, the light curve retrieved from our observations does not support greater than 30% variability in the planet’s Hα emission in six epochs over a five month timescale. We estimate a mass accretion rate of 1.4±0.2×10{sup −8}M{sub Jup} yr{sup −1}. Hα accounts for 36% of the total accretion luminosity. Such a high proportion of energy released in line emission suggests efficient production of Hα emission in planetary accretion, and motivates using the Hα band for searches of accreting planets. These results demonstrate HST/WFC3/UVIS’s excellent high-contrast imaging performance and highlight its potential for planet formation studies.