We study the properties of the local transverse deviations of magnetic field lines at a fixed moment in time. Those deviations "evolve" smoothly in a plane normal to the field-line direction as one moves that plane along the field line. Since the evolution can be described by a planar flow in the normal plane, we derive most of our results in the context of a toy model for planar fluid flow. We then generalize our results to include the effects of field-line curvature. We show that the type of flow is determined by the two non-zero eigenvalues of the gradient of the normalized magnetic field. The eigenvalue difference quantifies the local rate of squeezing or coiling of neighboring field lines, which we relate to standard notions of fluid vorticity and shear. The resulting squeezing rate can be used in the detection of null points, hyperbolic flux tubes and current sheets. Once integrated along field lines, that rate gives a squeeze factor, which is an approximation to the squashing factor, which is usually employed in locating quasi-separatrix layers (QSLs), which are possible sites for magnetic reconnection. Unlike the squeeze factor, the squashing factor can miss QSLs for which field lines are squeezed and then unsqueezed. In that regard, the squeeze factor is a better proxy for locating QSLs than the squashing factor. In another application of our analysis, we construct an approximation to the local rate of twist of neighboring field lines, which we refer to as the coiling rate. That rate can be integrated along a field line to give a coiling number, $\mathrm{N_c}$. We show that unlike the standard local twist number, $\mathrm{N_c}$ gives an unbiased approximation to the number of twists neighboring field lines make around one another. $\mathrm{N_c}$ can be useful for the study of flux rope instabilities, such as the kink instability, and can be used in the detection of flux ropes.

Quasi-Separatrix Layers (QSLs) are a useful proxy for the locations where current sheets can develop in the solar corona, and give valuable information about the connectivity in complicated magnetic field configurations. However, calculating QSL maps even for 2-dimensional slices through 3-dimensional models of coronal magnetic fields is a non-trivial task as it usually involves tracing out millions of magnetic field lines with immense precision. Thus, extending QSL calculations to three dimensions has rarely been done until now. In order to address this challenge, we present QSL Squasher -- a public, open-source code, which is optimized for calculating QSL maps in both two and three dimensions on GPUs. The code achieves large processing speeds for three reasons, each of which results in an order-of-magnitude speed-up. 1) The code is parallelized using OpenCL. 2) The precision requirements for the QSL calculation are drastically reduced by using perturbation theory. 3) A new boundary detection criterion between quasi-connectivity domains is used, which quickly identifies possible QSL locations which need to be finely sampled by the code. That boundary detection criterion relies on finding the locations of abrupt field-line length changes, which we do by introducing a new Field-line Length Edge (FLEDGE) map. We find FLEDGE maps useful on their own as a quick-and-dirty substitute for QSL maps. QSL Squasher allows constructing high-resolution 3D FLEDGE maps in a matter of minutes, which is two orders of magnitude faster than calculating the corresponding 3D QSL maps. We include a sample of calculations done using QSL Squasher to demonstrate its capabilities as a QSL calculator, as well as to compare QSL and FLEDGE maps.

Magnetic reconnection, the central engine that powers explosive phenomena throughout the Universe, is also perceived as one of the principal mechanisms for accelerating particles to high energies. Although various signatures of magnetic reconnection have been frequently reported, observational evidence that links particle acceleration directly to the reconnection site has been rare, especially for space plasma environments currently inaccessible to $\textit{in situ}$ measurements. Here we utilize broadband radio dynamic imaging spectroscopy available from the Karl G. Jansky Very Large Array to observe decimetric type III radio bursts in a solar jet with high angular ($\sim$20$''$), spectral ($\sim$1 %), and temporal resolution (50 milliseconds). These observations allow us to derive detailed trajectories of semi-relativistic (tens of keV) electron beams in the low solar corona with unprecedentedly high angular precision ($<0''.65$). We found that each group of electron beams, which corresponds to a cluster of type III bursts with 1-2-second duration, diverges from an extremely compact region ($\sim$600 km$^2$) in the low solar corona. The beam-diverging sites are located behind the erupting jet spire and above the closed arcades, coinciding with the presumed location of magnetic reconnection in the jet eruption picture supported by extreme ultraviolet/X-ray data and magnetic modeling. We interpret each beam-diverging site as a reconnection null point where multitudes of magnetic flux tubes join and reconnect. Our data suggest that the null points likely consist of a high level of density inhomogeneities possibly down to 10-km scales. These results, at least in the present case, strongly favor a reconnection-driven electron acceleration scenario.

