Disk Photosphere Research Articles

A Far-ultraviolet-detected Accretion Shock at the Star–Disk Boundary of FU Ori

FU Ori objects are the most extreme eruptive young stars known. Their 4–5 mag photometric outbursts last for decades and are attributed to a factor of up to 10,000 increase in the stellar accretion rate. The nature of the accretion disk-to-star interface in FU Ori objects has remained a mystery for decades. To date, attempts to directly observe a shock or boundary layer have been thwarted by the apparent lack of emission in excess of the accretion disk photosphere down to λ = 2300 Å. We present a new near-ultraviolet and the first high-sensitivity far-ultraviolet (FUV) spectrum of FU Ori. The FUV continuum is detected for the first time and, at λ = 1400 Å, is more than 104 times brighter than predicted by a viscous accretion disk. We interpret the excess as arising from a shock at the boundary between the disk and the stellar surface. We model the shock emission as a blackbody and find that the temperature of the shocked material is T FUV ≈ 16,000 ± 2000 K. The shock temperature corresponds to an accretion flow along the surface of the disk that reaches a velocity of 40 km s−1 at the boundary, consistent with predictions from simulations.

Common origin for black holes in both high mass X-ray binaries and gravitational-wave sources

Black-hole (BH) high-mass X-ray binary (HMXB) systems are likely to be the progenitors of BH-BH mergers detected in gravitational waves by LIGO/Virgo/KAGRA (LVK). Yet merging BHs reach higher masses ($ than BHs in HMXBs ($ and typically exhibit lower spins BH 0.25$ with a larger values tail) than what is often claimed for BHs in HMXBs BH 0.9$). This could suggest that these two classes of systems belong to different populations, but here we show that this may not necessarily be the case. The difference in masses is easily explained as the known HMXB-BHs are in galaxies with relatively high metallicity, so their progenitor stars are subject to strong mass loss from winds, leading to relatively low-mass BH at core collapse. Conversely, LVK is also able to detect BHs from low-metallicity galaxies that are known to naturally produce more massive stellar-origin BHs. However, the difference in spin is more difficult to explain. Models with efficient angular momentum transport in stellar interiors produce slowly spinning progenitors for both LVK and HMXB BHs. Known HMXBs have orbital periods that are too long for efficient tidal spin-up and are also unlikely to have undergone significant accretion spin-up. Instead, we show that the derived value of the BH spin depends strongly on how the HMXB accretion disc emission is modelled. We argue that since Cyg X-1 is never observed to be in a soft spectral state, the appropriate spectral models must take into account the Comptonisation of the disc photosphere. We show that such models are consistent with low spin values, namely: $a_ BH This was recently confirmed by other teams for both Cyg X-1 and LMC X-1 and here we show this is also the case for M33 X-7. We conclude that all known HMXB BHs can exhibit a low spin in accordance with the results of stellar evolution models. Hence, the observations presented in this work are consistent with the scenario where LVK BHs and HMXB BHs belong to the same population.

Turbulent processing of PAHs in protoplanetary discs

Context. Polycyclic aromatic hydrocarbons (PAHs) have been detected in numerous circumstellar discs. Despite the correlation between stellar temperature and low PAH detections rates, the diversity of PAH detections and non-detections at similar stellar properties is not well understood. Aims. We propose the continuous processing of PAHs through clustering, adsorption on dust grains, and their reverse-processes as key mechanisms to reduce the emission-capable PAH abundance in protoplanetary discs. This cycle of processing is driven by vertical turbulence in the disc mixing PAHs between the disc midplane and the photosphere. Methods. We used a theoretical Monte Carlo model for photodesorption in the photosphere and a coagulation code in the disc midplane to estimate the relevance and timescale of these processes in a Herbig Ae/Be disc environment. By combining these components in a 1D vertical model, we calculated the gas-phase depletion of PAHs that stick as clusters on dust grains. Results. Our results show that the clustering of gas-phase PAHs is very efficient, and that clusters with more than 100 monomers can grow for years before they are able to freeze out in the disc midplane. Once a PAH cluster is frozen on the dust grain surface, the large heat capacity of these clusters prevents them from evaporating off the grains in UV-rich environments such as the photosphere. Therefore, the clustering of PAHs followed by freeze-out can lead to a depletion of gas-phase PAHs in protoplanetary discs. We find that this mechanism is more efficient when the PAH species has fewer carbon atoms. In contrast, PAH monomers and very small clusters consisting of a few monomers can easily detach from the grain by absorption of a single UV photon. Evaluated over the lifetime of protoplanetary discs, we find a depletion of PAHs by a factor that ranges between 50 and 1000 compared to the standard ISM abundance of PAHs in the inner disc through turbulent processing. Conclusions. Through these processes, we favour PAHs smaller than circumovalene (C66H20) as the major gas-phase emitters of the disc photosphere as larger PAH monomers cannot photodesorb from the grain surface. These gas-phase PAHs co-exist with large PAH clusters sticking on dust grains. We find a close relation between the amount of PAHs frozen out on dust grains and the dust population, as well as the strength of the vertical turbulence.

