Stochastic accretion of planetesimals on to white dwarfs: Constraints on the mass distribution of accreted material from atmospheric pollution

This paper explores how the stochastic accretion of planetesimals on to white dwarfs would be manifested in observations of their atmospheric pollution. Archival observations of pollution levels for unbiased samples of DA and non-DA white dwarfs are used to derive the distribution of inferred accretion rates, confirming that rates become systematically lower as sinking time (assumed here to be dominated by gravitational settling) is decreased, with no discernable dependence on cooling age. The accretion rates expected from planetesimals that are all the same mass (i.e., a mono-mass distribution) are explored both analytically and using a Monte Carlo model, quantifying how measured accretion rates inevitably depend on sinking time, since different sinking times probe different times since the last accretion event. However, that dependence is so dramatic that a mono-mass distribution can be excluded within the context of this model. Consideration of accretion from a broad distribution of planetesimal masses uncovers an important conceptual difference: accretion is continuous (rather than stochastic) for planetesimals below a certain mass, and the accretion of such planetesimals determines the rate typically inferred from observations; smaller planetesimals dominate the rates for shorter sinking times. A reasonable fit to the observationally inferred accretion rate distributions is found with model parameters consistent with a collisionally evolved mass distribution up to Pluto-mass, and an underlying accretion rate distribution consistent with that expected from descendants of debris discs of main-sequence A stars. With these parameters, while both DA and non-DA white dwarfs accrete from the same broad planetesimal distribution, this model predicts that the pollution seen in DAs is dominated by the continuous accretion of 35 km objects (though the dominant size varies between stars by around an order of magnitude from this reference value). Furthermore, observations that characterize the dependence of inferred accretion rates on sinking time and cooling age (including a consideration of the effect of thermohaline convection on models used to derive those rates), and the decadal variability of DA accretion signatures, will improve constraints on the mass distribution of accreted material and the lifetime of the disc through which it is accreted

A self-consistent model for the evolution of the gas produced in the debris disc of β Pictoris

This paper presents a self-consistent model for the evolution of gas produced in the debris disc of β Pictoris. Our model proposes that atomic carbon and oxygen are created from the photodissociation of CO, which is itself released from volatile-rich bodies in the debris disc due to grain–grain collisions or photodesorption. While the CO lasts less than one orbit, the atomic gas evolves by viscous spreading resulting in an accretion disc inside the parent belt and a decretion disc outside. The temperature, ionization fraction and population levels of carbon and oxygen are followed with the photodissociation region model CLOUDY, which is coupled to a dynamical viscous α model. We present new gas observations of β Pic, of C I observed with Atacama Pathfinder EXperiment and O I observed with $\textit{Herschel}$, and show that these along with published CII and CO observations can all be explained with this new model. Our model requires a viscosity α > 0.1, similar to that found in sufficiently ionized discs of other astronomical objects; we propose that the magnetorotational instability is at play in this highly ionized and dilute medium. This new model can be tested from its predictions for high-resolution ALMA observations of C I. We also constrain the water content of the planetesimals in β Pic. The scenario proposed here might be at play in all debris discs and this model could be used more generally on all discs with C, O or CO detections.QK, MW and LM acknowledge support from the European Union through ERC grant number 279973. AJ acknowledges the support of the DISCSIM project, grant agreement 341137, funded by the European Research Council under ERC-2013-ADG.This is the final version of the article. It first appeared from Oxford University Press via http://dx.doi.org/10.1093/mnras/stw136

Scars of intense accretion episodes at metal-rich white dwarfs