We present the analysis of an unusual failed eruption captured in high cadence and in many wavelengths during the observing campaign in support of the VAULT2.0 sounding rocket launch. The refurbished Very high Angular resolution Ultraviolet Telescope (VAULT2.0) is a Ly$\alpha$ ($\lambda$ 1216 {\AA}) spectroheliograph launched on September 30, 2014. The campaign targeted active region NOAA AR 12172 and was closely coordinated with the Hinode and IRIS missions and several ground-based observatories (NSO/IBIS, SOLIS, and BBSO). A filament eruption accompanied by a low level flaring event (at the GOES C-class level) occurred around the VAULT2.0 launch. No Coronal Mass Ejection (CME) was observed. The eruption and its source region, however, were recorded by the campaign instruments in many atmospheric heights ranging from the photosphere to the corona in high cadence and spatial resolution. This is a rare occasion which enables us to perform a comprehensive investigation on a failed eruption. We find that a rising Magnetic Flux Rope-like (MFR) structure was destroyed during its interaction with the ambient magnetic field creating downflows of cool plasma and diffuse hot coronal structures reminiscent of "cusps". We employ magnetofrictional simulations to show that the magnetic topology of the ambient field is responsible for the destruction of the MFR. Our unique observations suggest that the magnetic topology of the corona is a key ingredient for a successful eruption.

The standard model for eruptive flares has in the past few years been extended to 3D. It predicts typical J-shaped photospheric footprints of the coronal current layer, forming at similar locations as the Quasi-Separatrix Layers (QSLs). Such a morphology is also found for flare ribbons observed in the EUV band, as well as in non-linear force-free field (NLFFF) magnetic field extrapolations and models. We study the evolution of the photospheric traces of the current density and flare ribbons, both obtained with the SDO instruments. We investigate the photospheric current evolution during the 6 September 2011 X-class flare (SOL2011-09-06T22:20) from observational data of the magnetic field obtained with HMI. This evolution is compared with that of the flare ribbons observed in the EUV filters of the AIA. We also compare the observed electric current density and the flare ribbon morphology with that of the QSLs computed from the flux rope insertion method/NLFFF model. The NLFFF model shows the presence of a fan-spine configuration of overlying field lines, due to the presence of a parasitic polarity, embedding an elongated flux rope that appears in the observations as two parts of a filament. The QSLs, evolved via a magnetofrictional method, also show similar morphology and evolution as both the current ribbons and the EUV flare ribbons obtained at several times during the flare. For the first time, we propose a combined analysis of the photospheric traces of an eruptive flare, in a complex topology, with direct measurements of electric currents and QSLs from observational data and a magnetic field model. The results, obtained by two different and independent approaches, 1) confirm previous results of current increase during the impulsive phase of the flare, 2) show how NLFFF models can capture the essential physical signatures of flares even in a complex magnetic field topology.

We present sCOLA -- an extension of the N-body COmoving Lagrangian Acceleration (COLA) method to the spatial domain. Similar to the original temporal-domain COLA, sCOLA is an N-body method for solving for large-scale structure in a frame that is comoving with observers following trajectories calculated in Lagrangian Perturbation Theory. Incorporating the sCOLA method in an N-body code allows one to gain computational speed by capturing the gravitational potential from the far field using perturbative techniques, while letting the N-body code solve only for the near field. The far and near fields are completely decoupled, effectively localizing gravity for the N-body side of the code. Thus, running an N-body code for a small simulation volume using sCOLA can reproduce the results of a standard N-body run for the same small volume embedded inside a much larger simulation. We demonstrate that sCOLA can be safely combined with the original temporal-domain COLA. sCOLA can be used as a method for performing zoom-in simulations. It also allows N-body codes to be made embarrassingly parallel, thus allowing for efficiently tiling a volume of interest using grid computing. Moreover, sCOLA can be useful for cheaply generating large ensembles of accurate mock halo catalogs required to study galaxy clustering. Surveys that will benefit the most are ones with large aspect ratios, such as pencil-beam surveys, where sCOLA can easily capture the effects of large-scale transverse modes without the need to substantially increase the simulated volume. As an illustration of the method, we present proof-of-concept zoom-in simulations using a freely available sCOLA-based N-body code.

