Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “accretion”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

River bank erosion and accretion rates, and planform metrics of the East River, downstream of Gothic Colorado over intervals between 1955 and 2015

This dataset provides the tabular summary of analysis of an alluvial floodplain reach of the East River, downstream of Gothic, CO near Crested Butte. The measurements include bank erosion and accretion rates, channel width, bank curvature, and the aspect/orientation of the river banks. The measurements were derived from binary masks of the location of the river channel from aerial and satellite imagery collected on seven dates over a sixty-year timespan: 1955, 1973, 1983, 1990, 2001, 2011, and 2015. The masks were analyzed using the Spatially Continuous Riverbank Erosion and Accretion Measurements (SCREAM) software detailed in Rowland et al. 2016. The masks used in this analysis can be found an accompanying dataset (DOI: ). Rates of change along the East River were measured over a total of 8-time intervals: 1955-1973; 1955-2015; 1973-1983; 1983-1990; 1990-2001; 2001-2011; 2001-2015; and 2011-2015.In files with “summary” in the name, the data is provided at a pixel level, where each mapped bank pixel has an associated erosion or accretion value, a channel width, a curvature value, and an aspect each river and time period will have an individual file. Files with “Segments” in name provide data that is averaged along segments of the rivers. These segments are approximately 10 channel widths in length. In addition to erosion and accretion rates, the segment-based results include area measurements of erosion and accretion, islands, and channels. The number of islands is also included.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

The criterion of supernova explosion revisited: The mass accretion history

By performing neutrino-radiation hydrodynamic simulations in spherical symmetry (1D) and axial symmetry (2D) with different progenitor models by Woosley and Heger from 12 to 100 M{sub ⊙}, we find that all 1D runs fail to produce an explosion and several 2D runs succeed. The difference in the shock evolutions for different progenitors can be interpreted by the difference in their mass accretion histories, which are in turn determined by the density structures of progenitors. The mass accretion history has two phases in the majority of the models: the earlier phase, in which the mass accretion rate is high and rapidly decreasing, and the later phase, with a low and almost constant accretion rate. They are separated by the so-called turning point, the origin of which is a change of the accreting layer. We argue that shock revival will most likely occur around the turning point and hence that its location in the M-dot --L{sub ν} plane will be a good measure for the possibility of shock revival: if the turning point lies above the critical curve and the system stays there for a long time, shock revival will obtain. In addition, we develop a phenomenological model to approximately evaluate the trajectories in the M-dot --L{sub ν} plane, which, after calibrating free parameters by a small number of 1D simulations, reproduces the location of the turning point reasonably well by using the initial density structure of progenitor alone. We suggest the application of the phenomenological model to a large collection of progenitors in order to infer without simulations which ones are more likely to explode.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Upper limit of the viscosity parameter in accretion flows around a black hole with shock waves

Black hole accretion is necessarily transonic; thus, flows must become supersonic and, therefore, sub-Keplerian before they enter into the black hole. The viscous timescale is much longer than the infall timescale close to a black hole. Hence, the angular momentum remains almost constant and the centrifugal force ∼l{sup 2}/r{sup 3} becomes increasingly dominant over the gravitational force ∼1/r{sup 2}. The slowed down matter piles creating an accretion shock. The flow between shock and inner sonic point is puffed up and behaves like a boundary layer. This so-called Comptonizing cloud/corona produces hard X-rays and jets/outflows and, therefore, is an important component of black hole astrophysics. In this paper, we study steady state viscous, axisymmetric, transonic accretion flows around a Schwarzschild black hole. We adopt a viscosity parameter α and compute the highest possible value of α (namely, α{sub cr}) for each pair of two inner boundary parameters (namely, specific angular momentum carried to horizon, l{sub in} and specific energy at inner sonic point, E(x{sub in})) which is still capable of producing a standing or oscillating shock. We find that while such possibilities exist for α as high as α{sub cr}=0.3 in very small regions of the flow parameter space, typical α{sub cr} appears to be about ∼0.05–0.1. Coincidentally, this also happens to be the typical viscosity parameter achieved by simulations of magnetorotational instabilities in accretion flows. We therefore believe that all realistic accretion flows are likely to have centrifugal pressure supported shocks unless the viscosity parameter everywhere is higher than α{sub cr}.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

