Normalized to: Terquem, C.
[1]
oai:arXiv.org:1909.08476 [pdf] - 1963800
Non-adiabatic tidal oscillations induced by a planetary companion
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
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
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
Gillen, Edward;
Aigrain, Suzanne;
Terquem, Caroline;
Bouvier, Jerome;
Alencar, Silvia H. P.;
Gandolfi, Davide;
Stauffer, John;
Cody, Ann Marie;
Venuti, Laura;
Almeida, Pedro Viana;
Micela, Giuseppina;
Favata, Fabio;
Deeg, Hans J.
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
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
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
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
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
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
Gillen, Edward;
Aigrain, Suzanne;
McQuillan, Amy;
Bouvier, Jerome;
Hodgkin, Simon;
Alencar, Silvia H. P.;
Terquem, Caroline;
Southworth, John;
Gibson, Neale P.;
Cody, Ann Marie;
Lendl, Monika;
Morales-Calderón, Maria;
Favata, Fabio;
Stauffer, John;
Micela, Giuseppina
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
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
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
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
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
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
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
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
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
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
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
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?
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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.