Motivated by the results presented in a companion paper, here we give a simple analytical expression for the matter n-point functions in the Zel'dovich approximation (ZA) both in real and in redshift space (including the angular case). We present numerical results for the 2-dimensional redshift-space correlation function, as well as for the equilateral configuration for the real-space 3-point function. We compare those to the tree-level results. Our analysis is easily extendable to include Lagrangian bias, as well as higher-order perturbative corrections to the ZA. The results should be especially useful for modelling probes of large-scale structure in the linear regime, such as the Baryon Acoustic Oscillations. We make the numerical code used in this paper freely available.

We present a pedagogical systematic investigation of the accuracy of Eulerian and Lagrangian perturbation theories of large-scale structure. We show that significant differences exist between them especially when trying to model the Baryon Acoustic Oscillations (BAO). We find that the best available model of the BAO in real space is the Zel'dovich Approximation (ZA), giving an accuracy of <~3% at redshift of z=0 in modelling the matter 2-pt function around the acoustic peak. All corrections to the ZA around the BAO scale are perfectly perturbative in real space. Any attempt to achieve better precision requires calibrating the theory to simulations because of the need to renormalize those corrections. In contrast, theories which do not fully preserve the ZA as their solution, receive O(1) corrections around the acoustic peak in real space at z=0, and are thus of suspicious convergence at low redshift around the BAO. As an example, we find that a similar accuracy of 3% for the acoustic peak is achieved by Eulerian Standard Perturbation Theory (SPT) at linear order only at z~4. Thus even when SPT is perturbative, one needs to include loop corrections for z<~4 in real space. In Fourier space, all models perform similarly, and are controlled by the overdensity amplitude, thus recovering standard results. However, that comes at a price. Real space cleanly separates the BAO signal from non-linear dynamics. In contrast, Fourier space mixes signal from short mildly non-linear scales with the linear signal from the BAO to the level that non-linear contributions from short scales dominate. Therefore, one has little hope in constructing a systematic theory for the BAO in Fourier space.

We present the COmoving Lagrangian Acceleration (COLA) method: an N-body method for solving for Large Scale Structure (LSS) in a frame that is comoving with observers following trajectories calculated in Lagrangian Perturbation Theory (LPT). Unlike standard N-body methods, the COLA method can straightforwardly trade accuracy at small-scales in order to gain computational speed without sacrificing accuracy at large scales. This is especially useful for cheaply generating large ensembles of accurate mock halo catalogs required to study galaxy clustering and weak lensing, as those catalogs are essential for performing detailed error analysis for ongoing and future surveys of LSS. As an illustration, we ran a COLA-based N-body code on a box of size 100Mpc/h with particles of mass ~5*10^9Msolar/h. Running the code with only 10 timesteps was sufficient to obtain an accurate description of halo statistics down to halo masses of at least 10^11Msolar/h. This is only at a modest speed penalty when compared to mocks obtained with LPT. A standard detailed N-body run is orders of magnitude slower than our COLA-based code. The speed-up we obtain with COLA is due to the fact that we calculate the large-scale dynamics exactly using LPT, while letting the N-body code solve for the small scales, without requiring it to capture exactly the internal dynamics of halos. Achieving a similar level of accuracy in halo statistics without the COLA method requires at least 3 times more timesteps than when COLA is employed.

We obtain approximations for the CDM particle trajectories starting from Lagrangian Perturbation Theory. These estimates for the CDM trajectories result in approximations for the density in real and redshift space, as well as for the momentum density that are better than what standard Eulerian and Lagrangian perturbation theory give. For the real space density, we find that our proposed approximation gives a good cross-correlation (>95%) with the non-linear density down to scales almost twice smaller than the non-linear scale, and six times smaller than the corresponding scale obtained using linear theory. This allows for a speed-up of an order of magnitude or more in the scanning of the cosmological parameter space with N-body simulations for the scales relevant for the baryon acoustic oscillations. Possible future applications of our method include baryon acoustic peak reconstruction, building mock galaxy catalogs, momentum field reconstruction.

