New, efficient, and accurate high order derivative and dissipation operators satisfying summation by parts, and applications in three-dimensional multi-block evolutions

We construct new, efficient, and accurate high-order finite differencing operators which satisfy summation by parts. Since these operators are not uniquely defined, we consider several optimization criteria: minimizing the bandwidth, the truncation error on the boundary points, the spectral radius, or a combination of these. We examine in detail a set of operators that are up to tenth order accurate in the interior, and we surprisingly find that a combination of these optimizations can improve the operators' spectral radius and accuracy by orders of magnitude in certain cases. We also construct high-order dissipation operators that are compatible with these new finite difference operators and which are semi-definite with respect to the appropriate summation by parts scalar product. We test the stability and accuracy of these new difference and dissipation operators by evolving a three-dimensional scalar wave equation on a spherical domain consisting of seven blocks, each discretized with a structured grid, and connected through penalty boundary conditions.Comment: 16 pages, 9 figures. The files with the coefficients for the derivative and dissipation operators can be accessed by downloading the source code for the document. The files are located in the "coeffs" subdirector

On the inviscid and non-resistive limit for the equations of incompressible magnetohydrodynamics

We prove the convergence of the solutions for the incompressible homogeneous magnetohydrodynamics (MHD) system to the solutions to ideal MHD one in the inviscid and non-resistive limit, detailing the explicit convergence rates. For this study we consider a fluid occupying the whole space R3 and we assume that the viscosity effects in this fluid can be described by two different operators: the usual Laplacian operator affected by the inverse of the Reynolds number or by a viscosity operator introduced by S. I. Braginskii in 1965

On the well posedness of the Baumgarte-Shapiro-Shibata-Nakamura formulation of Einstein's field equations

We give a well posed initial value formulation of the Baumgarte-Shapiro-Shibata-Nakamura form of Einstein's equations with gauge conditions given by a Bona-Masso like slicing condition for the lapse and a frozen shift. This is achieved by introducing extra variables and recasting the evolution equations into a first order symmetric hyperbolic system. We also consider the presence of artificial boundaries and derive a set of boundary conditions that guarantee that the resulting initial-boundary value problem is well posed, though not necessarily compatible with the constraints. In the case of dynamical gauge conditions for the lapse and shift we obtain a class of evolution equations which are strongly hyperbolic and so yield well posed initial value formulations

Stability of general-relativistic accretion disks

Self-gravitating relativistic disks around black holes can form as transient structures in a number of astrophysical scenarios such as binary neutron star and black hole-neutron star coalescences, as well as the core-collapse of massive stars. We explore the stability of such disks against runaway and non-axisymmetric instabilities using three-dimensional hydrodynamics simulations in full general relativity using the THOR code. We model the disk matter using the ideal fluid approximation with a $\Gamma$-law equation of state with $\Gamma=4/3$. We explore three disk models around non-rotating black holes with disk-to-black hole mass ratios of 0.24, 0.17 and 0.11. Due to metric blending in our initial data, all of our initial models contain an initial axisymmetric perturbation which induces radial disk oscillations. Despite these oscillations, our models do not develop the runaway instability during the first several orbital periods. Instead, all of the models develop unstable non-axisymmetric modes on a dynamical timescale. We observe two distinct types of instabilities: the Papaloizou-Pringle and the so-called intermediate type instabilities. The development of the non-axisymmetric mode with azimuthal number m = 1 is accompanied by an outspiraling motion of the black hole, which significantly amplifies the growth rate of the m = 1 mode in some cases. Overall, our simulations show that the properties of the unstable non-axisymmetric modes in our disk models are qualitatively similar to those in Newtonian theory.Comment: 30 pages, 21 figure

Well-Posed Initial-Boundary Evolution in General Relativity

Maximally dissipative boundary conditions are applied to the initial-boundary value problem for Einstein's equations in harmonic coordinates to show that it is well-posed for homogeneous boundary data and for boundary data that is small in a linearized sense. The method is implemented as a nonlinear evolution code which satisfies convergence tests in the nonlinear regime and is robustly stable in the weak field regime. A linearized version has been stably matched to a characteristic code to compute the gravitational waveform radiated to infinity.Comment: 5 pages, 6 figures; added another convergence plot to Fig. 2 + minor change

Oscillations and waves in solar spicules

Since their discovery, spicules have attracted increased attention as energy/mass bridges between the dense and dynamic photosphere and the tenuous hot solar corona. Mechanical energy of photospheric random and coherent motions can be guided by magnetic field lines, spanning from the interior to the upper parts of the solar atmosphere, in the form of waves and oscillations. Since spicules are one of the most pronounced features of the chromosphere, the energy transport they participate in can be traced by the observations of their oscillatory motions. Oscillations in spicules have been observed for a long time. However the recent high-resolutions and high-cadence space and ground based facilities with superb spatial, temporal and spectral capacities brought new aspects in the research of spicule dynamics. Here we review the progress made in imaging and spectroscopic observations of waves and oscillations in spicules. The observations are accompanied by a discussion on theoretical modelling and interpretations of these oscillations. Finally, we embark on the recent developments made on the presence and role of Alfven and kink waves in spicules. We also address the extensive debate made on the Alfven versus kink waves in the context of the explanation of the observed transverse oscillations of spicule axes

Testing outer boundary treatments for the Einstein equations

Various methods of treating outer boundaries in numerical relativity are compared using a simple test problem: a Schwarzschild black hole with an outgoing gravitational wave perturbation. Numerical solutions computed using different boundary treatments are compared to a `reference' numerical solution obtained by placing the outer boundary at a very large radius. For each boundary treatment, the full solutions including constraint violations and extracted gravitational waves are compared to those of the reference solution, thereby assessing the reflections caused by the artificial boundary. These tests use a first-order generalized harmonic formulation of the Einstein equations. Constraint-preserving boundary conditions for this system are reviewed, and an improved boundary condition on the gauge degrees of freedom is presented. Alternate boundary conditions evaluated here include freezing the incoming characteristic fields, Sommerfeld boundary conditions, and the constraint-preserving boundary conditions of Kreiss and Winicour. Rather different approaches to boundary treatments, such as sponge layers and spatial compactification, are also tested. Overall the best treatment found here combines boundary conditions that preserve the constraints, freeze the Newman-Penrose scalar Psi_0, and control gauge reflections.Comment: Modified to agree with version accepted for publication in Class. Quantum Gra

Awareness, Co-operation, Tackling to Stop Sexual Bullying: An empowerment pack for young people and the people working with them (The ACT Pack).

Available in English, Bulgarian, Slovenian and Italia

Telomere length shortening is associated with treatment-free remission in chronic myeloid leukemia patients

We studied telomere length in 32 CML patients who discontinued imatinib after achieving complete molecular remission and 32 age-sex-matched controls. The relative telomere length (RTL) was determined by q-PCR as the telomere to single copy gene (36B4) ratio normalized to a reference sample (K-562 DNA). Age-corrected RTL (acRTL) was also obtained. The 36-month probability of treatment-free remission (TFR) was 59.4 %. TFR patients showed shorter acRTL compared to relapsed (mean ± SD = 0.01 ± 0.14 vs 0.20 ± 0.21; p = 0.01). TFR was significantly higher in CML patients with acRTL ≤0.09 (78.9 vs 30.8 %, p = 0.002). CML stem cells harboring longer telomeres possibly maintain a proliferative potential after treatment discontinuation