Showing 1–50 of 110 results for author: Romanova, M

Eccentricity and Inclination of Massive Planets Inside Low-density Cavities: Results of 3D Simulations

Authors: M. M. Romanova, A. V. Koldoba, G. V. Ustyugova, C. Espaillat, R. V. E. Lovelace

Abstract: We study the evolution of eccentricity and inclination of massive planets in low-density cavities of protoplanetary discs using three-dimensional (3D) simulations. When the planet's orbit is aligned with the equatorial plane of the disc, the eccentricity increases to high values of 0.7-0.9 due to the resonant interaction with the inner parts of the disc. For planets on inclined orbits, the eccentr… ▽ More We study the evolution of eccentricity and inclination of massive planets in low-density cavities of protoplanetary discs using three-dimensional (3D) simulations. When the planet's orbit is aligned with the equatorial plane of the disc, the eccentricity increases to high values of 0.7-0.9 due to the resonant interaction with the inner parts of the disc. For planets on inclined orbits, the eccentricity increases due to the Kozai-Lidov mechanism, where the disc acts as an external massive body that perturbs the planet's orbit. At small inclination angles, < 30 degrees, the resonant interaction with the inner disc strongly contributes to the eccentricity growth, while at larger angles, eccentricity growth is mainly due to the Kozai-Lidov mechanism. We conclude that planets inside low-density cavities tend to acquire high eccentricity if favorable conditions give sufficient time for growth. The final value of the planet's eccentricity, after the disc dispersal depends on the planet's mass and properties of the cavity and protoplanetary disc. △ Less

Submitted 26 June, 2024; originally announced June 2024.

Comments: 19 pages, 17 figures, submitted to MNRAS

Accretion Funnel Reconfiguration during an Outburst in a Young Stellar Object: EX Lupi

Authors: Koshvendra Singh, Joe P. Ninan, Marina M. Romanova, David A. H. Buckley, Devendra K. Ojha, Arpan Ghosh, Andrew Monson, Malte Schramm, Saurabh Sharma, Daniel E. Reichart, Joanna Mikolajewska, Juan Carlos Beamin, J. Borissova, Valentin D. Ivanov, Vladimir V. Kouprianov, Franz-Josef Hambsch, Andrew Pearce

Abstract: EX Lupi, a low-mass young stellar object, went into an accretion-driven outburst in March of 2022. The outburst caused a sudden phase change of ~ 112$^{\circ}$ $\pm$ 5$^{\circ}$ in periodically oscillating multiband lightcurves. Our high resolution spectra obtained with HRS on SALT also revealed a consistent phase change in the periodically varying radial velocities, along with an increase in the… ▽ More EX Lupi, a low-mass young stellar object, went into an accretion-driven outburst in March of 2022. The outburst caused a sudden phase change of ~ 112$^{\circ}$ $\pm$ 5$^{\circ}$ in periodically oscillating multiband lightcurves. Our high resolution spectra obtained with HRS on SALT also revealed a consistent phase change in the periodically varying radial velocities, along with an increase in the radial velocity amplitude of various emission lines. The phase change and increase of radial velocity amplitude morphologically translates to a change in the azimuthal and latitudinal location of the accretion hotspot over the stellar surface, which indicates a reconfiguration of the accretion funnel geometry. Our 3D MHD simulations reproduce the phase change for EX Lupi. To explain the observations we explored the possibility of forward shifting of the dipolar accretion funnel as well as the possibility of an emergence of a new accretion funnel. During the outburst, we also found evidence of the hotspot's morphology extending azimuthally, asymmetrically with a leading hot edge and cold tail along the stellar rotation. Our high cadence photometry showed that the accretion flow has clumps. We also detected possible clumpy accretion events in the HRS spectra, that showed episodically highly blue-shifted wings in the Ca II IRT and Balmer H lines. △ Less

Submitted 8 April, 2024; originally announced April 2024.

Comments: Accepted for publication in The Astrophysical Journal

Density streams in the disc winds of Classical T Tauri stars

Authors: P. P. Petrov, K. N. Grankin, E. V. Babina, S. A. Artemenko, M. M. Romanova, S. Yu. Gorda, A. A. Djupvik, J. F. Gameiro

Abstract: Spectral and photometric variability of the Classical T Tauri stars RY Tau and SU Aur from 2013 to 2022 is analyzed. We find that in SU Aur the H-alpha line's flux at radial velocity RV = -50 +- 7 km/s varies with a period P = 255 +- 5 days. A similar effect previously discovered in RY Tau is confirmed with these new data: P = 21.6 days at RV = -95 +- 5 km/s. In both stars, the radial velocity of… ▽ More Spectral and photometric variability of the Classical T Tauri stars RY Tau and SU Aur from 2013 to 2022 is analyzed. We find that in SU Aur the H-alpha line's flux at radial velocity RV = -50 +- 7 km/s varies with a period P = 255 +- 5 days. A similar effect previously discovered in RY Tau is confirmed with these new data: P = 21.6 days at RV = -95 +- 5 km/s. In both stars, the radial velocity of these variations, the period, and the mass of the star turn out to be related by Kepler's law, suggesting structural features on the disc plane orbiting at radii of 0.2 AU in RY Tau and 0.9 AU in SU Aur, respectively. Both stars have a large inclination of the accretion disc to the line of sight - so that the line of sight passes through the region of the disc wind. We propose there is an azimuthal asymmetry in the disc wind, presumably in the form of 'density streams', caused by substructures of the accretion disc surface. These streams cannot dissipate until they go beyond the Alfven surface in the disc's magnetic field. These findings open up the possibility to learn about the structure of the inner accretion disc of CTTS on scales less than 1 AU and to reveal the orbital distances related to the planet's formation. △ Less

Submitted 17 August, 2023; originally announced August 2023.

Comments: 9 pages, 10 figures

Journal ref: Monthly Notices of the Royal Astronomical Society, 2023, V.524, pp.5944-5953

Eccentricity Growth of Massive Planets inside Cavities of Protoplanetary Discs

Authors: Marina M. Romanova, Alexander V. Koldoba, Galina V. Ustyugova, Dong Lai, Richard V. E. Lovelace

Abstract: We carry out hydrodynamical simulations to study the eccentricity growth of a 1-30 Jupiter mass planet located inside the fixed cavity of a protoplanetary disc. The planet exchanges energy and angular momentum with the disc at resonant locations, and its eccentricity grows due to Lindblad resonances. We observe several phases of eccentricity growth where different eccentric Lindblad resonances dom… ▽ More We carry out hydrodynamical simulations to study the eccentricity growth of a 1-30 Jupiter mass planet located inside the fixed cavity of a protoplanetary disc. The planet exchanges energy and angular momentum with the disc at resonant locations, and its eccentricity grows due to Lindblad resonances. We observe several phases of eccentricity growth where different eccentric Lindblad resonances dominate from 1:3 up to 3:5. The maximum values of eccentricity reached in our simulations are 0.65-0.75. We calculate the eccentricity growth rate for different planet masses and disc parameters and derive analytical dependencies on these parameters. We observe that the growth rate is proportional to both the planet's mass and the characteristic disc mass for a wide range of parameters. In a separate set of simulations, we derived the width of the 1:3 Lindblad resonance. △ Less

Submitted 30 March, 2023; originally announced March 2023.

Comments: 19 pages, 18 figures, accepted by MNRAS

Measuring the Density Structure of an Accretion Hot Spot

Authors: C. C. Espaillat, C. E. Robinson, M. M. Romanova, T. Thanathibodee, J. Wendeborn, N. Calvet, M. Reynolds, J. Muzerolle

Abstract: Magnetospheric accretion models predict that matter from protoplanetary disks accretes onto the star via funnel flows which follow the stellar field lines and shock on the stellar surface leaving a hot spot with a density gradient. Previous work has inferred different densities in the hot spot, but has not been sensitive to the radial density distribution. Attempts have been made to measure this w… ▽ More Magnetospheric accretion models predict that matter from protoplanetary disks accretes onto the star via funnel flows which follow the stellar field lines and shock on the stellar surface leaving a hot spot with a density gradient. Previous work has inferred different densities in the hot spot, but has not been sensitive to the radial density distribution. Attempts have been made to measure this with X-ray observations, but X-ray emission only traces a fraction of the hot spot and also coronal emission. Here we report periodic ultraviolet and optical light curves of the accreting star GM Aur that display a time lag of about 1 day between their peaks. The periodicity arises as the source of the ultraviolet and optical emission moves into and out of view as it rotates along with the star. The time-lag indicates a difference in the spatial distribution of ultraviolet and optical brightness over the stellar surface. Within the framework of the magnetospheric accretion model, this indicates a radial density gradient in a hot spot on the stellar surface since different density parts of the hot spot are expected to emit radiation at different wavelengths. These results are the first observational confirmation of the magnetospheric accretion model's prediction of a density gradient in the hot spot and demonstrate the insights gained from focusing on the wavelengths where the bulk of the accretion energy can be observed. △ Less

Submitted 1 September, 2021; originally announced September 2021.

Comments: published in Nature; open access link: https://www.nature.com/articles/s41586-021-03751-5

Journal ref: Nature volume 597, pages 41-44 (2021)

Modulated accretion in T Tauri star RY Tau -- stable MHD propeller or a planet at 0.2 AU?

Authors: P. P. Petrov, M. M. Romanova, K. N. Grankin, S. A. Artemenko, E. V. Babina, S. Yu. Gorda

Abstract: Planets are thought to form at the early stage of stellar evolution when the mass accretion is still ongoing. RY Tau is a T Tauri type star at the age of a few Myr, with accretion disc seen at high inclination, so that line of sight crosses both the wind and the accretion gas flows. In a long series of spectroscopic monitoring of the star in 2013-2020, we detected variations in H-alpha and NaI D a… ▽ More Planets are thought to form at the early stage of stellar evolution when the mass accretion is still ongoing. RY Tau is a T Tauri type star at the age of a few Myr, with accretion disc seen at high inclination, so that line of sight crosses both the wind and the accretion gas flows. In a long series of spectroscopic monitoring of the star in 2013-2020, we detected variations in H-alpha and NaI D absorptions at radial velocities of infall (accretion) and outflow (wind) with a period of about 22 days. The absorptions in the infalling and the outflowing gas streams vary in anti-phase: an increase of infall is accompanied by a decrease of outflow, and vice versa. These flip-flop oscillations retain phase over several years of observations. We suggest that this may result from the MHD processes at the disk-magnetosphere boundary in the propeller mode. Another possibility is that a massive planet modulates some processes in the disc and provides the observed effects. The period, if Keplerian, corresponds to a distance of 0.2 AU, which is close to the dust sublimation radius in this star. The presence of the putative planet may be confirmed by radial velocity measurements: expected amplitude is > 90 m/s if a planet mass is > 2 Mj. △ Less