An accretion disk-outflow model for hysteretic state transition in X-ray binaries

We suggest a model of the advection-dominated accretion flow (ADAF) with magnetically driven outflows to explain the hysteretic state transition observed in X-ray binaries (XRBs). The transition from a thin disk to an ADAF occurs when the mass accretion rate is below a critical value. The critical mass accretion rate for the ADAF can be estimated by equating the equilibration timescale to the accretion timescale of the ADAF, which is sensitive to its radial velocity. The radial velocity of thin disks is very small, which leads to the advection of the external field in thin disks becoming very inefficient. ADAFs are present in the low/hard states of XRBs, and their radial velocity is large compared with the thin disk. The external field can be dragged inward efficiently by the ADAF, so a strong large-scale magnetic field threading the ADAF can be formed, which may accelerate a fraction of gas in the ADAF into the outflows. Such outflows may carry away a large amount of angular momentum from the ADAF, which significantly increases the radial velocity of the ADAF. This leads to a high critical mass accretion rate, below which an ADAF with magnetic outflows can survive. Our calculations show that the critical luminosity of the ADAF with magnetic outflows can be one order of magnitude higher than that for a conventional ADAF, if the ratio of gas to magnetic pressure β∼4 in the disk. This can naturally explain the hysteretic state transition observed in XRBs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

ELECTROMAGNETIC SPINDOWN OF A TRANSIENT ACCRETING MILLISECOND PULSAR DURING QUIESCENCE

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.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

A HOT AND MASSIVE ACCRETION DISK AROUND THE HIGH-MASS PROTOSTAR IRAS 20126+4104

We present new spectral line observations of the CH{sub 3}CN molecule in the accretion disk around the massive protostar IRAS 20126+4104 with the Submillimeter Array, which, for the first time, measure the disk density, temperature, and rotational velocity with sufficient resolution (0.″37, equivalent to ∼600 au) to assess the gravitational stability of the disk through the Toomre- Q parameter. Our observations resolve the central 2000 au region that shows steeper velocity gradients with increasing upper state energy, indicating an increase in the rotational velocity of the hotter gas nearer the star. Such spin-up motions are characteristics of an accretion flow in a rotationally supported disk. We compare the observed data with synthetic image cubes produced by three-dimensional radiative transfer models describing a thin flared disk in Keplerian motion enveloped within the centrifugal radius of an angular-momentum-conserving accretion flow. Given a luminosity of 1.3 × 10{sup 4} L {sub ⊙}, the optimized model gives a disk mass of 1.5 M {sub ⊙} and a radius of 858 au rotating about a 12.0 M {sub ⊙} protostar with a disk mass accretion rate of 3.9 × 10{sup −5} M {sub ⊙} yr{sup −1}. Our study finds that, in contrast to some theoretical expectations, the disk is hot and stable to fragmentation with Q > 2.8 at all radii which permits a smooth accretion flow. These results put forward the first constraints on gravitational instabilities in massive protostellar disks, which are closely connected to the formation of companion stars and planetary systems by fragmentation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Global MHD simulations of accretion disks in cataclysmic variables. I. The importance of spiral shocks

We present results from the first global 3D MHD simulations of accretion disks in cataclysmic variable (CV) systems in order to investigate the relative importance of angular momentum transport via turbulence driven by the magnetorotational instability (MRI) compared with that driven by spiral shock waves. Remarkably, we find that even with vigorous MRI turbulence, spiral shocks are an important component of the overall angular momentum budget, at least when temperatures in the disk are high (so that Mach numbers are low). In order to understand the excitation, propagation, and damping of spiral density waves in our simulations more carefully, we perform a series of 2D global hydrodynamical simulations with various equation of states, both with and without mass inflow via the Lagrangian point (L1). Compared with previous similar studies, we find the following new results. (1) The linear wave dispersion relation fits the pitch angles of spiral density waves very well. (2) We demonstrate explicitly that mass accretion is driven by the deposition of negative angular momentum carried by the waves when they dissipate in shocks. (3) Using Reynolds stress scaled by gas pressure to represent the effective angular momentum transport rate α{sub eff} is not accurate when mass accretion is driven by non-axisymmetric shocks. (4) Using the mass accretion rate measured in our simulations to directly measure α defined in standard thin-disk theory, we find 0.02≲α{sub eff}≲0.05 for CV disks, consistent with observed values in quiescent states of dwarf novae. In this regime, the disk may be too cool and neutral for the MRI to operate and spiral shocks are a possible accretion mechanism. However, we caution that our simulations use unrealistically low Mach numbers in this regime and, therefore, future models with more realistic thermodynamics and non-ideal MHD are warranted.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

