Magnetohydrodynamic turbulence in warped accretion discs

Warped, precessing accretion discs appear in a range of astrophysical systems, for instance the X-ray binary Her X-1 and in the active nucleus of NGC4258. In a warped accretion disc there are horizontal pressure gradients that drive an epicyclic motion. We have studied the interaction of this epicyclic motion with the magnetohydrodynamic turbulence in numerical simulations. We find that the turbulent stress acting on the epicyclic motion is comparable in size to the stress that drives the accretion, however an important ingredient in the damping of the epicyclic motion is its parametric decay into inertial waves.Comment: to appear in the proceedings of the 20th Texas Symposium on Relativistic Astrophysics, J. C. Wheeler & H. Martel (eds.

The response of a turbulent accretion disc to an imposed epicyclic shearing motion

We excite an epicyclic motion, whose amplitude depends on the vertical position, $z$, in a simulation of a turbulent accretion disc. An epicyclic motion of this kind may be caused by a warping of the disc. By studying how the epicyclic motion decays we can obtain information about the interaction between the warp and the disc turbulence. A high amplitude epicyclic motion decays first by exciting inertial waves through a parametric instability, but its subsequent exponential damping may be reproduced by a turbulent viscosity. We estimate the effective viscosity parameter, $\alpha_{\rm v}$, pertaining to such a vertical shear. We also gain new information on the properties of the disc turbulence in general, and measure the usual viscosity parameter, $\alpha_{\rm h}$, pertaining to a horizontal (Keplerian) shear. We find that, as is often assumed in theoretical studies, $\alpha_{\rm v}$ is approximately equal to $\alpha_{\rm h}$ and both are much less than unity, for the field strengths achieved in our local box calculations of turbulence. In view of the smallness ($\sim 0.01$) of $\alpha_{\rm v}$ and $\alpha_{\rm h}$ we conclude that for $\beta = p_{\rm gas}/p_{\rm mag} \sim 10$ the timescale for diffusion or damping of a warp is much shorter than the usual viscous timescale. Finally, we review the astrophysical implications.Comment: 12 pages, 18 figures, MNRAS accepte

Convective intensification of magnetic fields in the quiet Sun

Kilogauss-strength magnetic fields are often observed in intergranular lanes at the photosphere in the quiet Sun. Such fields are stronger than the equipartition field B_e, corresponding to a magnetic energy density that matches the kinetic energy density of photospheric convection, and comparable with the field B_p that exerts a magnetic pressure equal to the ambient gas pressure. We present an idealised numerical model of three-dimensional compressible magnetoconvection at the photosphere, for a range of values of the magnetic Reynolds number. In the absence of a magnetic field, the convection is highly supercritical and is characterised by a pattern of vigorous, time-dependent, “granular” motions. When a weak magnetic field is imposed upon the convection, magnetic flux is swept into the convective downflows where it forms localised concentrations. Unless this process is significantly inhibited by magnetic diffusion, the resulting fields are often much greater than B_e, and the high magnetic pressure in these flux elements leads to their being partially evacuated. Some of these flux elements contain ultra-intense magnetic fields that are significantly greater than B_p. Such fields are contained by a combination of the thermal pressure of the gas and the dynamic pressure of the convective motion, and they are constantly evolving. These ultra-intense fields develop owing to nonlinear interactions between magnetic fields and convection; they cannot be explained in terms of “convective collapse” within a thin flux tube that remains in overall pressure equilibrium with its surroundings

On the classification and quantification of crystal defects after energetic bombardment by machine learned molecular dynamics simulations

The analysis of the damage on plasma facing materials (PFM), due to their direct interaction with the plasma environment, is needed to build the next generation of nuclear fusion reactors. After systematic analyses of numerous materials over the last decades, tungsten has become the most promising candidate for a nuclear fusion reactor. In this work, we perform molecular dynamics (MD) simulations using a machine learned interatomic potential, based on the Gaussian Approximation Potential framework, to model better neutron bombardment mechanisms in pristine W lattices. The MD potential is trained to reproduce realistic short-range dynamics, the liquid phase, and the material recrystallization, which are important for collision cascades. The formation of point defects is quantified and classified by a descriptor vector (DV) based method, which is independent of the sample temperature and its constituents, requiring only modest computational resources. The locations of vacancies are calculated by the k-d-tree algorithm. The analysis of the damage in the W samples is compared to results obtained by Finnis–Sinclair and Tersoff–Ziegler–Biersack–Littmark potentials, at a sample temperature of 300 K and a primary knock-on atom (PKA) energy range of 0.5–10 keV, where a good agreement with the reported number of Frenkel pair is observed. Our results provide information about the advantages and limits of the machine learned MD simulations with respect to the standard ones. The formation of dumbbell and crowdion defects as a function of PKA energy were identified and distinguished by our DV method.Peer reviewe

Dynamics of Magnetic Flux Elements in the Solar Photosphere

The interaction of magnetic fields and convection is investigated in the context of the coronal heating problem. We study the motions of photospheric magnetic elements using filtergrams obtained at the Swedish Vacuum Solar Telescope at La Palma. We use potential-field modeling to extrapolate the magnetic and velocity fields to larger height. We find that the velocity in the chromosphere can be locally enhanced at the separatrix surfaces between neighboring flux tubes. The predicted velocities are several km/s, significantly larger than those of the photospheric flux tubes, which may have important implications for coronal heating. sComment: submitted to ApJ, 21 pages, 10 figure

GALEX, Optical and IR Light Curves of MQ Dra: UV Excesses at Low Accretion Rates

Ultraviolet light curves constructed from NUV and FUV detectors on GALEX reveal large amplitude variations during the orbital period of the Low Accretion Rate Polar MQ Dra (SDSSJ1553+55). This unexpected variation from a UV source is similar to that seen and discussed in the Polar EF Eri during its low state of accretion, even though the accretion rate in MQ Dra is an order of magnitude lower than even the low state of EF Eri. The similarity in phasing of the UV and optical light curves in MQ Dra imply a similar location for the source of light. We explore the possibilities of hot spots and cyclotron emission with simple models fit to the UV, optical and IR light curves of MQ Dra. To match the GALEX light curves with a single temperature circular hot spot requires different sizes of spots for the NUV and FUV, while a cyclotron model that can produce the optical harmonics with a magnetic field near 60 MG requires multipoles with fields > 200 MG to match the UV fluxes.Comment: accepted for ApJ; 15 pages, 7 tables, 8 fig

The Quiet-Sun Photosphere and Chromosphere

The overall structure and the fine structure of the solar photosphere outside active regions are largely understood, except possibly important roles of a turbulent near-surface dynamo at its bottom, internal gravity waves at its top, and small-scale vorticity. Classical 1D static radiation-escape modelling has been replaced by 3D time-dependent MHD simulations that come closer to reality. The solar chromosphere, in contrast, remains ill-understood although its pivotal role in coronal mass and energy loading makes it a principal research area. Its fine structure defines its overall structure, so that hard-to-observe and hard-to-model small-scale dynamical processes are the key to understanding. However, both chromospheric observation and chromospheric simulation presently mature towards the required sophistication. The open-field features seem of greater interest than the easier-to-see closed-field features.Comment: Accepted for special issue "Astrophysical Processes on the Sun" of Phil. Trans. Royal Soc. A, ed. C. Parnell. Note: clicking on the year in a citation opens the corresponding ADS abstract page in the browse

Optimization of single crystal mirrors for ITER diagnostics

Diagnostic mirrors are planned to be used in all optical diagnostics in ITER. Degradation of mirrors due to e.g. deposition of plasma impurities will hamper the entire performance of affected diagnostics. in situ mirror cleaning by plasma sputtering is presently envisaged for the recovery of contaminated mirrors. There are observations showing a signature of sputtering dependence on crystal orientation. Should such a dependence exist, the sputtering of single crystal mirrors could be minimized, thus prolonging a mirror lifetime. Four single crystal molybdenum mirrors with different orientations were produced to study the effect of crystal orientation on sputtering. Mirrors were exposed to argon plasma under identical plasma conditions relevant to those expected in the mirror cleaning systems of ITER. The energy of impinging ions was about 60 eV. The amount of sputtered material corresponded to about a hundred mirror cleaning cycles in argon. Plasma exposures did not affect the mirror reflectivity. The maximum decrease of specular reflectivity did not exceed 5% at 250 nm. The mirrors with orientations [110]/[101] demonstrated up to 42% less sputtering than the mirrors with other crystal orientations. These findings outline the advantage of a favorable crystal orientation for a cleaning of heavy contaminants from ITER mirrors.Peer reviewe

Atomic-scale effects behind structural instabilities in Si lamellae during ion beam thinning

Peer reviewe