WhiskyMHD: a new numerical code for general relativistic magnetohydrodynamics

The accurate modelling of astrophysical scenarios involving compact objects and magnetic fields, such as the collapse of rotating magnetized stars to black holes or the phenomenology of gamma-ray bursts, requires the solution of the Einstein equations together with those of general-relativistic magnetohydrodynamics. We present a new numerical code developed to solve the full set of general-relativistic magnetohydrodynamics equations in a dynamical and arbitrary spacetime with high-resolution shock-capturing techniques on domains with adaptive mesh refinements. After a discussion of the equations solved and of the techniques employed, we present a series of testbeds carried out to validate the code and assess its accuracy. Such tests range from the solution of relativistic Riemann problems in flat spacetime, over to the stationary accretion onto a Schwarzschild black hole and up to the evolution of oscillating magnetized stars in equilibrium and constructed as consistent solutions of the coupled Einstein-Maxwell equations.Comment: minor changes to match the published versio

Numerical 3+1 general relativistic magnetohydrodynamics: a local characteristic approach

We present a general procedure to solve numerically the general relativistic magnetohydrodynamics (GRMHD) equations within the framework of the 3+1 formalism. The work reported here extends our previous investigation in general relativistic hydrodynamics (Banyuls et al. 1997) where magnetic fields were not considered. The GRMHD equations are written in conservative form to exploit their hyperbolic character in the solution procedure. All theoretical ingredients necessary to build up high-resolution shock-capturing schemes based on the solution of local Riemann problems (i.e. Godunov-type schemes) are described. In particular, we use a renormalized set of regular eigenvectors of the flux Jacobians of the relativistic magnetohydrodynamics equations. In addition, the paper describes a procedure based on the equivalence principle of general relativity that allows the use of Riemann solvers designed for special relativistic magnetohydrodynamics in GRMHD. Our formulation and numerical methodology are assessed by performing various test simulations recently considered by different authors. These include magnetized shock tubes, spherical accretion onto a Schwarzschild black hole, equatorial accretion onto a Kerr black hole, and magnetized thick accretion disks around a black hole prone to the magnetorotational instability.Comment: 18 pages, 8 figures, submitted to Ap

Non-axisymmetric relativistic Bondi-Hoyle accretion onto a Kerr black hole

In our program of studying numerically the so-called Bondi-Hoyle accretion in the fully relativistic regime, we present here first results concerning the evolution of matter accreting supersonically onto a rotating (Kerr) black hole. These computations generalize previous results where the non-rotating (Schwarzschild) case was extensively considered. We parametrize our initial data by the asymptotic conditions for the fluid and explore the dependence of the solution on the angular momentum of the black hole. Towards quantifying the robustness of our numerical results, we use two different geometrical foliations of the black hole spacetime, the standard form of the Kerr metric in Boyer-Lindquist coordinates as well as its Kerr-Schild form, which is free of coordinate singularities at the black hole horizon. We demonstrate some important advantages of using such horizon adapted coordinate systems. Our numerical study indicates that regardless of the value of the black hole spin the final accretion pattern is always stable, leading to constant accretion rates of mass and momentum. The flow is characterized by a strong tail shock, which, unlike the Schwarzschild case, is increasingly wrapped around the central black hole as the hole angular momentum increases. The rotation induced asymmetry in the pressure field implies that besides the well known drag, the black hole will experience also a lift normal to the flow direction. This situation exhibits some analogies with the Magnus effect of classical fluid dynamics.Comment: 33 pages, 20 figures, submited to MNRA

Magnetohydrodynamics in full general relativity: Formulation and tests

A new implementation for magnetohydrodynamics (MHD) simulations in full general relativity (involving dynamical spacetimes) is presented. In our implementation, Einstein's evolution equations are evolved by a BSSN formalism, MHD equations by a high-resolution central scheme, and induction equation by a constraint transport method. We perform numerical simulations for standard test problems in relativistic MHD, including special relativistic magnetized shocks, general relativistic magnetized Bondi flow in stationary spacetime, and a longterm evolution for self-gravitating system composed of a neutron star and a magnetized disk in full general relativity. In the final test, we illustrate that our implementation can follow winding-up of the magnetic field lines of magnetized and differentially rotating accretion disks around a compact object until saturation, after which magnetically driven wind and angular momentum transport inside the disk turn on.Comment: 28 pages, to be published in Phys. Rev.

Production mechanisms and single-spin asymmetry for kaons in high energy hadron-hadron collisions

Direct consequences on kaon production of the picture proposed in a recent Letter and subsequent publications are discussed. Further evidence supporting the proposed picture is obtained. Comparison with the data for the inclusive cross sections in unpolarized reactions is made. Quantitative results for the left-right asymmetry in single-spin processes are presented.Comment: 10 pages, 2 Postscript figure

Amorphous NiCu Thin Films Sputtered on TiO2 Nanotube Arrays: A Noble-Metal Free Photocatalyst for Hydrogen Evolution

In this work, NiCu co-catalysts on TiO2 are studied for photocatalytic hydrogen evolution. NiCu co-catalyst films are deposited at room temperature by argon plasma sputtering on high aspect-ratio anodic TiO2 nanotubes. To tune the Ni : Cu atomic ratio, alloys of various compositions were used as sputtering targets. Such co-catalyst films are found to be amorphous with small nanocrystalline domains. A series of parameters is investigated, i. e., i) Ni : Cu relative ratio in the sputtered films, ii) NiCu film thickness, and iii) thickness of the TiO2 nanotube layers. The highest photocatalytic activity is obtained with 8 μm long TiO2 nanotubes, sputter-coated with a 10 nm-thick NiCu films with a 1 : 1 Ni : Cu atomic ratio. This photocatalyst reaches a stable hydrogen evolution rate of 186 μL h−1 cm−2, 4.6 and 3 times higher than that of Ni- and Cu-TiO2, respectively, demonstrating a synergistic co-catalytic effect of Ni and Cu in the alloy co-catalyst film

Nonlocal calculation for nonstrange dibaryons and tribaryons

We study the possible existence of nonstrange dibaryons and tribaryons by solving the bound-state problem of the two- and three-body systems composed of nucleons and deltas. The two-body systems are $NN$, $N\Delta$, and $\Delta\Delta$, while the three-body systems are $NNN$, $NN\Delta$, $N\Delta\Delta$, and $\Delta\Delta\Delta$. We use as input the nonlocal $NN$, $N\Delta$, and $\Delta\Delta$ potentials derived from the chiral quark cluster model by means of the resonating group method. We compare with previous results obtained from the local version based on the Born-Oppenheimer approximation.Comment: 19 pages. To be published in Physical Review