ACCRETION DISK SIGNATURES IN TYPE I X-RAY BURSTS: PROSPECTS FOR FUTURE MISSIONS

Type I X-ray bursts and superbursts from accreting neutron stars illuminate the accretion disk and produce a reflection signal that evolves as the burst fades. Examining the evolution of reflection features in the spectra will provide insight into the burst–disk interaction, a potentially powerful probe of accretion disk physics. At present, reflection has been observed during only two bursts of exceptional duration. We investigate the detectability of reflection signatures with four of the latest well-studied X-ray observatory concepts: Hitomi , Neutron Star Interior Composition Explorer ( NICER ), Athena , and Large Observatory For X-ray Timing ( LOFT ). Burst spectra are modeled for different values for the flux, temperature, and the disk ionization parameter, which are representative for most known bursts and sources. The effective area and throughput of a Hitomi -like telescope are insufficient for characterizing burst reflection features. NICER and Athena will detect reflection signatures in Type I bursts with peak fluxes ≳10{sup 7.5} erg cm{sup 2} s{sup 1} and also effectively constrain the reflection parameters for bright bursts with fluxes of ∼10{sup 7} erg cm{sup 2} s{sup 1} in exposures of several seconds. Thus, these observatories will provide crucial new insight into the interaction of accretion flows and X-ray bursts. For sources with low line-of-sight absorption, the wide bandpass of these instruments allows for the detection of soft X-ray reflection features, which are sensitive to the disk metallicity and density. The large collecting area that is part of the LOFT design would revolutionize the field by tracing the evolution of the accretion geometry in detail throughout short bursts.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

SOFT X-RAY EXCESS FROM SHOCKED ACCRETING PLASMA IN ACTIVE GALACTIC NUCLEI

We propose a novel theoretical model to describe the physical identity of the soft X-ray excess that is ubiquitously detected in many Seyfert galaxies, by considering a steady-state, axisymmetric plasma accretion within the innermost stable circular orbit around a black hole (BH) accretion disk. We extend our earlier theoretical investigations on general relativistic magnetohydrodynamic accretion, which implied that the accreting plasma can develop into a standing shock under suitable physical conditions, causing the downstream flow to be sufficiently hot due to shock compression. We perform numerical calculations to examine, for sets of fiducial plasma parameters, the physical nature of fast magnetohydrodynamic shocks under strong gravity for different BH spins. We show that thermal seed photons from the standard accretion disk can be effectively Compton up-scattered by the energized sub-relativistic electrons in the hot downstream plasma to produce the soft excess feature in X-rays. As a case study, we construct a three-parameter Comptonization model of inclination angle θ {sub obs}, disk photon temperature kT {sub in}, and downstream electron energy kT {sub e}to calculate the predicted spectra in comparison with a 60 ks XMM-Newton /EPIC-pn spectrum of a typical radio-quiet Seyfert 1 active galactic nucleus, Ark 120. Our χ {sup 2}-analyses demonstrate that the model is plausible for successfully describing data for both non-spinning and spinning BHs with derived ranges of 61.3 keV ≲ kT {sub e}≲ 144.3 keV, 21.6 eV ≲ kT {sub in} ≲ 34.0 eV, and 17.°5 ≲ θ {sub obs} ≲ 42.°6, indicating a compact Comptonizing region of three to four gravitational radii that resembles the putative X-ray coronae.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

