Showing 1–50 of 119 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

Dynamical downfolding for localized quantum states

Authors: Mariya Romanova, Guorong Weng, Arsineh Apelian, Vojtech Vlcek

Abstract: We introduce an approach to treat localized correlated electronic states in the otherwise weakly correlated host medium. Here, the environment is dynamically downfolded on the correlated subspace. It is captured via renormalization of one and two quasiparticle interaction terms which are evaluated using many-body perturbation theory. We outline the strategy on how to take the dynamical effects int… ▽ More We introduce an approach to treat localized correlated electronic states in the otherwise weakly correlated host medium. Here, the environment is dynamically downfolded on the correlated subspace. It is captured via renormalization of one and two quasiparticle interaction terms which are evaluated using many-body perturbation theory. We outline the strategy on how to take the dynamical effects into account by going beyond the static limit approximation. Further, we introduce an efficient stochastic implementation that enables treating the host environment with a large number of electrons at a minimal computational cost. For small explicitly correlated subspace, the dynamical effects are critical. We demonstrate the methodology by reproducing optical excitations in the negatively charged NV center defect in diamond, that are in excellent agreement with experimental results. △ Less

Submitted 10 November, 2022; originally announced November 2022.

Reduced scaling of optimal regional orbital localization via sequential exhaustion of the single-particle space

Authors: Guorong Weng, Mariya Romanova, Arsineh Apelian, Hanbin Song, Vojtěch Vlček

Abstract: Wannier functions have become a powerful tool in the electronic structure calculations of extended systems. The generalized Pipek-Mezey Wannier functions exhibit appealing characteristics (e.g., reaching an optimal localization and the separation of the $σしぐま$-$πぱい$ orbitals) when compared with other schemes. However, when applied to giant nanoscale systems, the orbital localization suffers from a larg… ▽ More Wannier functions have become a powerful tool in the electronic structure calculations of extended systems. The generalized Pipek-Mezey Wannier functions exhibit appealing characteristics (e.g., reaching an optimal localization and the separation of the $σしぐま$-$πぱい$ orbitals) when compared with other schemes. However, when applied to giant nanoscale systems, the orbital localization suffers from a large computational cost overhead when one is interested in localized states in a small fragment of the system. Herein we present a swift, efficient, and robust approach for obtaining regionally localized orbitals of a subsystem within the generalized Pipek-Mezey scheme. The proposed algorithm introduces a reduced workspace and sequentially exhausts the entire orbital space until the convergence of the localization functional. It tackles systems with $\sim$10000 electrons within 0.5 hours with no loss in localization quality compared to the traditional approach. Regionally localized orbitals with a higher extent of localization are obtained via judiciously extending the subsystem's size. Exemplifying on large bulk and a 4-nm wide slab of diamond with NV$^-$ center, we demonstrate the methodology and discuss how the choice of the localization region affects the excitation energy of the defect. Furthermore, we show how the sequential algorithm is easily extended to stochastic methodologies that do not provide individual single-particle eigenstates. It is thus a promising tool to obtain regionally localized states for solving the electronic structure problems of a subsystem embedded in giant condensed systems. △ Less

Submitted 1 April, 2022; originally announced April 2022.

Comments: 36 pages, 7 figures, 2 tables

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)

Stochastic many-body calculations of moiré states in twisted bilayer graphene at high pressures

Authors: Mariya Romanova, Vojtech Vlcek

Abstract: Atomistic first principles many-body computational studies were so far limited by the system size. In this work, we apply and expand the stochastic GW method allowing calculations of quasiparticle energies of giant systems. We introduce three new developments to the stochastic many-body formalism: efficient evaluation of off-diagonal self-energy terms, construction of Dyson orbitals, and introduce… ▽ More Atomistic first principles many-body computational studies were so far limited by the system size. In this work, we apply and expand the stochastic GW method allowing calculations of quasiparticle energies of giant systems. We introduce three new developments to the stochastic many-body formalism: efficient evaluation of off-diagonal self-energy terms, construction of Dyson orbitals, and introduce stochastic cRPA, enabling efficient downfolding onto model Hamiltonians. We exemplify the new methodology by exploring twisted bilayer graphene (tBLG); its moire periodicity is associated with giant unit-cells hosting correlated electrons in flat bands. The resulting behavior of tBLG is governed by the coupling between the weakly and strongly correlated electrons in individual monolayers, controlled by twist or applied pressure. Here we study the latter scenario. The calculations document the formation of charge carrier localization at the Fermi level while most weakly correlated states are merely shifted in energy. This contrasts with the mean-field results. We show how to efficiently downfold the correlated subspace on a model Hamiltonian, namely the Hubbard model with a screened frequency dependent on-site interaction, computed with newly developed stochastic cRPA. For the $6^{\circ}$ tBLG system, the onsite interactions range between 200 and 300 meV under compression. The Dyson orbitals exhibit spatial distribution similar to the mean-field single particle states. Under pressure, the electron-electron interactions in the localized states increase; the dynamical screening does not fully compensate the dominant bare Coulomb interaction. The correlated moiré states thus appear only moderately screened at high pressures suggesting the presence of insulating states in compressed tBLG. △ Less