A re-evaluation of time-averaged accretion rates at DBZ-type white dwarfs points to historical, time-averaged rates significantly higher than the currently observed episodes at their DAZ counterparts. The difference between the ongoing, instantaneous accretion rates witnessed at DAZ white dwarfs, which often exceed 10 8gs -1, and those inferred over the past 10 5-10 6yr for the DBZ stars can be of a few orders of magnitude, and therefore must result from high-rate episodes of tens to hundreds of years so that they remain undetected to date. This paper explores the likelihood that such brief, intense accretion episodes of gas-phase material can account for existing data. For reasonable assumptions about the circumstellar gas, accretion rates approaching or exceeding 10 15gs -1 are possible, similar to rates observed in quiescent cataclysmic variables, and potentially detectable with future X-ray missions or wide-field monitoring facilities. Gaseous debris that is prone to such rapid accretion may be abundant immediately following a tidal disruption event via collisions and sublimation, or if additional bodies impinge upon an extant disc. Particulate disc matter accretes at or near the Poynting-Robertson drag rate for long periods between gas-producing events, consistent with rates inferred for dusty DAZ white dwarfs. In this picture, warm DAZ stars without infrared excesses have rates consistent with accretion from particulate discs that remain undetected. This overall picture has implications for quasi-steady state models of accretion and the derived chemical composition of asteroidal debris in DBZ white dwarfs

Photodissociation and the Morphology of HI in Galaxies

Young massive stars produce Far-UV photons which dissociate the molecular gas on the surfaces of their parent molecular clouds. Of the many dissociation products which result from this ``back-reaction'', atomic hydrogen \HI is one of the easiest to observe through its radio 21-cm hyperfine line emission. In this paper I first review the physics of this process and describe a simplified model which has been developed to permit an approximate computation of the column density of photodissociated \HI which appears on the surfaces of molecular clouds. I then review several features of the \HI morphology of galaxies on a variety of length scales and describe how photodissociation might account for some of these observations. Finally, I discuss several consequences which follow if this view of the origin of HI in galaxies continues to be successful.Comment: 18 pages, 7 figures in 8 files, invited review paper for the conference "Penetrating Bars Through Masks of Cosmic Dust: The Hubble Tuning Fork Strikes a New Note", South Africa, June 2004. Proceedings to be published by Kluwer, eds. D.L. Block, K.C. Freeman, I. Puerari, R. Groess, & E.K. Bloc

Accretion Disc Theory: From the Standard Model Until Advection

Accretion disc theory was first developed as a theory with the local heat balance, where the whole energy produced by a viscous heating was emitted to the sides of the disc. One of the most important new invention of this theory was a phenomenological treatment of the turbulent viscosity, known as ''alpha'' prescription, when the (r$\phi$) component of the stress tensor was approximated by ($\alpha$ P) with a unknown constant $\alpha$. This prescription played the role in the accretion disc theory as well important as the mixing-length theory of convection for stellar evolution. Sources of turbulence in the accretion disc are discussed, including nonlinear hydrodynamical turbulence, convection and magnetic field role. In parallel to the optically thick geometrically thin accretion disc models, a new branch of the optically thin accretion disc models was discovered, with a larger thickness for the same total luminosity. The choice between these solutions should be done of the base of a stability analysis. The ideas underlying the necessity to include advection into the accretion disc theory are presented and first models with advection are reviewed. The present status of the solution for a low-luminous optically thin accretion disc model with advection is discussed and the limits for an advection dominated accretion flows (ADAF) imposed by the presence of magnetic field are analysed.Comment: Roceeding of the Int. Workshop "Observational Evidence for Black Holes in the Universe". Calcutta, 11-17 January 1998. Kluwer Acad. Pu

Modeling Molecular-Line Emission from Circumstellar Disks

Molecular lines hold valuable information on the physical and chemical composition of disks around young stars, the likely progenitors of planetary systems. This invited contribution discusses techniques to calculate the molecular emission (and absorption) line spectrum based on models for the physical and chemical structure of protoplanetary disks. Four examples of recent research illutrate these techniques in practice: matching resolved molecular-line emission from the disk around LkCa15 with theoertical models for the chemistry; evaluating the two-dimensional transfer of ultraviolet radiation into the disk, and the effect on the HCN/CN ratio; far-infrared CO line emission from a superheated disk surface layer; and inward motions in the disk around L1489 IRS.Comment: 6 pages, no figures. To appear in "The Dense Interstellar Medium in Galaxies", Procs. Fourth Cologne-Bonn-Zermatt-Symposiu

Magnetospherically-trapped dust and a possible model for the unusual transits at WD 1145+017