Submitted 1 April, 2021; originally announced April 2021.

Comments: 7 pages, 12 figures, accepted for publication in MNRAS

3D MHD Simulations of Accretion onto Stars with Tilted Magnetic and Rotational Axes

Authors: M. M. Romanova, A. V. Koldoba, G. V. Ustyugova, A. A. Blinova, D. Lai, R. V. E. Lovelace

Abstract: We present results of global three-dimensional (3D) magnetohydrodynamic (MHD) simulations of accretion onto magnetized stars where both the magnetic and rotational axes of the star are tilted about the rotational axis of the disc. We observed that initially the inner parts of the disc are warped, tilted, and process due to the magnetic interaction between the magnetosphere and the disc. Later, lar… ▽ More We present results of global three-dimensional (3D) magnetohydrodynamic (MHD) simulations of accretion onto magnetized stars where both the magnetic and rotational axes of the star are tilted about the rotational axis of the disc. We observed that initially the inner parts of the disc are warped, tilted, and process due to the magnetic interaction between the magnetosphere and the disc. Later, larger tilted discs form with the size increasing with the magnetic moment of the star. The normal vector to the discs are tilted at different angles, from 5-10 degrees up to 30-40 degrees. Small tilts may result from the winding of the magnetic field lines about the rotational axis of the star and the action of the magnetic force which tends to align the disc. Another possible explanation is the magnetic Bardeen-Petterson effect in which the disc settles in the equatorial plane of the star due to precessional and viscous torques in the disc. Tilted discs slowly precess with the time scale of the order of 50 Keplerian periods at the reference radius (approx. 3 stellar radii). Our results can be applied to different types of stars where signs of tilted discs and/or slow precession have been observed. △ Less

Submitted 4 July, 2021; v1 submitted 19 December, 2020; originally announced December 2020.

Comments: 14 pages, 13 figures, MNRAS, accepted. The new version of the paper is more complete compared with the earlier version

3D Simulations of Planet Trapping at Disc-Cavity Boundaries

Authors: M. M. Romanova, P. S. Lii, A. V. Koldoba, G. V. Ustyugova, A. A. Blinova, R. V. E. Lovelace, L. Kaltenegger

Abstract: Inward migration of low-mass planets and embryos of giant planets can be stopped at the disc-cavity boundaries due to co-orbital corotation torque. We performed the first global three-dimensional (3D) simulations of planet migration at the disc-cavity boundary, and have shown that the boundary is a robust trap for low-mass planets and embryos. A protoplanetary disc may have several such trapping r… ▽ More Inward migration of low-mass planets and embryos of giant planets can be stopped at the disc-cavity boundaries due to co-orbital corotation torque. We performed the first global three-dimensional (3D) simulations of planet migration at the disc-cavity boundary, and have shown that the boundary is a robust trap for low-mass planets and embryos. A protoplanetary disc may have several such trapping regions at various distances from the star, such as at the edge of the stellar magnetosphere, the inner edge of the dead zone, the dust-sublimation radius and the snow lines. Corotation traps located at different distances from a star, and moving outward during the disc dispersal phase, may possibly explain the observed homogeneous distribution of low-mass planets with distance from their host stars. △ Less

Submitted 22 February, 2019; v1 submitted 11 September, 2018; originally announced September 2018.

Comments: 17 pages, 14 figures, accepted by MNRAS

Comparisons of MHD Propeller Model with Observations of Cataclysmic Variable AE Aqr

Authors: A. A. Blinova, M. M. Romanova, G. V. Ustyugova, A. V. Koldoba, R. V. E. Lovelace

Abstract: We have developed a numerical MHD model of the propeller candidate star AE Aqr using axisymmetric magneto-hydrodynamic (MHD) simulations. We suggest that AE Aqr is an intermediate polar-type star, where the magnetic field is relatively weak and an accretion disc may form around the white dwarf. The star is in the propeller regime, and many of its observational properties are determined by the disc… ▽ More We have developed a numerical MHD model of the propeller candidate star AE Aqr using axisymmetric magneto-hydrodynamic (MHD) simulations. We suggest that AE Aqr is an intermediate polar-type star, where the magnetic field is relatively weak and an accretion disc may form around the white dwarf. The star is in the propeller regime, and many of its observational properties are determined by the disc-magnetosphere interaction. Comparisons of the characteristics of the observed versus modelled AE Aqr star show that the model can explain many observational properties of AE Aqr. In a representative model, the magnetic field of the star is B\approx 3.3x10^5 G and the time averaged accretion rate in the disc is 5.5\times 10^{16} g/s. Most of this matter is ejected into conically-shaped winds. The numerical model explains the rapid spin-down of AE Aqr through the outflow of angular momentum from the surface of the star to the wind, corona and disc. The energy budget in the outflows, 9x10^{33} erg/s, is sufficient for explaining the observed flaring radiation in different wavebands. The time scale of ejections into the wind matches the short time scale variability in the light curves of AE Aqr. △ Less

Submitted 24 April, 2019; v1 submitted 26 April, 2018; originally announced April 2018.

Comments: 11 pages, 4 figures, accepted by MNRAS

Modelling the bow shock Pulsar Wind Nebulae propagating through a non-uniform ISM

Authors: O. D. Toropina, M. M. Romanova, R. V. E. Lovelace

Abstract: Many pulsars propagate through the interstellar medium (ISM) with supersonic velocities, and their pulsar winds interact with the interstellar medium (ISM), forming bow shocks and magnetotails (PWN). We model the propagation of pulsars through the inhomogeneous ISM using non-relativistic axisymmetric magneto-hydrodynamic (MHD) simulations. We take into account the wind from the star, and the azimu… ▽ More Many pulsars propagate through the interstellar medium (ISM) with supersonic velocities, and their pulsar winds interact with the interstellar medium (ISM), forming bow shocks and magnetotails (PWN). We model the propagation of pulsars through the inhomogeneous ISM using non-relativistic axisymmetric magneto-hydrodynamic (MHD) simulations. We take into account the wind from the star, and the azimuthal and poloidal components of the magnetic field, and investigate the PWN at different levels of magnetization (the ratio of magnetic to matter energy-densities) in the wind. We consider the interaction of PWN with small-scale and large-scale imhomogeneities in the ISM at different values of magnetization. We conclude that the inhomogeneities in the ISM can change the shapes of the bow shocks and magnetotails at different values of the magnetization.We compare the results of our simulations with the images of the Guitar Nebula and other PWN that show irregularities in the shapes of their bow shocks and magnetotails. We conclude that these irregularities may be caused by the interaction of PWN with the inhomogeneities in the ISM. △ Less

Submitted 15 March, 2018; originally announced March 2018.

Comments: 11 pages, 12 figures, submitted to MNRAS

Properties of Strong and Weak Propellers from MHD Simulations

Authors: Marina M. Romanova, Alisa A. Blinova, Galina V. Ustyugova, Alexander V. Koldoba, Richard V. E. Lovelace

