Researcher profile

Marina M. Romanova

Marina M. Romanova contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
13works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

13 published item(s)

preprint2015arXiv

Asymmetric MHD Outflows/Jets from Accreting T Tauri Stars

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.

preprint2014arXiv

Counter-Rotating Accretion Discs

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.

preprint2013arXiv

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

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 &#34;drift&#34; 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.

preprint2013arXiv

Propeller outflows from an MRI disc

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.

preprint2013arXiv

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

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&#39; 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.

preprint2012arXiv

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

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&#39; 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.

preprint2011arXiv

Accretion onto Stars with Octupole Magnetic Fields: Matter Flow, Hot Spots and Phase Shifts

We show results of global 3D MHD simulations of accretion onto stars with superposition of the dipole and octupole fields, where we vary the ratio between components. Simulations show that if octupolar field strongly dominates at the disc-magnetosphere boundary, then matter flows into the ring-like octupolar poles, forming ring-shape spots at the surface of the star above and below equator. The light-curves are complex and may have two peaks per period. In case where the dipole field dominates, matter accretes in two ordered funnel streams towards poles, however the polar spots are meridionally-elongated due to the action of the octupolar component. In the case when the fields are of similar strengths, both, polar and belt-like spots are present. In many cases the light-curves show the evidence of complex fields, excluding the cases of small inclinations angles, where sinusoidal light-curve %is observed and `hides&#39; the information about the field complexity. We also propose new mechanisms of phase shift in stars with complex magnetic fields. We suggest that the phase shifts can be connected with: (1) temporal variation of the star&#39;s intrinsic magnetic field and subsequent redistribution of main magnetic poles; (2) variation of the accretion rate, which causes the disc to interact with the magnetic fields associated with different magnetic moments. We use our model to demonstrate these phase shift mechanisms, and we discuss possible applications of these mechanisms to accreting millisecond pulsars and young stars.

preprint2011arXiv

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

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

preprint2011arXiv

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

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

preprint2011arXiv

MRI-driven Accretion onto Magnetized stars: Axisymmetric MHD Simulations

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

preprint2011arXiv

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

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

preprint2010arXiv

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

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

preprint2009arXiv

Possible Quasi-Periodic Oscillations from Unstable Accretion: 3D MHD Simulations

We investigate the photometric variability of magnetized stars, particularly neutron stars, accreting through a magnetic Rayleigh-Taylor-type instability at the disk-magnetosphere interface, and compare it with the variability during stable accretion, with the goal of looking for possible quasi-periodic oscillations. The lightcurves during stable accretion show periodicity at the star&#39;s frequency and sometimes twice that, due to the presence of two funnel streams that produce antipodal hotspots near the magnetic poles. On the other hand, lightcurves during unstable accretion through tongues penetrating the magnetosphere are more chaotic due to the stochastic behaviour of the tongues, and produce noisier power spectra. However, the power spectra do show some signs of quasi-periodic variability. Most importantly, the rotation frequency of the tongues and the resulting hotspots is close to the inner-disk orbital frequency, except in the most strongly unstable cases. There is therefore a high probability of observing QPOs at that frequency in longer simulations. In addition, the lightcurves in the unstable regime show periodicity at the star&#39;s rotation frequency in many of the cases investigated here, again except in the most strongly unstable cases which lack funnel flows and the resulting antipodal hotspots. The noisier power spectra result in the fractional rms amplitudes of the Fourier peaks being smaller. We also study in detail the effect of the misalignment angle between the rotation and magnetic axes of the star on the variability, and find that at misalignment angles $\gtrsim 25^\circ$, the star&#39;s period always appears in the lightcurves.