

CoCoNuT Code



Program:
The meeting will take place at the salon de grados
of the Mathematics Faculty in
the Burjassot campus of the University
of Valencia.
You can click on the title of the talks to see the abstract.
Invited talks: 45' + questions/discussion
Contributed talks: 20' + questions/discussion

Morning session (9:30  13:00) 


Antonio Marquina 
Anomalous wave structure in magnetized materials described by nonconvex equations of state We analyze the anomalous wave structure appearing in flow dynamics under the influence of magnetic field in materials described by nonideal equations of state. We consider the hyperbolic system of magnetohydrodynamics equations closed by a general equation of state (EOS) and propose a complete spectral decomposition of the fluxes that allows us to derive an expression of the nonlinearity factor as the mathematical tool to determine the nature of the wave phenomena. We prove that the possible formation of nonclassical wave structure is determined by the thermodynamic properties of the material. The analytical expression of the nonlinearity factor allows us to determine the specific amount of magnetic field necessary to prevent formation of complex structure induced by phase transition in the material. We illustrate our analytical approach by considering two nonconvex EOS that exhibit phase transitions and anomalous behavior in the evolution. We present numerical experiments validating the analysis performed through a set of onedimensional Riemann problems. (Work in collaboration with Dr. Susana Serna, Phys. Fluids, vol. 26, 016101, (2014))

45'+ 
Yuichiro Sekiguchi 
On Compact binary mergers, Short GRBs, and rprocess nucleosynthesis 
20'+ 

Coffee break (10:5011:30) 


Oliver Just 
Combined Analysis of Outflows from NeutronStar Mergers The astrophysical site for the production of rprocess elements is still mysterious. Recent studies suggest that instead of the long favored corecollapse supernovae the outflows from compact binary (i.e. double neutronstar, NSNS, or neutronstar blackhole, NSBH) mergers could provide suitable conditions to allow for the rprocess. In general, during each stage of the merger different types of outflows with individual nucleosynthesis signatures can be launched. I will present the results of our recent study in which we conducted a first comprehensive analysis of compact binary mergers that takes into account the prompt ejecta from both NSNS and NSBH mergers consistently together with their relic BHtorus systems. 
20'+ 
Riccardo Ciolfi 
Magnetically driven winds from differentially rotating neutron stars and Xray afterglows of short gammaray bursts Besides being among the most promising sources of gravitational waves, merging neutron star binaries also represent a leading scenario to explain the phenomenology of short gammaray bursts (SGRBs). Recent observations have revealed a large subclass of SGRBs with roughly constant luminosity in their Xray afterglows, lasting 10–10^4 s. These features are generally taken as evidence of a longlived central engine powered by the magnetic spindown of a uniformly rotating, magnetized object. We propose a scenario in which the central engine powering the early Xray afterglow emission is a differentially rotating hypermassive or supramassive neutron star (HMNS or SMNS). This emission is associated with a quasiisotropic and baryonloaded wind driven by the magnetic field, which is builtup through differential rotation. Our model is supported by longterm, threedimensional, generalrelativistic, and ideal magnetohydrodynamic simulations, showing that this isotropic emission is a very robust feature. We show that our results are compatible with the timescales and luminosities of the observed Xray afterglows. 
20'+ 
Daniel Siegel 
Short gammaray bursts in the "timereversal" scenario Short gammaray bursts (SGRBs) are among the most luminous explosions in the Universe and their origin still remains uncertain. Observational evidence favors the association with binary neutron star or neutron starblack hole (NSBH) binary mergers. Leading models relate SGRBs to a relativistic jet launched by the BHtorus system resulting from the merger. However, recent observations have revealed a large fraction of SGRB events accompanied by Xray afterglows with durations $10^210^5$ s, suggesting continuous energy injection from a longlived central engine, which is incompatible with the short ($lt 1$ s) accretion timescale of a BHtorus system. The formation of a supramassive NS, resisting the collapse on much longer spindown timescales, can explain these afterglow durations, but leaves serious doubts on whether a relativistic jet can be launched at merger. Here we present a novel scenario accommodating both aspects, where the SGRB is produced after the collapse of a supramassive NS. Early differential rotation and subsequent spindown emission generate an optically thick environment around the NS consisting of a photonpair nebula and an outer shell of baryonloaded ejecta. While the jet easily drills through this environment, spindown radiation diffuses outwards on much longer timescales and accumulates a delay that allows the SGRB to be observed before (part of) the longlasting Xray signal. By analyzing diffusion timescales for a wide range of physical parameters, we find delays that can generally reach $10^5$ s, compatible with observations. The success of this fundamental test makes this "timereversal" scenario an attractive alternative to current SGRB models. 
20'+ 

Lunch break (13:0015:00) 



Afternoon session (15:00  18:00) 


Vassilios Mewes 
Measuring BH spin in numerical relativity simulations We use the LandauLifshitz Pseudotensor to measure the spin of Kerr black holes with arbitrary initial spin orientation in BBH simulations. We compare the method to measuring the spin in the quasilocal framework using an approximate Killing vector on the horizon and a method using flat space coordinate rotational Killing vectors for the calculation of the spin. 
20'+ 
David Hilditch 
Recent progress in the numerical treatment of collapsing gravitational waves I will present simulations of collapsing gravitational waves obtained using two methods. First, using the moving puncture gauge, I show that it is possible to evolve certain types of initial data through apparent horizon formation. Secondly I will present a new pseudospectral code which we are using to tackle the same initial data. 
20'+ 

Coffee break (16:0016:30) 



Morning session (9:30  13:00) 