Surveying the inner structure of massive young stellar objects using L-band spectroscopy

We present results from a high spectral resolution (6 km s-1) survey of five massive protostars in the wavelength range of 2.95 and 3.25 µm, conducted with iSHELL at the InfraRed Telescope Facility (IRTF). Our targets are Mon R2 IRS 2, Mon R2 IRS 3, AFGL 2136, Orion BN and S140 IRS 1. Two of our five targets (Mon R2 IRS 3 and AFGL 2136) show transitions from organic species, with MonR2 IRS 3 showing HCN lines in emission, and AFGL 2136 showing HCN and C2H2 lines in absorption. The velocity of the emission lines of HCN of MonR2 IRS 3A are consistent with CO emission features in lines up to J = 26, as both are red-shifted with respect to the systemic velocity. Carbon monoxide lines also show blue-shifted absorption. This P-Cygni line profile, commonly observed towards massive young stellar objects, is likely due to an expanding shell, which is supported by sub-millimetre velocity maps of HCN. Alternatively HCN emission may arise from the upper layers of a disk photosphere, as has been suggested for the massive protostar AFGL 2591. Absorption lines in AFGL 2136 may either originate in foreground cloud or in the disk photosphere. For a foreground cloud, the data require that the foreground gas only covers the source partially (0.3) at 13 µm. In contrast, absorption lines at 3 and 7 µm require a covering factor of >0.9. Analysing the 13 µm HCN absorption lines in terms of absorption by gas in the photosphere of a disk, results in physical conditions that are consistent over all three vibrational modes. C2H2 absorption lines reveal an increasing temperature and abundance with decreasing wavelength, indicative of a radial abundance gradient. We conclude that the disk model is the best interpretation of the absorption lines of AFGL 2136.

Electromagnetic Signatures of Relativistic Explosions in the Disks of Active Galactic Nuclei

Abstract The disks of active galactic nuclei (AGNs), traditionally studied as feeders of the supermassive black holes (SMBHs) at their centers, are also hosts to massive stars and hence their neutron star (NS) and black hole (BH) remnants. Migration traps and gas torques in these disks favor binary formation and enhance the rate of compact object mergers. In these environments both long gamma-ray bursts (GRBs) from the death of massive stars and short GRBs from NS–NS to NS–BH mergers are expected. However, their properties in the environment of AGN disks have never been studied. Here we show that GRBs in AGN disks can display unique features, owing to the unusual relative position of the shocks that characterize the burst evolution and the Thomson photosphere of the AGN disk. In dense environments, for example, a relativistic reverse shock develops early, likely powering the prompt emission instead of internal shocks. The transient’s time evolution is also compressed, yielding afterglow emission that is brighter and may peak earlier than for GRBs in the interstellar medium. Additionally, in regions of the disk that are sufficiently dense and extended, the light curves are dominated by diffusion, since the fireball remains inside the disk photosphere throughout the entire evolution. These diffusion-dominated transients emerge on timescales of days in disks around SMBHs of ∼ 106 M ⊙ to years for SMBHs of ∼ 108 M ⊙. Finally, a large fraction of events, especially in AGNs with SMBHs ≲ 107 M ⊙, display time-variable absorption in the X-ray band.

Disc light variability in the FUor star V646 Puppis as observed by TESS and from the ground

