Source author record

C. J. Clarke

C. J. Clarke appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

24works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

24 published item(s)

preprint2025arXiv

Hot Jupiter formation in dense stellar clusters: A Monte Carlo model applied to 47 Tucanae

We study the efficiency of high-e migration as a pathway for Hot Jupiter formation in the dense globular cluster 47 Tuc. Gravitational N-body simulations are performed to investigate the orbital evolution of star-planet systems due to dynamical stellar perturbations. Planetary systems that have been scattered into orbits of sufficiently high eccentricity can undergo tidal circularisation, with Hot Jupiter formation being one possible stopping condition. We also account for the possibility of (i) ionisation due to high-energy encounters, (ii) tidal disruption of the planet by tidal forces inside the Roche limit and (iii) Warm Jupiter formation. The orbital evolution of a population of cold Jupiter progenitors, with initial semi-major axes between 1-30 au, is simulated over 12 Gyr using a simplified dynamical model of 47 Tuc. Our computational treatment of dynamical encounters yields an overall HJ occurrence rate of F_HJ = 5.9 x 10^-4 per cluster star (a 51 per cent enhancement relative to the analytic baseline). The probability of Hot Jupiter formation is highest in the core and falls off steeply beyond a few parsecs from the centre of the cluster, where the stellar density is too low to drive efficient eccentricity diffusion. The code can be found here: https://github.com/James-Wirth/HotJupiter.

preprint2022arXiv

Inward and outward migration of massive planets: moving towards a stalling radius

