sort results by

Use logical operators AND, OR, NOT and round brackets to construct complex queries. Whitespace-separated words are treated as ANDed.

Show articles per page in mode

Terquem, Caroline E. J. M. L. J.

Normalized to: Terquem, C.

43 article(s) in total. 35 co-authors, from 1 to 23 common article(s). Median position in authors list is 1,0.

[1]  oai:arXiv.org:1909.08476  [pdf] - 1963800
Non-adiabatic tidal oscillations induced by a planetary companion
Comments:
Submitted: 2019-09-18
We calculate the dynamical tides raised by a close planetary companion on non-rotating stars of $1$ $\text{M}_{\odot}$ and $1.4$ $\text{M}_{\odot}$. Using the Henyey method, we solve the fully non-adiabatic equations throughout the star. The horizontal Lagrangian displacement is found to be 10 to 100 times larger than the equilibrium tide value in a thin region near the surface of the star. This is because non--adiabatic effects dominate in a region that extends from below the outer edge of the convection zone up to the stellar surface, and the equilibrium tide approximation is inconsistent with non--adiabaticity. Although this approximation generally applies in the low frequency limit, it also fails in the parts of the convection zone where the forcing frequency is small but larger than the Brunt-V\"ais\"al\"a frequency. We derive analytical estimates which give a good approximation to the numerical values of the magnitude of the ratio of the horizontal and radial displacements at the surface. The relative surface flux perturbation is also significant, on the order of 0.1 % for a system modelled on 51 Pegasi b. Observations affected by the horizontal displacement may therefore be more achievable than previously thought, and brightness perturbations may be the result of flux perturbations rather than due to the radial displacement. We discuss the implication of this on the possibility of detecting such tidally excited oscillations, including the prospect of utilising the large horizontal motion for observations of systems such as 51 Pegasi.
[2]  oai:arXiv.org:1809.10042  [pdf] - 1767693
First-order mean motion resonances in two-planet systems: general analysis and observed systems
Comments: Accepted for publication in MNRAS
Submitted: 2018-09-26
This paper focuses on two-planet systems in a first-order $(q+1):q$ mean motion resonance and undergoing type-I migration in a disc. We present a detailed analysis of the resonance valid for any value of $q$. Expressions for the equilibrium eccentricities, mean motions and departure from exact resonance are derived in the case of smooth convergent migration. We show that this departure, not assumed to be small, is such that period ratio normally exceeds, but can also be less than, $ (q+1)/q.$ Departure from exact resonance as a function of time for systems starting in resonance and undergoing divergent migration is also calculated. We discuss observed systems in which two low mass planets are close to a first-order resonance. We argue that the data are consistent with only a small fraction of the systems having been captured in resonance. Furthermore, when capture does happen, it is not in general during smooth convergent migration through the disc but after the planets reach the disc inner parts. We show that although resonances may be disrupted when the inner planet enters a central cavity, this alone cannot explain the spread of observed separations. Disruption is found to result in either the system moving interior to the resonance by a few percent, or attaining another resonance. We postulate two populations of low mass planets: a small one for which extensive smooth migration has occurred, and a larger one that formed approximately in-situ with very limited migration.
[3]  oai:arXiv.org:1711.07932  [pdf] - 1670522
The TRAPPIST-1 system: Orbital evolution, tidal dissipation, formation and habitability
Comments: 27 pages, 15 figures, accepted for publication in MNRAS
Submitted: 2017-11-21
We study the dynamical evolution of the TRAPPIST-1 system under the influence of orbital circularization through tidal interaction with the central star. We find that systems with parameters close to the observed one evolve into a state where consecutive planets are linked by first order resonances and consecutive triples, apart from planets c, d and e, by connected three body Laplace resonances. The system expands with period ratios increasing and mean eccentricities decreasing with time. This evolution is largely driven by tides acting on the innermost planets which then influence the outer ones. In order that deviations from commensurability become significant only on $Gy$ time scales or longer, we require that the tidal parameter associated with the planets has to be such that $Q' > \sim 10^{2-3}.$ At the same time, if we start with two subsystems, with the inner three planets comprising the inner one, $Q'$ associated with the planets has to be on the order (and not significantly exceeding) $10^{2-3}$ for the two subsystems to interact and end up in the observed configuration. This scenario is also supported by modelling of the evolution through disk migration which indicates that the whole system cannot have migrated inwards together. Also in order to avoid large departures from commensurabilities, the system cannot have stalled at a disk inner edge for significant time periods. We discuss the habitability consequences of the tidal dissipation implied by our modelling, concluding that planets d, e and f are potentially in habitable zones.
[4]  oai:arXiv.org:1611.05461  [pdf] - 1542803
CoRoT 223992193: Investigating the variability in a low-mass, pre-main sequence eclipsing binary with evidence of a circumbinary disk
Comments: 18 pages, 14 figures, 2 tables. Figures 1 & 2 reproduced from Gillen et al. (2014) for clarity. Accepted for publication in A&A
Submitted: 2016-11-16
CoRoT 223992193 is the only known low-mass, pre-main sequence eclipsing binary that shows evidence of a circumbinary disk. The system displays complex photometric and spectroscopic variability over a range of timescales and wavelengths. Using two optical CoRoT runs, and infrared Spitzer 3.6 and 4.5 $\mu$m observations (simultaneous with the second CoRoT run), we model the out-of-eclipse light curves. The large scale structure in both CoRoT light curves is consistent with the constructive and destructive interference of starspot signals at two slightly different periods. Using the stellar $v\sin i$ 's, we infer different rotation periods: the primary is consistent with synchronisation and the secondary is slightly supersynchronous. Comparison of the raw data to the residuals of our spot model in colour-magnitude space indicates additional contributions consistent with variable dust emission and obscuration. We also identify short-duration flux dips preceding secondary eclipse in all three CoRoT and Spitzer bands. We construct a model of the inner regions of the binary and propose that these dips could be caused by partial occultation of the central binary by the accretion stream onto the primary star. Analysis of 15 VLT/FLAMES H$\alpha$ profiles reveal an emission profile associated with each star: the majority is consistent with chromospheric emission but higher velocity emission is also seen, which could be due to prominences or accretion-related phenomena. In addition, simultaneous u and r-band observations from CFHT/MEGACam reveal a short-lived u-band excess consistent with either an accretion hot spot or stellar flare. The photometric and spectroscopic variations are complex but are consistent with the picture of two active stars possibly undergoing non-steady, low-level accretion; the system's very high inclination provides a new view of such variability. [Abridged]
[5]  oai:arXiv.org:1609.08465  [pdf] - 1494839
On the energy dissipation rate at the inner edge of circumbinary discs
Comments: 14 pages, 11 figures, accepted for publication in MNRAS
Submitted: 2016-09-27
We study, by means of numerical simulations and analysis, the details of the accretion process from a disc onto a binary system. We show that energy is dissipated at the edge of a circumbinary disc and this is associated with the tidal torque that maintains the cavity: angular momentum is transferred from the binary to the disc through the action of compressional shocks and viscous friction. These shocks can be viewed as being produced by fluid elements which drift into the cavity and, before being accreted, are accelerated onto trajectories that send them back to impact the disc. The rate of energy dissipation is approximately equal to the product of potential energy per unit mass at the disc's inner edge and the accretion rate, estimated from the disc parameters just beyond the cavity edge, that would occur without the binary. For very thin discs, the actual accretion rate onto the binary may be significantly less. We calculate the energy emitted by a circumbinary disc taking into account energy dissipation at the inner edge and also irradiation arising there from reprocessing of light from the stars. We find that, for tight PMS binaries, the SED is dominated by emission from the inner edge at wavelengths between 1-4 and 10 $\mu$m. This may apply to systems like CoRoT 223992193 and V1481 Ori.
[6]  oai:arXiv.org:1609.06056  [pdf] - 1483629
On the formation of planetary systems in photoevaporating transition discs
Comments: 8 pages, 5 figures, accepted for publication in MNRAS
Submitted: 2016-09-20
In protoplanetary discs, planetary cores must be at least 0.1 earth mass at 1 au for migration to be significant; this mass rises to 1 earth mass at 5 au. Planet formation models indicate that these cores form on million year timescales. We report here a study of the evolution of 0.1 earth mass and 1 earth mass cores, migrating from about 2 and 5 au respectively, in million year old photoevaporating discs. In such a disc, a gap opens up at around 2 au after a few million years. The inner region subsequently accrete onto the star on a smaller timescale. We find that, typically, the smallest cores form systems of non-resonant planets beyond 0.5 au with masses up to about 1.5 earth mass. In low mass discs, the same cores may evolve in situ. More massive cores form systems of a few earth masses planets. They migrate within the inner edge of the disc gap only in the most massive discs. Delivery of material to the inner parts of the disc ceases with opening of the gap. Interestingly, when the heavy cores do not migrate significantly, the type of systems that are produced resembles our solar system. This study suggests that low mm flux transition discs may not form systems of planets on short orbits but may instead harbour earth mass planets in the habitable zone.
[7]  oai:arXiv.org:1509.08509  [pdf] - 1300364
A circumbinary disc model for the variability of the eclipsing binary CoRoT 223992193
Comments: 17 pages, 6 figures, accepted for publication in MNRAS
Submitted: 2015-09-28
We calculate the flux received from a binary system obscured by a circumbinary disc. The disc is modelled using two dimensional hydrodynamical simulations, and the vertical structure is derived by assuming it is isothermal. The gravitational torque from the binary creates a cavity in the disc's inner parts. If the line of sight along which the system is observed has a high inclination $I$, it intersects the disc and some absorption is produced. As the system is not axisymmetric, the resulting light curve displays variability. We calculate the absorption and produce light curves for different values of the dust disc aspect ratio $H/r$ and mass of dust in the cavity $M_{\rm dust}$. This model is applied to the high inclination ($I=85^{\circ}$) eclipsing binary CoRoT 223992193, which shows 5-10% residual photometric variability after the eclipses and a spot model are subtracted. We find that such variations for $I \sim 85^{\circ}$ can be obtained for $H/r=10^{-3}$ and $M_{\rm dust} \ge 10^{-12}$ M$_{\odot}$. For higher $H/r$, $M_{\rm dust}$ would have to be close to this lower value and $I$ somewhat less than $85^{\circ}$. Our results show that such variability in a system where the stars are at least 90% visible at all phases can be obtained only if absorption is produced by dust located inside the cavity. If absorption is dominated by the parts of the disc located close to or beyond the edge of the cavity, the stars are significantly obscured.
[8]  oai:arXiv.org:1407.7682  [pdf] - 1215978
On the formation of the Kepler-10 planetary system
Comments: 12 pages, accepted for publication in MNRAS
Submitted: 2014-07-29
In this paper, we investigate the conditions required for the 3 and 17 Earth mass solid planets in the Kepler-10 system to have formed through collisions and mergers within an initial population of embryos. By performing a large number of N-body simulations, we show that the total mass of the initial population had to be significantly larger than the masses of the two planets, and that the two planets must have built-up farther away than their present location, at a distance of at least a few au from the central star. The planets had to grow fast enough so that they would detach themselves from the population of remaining, less massive, cores and migrate in to their present location. By the time the other cores migrated in, the disc's inner edge would have moved out so that these cores cannot be detected today. We also compute the critical core mass beyond which a massive gaseous envelope would be accreted and show that it is larger than 17 Earth masses if the planetesimal accretion rate onto the core is larger than 10^{-6} Earth mass per year. For a planetesimal accretion rate between 10^{-6} and 10^{-5} Earth mass per year, the 17 Earth mass core would not be expected to have accreted more than about 1 Earth mass of gas. The results presented in this paper suggest that a planetary system like Kepler-10 may not be unusual, although it has probably formed in a rather massive disc.
[9]  oai:arXiv.org:1406.2189  [pdf] - 1210023
Evolution of eccentricity and orbital inclination of migrating planets in 2:1 mean motion resonance
Comments: 22 pages, 10 figures, accepted for publication in MNRAS
Submitted: 2014-06-09
We determine, analytically and numerically, the conditions needed for a system of two migrating planets trapped in a 2:1 mean motion resonance to enter an inclination-type resonance. We provide an expression for the asymptotic equilibrium value that the eccentricity $e_{\rm i}$ of the inner planet reaches under the combined effects of migration and eccentricity damping. We also show that, for a ratio $q$ of inner to outer masses below unity, $e_{\rm i}$ has to pass through a value $e_{\rm i,res}$ of order 0.3 for the system to enter an inclination-type resonance. Numerically, we confirm that such a resonance may also be excited at another, larger, value $e_{\rm i, res} \simeq 0.6$, as found by previous authors. A necessary condition for onset of an inclination-type resonance is that the asymptotic equilibrium value of $e_{\rm i}$ is larger than $e_{\rm i,res}$. We find that, for $q \le 1$, the system cannot enter an inclination-type resonance if the ratio of eccentricity to semimajor axis damping timescales $t_e/t_a$ is smaller than 0.2. This result still holds if only the eccentricity of the outer planet is damped and $q \lesssim 1$. As the disc/planet interaction is characterized by $t_e/t_a \sim 10^{-2}$, we conclude that excitation of inclination through the type of resonance described here is very unlikely to happen in a system of two planets migrating in a disc.
[10]  oai:arXiv.org:1311.3990  [pdf] - 1180758
CoRoT 223992193: A new, low-mass, pre-main sequence eclipsing binary with evidence of a circumbinary disk
Comments: 19 pages, 14 figures, 5 tables. Accepted for publication in A&A
Submitted: 2013-11-15
We present the discovery of CoRoT 223992193, a double-lined, detached eclipsing binary, comprising two pre-main sequence M dwarfs, discovered by the CoRoT space mission during a 23-day observation of the 3 Myr old NGC 2264 star-forming region. Using multi-epoch optical and near-IR follow-up spectroscopy with FLAMES on the Very Large Telescope and ISIS on the William Herschel Telescope we obtain a full orbital solution and derive the fundamental parameters of both stars by modelling the light curve and radial velocity data. The orbit is circular and has a period of $3.8745745 \pm 0.0000014$ days. The masses and radii of the two stars are $0.67 \pm 0.01$ and $0.495 \pm 0.007$ $M_{\odot}$ and $1.30 \pm 0.04$ and $1.11 ~^{+0.04}_{-0.05}$ $R_{\odot}$, respectively. This system is a useful test of evolutionary models of young low-mass stars, as it lies in a region of parameter space where observational constraints are scarce; comparison with these models indicates an apparent age of $\sim$3.5-6 Myr. The systemic velocity is within $1\sigma$ of the cluster value which, along with the presence of lithium absorption, strongly indicates cluster membership. The CoRoT light curve also contains large-amplitude, rapidly evolving out-of-eclipse variations, which are difficult to explain using starspots alone. The system's spectral energy distribution reveals a mid-infrared excess, which we model as thermal emission from a small amount of dust located in the inner cavity of a circumbinary disk. In turn, this opens up the possibility that some of the out-of-eclipse variability could be due to occultations of the central stars by material located at the inner edge or in the central cavity of the circumbinary disk.
[11]  oai:arXiv.org:1309.1025  [pdf] - 1178995
The effects of disc warping on the inclination of planetary orbits
Comments: 25 pages, 4 figures, accepted for publication in MNRAS
Submitted: 2013-09-04
The interaction between a planet located in the inner region of a disc and the warped outer region is studied. We consider the stage of evolution after the planet has cleared-out a gap, so that the planetary orbit evolves only under the gravitational potential from the disc. We develop a secular analysis and compute the evolution of the orbital elements by solving Lagrange's equations valid to second order in the eccentricity. We also perform numerical simulations with the full disc potential. In general, the interaction between the disc and the planet leads to the precession of the orbit. The orbital plane therefore becomes tilted relative to the disc's inner parts, with no change in the eccentricity. When the inclination approaches 90 degrees, there is an instability and the eccentricity increases. In this case, both the inclination and the eccentricity develop large variations, with the orbit becoming retrograde. As the eccentricity reaches high values, we would expect tidal capture on a short orbit of the planet by the star to occur. This instability happens when the disc is severely warped, or if there is a significant amount of mass in a ring inclined by at least 45 degrees relative to the initial orbital plane. The inclination of the orbit does not depend on the semimajor axis nor on the planet's mass. However, for a significant inclination to be generated on a timescale of at most a few Myr, the planet should be beyond the snow line. The process described here would therefore produce two distinct populations of inclined planets: one with objects beyond the snow line with at most moderate eccentricities, and another with objects on short circularized orbits.
[12]  oai:arXiv.org:1209.4843  [pdf] - 1151561
Orbital evolution of a planet on an inclined orbit interacting with a disc
Comments: 28 pages, 8 figures, accepted for publication in MNRAS
Submitted: 2012-09-21
We study the dynamics of a planet on an orbit inclined with respect to a disc. If the initial inclination of the orbit is larger than some critical value, the gravitational force exerted by the disc on the planet leads to a Kozai cycle in which the eccentricity of the orbit is pumped up to large values and oscillates with time in antiphase with the inclination. On the other hand, both the inclination and the eccentricity are damped by the frictional force that the planet is subject to when it crosses the disc. We show that, by maintaining either the inclination or the eccentricity at large values, the Kozai effect provides a way of delaying alignment with the disc and circularization of the orbit. We find the critical value to be characteristically as small as about 20 degrees. Typically, Neptune or lower mass planets would remain on inclined and eccentric orbits over the disc lifetime, whereas orbits of Jupiter or higher mass planets would align and circularize. This could play a significant role in planet formation scenarios.
[13]  oai:arXiv.org:1109.1704  [pdf] - 1083922
Protoplanets with core masses below the critical mass fill in their Roche lobe
Comments: 17 pages, 1 figure, accepted by MNRAS
Submitted: 2011-09-08
We study the evolution of a protoplanet of a few earth masses embedded in a protoplanetary disc. If we assume that the atmosphere of the protoplanet, i.e. the volume of gas in hydrostatic equilibrium bound to the core, has a surface radius smaller than the Roche lobe radius, we show that it expands as it accretes both planetesimals and gas at a fixed rate from the nebula until it fills in the Roche lobe. The evolution occurs on a timescale shorter than the formation or migration timescales. Therefore, we conclude that protoplanets of a few earth masses have an atmosphere that extends to the Roche lobe surface, where it matches to the nebula. This is true even when the Bondi radius is smaller than the Roche lobe radius. This is in contrast to the commonly used models in which the static atmosphere extends up to the Bondi radius and is surrounded by a cold accretion flow. As a result, any calculation of the tidal torque exerted by the disc onto the protoplanet should exclude the material present in the Roche lobe, since it is bound to the protoplanet.
[14]  oai:arXiv.org:1002.2125  [pdf] - 1025107
On the dynamics of multiple systems of hot super-Earths and Neptunes: Tidal circularization, resonance and the HD 40307 system
Comments: 38 pages, 12 figures, accepted for publication in MNRAS
Submitted: 2010-02-10
[Abridged] We consider the dynamics of a system of hot super-Earths or Neptunes such as HD 40307. We show that, as tidal interaction with the central star leads to small eccentricities, the planets in this system could be undergoing resonant coupling even though the period ratios depart significantly from very precise commensurability. In a three planet system, this is indicated by the fact that resonant angles librate or are associated with long term changes to the orbital elements. We propose that the planets in HD 40307 were in a strict Laplace resonance while they migrated through the disc. After entering the disc inner cavity, tidal interaction would cause the period ratios to increase from two but with the inner pair deviating less than the outer pair, counter to what occurs in HD 40307. However, the relationship between these pairs that occurs in HD 40307 might be produced if the resonance is impulsively modified by an event like a close encounter shortly after the planetary system decouples from the disc. We find this to be in principle possible for a small relative perturbation on the order of a few 1.d-3 but then we find that evolution to the present system in a reasonable time is possible only if the masses are significantly larger than the minimum masses and tidal dissipation is very effective. On the other hand we found that a system like HD 40307 with minimum masses and more realistic tidal dissipation could be produced if the eccentricity of the outermost planet was impulsively increased to about 0.15. The form of resonantly coupled tidal evolution we consider here is quite general and could be of greater significance for systems with inner planets on significantly shorter orbital periods characteristic of for example CoRoT 7 b.
[15]  oai:arXiv.org:1001.0657  [pdf] - 1018998
Eccentricity pumping of a planet on an inclined orbit by a disc
Comments: 13 pages, 4 figures, accepted for publication in MNRAS
Submitted: 2010-01-05
In this paper, we show that the eccentricity of a planet on an inclined orbit with respect to a disc can be pumped up to high values by the gravitational potential of the disc, even when the orbit of the planet crosses the disc plane. This process is an extension of the Kozai effect. If the orbit of the planet is well inside the disc inner cavity, the process is formally identical to the classical Kozai effect. If the planet's orbit crosses the disc but most of the disc mass is beyond the orbit, the eccentricity of the planet grows when the initial angle between the orbit and the disc is larger than some critical value which may be significantly smaller than the classical value of 39 degrees. Both the eccentricity and the inclination angle then vary periodically with time. When the period of the oscillations of the eccentricity is smaller than the disc lifetime, the planet may be left on an eccentric orbit as the disc dissipates.
[16]  oai:arXiv.org:0808.3897  [pdf] - 15750
New composite models of partially ionized protoplanetary disks
Comments: 13 pages, 5 figures, accepted for publication in ApJ
Submitted: 2008-08-28
We study an accretion disk in which three different regions may coexist: MHD turbulent regions, dead zones and gravitationally unstable regions. Although the dead zones are stable, there is some transport due to the Reynolds stress associated with waves emitted from the turbulent layers. We model the transport in each of the different regions by its own $\alpha$ parameter, this being 10 to $10^{3}$ times smaller in dead zones than in active layers. In gravitationally unstable regions, $\alpha$ is determined by the fact that the disk self-adjusts to a state of marginal stability. We construct steady-state models of such disks. We find that for uniform mass flow, the disk has to be more massive, hotter and thicker at the radii where there is a dead zone. In disks in which the dead zone is very massive, gravitational instabilities are present. Whether such models are realistic or not depends on whether hydrodynamical fluctuations driven by the turbulent layers can penetrate all the way inside the dead zone. This may be more easily achieved when the ratio of the mass of the active layer to that of the dead zone is relatively large, which in our models corresponds to $\alpha$ in the dead zone being about 10% of $\alpha$ in the active layers. If the disk is at some stage of its evolution not in steady-state, then the surface density will evolve toward the steady-state solution. However, if $\alpha$ in the dead zone is much smaller than in the active zone, the timescale for the parts of the disk beyond a few AU to reach steady-state may become longer than the disk lifetime. Steady-state disks with dead zones are a more favorable environment for planet formation than standard disks, since the dead zone is typically 10 times more massive than a corresponding turbulent zone at the same location.
[17]  oai:arXiv.org:astro-ph/0609779  [pdf] - 85371
Migration and the formation of systems of hot super-Earths and Neptunes
Comments: 29 pages, to be published in ApJ
Submitted: 2006-09-28
The existence of extrasolar planets with short orbital periods suggests that planetary migration induced by tidal interaction with the protoplanetary disk is important. Cores and terrestrial planets may undergo migration as they form. In this paper we investigate the evolution of a population of cores with initial masses in the range 0.1-1 earth mass embedded in a disk. Mutual interactions lead to orbit crossing and mergers, so that the cores grow during their evolution. Interaction with the disk leads to orbital migration, which results in the cores capturing each other in mean motion resonances. As the cores migrate inside the disk inner edge, scatterings and mergers of planets on unstable orbits together with orbital circularization causes strict commensurability to be lost. Near commensurability however is usually maintained. All the simulations end with a population of typically between two and five planets, with masses depending on the initial mass. These results indicate that if hot super-Earths or Neptunes form by mergers of inwardly migrating cores, then such planets are most likely not isolated. We would expect to always find at least one, more likely a few, companions on close and often near-commensurable orbits. To test this hypothesis, it would be of interest to look for planets of a few to about 10 earth masses in systems where hot super-Earths or Neptunes have already been found.
[18]  oai:arXiv.org:astro-ph/0510487  [pdf] - 76983
Planet formation and migration
Comments: Revised version. Typos corrected in equations (36)-(48). Article available free at www.iop.org/journals/thismonth until December 8
Submitted: 2005-10-17, last modified: 2005-11-28
We review the observations of extrasolar planets, ongoing developments in theories of planet formation, orbital migration, and the evolution of multiplanet systems.
[19]  oai:arXiv.org:astro-ph/0508671  [pdf] - 75560
Numerical simulations of type I planetary migration in nonturbulent magnetized discs
Comments: 17 pages, 12 figures, accepted in MNRAS. A version with full resolution, colour figures is available at http://www.maths.qmul.ac.uk/~rpn/preprints/
Submitted: 2005-08-31
Using 2D MHD numerical simulations performed with two different finite difference Eulerian codes, we analyze the effect that a toroidal magnetic field has on low mass planet migration in nonturbulent protoplanetary discs. The presence of the magnetic field modifies the waves that can propagate in the disc. In agreement with a recent linear analysis (Terquem 2003), we find that two magnetic resonances develop on both sides of the planet orbit, which contribute to a significant global torque. In order to measure the torque exerted by the disc on the planet, we perform simulations in which the latter is either fixed on a circular orbit or allowed to migrate. For a 5 earth mass planet, when the ratio \beta between the square of the sound speed and that of the Alfven speed at the location of the planet is equal to 2, we find inward migration when the magnetic field B_{\phi} is uniform in the disc, reduced migration when B_{\phi} decreases as r^{-1} and outward migration when B_{\phi} decreases as r^{-2}. These results are in agreement with predictions from the linear analysis. Taken as a whole, our results confirm that even a subthermal stable field can stop inward migration of an earth--like planet.
[20]  oai:arXiv.org:astro-ph/0409404  [pdf] - 67485
Evolution of self-gravitating magnetized disks. II- Interaction between MHD turbulence and gravitational instabilities
Comments: 31 pages, 19 figures, accepted for publication in ApJ, animation avalaible at http://www2.iap.fr/users/fromang/simu3d/simu3d.html
Submitted: 2004-09-16
We present 3D magnetohydrodynamic (MHD) numerical simulations of the evolution of self--gravitating and weakly magnetized disks with an adiabatic equation of state. Such disks are subject to the development of both the magnetorotational and gravitational instabilities, which transport angular momentum outward. As in previous studies, our hydrodynamical simulations show the growth of strong m=2 spiral structure. This spiral disturbance drives matter toward the central object and disappears when the Toomre parameter Q has increased well above unity. When a weak magnetic field is present as well, the magnetorotational instability grows and leads to turbulence. In that case, the strength of the gravitational stress tensor is lowered by a factor of about~2 compared to the hydrodynamical run and oscillates periodically, reaching very small values at its minimum. We attribute this behavior to the presence of a second spiral mode with higher pattern speed than the one which dominates in the hydrodynamical simulations. It is apparently excited by the high frequency motions associated with MHD turbulence. The nonlinear coupling between these two spiral modes gives rise to a stress tensor that oscillates with a frequency which is a combination of the frequencies of each of the modes. This interaction between MHD turbulence and gravitational instabilities therefore results in a smaller mass accretion rate onto the central object.
[21]  oai:arXiv.org:astro-ph/0402373  [pdf] - 62887
Evolution of massive and magnetized protoplanetary disks
Comments: 8 pages, 7 figures, to appear in the proceedings of the XIXth IAP colloquium "Extrasolar Planets: Today and Tomorrow" held in Paris, France, 2003, June 30 - July 4, ASP Conf. Series
Submitted: 2004-02-16
We present global 2D and 3D simulations of self-gravitating magnetized tori. We used the 2D calculations to demonstrate that the properties of the MRI are not affected by the presence of self-gravity: MHD turbulence and enhanced angular momentum transport follow the linear growth of the instability. In 3D, we have studied the interaction between an $m=2$ gravitational instability and MHD turbulence. We found its strength to be significantly decreased by the presence of the latter, showing that both instabilities strongly interact in their non-linear phases. We discuss the consequences of these results.
[22]  oai:arXiv.org:astro-ph/0309175  [pdf] - 59002
Is planetary migration inevitable?
Comments: 8 pages, submitted for publication in 'Planetary Systems and Planets in Systems', eds. S. Udry, W. Benz and R. von Steiger (ISSI Space Sciences Series, 19, reprinted from Space Science Reviews - Kluwer)
Submitted: 2003-09-05
According to current theories, tidal interactions between a disk and an embedded planet may lead to the rapid migration of the protoplanet on a timescale shorter than the disk lifetime or estimated planetary formation timescales. Therefore, planets can form only if there is a mechanism to hold at least some of the cores back on their way in. Once a giant planet has assembled, there also has to be a mechanism to prevent it from migrating down to the disk center. This paper reviews the different mechanisms that have been proposed to stop or slow down migration.
[23]  oai:arXiv.org:astro-ph/0301556  [pdf] - 54527
Stopping inward planetary migration by a toroidal magnetic field
Comments: 25 pages including 10 figures, Latex in the MN style - Accepted by MNRAS
Submitted: 2003-01-28
We calculate the linear torque exerted by a planet on a circular orbit on a disc containing a toroidal magnetic field. All fluid perturbations are singular at the so--called magnetic resonances, where the Doppler shifted frequency of the perturbation matches that of a slow MHD wave propagating along the field line. These lie on both sides of the corotation radius. Waves propagate outside the Lindblad resonances, and also in a restricted region around the magnetic resonances. The magnetic resonances contribute to a significant global torque which, like the Lindblad torque, is negative (positive) inside (outside) the planet's orbit. Since these resonances are closer to the planet than the Lindblad resonances, the torque they contribute dominates over the Lindblad torque if the magnetic field is large enough. In addition, if beta=c^2/v_A^2 increases fast enough with radius, the outer magnetic resonance becomes less important and the total torque is then negative, dominated by the inner magnetic resonance. This leads to outward migration of the planet. Even for beta=100 at corotation, a negative torque may be obtained. A planet migrating inward through a nonmagnetized region of a disc would then stall when reaching a magnetized region. It would then be able to grow to become a terrestrial planet or the core of a giant planet. In a turbulent magnetized disc in which the large scale field structure changes sufficiently slowly, a planet may alternate between inward and outward migration, depending on the gradients of the field encountered. Its migration could then become diffusive, or be limited only to small scales.
[24]  oai:arXiv.org:astro-ph/0203187  [pdf] - 48208
Dynamical relaxation and the orbits of low-mass extrasolar planets
Comments: 6 pages including 6 figures, Latex in the MN style - Accepted by MNRAS Letters
Submitted: 2002-03-12
We consider the evolution of a system containing a population of massive planets formed rapidly through a fragmentation process occurring on a scale on the order of 100 au and a lower mass planet that assembles in a disc on a much longer timescale. During the formation phase, the inner planet is kept on a circular orbit due to tidal interaction with the disc, while the outer planets undergo dynamical relaxation. Interaction with the massive planets left in the system after the inner planet forms may increase the eccentricity of the inner orbit to high values, producing systems similar to those observed.
[25]  oai:arXiv.org:astro-ph/0109008  [pdf] - 44507
The ionization fraction in alpha-models of protoplanetary disks
Comments: 25 pages including 8 figures, Latex in the MN style - Accepted by MNRAS
Submitted: 2001-09-02
We calculate the ionization fraction of protostellar alpha disks, taking into account vertical temperature structure, and the possible presence of trace metal atoms. Both thermal and X-ray ionization are considered. Previous investigations of layered disks used radial power-law models with isothermal vertical structure. But alpha models are used to model accretion, and the present work is a step towards a self-consistent treatment. The extent of the magnetically uncoupled (``dead'') zone depends sensitively on alpha, on the assumed accretion rate, and on the critical magnetic Reynolds number, below which MHD turbulence cannot be self-sustained. Its extent is extremely model-dependent. It is also shown that a tiny fraction of the cosmic abundance of metal atoms can dramatically affect the ionization balance. Gravitational instabilities are an unpromising source of transport, except in the early stages of disk formation.
[26]  oai:arXiv.org:astro-ph/0107408  [pdf] - 1232714
Theory of Turbulent Accretion Disks
Comments: 25 pages, LaTex, To appear in the proceedings of the Aussois 2000 summer school ``Formation Stellaire et Physique des Etoiles Jeunes''
Submitted: 2001-07-20
In low-mass disks, turbulent torques are probably the most important way of redistributing angular momentum. Here we present the theory of turbulent accretion disks. We show the molecular viscosity is far too small to account for the evolutionary timescale of disks, and we describe how turbulence may result in enhanced transport of (angular) momentum. We then turn to the magnetorotational instability, which thus far is the only mechanism that has been shown to initiate and sustain turbulence in disks. Finally, we present both the basis and the structure of alpha disk models.
[27]  oai:arXiv.org:astro-ph/0012391  [pdf] - 39980
Dynamical Relaxation and Massive Extrasolar Planets
Comments: Submitted to MNRAS on December 18, 2000
Submitted: 2000-12-18
Following the suggestion of Black (1997) that some massive extrasolar planets may be associated with the tail of the distribution of stellar companions, we investigate a scenario in which 5 < N < 100 planetary mass objects are assumed to form rapidly through a fragmentation process occuring in a disc or protostellar envelope on a scale of 100 au. These are assumed to have formed rapidly enough through gravitational instability or fragmentation that their orbits can undergo dynamical relaxation on a timescale of about 100 orbits. Under a wide range of initial conditions and assumptions the relaxation process ends with either (i) one potential 'hot Jupiter' plus up to two 'external' companions, i.e. planets orbiting near the outer edge of the initial distribution; (ii) one or two 'external' planets or even none at all; (iii) one planet on an orbit with a semi--major axis 10 to a 100 times smaller than the outer boundary radius of the inital distribution together with an 'external' companion. Most of the other objects are ejected and could contribute to a population of free floating planets. Apart from the potential 'hot Jupiters', all the bound objects are on highly eccentric orbits. We found that, apart from the close orbiters, the probability of ending up with a planet orbiting at a given distance from the central star increases with the distance. This is because of the tendency of the relaxation process to lead to collisions with the central star. We discuss the application of these results to some of the more massive extrasolar planets.
[28]  oai:arXiv.org:astro-ph/0010229  [pdf] - 38576
Linear Analysis of the Hall Effect in Protostellar Disks
Comments: 34 pages, 10 figures, AAS LaTEx, v.4.0. Submitted to ApJ
Submitted: 2000-10-11
The effects of Hall electromotive forces (HEMFs) on the linear stability of protostellar disks are examined. Earlier work on this topic focused on axial field and perturbation wavenumbers. Here we treat the problem more generally. Both axisymmetric and nonaxisymmetric cases are investigated. Though seldom explicitly included in calculations, HEMFs appear to be important whenever Ohmic dissipation is. They allow for the appearance of electron whistler waves, and since these have right-handed polarization, a helicity factor is also introduced into the stability problem. This factor is the product of the components of the angular velocity and magnetic field along the perturbation wavenumber, and it is destabilizing when negative. Unless the field and angular velocity are exactly aligned, it is always possible to find destabilizing wavenumbers. HEMFs can destabilize any differential rotation law, even those with angular velocity increasing outward. Regardless of the sign of the angular velocity gradient, the maximum growth rate is always given in magnitude by the local Oort A value of the disk, as in the standard magnetorotational instability. The role of Hall EMFs may prove crucial to understanding how turbulence is maintained in the ``low state'' of eruptive disk systems.
[29]  oai:arXiv.org:astro-ph/0008514  [pdf] - 37830
Tidally-induced angular momentum transport in disks
Comments: 4 pages, LaTeX, uses newpasp.sty, to be published in 'The Formation of Binary Stars', eds. R.D. Mathieu and H. Zinnecker (ASP Conference Series)
Submitted: 2000-08-31
We discuss the transport of angular momentum induced by tidal effects in a disk surrounding a star in a pre-main sequence binary system. We consider the effect of both density and bending waves. Although tidal effects are important for truncating protostellar disks and for determining their size, it is unlikely that tidally-induced angular momentum transport plays a dominant role in the evolution of protostellar disks. Where the disk is magnetized, transport of angular momentum is probably governed by MHD turbulence. In a non self-gravitating laminar disk, the amount of transport provided by tidal waves is probably too small to account for the lifetime of protostellar disks. In addition, tidal effects tend to be localized in the disk outer regions.
[30]  oai:arXiv.org:astro-ph/0006113  [pdf] - 36493
The response of an accretion disc to an inclined dipole with application to AA Tau
Comments: 18 pages (including 7 figures), Latex, to be published in A&A, also avalaible at http://www.iap.fr/users/terquem/index.html
Submitted: 2000-06-08
We compute the warping of a disc induced by an inclined dipole. We consider a magnetised star surrounded by a thin Keplerian diamagnetic disc with an inner edge that corotates with the star. We suppose the stellar field is a dipole with an axis that is slightly misaligned with the stellar rotation axis. The rotation axes of the disc material orbiting at large distances from the star and that of the star are supposed to coincide. The misalignment of the magnetic and rotation axes results in the magnetic pressure not being the same on the upper and lower surfaces of the disc. The resultant net vertical force produces a warp which appears stationary in a frame corotating with the star. We find that, if viscosity is large enough (alpha on the order of 0.01-0.1) to damp bending waves as they propagate away, a smoothly varying warp of the inner region of the disc is produced. The amplitude of the warp can easily be on the order of ten percent of the disc inner radius for reasonably small misalignment angles (less than 30 degrees). Viscous damping also introduces a phase shift between the warp and the forcing torque, which results in the locations of maximum elevation above the disc forming a trailing spiral pattern. We apply these results to recent observations of AA Tau, and show that the variability of its light curve, which occurs with a period comparable to the expected stellar rotation period, could be due to obscuration produced by a warp configuration of the type we obtain.
[31]  oai:arXiv.org:astro-ph/9905229  [pdf] - 106555
Disks, Extrasolar Planets and Migration
Comments: 22 pages (including 3 figures), LaTeX, uses kluwer.cls, to be published in 'From Dust to Terrestrial Planets', eds W. Benz, R. Kallenbach, G. Lugmair and F. Podosek, ISSI Space Sciences Series No. 9, reprinted in Space Science Reviews (January 1, 2000)
Submitted: 1999-05-18
We review results about protoplanetary disk models, protoplanet migration and formation of giant planets with migrating cores. We first model the protoplanetary nebula as an \alpha-accretion disk and present steady state calculations for different values of \alpha and gas accretion rate through the disk. We then review the current theories of protoplanet migration in the context of these models, focusing on the gaseous disk-protoplanet tidal interaction. According to these theories, the migration timescale may be shorter than the planetary formation timescale. Therefore we investigate planet formation in the context of a migrating core, considering both the growth of the core and the build-up of the envelope in the course of the migration.
[32]  oai:arXiv.org:astro-ph/9903310  [pdf] - 105714
Critical Protoplanetary Core Masses in Protoplanetary Disks and the Formation of Short-Period Giant Planets
Comments: 24 pages (including 9 figures), LaTeX, uses emulateapj.sty, to be published in ApJ, also available at http://www.ucolick.org/~ct/home.html
Submitted: 1999-03-19
We study a solid protoplanetary core of 1-10 earth masses migrating through a disk. We suppose the core luminosity is generated as a result of planetesimal accretion and calculate the structure of the gaseous envelope assuming equilibrium. This is a good approximation when the core mass is less than the critical value, M_{crit}, above which rapid gas accretion begins. We model the structure of the protoplanetary nebula as an accretion disk with constant \alpha. We present analytic fits for the steady state relation between disk surface density and mass accretion rate as a function of radius r. We calculate M_{crit} as a function of r, gas accretion rate through the disk, and planetesimal accretion rate onto the core \dot{M}. For a fixed \dot{M}, M_{crit} increases inwards, and it decreases with \dot{M}. We find that \dot{M} onto cores migrating inwards in a time 10^3-10^5 yr at 1 AU is sufficient to prevent the attainment of M_{crit} during the migration process. Only at small radii where planetesimals no longer exist can M_{crit} be attained. At small radii, the runaway gas accretion phase may become longer than the disk lifetime if the core mass is too small. However, massive cores can be built-up through the merger of additional incoming cores on a timescale shorter than for in situ formation. Therefore, feeding zone depletion in the neighborhood of a fixed orbit may be avoided. Accordingly, we suggest that giant planets may begin to form early in the life of the protostellar disk at small radii, on a timescale that may be significantly shorter than for in situ formation. (abridged)
[33]  oai:arXiv.org:astro-ph/9812155  [pdf] - 104262
Precession of collimated outflows from young stellar objects
Comments: 5 pages (including 1 figure), LaTeX, uses emulateapj.sty, to be published in ApJ Letters, also available at http://www.ucolick.org/~ct/home.html and http://www.tls-tautenburg.de/research/research.html
Submitted: 1998-12-08
We consider several protostellar systems where either a precessing jet or at least two misaligned jets have been observed. We assume that the precession of jets is caused by tidal interactions in noncoplanar binary systems. For Cep E, V1331 Cyg and RNO 15-FIR the inferred orbital separations and disk radii are in the range 4-160 AU and 1-80 AU, respectively, consistent with those expected for pre-main sequence stars. Furthermore, we assume or use the fact that the source of misaligned outflows is a binary, and evaluate the lengthscale over which the jets should precess as a result of tidal interactions. For T Tau, HH1 VLA 1/2 and HH 24 SVS63, it may be possible to detect a bending of the jets rather than 'wiggling'. In HH 111 IRS and L1551 IRS5, 'wiggling' may be detected on the current observed scale. Our results are consistent with the existence of noncoplanar binary systems in which tidal interactions induce jets to precess.
[34]  oai:arXiv.org:astro-ph/9810172  [pdf] - 103288
The response of accretion disks to bending waves: angular momentum transport and resonances
Comments: 23 pages (including 7 figures), LaTeX, uses emulateapj.sty, to be published in ApJ, also available at http://www.ucolick.org/~ct/home.html
Submitted: 1998-10-11
We investigate the linear tidal perturbation of a viscous Keplerian disk by a companion star orbiting in a plane inclined to the disk. We consider m=1 perturbations with odd symmetry with respect to the z=0 midplane. Since the response of a viscous disk is not in phase with the perturbing potential, a tidal torque is exerted on the disk, resulting in a decrease of its angular momentum. This tidal torque is found to be comparable to the horizontal viscous stress acting on the background flow when the perturbed velocities in the disk are on the order of the sound speed. If these velocities remain subsonic, the tidal torque can exceed the horizontal viscous stress only if the viscous parameter \alpha which couples to the vertical shear is larger than that coupled to the horizontal shear. In protostellar disks, bending waves are found to propagate deep into the disk inner parts. If the waves are reflected at the center, resonances occur when the frequency of the tidal waves is equal to that of some free normal global bending mode of the disk. If such resonances exist, tidal interactions may then be important even when the binary separation is large. Out of resonance, the torque associated with the secular perturbation is generally much larger than that associated with the finite frequency perturbations. As long as the waves are damped before they reach the center, the torque associated with the finite frequency perturbations does not depend on the viscosity. These calculations are relevant to disks around young stars and maybe also to disks in X-ray binary systems. (abridged)
[35]  oai:arXiv.org:astro-ph/9810014  [pdf] - 103130
Tidally-induced warps in protostellar discs
Comments: 10 pages, LaTeX (including 1 figure), uses paspconf.sty and epsf.sty, to be published in Proceedings of the EC Summer School on 'Astrophysical Discs', eds J. A. Sellwood and J. Goodman, ASP Conf. Series
Submitted: 1998-10-01
We review results on the dynamics of warped gaseous discs. We consider tidal perturbation of a Keplerian disc by a companion star orbiting in a plane inclined to the disc. The perturbation induces the precession of the disc, and thus of any jet it could drive. In some conditions the precession rate is uniform, and as a result the disc settles into a warp mode. The tidal torque also leads to the truncation of the disc, to the evolution of the inclination angle (not necessarily towards alignment of the disc and orbital planes) and to a transport of angular momentum in the disc. We note that the spectral energy distribution of such a warped disc is different from that of a flat disc. We conclude by listing observational effects of warps in protostellar discs.
[36]  oai:arXiv.org:astro-ph/9810027  [pdf] - 103143
Discs and Planetary Formation
Comments: 22 pages, LaTeX (including 6 figures), uses paspconf.sty, epsf.sty and rotate.sty, to be published in Proceedings of the EC Summer School on 'Astrophysical Discs', eds J. A. Sellwood and J. Goodman, ASP Conf. Series
Submitted: 1998-10-01
The formation, structure and evolution of protoplanetary discs is considered. The formation of giant planets within the environment of these models is also discussed.
[37]  oai:arXiv.org:astro-ph/9809200  [pdf] - 102904
Planet formation, orbital evolution and planet-star tidal interaction
Comments: 26 pages, TeX, 3 figures (2 in color), 2 style files (ppiv-style.tex and epsf.sty), to be published in "Protostars and Planets IV," eds. A. Boss, V. Mannings, and S. Russell
Submitted: 1998-09-15
We consider several processes operating during the late stages of planet formation that can affect observed orbital elements. Disk-planet interactions, tidal interactions with the central star, long term orbital instability and the Kozai mechanism are discussed.
[38]  oai:arXiv.org:astro-ph/9801280  [pdf] - 100129
On the tidal interaction of a solar-type star with an orbiting companion: Excitation of g mode oscillation and orbital evolution
Comments: 36 pages including 6 PostScript figures, LaTex -- To be published in ApJ -- Also available at http://www.ucolick.org/~ct/home.html
Submitted: 1998-01-28
We calculate the dynamical tides raised on a non-rotating solar-type star by a close stellar or planetary companion. Dissipation arising from a turbulent viscosity operating in the convection zone and radiative damping in the radiative core are considered. We compute the torque exerted on the star by a companion in circular orbit, and determine the potentially observable magnitude of the tidally induced velocity at the stellar photosphere. These calculations are compared with the results obtained by assuming that a very small frequency limit can be taken in order to calculate the tidal response (equilibrium tide). For a standard solar model, the latter is found to give a relatively poor approximation at the periods of interest of several days, even when the system is far from resonance with a normal mode. It is shown that although the companion may go through a succession of resonances as it spirals in under the action of the tides, for a fixed spectrum of normal modes its migration is controlled essentially by the non-resonant interaction. We find that the turbulent viscosity that is required to provide the observed circularization rates of main sequence solar-type binaries is about fifty times larger than that simply estimated from mixing length theory for non-rotating stars. We discuss the means by which this enhanced viscosity might be realized. These calculations are applied to 51 Pegasi. We show that the perturbed velocity induced by the tides at the stellar surface is too small to be observed.
[39]  oai:arXiv.org:astro-ph/9711213  [pdf] - 99346
Oscillations in solar-type stars tidally induced by orbiting planets
Comments: 5 pages with 1 Postscript figure, LaTeX. To be published in the Proceedings of the IXemes Rencontres de Blois "Planetary Systems: The Long View" (Blois, June 1997), Editions Frontieres
Submitted: 1997-11-18
We examine the effect of dynamical tides raised by a companion on a solar-type star. In these binaries, gravity or g mode oscillations are excited by the companion in the radiative region beneath the convective envelope of the star. They become evanescent in the convection zone. This is of particular interest in connection with the newly discovered planets, some of which are found to orbit around solar-type stars with a period comparable to that of the high order g modes of the star. One such example is 51 Pegasi. Here, we determine the magnitude of the perturbed velocity induced by the tides at the stellar surface. We show that, in the case of 51 Pegasi, this velocity is too small to be observed. This result is insensitive to the magnitude of the stellar turbulent viscosity assumed and is not affected by the possibility of resonance, which occurs when the frequency of the tidal disturbance is close to that of some normal mode of the star. We also discuss the orbital evolution and synchronization timescales associated with the tidal interaction.
[40]  oai:arXiv.org:astro-ph/9710047  [pdf] - 1469565
On the global warping of a thin self-gravitating near Keplerian gaseous disk with application to the disk in NGC 4258
Comments: 36 pages (including 4 figures), Latex, to appear in ApJ
Submitted: 1997-10-04
On the global warping of a thin self-gravitating near Keplerian gaseous disk with application to the disk in NGC 4258
[41]  oai:arXiv.org:astro-ph/9707316  [pdf] - 98140
Bending Instabilities in Magnetized Accretion Discs
Comments: Accepted by MNRAS. 18 pages, 11 figures. LaTeX2e in the MN style. PostScript and HTML files are also available from http://www-star.qmw.ac.uk/~va/ or by e-mail: V.Agapitou@qmw.ac.uk
Submitted: 1997-07-29
We study the global bending modes of a thin annular disc subject to both an internally generated magnetic field and a magnetic field due to a dipole embedded in the central star with axis aligned with the disc rotation axis. When there is a significant inner region of the disc corotating with the star, we find spectra of unstable bending modes. These may lead to elevation of the disc above the original symmetry plane facilitating accretion along the magnetospheric field lines. The resulting non-axisymmetric disc configuration may result in the creation of hot spots on the stellar surface and the periodic photometric variations observed in many classical T Tauri stars (CTTS). Time-dependent behaviour may occur including the shadowing of the central source in magnetic accretors even when the dipole and rotation axes are aligned.
[42]  oai:arXiv.org:astro-ph/9701106  [pdf] - 96379
Precessing warped discs in close binary systems
Comments: 14 pages (LaTeX, lamuphys.sty), 6 figures available from J.D.Larwood@qmw.ac.uk. Full PostScript available from http://www.maths.qmw.ac.uk/~jdl/preprint.html . To appear in the Proceedings of the EARA Workshop on Accretion Disks (Garching, Oct. 1996), Lecture Notes in Physics
Submitted: 1997-01-16
We describe some recent nonlinear three dimensional hydrodynamic simulations of accretion discs in binary systems where the orbit is circular and not necessarily coplanar with the disc midplane. The calculations are relevant to a number of observed astrophysical phenomena, including the precession of jets associated with young stars, the high spectral index of some T Tauri stars, and the light curves of X-ray binaries such as Hercules X-1 which suggest the presence of precessing accretion discs.
[43]  oai:arXiv.org:astro-ph/9604013  [pdf] - 1234627
The tidally induced warping, precession and truncation of accretion discs in binary systems: three-dimensional simulations
Comments: Accepted by MNRAS, received in original form 23 January 1996. LaTeX, 22 pages. 15 figures available from J.D.Larwood@qmw.ac.uk, or by anonymous ftp (see http://www.maths.qmw.ac.uk/~jdl/preprint.html for instructions)
Submitted: 1996-04-02, last modified: 1996-07-31
We present the results of non linear, hydrodynamic simulations, in three dimensions, of the tidal perturbation of accretion discs in binary systems where the orbit is circular and not necessarily coplanar with the disc mid-plane. The accretion discs are assumed to be geometrically thin, and of low mass relative to the stellar mass so that they are governed by thermal pressure and viscosity, but not self-gravity. The parameters that we consider in our models are the ratio of the orbital distance to the disc radius, D/R, the binary mass ratio, the initial inclination angle between the orbit and disc planes and the Mach number in the outer parts of the unperturbed disc. For binary mass ratios of around unity and D/R in the range 3 to 4, we find that the global evolution of the discs is governed primarily by the value of the Mach number. For relatively low Mach numbers (i.e. 10 to 20) we find that the discs develop a mildly warped structure, are tidally truncated, and undergo a near rigid body precession at a rate which is in close agreement with analytical arguments. For higher Mach numbers (approx. 30), the evolution is towards a considerably more warped structure, but the disc none the less maintains itself as a long-lived, coherent entity. Additionally, it is found that the inclination angle between the disc and orbital angular momentum vectors evolves on a longer timescale which is probably the viscous evolution timescale of the disc. The calculations are relevant to a number of observed astrophysical phenomena, including the precession of jets associated with young stars, the high spectral index of some T-Tauri stars, and the light curves of X-ray binaries such as Hercules X-1 which suggest the presence of precessing accretion discs.