An accretion-jet model for M87: interpreting the spectral energy distribution and faraday rotation measure

M87 is arguably the best supermassive black hole (BH) to explore jet and/or accretion physics, due to its proximity and fruitful high-resolution multi-waveband observations. We model the multi-wavelength spectral energy distribution (SED) of the M87 core that observed at a scale of 0.4 arcsec (∼10{sup 5} R {sub g}, R {sub g} is gravitational radius), as recently presented by Prieto et al. Similar to Sgr A*, we find that the millimeter bump as observed by the Atacama Large Millimeter/submillimeter Array can be modeled by the synchrotron emission of the thermal electrons in an advection-dominated accretion flow (ADAF), while the low-frequency radio emission and X-ray emission may predominantly come from the jet. The millimeter radiation from ADAF predominantly comes from the region within 10R {sub g}, which is roughly consistent with the recent very long baseline interferometry observations at 230 GHz. We further calculate the Faraday rotation measure (RM) from both ADAF and jet models, and find that the RM predicted from the ADAF is roughly consistent with the measured value, while the RM predicted from the jet is much higher if jet velocity close to the BH is low or moderate (e.g., v {sub jet} ≲ 0.6 c). With the constraints from the SED modeling and RM, we find that the accretion rate close to the BH horizon is ∼(0.2--1) ×10{sup −3} M{sub ⊙} yr{sup −1}≪ M-dot {sub B}∼0.2 M{sub ⊙} yr{sup −1} ( M-dot {sub B} is Bondi accretion rate), where the electron density profile, n {sub e} ∝ r {sup ∼−1}, in the accretion flow, is consistent with that determined from X-ray observation inside the Bondi radius and recent numerical simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

SHOCK-DRIVEN ACCRETION IN CIRCUMPLANETARY DISKS: OBSERVABLES AND SATELLITE FORMATION

Circumplanetary disks (CPDs) control the growth of planets, supply material for satellites to form, and provide observational signatures of young forming planets. We have carried out two-dimensional hydrodynamical simulations with radiative cooling to study CPDs and suggested a new mechanism to drive the disk accretion. Two spiral shocks are present in CPDs, excited by the central star. We find that spiral shocks can at least contribute to, if not dominate, the angular momentum transport and energy dissipation in CPDs. Meanwhile, dissipation and heating by spiral shocks have a positive feedback on shock-driven accretion itself. As the disk is heated up by spiral shocks, the shocks become more open, leading to more efficient angular momentum transport. This shock-driven accretion is, on the other hand, unsteady due to production and destruction of vortices in disks. After being averaged over time, a quasi-steady accretion is reached from the planet’s Hill radius all the way to the planet surface, and the disk α coefficient characterizing angular momentum transport is ∼0.001–0.02. The disk surface density ranges from 10 to 1000 g cm{sup −2} in our simulations, which is at least three orders of magnitude smaller than the “minimum-mass subnebula” model used to study satellite formation; instead it is more consistent with the “gas-starved” satellite formation model. Finally, we calculate the millimeter flux emitted by CPDs at ALMA and EVLA wavelength bands and predict the flux for several recently discovered CPD candidates, which suggests that ALMA is capable of discovering these accreting CPDs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Accretion Properties of PDS 70b with MUSE

We report a new evaluation of the accretion properties of PDS 70b obtained with the Very Large Telescope/Multi Unit Spectroscopic Explorer. The main difference from the previous studies of Haffert et al. and Aoyama & Ikoma is in the mass accretion rate. Simultaneous multiple line observations, such as Hα and Hβ, can better constrain the physical properties of an accreting planet. While we clearly detected Hα emissions from PDS 70b, no Hβ emissions were detected. We estimate the line flux of Hβ with a 3σ upper limit to be 2.3 × 10{sup −16} erg s{sup −1} cm{sup −2}. The flux ratio F {sub Hβ}/F {sub Hα} for PDS 70b is <0.28. Numerical investigations by Aoyama et al. suggest that F {sub Hβ}/F {sub Hα} should be close to unity if the extinction is negligible. We attribute the reduction of the flux ratio to the extinction, and estimate the extinction of Hα (A {sub Hα}) for PDS 70b to be >2.0 mag using the interstellar extinction value. By combining with the Hα linewidth and the dereddening line luminosity of Hα, we derive the PDS 70b mass accretion rate to be ≳5 × 10{sup −7} M {sub Jup} yr{sup −1}. The PDS 70b mass accretion rate is an order of magnitude larger than that of PDS 70. We found that the filling factor f {sub f} (the fractional area of the planetary surface emitting Hα) is ≳0.01, which is similar to the typical stellar value. The small value of f {sub f} indicates that the Hα emitting areas are localized at the surface of PDS 70b.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Keck/NIRC2 L’-band Imaging of Jovian-mass Accreting Protoplanets around PDS 70