Recent studies on the planet-dominated regime of Type II migration showed that, contrary to the conventional wisdom, massive planets can migrate outwards. Using `fixed-planet' simulations these studies found a correlation between the sign of the torques acting on the planet and the parameter $K'$ (which describes the depth of the gap carved by the planet in the disc). We perform `live-planet' simulations exploring a range of $K'$ and disc mass values to test and extend these results. The excitation of planet eccentricity in live-planet simulations breaks the direct dependence of migration rate (rate of change of semi-major axis) on the torques imposed, an effect that `fixed-planet' simulations cannot treat. By disentangling the contribution to the torque due to the semi-major axis evolution from that due to the eccentricity evolution, we recover the relation between the magnitude and sign of migration and $K'$ and argue that this relation may be better expressed in terms of the related gap depth parameter $K$. We present a toy model in which the sign of planetary migration changes at a limiting value of $K$, through which we explore planets' migration in viscously evolving discs. The existence of the torque reversal shapes the planetary system's architecture by accumulating planets either at the stalling radius or in a band around it (defined by the interplay between the planet migration and the disc evolution). In either case, planets pile up in the area $1-10$ au, disfavouring hot Jupiter formation through Type II migration in the planet-dominated regime.

preprint2020arXiv

Forbidden line diagnostics of photoevaporative disc winds

Photoevaporation driven by high energy radiation from the central star plays an important role in the evolution of protoplanetary discs. Photoevaporative winds have been unambiguously detected through blue-shifted emission lines, but their detailed properties remain uncertain. Here we present a new empirical approach to make observational predictions of these thermal winds, seeking to fill the gap between theory and observations. We use a self-similar model of an isothermal wind to compute line profiles of several characteristic emission lines (the [Ne${\rm{\scriptsize II}}$] line at 12.81 $μ$m, and optical forbidden lines such as [O${\rm{\scriptsize I}}$] 6300 $\mathring{A}$ and [S${\rm{\scriptsize II}}$] 4068/4076 $\mathring{A}$), studying how the lines are affected by parameters such as the gas temperature, disc inclinations, and density profile. Our model successfully reproduces blue-shifted lines with $v_{\rm peak} \lesssim 10$ km/s, which decrease with increasing disc inclination. The line widths increase with increasing disc inclinations and range from $Δv \sim 15-30$ km/s. The predicted blue-shifts are mostly sensitive to the gas sound speed. The observed [Ne${\rm{\scriptsize II}}$] line profiles are consistent with a thermal wind and point towards a relatively high sound speed, as expected for EUV photoevaporation. However, the observed [O${\rm{\scriptsize I}}$] line profiles require lower temperatures, as expected in X-ray photoevaporation, and show a wider scatter that is difficult to reconcile with a single wind model; it seems likely that these lines trace different components of a multi-phase wind. We also note that the spectral resolution of current observations remains an important limiting factor in these studies, and that higher resolution spectra are required if emission lines are to further our understanding of protoplanetary disc winds.

preprint2016arXiv

Violent environment of the inner disk of RW Aur A probed by the 2010 and 2015 dimming events

RW Aur is a young binary system showing strong signatures of a recent tidal encounter between the circumprimary disk and the secondary star. The primary star has recently undergone two major dimming events ($Δ$mag $\approx$ 2 in V-band), whose origin is still under debate. To shed light on the mechanism leading to the dimming events, we study the extinction properties, accretion variability, and gas kinematics using absorption lines from the material obscuring star RW Aur A. We compare our moderate resolution X-Shooter spectra of the dim state of RW Aur A with other spectral observations. In particular, we analyse archival high resolution UVES spectra obtained during the bright state of the system, in order to track the evolution of the spectral properties across the second dimming event. The spectrum obtained during the dim state shows narrow absorption lines in the Na and K optical doublets, where the former is saturated. With a velocity of -60 km/s these lines indicate that during the dim state the disk wind is either enhanced, or significantly displaced into the line of sight. The photometric evolution across the dimming event shows a gray extinction, and is correlated with a significant reduction of the EW of all photospheric lines. Emission lines tracing accretion do not vary significantly across the dimming. We conclude that the dimming event is related to a major perturbation on the inner disk. We suggest that the inner disk is occulting (most of) the star, and thus its photosphere, but is not occulting the accretion regions within a few stellar radii. Since observations of the outer disk indicate that the disk is modestly inclined (45 - 60 deg), we propose that the inner disk might be warped by a yet unseen (sub-)stellar companion, which may also explain the 2.77 day periodic variability of the spectral lines.

preprint2014arXiv

Probing X-ray photoevaporative winds through their interaction with ionising radiation in cluster environments: the case for X-ray proplyds

We show that if young low mass stars undergo vigorous X-ray driven disc winds, these may be detected in clusters through their interaction with ionising radiation from massive stars. We argue that in the ONC (Orion Nebula Cluster) one should see $\sim$ 10s of `X-ray proplyds' ( objects with optically imaged offset ionisation fronts) in the range $0.3-0.6$pc from $θ_1$C Ori (the dominant O star in the ONC). Such objects lie outside the central `FUV zone' in the ONC where proplyds are instead well explained by neutral winds driven by external Far Ultraviolet (FUV) emission from $θ_1$C. We show that the predicted numbers and sizes of X-ray proplyds are compatible with those observed and that this may also explain at least some of the far flung proplyds seen in the Carina nebula. We compare the sizes of observed proplyds outside the FUV region of the ONC with model predictions based on the current observed X-ray luminosities of these sources ( bearing in mind that the current size is actually set by the X-ray luminosity a few hundred years ago, $\sim$ the flow time to the ionisation front). We discuss whether variability on this timescale can plausibly explain proplyd size data on a case by case basis. We also show that the predicted radio free-free emission signature of X-ray proplyds is readily detectable. Monitoring is however required to distinguish such emission from non-thermal emission from active coronae. We also predict that it is only at $> 1$ pc from $θ_1$C that the free-free emission from such offset ionised structures would be clearly distinguishable from an externally driven ionised disc wind. We argue that the fortuitous proximity of massive stars in the ONC can be used as a beacon to light up internally driven X-ray winds and that this is a promising avenue for observationally testing X-ray photoevaporation.

preprint2013arXiv

Evolutionary constraints on the planetary hypothesis for transition discs