The rapidly evolving dust and gas extinction observed towards WD 1145+017 has opened a real-time window on to the mechanisms for destruction-accretion of planetary bodies on to white dwarf stars, and has served to underline the importance of considering the dynamics of dust particles around such objects. Here it is argued that the interaction between (charged) dust grains and the stellar magnetic field is an important ingredient in understanding the physical distribution of infrared emitting particles in the vicinity of such white dwarfs. These ideas are used to suggest a possible model for WD 1145+017 in which the unusual transit shapes are caused by opaque clouds of dust trapped in the stellar magnetosphere. The model can account for the observed transit periodicities if the stellar rotation is near 4.5 h, as the clouds of trapped dust are then located near or within the co-rotation radius. The model requires the surface magnetic field to be at least around some tens of kG. In contrast to the eccentric orbits expected for large planetesimals undergoing tidal disintegration, the orbits of magnetospherically-trapped dust clouds are essentially circular, consistent with the observations

Phase transitions in biological membranes

Native membranes of biological cells display melting transitions of their lipids at a temperature of 10-20 degrees below body temperature. Such transitions can be observed in various bacterial cells, in nerves, in cancer cells, but also in lung surfactant. It seems as if the presence of transitions slightly below physiological temperature is a generic property of most cells. They are important because they influence many physical properties of the membranes. At the transition temperature, membranes display a larger permeability that is accompanied by ion-channel-like phenomena even in the complete absence of proteins. Membranes are softer, which implies that phenomena such as endocytosis and exocytosis are facilitated. Mechanical signal propagation phenomena related to nerve pulses are strongly enhanced. The position of transitions can be affected by changes in temperature, pressure, pH and salt concentration or by the presence of anesthetics. Thus, even at physiological temperature, these transitions are of relevance. There position and thereby the physical properties of the membrane can be controlled by changes in the intensive thermodynamic variables. Here, we review some of the experimental findings and the thermodynamics that describes the control of the membrane function.Comment: 23 pages, 15 figure

Order-of-magnitude speedup for steady states and traveling waves via Stokes preconditioning in Channelflow and Openpipeflow

Steady states and traveling waves play a fundamental role in understanding hydrodynamic problems. Even when unstable, these states provide the bifurcation-theoretic explanation for the origin of the observed states. In turbulent wall-bounded shear flows, these states have been hypothesized to be saddle points organizing the trajectories within a chaotic attractor. These states must be computed with Newton's method or one of its generalizations, since time-integration cannot converge to unstable equilibria. The bottleneck is the solution of linear systems involving the Jacobian of the Navier-Stokes or Boussinesq equations. Originally such computations were carried out by constructing and directly inverting the Jacobian, but this is unfeasible for the matrices arising from three-dimensional hydrodynamic configurations in large domains. A popular method is to seek states that are invariant under numerical time integration. Surprisingly, equilibria may also be found by seeking flows that are invariant under a single very large Backwards-Euler Forwards-Euler timestep. We show that this method, called Stokes preconditioning, is 10 to 50 times faster at computing steady states in plane Couette flow and traveling waves in pipe flow. Moreover, it can be carried out using Channelflow (by Gibson) and Openpipeflow (by Willis) without any changes to these popular spectral codes. We explain the convergence rate as a function of the integration period and Reynolds number by computing the full spectra of the operators corresponding to the Jacobians of both methods.Comment: in Computational Modelling of Bifurcations and Instabilities in Fluid Dynamics, ed. Alexander Gelfgat (Springer, 2018

The evolution of pitch angles of spiral arms

This is the final version. Available from Oxford University Press via the DOI in this recordIn spiral galaxies, the pitch angle, α, of the spiral arms is often proposed as a discriminator between theories for the formation of the spiral structure. In Lin–Shu density wave theory, α stays constant in time, being simply a property of the underlying galaxy. In other theories (e.g. tidal interaction, and self-gravity), it is expected that the arms wind up in time, so that to a first approximation cot α ∝ t. For these theories, it would be expected that a sample of galaxies observed at random times should show a uniform distribution of cot α. We show that a recent set of measurements of spiral pitch angles (Yu & Ho) is broadly consistent with this expectation.European Union FP7Science and Technology Facilities Council (STFC