Context. We investigate small-scale light variations in V646 Pup occurring on timescales of days, weeks, and years. Aims. We aim to investigate whether this variability is similar to that observed in FU Ori. Methods. We observed V646 Pup on six occasions at the SAAO and CTIO between 2013 and 2018 with Johnson and Sloan filters, typically using a one-day cadence maintained for two to four weeks. We also utilised the public-domain 1512-day-long ASAS-SN light curve and TESS photometry obtained in 2019 over 24.1 days with a 30 min cadence. New SAAO low-resolution spectra assist in updating major disc parameters, while the archival high-resolution Keck spectra are used to search for temporal changes in the disc rotational profiles. Results. The ground-based observations confirm the constantly decreasing brightness of V646 Pup at the rate of 0.018 mag yr−1. Precise i-band sensitive TESS data show that the slight, 0.005−0.01 mag, light variations imposed on this general trend do consist of a few independent wave trains of an apparently time-coherent nature. Assuming that this is typical situation, based on an analysis of colour–magnitude diagrams obtained for earlier epochs, we were able to make a preliminarily inference that the bulk of the light changes observed could be due to the rotation of disc photosphere inhomogeneities, arising between 10–12 R⊙ from the star. We do not exclude the possibility that these inhomogeneities could also manifest themselves in the rotational profiles of the disc, as obtained from the high-resolution spectra. Assuming Keplerian rotation of these inhomogeneities, we give a preliminary determination of the stellar mass at 0.7–0.9 M⊙. Conclusions. Over certain weeks, at least, V646 Pup has shown time-coherent light variability pattern(s) that could be explained by the rotation of an inhomogeneous disc photosphere. These preliminary results are similar to those better established for FU Ori, which suggests a common driving mechanism(s).

High-resolution Infrared Spectroscopy of Hot Molecular Gas in AFGL 2591 and AFGL 2136: Accretion in the Inner Regions of Disks around Massive Young Stellar Objects

Abstract We have performed a high-resolution 4–13 μm spectral survey of the hot molecular gas associated with the massive protostars AFGL 2591 and AFGL 2136, utilizing the Echelon Cross Echelle Spectrograph (EXES) on board the Stratospheric Observatory for Infrared Astronomy, and the iSHELL instrument and Texas EchelonCross Echelle Spectrograph (TEXES) on the NASA Infrared Telescope Facility (IRTF). Here we present the results of this survey with analysis of CO, HCN, C2H2, NH3, and CS, deriving the physical conditions for each species. Also from the IRTF, iSHELL data at 3 μm for AFGL 2591 are presented that show HCN and C2H2 in emission. In the EXES and TEXES data, all species are detected in absorption, and temperatures and abundances are found to be high (600 K and 10−6, respectively). Differences of up to an order of magnitude in the abundances of transitions that trace the same ground-state level are measured for HCN and C2H2. The mid-infrared continuum is known to originate in a disk, hence we attribute the infrared absorption to arise in the photosphere of the disk. As absorption lines require an outwardly decreasing temperature gradient, we conclude that the disk is heated in the midplane by viscous heating due to accretion. We attribute the near-IR emission lines to scattering by molecules in the upper layers of the disk photosphere. The absorption lines trace the disk properties at 50 au where high-temperature gas-phase chemistry is taking place. Abundances are consistent with chemical models of the inner disk of Herbig disks.

Global 3D radiation magnetohydrodynamic simulations for FU Ori’s accretion disc and observational signatures of magnetic fields