We assume a scenario in which transition discs (i.e. discs around young stars that have signatures of cool dust but lack significant near infra-red emission from warm dust) are associated with the presence of planets (or brown dwarfs). These are assumed to filter the dust content of any gas flow within the planetary orbit and produce an inner `opacity hole'. In order to match the properties of transition discs with the largest (~50 A.U. scale) holes, we place such `planets' at large radii in massive discs and then follow the evolution of the tidally coupled disc-planet system, comparing the system's evolution in the plane of mm flux against hole radius with the properties of observed transition discs. We find that, on account of the high disc masses in these systems, all but the most massive `planets' (100 Jupiter masses) are conveyed to small radii by Type II migration without significant fading at millimetre wavelengths. Such behaviour would contradict the observed lack of mm bright transition discs with small (<10 A.U.) holes. On the other hand, imaging surveys clearly rule out the presence of such massive companions in transition discs. We conclude that this is a serious problem for models that seek to explain transition discs in terms of planetary companions unless some mechanism can be found to halt inward migration and/or suppress mm flux production. We suggest that the dynamical effects of substantial accretion on to the planet/through the gap may offer the best prospect for halting such migration but that further long term simulations are required to clarify this issue.

preprint2013arXiv

The relation between accretion rates and the initial mass function in hydrodynamical simulations of star formation

We analyse a hydrodynamical simulation of star formation. Sink particles in the simulations which represent stars show episodic growth, which is presumably accretion from a core that can be regularly replenished in response to the fluctuating conditions in the local environment. The accretion rates follow $\dot{m} \propto m^{2/3}$, as expected from accretion in a gas-dominated potential, but with substantial variations over-laid on this. The growth times follow an exponential distribution which is tapered at long times due to the finite length of the simulation. The initial collapse masses have an approximately lognormal distribution with already an onset of a power-law at large masses. The sink particle mass function can be reproduced with a non-linear stochastic process, with fluctuating accretion rates $\propto m^{2/3}$, a distribution of seed masses and a distribution of growth times. All three factors contribute equally to the form of the final sink mass function. We find that the upper power law tail of the IMF is unrelated to Bondi-Hoyle accretion.

preprint2011arXiv

Global mass segregation in hydrodynamical simulations of star formation

Recent analyses of mass segregation diagnostics in star forming regions invite a comparison with the output of hydrodynamic simulations of star formation. In this work we investigate the state of mass segregation of 'stars' (i.e. sink particles in the simulations) in the case of hydrodynamical simulations which omit feedback. We first discuss methods to quantify mass segregation in substructured regions, either based on the minimum spanning tree (Allison's Lambda), or through analysis of correlations between stellar mass and local stellar surface number densities. We find that the presence of even a single 'outlier' (i.e. a massive object far from other stars) can cause the Allison Lambda method to describe the system as inversely mass segregated, even where in reality the most massive sink particles are overwhelmingly in the centres of the subclusters. We demonstrate that a variant of the Lambda method is less susceptible to this tendency but also argue for an alternative representation of the data in the plane of stellar mass versus local surface number density. The hydrodynamical simulations show global mass segregation from very early times which continues throughout the simulation, being only mildly influenced during sub-cluster merging. We find that up to approx. 2-3% of the "massive" sink particles (m > 2.5 Msun) are in relative isolation because they have formed there, although other sink particles can form later in their vicinity. Ejections of massive sinks from subclusters do not contribute to the number of isolated massive sink particles, as the gravitational softening in the calculation suppresses this process.

preprint2011arXiv

Stability of self-gravitating discs under irradiation

Self-gravity becomes competitive as an angular momentum transport process in accretion discs at large radii, where the temperature is low enough that external irradiation likely contributes to the thermal balance. Irradiation is known to weaken the strength of disc self-gravity, and can suppress it entirely if the disc is maintained above the threshold for linear instability. However, its impact on the susceptibility of the disc to fragmentation is less clear. We use two-dimensional numerical simulations to investigate the evolution of self-gravitating discs as a function of the local cooling time and strength of irradiation. In the regime where the disc does not fragment, we show that local thermal equilibrium continues to determine the stress - which can be represented as an effective viscous alpha - out to very long cooling times (at least 240 dynamical times). In this regime, the power spectrum of the perturbations is uniquely set by the effective viscous alpha and not by the cooling rate. Fragmentation occurs for cooling times tau < beta_crit / Omega, where beta_crit is a weak function of the level of irradiation. We find that beta_crit declines by approximately a factor of two, as irradiation is increased from zero up to the level where instability is almost quenched. The numerical results imply that irradiation cannot generally avert fragmentation of self-gravitating discs at large radii; if other angular momentum transport sources are weak mass will build up until self-gravity sets in, and fragmentation will ensue.

preprint2010arXiv

OB Stars in Stochastic Regimes

The highest-mass stars have the lowest frequency in the stellar IMF, and they are also the most easily observed stars. Thus, the counting statistics for OB stars provide important tests for the fundamental nature and quantitative parameters of the IMF. We first examine some local statistics for the stellar upper-mass limit itself. Then, we examine the parameter space and statistics for extremely sparse clusters that contain OB stars, in the SMC. We find that thus far, these locally observed counting statistics are consistent with a constant stellar upper-mass limit. The sparse OB star clusters easily fall within the parameter space of Monte Carlo simulations of cluster populations. If the observed objects are representative of their cluster birth masses, their existence implies that the maximum stellar mass is largely independent of the parent cluster mass.

preprint2010arXiv

Properties of hierarchically forming star clusters

We undertake a systematic analysis of the early (< 0.5 Myr) evolution of clustering and the stellar initial mass function in turbulent fragmentation simulations. These large scale simulations for the first time offer the opportunity for a statistical analysis of IMF variations and correlations between stellar properties and cluster richness. The typical evolutionary scenario involves star formation in small-n clusters which then progressively merge; the first stars to form are seeds of massive stars and achieve a headstart in mass acquisition. These massive seeds end up in the cores of clusters and a large fraction of new stars of lower mass is formed in the outer parts of the clusters. The resulting clusters are therefore mass segregated at an age of 0.5 Myr, although the signature of mass segregation is weakened during mergers. We find that the resulting IMF has a smaller exponent (alpha=1.8-2.2) than the Salpeter value (alpha=2.35). The IMFs in subclusters are truncated at masses only somewhat larger than the most massive stars (which depends on the richness of the cluster) and an universal upper mass limit of 150 Msun is ruled out. We also find that the simulations show signs of the IGIMF effect proposed by Weidner & Kroupa, where the frequency of massive stars is suppressed in the integrated IMF compared to the IMF in individual clusters. We identify clusters through the use of a minimum spanning tree algorithm which allows easy comparison between observational survey data and the predictions of turbulent fragmentation models. In particular we present quantitative predictions regarding properties such as cluster morphology, degree of mass segregation, upper slope of the IMF and the relation between cluster richness and maximum stellar mass. [abridged]

preprint2009arXiv

Radiation-Hydrodynamic Models of X-Ray & EUV Photoevaporating Protoplanetary Discs

(Abridged) We present the first radiation-hydrodynamic model of a protoplanetary disc irradiated with an X-EUV spectrum. In a model where the total ionizing luminosity is divided equally between X-ray and EUV luminosity, we find a photoevaporation rate of 1.4e-8 M_sun/yr, which is two orders of magnitude greater than the case of EUV photoevaporation alone. Thus it is clear that the X-rays are the dominant driving mechanism for photoevaporation. This can be understood inasmuch as X-rays are capable of penetrating much larger columns (~1e22 cm^-2) and can thus effect heating in denser regions and at larger radius than the EUV can. The radial extent of the launching region of the X-ray heated wind is 1-70AU compared with the pure EUV case where the launch region is concentrated around a few AU. When we couple our wind mass-loss rates with models for the disc's viscous evolution, we find that, as in the pure EUV case, there is a photoevaporative switch, such that an inner hole develops at ~ 1 AU at the point that the accretion rate in the disc drops below the wind mass loss rate. At this point, the remaining disc material is quickly removed in the final 15-20% of the disc's lifetime. This is consistent with the 1e5 yr transitional timescale estimated from observations of T-Tauri stars. We also caution that although our mass-loss rates are high compared to some accretion rates observed in young stars, our model has a rather large X-ray luminosity of 2e30 erg/s; further modeling is required in order to investigate the evolutionary implications of the large observed spread of X-ray luminosities in T-Tauri stars.

preprint2008arXiv

Gas Phase diagnostics of Protoplanetary disk extension

We investigate the potential of using ratios of fine structure and near-infrared forbidden line transitions of atomic carbon to diagnose protoplanetary disk extension. Using results from 2D photoionisation and radiative transfer modeling of a realistic protoplanetary disk structure irradiated by X-rays from a T Tauri star, we obtain theoretical emission maps from which we construct radial distributions of the strongest emission lines produced in the disk. We show that ratios of fine structure to near-infrared forbidden line emission of atomic carbon are especially promising to constrain the minimum size of gaseous protoplanetary disks. While theoretically viable, the method presents a number of observational difficulties that are also discussed here.

preprint2006arXiv

Disk Evolution in Young Binaries: from Observations to Theory

The formation of a binary system surrounded by disks is the most common outcome of stellar formation. Hence studying and understanding the formation and the evolution of binary systems and associated disks is a cornerstone of star formation science. Moreover, since the components within binary systems are coeval and the sizes of their disks are fixed by the tidal truncation of their companion, binary systems provide an ideal "laboratory" in which to study disk evolution under well defined boundary conditions. In this paper, we review observations of several inner disk diagnostics in multiple systems, including hydrogen emission lines (indicative of ongoing accretion), $K-L$ and $K-N$ color excesses (evidence of warm inner disks), and polarization (indicative of the relative orientations of the disks around each component). We examine to what degree these properties are correlated within binary systems and how this degree of correlation depends on parameters such as separation and binary mass ratio. These findings will be interpreted both in terms of models that treat each disk as an isolated reservoir and those in which the disks are subject to re-supply from some form of circumbinary reservoir, the observational evidence for which we will also critically review. The planet forming potential of multiple star systems is discussed in terms of the relative lifetimes of disks around single stars, binary primaries and binary secondaries. Finally, we summarize several potentially revealing observational problems and future projects that could provide further insight into disk evolution in the coming decade

preprint2005arXiv

Photoionising feedback in star cluster formation

We present the first ever hydrodynamic calculations of star cluster formation that incorporate the effect of feedback from ionising radiation. In our simulations, the ionising source forms in the cluster core at the intersection of several dense filaments of inflowing gas. We show that these filaments collimate ionised outflows and suggest such an environmental origin for at least some observed outflows in regions of massive star formation. Our simulations show both positive feedback (i.e. promotion of star formation in neutral gas compressed by expanding HII regions) and negative feedback (i.e. suppression of the accretion flow in to the central regions). We show that the volume filling factor of ionised gas is very different in our simulations than would result from the case where the central source interacted with an azimuthally smoothed gas density distribution. As expected, gas density is the key parameter in determining whether clusters are unbound by photoionising radiation. Nevertheless, we find - on account of the acceleration of a small fraction of the gas to high velocities in the outflows - that the deposition in the gas of an energy that exceeds the binding energy of the cluster is not a sufficient criterion for unbinding the bulk of the cluster mass.

preprint2005arXiv

Statistical Confirmation of a Stellar Upper Mass Limit

We derive the expectation value for the maximum stellar mass (m_max) in an ensemble of N stars, as a function of the IMF upper-mass cutoff (m_up) and N. We statistically demonstrate that the upper IMF of the local massive star census observed thus far in the Milky Way and Magellanic Clouds clearly exhibits a universal upper mass cutoff around 120 - 200 M_sun for a Salpeter IMF, although the result is more ambiguous for a steeper IMF.

preprint2001arXiv

Competitive accretion in embedded stellar cluster

We investigate the physics of gas accretion in young stellar clusters. Accretion in clusters is a dynamic phenomenon as both the stars and the gas respond to the same gravitational potential. Accretion rates are highly non-uniform with stars nearer the centre of the cluster, where gas densities are higher, accreting more than others. This competitive accretion naturally results in both initial mass segregation and a spectrum of stellar masses. Accretion in gas-dominated clusters is well modelled using a tidal-lobe radius instead of the commonly used Bondi-Hoyle accretion radius. This works as both the stellar and gas velocities are under the influence of the same gravitational potential and are thus comparable. The low relative velocity that results means that the tidal radius is smaller than the Bondi-Hoyle radius in these systems. In contrast, when the stars dominate the potential and are virialised, the Bondi-Hoyle radius is smaller than the tidal radius and thus Bondi-Hoyle accretion is a better fit to the accretion rates.

preprint2000arXiv

Jet directions in Seyfert galaxies

We present the study of the relative angle between the accretion disk (or radio jet) and the galaxy disk for a sample of Seyfert galaxies selected from a mostly isotropic property, the 60$μ$m flux. For each galaxy we have a pair of points ($i$,$δ$), which are the inclination of the galaxy relative to the line of sight and the angle between the jet projected into the plane of the sky and the host galaxy major axis, respectively. This data is combined with a statistical technique to determine the distribution of $β$ angles {\it in 3 dimensions}, the angle between the jet and the host galaxy plane axis. The analysis of the data, not differentiating between Seyfert 1's and 2's, showed that the observed distribution of $i$ and $δ$ values can be well represented by a homogeneous $β-$distribution in the range $0^{\circ}\leqβ\leq90^{\circ}$. However, when we distinguished between Seyfert 1's and 2's, the models could not represent Seyfert 1's adequately. It was necessary to introduce viewing angle restrictions, that a galaxy can only be recognized as a Seyfert 1 if the angle between the jet and the line of sight ($|ϕ|$) is smaller than a given angle $ϕ_c$ and that the galaxy inclination $i$ is smaller than an angle $i_c$, in order to have statistically acceptable models. This is a direct and independent confirmation of the underlying concepts of the Unified Model. We discuss several ways to explain the misalignment between the accretion disk axis and the host galaxy disk axis.

preprint1997arXiv

The Superbubble Size Distribution in the Interstellar Medium of Galaxies

We use the standard, adiabatic shell evolution to predict the size distribution N(R) for populations of OB superbubbles in a uniform ISM. We derive N(R) for simple cases of superbubble creation rate and mechanical luminosity function (MLF). For R < the characteristic radius R_e, N(R) is dominated by stalled objects, while for R>R_e it is dominated by growing objects. We also briefly investigate N(R) resulting from momentum-conserving shell evolution. We predict a peak in N(R) corresponding to individual SNRs. To estimate the MLF, we also examine evolutionary effects on the HII region luminosity function (HII LF), finding that for nebular luminosity fading as a power law in time, there is a minimum observed slope for the HII LFs. Comparison with the largely complete HI hole catalog for the SMC shows surprising agreement in the predicted and observed slope of N(R), suggesting that no other fundamental process is needed to explain the size distribution of shells in the SMC. Further comparison with largely incomplete HI data for M31, M33, and Holmberg II is also encouraging. We present expressions for the ISM porosity parameters, and estimate that they are substantially <1 for all of the galaxies except Holmberg II. Most of these galaxies therefore may not be strongly dominated by a hot interstellar component. However, porosity results for the Galaxy remain inconclusive.

preprint1996arXiv

The ejection of T Tauri stars from molecular clouds and the fate of circumstellar discs

We investigate the evolution of circumstellar discs around T Tauri stars that are ejected from small stellar clusters within molecular clouds. In particular, we study how the interaction that leads to ejection may hasten the transition between Classical and Weak-lined T Tauri status. In our models, ejections of T Tauri stars at velocities of 3-10 km/s truncate the accretion disc at radii between 1 and 10 a.u., reducing the viscous evolution time of the disc so that accretion rapidly ceases. The observational appearance of the resulting systems is then dependent on the presence or absence of a stellar magnetic field. For non-magnetic stars we find that a near-infra red excess should persist due to reprocessing of stellar radiation, but that this is greatly diminished for magnetic T Tauri stars by the presence of a magnetosphere extending to corotation. In either case, there is a period when ejected stars should appear as non-accreting systems with detectable circumstellar material at wavelengths of 5 microns and beyond. We discuss the implications of these results for models in which ejected stars contribute to the halo of pre-main-sequence objects discovered from ROSAT observations of star forming regions and the All-Sky Survey.

preprint1995arXiv

Energetics of star-disc encounters in the non-linear regime

We investigate the response of a circumstellar accretion disc to the fly-by of a perturbing mass on a parabolic orbit. The energy and angular momentum transferred during the encounter are calculated using a reduced three-body method. In almost all close encounters the energy and angular momentum transfer is dominated by disc material becoming unbound from the system, with the contributions from close disc particle -- star encounters being significant. For more distant encounters with some prograde element to the motion the disc material loses energy and angular momentum to the perturber's orbit through a resonance feature. The magnitude of the energy transfer calculated in our simulations is greater than that of the binding energy of material exterior to periastron by a factor of two in the prograde case, and up to a factor of five in the case of the retrograde encounter. The destructive nature of the encounters indicates that a non-linear treatment is essential in all but the most distant encounters.

preprint1995arXiv

Magnetic braking of T Tauri stars

We construct models for the rotation rates of T Tauri stars whose spin is regulated by magnetic linkage between the star and a surrounding accretion disc. Our models utilise a time-dependent disc code to follow the accretion process and include the effects of pre-main-sequence stellar evolution. We find that the initial disc mass controls the evolution of the star-disc system. For sufficiently massive discs, a stellar field of $\sim$ 1 kG is able to regulate the spin rate to the observed values during the classical T Tauri phase. The field then acts to expel the disc and the star spins up at constant angular momentum as a weak-line system. Lower mass discs are ejected at an early epoch and fail to brake the star significantly. We extend the model to close binary systems, and find that the removal of angular momentum from the disc by the secondary significantly prolongs the inner disc lifetime. Such systems should therefore be relatively slow rotators. We also discuss the implications of our model for the spectral energy distributions and variability of T Tauri stars.

preprint1995arXiv

Magnetically modulated accretion in T Tauri stars

We examine how accretion on to T Tauri stars may be modulated by a time-dependent `magnetic gate' where the inner edge of the accretion disc is disrupted by a varying stellar field. We show that magnetic field variations on time-scales shorter than 10^5 yr can modulate the accretion flow, thus providing a possible mechanism both for the marked photometric variability of T Tauri stars and for the possible conversion of T Tauri stars between classical and weak line status. We thus suggest that archival data relating to the spectrophotometric variability of T Tauri stars may provide an indirect record of magnetic activity cycles in low-mass pre-main-sequence stars.

preprint1993arXiv

N-Body Simulations of Star-Disc Captures in Globular Clusters

The presence of protostellar disks can greatly increase the dissipation during close stellar encounters, leading to the formation of a significant population of binaries during the initial collapse and virialization of a cluster. We have used N-body simulations of collapsing globular clusters to find the major factors that determine the efficiency of binary formation through star-disk captures. This work serves the dual purpose of verifying the results of earlier analytic work as well as examining parameters not testable by that work. As in the earlier work, typical binary fractions of a few percent are found. For the parameters studied, the results are found to depend remarkably little upon disk evolution, the mass distribution of the stars, or their spatial distribution, though distributions in which the stars are highly clumped yield binary fractions larger by a factor of a few. The direct N-body integrations limit the models to relatively small values of N. Semiempirical relations are derived, however, which allow the results to be extrapolated to values of N appropriate to globular clusters.