New Measure of the Dissipation Region in Collisionless Magnetic Reconnection

A new measure to identify a small-scale dissipation region in collisionless magnetic reconnection is proposed. The energy transfer from the electromagnetic field to plasmas in the electron's rest frame is formulated as a Lorentz-invariant scalar quantity. The measure is tested by two-dimensional particle-in-cell simulations in typical configurations: symmetric and asymmetric reconnection, with and without the guide field. The innermost region surrounding the reconnection site is accurately located in all cases. We further discuss implications for nonideal MHD dissipation

On Electrostatic Positron Acceleration In The Accretion Flow Onto Neutron Stars

As first shown by Shvartsman (1970), a neutron star accreting close to the Eddington limit must acquire a positive charge in order for electrons and protons to move at the same speed. The resulting electrostatic field may contribute to accelerating positrons produced near the star surface in conjunction with the radiative force. We reconsider the balance between energy gains and losses, including inverse Compton (IC), bremsstrahlung and non--radiative scatterings. It is found that, even accounting for IC losses only, the maximum positron energy never exceeds $\approx 400$ keV. The electrostatic field alone may produce energies $\approx 50$ keV at most. We also show that Coulomb collisions and annihilation with accreting electrons severely limit the number of positrons that escape to infinity.Comment: 9 pages plus 3 postscript figures, to be published in Ap

Equation of State in Relativistic Magnetohydrodynamics: variable versus constant adiabatic index

The role of the equation of state for a perfectly conducting, relativistic magnetized fluid is the main subject of this work. The ideal constant $\Gamma$-law equation of state, commonly adopted in a wide range of astrophysical applications, is compared with a more realistic equation of state that better approximates the single-specie relativistic gas. The paper focus on three different topics. First, the influence of a more realistic equation of state on the propagation of fast magneto-sonic shocks is investigated. This calls into question the validity of the constant $\Gamma$-law equation of state in problems where the temperature of the gas substantially changes across hydromagnetic waves. Second, we present a new inversion scheme to recover primitive variables (such as rest-mass density and pressure) from conservative ones that allows for a general equation of state and avoids catastrophic numerical cancellations in the non-relativistic and ultrarelativistic limits. Finally, selected numerical tests of astrophysical relevance (including magnetized accretion flows around Kerr black holes) are compared using different equations of state. Our main conclusion is that the choice of a realistic equation of state can considerably bear upon the solution when transitions from cold to hot gas (or viceversa) are present. Under these circumstances, a polytropic equation of state can significantly endanger the solution.Comment: 14 pages, 14 figure

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.

Nonequilibrium corrections in the pressure tensor due to an energy flux

The physical interpretation of the nonequilibrium corrections in the pressure tensor for radiation submitted to an energy flux obtained in some previous works is revisited. Such pressure tensor is shown to describe a moving equilibrium system but not a real nonequilibrium situation.Comment: 4 pages, REVTeX, Brief Report to appear in PRE Dec 9

Mechanisms of endothelial cell dysfunction in cystic fibrosis

Although cystic fibrosis (CF) patients exhibit signs of endothelial perturbation, the functions of the cystic fibrosis conductance regulator (CFTR) in vascular endothelial cells (EC) are poorly defined. We sought to uncover biological activities of endothelial CFTR, relevant for vascular homeostasis and inflammation. We examined cells from human umbilical cords (HUVEC) and pulmonary artery isolated from non-cystic fibrosis (PAEC) and CF human lungs (CF-PAEC), under static conditions or physiological shear. CFTR activity, clearly detected in HUVEC and PAEC, was markedly reduced in CF-PAEC. CFTR blockade increased endothelial permeability to macromolecules and reduced trans‑endothelial electrical resistance (TEER). Consistent with this, CF-PAEC displayed lower TEER compared to PAEC. Under shear, CFTR blockade reduced VE-cadherin and p120 catenin membrane expression and triggered the formation of paxillin- and vinculin-enriched membrane blebs that evolved in shrinking of the cell body and disruption of cell-cell contacts. These changes were accompanied by enhanced release of microvesicles, which displayed reduced capability to stimulate proliferation in recipient EC. CFTR blockade also suppressed insulin-induced NO generation by EC, likely by inhibiting eNOS and AKT phosphorylation, whereas it enhanced IL-8 release. Remarkably, phosphodiesterase inhibitors in combination with a β2 adrenergic receptor agonist corrected functional and morphological changes triggered by CFTR dysfunction in EC. Our results uncover regulatory functions of CFTR in EC, suggesting a physiological role of CFTR in the maintenance EC homeostasis and its involvement in pathogenetic aspects of CF. Moreover, our findings open avenues for novel pharmacology to control endothelial dysfunction and its consequences in CF