ABSTRACT FU Ori is the prototype of FU Orionis systems that are outbursting protoplanetary discs. Magnetic fields in FU Ori’s accretion discs have previously been detected using spectropolarimetry observations for Zeeman effects. We carry out global radiation ideal MHD simulations to study FU Ori’s inner accretion disc. We find that (1) when the disc is threaded by vertical magnetic fields, most accretion occurs in the magnetically dominated atmosphere at z ∼ R, similar to the ‘surface accretion’ mechanism in previous locally isothermal MHD simulations. (2) A moderate disc wind is launched in the vertical field simulations with a terminal speed of ∼300–500 km s−1 and a mass-loss rate of 1–10 per cent the disc accretion rate, which is consistent with observations. Disc wind fails to be launched in simulations with net toroidal magnetic fields. (3) The disc photosphere at the unit optical depth can be either in the wind launching region or the accreting surface region. Magnetic fields have drastically different directions and magnitudes between these two regions. Our fiducial model agrees with previous optical Zeeman observations regarding both the field directions and magnitudes. On the other hand, simulations indicate that future Zeeman observations at near-IR wavelengths or towards other FU Orionis systems may reveal very different magnetic field structures. (4) Due to energy loss by the disc wind, the disc photosphere temperature is lower than that predicted by the thin disc theory, and the previously inferred disc accretion rate may be lower than the real accretion rate by a factor of ∼2–3.

Medium-resolution échelle spectroscopy of the Red Square Nebula, MWC 922

Context. Medium-resolution echelle spectra of the Red Square Nebula surrounding the star MWC 922 are presented. The spectra have been obtained in 2010 and 2012 using the X-shooter spectrograph mounted on the Very Large Telescope (VLT) in Paranal, Chile. The spectrum covers a wavelength range between 300 nm−2.5 μ m and shows that the nebula is rich in emission lines. Aims. We aim to identify the emission lines and use them as a tool to determine the physical and chemical characteristics of the nebula. The emission lines are also used to put constraints on the structure of the nebula and on the nature of the central stars.Methods. We analyzed and identified emission lines that indicated that the Red Square Nebula consists of a low density bipolar outflow, eminent in the broad emission component seen in [Fe ii], as well as in P Cygni line profiles indicative of fast outflowing material. The narrow component in the [Fe ii] lines is most likely formed in the photosphere of a surrounding disk. Some of the emission lines show a pronounced double peaked profile, such as Ca ii, indicating an accretion disk in Keplerian rotation around the central star. [O i] emission lines are formed in the neutral atomic zone separating the ionized disk photosphere from the molecular gas in the interior of the disk, which is prominent in molecular CO emission in the near-IR. [N ii] and [S ii] emission clearly originates in a low density but fairly hot (7 000−10 000 K) nebular environment. H i recombination lines trace the extended nebula as well as the photosphere of the disk.Results. These findings put constraints on the evolution of the central objects in MWC 922. The Red Square shows strong similarities to the Red Rectangle Nebula, both in morphology and in its mid-IR spectroscopic characteristics. As for the Red Rectangle, the observed morphology of the nebula reflects mass-loss in a binary system. Specifically, we attribute the biconical morphology and the associated rung-like structure to the action of intermittent jets blown by the accreting companion in a dense shell, which has been created by the primary. We stress, though, that despite the morphological similarities, these two objects represent very different classes of stellar objects.

Investigating the mass of the intermediate mass black hole candidate HLX-1 with the slimbh model

In this paper we present a comprehensive study of the mass of the intermediate mass black hole candidate HLX-1 in the galaxy ESO 243-49. We analyse the continuum X-ray spectra collected by Swift, XMM-Newton, and Chandra with the slim disc model, SLIMBH, and estimate the black hole mass for the full range of inclination (inc = 0{\deg} - 85{\deg}) and spin (a* = 0 - 0.998). The relativistic SLIMBH model is particularly suited to study high luminosity disc spectra as it incorporates the effects of advection, such as the shift of the inner disc edge towards smaller radii and the increasing height of the disc photosphere (including relativistic ray-tracing from its proper location rather than the mid-plane of the disc). We find for increasing values of inclination that a zero spin black hole has a mass range of 6,300 - 50,900 M_sun and a maximally spinning black hole has a mass between 16,900 - 191,700 M_sun. This is consistent with previous estimates and reinforces the idea that HLX-1 contains an intermediate mass black hole.

TIME EVOLUTION OF A VISCOUS PROTOPLANETARY DISK WITH A FREE GEOMETRY: TOWARD A MORE SELF-CONSISTENT PICTURE