We present L’-band imaging of the PDS 70 planetary system with Keck/NIRC2 using the new infrared pyramid wave front sensor. We detected both PDS 70 b and c in our images, as well as the front rim of the circumstellar disk. After subtracting off a model of the disk, we measured the astrometry and photometry of both planets. Placing priors based on the dynamics of the system, we estimated PDS 70 b to have a semimajor axis of 20{sub −4}{sup +3} au and PDS 70 c to have a semimajor axis of 34{sub −6}{sup +12} au (95% credible interval). We fit the spectral energy distribution (SED) of both planets. For PDS 70 b, we were able to place better constraints on the red half of its SED than previous studies and inferred the radius of the photosphere to be 2–3 R {sub Jup}. The SED of PDS 70 c is less well constrained, with a range of total luminosities spanning an order of magnitude. With our inferred radii and luminosities, we used evolutionary models of accreting protoplanets to derive a mass of PDS 70 b between 2 and 4 M {sub Jup} and a mean mass accretion rate between 3 × 10{sup −7} and 8 × 10{sup −7} M {sub Jup}/yr. For PDS 70 c, we computed a mass between 1 and 3 M {sub Jup} and mean mass accretion rate between 1 × 10{sup −7} and 5 × 10{sup −7} M {sub Jup}/yr. The mass accretion rates imply dust accretion timescales short enough to hide strong molecular absorption features in both planets’ SEDs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

HIGH ENERGY NEUTRINOS PRODUCED IN THE ACCRETION DISKS BY NEUTRONS FROM NUCLEI DISINTEGRATED IN THE AGN JETS

We investigate the consequences of acceleration of nuclei in jets of active galaxies not far from the surface of an accretion disk. The nuclei can be accelerated in the re-connection regions in the jet and/or at the jet boundary, between the relativistic jet and its cocoon. It is shown that the relativistic nuclei can efficiently fragment onto specific nucleons in collisions with the disk radiation. Neutrons, directed toward the accretion disk, take a significant part of energy from the relativistic nuclei. These neutrons develop a cascade in the dense accretion disk. We calculate the neutrino spectra produced in such a hadronic cascade within the accretion disk. We propose that the neutrinos produced in such a scenario, from the whole population of super-massive black holes in active galaxies, can explain the extragalactic neutrino background recently measured by the IceCube neutrino detector, provided that a 5% fraction of galaxies have an active galactic nucleus and a few percent of neutrons reach the accretion disk. We predict that the neutrino signals in the present neutrino detectors, produced in terms of such a model, will not be detectable even from the nearby radio galaxies similar to M87.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

rHARM: ACCRETION AND EJECTION IN RESISTIVE GR-MHD

Turbulent magnetic diffusivity plays an important role for accretion disks and the launching of disk winds. We have implemented magnetic diffusivity and respective resistivity in the general relativistic MHD code HARM. This paper describes the theoretical background of our implementation, its numerical realization, our numerical tests, and preliminary applications. The test simulations of the new code rHARM are compared to an analytic solution of the diffusion equation and a classical shock tube problem. We have further investigated the evolution of the magnetorotational instability (MRI) in tori around black holes (BHs) for a range of magnetic diffusivities. We find an indication for a critical magnetic diffusivity (for our setup) beyond which no MRI develops in the linear regime and for which accretion of torus material to the BH is delayed. Preliminary simulations of magnetically diffusive thin accretion disks around Schwarzschild BHs that are threaded by a large-scale poloidal magnetic field show the launching of disk winds with mass fluxes of about 50% of the accretion rate. The disk magnetic diffusivity allows for efficient disk accretion that replenishes the mass reservoir of the inner disk area and thus allows for long-term simulations of wind launching for more than 5000 time units.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Vertical Advection Effects on Hyper-accretion Disks and Potential Link between Gamma-Ray Bursts and Kilonovae

