Magnetic field are transported and tangled by turbulence, even as they lose identity due to nonideal or resistive effects. On balance field lines undergo stretch-twist-fold processes. The curvature field, a scalar that measures the tangling of the magnetic field lines, is studied in detail here, in the context of magnetohydrodynamic turbulence. A central finding is that the magnitudes of the curvature and the magnetic field are anti-correlated. High curvature co-locates with low magnetic field, which gives rise to power-law tails of the probability density function of the curvature field. The curvature drift term that converts magnetic energy into flow and thermal energy, largely depends on the curvature field behavior, a relationship that helps to explain particle acceleration due to curvature drift. This adds as well to evidence that turbulent effects most likely play an essential role in particle energization since turbulence drives stronger tangled field configurations, and therefore curvature.

In the context of space and astrophysical plasma turbulence and particle heating, several vocabularies emerge for estimating turbulent energy dissipation rate, including Kolmogorov-Yaglom third-order law and, in its various forms, $\boldsymbol{j}\cdot\boldsymbol{E}$ (work done by the electromagnetic field on particles), and $-\left( \boldsymbol{P} \cdot \nabla \right) \cdot \boldsymbol{u}$ (pressure-strain interaction), to name a couple. It is now understood that these energy transfer channels, to some extent, are correlated with coherent structures. In particular, we find that different energy dissipation proxies, although not point-wise correlated, are concentrated in proximity to each other, for which they decorrelate in a few $d_i$(s). However, the energy dissipation proxies dominate at different scales. For example, there is an inertial range over which the third-order law is meaningful. Contributions from scale bands stemming from scale-dependent spatial filtering show that, the energy exchange through $\boldsymbol{j}\cdot\boldsymbol{E}$ mainly results from large scales, while the energy conversion from fluid flow to internal through $-\left( \boldsymbol{P} \cdot \nabla \right) \cdot \boldsymbol{u}$ dominates at small scales.

We numerically study the gravitational collapse of a massless scalar field with spherical symmetry in Einstein-aether theory, and show that apparent, spin-0 and dynamical universal horizons (dUHs) can be all formed. The spacetime and the aether field are well-behaved and regular, including regions nearby these horizons (but away from the center of spherical symmetry). The spacetime outside the apparent and spin-0 horizons settles down to a static configuration, and some of such resulting static black holes were already found numerically in the literature. On the other hand, the proper distance of the outermost dUH from the apparent (or spin-0) horizon keeps increasing on $\mbox{\ae}$ther-orthogonal time slices. This indicates that the outermost dUH is evolving into the causal boundary, even for excitations with large speeds of propagation.

We present GiRaFFE, the first open-source general relativistic force-free electrodynamics (GRFFE) code for dynamical, numerical-relativity generated spacetimes. GiRaFFE adopts the strategy pioneered by McKinney and modified by Paschalidis and Shapiro to convert a GR magnetohydrodynamic (GRMHD) code into a GRFFE code. In short, GiRaFFE exists as a modification of IllinoisGRMHD, a user-friendly, open-source, dynamical-spacetime GRMHD code. Both GiRaFFE and IllinoisGRMHD leverage the Einstein Toolkit's highly-scalable infrastructure to make possible large-scale simulations of magnetized plasmas in strong, dynamical spacetimes on adaptive-mesh refinement (AMR) grids. We demonstrate that GiRaFFE passes a large suite of both flat and curved-spacetime code tests passed by a number of other state-of-the-art GRFFE codes, and is thus ready for production-scale simulations of GRFFE phenomena of key interest to relativistic astrophysics.

Kinetic plasma turbulence cascade spans multiple scales ranging from macroscopic fluid flow to sub-electron scales. Mechanisms that dissipate large scale energy, terminate the inertial range cascade and convert kinetic energy into heat are hotly debated. Here we revisit these puzzles using fully kinetic simulation. By performing scale-dependent spatial filtering on the Vlasov equation, we extract information at prescribed scales and introduce several energy transfer functions. This approach allows highly inhomogeneous energy cascade to be quantified as it proceeds down to kinetic scales. The pressure work, $-\left( \boldsymbol{P} \cdot \nabla \right) \cdot \boldsymbol{u}$, can trigger a channel of the energy conversion between fluid flow and random motions, which is a collision-free generalization of the viscous dissipation in collisional fluid. Both the energy transfer and the pressure work are strongly correlated with velocity gradients.

