We present monitoring observations of the active T Tauri star RW Aur, from 2010 October to 2015 January, using optical high-resolution (R>10000) spectroscopy with CFHT-ESPaDOnS. Optical photometry in the literature shows bright, stable fluxes over most of this period, with lower fluxes (by 2-3 mag.) in 2010 and 2014. In the bright period our spectra show clear photospheric absorption, complicated variation in the Ca II 8542 A emission}profile shapes, and a large variation in redshifted absorption in the O I 7772 and 8446 A and He I 5876 A lines, suggesting unstable mass accretion during this period. In contrast, these line profiles are relatively uniform during the faint periods, suggesting stable mass accretion. During the faint periods the photospheric absorption lines are absent or marginal, and the averaged Li I profile shows redshifted absorption due to an inflow. We discuss (1) occultation by circumstellar material or a companion and (2) changes in the activity of mass accretion to explain the above results, together with near-infrared and X-ray observations from 2011-2015. Neither scenario can simply explain all the observed trends, and more theoretical work is needed to further investigate their feasibilities.

We present a study of the wind launching region of the Herbig Be star HD 58647 using high angular (lambda/2B=0.003") and high spectral (R=12000) resolution interferometric VLTI-AMBER observations of the near-infrared hydrogen emission line, Br-gamma. The star displays double peaks in both Br-gamma line profile and wavelength-dependent visibilities. The wavelength-dependent differential phases show S-shaped variations around the line centre. The visibility level increases in the line (by ~0.1) at the longest projected baseline (88 m), indicating that the size of the line emission region is smaller than the size of the K-band continuum-emitting region, which is expected to arise near the dust sublimation radius of the accretion disc. The data have been analysed using radiative transfer models to probe the geometry, size and physical properties of the wind that is emitting Br-gamma. We find that a model with a small magnetosphere and a disc wind with its inner radius located just outside of the magnetosphere can well reproduce the observed Br-gamma profile, wavelength-dependent visibilities, differential and closure phases, simultaneously. The mass-accretion and mass-loss rates adopted for the model are Mdot_a = 3.5 x 10^{-7} Msun/yr and Mdot_dw = 4.5 x 10^{-8} Msun/yr, respectively (Mdot_dw/Mdot_a =0.13). Consequently, about 60 per cent of the angular momentum loss rate required for a steady accretion with the measured accretion rate is provide by the disc wind. The small magnetosphere in HD 58647 does not contribute to the Br-gamma line emission significantly.

The origin of the near-infrared (NIR) HI emission lines in young stellar objects are not yet understood. To probe it, we present multi-epoch LBT-LUCIFER spectroscopic observations of the Pa{\delta}, Pa{\beta}, and Br{\gamma} lines observed in the Herbig star VVSer, along with VLTI-AMBER Br{\gamma} spectro-interferometric observations at medium resolution. Our spectroscopic observations show line profile variability in all the HI lines. The strongest variability is observed in the redshifted part of the line profiles. The Br{\gamma} spectro-interferometric observations indicate that the Br{\gamma} line emitting region is smaller than the continuum emitting region. To interpret our results, we employed radiative transfer models with three different flow configurations: magnetospheric accretion, a magneto-centrifugally driven disc wind, and a schematic bipolar outflow. Our models suggest that the HI line emission in VVSer is dominated by the contribution of an extended wind, perhaps a bipolar outflow. Although the exact physical process for producing such outflow is not known, this model is capable of reproducing the averaged single-peaked line profiles of the HI lines. Additionally, the observed visibilities, differential and closure phases are best reproduced when a wind is considered. Nevertheless, the complex line profiles and variability could be explained by changes in the relative contribution of the magnetosphere and/or winds to the line emission. This might indicate that the NIR HI lines are formed in a complex inner disc region where inflow and outflow components might coexist. Furthermore, the contribution of each of these mechanisms to the line appears time variable, suggesting a non-steady accretion/ejection flow.

The mass outflows in T Tauri stars (TTS) are thought to be an effective mechanism to remove angular momentum during the pre-main-sequence contraction of a low-mass star. The most powerful winds are observed at the FUor stage of stellar evolution. V1331 Cyg has been considered as a TTS at the pre-FUor stage. We analyse high-resolution spectra of V1331 Cyg collected in 1998-2007 and 20-d series of spectra taken in 2012. For the first time the photospheric spectrum of the star is detected and stellar parameters are derived: spectral type G7-K0 IV, mass 2.8 Msun, radius 5 Rsun, vsini < 6 km/s. The photospheric spectrum is highly veiled, but the amount of veiling is not the same in different spectral lines, being lower in weak transitions and much higher in strong transitions. The Fe II 5018, Mg I 5183, K I 7699 and some other lines of metals are accompanied by a `shell' absorption at radial velocity of about -240 km/s. We show that these absorptions form in the post-shock gas in the jet, i.e. the star is seen though its jet. The P Cyg profiles of H-alpha and H-beta indicate the terminal wind velocity of about 500 km/s, which vary on time-scales from several days to years. A model of the stellar wind is developed to interpret the observations. The model is based on calculation of hydrogen spectral lines using the radiative transfer code TORUS. The observed H-alpha and H-beta line profiles and their variability can be well reproduced with a stellar wind model, where the mass-loss rate and collimation (opening angle) of the wind are variable. The changes of the opening angle may be induced by small variability in magetization of the inner disc wind. The mass-loss rate is found to vary within (6-11)x10^{-8} Msun/yr, with the accretion rate of 2.0x10^{-6} Msun/yr.

We discuss the results of modelling of young magnetized stars, where matter flow is calculated using the three-dimensional (3D) magneto-hydrodynamic (MHD) "Cubed Sphere" code, and the spectra are calculated using the 3D Monte Carlo radiative transfer code "TORUS". Two examples of modelling are shown: (1) accretion onto stars in stable and unstable regimes, and (2) accretion to a young star V2129 Oph, modelled with realistic parameters.

We review recent axisymmetric and three-dimensional (3D) magnetohydrodynamic (MHD) numerical simulations of magnetospheric accretion, plasma-field interaction and outflows from the disk-magnetosphere boundary.

Recent spectropolarimetric observations suggest that young low-mass stars such as classical T Tauri stars (CTTSs) possess relatively strong (~kG) magnetic field. This supports a scenario in which the final accretion onto the stellar surface proceeds through a magnetosphere, and the winds are formed in magnetohydrodynamics (MHD) processes. We examine recent numerical simulations of magnetospheric accretions via an inclined dipole and a complex magnetic fields. The difference between a stable accretion regime, in which accretion occurs in ordered funnel streams, and an unstable regime, in which gas penetrates through the magnetosphere in several unstable streams due to the magnetic Rayleigh-Taylor instability, will be discussed. We describe how MHD simulation results can be used in separate radiative transfer (RT) models to predict observable quantiles such as line profiles and light curves. The plausibility of the accretion flows and outflows predicted by MHD simulations (via RT models) can be tested against observations. We also address the issue of outflows/winds that arise from the innermost part of CTTSs. First, we discuss the line formations in a simple disk wind and a stellar wind models. We then discuss the formation of the conically shaped magnetically driven outflow that arises from the disk-magnetosphere boundary when the magnetosphere is compressed into an X-type configuration.

We discuss key observational signatures of Classical T Tauri stars (CTTSs) accreting through Rayleigh-Taylor instability, which occurs at the interface between an accretion disk and a stellar magnetosphere. In this study, the results of global 3-D MHD simulations of accretion flows, in both stable and unstable regimes, are used to predict the variability of hydrogen emission lines and light curves associated with those two distinctive accretion flow patterns. In the stable regime, a redshifted absorption component (RAC) periodically appears in some hydrogen lines, but only during a fraction of a stellar rotation period. In the unstable regime, the RAC is present rather persistently during a whole stellar rotation period, and its strength varies non-periodically. The latter is caused by multiple accreting streams, formed randomly due to the instability, passing across the line of sight to an observer during one stellar rotation. This results in the quasi-stationarity appearance of the RAC because at least one of the accretion stream is almost always in the line of sight to an observer. In the stable regime, two stable hot spots produce a smooth and periodic light curve that shows only one or two peaks per stellar rotation. In the unstable regime, multiple hot spots formed on the surface of the star, produce the stochastic light curve with several peaks per rotation period.

Classical T Tauri stars (CTTSs) are variable in different time-scales. One type of variability is possibly connected with the accretion of matter through the Rayleigh-Taylor instability that occurs at the interface between an accretion disc and a stellar magnetosphere. In this regime, matter accretes in several temporarily formed accretion streams or `tongues' which appear in random locations, and produce stochastic photometric and line variability. We use the results of global three-dimensional magnetohydrodynamic simulations of matter flows in both stable and unstable accretion regimes to calculate time-dependent hydrogen line profiles and study their variability behaviours. In the stable regime, some hydrogen lines (e.g. H-beta, H-gamma, H-delta, Pa-beta and Br-gamma) show a redshifted absorption component only during a fraction of a stellar rotation period, and its occurrence is periodic. However, in the unstable regime, the redshifted absorption component is present rather persistently during a whole stellar rotation cycle, and its strength varies non-periodically. In the stable regime, an ordered accretion funnel stream passes across the line of sight to an observer only once per stellar rotation period while in the unstable regime, several accreting streams/tongues, which are formed randomly, pass across the line of sight to an observer. The latter results in the quasi-stationarity appearance of the redshifted absorption despite the strongly unstable nature of the accretion. In the unstable regime, multiple hot spots form on the surface of the star, producing the stochastic light curve with several peaks per rotation period. This study suggests a CTTS that exhibits a stochastic light curve and a stochastic line variability, with a rather persistent redshifted absorption component, may be accreting in the unstable accretion regime.

We present the emission line profile models of hydrogen and helium based on the results from axisymmetric magnetohydrodynamics (MHD) simulations of the wind formed near the disk-magnetosphere boundary of classical T Tauri stars (CTTSs). We extend the previous outflow models of `the conical-shell wind' by Romanova et al. to include a well defined magnetospheric accretion funnel flow which is essential for modelling the optical and near-infrared hydrogen and helium lines of CTTSs. The MHD model with an intermediate mass-accretion rate shows outflows in conical-shell shape with a half opening angle about 35 degrees. The flow properties such as the maximum outflow speed in the conical-shell wind, maximum inflow speed in the accretion funnel, mass-accretion and mass-loss rates are comparable to those found in a typical CTTS. The density, velocity and modified temperature from the MHD simulations are used in a separate radiative transfer model to predict the line profiles and test the consistency of the MHD models with observations. The line profiles are computed with various combinations of X-ray luminosities, temperatures of X-ray emitting plasma, and inclination angles. A rich diversity of line profile morphology is found, and many of the model profiles are very similar to those found in observations. We find that the conical-shell wind may contribute to the emission in some hydrogen lines (e.g. H-alpha, H-beta, Pa-beta and Pa-gamma) significantly when the temperature in the wind is relatively high (e.g. \sim 10^{4} K); however, the wind contribution decreases rapidly when a lower wind temperature is adopted. The model well reproduces a relatively narrow and low-velocity blueshifted absorption component in He I (10830), which are often seen in observations.

Ultraviolet and X-ray observations show evidence of outflowing gas around many active galactic nuclei. Some of these outflows may be driven off gas infalling towards the central black hole. We perform radiative transfer calculations to compute the gas ionization state and X-ray spectra for two- and three-dimensional (3D) hydrodynamical simulations of this outflow-from-inflow scenario. By comparison with observations, our results can be used to test the theoretical models and guide future numerical simulations. We predict both absorption and emission features, most of which are formed in a polar funnel of outflowing gas. This outflow causes strong absorption for observer orientation angles of < 35 degrees. Particularly in 3D, the strength of this absorption varies significantly for different lines-of-sight owing to the fragmentary structure of the gas flow. Although infalling material occupies a large fraction of the simulation volume, we do not find that it imprints strong absorption features since the ionization state is very high. Thus, an absence of observed inflow absorption features does not exclude the models. The main spectroscopic consequence of the infalling gas is a scattered continuum component that partially re-fills the absorption features caused by the outflowing polar funnel. Fluorescence and scattering in the outflow is predicted to give rise to several emission features for all observer orientations. For the hydrodynamical simulations considered we find both ionization states and column densities for the outflowing gas that are too high to be quantitatively consistent with well-observed X-ray absorption systems. Nevertheless, our results are qualitatively encouraging and further exploration of the model parameter space is warranted. (Abridged.)

We analyze the photometric and spectroscopic variability of the classical T Tauri star V2129 Oph over several rotational cycles to test the dynamical predictions of magnetospheric accretion models. The photometric variability and the radial velocity variations in the photospheric lines can be explained by rotational modulation due to cold spots, while the radial velocity variations of the He I (5876 \AA) line and the veiling variability are due to hot spot rotational modulation. The hot and cold spots are located at high latitudes and about the same phase, but the hot spot is expected to sit at the chromospheric level, while the cold spot is at the photospheric level. Using the dipole+octupole magnetic-field configuration previously proposed in the literature for the system, we compute 3D MHD magnetospheric simulations of the star-disk system. We use the simulation's density, velocity and scaled temperature structures as input to a radiative transfer code, from which we calculate theoretical line profiles at all rotational phases. The theoretical profiles tend to be narrower than the observed ones, but the qualitative behavior and the observed rotational modulation of the H\alpha and H\beta emission lines are well reproduced by the theoretical profiles. The spectroscopic and photometric variability observed in V2129 Oph support the general predictions of complex magnetospheric accretion models with non-axisymmetric, multipolar fields.

We present multidimensional non-LTE radiative transfer models of hydrogen and helium line profiles formed in the accretion flows and the outflows near the star-disk interaction regions of classical T Tauri stars (CTTSs). The statistical equilibrium calculations, performed under the assumption of the Sobolev approximation using the radiative transfer code TORUS, has been improved to include He I and He II energy levels. This allows us to probe the physical conditions of the inner wind of CTTSs by simultaneously modelling the robust wind diagnostic line He I (10830) and the accretion diagnostic lines such as Pa-beta, Br-gamma and He I (5876). The code has been tested in 1 and 2-D problems, and we have shown that the results are in agreement with established codes. We apply the model to the complex flow geometries of CTTSs. Example model profiles are computed using the combinations of (1) magnetospheric accretion and disc wind, and (2) magnetospheric accretion and the stellar wind. In both cases, the model produces line profiles which are qualitatively similar to those found in observations. Our models are consistent with the scenario in which the narrow blueshifted absorption component of He I (10830) seen in observations is caused by a disc wind, and the wider blueshifted absorption component (the P-Cygni profile) is caused by a bipolar stellar wind. However, we do not have a strong constraint on the relative importance of the wind and the magnetosphere for the `emission' component. Our preliminary calculations suggest that the temperatures of the disc wind, stellar wind and the magnetosphere cannot be much higher than ~10,000 K, on the basis of the strengths of hydrogen lines. With these low temperatures, we find that the photoionzation by high energy photons (e.g. X-rays) is necessary to produce He I (10830) in emission and to produce the blueshifted absorption components.

We summarize the results from numerical simulations of mass outflows from AGN. We focus on simulations of outflows driven by radiation from large-scale inflows. We discuss the properties of these outflows in the context of the so-called AGN feedback problem. Our main conclusion is that this type of outflows are efficient in removing matter but inefficient in removing energy.

We present a brief summary of the main results from our multi-dimensional, time-dependent simulations of gas dynamics in AGN. We focus on two types of outflows powered by radiation emitted from the AGN: disk winds and winds driven from large-scale inflows. We show spectra predicted by the simulations and discuss their relevance to observations of broad- and narrow-line regions of the AGN. We finish with a few remarks on whether these outflows can have a significant impact on their environment and host galaxy.

We review the results from multi-dimensional, time-dependent simulations of gas dynamics in AGN. We will focus on two types of outflows powered by radiation emitted from the AGN central engine: (i) outflows driven from the innermost part of an accretion disk and (2) outflows driven from a large-scale inflow that is likely the main supplier of material to the central engine. We discuss the relevance of both types of outflows to the so-called AGN feedback problem. However, the AGN feedback should not be considered separately from the AGN physics. Therefore, we also discuss the issue whether the properties of the same outflows are consistent with the gas properties in broad- and narrow-line regions.

(Abridged) We measure and analyze the energy, momentum, and mass feedback efficiencies due to radiation from AGN in relatively large scale outflows. Our measurements are based on the axisymmetric and time-dependent radiation-hydrodynamical simulations recently presented in Kurosawa & Proga. These simulations follow dynamics of gas under the influences of the AGN radiation and the gravity of the central 10^8 Msun black hole on scales from ~0.01 to ~10pc. We compare our Mdot_a-rho_o relation with that predicted by the Bondi accretion model. For high luminosities comparable to the Eddington limit, the power-law fit (Mdot_a ~ rho_o^q) to our models yields q=~0.5 instead of q=1.0 which is predicted by the Bondi model. This difference is caused by the outflows which are important for the overall mass budget at high luminosities. The maximum momentum and mass feedback efficiencies found in our models are ~10^-2 and ~10^-1, respectively. However, the outflows are much less important energetically: their thermal and kinetic powers in units of the radiative luminosity are ~10^-5 and ~10^-4, respectively. The efficiencies do not increase monotonically with the accretion luminosity but rather peak around the Eddington limit beyond which a steady state disk-wind-like solution exists. Our energy feedback efficiencies are significantly lower than 0.05, which is required in some cosmological and galaxy merger simulations. The low feedback efficiencies found here could have significant implications on the mass growth of super massive black holes in the early universe. We stress however that we have not considered the innermost parts of the accretion and outflow where radiation and matter interact most strongly. The feedback from this region could have efficiencies significantly above the low values found here.

(Abridged) We present two-dimensional hydrodynamical simulations of slowly rotating gas that is under the influence of the gravity of a super massive black hole and is irradiated by a thin UV accretion disc and a spherical X-ray corona. We calculate the accretion luminosity of a system based on the accretion-rate which is assumed to be equal to the mass-supply rate at the radius of ~10^-2 pc. For the models with high temperature gas at large radii (~10 pc) and high luminosities, we find a strong correlation between the mass-outflow rate (Mdot_out) and the luminosity (L). The power law index (q) describing the Mdot_out-L relation is q=2.0(+/-0.1), which is very similar to that for radiation-driven stellar and disc wind models. More surprisingly, for high density at large radii, we find steady state solutions with the accretion luminosity exceeding the Eddington limit. The super-Eddington accretion proceeds in the equatorial region and is possible because the radiation flux from the disc is significantly reduced in the equatorial direction due to the geometrical foreshortening effect. In all models, an outflow is driven from an inflow with sub-Keplerian rotation. For high temperature at large radii, the inflow occurs over a wide range of the polar angles, and the outflow occurs in a relatively narrow polar cone. However, for the super-Eddington cases with low temperature at large radii, the inflow persists only very close to the equatorial plane, resembling a thin accretion disc, while the outflow arises in a wide range of radii and polar angles. The geometry of this extreme inflow-outflow solution is very similar to a radiation-driven wind from a luminous Keplerian accretion disc.

We study the axisymmetric and non-axisymmetric, time-dependent hydrodynamics of gas that is under the influence of the gravity of a super massive black hole (SMBH) and the radiation force produced by a radiatively efficient flow accreting onto the SMBH. We have considered two cases: (1) the formation of an outflow from the accretion of the ambient gas without rotation and (2) that with weak rotation. The main goals of this study are: (1) to examine if there is a significant difference between the models with identical initial and boundary conditions but in different dimensionality (2-D and 3-D), and (2) to understand the gas dynamics in AGN. Our 3-D simulations of a non-rotating gas show small yet noticeable non-axisymmetric small-scale features inside the outflow. The outflow as a whole and the inflow do not seem to suffer from any large-scale instability. In the rotating case, the non-axisymmetric features are very prominent, especially in the outflow which consists of many cold dense clouds entrained in a smoother hot flow. The 3-D outflow is non-axisymmetric due to the shear and thermal instabilities. In both 2-D and 3-D simulations, gas rotation increases the outflow thermal energy flux, but reduces the outflow mass and kinetic energy fluxes. Rotation also leads to time variability and fragmentation of the outflow in the radial and latitudinal directions. The collimation of the outflow is reduced in the models with gas rotation. The time variability in the mass and energy fluxes is reduced in the 3-D case because of the outflow fragmentation in the azimuthal direction. The virial mass estimated from the kinematics of the dense cold clouds found in our 3-D simulations of rotating gas underestimates the actual mass used in the simulations by about 40 %. (Abbreviated)

We discuss a number of topics relevant to disk-magnetosphere interaction and how numerical simulations illuminate them. The topics include: (1) disk-magnetosphere interaction and the problem of disk-locking; (2) the wind problem; (3) structure of the magnetospheric flow, hot spots at the star's surface, and the inner disk region; (4) modeling of spectra from 3D funnel streams; (5) accretion to a star with a complex magnetic field; (6) accretion through 3D instabilities; (7) magnetospheric gap and survival of protoplanets. Results of both 2D and 3D simulations are discussed.

We report on the fourth phase of our study of slightly rotating accretion flows onto black holes. The main new element of this study is that we used fully three dimensional (3-D) numerical simulations. We consider hydrodynamics of inviscid accretion flows. We assume a spherically symmetric density distribution at the outer boundary, but brake the flow symmetry by introducing a small, latitude-dependent angular momentum. We also consider cases where angular momentum at large radii is latitude- and azimuth-dependent. For the latitude-dependent angular momentum, 3-D simulations confirm axisymmetric results: the material that has too much angular momentum to be accreted forms a thick torus near the equator. Consequently, accretion proceeds only through the polar funnel, and the mass accretion rate through the funnel is constrained by the size and shape of the torus, not by the outer conditions. In 3-D simulations, we found that the torus precesses, even for axisymmetric conditions at large radii. For the latitude and azimuth-dependent angular momentum, the non-rotating gas near the equator can also significantly affect the evolution of the rotating gas. In particular, it may prevent the formation of a proper torus (i.e. its closing, in the azimuthal direction). In such models, the mass accretion rate is only slightly less than the corresponding Bondi rate.

We present three-dimensional (3-D) simulations of rotationally induced line variability arising from complex circumstellar environment of classical T Tauri stars (CTTS) using the results of the 3-D magnetohydrodynamic (MHD) simulations of Romanova et al., who considered accretion onto a CTTS with a misaligned dipole magnetic axis with respect to the rotational axis. The density, velocity and temperature structures of the MHD simulations are mapped on to the radiative transfer grid, and corresponding line source function and the observed profiles of neutral hydrogen lines (H-beta, Pa-beta and Br-gamma) are computed using the Sobolev escape probability method. We study the dependency of line variability on inclination angles (i) and magnetic axis misalignment angles (Theta). By comparing our models with the Pa-beta profiles of 42 CTTS observed by Folha & Emerson, we find that models with a smaller misaligngment angle (Theta<~15 deg.) are more consistent with the observations which show that majority of Pa-beta are rather symmetric around the line centre. For a high inclination system with a small dipole misalignment angle (Theta ~ 15 deg.), only one accretion funnel (on the upper hemisphere) is visible to an observer at any given rotational phase. This can cause an anti-correlation of the line equivalent width in the blue wing (v<0) and that in the red wing (v>0) over a half of a rotational period, and a positive correlation over other half. We find a good overall agreement of the line variability behaviour predicted by our model and those from observations. (Abridged)

We study the axisymmetric, time-dependent hydrodynamics of rotating flows that are under the influence of supermassive black hole gravity and radiation from an accretion disk surrounding the black hole. This work is an extension of the earlier work presented by Proga, where nonrotating flows were studied. Here, we consider effects of rotation, a position-dependent radiation temperature, density at large radii, and uniform X-ray background radiation. As in the non-rotating case, the rotating flow settles into a configuration with two components (1) an equatorial inflow and (2) a bipolar inflow/outflow with the outflow leaving the system along the pole. However, with rotation the flow does not always reach a steady state. In addition, rotation reduces the outflow collimation and the outward flux of mass and kinetic energy. Moreover rotation increases the outward flux of the thermal energy and can lead to fragmentation and time-variability of the outflow. We also show that a position-dependent radiation temperature can significantly change the flow solution. In particular, the inflow in the equatorial region can be replaced by a thermally driven outflow. Generally, as it have been discussed and shown in the past, we find that self-consistently determined preheating/cooling from the quasar radiation can significantly reduce the rate at which the central BH is fed with matter. However, our results emphasize also a little appreciated feature. Namely, quasar radiation drives a non-spherical, multi-temperature and very dynamic flow. These effects become dominant for luminosities in excess of 0.01 of the Eddington luminosity.

We present three-dimensional (3-D) hydrodynamical simulations of gas flows in the vicinity of an active galactic nucleus (AGN) powered by a precessing accretion disk. We consider the effects of the radiation force from such a disk on its environment on a relatively large scale (up to ~10 pc. We implicitly include the precessing disk by forcing the disk radiation field to precess around a symmetry axis with a given period ($P$) and a tilt angle ($\Theta$). We study time evolution of the flows irradiated by the disk, and investigate basic dependencies of the flow morphology, mass flux, angular momentum on different combinations of $\Theta$ and $P$. We find the gas flow settles into a configuration with two components, (1) an equatorial inflow and (2) a bipolar inflow/outflow with the outflow leaving the system along the poles (the directions of disk normals). However, the flow does not always reach a steady state. We find that the maximum outflow velocity and the kinetic outflow power at the outer boundary can be reduced significantly with increasing $\Theta$. We also find that of the mass inflow rate across the inner boundary does not change significantly with increasing $\Theta$. (Abbreviated)

We present the results of a radial velocity (RV) survey of 14 brown dwarfs (BDs) and very low-mass (VLM) stars in the Upper Scorpius OB association (UScoOB) and 3 BD candidates in the rho Ophiuchi dark cloud core. We obtained high-resolution echelle spectra at the Very Large Telescope using Ultraviolet and Visual Echelle Spectrograph (UVES) at two different epochs for each object, and measured the shifts in their RVs to identify candidates for binary/multiple systems in the sample. The average time separation of the RV measurements is 21.6d, and our survey is sensitive to the binaries with separation < 0.1 au. We found that 4 out of 17 objects (or 24^{+16}_{-13} per cent by fraction) show a significant RV change in 4-33d time scale, and are considered as binary/multiple `candidates.' We found no double-lined spectroscopic binaries in our sample, based on the shape of cross-correlation curves. The RV dispersion of the objects in UScoOB is found to be very similar to that of the BD and VLM stars in Chamaeleon I (Cha I). We also found the distribution of the mean rotational velocities (v sin i) of the UScoOB objects is similar to that of the Cha I, but the dispersion of v sin i is much larger than that of the Cha I objects.

We present radiative transfer models of the circumstellar environment of classical T Tauri stars, concentrating on the formation of the H-alpha emission. The wide variety of line profiles seen in observations are indicative of both inflow and outflow, and we therefore employ a circumstellar structure that includes both magnetospheric accretion and a disc wind. We perform systematic investigations of the model parameters for the wind and the magnetosphere to search for possible geometrical and physical conditions which lead to the types of profiles seen in observations. We find that the hybrid models can reproduce the wide range of profile types seen in observations, and that the most common profile types observed occupy a large volume of parameter space. Conversely, the most infrequently observed profile morphologies require a very specific set of models parameters. We find our model profiles are consistent with the canonical value of the mass-loss rate to mass-accretion rate ratio (mu=0.1) found in earlier magneto-hydrodynamic calculations and observations, but the models with 0.05<mu<0.2 are still in accord with observed H-alpha profiles. We investigate the wind contribution to the line profile as a function of model parameters, and examine the reliability of H-alpha as a mass accretion diagnostic. Finally, we examine the H-alpha spectroscopic classification used by Reipurth et. al, and discuss the basic physical conditions that are required to reproduce the profiles in each classified type.

We present measurements of the longitudinal magnetic field in the circumstellar environment of seven classical T Tauri stars. The measurements are based on high-resolution circular spectropolarimetry of the He I 5876 emission line, which is thought to form in accretion streams controlled by a stellar magnetosphere. We detect magnetic fields in BP Tau, DF Tau and DN Tau, and detect statistically significant fields in GM Aur and RW Aur A at one epoch but not at others. We detect no field for DG Tau and GG Tau, with the caveat that these objects were observed at one epoch only. Our measurements for BP Tau and DF Tau are consistent, both in terms of sign and magnitude, with previous studies, suggesting that the characteristics of T Tauri magnetospheres are persistent over several years. We observed the magnetic field of BP Tau to decline monotonically over three nights, and have detected a peak field of 4kG in this object, the highest magnetic field yet observed in a T Tauri star. We combine our observations with results from the literature in order to perform a statistical analysis of the magnetospheric fields in BP Tau and DF Tau. Assuming a dipolar field, we determine a polar field of ~3kG and a dipole offset of 40deg for BP Tau, while DF Tau's field is consistent with a polar field of ~-4.5kG and a dipole offset of 10deg. We conclude that many classical T Tauri stars have circumstellar magnetic fields that are both strong enough and sufficiently globally-ordered to sustain large-scale magnetospheric accretion flows.

We present time-series echelle spectra of the Pa-beta line of the T Tauri star SU Aur, observed over three consecutive nights. The line shows strong variability (~10 per cent) over the velocity range (100 km/s, 420 km/s) in the red broad absorption component, and weaker variability (~2 per cent) over the velocity range (-200 km/s, 0 km/s) in the blue wing. The variability in the velocity range (-200 km/s, 0 km/s) is correlated with that in (200 km/s, 400 km/s), and the variability in these velocity ranges anti-correlates with that in (0 km/s, 100 km/s). The mean spectrum from the second night shows a suggestion of a blue-shifted absorption component at about -150 km/s, similar to that found in the H-alpha and H-beta lines. We find the position of the subpeak in the red absorption component changes steadily with time, and its motion modulates on half the rotational period. We also find that the modulation of the line equivalent width is associated with a half and a third of the rotational period, which is consistent with the surface Doppler images of SU Aur. Radiative transfer models of a rotationally modulated Pa-beta line, produced in the shock-heated magnetospheric accretion flow, are also presented. Models with a magnetic dipole offset reproduce the overall characteristics of the observed line variability, including the line equivalent width and the motion of the subpeak in the red absorption trough.

We present hydrogen emission line profile models of magnetospheric accretion onto Classical T Tauri stars. The models are computed under the Sobolev approximation using the three-dimensional Monte Carlo radiative-transfer code TORUS. We have calculated four illustrative models in which the accretion flows are confined to azimuthal curtains - a geometry predicted by magneto-hydrodynamical simulations. Properties of the line profile variability of our models are discussed, with reference to dynamic spectra and cross-correlation images. We find that some gross characteristics of observed line profile variability are reproduced by our models, although in general the level of variability predicted is larger than that observed. We conclude that this excessive variability probably excludes dynamical simulations that predict accretion flows with low degrees of axisymmetry.

The aim of this article is to demonstrate the useful role that can be played by spectropolarimetric observations of young and evolved emission line stars that analyse the linearly polarized component in their spectra. At the time of writing, this demonstration has to be made on the basis of optical data since there is no common-user infrared facility, in operation, that offers the desired combination of spectral resolution and sensitivity. Here we focus on what can be learned from linear spectropolarimetry alone at reasonably high spectral resolution and at $10^3 < $S/N$ < 10^4$. And we remind that the near infrared (1--2 micron) has the potential to out-perform the optical as a domain to work in because of the greatly reduced interstellar obscuration at these wavelengths. This point has been reached at a time when theory, exploiting flexible Monte Carlo methods, is fast becoming a powerful tool. In short we have the complex phenomena, and the rise of the modelling capability to match -- good data are the missing link.

We present three-dimensional Monte Carlo radiative transfer models of a very young (<10^5 years old) low mass (50 Msun) stellar cluster containing 23 stars and 27 brown dwarfs. The models use the density and the stellar mass distributions from the large-scale smoothed particle hydrodynamics (SPH) simulation of the formation of a low-mass stellar cluster by Bate, Bonnell and Bromm. We compute the spectral energy distributions of the whole cluster and the individual objects, simulated far-infrared Spitzer Space Telescope (SST) images and construct colour-colour diagrams (near-infrared HKL and SST mid-infrared bands). The presence of accretion discs around the light sources influences the morphology of the dust temperature structure on a large scale (up to a several 10^4 au). A considerable fraction of the interstellar dust is underheated compared to a model without the accretion discs because the radiation from the light sources is blocked/shadowed by the discs. The spectral energy distribution (SED) of the model cluster with accretion discs shows excess emission at lambda = 3-30 micron and lambda > 500 micron compared to that without accretion discs. Our model with accretion discs around each object shows a similar distribution of spectral index (2.2-20 micron) values (i.e. Class 0--III sources) as seen in the Rho Ophiuchus cloud. We confirm that the best diagnostics for identifying objects with accretion discs are mid-infrared (lambda = 3-10 micron) colours (e.g. SST IRAC bands) rather than HKL colours.

We present radiative-transfer modelling of the dusty spiral Pinwheel Nebula observed around the Wolf-Rayet/OB-star binary WR104. The models are based on the three-dimensional radiative-transfer code TORUS, modified to include an adaptive mesh that allows us to adequately resolve both the inner spiral turns (sub-AU scales) and the outer regions of the nebula (distances of 10^4 AU from the central source). The spiral model provides a good fit to both the spectral energy distribution and Keck aperture masking interferometry, reproducing both the maximum entropy recovered images and the visibility curves. We deduce a dust creation rate of 8+-1 x 10^{-7} solar masses per year, corresponding to approximately 2% by mass of the carbon produced by the Wolf-Rayet star. Simultaneous modelling of the imaging and spectral data enables us to constrain both the opening-angle of the wind-wind collision interface and the dust grain size. We conclude that the dust grains in the inner part of the Pinwheel nebula are small (~100A), in agreement with theoretical predictions, although we cannot rule out the presence of larger grains (~1 micron) further from the central binary. The opening angle of the wind-wind collision interface appears to be about 40 degrees, in broad agreement with the wind parameters estimated for the central binary. We discuss the success and deficiencies of the model, and the likely benefits of applying similar techniques to the more the more complex nebulae observed around other WR/O star binaries.

A newly developed 3-D Monte Carlo model is used, in conjunction with a multi-line non-LTE radiative transfer model, to determine the mass-loss rate of the Wolf-Rayet (W-R) star in the massive binary \object{V444 Cyg} (WN5+O6). This independent estimate of mass-loss rate is attained by fitting the observed \HeI (5876) \AA and \HeII (5412) \AA line profiles, and the continuum light curves of three Stokes parameters ((I, Q, U)) in the (V) band simultaneously. The high accuracy of our determination arises from the use of many observational constraints, and the sensitivity of the continuum polarization to the mass-loss rate. Our best fit model suggests that the mass-loss rate of the system is (\dot{M}_{\WR}=0.6(\pm 0.2) \times 10^{-5} M_{\sun} \mathrm{yr}^{-1} ), and is independent of the assumed distance to \object{V444 Cyg}. The fits did not allow a unique value for the radius of the W-R star to be derived. The range of the volume filling factor for the W-R star atmosphere is estimated to be in the range of 0.050 (for $R_{\WR}=5.0 R_{\sun}$) to 0.075 (for $R_{\WR}=2.5 R_{\sun}$). We also found that the blue-side of \HeI (5876 ) \AA and \HeII (5412) \AA lines at phase 0.8 is relatively unaffected by the emission from the wind-wind interaction zone and the absorption by the O-star atmosphere; hence, the profiles at this phase are suitable for spectral line fittings using a spherical radiative transfer model.

We have developed a 3-D Monte Carlo radiative transfer model which computes line and continuum polarization variability for a binary system with an optically thick non-axisymmetric envelope. This allows us to investigate the complex (phase-locked) line and continuum polarization variability features displayed by many massive binaries: W-R+O, O+O, etc. An 8-way tree data structure constructed via a ``cell-splitting'' method allows for high precision with efficient use of computer resources. The model is not restricted to binary systems; it can easily be adapted to a system with an arbitrary density distribution and large density gradients. As an application to a real system, the phase dependent Stokes parameters (I, Q, U) and the phase dependent He I (5876) profiles of the massive binary system V444 Cyg (WN5+O6 III-V) are computed.