Submitted 3 November, 2021; v1 submitted 3 June, 2021; originally announced June 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

Role of precursor composition in the polymorph transformations, morphology control and ferromagnetic properties of nanosized TiO$_2$

Authors: Dmitry Zablotsky, Mikhail M. Maiorov, Aija Krumina, Marina Romanova, Elmars Blums

Abstract: Pure phase and mixed phase TiO$_2$ nanoparticles have been produced using a pyrolytic method from a non-aqueous carboxylate precursor. The precursor was prepared by a multiphase cation exchange using pentanoic acid (C$_4$H$_9$COOH). The thermal stability, polymorph content, morphology, size distribution and surface region of the produced nanoparticles were studied by TGA/DSC, XRD, FTIR and TEM. Hi… ▽ More Pure phase and mixed phase TiO$_2$ nanoparticles have been produced using a pyrolytic method from a non-aqueous carboxylate precursor. The precursor was prepared by a multiphase cation exchange using pentanoic acid (C$_4$H$_9$COOH). The thermal stability, polymorph content, morphology, size distribution and surface region of the produced nanoparticles were studied by TGA/DSC, XRD, FTIR and TEM. High quality monodisperse nanoparticles have been produced in the size range from 7 to 27 nm. The nanoparticles showed room temperature ferromagnetism revealed by VSM within bound polaron model. The carboxylate precursor is a good alternative to standard sol-gel to produce nanoparticles free from impurities. △ Less

Submitted 5 February, 2021; originally announced February 2021.

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

Laser scattering by submicron droplets originated during the electrical explosion of thin metal wires

Authors: V. M. Romanova, G. V. Ivanenkov, E. V. Parkevich, I. N. Tilikin, M. A. Medvedev, T. A. Shelkovenko, S. A. Pikuz, A. S. Selyukov

Abstract: This paper presents the results of studying of dispersed media formation during the electrical explosion of thin metal wires in vacuum by using low-current generators ($\sim 1$-$10$~kA). Particular attention is paid to the analysis of the composition and structure of the corresponding explosion products as well as to the problem of their visualization using simultaneous laser interferometry and sh… ▽ More This paper presents the results of studying of dispersed media formation during the electrical explosion of thin metal wires in vacuum by using low-current generators ($\sim 1$-$10$~kA). Particular attention is paid to the analysis of the composition and structure of the corresponding explosion products as well as to the problem of their visualization using simultaneous laser interferometry and shadow imaging at two wavelengths (1.064~$μみゅー$m and 0.532~$μみゅー$m). Our findings point to the fact that the important role in the visualization of the explosion products belongs to multiple scattering by submicron droplets of dense condensed matter, which are mixed with metal vapor. The hypothesis on the existence of submicron droplets in the products of exploding metal wires correlates with the results obtained by soft x-ray radiography combined with a laser probing technique. Taking into account the multiple scattering by submicron droplets, it is possible to significantly clarify the parameters of the explosion products visualized via laser probing techniques as well as to gain a deeper insight into the physics behind the electrical wire explosion. △ Less

Submitted 3 December, 2020; v1 submitted 27 October, 2020; originally announced November 2020.

Comments: 8 pages, 5 figures Changes: 1. 1st page: text was removed "(Dated: December 1, 2020)" - this (compilation) date do not relate to the submission date; 2. 7th page: the acknolegment section was changed. 3. 3rd page: width of the 1st figure was changed from 15.5 cm to 15.6 cm (corrected structure of preprint)

Are multi-quasiparticle interactions important in molecular ionization?

Authors: Carlos Mejuto-Zaera, Guorong Weng, Mariya Romanova, Stephen J. Cotton, K. Birgitta Whaley, Norm M. Tubman, Vojtěch Vlček

Abstract: Photo-emission spectroscopy directly probes individual electronic states, ranging from single excitations to high-energy satellites, which simultaneously represent multiple quasiparticles (QPs) and encode information about electronic correlation. First-principles description of the spectra requires an efficient and accurate treatment of all many-body effects. This is especially challenging for inn… ▽ More Photo-emission spectroscopy directly probes individual electronic states, ranging from single excitations to high-energy satellites, which simultaneously represent multiple quasiparticles (QPs) and encode information about electronic correlation. First-principles description of the spectra requires an efficient and accurate treatment of all many-body effects. This is especially challenging for inner valence excitations where the single QP picture breaks down. Here, we provide the full valence spectra of small closed-shell molecules, exploring the independent and interacting quasiparticle regimes, computed with the fully-correlated adaptive sampling configuration interaction (ASCI) method. We critically compare these results to calculations with the many-body perturbation theory, based on the $GW$ and vertex corrected $GWΓがんま$ approaches. The latter explicitly accounts for two-QP quantum interactions, which have been often neglected. We demonstrate that for molecular systems, the vertex correction universally improves the theoretical spectra, and it is crucial for accurate prediction of QPs as well as capturing the rich satellite structures of high-energy excitations. $GWΓがんま$ offers a unified description across all relevant energy scales. Our results suggest that the multi-QP regime corresponds to dynamical correlations, which can be described via perturbation theory. △ Less