Abstract: We investigate the properties of magnetized stars in the propeller regime using axisymmetric numerical simulations. We modelled the propeller regime for stars with realistically large magnetospheres (5-7 stellar radii), so that our results could be applied to different types of magnetized stars, including Classical T Tauri stars, cataclysmic variables, and accreting millisecond pulsars. A wide ran… ▽ More We investigate the properties of magnetized stars in the propeller regime using axisymmetric numerical simulations. We modelled the propeller regime for stars with realistically large magnetospheres (5-7 stellar radii), so that our results could be applied to different types of magnetized stars, including Classical T Tauri stars, cataclysmic variables, and accreting millisecond pulsars. A wide range of propeller strengths has been studied, from very strong propellers (where the magnetosphere rotates much more rapidly than the inner disk) to very weak propellers (where the magnetosphere rotates only slightly faster than the inner disc. In both propellers, matter is accumulated at the inner disc for the majority of the time, while episodes of accretion onto the star and ejection into the wind are relatively brief. The efficiency of the propeller, which characterizes the part of inner disk matter flowing into the wind, strongly depends on the fastness parameter omega_s: propeller efficiency increases with omega_s. The properties of the winds are different in strong and weak propellers. In the strong propellers, matter is accelerated rapidly above the escape velocity and flows at a relatively small opening angle of 40-45 degrees. In the weak propellers, matter may flow faster or slower than the escape velocity and at a large opening angle of 60-70 degrees. A star-disk system loses energy and angular momentum. A part of the rotational energy of the star is ejected to the magnetically-dominated (Poynting flux) jet, which is only present in the strong propellers. The other part of the energy flows from the inner disk into a propeller-driven wind. A star spins down partly due to the flow of angular momentum from the star to the corona or to the Poynting flux jet along the open field lines, and partly due to the flow of angular momentum to the inner disk along the closed field lines. △ Less

Submitted 26 April, 2017; originally announced April 2017.

Comments: 31 pages, 23 figures, submitted to New Astronomy

Accretion, Outflows, and Winds of Magnetized Stars

Authors: M. M. Romanova, S. P. Owocki

Abstract: Many types of stars have strong magnetic fields that can dynamically influence the flow of circumstellar matter. In stars with accretion disks, the stellar magnetic field can truncate the inner disk and determine the paths that matter can take to flow onto the star. These paths are different in stars with different magnetospheres and periods of rotation. External field lines of the magnetosphere m… ▽ More Many types of stars have strong magnetic fields that can dynamically influence the flow of circumstellar matter. In stars with accretion disks, the stellar magnetic field can truncate the inner disk and determine the paths that matter can take to flow onto the star. These paths are different in stars with different magnetospheres and periods of rotation. External field lines of the magnetosphere may inflate and produce favorable conditions for outflows from the disk-magnetosphere boundary. Outflows can be particularly strong in the propeller regime, wherein a star rotates more rapidly than the inner disk. Outflows may also form at the disk-magnetosphere boundary of slowly rotating stars, if the magnetosphere is compressed by the accreting matter. In isolated, strongly magnetized stars, the magnetic field can influence formation and/or propagation of stellar wind outflows. Winds from low-mass, solar-type stars may be either thermally or magnetically driven, while winds from massive, luminous O and B type stars are radiatively driven. In all of these cases, the magnetic field influences matter flow from the stars and determines many observational properties. In this chapter we review recent studies of accretion, outflows, and winds of magnetized stars with a focus on three main topics: (1) accretion onto magnetized stars; (2) outflows from the disk-magnetosphere boundary; and (3) winds from isolated massive magnetized stars. We show results obtained from global magnetohydrodynamic simulations and, in a number of cases compare global simulations with observations. △ Less

Submitted 16 May, 2016; originally announced May 2016.

Comments: 60 pages, 44 figures

Journal ref: Space Science Reviews, Volume 191, Issue 1-4, pp. 339-389, 2015

The effects of a magnetic field on planetary migration in laminar and turbulent discs

Authors: M. L. Comins, M. M. Romanova, A. V. Koldoba, G. V. Ustyugova, R. V. E. Lovelace

Abstract: We investigate the migration of low-mass planets ($5 M_{\oplus}$ and $20 M_{\oplus}$) in accretion discs threaded with a magnetic field using 2D MHD code in polar coordinates. We observed that, in the case of a strong azimuthal magnetic field where the plasma parameter is $β\sim 1-2$, density waves at the magnetic resonances exert a positive torque on the planet and may slow down or reverse its mi… ▽ More We investigate the migration of low-mass planets ($5 M_{\oplus}$ and $20 M_{\oplus}$) in accretion discs threaded with a magnetic field using 2D MHD code in polar coordinates. We observed that, in the case of a strong azimuthal magnetic field where the plasma parameter is $β\sim 1-2$, density waves at the magnetic resonances exert a positive torque on the planet and may slow down or reverse its migration. However, when the magnetic field is weaker (i.e., the plasma parameter $β$ is relatively large), then non-axisymmetric density waves excited by the planet lead to growth of the radial component of the field and, subsequently, to development of the magneto-rotational instability, such that the disc becomes turbulent. Migration in a turbulent disc is stochastic, and the migration direction may change as such. To understand migration in a turbulent disc, both the interaction between a planet and individual turbulent cells, as well as the interaction between a planet and ordered density waves, have been investigated. △ Less

Submitted 13 August, 2015; originally announced August 2015.

Comments: 16 pages, 12 figures

Numerical MHD Codes for Modeling Astrophysical Flows

Authors: A. V. Koldoba, G. V. Ustyugova, P. S. Lii, M. L. Comins, S. Dyda, M. M. Romanova, R. V. E. Lovelace

Abstract: We describe a Godunov-type magnetohydrodynamic (MHD) code based on the Miyoshi and Kusano (2005) solver which can be used to solve various astrophysical hydrodynamic and MHD problems. The energy equation is in the form of entropy conservation. The code has been implemented on several different coordinate systems: 2.5D axisymmetric cylindrical coordinates, 2D Cartesian coordinates, 2D plane polar c… ▽ More We describe a Godunov-type magnetohydrodynamic (MHD) code based on the Miyoshi and Kusano (2005) solver which can be used to solve various astrophysical hydrodynamic and MHD problems. The energy equation is in the form of entropy conservation. The code has been implemented on several different coordinate systems: 2.5D axisymmetric cylindrical coordinates, 2D Cartesian coordinates, 2D plane polar coordinates, and fully 3D cylindrical coordinates. Viscosity and diffusivity are implemented in the code to control the accretion rate in the disk and the rate of penetration of the disk matter through the magnetic field lines. The code has been utilized for the numerical investigations of a number of different astrophysical problems, several examples of which are shown. △ Less

Submitted 31 July, 2015; originally announced August 2015.

Comments: 23 pages, 17 figures, submitted to New Astronomy

Asymmetric MHD Outflows/Jets from Accreting T Tauri Stars

Authors: Sergei Dyda, Richard V. E. Lovelace, Galina V. Ustyugova, Patrick S. Lii, Marina M. Romanova, Alexander V. Koldoba

Abstract: Observations of jets from young stellar objects reveal the asymmetric outflows from some sources. A large set of $2.5$D MHD simulations has been carried out for axisymmetric viscous/diffusive disc accretion to rotating magnetized stars for the purpose of assessing the conditions where the outflows or jets are asymmetric relative to the equatorial plane. We consider initial magnetic fields that are… ▽ More Observations of jets from young stellar objects reveal the asymmetric outflows from some sources. A large set of $2.5$D MHD simulations has been carried out for axisymmetric viscous/diffusive disc accretion to rotating magnetized stars for the purpose of assessing the conditions where the outflows or jets are asymmetric relative to the equatorial plane. We consider initial magnetic fields that are symmetric about the equatorial plane and consist of a radially distributed field threading the disc (disc-field) and a stellar dipole field.({\bf 1}). For pure disc-fields the symmetry or asymmetry of the outflows is affected by the midplane plasma $β$ of the disc (where $β$ is the ratio of the plasma pressure to the magnetic pressure). For the low density discs with small plasma $β$ values, outflows are observed to be symmetric to within $10\%$ over timescales of hundreds of inner disc orbits. For the denser higher $β$ discs, the coupling of the upper and lower coronal plasmas is broken, and quasi-periodic field motion in the two hemispheres becomes different. This asymmetry leads to asymmetric episodic outflows. ({\bf 2.}) Accreting stars with a stellar dipole field and no disc-field exhibit episodic, two component outflows - a magnetospheric wind and an inner disc wind from somewhat larger radial distances. Both are characterized by similar velocity profiles but the magnetospheric wind has densities $\gtrsim 10$ times that of the disc wind. ({\bf 3}.)Adding a disc-field parallel to the stellar dipole field acts to enhance the magnetospheric winds but suppress the disc wind. ({\bf 4}.) In contrast, adding a disc-field which is anti-parallel to the stellar dipole field in the disc acts to suppress the magnetospheric and disc winds. Our simulations reproduce some key features of observations of asymmetric outflows of T Tauri stars. △ Less

Submitted 16 March, 2015; v1 submitted 4 February, 2015; originally announced February 2015.

Comments: 16 pages, 17 figures, v2 has additional discussion regarding observations

Boundary between Stable and Unstable Regimes of Accretion. Ordered and Chaotic Unstable Regimes

Authors: A. A. Blinova, M. M. Romanova, R. V. E. Lovelace

Abstract: We present a study of the Rayleigh-Taylor unstable regime of accretion onto rotating magnetized stars in a set of high grid resolution three-dimensional (3D) magnetohydrodynamic (MHD) simulations performed in low-viscosity discs. We find that the boundary between the stable and unstable regimes is determined almost entirely by the fastness parameter omega_s=Omega_star/Omega_K(r_m), where Omega_sta… ▽ More We present a study of the Rayleigh-Taylor unstable regime of accretion onto rotating magnetized stars in a set of high grid resolution three-dimensional (3D) magnetohydrodynamic (MHD) simulations performed in low-viscosity discs. We find that the boundary between the stable and unstable regimes is determined almost entirely by the fastness parameter omega_s=Omega_star/Omega_K(r_m), where Omega_star is the angular velocity of the star and Omega_K(r_m) is the angular velocity of the Keplerian disc at the disc-magnetosphere boundary r=r_m. We found that accretion is unstable if omega_s < 0.6. Accretion through instabilities is present in stars with different magnetospheric sizes. However, only in stars with relatively small magnetospheres, r_m/R_star < 7, do the unstable tongues produce chaotic hot spots on the stellar surface and irregular light-curves. At even smaller values of the fastness parameter, omega_s < 0.45, multiple irregular tongues merge, forming one or two ordered unstable tongues that rotate with the angular frequency of the inner disc. This transition occurs in stars with even smaller magnetospheres, r_m/R_star < 4.2. Most of our simulations were performed at a small tilt of the dipole magnetosphere, Theta=5 degrees, and a small viscosity parameter alpha=0.02. Test simulations at higher alpha values show that many more cases become unstable, and the light-curves become even more irregular. Test simulations at larger tilts of the dipole Theta show that instability is present, however, accretion in two funnel streams dominates if Theta > 15 degrees. The results of these simulations can be applied to accreting magnetized stars with relatively small magnetospheres: Classical T Tauri stars, accreting millisecond X-ray pulsars, and cataclysmics variables. △ Less

Submitted 16 May, 2016; v1 submitted 8 January, 2015; originally announced January 2015.

Comments: 17 pages, 14 figures, published in MNRAS

Journal ref: MNRAS, Volume 459, p.2354, 2016

Counter-Rotating Accretion Discs

Authors: Sergei Dyda, Richard V. E. Lovelace, Galina V. Ustyugova, Marina M. Romanova, Alexander V. Koldoba

Abstract: Counter-rotating discs can arise from the accretion of a counter-rotating gas cloud onto the surface of an existing co-rotating disc or from the counter-rotating gas moving radially inward to the outer edge of an existing disc. At the interface, the two components mix to produce gas or plasma with zero net angular momentum which tends to free-fall towards the disc center. We discuss high-resolutio… ▽ More Counter-rotating discs can arise from the accretion of a counter-rotating gas cloud onto the surface of an existing co-rotating disc or from the counter-rotating gas moving radially inward to the outer edge of an existing disc. At the interface, the two components mix to produce gas or plasma with zero net angular momentum which tends to free-fall towards the disc center. We discuss high-resolution axisymmetric hydrodynamic simulations of a viscous counter-rotating disc for cases where the two components are vertically separated and radially separated. The viscosity is described by an isotropic $α-$viscosity including all terms in the viscous stress tensor. For the vertically separated components a shear layer forms between them. The middle of this layer free-falls to the disk center. The accretion rates are increased by factors $\sim 10^2-10^4$ over that of a conventional disc rotating in one direction with the same viscosity. The vertical width of the shear layer and the accretion rate are strongly dependent on the viscosity and the mass fraction of the counter-rotating gas. In the case of radially separated components where the inner disc co-rotates and the outer disc rotates in the opposite direction, a gap between the two components opens and closes quasi-periodically. The accretion rates are $\gtrsim 25$ times larger than those for a disc rotating in one direction with the same viscosity. △ Less

Submitted 29 October, 2014; v1 submitted 24 August, 2014; originally announced August 2014.

Comments: 11 pages, 13 figures, 1 table

Facing the wind of the pre-FUor V1331 Cyg

Authors: P. P. Petrov, R. Kurosawa, M. M. Romanova, J. F. Gameiro, M. Fernandez, E. V. Babina, S. A. Artemenko

Abstract: 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 s… ▽ More 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. △ Less

Submitted 19 June, 2014; v1 submitted 6 June, 2014; originally announced June 2014.

Comments: 11 pages, 12 figures; accepted for publication in MNRAS. Typographical errors have been corrected after the proof stage

Stellar Wind Erosion of Protoplanetary Discs

Authors: Neesha R. Schnepf, Richard V. E. Lovelace, Marina Romanova, Vladimir Airapetian

Abstract: An analytic model is developed for the erosion of protoplanetary gas discs by high velocity magnetized stellar winds. The winds are centrifugally driven from the surface of rapidly rotating, strongly magnetized young stars. The presence of the magnetic field in the wind leads to Reynolds numbers sufficiently large to cause a strongly turbulent wind/disk boundary layer which entrains and carries aw… ▽ More An analytic model is developed for the erosion of protoplanetary gas discs by high velocity magnetized stellar winds. The winds are centrifugally driven from the surface of rapidly rotating, strongly magnetized young stars. The presence of the magnetic field in the wind leads to Reynolds numbers sufficiently large to cause a strongly turbulent wind/disk boundary layer which entrains and carries away the disc gas. The model uses the conservation of mass and momentum in the turbulent boundary layer. The time-scale for significant erosion depends on the disc accretion speed and accretion rate and on the wind mass loss rate and the wind velocity. The time-scale is estimated to be 2E6 yr. A significant contribution to the disc erosion can come from frequent powerful coronal mass ejections (CMEs) where the average mass loss rate in CMEs and velocities have values comparable to those for the steady wind. △ Less

Submitted 9 January, 2015; v1 submitted 22 May, 2014; originally announced May 2014.

Comments: 7 pages, 3 figures

Simulations of Accretion onto Magnetized Stars: Results of 3D MHD Simulations and 3D Radiative Transfer

Authors: Marina Romanova, Ryuichi Kurosawa

Abstract: 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 O… ▽ More 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. △ Less

Submitted 19 January, 2014; originally announced January 2014.

Comments: 6 pages, 5 figures, to appear in the proceedings of ASTRONUM-2013 meeting: "8th International Conference on Numerical Modeling of Space Plasma Flows", Biarritz, France, July 1-5, 2013

Episodic Accretion in Young Stars

Authors: Marc Audard, Péter Ábrahám, Michael M. Dunham, Joel D. Green, Nicolas Grosso, Kenji Hamaguchi, Joel H. Kastner, Ágnes Kóspál, Giuseppe Lodato, Marina Romanova, Stephen L. Skinner, Eduard I. Vorobyov, Zhaohuan Zhu

Abstract: In the last twenty years, the topic of episodic accretion has gained significant interest in the star formation community. It is now viewed as a common, though still poorly understood, phenomenon in low-mass star formation. The FU Orionis objects (FUors) are long-studied examples of this phenomenon. FUors are believed to undergo accretion outbursts during which the accretion rate rapidly increases… ▽ More In the last twenty years, the topic of episodic accretion has gained significant interest in the star formation community. It is now viewed as a common, though still poorly understood, phenomenon in low-mass star formation. The FU Orionis objects (FUors) are long-studied examples of this phenomenon. FUors are believed to undergo accretion outbursts during which the accretion rate rapidly increases from typically $10^{-7}$ to a few $10^{-4}$ $M_\odot$ yr$^{-1}$, and remains elevated over several decades or more. EXors, a loosely defined class of pre-main sequence stars, exhibit shorter and repetitive outbursts, associated with lower accretion rates. The relationship between the two classes, and their connection to the standard pre-main sequence evolutionary sequence, is an open question: do they represent two distinct classes, are they triggered by the same physical mechanism, and do they occur in the same evolutionary phases? Over the past couple of decades, many theoretical and numerical models have been developed to explain the origin of FUor and EXor outbursts. In parallel, such accretion bursts have been detected at an increasing rate, and as observing techniques improve each individual outburst is studied in increasing detail. We summarize key observations of pre-main sequence star outbursts, and review the latest thinking on outburst triggering mechanisms, the propagation of outbursts from star/disk to disk/jet systems, the relation between classical EXors and FUors, and newly discovered outbursting sources -- all of which shed new light on episodic accretion. We finally highlight some of the most promising directions for this field in the near- and long-term. △ Less

Submitted 14 January, 2014; originally announced January 2014.

Comments: 24 pages, 8 figures. Accepted for publication as a review chapter in Protostars and Planets VI, University of Arizona Press (2014), eds. H. Beuther, R. Klessen, C. Dullemond, Th. Henning

Rossby Wave Instability in Astrophysical Discs

Authors: R. V. E. Lovelace, M. M. Romanova

Abstract: A brief review is given of the Rossby wave instability (RWI) in astrophysical discs. In non-self-gravitating discs, around for example a newly forming stars, the instability can be triggered by an axisymmetric bump at some radius $r_0$ in the disc surface mass-density. It gives rise to exponentially growing non-axisymmetric perturbation [$\propto \exp({\rm i}mφ)$, $m=1,2..$] in the vicinity of… ▽ More A brief review is given of the Rossby wave instability (RWI) in astrophysical discs. In non-self-gravitating discs, around for example a newly forming stars, the instability can be triggered by an axisymmetric bump at some radius $r_0$ in the disc surface mass-density. It gives rise to exponentially growing non-axisymmetric perturbation [$\propto \exp({\rm i}mφ)$, $m=1,2..$] in the vicinity of $r_0$ consisting of {\it anticyclonic} vortices. These vortices are regions of high pressure and consequently act to trap dust particles which in turn can facilitate planetesimal growth in proto-planetary discs. The Rossby vortices in the discs around stars and black holes may cause the observed quasi-periodic modulations of the disc's thermal emission. △ Less

Submitted 10 April, 2014; v1 submitted 16 December, 2013; originally announced December 2013.

Comments: 11 pages, 4 figures

Journal ref: Fluid Dynamics Research 46 (2014) 041401 (IOP Publishing: The Japan Society of Fluid Mechanics)

MHD Simulations of Magnetospheric Accretion, Ejection and Plasma-field Interaction

Authors: M. M. Romanova, R. V. E. Lovelace, M. Bachetti, A. A. Blinova, A. V. Koldoba, R. Kurosawa, P. S. Lii, G. V. Ustyugova

Abstract: 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. 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. △ Less

Submitted 18 November, 2013; originally announced November 2013.

Comments: 11 pages, 8 figures, conference proceedings: "Physics at the Magnetospheric Boundary", Geneva, Switzerland, 25-28 June, 2013

The effect of an ordered azimuthal magnetic field on a migrating planet in a non-turbulent disc

Authors: M. L. Comins, M. M. Romanova, A. V. Koldoba, G. V. Ustyugova, R. V. E. Lovelace

Abstract: In this work, we consider the physics of the interaction between a planet and a magnetized gaseous protoplanetary disc. We investigate the migration of a planet in a disc that is threaded with an azimuthal magnetic field. We find that, for a larger magnetic field amplitude, there is an increasingly large positive torque on the planet from the disc, resulting in slowed and even outward migration. O… ▽ More In this work, we consider the physics of the interaction between a planet and a magnetized gaseous protoplanetary disc. We investigate the migration of a planet in a disc that is threaded with an azimuthal magnetic field. We find that, for a larger magnetic field amplitude, there is an increasingly large positive torque on the planet from the disc, resulting in slowed and even outward migration. Our results indicate that magnetic resonances due to a purely azimuthal, ordered magnetic field can slow or stop the inward migration of Jupiter-mass, Saturn-mass, and $5 M_{\oplus}$ planets. △ Less

Submitted 24 October, 2013; originally announced October 2013.

Comments: 14 pages, 11 figures, submitted to MNRAS

Magnetospheric Accretions and the Inner Winds of Classical T Tauri Stars

Authors: Ryuichi Kurosawa, M. M. Romanova

Abstract: 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 a… ▽ More 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. △ Less

Submitted 5 October, 2013; originally announced October 2013.

Comments: 10 pages, 4 figures, to appear in the proceedings of IAUS 302: "Magnetic Fields Throughout Stellar Evolution", Biarritz, France, Aug. 26-30, 2013

Oscillations of the Boundary Layer and High-frequency QPOs

Authors: A. A. Blinova, M. Bachetti, M. M. Romanova

Abstract: We observed persistent high-frequency oscillations of the boundary layer near an accreting, weakly-magnetized star in global 3D MHD simulations. The tilted dipole magnetic field is not strong enough to open a gap between the star and the disk. Instead, it forms a highly-wrapped azimuthal field near the surface of the star which slows down rotation of the disk matter, while a small tilt of the fiel… ▽ More We observed persistent high-frequency oscillations of the boundary layer near an accreting, weakly-magnetized star in global 3D MHD simulations. The tilted dipole magnetic field is not strong enough to open a gap between the star and the disk. Instead, it forms a highly-wrapped azimuthal field near the surface of the star which slows down rotation of the disk matter, while a small tilt of the field excites oscillations of the boundary layer with a frequency below the Keplerian frequency. This mechanism may be responsible for the high-frequency oscillations in accreting neutron stars, white dwarfs and classical T Tauri stars. △ Less

Submitted 17 September, 2013; originally announced September 2013.

Comments: 4 pages, 3 figures, conference proceedings: "Physics at the Magnetospheric Boundary", Geneva, Switzerland, 25-28 June, 2013

Boundary Between Stable and Unstable Regimes of Accretion

Authors: A. A. Blinova, R. V. E. Lovelace, M. M. Romanova

Abstract: We investigated the boundary between stable and unstable regimes of accretion and its dependence on different parameters. Simulations were performed using a "cubed sphere" code with high grid resolution (244 grid points in the azimuthal direction), which is twice as high as that used in our earlier studies. We chose a very low viscosity value, with alpha-parameter alpha=0.02. We observed from the… ▽ More We investigated the boundary between stable and unstable regimes of accretion and its dependence on different parameters. Simulations were performed using a "cubed sphere" code with high grid resolution (244 grid points in the azimuthal direction), which is twice as high as that used in our earlier studies. We chose a very low viscosity value, with alpha-parameter alpha=0.02. We observed from the simulations that the boundary strongly depends on the ratio between magnetospheric radius r_m (where the magnetic stress in the magnetosphere matches the matter stress in the disk) and corotation radius r_cor (where the Keplerian velocity in the disk is equal to the angular velocity of the star). For a small misalignment angle of the dipole field, Theta=5 degrees, accretion is unstable if r_cor/r_m>1.35, and is stable otherwise. In cases of a larger misalignment angle of the dipole, Theta=20 degrees, instability occurs at slightly larger values, r_cor/r_m>1.41. △ Less

Submitted 17 September, 2013; originally announced September 2013.

Comments: 4 pages, 4 figures, conference proceedings: "Physics at the Magnetospheric Boundary", Geneva, Switzerland, 25-28 June, 2013

Magnetospheric accretion on the fully-convective classical T Tauri star DN Tau

Authors: JF Donati, SG Gregory, SHP Alencar, G Hussain, J Bouvier, MM Jardine, F Menard, C Dougados, MM Romanova, the MaPP collaboration

Abstract: We report here results of spectropolarimetric observations of the classical T Tauri star DN Tau carried out (at 2 epochs) with ESPaDOnS at the Canada-France-Hawaii Telescope within the `Magnetic Protostars and Planets' programme. We infer that DN Tau, with a photospheric temperature of 3,950+-50 K, a luminosity of 0.8+-0.2 Lsun and a rotation period of 6.32 d, is a ~2Myr-old fully-convective 0.65+… ▽ More We report here results of spectropolarimetric observations of the classical T Tauri star DN Tau carried out (at 2 epochs) with ESPaDOnS at the Canada-France-Hawaii Telescope within the `Magnetic Protostars and Planets' programme. We infer that DN Tau, with a photospheric temperature of 3,950+-50 K, a luminosity of 0.8+-0.2 Lsun and a rotation period of 6.32 d, is a ~2Myr-old fully-convective 0.65+-0.05 Msun star with a radius of 1.9+-0.2 Dsun, viewed at an inclination of 35+-10degr. Clear circularly-polarized Zeeman signatures are detected in both photospheric and accretion-powered emission lines, probing longitudinal fields of up to 1.8 kG (in the He1 D3 accretion proxy). Rotational modulation of Zeeman signatures, detected both in photospheric and accretion lines, is different between our 2 runs, providing further evidence that fields of cTTSs are generated by non-stationary dynamos. Using tomographic imaging, we reconstruct maps of the large-scale field, of the photospheric brightness and of the accretion-powered emission at the surface of DN Tau at both epochs. We find that the magnetic topology is mostly poloidal, and largely axisymmetric, with an octupolar component (of polar strength 0.6-0.8 kG) 1.5-2.0x larger than the dipolar component (of polar strength 0.3-0.5 kG). DN Tau features dominantly poleward accretion at both epochs. The large-scale dipole component of DN Tau is however too weak to disrupt the surrounding accretion disc further than 65-90% of the corotation radius (at which the disc Keplerian period matches the stellar rotation period), suggesting that DN Tau is already spinning up despite being fully convective. △ Less

Submitted 23 August, 2013; originally announced August 2013.

Comments: MNRAS in press (18 pages, 11 figures, 1 table)

Observable Signatures of Classical T Tauri Stars Accreting in an Unstable Regime

Authors: Ryuichi Kurosawa, M. M. Romanova

Abstract: 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… ▽ More 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. △ Less

Submitted 13 July, 2013; originally announced July 2013.

Comments: 6 pages, 3 figures, to appear in the proceedings of the conference: "Physics at the Magnetspheric Boundary", Geneva, Switzerland, 25-28 June, 2013

Magnetically Driven and Collimated Jets from the Disc-Magnetosphere Boundary of Rotating Stars

Authors: R. V. E. Lovelace, M. M. Romanova, P. Lii

Abstract: We discuss recent progress in understanding the launching of outflows/jets from the disc-magnetosphere boundary of slowly and rapidly rotating magnetized stars. In most of the discussed models the interior of the disc is assumed to have a turbulent viscosity and magnetic diffusivity (as described by two ``alpha'' parameters), whereas the coronal region outside of the disc is treated using ideal ma… ▽ More We discuss recent progress in understanding the launching of outflows/jets from the disc-magnetosphere boundary of slowly and rapidly rotating magnetized stars. In most of the discussed models the interior of the disc is assumed to have a turbulent viscosity and magnetic diffusivity (as described by two ``alpha'' parameters), whereas the coronal region outside of the disc is treated using ideal magnetohydrodynamics (MHD). Extensive MHD simulations have established the occurrence of long-lasting outflows in both the cases of slowly and rapidly rotating stars. In the case of {\it slowly rotating stars}, a new type of outflow, {\it a conical wind}, is found and studied in simulations. The conical winds appear in cases where the magnetic flux of the star is bunched up by the inward motion of the accretion disc. Near their region of origin, the winds have the shape of a thin conical shell with a half-opening angle of $ \sim 30^\circ$. At large distances these outflows become magnetically collimated by their toroidal magnetic field and form matter dominated jets. That is, the jets are current carrying. The predominant driving force for the conical winds is the magnetic force proportional to the negative gradient of the square of the toroidal magnetic field and not the centrifugal force. In the case of {\it rapidly rotating stars} in the ``propeller regime,'' a two-component outflow is observed. The first component is similar to the matter dominated conical winds. A large fraction of the disc matter may be ejected into the winds in this regime. The second component is a high-velocity, low-density magnetically dominated {\it axial jet} where matter flows along the opened polar field lines of the star. △ Less

Submitted 5 June, 2013; originally announced June 2013.

Comments: 14 pages, 8 figures. arXiv admin note: substantial text overlap with arXiv:0907.3394

Propeller outflows from an MRI disc

Authors: Patrick S Lii, Marina M. Romanova, Galina V. Ustyugova, Alexander V. Koldoba, Richard V. E. Lovelace

Abstract: We present the results of axisymmetric simulations of MRI-driven accretion onto a rapidly rotating, magnetized star accreting in the propeller regime. The stellar magnetosphere corotates with the star, forming a centrifugal barrier at the disc-magnetosphere boundary which inhibits matter accretion onto the star. Instead, the disc matter accumulates at the disc-magnetosphere interface and slowly di… ▽ More We present the results of axisymmetric simulations of MRI-driven accretion onto a rapidly rotating, magnetized star accreting in the propeller regime. The stellar magnetosphere corotates with the star, forming a centrifugal barrier at the disc-magnetosphere boundary which inhibits matter accretion onto the star. Instead, the disc matter accumulates at the disc-magnetosphere interface and slowly diffuses into the inner magnetosphere where it picks up angular momentum and is quickly ejected from the system as an outflow. Due to the interaction of the matter with the magnetosphere, this wind is discontinuous and is launched as discrete plasmoids. If the ejection rate is lower than the disc accretion rate, the matter accumulates at the disc-magnetosphere boundary faster than it can be ejected. In this case, accretion onto the star proceeds through the episodic accretion instability in which episodes of matter accumulation are followed by simultaneous accretion and ejection. During the accretion phase of this instability in which matter flows onto the star in funnel streams, we observe a corresponding rise in the outflow rate. Both the accretion and ejection processes observed in our simulations are highly non-stationary. The stars undergo strong spin-down due to the coupling of the stellar field with the disc and corona and we measure the spin-down timescales of around 1 Myr for a typical CTTS in the propeller regime. △ Less

Submitted 9 April, 2013; originally announced April 2013.

Comments: 13 pages, 10 figures, submitted to MNRAS

Analytical Hot Spot Shapes and Magnetospheric Radius from 3D Simulations of Magnetospheric Accretion

Authors: Akshay K. Kulkarni, Marina M. Romanova

Abstract: We present an analytical formula for the position and shape of the spots on the surface of accreting magnetized stars in cases where a star has a dipole magnetic field tilted at a small misalignment angle Theta < 30 degrees about the rotational axis, and the magnetosphere is 2.5-5 times the radius of the star. We observed that the azimuthal position of the spots varies significantly when the posit… ▽ More We present an analytical formula for the position and shape of the spots on the surface of accreting magnetized stars in cases where a star has a dipole magnetic field tilted at a small misalignment angle Theta < 30 degrees about the rotational axis, and the magnetosphere is 2.5-5 times the radius of the star. We observed that the azimuthal position of the spots varies significantly when the position of the inner disc varies. In contrast, the polar position of the spots varies only slightly because of the compression of the magnetosphere. The azimuthal width of the spots strongly varies with Theta: spots have the shape of an arc at larger misalignment angles, and resemble a ring at very small misalignment angles. The polar width of the spots varies only slightly with changes in parameters. The motion of the spots in the azimuthal direction can provide phase-shifts in accreting millisecond pulsars, and the "drift" of the period in Classical T Tauri stars. The position and shape of the spots are determined by three parameters: misalignment angle Theta; normalized corotation radius, r_c/R_* and normalized magnetospheric radius, r_m/R_*, where R_* is the stellar radius. We also use our data to check the formula for the Alfvén radius, where the main dependencies are r_m \sim (μ^2/\dot M)^{2/7}, where μis the magnetic moment of the star, and \dot M is the accretion rate. We found that the dependence is more gradual, r_m \sim (μ^2/\dot M)^{1/5}, which can be explained by the compression of the magnetosphere by the disc matter and by the non-dipole shape of the magnetic field lines of the external magnetosphere. △ Less

Submitted 19 March, 2013; originally announced March 2013.

Comments: 13 pages, 12 figures, Submitted to MNRAS

Spectral variability of classical T Tauri stars accreting in an unstable regime

Authors: Ryuichi Kurosawa, Marina M. Romanova

Abstract: 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… ▽ More 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. △ Less

Submitted 25 February, 2013; v1 submitted 3 January, 2013; originally announced January 2013.

Comments: 20 pages, 11 figures, 1 table, accepted for publication in MNRAS

Advection of Matter and B-Fields in Alpha-Discs

Authors: S. Dyda, R. V. E. Lovelace, G. V. Ustyugova, P. S. Lii, M. M. Romanova, A. V. Koldoba

Abstract: We have carried out and analyzed a set of axisymmetric MHD simulations of the evolution of a turbulent/diffusive accretion disc around an initially unmagnetized star. The disc is initially threaded by a weak magnetic field where the magnetic pressure is significantly less than the kinetic pressure in the disc. The viscosity and magnetic diffusivity are modelled by two "alpha" parameters, while the… ▽ More We have carried out and analyzed a set of axisymmetric MHD simulations of the evolution of a turbulent/diffusive accretion disc around an initially unmagnetized star. The disc is initially threaded by a weak magnetic field where the magnetic pressure is significantly less than the kinetic pressure in the disc. The viscosity and magnetic diffusivity are modelled by two "alpha" parameters, while the coronal region above the disc is treated using ideal MHD. The initial magnetic field is taken to consist of three poloidal field loops threading the disc. The motivation for this study is to understand the advection of disc matter and magnetic field by the turbulent/diffusive disc. At early times the innermost field loop twists and its field lines become open. The twisting of the opened field lines leads to the formation of both an inner collimated, magnetically-dominated jet, and at larger distances from the axis a matter dominated uncollimated wind. For later times, the strength of the magnetic field decreases owing to field reconnection and annihilation in the disc. For the early times, we have derived from the simulations both the matter accretion speed in the disc and the accretion speed of the magnetic field. We show that the derived matter accretion speed agrees approximately with the predictions of a model where the accretion speed is the sum of two terms, one due to the disc's viscosity (which gives a radial outflow of angular momentum in the disc), and a second due to the twisted magnetic field at the disc's surface (which gives a vertical outflow of angular momentum). For early times we find that the magnetic contribution is roughly twice the viscous contribution for the case where the alpha parameters are both equal to 0.1. At later times the magnetic contribution to the matter speed becomes small compared to the viscous contribution. △ Less

Submitted 1 April, 2013; v1 submitted 3 December, 2012; originally announced December 2012.

Comments: 14 pages, 17 figures

Warps, bending and density waves excited by rotating magnetized stars: results of global 3D MHD simulations

Authors: M. M. Romanova, G. V. Ustyugova, A. V. Koldoba, R. V. E. Lovelace

Abstract: We report results of the first global three-dimensional magnetohydrodynamic simulations of the waves excited in an accretion disc by a rotating star with a dipole magnetic field misaligned from the star's rotation axis (which is aligned with the disc axis). The main results are the following: (1) If the magnetosphere of the star corotates approximately with the inner disc, then we observe a strong… ▽ More We report results of the first global three-dimensional magnetohydrodynamic simulations of the waves excited in an accretion disc by a rotating star with a dipole magnetic field misaligned from the star's rotation axis (which is aligned with the disc axis). The main results are the following: (1) If the magnetosphere of the star corotates approximately with the inner disc, then we observe a strong one-armed bending wave (a warp). This warp corotates with the star and has a maximum amplitude between corotation radius and the radius of the vertical resonance. The disc's center of mass can deviate from the equatorial plane up to the distance of z_w\approx 0.1 r. However, the effective height of the warp can be larger, h_w \approx 0.3 r due to the finite thickness of the disc. Stars with a range of misalignment angles excite warps. However, the amplitude of the warps is larger for misalignment angles between 15 and 60 degrees. (2) If the magnetosphere rotates slower, than the inner disc, then a bending wave is excited at the disc-magnetosphere boundary, but does not form a large-scale warp. Instead, high-frequency oscillations become strong at the inner region of the disc. These are (a) trapped density waves which form inside the radius where the disc angular velocity has a maximum, and (b) inner bending waves which appear in the case of accretion through magnetic Raleigh-Taylor instability. These two types of waves are connected with the inner disc and their frequencies will vary with accretion rate. Bending oscillations at lower frequencies are also excited including global oscillations of the disc. In cases where the simulation region is small, slowly-precessing warp forms. Simulations are applicable to young stars, cataclysmic variables, and accreting millisecond pulsars. △ Less

Submitted 19 March, 2013; v1 submitted 5 September, 2012; originally announced September 2012.

Comments: 26 pages, 25 figures

Journal ref: MNRAS, 2013, vol. 430, p. 699-724

Magnetometry of the classical T Tauri star GQ Lup: non-stationary dynamos & spin evolution of young Suns

Authors: J. F. Donati, S. G. Gregory, S. H. P. Alencar, G. Hussain, J. Bouvier, C. Dougados, M. M. Jardine, F. Menard, M. M. Romanova, the MaPP Collaboration

Abstract: We report here results of spectropolarimetric observations of the classical T Tauri star (cTTS) GQ Lup carried out with ESPaDOnS at the Canada-France-Hawaii Telescope (CFHT) in the framework of the "Magnetic Protostars and Planets" (MaPP) programme, and obtained at 2 different epochs (2009 July & 2011 June). From these observations, we first infer that GQ Lup has a photospheric temperature of 4,30… ▽ More We report here results of spectropolarimetric observations of the classical T Tauri star (cTTS) GQ Lup carried out with ESPaDOnS at the Canada-France-Hawaii Telescope (CFHT) in the framework of the "Magnetic Protostars and Planets" (MaPP) programme, and obtained at 2 different epochs (2009 July & 2011 June). From these observations, we first infer that GQ Lup has a photospheric temperature of 4,300+-50Â K and a rotation period of 8.4+-0.3 d; it implies that it is a 1.05+-0.07 Msun star viewed at an inclination of ~30deg, with an age of 2-5 Myr, a radius of 1.7+-0.2 Rsun, and has just started to develop a radiative core. Large Zeeman signatures are clearly detected at all times, both in photospheric lines & in accretion-powered emission lines, probing longitudinal fields of up to 6 kG and hence making GQ Lup the cTTS with the strongest large-scale fields known as of today. Rotational modulation of Zeeman signatures is clearly different between our 2 runs, demonstrating that large-scale fields of cTTSs are evolving with time and are likely produced by non-stationary dynamo processes. Using tomographic imaging, we reconstruct maps of the large-scale field, of the photospheric brightness & of the accretion-powered emission of GQ Lup. We find that the magnetic topology is mostly poloidal & axisymmetric; moreover, the octupolar component of the large-scale field (of strength 2.4 & 1.6 kG in 2009 & 2011) dominates the dipolar component (of strength ~1 kG) by a factor of ~2, consistent with the fact that GQ Lup is no longer fully-convective. GQ Lup also features dominantly poleward magnetospheric accretion at both epochs. The large-scale dipole of GQ Lup is however not strong enough to disrupt the surrounding accretion disc further than about half-way to the corotation radius, suggesting that GQ Lup should rapidly spin up like other similar partly-convective cTTSs (abridged). △ Less

Submitted 8 June, 2012; originally announced June 2012.

Comments: MNRAS, in press (17 pages, 10 figures, 1 table)

Accretion dynamics in the classical T Tauri star V2129 Oph

Authors: S. H. P. Alencar, J. Bouvier, F. M. Walter, C. Dougados, J. -F. Donati, R. Kurosawa, M. Romanova, X. Bonfils, G. H. R. A. Lima, S. Massaro, M. Ibrahimov, E. Poretti

Abstract: 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 (58… ▽ More 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 Å) 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αand Hβ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. △ Less

Submitted 28 March, 2012; originally announced March 2012.

Comments: Accepted by Astronomy and Astrophysics

Line formation in the inner winds of classical T Tauri stars: testing the conical-shell wind solution

Authors: Ryuichi Kurosawa, Marina M. Romanova

Abstract: 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 essent… ▽ More 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. △ Less

Submitted 2 August, 2012; v1 submitted 22 March, 2012; originally announced March 2012.

Comments: 16 pages, 10 figures, accepted for publication in MNRAS

MRI-driven Accretion on to Magnetized stars: Global 3D MHD Simulations of Magnetospheric and Boundary Layer Regimes

Authors: M. M. Romanova, G. V. Ustyugova, A. V. Koldoba, R. V. E. Lovelace

Abstract: We discuss results of global 3D MHD simulations of accretion on to a rotating magnetized star with a tilted dipole magnetic field, where the accretion is driven by the magneto-rotational instability (MRI). The simulations show that MRI-driven turbulence develops in the disc, and angular momentum is transported outwards due primarily to the magnetic stress. The turbulent flow is strongly inhomogene… ▽ More We discuss results of global 3D MHD simulations of accretion on to a rotating magnetized star with a tilted dipole magnetic field, where the accretion is driven by the magneto-rotational instability (MRI). The simulations show that MRI-driven turbulence develops in the disc, and angular momentum is transported outwards due primarily to the magnetic stress. The turbulent flow is strongly inhomogeneous and the densest matter is in azimuthally-stretched turbulent cells. We investigate two regimes of accretion: a magnetospheric regime and a boundary layer (BL) regime. In the magnetospheric regime, the accretion disc is truncated by the star's magnetic field within a few stellar radii from the star, and matter flows to the star in funnel streams. The funnel streams flowing towards the south and north magnetic poles but are not equal due to the inhomogeneity of the flow. In the BL regime, matter accretes to the surface of the star through the boundary layer. The magnetic field in the inner disc is strongly amplified by the shear of the accretion flow, and the matter and magnetic stresses become comparable. Accreting matter forms a belt-shaped region on the surface of the star. The belt has inhomogeneous density distribution which varies in time due to variable accretion rate. Results of simulations can be applied to classical T Tauri stars, accreting brown dwarfs, millisecond pulsars, dwarf novae cataclysmic variables, and other stars with magnetospheres smaller than several stellar radii. △ Less

Submitted 13 November, 2011; originally announced November 2011.

Comments: 15 pages, 13 figures, accepted by MNRAS

Bondi-Hoyle Accretion onto Magnetized Neutron Star

Authors: O. D. Toropina, M. M. Romanova, R. V. E. Lovelace

Abstract: Axisymmetric MHD simulations are used to investigate the Bondi-Hoyle accretion onto an isolated magnetized neutron star moving supersonically (with Mach number of 3) through the interstellar medium. The star is assumed to have a dipole magnetic field aligned with its motion and a magnetospheric radius R_m less then the accretion radius R_BH, so that the gravitational focusing is important. We find… ▽ More Axisymmetric MHD simulations are used to investigate the Bondi-Hoyle accretion onto an isolated magnetized neutron star moving supersonically (with Mach number of 3) through the interstellar medium. The star is assumed to have a dipole magnetic field aligned with its motion and a magnetospheric radius R_m less then the accretion radius R_BH, so that the gravitational focusing is important. We find that the accretion rate to a magnetized star is smaller than that to a non-magnetized star for the parameters considered. Close to the star the accreting matter falls to the star's surface along the magnetic poles with a larger mass flow to the leeward pole of the star. In the case of a relatively large stellar magnetic field, the star's magnetic field is stretched in the direction of the matter flow outside of R_m (towards the windward side of the star). For weaker magnetic fields we observed oscillations of the closed magnetosphere frontward and backward. These are accompanied by strong oscillations of the mass accretion rate which varies by factors ~ 3. Old slowly rotating neutron stars with no radio emission may be visible in the X-ray band due to accretion of interstellar matter. In general, the star's velocity, magnetic moment, and angular velocity vectors may all be in different directions so that the accretion luminosity will be modulated at the star's rotation rate. △ Less

Submitted 10 November, 2011; originally announced November 2011.

Comments: 9 pages, 9 figures, accepted for publication in MNRAS

The large-scale magnetic field and poleward mass accretion of the classical T Tauri star TW Hya

Authors: JF Donati, SG Gregory, SHP Alencar, J Bouvier, G Hussain, M Skelly, C Dougados, MM Jardine, F Menard, MM Romanova, YC Unruh, the MaPP collaboration

Abstract: We report here results of spectropolarimetric observations of the ~8Myr classical TTauri star (cTTS) TWHya carried out with ESPaDOnS at the Canada-France-Hawaii Telescope (CFHT) in the framework of the `Magnetic Protostars and Planets' (MaPP) programme, and obtained at 2 different epochs (2008 March and 2010 March). Obvious Zeeman signatures are detected at all times, both in photospheric lines an… ▽ More We report here results of spectropolarimetric observations of the ~8Myr classical TTauri star (cTTS) TWHya carried out with ESPaDOnS at the Canada-France-Hawaii Telescope (CFHT) in the framework of the `Magnetic Protostars and Planets' (MaPP) programme, and obtained at 2 different epochs (2008 March and 2010 March). Obvious Zeeman signatures are detected at all times, both in photospheric lines and in accretion-powered emission lines. Significant intrinsic variability and moderate rotational modulation is observed in both photospheric and accretion proxies. Using tomographic imaging, we reconstruct maps of the large-scale field, of the photospheric brightness and of the accretion-powered emission at the surface of TWHya at both epochs. We find that the magnetic topology is mostly poloidal and axisymmetric with respect to the rotation axis of the star, and that the octupolar component of the large-scale field (2.5-2.8kG at the pole) largely dominates the dipolar component. This large-scale field topology is characteristic of partly-convective stars, supporting the conclusion (from evolutionary models) that TWHya already hosts a radiative core. We also show that TWHya features a high-latitude photospheric cool spot overlapping with the main magnetic pole (and producing the observed radial velocity fluctuations); this is also where accretion concentrates most of the time, although accretion at lower latitudes is found to occur episodically. We propose that the relatively rapid rotation of TWHya (with respect to AATau-like cTTSs) directly reflects the weakness of the large-scale dipole, no longer capable of magnetically disrupting the accretion disc up to the corotation radius (at which the Keplerian period equals the stellar rotation period). We therefore conclude that TWHya is in a phase of rapid spin-up as its large-scale dipole field progressively vanishes. △ Less

Submitted 21 June, 2011; originally announced June 2011.

Comments: MNRAS in press

Are the outflows in FU Orionis systems driven by the stellar magnetic field?

Authors: Arieh Königl, Marina M. Romanova, Richard V. E. Lovelace

Abstract: FU Orionis (FUOR) outbursts are major optical brightening episodes in low-mass protostars that correspond to rapid mass-accretion events in the innermost region of a protostellar disc. The outbursts are accompanied by strong outflows, with the inferred mass outflow rates reaching ~10% of the mass inflow rates. Shu et al. proposed that the outflows represent accreted disc material that is driven ce… ▽ More FU Orionis (FUOR) outbursts are major optical brightening episodes in low-mass protostars that correspond to rapid mass-accretion events in the innermost region of a protostellar disc. The outbursts are accompanied by strong outflows, with the inferred mass outflow rates reaching ~10% of the mass inflow rates. Shu et al. proposed that the outflows represent accreted disc material that is driven centrifugally from the spun-up surface layers of the protostar by the stellar magnetic field. This model was critiqued by Calvet et al., who argued that it cannot reproduce the photospheric absorption-line shifts observed in the prototype object FU Ori. Calvet et al. proposed that the wind is launched, instead, from the surface of the disc on scales of a few stellar radii by a non-stellar magnetic field. In this paper we present results from numerical simulations of disc accretion on to a slowly rotating star with an aligned magnetic dipole moment that gives rise to a kilogauss-strength surface field. We demonstrate that, for parameters appropriate to FU Ori, such a system can develop a strong, collimated disc outflow of the type previously identified by Romanova et al. in simulations of protostars with low and moderate accretion rates. At the high accretion rate that characterizes FUOR outbursts, the radius at which the disc is truncated by the stellar magnetic field moves much closer to the stellar surface, but the basic properties of the outflow, which is launched from the vicinity of the truncation radius along opened-up stellar magnetic field lines, remain the same, and are distinct from those of the mechanism proposed by Shu et al. We show that the simulated outflow can in principle account for the main observed characteristics of FUOR winds, including the photospheric line shifts measured in FU Ori. △ Less

Submitted 12 June, 2011; originally announced June 2011.

Comments: 11 pages, 4 figures, 1 table, to appear in MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc. 416 (2011) 757-766

Magnetic launching and collimation of jets from the disk-magnetosphere boundary: 2.5D MHD simulations

Authors: Patrick S Lii, Marina M. Romanova, Richard V. Lovelace

Abstract: We use axisymmetric magnetohydrodynamic (MHD) simulations to investigate the launching and collimation of jets emerging from the disk-magnetosphere boundary of accreting magnetized stars. Our analysis shows that the matter flows into the jet from the inner edge of the accretion disk. It is magnetically accelerated along field lines extending up from the disk and simultaneously collimated by the ma… ▽ More We use axisymmetric magnetohydrodynamic (MHD) simulations to investigate the launching and collimation of jets emerging from the disk-magnetosphere boundary of accreting magnetized stars. Our analysis shows that the matter flows into the jet from the inner edge of the accretion disk. It is magnetically accelerated along field lines extending up from the disk and simultaneously collimated by the magnetic pinch force. In the reference run which we use for analysis, the matter in the jet crosses the Alfvén surface a few R_* above the disk and the fast magnetosonic surface ~13 R_* above the disk. At larger distances, the magnetic pressure is a few times smaller than the total matter pressure but the magnetic force continues to accelerate and collimate the jet. In steady state, we observe a matter ejection-to-accretion ratio of ~0.2. Across different simulation runs, we measure a range of half-opening angles between Θ \approx 4° and 20° at the top of the simulation region, depending on the degree of magnetization in the outflow. We consider the case of stars undergoing epochs of high accretion (such as EXors, FUORs, and CTTSs) where the stellar magnetosphere is strongly compressed by the incoming accretion disk. For a typical EXor (mass 0.8 M_sun, radius 2 R_sun) accreting at ~10^-5 M_sun/yr, we measure poloidal velocities in the jet ranging from 30 km/s on the outer edge of the jet to more than 260 km/s on the inner edge. In general, the models can be applied to a variety of magnetized stars--white dwarfs, neutron stars, and brown dwarfs--which exhibit periods of high accretion. △ Less

Submitted 3 November, 2011; v1 submitted 21 April, 2011; originally announced April 2011.

Comments: 15 pages, 13 figures, accepted for publication by MNRAS, supplementary animation available at http://astro.cornell.edu/~pslii/research.htm

MRI-driven Accretion onto Magnetized stars: Axisymmetric MHD Simulations

Authors: Marina M. Romanova, Galina V. Ustyugova, Alexander V. Koldoba, Richard V. E. Lovelace

Abstract: We present the first results of a global axisymmetric simulation of accretion onto rotating magnetized stars from a turbulent, MRI-driven disk. The angular momentum is transported outward by the magnetic stress of the turbulent flow with a rate corresponding to a Shakura-Sunyaev viscosity parameter alpha\approx 0.01-0.04. The result of the disk-magnetosphere interaction depends on the orientation… ▽ More We present the first results of a global axisymmetric simulation of accretion onto rotating magnetized stars from a turbulent, MRI-driven disk. The angular momentum is transported outward by the magnetic stress of the turbulent flow with a rate corresponding to a Shakura-Sunyaev viscosity parameter alpha\approx 0.01-0.04. The result of the disk-magnetosphere interaction depends on the orientation of the poloidal field in the disk relative to that of the star at the disk-magnetosphere boundary. If fields have the same polarity, then the magnetic flux is accumulated at the boundary and blocks the accretion which leads to the accumulation of matter at the boundary. Subsequently, this matter accretes to the star in outburst before accumulating again. Hence, the cycling, `bursty' accretion is observed. If the disc and stellar fields have opposite polarity, then the field reconnection enhances the penetration of the disk matter towards the deeper field lines of the magnetosphere. However, the magnetic stress at the boundary is lower due to the field reconnection. This decreases the accretion rate and leads to smoother accretion at a lower rate. Test simulations show that in the case of higher accretion rate corresponding to alpha=0.05-0.1, accretion is bursty in cases of both polarities. On the other hand, at much lower accretion rates corresponding to alpha < 0.01, accretion is not bursty in any of these cases. We conclude that the episodic, bursty accretion is expected during periods of higher accretion rates in the disc, and in some cases it may alternate between bursty and smooth accretion, if the disk brings the poloidal field of alternating polarity. We find that a rotating, magnetically-dominated corona forms above and below the disk, and that it slowly expands outward, driven by the magnetic force. △ Less

Submitted 5 February, 2011; originally announced February 2011.

Comments: 20 pages, 20 figures, 2 tables, submitted to MNRAS

Multidimensional models of hydrogen and helium emission line profiles for classical T Tauri Stars: method, tests and examples

Authors: Ryuichi Kurosawa, Marina M. Romanova, Tim J. Harries

Abstract: 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 an… ▽ More 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. △ Less

Submitted 9 June, 2011; v1 submitted 3 February, 2011; originally announced February 2011.

Comments: 18 pages, 12 figures, accepted for publication in MNRAS

Non-stationary dynamo & magnetospheric accretion processes of the classical T Tauri star V2129 Oph

Authors: JF Donati, J Bouvier, FM Walter, SG Gregory, MB Skelly, GAJ Hussain, E Flaccomio, C Argiroffi, KN Grankin, MM Jardine, F Menard, C Dougados, MM Romanova, the MaPP collaboration

Abstract: We report here the first results of a multi-wavelength campaign focussing on magnetospheric accretion processes of the classical TTauri star (cTTS) V2129Oph. In this paper, we present spectropolarimetric observations collected in 2009 July with ESPaDOnS at the Canada-France-Hawaii Telescope (CFHT). Circularly polarised Zeeman signatures are clearly detected, both in photospheric absorption and acc… ▽ More We report here the first results of a multi-wavelength campaign focussing on magnetospheric accretion processes of the classical TTauri star (cTTS) V2129Oph. In this paper, we present spectropolarimetric observations collected in 2009 July with ESPaDOnS at the Canada-France-Hawaii Telescope (CFHT). Circularly polarised Zeeman signatures are clearly detected, both in photospheric absorption and accretion-powered emission lines, from time-series of which we reconstruct new maps of the magnetic field, photospheric brightness and accretion-powered emission at the surface of V2129Oph using our newest tomographic imaging tool - to be compared with those derived from our old 2005 June data set, reanalyzed in the exact same way. We find that in 2009 July, V2129Oph hosts octupolar & dipolar field components of about 2.1 & 0.9kG respectively, both tilted by about 20deg with respect to the rotation axis; we conclude that the large-scale magnetic topology changed significantly since 2005 June (when the octupole and dipole components were about 1.5 and 3 times weaker respectively), demonstrating that the field of V2129Oph is generated by a non-stationary dynamo. We also show that V2129Oph features a dark photospheric spot and a localised area of accretion-powered emission, both close to the main surface magnetic region (hosting fields of up to about 4kG in 2009 July). We finally obtain that the surface shear of V2129Oph is about half as strong as solar. From the fluxes of accretion-powered emission lines, we estimate that the observed average logarithmic accretion rate (in Msun/yr) at the surface of V2129Oph is -9.2+-0.3 at both epochs, peaking at -9.0 at magnetic maximum. It implies in particular that the radius at which the magnetic field of V2129Oph truncates the inner accretion disc is 0.93x and 0.50x the corotation radius in 2009 July and 2005 June respectively. △ Less

Submitted 23 November, 2010; v1 submitted 22 November, 2010; originally announced November 2010.

Comments: MNRAS in press - 16 pages, 9 figures

The Moving Hotspots model for kHz QPOs in accreting neutron stars

Authors: Matteo Bachetti, Marina Romanova

Abstract: 3D MHD simulation of accretion onto neutron stars have shown in the last few years that the footprint (hotspot) of the accretion flow changes with time. Two different kinds of accretion, namely the funnel flow and the equatorial accretion produced by instabilities at the inner disk, produce different kinds of motion of the hotspot. The funnel flow produces hotspots that move around the magnetic po… ▽ More 3D MHD simulation of accretion onto neutron stars have shown in the last few years that the footprint (hotspot) of the accretion flow changes with time. Two different kinds of accretion, namely the funnel flow and the equatorial accretion produced by instabilities at the inner disk, produce different kinds of motion of the hotspot. The funnel flow produces hotspots that move around the magnetic pole, while instabilities produce other hotspots that appear randomly and move along the equator or slightly above. The angular velocities of the two hotspots are different, the equatorial one being higher and both close to the Keplerian velocity in the inner region. Modeling of the lightcurves of these hotspots with Monte Carlo simulations show that the signatures produced in power specra by them, if observed, are QPOs plus low frequency components. Their frequencies, general behavior and features describe correctly most of the properties of kHz QPOs, if we assume the funnel flow hotspots as the origin of the lower kHz QPO and instabilities as the origin of the upper kHz QPO. △ Less

Submitted 31 October, 2010; originally announced November 2010.

Comments: 4 pages, 2 figures, to appear in the Proceedings of the ASTRONS 2010 conference

Global 3D Simulations of Disc Accretion onto the classical T Tauri Star BP Tauri

Authors: M. Long, M. M. Romanova, A. K. Kulkarni, J. -F. Donati

Abstract: The magnetic field of the classical T Tauri star BP Tau can be approximated as a superposition of dipole and octupole moments with respective strengths of the polar magnetic fields of 1.2 kG and 1.6 kG (Donati et al. 2008). We adopt the measured properties of BP Tau and model the disc accretion onto the star. We observed in simulations that the disc is disrupted by the dipole component and matter… ▽ More The magnetic field of the classical T Tauri star BP Tau can be approximated as a superposition of dipole and octupole moments with respective strengths of the polar magnetic fields of 1.2 kG and 1.6 kG (Donati et al. 2008). We adopt the measured properties of BP Tau and model the disc accretion onto the star. We observed in simulations that the disc is disrupted by the dipole component and matter flows towards the star in two funnel streams which form two accretion spots below the dipole magnetic poles. The octupolar component becomes dynamically important very close to the star and it redirects the matter flow to higher latitudes. The spots are meridionally elongated and are located at higher latitudes, compared with the pure dipole case, where crescent-shaped, latitudinally elongated spots form at lower latitudes. The position and shape of the spots are in good agreement with observations. The disk-magnetosphere interaction leads to the inflation of the field lines and to the formation of magnetic towers above and below the disk. The magnetic field of BP Tau is close to the potential only near the star, inside the magnetospheric surface, where magnetic stress dominates over the matter stress. A series of simulation runs were performed for different accretion rates. They show that an accretion rate is lower than obtained in many observations, unless the disc is truncated close to the star. The torque acting on the star is about an order of magnitude lower than that which is required for the rotational equilibrium. We suggest that a star could lose most of its angular momentum at earlier stages of its evolution. △ Less

Submitted 16 September, 2010; originally announced September 2010.

Comments: 11 pages, 13 figures, submitted to MNRAS

One-sided Outflows/Jets from Rotating Stars with Complex Magnetic Fields

Authors: R. V. E. Lovelace, M. M. Romanova, G. V. Ustyugova, A. V. Koldoba

Abstract: We investigate the generation of intrinsically asymmetric or {\it one-sided} outflows or jets from disk accretion onto rotating stars with complex magnetic fields using axisymmetric (2.5D) magnetohydrodynamic simulations. The intrinsic magnetic field of the star is assumed to consist of a superposition of an aligned dipole and an aligned quadrupole in different proportions. The star is assume… ▽ More We investigate the generation of intrinsically asymmetric or {\it one-sided} outflows or jets from disk accretion onto rotating stars with complex magnetic fields using axisymmetric (2.5D) magnetohydrodynamic simulations. The intrinsic magnetic field of the star is assumed to consist of a superposition of an aligned dipole and an aligned quadrupole in different proportions. The star is assumed to be rapidly rotating in the sense that the star's magnetosphere is in the propeller regime where strong outflows occur. Our simulations show that for conditions where there is a significant quadrupole component in addition to the dipole component, then a dominantly {\it one-sided} conical wind tends to form on the side of the equatorial plane with the larger value of the intrinsic axial magnetic field at a given distance. For cases where the quadrupole component is absent or very small, we find that dominantly one-sided outflows also form, but the direction of the flow "flip-flops" between upward and downward on a time-scale of $\sim 30$ days for a protostar. The average outflow will thus be symmetrical. In the case of a pure quadrupole field we find symmetric outflows in the upward and downward directions. △ Less

Submitted 2 April, 2010; originally announced April 2010.

Comments: 10 pages, 11 figures

Global 3D Simulations of Disc Accretion onto the classical T Tauri Star V2129 Oph

Authors: Marina M. Romanova, Min Long, Frederick K. Lamb, Akshay K. Kulkarni, Jean-Francois Donati

Abstract: The magnetic field of the classical T Tauri star V2129 Oph can be modeled approximately by superposing slightly tilted dipole and octupole moments, with polar magnetic field strengths of 0.35kG and 1.2kG respectively (Donati et al. 2007). Here we construct a numerical model of V2129 Oph incorporating this result and simulate accretion onto the star. Simulations show that the disk is truncated by t… ▽ More The magnetic field of the classical T Tauri star V2129 Oph can be modeled approximately by superposing slightly tilted dipole and octupole moments, with polar magnetic field strengths of 0.35kG and 1.2kG respectively (Donati et al. 2007). Here we construct a numerical model of V2129 Oph incorporating this result and simulate accretion onto the star. Simulations show that the disk is truncated by the dipole component and matter flows towards the star in two funnel streams. Closer to the star, the flow is redirected by the octupolar component, with some of the matter flowing towards the high-latitude poles, and the rest into the octupolar belts. The shape and position of the spots differ from those in a pure dipole case, where crescent-shaped spots are observed at the intermediate latitudes. Simulations show that if the disk is truncated at the distance of 6.2 R_* which is comparable with the co-rotation radius, 6.8 R_*, then the high-latitude polar spots dominate, but the accretion rate obtained from the simulations is about an order of magnitude lower than the observed one. The accretion rate matches the observed one if the disk is disrupted much closer to the star, at 3.4 R_*. However, the octupolar belt spots strongly dominate. Better match has been obtained in experiments with a dipole field twice as strong. The torque on the star from the disk-magnetosphere interaction is small, and the time-scale of spin evolution, 2 x10^7-10^9 years is longer than the 2x10^6 years age of V2129 Oph. The external magnetic flux of the star is strongly influenced by the disk: the field lines connecting the disk and the star inflate and form magnetic towers above and below the disk. The potential (vacuum) approximation is still valid inside the Alfvén (magnetospheric) surface where the magnetic stress dominates over the matter stress. △ Less

Submitted 16 September, 2010; v1 submitted 9 December, 2009; originally announced December 2009.

Comments: 15 pages, 15 figures, after major revision, added 3 figures, 2 tables. Accepted to MNRAS