Manuel Castro 
An overview of pathconservative numerical schemes for balance laws and nonconservative hyperbolic systems: main difficulties and efficient Riemann solvers. In this talk we present a general overview of pathconservative numerical schemes for balance laws and nonconservative hyperbolic systems introduced by Parés. In particular, we show the main difficulties related to the presence of nonconservative products concerning the definition of weak solutions and the derivation of numerical schemes. Next, we present a class of incomplete Riemann solvers, whose viscosity matrix are defined in terms of an appropriate real function $R(x)$, that approximates the function $x$. The resulting approximate Riemann solvers are incomplete, in the sense that we do not use the complete spectral decomposition of the system, and only some information about the maximum in absolute value of the characteristic speeds are needed. Finally, some numerical tests will be presented. 
45'+ 
Michael Gabler 
RCS in Magnetar Magnetospheres We present a new numerical tool to study the resonant cyclotron scattering to calculate the modulation of the emission of highly magnetized neutron stars (magnetars). The scattering of the photons is calculated with a MonteCarlo approach. In order to obtain a realistic prescription for the scattering targets, i.e. electrons and positrons that are accelerated along the magnetic field lines, we need to follow their motion. We follow the particles with a particleinaline code that calculates their trajectories by taking the interaction with the photons into account. Our final aim is to apply the new code to describe the quasiperiodic oscillations (QPOs) observed in magnetar giant flares. 
20'+ 

Coffee break (10:5011:30) 


Alejandro TorresFornés 
Accretion onto neutron stars and the hidden magnetic field model We present 1D and 2D relativistic MHD simulations of the accretion of matter onto neutron stars. We study the consequences of the accretion process in the distribution of the magnetic field across the magnetosphere of the star. To this end we explore different accretion rates, fluid composition (through various equations of state) and magnetic field distributions, in order to test the viability of the hidden magnetic field model. Such model has been put forward to explain the evolution of Xray sources with low values ($\lt 10^11$ G) of the magnetic field. In this scenario, the reverse supernova shock falling back onto the neutron star compresses the magnetic field, a process that can eventually bury the field inside the neutron star crust. 
20'+ 
Debarati Chatterjee 
Magnetic field effects in neutron stars Previously there have been attempts to determine the structure of neutron stars in strong magnetic fields. However, there is no such study which simultaneously includes the effect of the magnetic field on the equation of state, general relativistic aspects as well as the pressure anisotropy caused by the breaking of spherical symmetry of the star by the field. We present a selfconsistent study of the structure of a neutron star in strong magnetic fields, including all the magnetic field effects. We solve Einstein’s equations in an axisymmetric metric, which is determined self consistently from the axisymmetric energymomentum tensor of the star, and solve it numerically. 
20'+ 
David Zorío 
High order boundary extrapolation technique for finite difference methods on complex domains with cartesian meshes The application of suitable numerical boundary conditions for hyperbolic conservation laws on domains with complex geometry has become a problem with certain difficulty that has been tackled in different ways according to the nature of the numerical methods and mesh type. We present an extrapolation technique on structured Cartesian meshes (which, as opposed to nonstructured meshes, cannot be adapted to the morphology of the domain boundary) of the information in the interior of the computational domain to ghost cells. This technique is based on the application of Lagrange interpolation with a previous filter for the detection of discontinuities that permits a data dependent extrapolation, with higher order at smooth regions and essentially nonoscillatory properties near discontinuities. 
20'+ 

Lunch break (13:0015:00) 



Afternoon session (15:00  18:00) 


Isabel CorderoCarrión 
Minimally Implicit RungeKutta methods I: Motivation and derivation of the schemes. In this talk I will present some RungeKutta methods, named as Minimally Implicit RungeKutta (MIRK) methods, designed to integrate in time the resistive (non ideal) relativistic magnetohydrodynamic equations. This system of equations contains stiff source terms in the evolution of the electric field. I will comment on the motivation behind these methods and the derived schemes for several orders of convergence. 
20'+ 
MiguelÁngel Aloy 
Minimally Implicit RungeKutta methods II: Numerical tests. Following the talk of I. CorderoCarrión, I will show the numerical performance of the Minimally Implicit RungeKutta (MIRK) methods on the system of resistive relativistic magnetohydrodynamics. 
20'+ 

Coffee break (16:0016:30) 



Conference dinner (~21:30) 



Morning session (9:30  13:00) 


Edouard Audit 
TBA 
45'+ 
Bruno Peres 
Neutrino treatment in CoCoNuT with an M1 scheme Recently, a general relativistic and energy dependent M1 scheme has been developed for CoCoNuT. I will report on its current state and give some technical details. 
20'+ 

Coffee break (10:5011:30) 


Martin Obergaulinger 
Amplification of magnetic fields in core collapse The precollapse magnetic fields of a supernova core can, apart from compression during collapse, be amplified by hydrodynamic instabilities developing in the protoneutron star and the postshock region. As a consequence, complex field geometries can develop and, for significantly strong field, an explosion can be facilitated. I will present simulations of the evolution of magnetised cores including neutrino transport in the twomoment approach. 
20'+ 
Miguel Marques 
Numerical model of stationary axisymmetric stars with hot matter Finite temperature effects are expected to be important in several astrophysical scenarios (e.g. corecollapse supernovae, cooling of protoneutron stars). In this talk, we will present a model for stationary axisymmetric stars with nonbarotropic equations of state. We propose a numerical scheme to find exact solutions of the equilibrium equations for stars with nonconstant entropy profiles. We will present preliminary results for a simple analytic equation of state. Finally, we will discuss the possible applications of this model. 
20'+ 
Pablo CerdaDuran 
Status of CoCoNuT I will report on the present status of CoCoNuT. I will focus in the development of a MPIparallel version and the possibility of going to an open source model. 
20'+ 

Lunch break (13:0015:00) 



Free afternoon (rooms available for smaller meetings) 