Recent simulations on super-Eddington accretion flows have shown that, apart from the diffusion process, the vertical advection based on magnetic buoyancy can be a more efficient process to release the trapped photons in the optically thick disk. As a consequence, the radiative luminosity from the accretion disk can be far beyond the Eddington value. Following this spirit, we revisit the structure and radiation of hyper-accretion disks with mass accretion rates in the range of 10{sup −3}∼10 M{sub ⊙} s{sup −1}. Our results show that, due to the strong cooling through the vertical advection, the disk temperature becomes lower than that in the classic model without the vertical advection process, and therefore the neutrino luminosity from the disk is lower. On the other hand, the gamma-ray photons released through the vertical advection can be extremely super-Eddington. We argue that the large amount of escaped gamma-ray photons may have more significant contribution to the primordial fireball than the neutrino annihilation, and may hint at a link between gamma-ray bursts and kilonovae in the black hole hyper-accretion scenario.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Accretion in Radiative Equipartition (AiRE) Disks

Standard accretion disk theory predicts that the total pressure in disks at typical (sub-)Eddington accretion rates becomes radiation pressure dominated. However, radiation pressure dominated disks are thermally unstable. Since these disks are observed in approximate steady state over the instability timescale, our accretion models in the radiation-pressure-dominated regime (i.e., inner disk) need to be modified. Here, we present a modification to the Shakura and Sunyaev model, where the radiation pressure is in equipartition with the gas pressure in the inner region. We call these flows accretion in radiative equipartition (AiRE) disks. We introduce the basic features of AiRE disks and show how they modify disk properties such as the Toomre parameter and the central temperature. We then show that the accretion rate of AiRE disks is limited from above and below, by Toomre and nodal sonic point instabilities, respectively. The former leads to a strict upper limit on the mass of supermassive black holes as a function of cosmic time (and spin), while the latter could explain the transition between hard and soft states of X-ray binaries.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Constraints from Dust Mass and Mass Accretion Rate Measurements on Angular Momentum Transport in Protoplanetary Disks

In this paper, we investigate the relation between disk mass and mass accretion rate to constrain the mechanism of angular momentum transport in protoplanetary disks. We find a correlation between dust disk mass and mass accretion rate in Chamaeleon I with a slope that is close to linear, similar to the one recently identified in Lupus. We investigate the effect of stellar mass and find that the intrinsic scatter around the best-fit M{sub dust}–M{sub ⋆} and M-dot {sub acc}–M{sub ⋆} relations is uncorrelated. We simulate synthetic observations of an ensemble of evolving disks using a Monte Carlo approach and find that disks with a constant α viscosity can fit the observed relations between dust mass, mass accretion rate, and stellar mass but overpredict the strength of the correlation between disk mass and mass accretion rate when using standard initial conditions. We find two possible solutions. In the first one, the observed scatter in M{sub dust} and M-dot {sub acc} is not primordial, but arises from additional physical processes or uncertainties in estimating the disk gas mass. Most likely grain growth and radial drift affect the observable dust mass, while variability on large timescales affects the mass accretion rates. In the second scenario, the observed scatter is primordial, but disks have not evolved substantially at the age of Lupus and Chamaeleon I owing to a low viscosity or a large initial disk radius. More accurate estimates of the disk mass and gas disk sizes in a large sample of protoplanetary disks, through either direct observations of the gas or spatially resolved multiwavelength observations of the dust with ALMA, are needed to discriminate between both scenarios or to constrain alternative angular momentum transport mechanisms such as MHD disk winds.

79 ASTRONOMY AND ASTROPHYSICS↗