The Baryon Acoustic Oscillations (BAO) in the large-scale structure of the universe leave a distinct peak in the two-point correlation function of the matter distribution. That acoustic peak is smeared and shifted by bulk flows and non-linear evolution. However, it has been shown that it is still possible to sharpen the peak and remove its shift by undoing the effects of the bulk flows. We propose an improvement to the standard acoustic peak reconstruction. Contrary to the standard approach, the new scheme has no free parameters, treats the large-scale modes consistently, and uses optimal filters to extract the BAO information. At redshift of zero, the reconstructed linear matter power spectrum leads to a markedly improved sharpening of the reconstructed acoustic peak compared to standard reconstruction.

We present a simple physically motivated picture for the mildly non-linear regime of structure formation, which captures the effects of the bulk flows. We apply this picture to develop a method to significantly reduce the sample variance in cosmological N-body simulations at the scales relevant to the Baryon Acoustic Oscillations (BAO). The results presented in this paper will allow for a speed-up of an order of magnitude (or more) in the scanning of the cosmological parameter space using N-body simulations for studies which require a good handle of the mildly non-linear regime, such as those targeting the BAO. Using this physical picture we develop a simple formula, which allows for the rapid calculation of the mildly non-linear matter power spectrum to percent level accuracy, and for robust estimation of the BAO scale.

We present a new formalism to study large-scale structure in the universe. The result is a hierarchy (which we call the "Helmholtz Hierarchy") of equations describing the phase space statistics of cold dark matter (CDM). The hierarchy features a physical ordering parameter which interpolates between the Zel'dovich approximation and fully-fledged gravitational interactions. The results incorporate the effects of stream crossing. We show that the Helmholtz hierarchy is self-consistent and obeys causality to all orders. We present an interpretation of the hierarchy in terms of effective particle trajectories.

We approximately compute the bispectrum induced on the CMB temperature by fluctuations in the standard recombination epoch. Of all the second order sources that can induce non-Gaussianity during recombination, we concentrate on those proportional to the perturbation in the free electron density, which is about a factor of 5 larger than the other first order perturbations. This term induces some non-Gaussianity by delaying the time of recombination and by changing the photon diffusion scale. We find that the signal is not scale invariant, peaked on squeezed triangles with the smaller multipole around the scale of the first acoustic peak, and that its size corresponds to an effective f_NL ~ -3.5, which could be marginally detected by Planck if both temperature and polarization are measured.

We derive the full set of second-order equations governing the evolution of cosmological perturbations, including the effects of the first-order electron number density perturbations, \delta_e. We provide a detailed analysis of the perturbations to the recombination history of the universe and show that a perturbed version of the Peebles effective 3-level atom is sufficient for obtaining the evolution of \delta_e for comoving wavenumbers smaller than 1Mpc^{-1}. We calculate rigorously the perturbations to the Ly\alpha escape probability and show that to a good approximation it is governed by the local baryon velocity divergence. For modes shorter than the photon diffusion scale, we find that \delta_e is enhanced during recombination by a factor of roughly 5 relative to other first-order quantities sourcing the CMB anisotropies at second order. Using these results, in a companion paper we calculate the CMB bispectrum generated during recombination.

When thin accretion disks around black holes are perturbed, the main restoring force is gravity. If gas pressure, magnetic stresses, and radiation pressure are neglected, the disk remains thin as long as orbits do not intersect. Intersections would result in pressure forces which limit the growth of perturbations. We find that a discrete set of perturbations is possible for which orbits remain non-intersecting for arbitrarily long times. These modes define a discrete set of frequencies. We classify all long-lived perturbations for arbitrary potentials and show how their mode frequencies are related to pattern speeds computed from the azimuthal and epicyclic frequencies. We show that modes are concentrated near radii where the pattern speed has vanishing radial derivative. We explore these modes around Kerr black holes as a possible explanation for the high-frequency quasi-periodic oscillations of black hole binaries such as GRO J1655-40. The long-lived modes are shown to coincide with diskoseismic waves in the limit of small sound speed. While the waves have long lifetime, they have the wrong frequencies to explain the pairs of high-frequency quasi-periodic oscillations observed in black hole binaries.