Submitted 1 December, 2020; v1 submitted 4 September, 2020; originally announced September 2020.

Comments: 14 pages, 3 figures, plus SI

Journal ref: J. Chem. Phys. 154, 121101 (2021)

Decomposition and embedding in the stochastic $GW$ self-energy

Authors: Mariya Romanova, Vojtěch Vlček

Abstract: We present two new developments for computing excited state energies within the $GW$ approximation. First, calculations of the Green's function and the screened Coulomb interaction are decomposed into two parts: one is deterministic while the other relies on stochastic sampling. Second, this separation allows constructing a subspace self-energy, which contains dynamic correlation from only a parti… ▽ More We present two new developments for computing excited state energies within the $GW$ approximation. First, calculations of the Green's function and the screened Coulomb interaction are decomposed into two parts: one is deterministic while the other relies on stochastic sampling. Second, this separation allows constructing a subspace self-energy, which contains dynamic correlation from only a particular (spatial or energetic) region of interest. The methodology is exemplified on large-scale simulations of nitrogen-vacancy states in a periodic hBN monolayer and hBN-graphene heterostructure. We demonstrate that the deterministic embedding of strongly localized states significantly reduces statistical errors, and the computational cost decreases by more than an order of magnitude. The computed subspace self-energy unveils how interfacial couplings affect electronic correlations and identifies contributions to excited-state lifetimes. While the embedding is necessary for the proper treatment of impurity states, the decomposition yields new physical insight into quantum phenomena in heterogeneous systems. △ Less

Submitted 30 June, 2020; originally announced July 2020.

Comments: 11 pages, 3 figures

Symmetry aspects of spin-filtering in molecular junctions: hybridization and quantum interference effects

Authors: Dongzhe Li, Rajdeep Banerjee, Sourav Mondal, Ivan Maliyov, Mariya Romanova, Yannick J. Dappe, Alexander Smogunov

Abstract: Control and manipulation of electric current and, especially, its degree of spin polarization (spin filtering) across single molecules are currently of great interest in the field of molecular spintronics. We explore one possible strategy based on the modification of nanojunction symmetry which can be realized, for example, by a mechanical strain. Such modification can activate new molecular orbit… ▽ More Control and manipulation of electric current and, especially, its degree of spin polarization (spin filtering) across single molecules are currently of great interest in the field of molecular spintronics. We explore one possible strategy based on the modification of nanojunction symmetry which can be realized, for example, by a mechanical strain. Such modification can activate new molecular orbitals which were inactive before due to their orbital mismatch with the electrode's conduction states. This can result in several important consequences such as (i) quantum interference effects appearing as Fano-like features in electron transmission and (ii) the change in molecular level hybridization with the electrode's states. We argue that the symmetry change can affect very differently two majority- and minority-spin conductances and thus alter significantly the resulting spin-filtering ratio as the junction symmetry is modified. We illustrate the idea for two basic molecular junctions: Ni/benzene/Ni (perpendicular vs tilted orientations) and Ni/Si chain/Ni (zigzag vs linear chains). In both cases, one highest occupied molecular orbital (HOMO) and one lowest unoccupied molecular orbital (LUMO) (out of HOMO and LUMO doublets) are important. In particular, their destructive interference with other orbitals leads to dramatic suppression of majority-spin conductance in low-symmetry configurations. For a minority-spin channel, on the contrary, the conductance is strongly enhanced when the symmetry is lowered due to an increase in hybridization strength. We believe that our results may offer a potential route for creating molecular devices with a large on-off ratio of spin polarization via quantum interference effects. △ Less

Submitted 18 August, 2019; originally announced August 2019.

Journal ref: Phys. Rev. B 99, 115403 (2019)

Interpolation sets of algebras of generalized analytic functions

Authors: A. R. Mirotin, M. A. Romanova

Abstract: The properties of the compactness of interpolation sets of algebras of generalized analytic functions are investigated and convenient sufficient conditions for interpolation are given. The properties of the compactness of interpolation sets of algebras of generalized analytic functions are investigated and convenient sufficient conditions for interpolation are given. △ Less

Submitted 28 February, 2019; originally announced March 2019.

Comments: in Russian

MSC Class: 43A

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