We have programmed a full simulation of the open field lines above the polar cap of a magnetized pulsar, using time dependent Particle In Cell (PIC) numerical code (Birdsall and Langdon 1991). We consider the case of free ejection of electrons from the star surface, a Space Charge Limited Flow (SCLF) model. We report here the first results of the simulation. Electrons are accelerating along the open field lines to the flat space-time SCLF maximum Lorentz factor prediction, with oscillation pattern. Than we add the General Relativistic Frame Dragging correction that provide the particles the high Lorentz factor (10^6 - 10^7) needed to initiate pair production. The electrons accelerate according to the analytic expressions given in Muslimov and Tsygan 1992, and Muslimov and Harding 1997, with oscillation pattern. Electron-positron pair production is now being programmed, using the cross sections appears in the literature, and Monte-Carlo code. After completing this stage, we will automatically get: a) the positron return current (thus we could calculate the polar cap heating and the X-ray emission). b). The photons and electrons observed spectrum (photons and electrons that escape the magnetosphere after the cascade). c). The pulsar death line (pulsars with not enough pair production). d). The PFF height (pair formation front location). Those results will be report in a different paper.

We compute the structure and degree of neutronization of general relativistic magnetohydrodynamic (GRMHD) outflows originating from the inner region of neutrino-cooled disks. We consider both, outflows expelled from a hydrostatic disk corona and outflows driven by disk turbulence. We show that in outflows driven thermally from a static disk the electron fraction quickly evolves to its equilibrium value which is dominated by neutrino capture. Those outflows are generally proton rich and, under certain conditions, can be magnetically dominated. They may also provide sites for effective production of 56Ni. Centrifugally driven outflows and outflows driven by disk turbulence, on the other hand, can preserve the large in-disk neutron excess. Those outflows are, quite generally, subrelativistic by virtue of the large mass flux driven by the additional forces.