Observations of protoplanetary disks show that some characteristics seem recurrent, even in star formation regions that are physically distant such as surface mass density profiles varying as $r^{-1}$, or aspect ratios about 0.03 to 0.23. Accretion rates are also recurrently found around $10^{-8} - 10^{-6}~\mathrm{M_{\odot}~yr^{-1}}$ for disks already evolved (Isella et al., 2009, Andrews et al., 2009, 2010). Several models have been developed in order to recover these properties. However, most of them usually simplify the disk geometry if not its mid-plane temperature. This has major consequences for modeling the disk evolution over million years and consequently planet migration. In the present paper, we develop a viscous evolution hydrodynamical numerical code that determines simultaneously the disk photosphere geometry and the mid-plane temperature. We then compare our results of long-term simulations with similar simulations of disks with a constrained geometry along the Chiang & Goldreich (1997) prescription (dlnH/dlnr = 9/7). We find that the constrained geometry models provide a good approximation of the disk surface density evolution. However, they differ significantly regarding the temperature time evolution. In addition, we find that shadowed regions naturally appear at the transition between viscously dominated and radiation dominated regions that falls in the region of planetary formation. We show that $\chi$ (photosphere height to pressure scale height ratio) cannot be considered as a constant, consistently with Watanabe et al. (2008). Comparisons with observations show that all disk naturally evolve toward a shallow surface density disk ($\Sigma \propto r^{-1}$). The mass flux across the disk stabilizes in about 1 million year typically.

VERTICAL STRUCTURE AND CORONAL POWER OF ACCRETION DISKS POWERED BY MAGNETOROTATIONAL-INSTABILITY TURBULENCE

In this paper we consider two outstanding intertwined problems in modern high-energy astrophysics: (1) the vertical thermal structure of an optically thick accretion disk heated by the dissipation of magnetohydrodynamic (MHD) turbulence driven by the magneto-rotational instability (MRI), and (2) determining the fraction of the accretion power released in the corona above the disk. For simplicity, we consider a gas-pressure-dominated disk and assume a constant opacity. We argue that the local turbulent dissipation rate due to the disruption of MRI channel flows by secondary parasitic instabilities should be uniform across most of the disk, almost up to the disk photosphere. We then obtain a self-consistent analytical solution for the vertical thermal structure of the disk, governed by the balance between the heating by MRI turbulence and the cooling by radiative diffusion. Next, we argue that the coronal power fraction is determined by the competition between the Parker instability, viewed as a parasitic instability feeding off of MRI channel flows, and other parasitic instabilities. We show that the Parker instability inevitably becomes important near the disk surface, leading to a certain lower limit on the coronal power. While most of the analysis in this paper focuses on the case of a disk threaded by an externally imposed vertical magnetic field, we also discuss the zero-net-flux case, in which the magnetic field is produced by the MRI dynamo itself, and show that most of our arguments and conclusions should be valid in this case as well.

The effect of advection at luminosities close to Eddington: The ULX in M 31

The transient, ultra-luminous X-ray source CXOM31 J004253.1+411422 in the Andromeda galaxy is most likely a 10 solar mass black hole, with super-Eddington luminosity at its peak. The XMM-Newton spectra taken during the decline then trace luminosities of 0.86-0.27 L_Edd. These spectra are all dominated by a hot disc component, which roughly follows a constant inner radius track in luminosity and temperature as the source declines. At the highest luminosity the disc structure should change due to advection of radiation through the disc. This advected flux can be partly released at lower radii thus modifying the spectral shape. To study the effect of advection at luminosities close to Eddington we employ a fully relativistic slim disc model, SLIMBH, that includes advective cooling and full radiative transfer through the photosphere based on TLUSTY. The model also incorporates relativistic photon ray-tracing from the proper location of the disc photosphere rather than the mid-plane as the slim disc is no longer geometrically thin. We find that these new models differ only slightly from the non-advective (standard) BHSPEC models even at the highest luminosities considered here. While both discs can fit the highest luminosity data, neither is a very good fit to the *lower* luminosities. This could indicate a missing physical process that acts in low luminosity discs and subsides as the disc luminosity approaches the Eddington limit.

The periodic spectroscopic variability of FU Orionis