We report long-term simulations of black hole-neutron star binary mergers where the neutron star possesses an asymmetric magnetic field dipole. Focusing on the scenario where the neutron star is tidally disrupted by the black hole, we track the evolution of the binary up to $\approx 100$ms after the merger. We uncover more than one episode of thermally driven winds being launched along a funnel wall in all these cases beginning from $\approx 25$ms after the merger. On the other hand, we are unable to conclude presently whether the amount of ejected mass increases with the \it{degree} of asymmetry. A large-scale magnetic field configuration in the poloidal direction is formed along the funnel wall accompanied by the generation of a large Poynting flux. The magnetic field in the accretion disk around the black hole remnant is amplified by the nonaxisymmetric magnetorotational instability (MRI). The MRI growth is estimated to be in the ideal magnetohydrodynamics (MHD) regime and thus would be free from significant effects induced by potential neutrino radiation. However, the asymmetry in the magnetic field leads to increased turbulence, which causes the vertical magnetic field in the accretion disk to grow largely in a nonlinear manner.

The investigation of the diffusive transport of charged particles in a turbulent magnetic field remains a subject of considerable interest. Research has most frequently concentrated on determining the diffusion coefficient in the presence of a mean magnetic field. Here we consider diffusion of charged particles in fully three-dimensional isotropic turbulent magnetic fields with no mean field, which may be pertinent to many astrophysical situations. We identify different ranges of particle energy depending upon the ratio of the Larmor radius of the charged particle to the characteristic outer length scale of the turbulence. Two different theoretical models are proposed to calculate the diffusion coefficient, each applicable to a distinct range of particle energies. The theoretical results are compared with those from computer simulations, showing good agreement.

Evidence for inhomogeneous heating in the interplanetary plasma near current sheets dynamically generated by magnetohydrodynamic (MHD) turbulence is obtained using measurements from the ACE spacecraft. These coherent structures only constitute 19% of the data, but contribute 50% of the total plasma internal energy. Intermittent heating manifests as elevations in proton temperature near current sheets, resulting in regional heating and temperature enhancements extending over several hours. The number density of non-Gaussian structures is found to be proportional to the mean proton temperature and solar wind speed. These results suggest magnetofluid turbulence drives intermittent dissipation through a hierarchy of coherent structures, which collectively could be a significant source of coronal and solar wind heating.

Correlation anisotropy emerges dynamically in magnetohydrodynamics (MHD), producing stronger gradients across the large-scale mean magnetic field than along it. This occurs both globally and locally, and has significant implications in space and astrophysical plasmas, including particle scattering and transport, and theories of turbulence. Properties of local correlation anisotropy are further documented here by showing through numerical experiments that the effect is intensified in more localized estimates of the mean field. The mathematical formulation of this property shows that local anisotropy mixes second-order with higher order correlations. Sensitivity of local statistical estimates to higher order correlations can be understood in connection with the stochastic coordinate system inherent in such formulations. We demonstrate this in specific cases, illustrate the connection to higher order statistics by showing the sensitivity of local anisotropy to phase randomization, and thus establish that the local structure function is not a measure of the energy spectrum. Evidently the local enhancement of correlation anisotropy is of substantial fundamental interest, and this phenomenon must be understood in terms of higher order correlations, fourth-order and above.

We solve the Tolman-Oppenheimer-Volkoff equation using an equation of state (EoS) calculated in holographic QCD. The aim is to use compact astrophysical objects like neutron stars as an indicator to test holographic equations of state. We first try an EoS from a dense D4/D8/\textoverline {D8} model. In this case, however, we could not find a stable compact star, a star satisfying pressure-zero condition with a radius $R$, $p(R)=0$, within a reasonable value of the radius. This means that the EoS from the D4/D8/\textoverline {D8} model may not support any stable compact stars or may support one whose radius is very large. This might be due to a deficit of attractive force from a scalar field or two-pion exchange in the D4/D8/\textoverline {D8} model. Then, we consider D4/D6 type models with different number of quark flavors, $N_f=1,2,3$. Though the mass and radius of a holographic star is larger than those of normal neutron stars, the D4/D6 type EoS renders a stable compact star.