Conservative 3+1 General Relativistic Variable Eddington Tensor Radiation Transport Equations

We present conservative 3+1 general relativistic variable Eddington tensor radiation transport equations, including greater elaboration of the momentum space divergence (that is, the energy derivative term) than in previous work. These equations are intended for use in simulations involving numerical relativity, particularly in the absence of spherical symmetry. The independent variables are the lab frame coordinate basis spacetime position coordinates and the particle energy measured in the comoving frame. With an eye towards astrophysical applications---such as core-collapse supernovae and compact object mergers---in which the fluid includes nuclei and/or nuclear matter at finite temperature, and in which the transported particles are neutrinos, we pay special attention to the consistency of four-momentum and lepton number exchange between neutrinos and the fluid, showing the term-by-term cancellations that must occur for this consistency to be achieved.Comment: Version accepted by Phys. Rev.

Relativistic MHD with Adaptive Mesh Refinement

This paper presents a new computer code to solve the general relativistic magnetohydrodynamics (GRMHD) equations using distributed parallel adaptive mesh refinement (AMR). The fluid equations are solved using a finite difference Convex ENO method (CENO) in 3+1 dimensions, and the AMR is Berger-Oliger. Hyperbolic divergence cleaning is used to control the $\nabla\cdot {\bf B}=0$ constraint. We present results from three flat space tests, and examine the accretion of a fluid onto a Schwarzschild black hole, reproducing the Michel solution. The AMR simulations substantially improve performance while reproducing the resolution equivalent unigrid simulation results. Finally, we discuss strong scaling results for parallel unigrid and AMR runs.Comment: 24 pages, 14 figures, 3 table

Vortical and Wave Modes in 3D Rotating Stratified Flows: Random Large Scale Forcing

Utilizing an eigenfunction decomposition, we study the growth and spectra of energy in the vortical and wave modes of a 3D rotating stratified fluid as a function of $\epsilon = f/N$. Working in regimes characterized by moderate Burger numbers, i.e. $Bu = 1/\epsilon^2 < 1$ or $Bu \ge 1$, our results indicate profound change in the character of vortical and wave mode interactions with respect to $Bu = 1$. As with the reference state of $\epsilon=1$, for $\epsilon < 1$ the wave mode energy saturates quite quickly and the ensuing forward cascade continues to act as an efficient means of dissipating ageostrophic energy. Further, these saturated spectra steepen as $\epsilon$ decreases: we see a shift from $k^{-1}$ to $k^{-5/3}$ scaling for $k_f < k < k_d$ (where $k_f$ and $k_d$ are the forcing and dissipation scales, respectively). On the other hand, when $\epsilon > 1$ the wave mode energy never saturates and comes to dominate the total energy in the system. In fact, in a sense the wave modes behave in an asymmetric manner about $\epsilon = 1$. With regard to the vortical modes, for $\epsilon \le 1$, the signatures of 3D quasigeostrophy are clearly evident. Specifically, we see a $k^{-3}$ scaling for $k_f < k < k_d$ and, in accord with an inverse transfer of energy, the vortical mode energy never saturates but rather increases for all $k < k_f$. In contrast, for $\epsilon > 1$ and increasing, the vortical modes contain a progressively smaller fraction of the total energy indicating that the 3D quasigeostrophic subsystem plays an energetically smaller role in the overall dynamics.Comment: 18 pages, 6 figs. (abbreviated abstract

Unconstrained Hamiltonian formulation of General Relativity with thermo-elastic sources

A new formulation of the Hamiltonian dynamics of the gravitational field interacting with(non-dissipative) thermo-elastic matter is discussed. It is based on a gauge condition which allows us to encode the six degrees of freedom of the ``gravity + matter''-system (two gravitational and four thermo-mechanical ones), together with their conjugate momenta, in the Riemannian metric q_{ij} and its conjugate ADM momentum P^{ij}. These variables are not subject to constraints. We prove that the Hamiltonian of this system is equal to the total matter entropy. It generates uniquely the dynamics once expressed as a function of the canonical variables. Any function U obtained in this way must fulfil a system of three, first order, partial differential equations of the Hamilton-Jacobi type in the variables (q_{ij},P^{ij}). These equations are universal and do not depend upon the properties of the material: its equation of state enters only as a boundary condition. The well posedness of this problem is proved. Finally, we prove that for vanishing matter density, the value of U goes to infinity almost everywhere and remains bounded only on the vacuum constraints. Therefore the constrained, vacuum Hamiltonian (zero on constraints and infinity elsewhere) can be obtained as the limit of a ``deep potential well'' corresponding to non-vanishing matter. This unconstrained description of Hamiltonian General Relativity can be useful in numerical calculations as well as in the canonical approach to Quantum Gravity.Comment: 29 pages, TeX forma