FU Orionis systems are young stars undergoing outbursts of disc accretion and where the optical spectrum contains lines associated with both the disc photosphere and a wind component. Previous observations of the prototype FU Orionis have suggested that the wind lines and the photospheric lines are modulated with periods of 14.54 and 3.54 days respectively (Herbig et al. 2003). We have re-observed the system at higher spectral resolution, by monitoring variations of optical line profiles over 21 nights in 2007 and have found periods of 13.48 and 3.6 days in the wind and disc components consistent with the above: this implies variability mechanisms that are stable over at least a decade. In addition we have found: i) that the variations in the photospheric absorption lines are confined to the blue wing of the line (around -9km/s): we tentatively ascribe this to an orbiting hotspot in the disc which is obscured by a disc warp during its receding phase. ii) The wind period is manifested not only in blue-shifted Halpha absorption, but also in red-shifted emission of Halpha and Hbeta, as well as in blue-shifted absorption of Na I D, Li I and Fe II. iii) We find that the periodic modulation of blue-shifted Halpha absorption at around -100km/s, is phase lagged with respect to variations in the other lines by ~1.8days. This is consistent with a picture in which variations at the wind base first affect chromospheric emission and then low velocity blue-shifted absorption, followed - after a lag equal to the propagation time of disturbances across the wind's acceleration region - by a response in high velocity blue-shifted absorption. Such arguments constrain the size of the acceleration region to ~10^12cm. We discuss possible mechanisms for periodic variations within the innermost 0.1AU of the disc, including the possibility that these variations indicate the presence of an embedded hot Jupiter.

Testing slim-disk models on the thermal spectra of LMC X-3

Slim-disk models describe advective accretion flows at high luminosities, while reducing to the standard thin disk form in the low luminosity limit. We have developed a new spectral model, slimbb, within the framework of XSPEC, which describes fully relativistic slim-disk accretion and includes photon ray-tracing that starts from the disk photosphere, rather than the equatorial plane. We demonstrate the features of this model by applying it to RXTE spectra of the persistent black-hole X-ray binary LMC X-3. LMC X-3 has the virtues of exhibiting large intensity variations while maintaining itself in soft spectral states which are well described using accretion-disk models, making it an ideal candidate to test the aptness of slimbb. Our results demonstrate consistency between the low-luminosity (thin-disk) and high luminosity (slim-disk) regimes. The results also illustrate that advection alone does not solve the problem of the origin of the surprisingly soft high-luminosity spectra in LMC X-3. We show that X-ray continuum-fitting in the high accretion rate regime can powerfully test black-hole accretion disk models.

The imprint of photoevaporation on edge-on discs

We have performed hydrodynamic and radiative transfer calculations of a photoevaporating disc around a Herbig Ae/Be star to determine the evolution and observational impact of dust entrained in the wind. We find that the wind selectively entrains grains of different sizes at different radii resulting in a dust population that varies spatially and increases with height above the disc at radii > 10 AU. This variable grain population results in a 'wingnut' morphology to the dust density distribution. We calculate images of this dust distribution at NIR wavelengths that also show a wingnut morphology at all wavelengths considered. We have also considered the contribution that crystalline dust grains will have in the wind and show that a photoevaporative wind can result in a significant crystallinity fraction at all radii, when the disc is edge-on. However, when the disc's photosphere is unobscured, a photoevaporative wind makes no contribution to the observable crystallinity fraction in the disc. Finally, we conclude that the analysis of extended emission around edge-on discs could provide a new and independent method of testing photoevaporation models.

ON THE PROPERTIES OF THERMAL DISK WINDS IN X-RAY TRANSIENT SOURCES: A CASE STUDY OF GRO J1655–40

We present the results of hydrodynamical simulations of the disk photosphere irradiated by strong X-rays produced in the inner most part of the disk. As expected, the irradiation heats the photosphere and drives a thermal wind. To apply our results to the well-studied X-ray transient source GRO J1655-40, we adopted the observed mass of its black hole, and the observed properties of its X-ray radiation. To compare the results with the observations, we also computed transmitted X-ray spectra based on the wind solution. Our main finding is: the density of the fast moving part of the wind is more than one order of magnitude lower than that inferred from the observations. Consequently, the model fails to predict spectra with line absorption as strong and as blueshifted as those observed. However, despite the thermal wind being weak and Compton thin, the ratio between the mass-loss rate and the mass accretion rate is about seven. This high ratio is insensitive to the accretion luminosity, in the limit of lower luminosities. Most of the mass is lost from the disk between 0.07 and 0.2 of the Compton radius. We discovered that beyond this range the wind solution is self-similar. In particular, soon after it leaves the disk, the wind flows at a constant angle with respect to the disk. Overall, the thermal winds generated in our comprehensive simulations do not match the wind spectra observed in GRO J1655-40. This supports the conclusion of Miller et al. and Kallman et al. that the wind in GRO J1655-40, and possibly other X-ray transients, may be driven by magnetic processes. This in turn implies that the disk wind carries even more material than our simulations predict and as such has a very significant impact on the accretion disk structure and dynamics.

