-
SPH simulations of irradiation-driven warped accretion discs and the long periods in X-ray binaries
Authors:
Stephen B. Foulkes,
Carole A. Haswell,
James R. Murray
Abstract:
We present three dimensional smoothed particle hydrodynamics (SPH) calculations of irradiation-driven warping of accretion discs. Initially unwarped planar discs are unstable to the radiation reaction when the disc is illuminated by a central radiation source. The disc warps and tilts and precesses slowly in a retrograde direction; its shape continuously flexes in response to the changing orient…
▽ More
We present three dimensional smoothed particle hydrodynamics (SPH) calculations of irradiation-driven warping of accretion discs. Initially unwarped planar discs are unstable to the radiation reaction when the disc is illuminated by a central radiation source. The disc warps and tilts and precesses slowly in a retrograde direction; its shape continuously flexes in response to the changing orientation of the Roche potential. We simulate ten systems: eight X-ray binaries, one cataclysmic variable (CV), and a `generic' low mass X-ray binary (LMXB). We adopt system parameters from observations and tune a single parameter: our model X-ray luminosity ($L_{*}$) to reproduce the observed or inferred super-orbital periods. Without exception, across a wide range of parameter space, we find an astonishingly good match between the observed $L_{X}$ and the model $L_{*}$. We conclude irradiation-driven warping is the mechanism underlying the long periods in X-ray binaries. Our Her X-1 simulation simultaneously reproduces the observed $L_{X}$, the "main-" and "short-high" X-ray states and the orbital inclination. Our simulations of SS 433 give a maximum warp angle of $18.6^{\circ}$, a good match to the cone traced by the jets, but this angle is reached only in the outer disc. In all cases, the overall disc tilt is less than $\degrees{13}$ and the maximum disc warp is less than and or equal to $\degrees{21}$.
△ Less
Submitted 10 January, 2010;
originally announced January 2010.
-
The effect of a planet on the dust distribution in a 3D protoplanetary disk
Authors:
L. Fouchet,
S. T. Maddison,
J. -F. Gonzalez,
J. R. Murray
Abstract:
Aims: We investigate the behaviour of dust in protoplanetary disks under the action of gas drag in the presence of a planet. Our goal is twofold: to determine the spatial distribution of dust depending on grain size and planet mass, and therefore to provide a framework for interpretation of coming observations and future studies of planetesimal growth. Method: We numerically model the evolution…
▽ More
Aims: We investigate the behaviour of dust in protoplanetary disks under the action of gas drag in the presence of a planet. Our goal is twofold: to determine the spatial distribution of dust depending on grain size and planet mass, and therefore to provide a framework for interpretation of coming observations and future studies of planetesimal growth. Method: We numerically model the evolution of dust in a protoplanetary disk using a two-fluid (gas + dust) Smoothed Particle Hydrodynamics (SPH) code, which is non-self-gravitating and locally isothermal. The code follows the three dimensional distribution of dust in a protoplanetary disk as it interacts with the gas via aerodynamic drag. In this work, we present the evolution of a minimum mass solar nebula (MMSN) disk comprising 1% dust by mass in the presence of an embedded planet. We run a series of simulations which vary the grain size and planetary mass to see how they affect the resulting disk structure. Results: We find that gap formation is much more rapid and striking in the dust layer than in the gaseous disk and that a system with a given stellar, disk and planetary mass will have a completely different appearance depending on the grain size. For low mass planets in our MMSN disk, a gap can open in the dust disk while not in the gas disk. We also note that dust accumulates at the external edge of the planetary gap and speculate that the presence of a planet in the disk may enhance the formation of a second planet by facilitating the growth of planetesimals in this high density region.
△ Less
Submitted 30 August, 2007;
originally announced August 2007.
-
Comprehensive simulations of superhumps
Authors:
Amanda J. Smith,
Carole A. Haswell,
James R. Murray,
Michael R. Truss,
Stephen B. Foulkes
Abstract:
(Abridged) We use 3D SPH calculations with higher resolution, as well as with more realistic viscosity and sound-speed prescriptions than previous work to examine the eccentric instability which underlies the superhump phenomenon in semi-detached binaries. We illustrate the importance of the two-armed spiral mode in the generation of superhumps. Differential motions in the fluid disc cause conve…
▽ More
(Abridged) We use 3D SPH calculations with higher resolution, as well as with more realistic viscosity and sound-speed prescriptions than previous work to examine the eccentric instability which underlies the superhump phenomenon in semi-detached binaries. We illustrate the importance of the two-armed spiral mode in the generation of superhumps. Differential motions in the fluid disc cause converging flows which lead to strong spiral shocks once each superhump cycle. The dissipation associated with these shocks powers the superhump. We compare 2D and 3D results, and conclude that 3D simulations are necessary to faithfully simulate the disc dynamics. We ran our simulations for unprecedented durations, so that an eccentric equilibrium is established except at high mass ratios where the growth rate of the instability is very low. Our improved simulations give a closer match to the observed relationship between superhump period excess and binary mass ratio than previous numerical work. The observed black hole X-ray transient superhumpers appear to have systematically lower disc precession rates than the cataclysmic variables. This could be due to higher disc temperatures and thicknesses. The modulation in total viscous dissipation on the superhump period is overwhelmingly from the region of the disc within the 3:1 resonance radius. As the eccentric instability develops, the viscous torques are enhanced, and the disc consequently adjusts to a new equilibrium state, as suggested in the thermal-tidal instability model. We quantify this enhancement in the viscosity, which is ~10 per cent for q=0.08. We characterise the eccentricity distributions in our accretion discs, and show that the entire body of the disc partakes in the eccentricity.
△ Less
Submitted 11 April, 2007;
originally announced April 2007.
-
Populating the Galaxy with close Be binaries
Authors:
P. D. Kiel,
J. R. Hurley,
J. R. Murray,
K. Hayasaki
Abstract:
Be/X-ray binaries comprise roughly two-thirds of the high-mass X-ray binaries (HMXBs), which is a class of X-ray binaries that results from the high mass of the companion or donor star (> 10 solar masses). Currently the formation and evolution of X-ray producing Be binaries is a matter of great debate. Modelling of these systems requires knowledge of Be star evolution and also consideration of h…
▽ More
Be/X-ray binaries comprise roughly two-thirds of the high-mass X-ray binaries (HMXBs), which is a class of X-ray binaries that results from the high mass of the companion or donor star (> 10 solar masses). Currently the formation and evolution of X-ray producing Be binaries is a matter of great debate. Modelling of these systems requires knowledge of Be star evolution and also consideration of how the evolution changes when the star is in close proximity to a companion. Within this work we complete a full population synthesis study of Be binaries for the Galaxy. The results for the first time match aspects of the observational data, most notably the observed upper limit to the period distribution. We conclude that greater detailed studies on the evolution of Be stars within X-ray binaries needs to be completed, so that rapid binary evolution population synthesis packages may best evolve these systems.
△ Less
Submitted 7 March, 2007;
originally announced March 2007.
-
Three dimensional SPH simulations of radiation-driven warped accretion discs
Authors:
Stephen B. Foulkes,
Carole A. Haswell,
James R. Murray
Abstract:
We present three dimensional smoothed particle hydrodynamics (SPH) calculations of warped accretion discs in X-ray binary systems. Geometrically thin, optically thick accretion discs are illuminated by a central radiation source. This illumination exerts a non-axisymmetric radiation pressure on the surface of the disc resulting in a torque that acts on the disc to induce a twist or warp. Initial…
▽ More
We present three dimensional smoothed particle hydrodynamics (SPH) calculations of warped accretion discs in X-ray binary systems. Geometrically thin, optically thick accretion discs are illuminated by a central radiation source. This illumination exerts a non-axisymmetric radiation pressure on the surface of the disc resulting in a torque that acts on the disc to induce a twist or warp. Initially planar discs are unstable to warping driven by the radiation torque and in general the warps also precess in a retrograde direction relative to the orbital flow. We simulate a number of X-ray binary systems which have different mass ratios using a number of different luminosities for each. Radiation-driven warping occurs for all systems simulated. For mass ratios q ~ 0.1 a moderate warp occurs in the inner disc while the outer disc remains in the orbital plane (c.f. X 1916-053). For less extreme mass ratios the entire disc tilts out of the orbital plane (c.f. Her X-1). For discs that are tilted out of the orbital plane in which the outer edge material of the disc is precessing in a prograde direction we obtain both positive and negative superhumps simultaneously in the dissipation light curve (c.f. V603 Aql).
△ Less
Submitted 6 January, 2006;
originally announced January 2006.
-
SPH simulations of accretion flow via Roche lobe overflow and via mass transfer from Be disk
Authors:
Kimitake Hayasaki,
Atsuo T. Okazaki,
James. R. Murray
Abstract:
We compare the accretion flow onto the neutron star induced by Roche lobe overflow with that by the overflow from the Be disk, in a zero eccentricity, short period binary with the same mass transfer rate, performing three-dimensional Smoothed Particle Hydrodynamics simulations. We find that a persistent accretion disk is formed around the neutron star in both cases. The circularization radius of…
▽ More
We compare the accretion flow onto the neutron star induced by Roche lobe overflow with that by the overflow from the Be disk, in a zero eccentricity, short period binary with the same mass transfer rate, performing three-dimensional Smoothed Particle Hydrodynamics simulations. We find that a persistent accretion disk is formed around the neutron star in both cases. The circularization radius of the material transferred via Roche lobe overflow is larger than that of the material transfered from the Be disk. Thus, the growth of the accretion disk in the former case becomes significantly slower than in the latter case. In both cases, the mass accretion rate is very small and varies little with orbital phase, which is consistent with the observed X-ray behaviour of Be/X-ray binaries with circular orbits (e.g. XTE J1543-568).
△ Less
Submitted 17 November, 2005;
originally announced November 2005.
-
Discovery and modelling of disc precession in the M31 X-ray binary Bo 158?
Authors:
R. Barnard,
S. B. Foulkes,
C. A. Haswell,
U. Kolb,
J. P. Osborne,
J. R. Murray
Abstract:
The low mass X-ray binary (LMXB) associated with the M31 globular cluster Bo 158 is known to exhibit intensity dips on a ~2.78 hr period. This is due to obscuration of the X-ray source on the orbital period by material on the outer edge of the accretion disc. However, the depth of dipping varied from <10% to \~83% in three archival XMM-Newton observations of Bo 158. Previous work suggested that…
▽ More
The low mass X-ray binary (LMXB) associated with the M31 globular cluster Bo 158 is known to exhibit intensity dips on a ~2.78 hr period. This is due to obscuration of the X-ray source on the orbital period by material on the outer edge of the accretion disc. However, the depth of dipping varied from <10% to \~83% in three archival XMM-Newton observations of Bo 158. Previous work suggested that the dip depth was anticorrelated with the X-ray luminosity. However, we present results from three new XMM-Newton observations that suggest that the evolution of dipping is instead due to precession of the accretion disc. Such precession is expected in neutron star LMXBs with mass ratios <0.3 (i.e. with orbital periods <4 hr), such as the Galactic dipping LMXB 4U 1916-053. We simulated the accretion disc of Bo 158 using cutting-edge 3D smoothed particle hydrodynamics (SPH), and using the observed parameters. Our results show disc variability on two time-scales. The disc precesses in a prograde direction on a period of 81+/-3 hr. Also, a radiatively-driven disc warp is present in the inner disc, which undergoes retrograde precesson on a \~31 hr period. From the system geometry, we conclude that the dipping evolution is driven by the disc precession. Hence we predict that the dipping behaviour repeats on a ~81 hr cycle.
△ Less
Submitted 2 November, 2005;
originally announced November 2005.
-
Dust distribution in protoplanetary disks - Vertical settling and radial migration
Authors:
Laure Barrière-Fouchet,
Jean-François Gonzalez,
James R. Murray,
Robin J. Humble,
Sarah T. Maddison
Abstract:
We present the results of a three dimensional, locally isothermal, non-self-gravitating SPH code which models protoplanetary disks with two fluids: gas and dust. We ran simulations of a 1 Msun star surrounded by a 0.01 Msun disk comprising 99% gas and 1% dust in mass and extending from 0.5 to ~300 AU. The grain size ranges from 0.001 mm to 10 m for the low resolution (~25 000 SPH particles) simu…
▽ More
We present the results of a three dimensional, locally isothermal, non-self-gravitating SPH code which models protoplanetary disks with two fluids: gas and dust. We ran simulations of a 1 Msun star surrounded by a 0.01 Msun disk comprising 99% gas and 1% dust in mass and extending from 0.5 to ~300 AU. The grain size ranges from 0.001 mm to 10 m for the low resolution (~25 000 SPH particles) simulations and from 0.1 mm to 10 cm for the high resolution (~160 000 SPH particles) simulations. Dust grains are slowed down by the sub-Keplerian gas and lose angular momentum, forcing them to migrate towards the central star and settle to the midplane. The gas drag efficiency varies according to the grain size, with the larger bodies being weakly influenced and following marginally perturbed Keplerian orbits, while smaller grains are strongly coupled to the gas. For intermediate sized grains, the drag force decouples the dust and gas, allowing the dust to preferentially migrate radially and efficiently settle to the midplane. The resulting dust distributions for each grain size will indicate, when grain growth is added, the regions when planets are likely to form.
△ Less
Submitted 22 August, 2005;
originally announced August 2005.
-
Spiral waves and the secondary in the novalike variable, V3885 Sgr
Authors:
L. E. Hartley,
J. R. Murray,
J. E. Drew,
K. S. Long
Abstract:
We present seven nights' blue spectroscopy of the novalike variable star, V3885 Sgr. The line spectrum shows a typical combination of broad absorption and emission in H gamma and H beta and He I, which is associated with the accretion disk. We also observe anti-phased narrow emission, which we attribute to irradiation of the secondary star. The He II 4686 and N III--C III--C IV emission lines ar…
▽ More
We present seven nights' blue spectroscopy of the novalike variable star, V3885 Sgr. The line spectrum shows a typical combination of broad absorption and emission in H gamma and H beta and He I, which is associated with the accretion disk. We also observe anti-phased narrow emission, which we attribute to irradiation of the secondary star. The He II 4686 and N III--C III--C IV emission lines are devoid of structure and are most likely formed in an outflow. We measure radial velocity shifts in the absorption and emission lines, from which we fit an orbital period of 4.97126+/-0.00036h. From the velocity semi-amplitudes of the disk and companion star, we are able to constrain the binary mass ratio to q>0.7.
The phase-folded spectra provide dense coverage of the entire orbital cycle. Doppler tomograms of the hydrogen and He I lines reveal spiral structure in the accretion disk and the irradiated donor star. We believe this is the first unambiguous detection of spiral waves in a novalike variable.
△ Less
Submitted 20 July, 2005;
originally announced July 2005.
-
Discovery of disc precession in the M31 dipping X-ray binary Bo 158
Authors:
R. Barnard,
S. B. Foulkes,
C. A. Haswell,
U. Kolb,
J. P. Osborne,
J. R. Murray
Abstract:
We present results from three XMM-Newton observations of the M31 low mass X-ray binary XMMU J004314.4+410726.3 (Bo 158), spaced over 3 days in 2004, July. Bo 158 was the first dipping LMXB to be discovered in M31. Periodic intensity dips were previously seen to occur on a 2.78-hr period, due to absorption in material that is raised out of the plane of the accretion disc. The report of these obse…
▽ More
We present results from three XMM-Newton observations of the M31 low mass X-ray binary XMMU J004314.4+410726.3 (Bo 158), spaced over 3 days in 2004, July. Bo 158 was the first dipping LMXB to be discovered in M31. Periodic intensity dips were previously seen to occur on a 2.78-hr period, due to absorption in material that is raised out of the plane of the accretion disc. The report of these observations stated that the dip depth was anti-correlated with source intensity. However, our new observations do not favour a strict intensity dependance, but rather suggest that the dip variation is due to precession of the accretion disc. This is to be expected in LMXBs with a mass ratio <~ 0.3 (period <~ 4 hr), as the disc reaches the 3:1 resonance with the binary companion, causing elongation and precession of the disc. A smoothed particle hydrodynamics simulation of the disc in this system shows retrograde rotation of a disc warp on a period of ~11 P_orb, and prograde disc precession on a period of ~29 P_orb. This is consistent with the observed variation in the depth of the dips. We find that the dipping behaviour is most likely to be modified by the disc precession, hence we predict that the dipping behaviour repeats on a 81+/-3 hr cycle.
△ Less
Submitted 8 November, 2005; v1 submitted 11 March, 2005;
originally announced March 2005.
-
Hydrodynamic modelling of accretion flows
Authors:
J. R. Murray,
M. R. Truss,
S. B. Foulkes,
C. A. Haswell,
K. Manson
Abstract:
In the proceedings of this, and of several recent close binary conferences, there have been several contributions describing smoothed particle hydrodynamics simulations of accretion disks. It is apposite therefore to review the numerical scheme itself with emphasis on its advantages for disk modelling, and the methods used for modelling viscous processes.
In the proceedings of this, and of several recent close binary conferences, there have been several contributions describing smoothed particle hydrodynamics simulations of accretion disks. It is apposite therefore to review the numerical scheme itself with emphasis on its advantages for disk modelling, and the methods used for modelling viscous processes.
△ Less
Submitted 18 March, 2004;
originally announced March 2004.
-
Simulations of spectral lines from an eccentric precessing accretion disc
Authors:
S. B. Foulkes,
C. A. Haswell,
J. R. Murray,
D. J. Rolfe
Abstract:
Two dimensional SPH simulations of a precessing accretion disc in a q=0.1 binary system (such as XTE J1118+480) reveal complex and continuously varying shape, kinematics, and dissipation. The stream-disc impact region and disc spiral density waves are prominent sources of energy dissipation.The dissipated energy is modulated on the period P_{sh} = ({P_{orb}}^{-1}-{P_{prec}}^{-1}^{-1} with which…
▽ More
Two dimensional SPH simulations of a precessing accretion disc in a q=0.1 binary system (such as XTE J1118+480) reveal complex and continuously varying shape, kinematics, and dissipation. The stream-disc impact region and disc spiral density waves are prominent sources of energy dissipation.The dissipated energy is modulated on the period P_{sh} = ({P_{orb}}^{-1}-{P_{prec}}^{-1}^{-1} with which the orientation of the disc relative to the mass donor repeats. This superhump modulation in dissipation energy has a variation in amplitude of ~10% relative to the total dissipation energy and evolves, repeating exactly only after a full disc precession cycle. A sharp component in the light curve is associated with centrifugally expelled material falling back and impacting the disc. Synthetic trailed spectrograms reveal two distinct "S-wave" features, produced respectively by the stream gas and the disc gas at the stream-disc impact shock. These S-waves are non-sinusoidal, and evolve with disc precession phase. We identify the spiral density wave emission in the trailed spectrogram. Instantaneous Doppler maps show how the stream impact moves in velocity space during an orbit. In our maximum entropy Doppler tomogram the stream impact region emission is distorted, and the spiral density wave emission is uppressed. A significant radial velocity modulation of the whole line profile occurs on the disc precession period. We compare our SPH simulation with a simple 3D model: the former is appropriate for comparison with emission lines while the latter is preferable for skewed absorption lines from precessing discs.
△ Less
Submitted 19 January, 2004;
originally announced January 2004.
-
Predicting Dust Distribution in Protoplanetary Discs
Authors:
Laure Barriere-Fouchet,
Jean-Francois Gonzalez,
Robin J. Humble,
Sarah T. Maddison,
James R. Murray
Abstract:
We present the results of three-dimensional numerical simulations that include the effects of hydrodynamical forces and gas drag upon an evolving dusty gas disk. We briefly describe a new parallel, two phase numerical code based upon the smoothed particle hydrodynamics (SPH) technique in which the gas and dust phases are represented by two distinct types of particles. We use the code to follow t…
▽ More
We present the results of three-dimensional numerical simulations that include the effects of hydrodynamical forces and gas drag upon an evolving dusty gas disk. We briefly describe a new parallel, two phase numerical code based upon the smoothed particle hydrodynamics (SPH) technique in which the gas and dust phases are represented by two distinct types of particles. We use the code to follow the dynamical evolution of a population of grains in a gaseous protoplanetary disk in order to understand the distribution of grains of different sizes within the disk. Our ``grains'' range from metre to submillimetre in size.
△ Less
Submitted 31 October, 2003;
originally announced October 2003.
-
Building Planets with Dusty Gas
Authors:
S. T. Maddison,
R. J. Humble,
J. R. Murray
Abstract:
We have developed a new numerical technique for simulating dusty-gas flows. Our unique code incorporates gas hydrodynamics, self-gravity and dust drag to follow the dynamical evolution of a dusty-gas medium. We have incorporated several descriptions for the drag between gas and dust phases and can model flows with submillimetre, centimetre and metre size "dust". We present calculations run on th…
▽ More
We have developed a new numerical technique for simulating dusty-gas flows. Our unique code incorporates gas hydrodynamics, self-gravity and dust drag to follow the dynamical evolution of a dusty-gas medium. We have incorporated several descriptions for the drag between gas and dust phases and can model flows with submillimetre, centimetre and metre size "dust". We present calculations run on the APAC* supercomputer following the evolution of the dust distribution in the pre-solar nebula. *Australian Partnership for Advanced Computing
△ Less
Submitted 9 September, 2002;
originally announced September 2002.
-
Superhumps in Low-Mass X-Ray Binaries
Authors:
Carole A. Haswell,
Andrew R. King,
James R. Murray,
Philip A. Charles
Abstract:
We propose a mechanism for the superhump modulations observed in optical photometry of at least two black hole X-ray transients (SXTs). As in extreme mass-ratio cataclysmic variables (CVs), superhumps are assumed to result from the presence of the 3:1 orbital resonance in the accretion disc. This causes the disc to become non-axisymmetric and precess. However the mechanism for superhump luminosi…
▽ More
We propose a mechanism for the superhump modulations observed in optical photometry of at least two black hole X-ray transients (SXTs). As in extreme mass-ratio cataclysmic variables (CVs), superhumps are assumed to result from the presence of the 3:1 orbital resonance in the accretion disc. This causes the disc to become non-axisymmetric and precess. However the mechanism for superhump luminosity variations in low mass X-ray binaries (LMXBs) must differ from that in CVs, where it is attributed to a tidally-driven modulation of the disc's viscous dissipation, varying on the beat between the orbital and disc precession period. By contrast in LMXBs, tidal dissipation in the outer accretion disc is negligible: the optical emission is overwhelming dominated by reprocessing of intercepted central X-rays. Thus a different origin for the superhump modulation is required. Recent observations and numerical simulations indicate that in an extreme mass-ratio system the disc area changes on the superhump period. We deduce that the superhumps observed in SXTs arise from a modulation of the reprocessed flux by the changing area. Therefore, unlike the situation in CVs, where the superhump amplitude is inclination-independent, superhumps should be best seen in low-inclination LMXBs, whereas an orbital modulation from the heated face of the secondary star should be more prominent at high inclinations. Modulation at the disc precession period (10s of days) may indicate disc asymmetries such as warping. We comment on the orbital period determinations of LMXBs, and the possibility and significance of possible permanent superhump LMXBs.
△ Less
Submitted 21 September, 2000; v1 submitted 23 August, 2000;
originally announced August 2000.
-
Eccentric discs in binaries with intermediate mass ratios: Superhumps in the VY Sculptoris stars
Authors:
J. R. Murray,
B. Warner,
D. T. Wickramasinghe
Abstract:
We investigate the role of the eccentric disc resonance in systems with mass ratios q greater than 1/4, and demonstrate the effects that changes in the mass flux from the secondary star have upon the disc radius and structure. The addition of material with low specific angular momentum to its outer edge restricts a disc radially. Should the mass flux from the secondary be reduced, it is possible…
▽ More
We investigate the role of the eccentric disc resonance in systems with mass ratios q greater than 1/4, and demonstrate the effects that changes in the mass flux from the secondary star have upon the disc radius and structure. The addition of material with low specific angular momentum to its outer edge restricts a disc radially. Should the mass flux from the secondary be reduced, it is possible for the disc in a system with mass ratio as large as 1/3 to expand to the 3:1 eccentric inner Lindblad resonance and for superhumps to be excited.
△ Less
Submitted 23 February, 2000;
originally announced February 2000.
-
The precession of eccentric discs in close binaries
Authors:
James R. Murray
Abstract:
We consider the precession rates of eccentric discs in close binaries, and compare theoretical predictions with the results of numerical disc simulations and with observed superhump periods. A simple dynamical model for precession is found to be inadequate. For mass ratios less than approximately 1/4 a linear dynamical model does provide an upper limit for disc precession rates. Theory suggests…
▽ More
We consider the precession rates of eccentric discs in close binaries, and compare theoretical predictions with the results of numerical disc simulations and with observed superhump periods. A simple dynamical model for precession is found to be inadequate. For mass ratios less than approximately 1/4 a linear dynamical model does provide an upper limit for disc precession rates. Theory suggests that pressure forces have a significant retrograde impact upon the precession rate (Lubow 1992). We find that the disc precession rates for three systems with accurately known mass ratios are significantly slower than predicted by the dynamical theory, and we attribute the difference to pressure forces. By assuming that pressure forces of similar magnitude occur in all superhumping systems, we obtain an improved fit to superhump observations.
△ Less
Submitted 25 November, 1999;
originally announced November 1999.
-
Angular Momentum Transfer in the Binary X-ray Pulsar GX 1+4
Authors:
J. G. Greenhill,
D. K. Galloway,
J. R. Murray
Abstract:
We describe three presentations relating to the X-ray pulsar GX 1+4 at a workshop on magnetic fields and accretion at the Astrophysical Theory Centre, Australian National University on 1998, November 12-13. Optical and X-ray spectroscopy indicate that GX 1+4 is seen through a cloud of gravitationaly bound matter. We discuss an unstable negative feedback mechanism (originally proposed by Kotani e…
▽ More
We describe three presentations relating to the X-ray pulsar GX 1+4 at a workshop on magnetic fields and accretion at the Astrophysical Theory Centre, Australian National University on 1998, November 12-13. Optical and X-ray spectroscopy indicate that GX 1+4 is seen through a cloud of gravitationaly bound matter. We discuss an unstable negative feedback mechanism (originally proposed by Kotani et al, 1999), based on X-ray heating of this matter which controls the accretion rate when the source is in a low X-ray luminosity state. A deep minimum lasting ~6 hours occurred during observations with the RXTE satellite over 1996, July 19-21. The shape of the X-ray pulses changed remarkably from before to after the minimum. These changes may be related to the transition from neutron star spin-down to spin-up which occurred at about the same time. Smoothed particle hydrodynamic simulations of the effect of adding matter with opposite angular momentum to an existing disc, show that it is possible for a number of concentric rings with alternating senses of rotation to co-exist in a disc. This could provide an explanation for the step-like changes in Pdot which are observed in GX 1+4. Changes at the inner boundary of the disc occur at the same timescale as that imposed at the outer boundary. Reversals of material torque on the neutron star occur at a minimum in L_X.
△ Less
Submitted 23 September, 1999;
originally announced September 1999.
-
Accretion disk reversal and the spin-up/spin-down of accreting pulsars
Authors:
James R. Murray,
Martijn de Kool,
Jianke Li
Abstract:
We numerically investigate the hydrodynamics of accretion disk reversal and relate our findings to the observed spin-rate changes in the accreting X-ray pulsar GX~1+4. In this system, which accretes from a slow wind, the accretion disk contains two dynamically distinct regions. In the inner part viscous forces are dominant and disk evolution occurs on a viscous timescale. In the outer part dynam…
▽ More
We numerically investigate the hydrodynamics of accretion disk reversal and relate our findings to the observed spin-rate changes in the accreting X-ray pulsar GX~1+4. In this system, which accretes from a slow wind, the accretion disk contains two dynamically distinct regions. In the inner part viscous forces are dominant and disk evolution occurs on a viscous timescale. In the outer part dynamical mixing of material with opposite angular momentum is more important, and the externally imposed angular momentum reversal timescale governs the flow. In this outer region the disk is split into concentric rings of material with opposite senses of rotation that do not mix completely but instead remain distinct, with a clear gap between them. We thus predict that torque reversals resulting from accretion disk reversals will be accompanied by minima in accretion luminosity.
△ Less
Submitted 7 October, 1998;
originally announced October 1998.
-
Tilted accretion discs in cataclysmic variables: tidal instabilities and superhumps
Authors:
J. R. Murray,
P. J. Armitage
Abstract:
We investigate the growth of tidal instabilities in accretion discs in a binary star potential, using three dimensional numerical simulations. As expected from analytic work, the disc is prone to an eccentric instability provided that it is large enough to extend to the 3:1 resonance. The eccentric disc leads to positive superhumps in the light curve. It has been proposed that negative superhump…
▽ More
We investigate the growth of tidal instabilities in accretion discs in a binary star potential, using three dimensional numerical simulations. As expected from analytic work, the disc is prone to an eccentric instability provided that it is large enough to extend to the 3:1 resonance. The eccentric disc leads to positive superhumps in the light curve. It has been proposed that negative superhumps might arise from a tilted disc, but we find no evidence that the companion gravitational tilt instability can grow fast enough in a fluid disc to create a measurable inclination. The origin of negative superhumps in the light curves of cataclysmic variables remains a puzzle.
△ Less
Submitted 10 September, 1998;
originally announced September 1998.
-
The effects of tidally induced disc structure on white dwarf accretion in intermediate polars
Authors:
J. R. Murray,
P. J. Armitage,
L. Ferrario,
D. T Wickramasinghe
Abstract:
We investigate the effects of tidally induced asymmetric disc structure on accretion onto the white dwarf in intermediate polars. Using numerical simulation, we show that it is possible for tidally induced spiral waves to propagate sufficiently far into the disc of an intermediate polar that accretion onto the central white dwarf could be modulated as a result. We suggest that accretion from the…
▽ More
We investigate the effects of tidally induced asymmetric disc structure on accretion onto the white dwarf in intermediate polars. Using numerical simulation, we show that it is possible for tidally induced spiral waves to propagate sufficiently far into the disc of an intermediate polar that accretion onto the central white dwarf could be modulated as a result. We suggest that accretion from the resulting asymmetric inner disc may contribute to the observed X-ray and optical periodicities in the light curves of these systems. In contrast to the stream-fed accretion model for these periodicities, the tidal picture predicts that modulation can exist even for systems with weaker magnetic fields where the magnetospheric radius is smaller than the radius of periastron of the mass transfer stream. We also predict that additional periodic components should exist in the emission from low mass ratio intermediate polars displaying superhumps.
△ Less
Submitted 10 September, 1998;
originally announced September 1998.
-
Simulations of spiral structure in the accretion disc of IP Pegasi during outburst
Authors:
P. J. Armitage,
J. R. Murray
Abstract:
We consider the implications of the detection of spiral structure in the accretion disc of the binary IP Pegasi. We use numerical simulations of the development of a disc outburst to construct predicted Doppler tomograms, which are found to be in close agreement with the observations if the spiral pattern arises as a transient feature when the disc expands viscously at the start of the outburst.…
▽ More
We consider the implications of the detection of spiral structure in the accretion disc of the binary IP Pegasi. We use numerical simulations of the development of a disc outburst to construct predicted Doppler tomograms, which are found to be in close agreement with the observations if the spiral pattern arises as a transient feature when the disc expands viscously at the start of the outburst. The good agreement of such viscous disc simulations with the data is consistent with models in which most of the angular momentum transport in the disc originates in internal stresses rather than globally excited waves or shocks. Future detailed observations of the development of transient spiral features offer the potential to measure the dependence of the disc viscosity on the local physical conditions in the disc.
△ Less
Submitted 28 April, 1998;
originally announced April 1998.
-
Angular Momentum Transport in Simulations of Accretion Disks
Authors:
James Rhys Murray
Abstract:
In this paper we briefly discuss the ways in which angular momentum transport is included in simulations of non-self-gravitating accretion disks, concentrating on disks in close binaries. Numerical approaches fall in two basic categories; particle based Lagrangian schemes, and grid based Eulerian techniques. Underlying the choice of numerical technique are assumptions that are made about disk ph…
▽ More
In this paper we briefly discuss the ways in which angular momentum transport is included in simulations of non-self-gravitating accretion disks, concentrating on disks in close binaries. Numerical approaches fall in two basic categories; particle based Lagrangian schemes, and grid based Eulerian techniques. Underlying the choice of numerical technique are assumptions that are made about disk physics, in particular about the angular momentum transport mechanism. Grid-based simulations have generally been of hot, relatively inviscid disks whereas particle-based simulations are more commonly of cool, viscous disks. Calculations of the latter type have been instrumental in developing a model for the superhump phenomenon. We describe how we use an artificial viscosity term to introduce angular momentum transport into our smoothed particle hydrodynamics (SPH) disk code.
△ Less
Submitted 10 March, 1997;
originally announced March 1997.
-
SPH simulations of tidally unstable accretion disks in cataclysmic variables
Authors:
James R. Murray
Abstract:
We numerically study the precessing disk model for superhump in the SU~UMa subclass of cataclysmic variables, using a two dimensional SPH code specifically designed for thin disk problems.
Two disk simulations for a binary with mass ratio $q=\frac{3}{17}$ (similar to OY~Car) are performed, in order to investigate the Lubow (1991 a,b) tidal resonance instability mechanism. In the first calculati…
▽ More
We numerically study the precessing disk model for superhump in the SU~UMa subclass of cataclysmic variables, using a two dimensional SPH code specifically designed for thin disk problems.
Two disk simulations for a binary with mass ratio $q=\frac{3}{17}$ (similar to OY~Car) are performed, in order to investigate the Lubow (1991 a,b) tidal resonance instability mechanism. In the first calculation, a disk evolves under steady mass transfer from $L_1$. In the second simulation, mass is added in Keplerian orbit to the inner disk.
The two disks follow similar evolutionary paths. However the $L_1$ stream-disk interaction is found to slow the disk's radial expansion and to circularise gas orbits.
The initial eccentricity growth in our simulations is exponential at a rate slightly less than predicted by Lubow (1991a). We do not observe a clearing of material from the resonance region via the disk's tidal response to the $m=2$ component of the binary potential as was described in Lubow (1992). Instead the $m=2$ response weakens as the disk eccentricty increases.
Both disks reach an eccentric equilibrium state, in which they undergo prograde precession. The rate of viscous energy dissipation in the disks has a periodic excess with a period matching the disk's rotation. The source is identified as a large region in the outer disk, and the mechanism by which it is produced is identified. The time taken for the periodic excess to develop is consistent with the first appearance of superhumps in a superoutburst.
△ Less
Submitted 8 November, 1995;
originally announced November 1995.
-
SPH Simulations of Accretion Disks and Narrow Rings
Authors:
Sarah T. Maddison,
James R. Murray,
Joe J. Monaghan
Abstract:
We model a massless viscous disk using Smoothed Particle Hydrodynamics (SPH) and note that it evolves according to the Lynden-Bell \& Pringle theory (1974) until a non-axisymmetric instability develops at the inner edge of the disk. This instability may have the same origin as the instability of initially axisymmetric viscous disks discussed by Lyubarskij et al. (1994). To clarify the evolution…
▽ More
We model a massless viscous disk using Smoothed Particle Hydrodynamics (SPH) and note that it evolves according to the Lynden-Bell \& Pringle theory (1974) until a non-axisymmetric instability develops at the inner edge of the disk. This instability may have the same origin as the instability of initially axisymmetric viscous disks discussed by Lyubarskij et al. (1994). To clarify the evolution we evolved single and double rings of particles. It is actually inconsistent with the SPH scheme to set up a single ring as an initial condition because SPH assumes a smoothed initial state. As would be expected from an SPH simulation, the ring rapidly breaks up into a band. We analyse the stability of the ring and show that the predictions are confirmed by the simulation.
△ Less
Submitted 12 July, 1995;
originally announced July 1995.