GIANT PLANET FORMATION BY DISK INSTABILITY: A COMPARISON SIMULATION WITH AN IMPROVED RADIATIVE SCHEME

There has been disagreement currently about whether cooling in protoplanetary disks can be sufficiently fast to induce the formation of gas giant protoplanets via gravitational instabilities. Simulations by our own group and others indicate that this method of planet formation does not work for disks around young, low- mass stars inside several tens of AU, while simulations by other groups show fragmentation into protoplanetary clumps in this region. To allow direct comparison in hopes of isolating the cause of the differences, we here present a high resolution three-dimensional hydrodynamics simulation of a protoplanetary disk, where the disk model, initial perturbation, and simulation conditions are essentially identical to those used in a set of simulations by Boss. As in earlier papers by the same author, Boss (2007, hereafter B07) purports to show that cooling is fast enough to produce protoplanetary clumps. Here, we evolve the same B07 disk using an improved version of one of our own radiative schemes and find that the disk does not fragment in our code but instead quickly settles into a state with only low amplitude nonaxisymmetric structure, which persists for at least several outer disk rotations. We see no rapid radiative or convective cooling. We conclude that the differences in results are due to different treatments of regions at and above the disk photosphere, and we explain at least one way in which the scheme in B07 may lead to artificially fast cooling.

TESTING ACCRETION DISK STRUCTURE WITHSUZAKUDATA OF LMC X-3

The Suzaku observation of LMC X-3 gives the best data to date on the shape of the accretion disk spectrum. This is due to the combination of very low absorbing column density along this line of sight which allows the shape of the disk emisison to be constrained at low energies by the CCD's, while the tail can be simultaneously determined up to 30 keV by the high energy detectors. These data clearly demonstrate that the observed disk spectrum is broader than a simple 'sum of blackbodies', and relativistic smearing of the emission is strongly required. However, the intrinsic emission should be more complex than a (color-corrected) sum of blackbodies as it should also contain photo-electric absorption edges from the partially ionised disk photosphere. These are broadened by the relativistic smearing, but the models predict ~ 3-5 per cent deviations for 1/3- 1 solar abundance around the edge energies, significantly stronger than observed. This indicate that the models need to include more physical processes such as self-irradiation, bound-bound (line) absorption, and/or emission from recombination continuua and/or lines. Alternatively, if none of these match the data, it may instead require that the accretion disk density and/or emissivity profile with height is different to that assumed. Thus these data demonstrate the feasibility of observational tests of our fundamental understanding of the vertical structure of accretion disks.

TWO-DIMENSIONAL SIMULATIONS OF FU ORIONIS DISK OUTBURSTS

We have developed time-dependent models of FU Ori accretion outbursts to explore the physical properties of protostellar disks. Our two-dimensional, axisymmetric models incorporate full vertical structure with a new treatment of the radiative boundary condition for the disk photosphere. We find that FU Ori-type outbursts can be explained by a slow accumulation of matter due to gravitational instability. Eventually this triggers the magnetorotational instability, which leads to rapid accretion. The thermal instability is triggered in the inner disk but this instability is not necessary for the outburst. An accurate disk vertical structure, including convection, is important for understanding the outburst behavior. Large convective eddies develop during the high state in the inner disk. The models are in agreement with Spitzer IRS spectra and also with peak accretion rates and decay timescales of observed outbursts, though some objects show faster rise timescale. We also propose that convection may account for the observed mild-supersonic turbulence and the short-timescale variations of FU Orionis objects.