The Galactic microquasar GRS 1915+105 shows rich variability which is categorized into different classes. In this paper, we report temporal and spectral analysis of GRS 1915+105 to study the properties of the accretion flow when the light curve is showing $\theta$ class variability. For this purpose, we use the LAXPC data from the Target of Opportunity observations of India's first multi-wavelength astronomy satellite AstroSat. The $\theta$ class is marked by the recurrent appearance of U-shaped regions in the light curve, where the photon count rate first decreases rapidly and then increases slowly. For our analysis, we use U-shaped regions of first two orbits (02345 and 02346) on 2016 March 04. In both of the cases, the dynamic Power Density Spectra (PDS) showed significant power at around $4-5$ Hz, suggesting the presence of a low-frequency Quasi-Periodic Oscillation (QPO) around that frequency interval. The QPO frequency is found to increase with time when the enhancement of the energy flux also takes place. From the evolution of the spectra, we determine the evolution of the accretion flow parameters in both of these observations. Fitting the spectra with the transonic flow solution based Two-Component Advective Flow (TCAF) model in $4-25$ keV energy band shows that the Keplerian disk accretion rate rises with the rise in the radiation intensity, while the location of the centrifugal pressure driven shock front decreases. In both these data, a gradual increment of power-law photon index with intensity is observed suggesting the progressive softening of the source.

Galactic transient black hole candidate (BHC) MAXI~J1535-571 was discovered on 2017 September 02 simultaneously by {\it MAXI}/GSC and {\it Swift}/BAT instruments. It has also been observed by India's first multi-wavelength astronomy-mission satellite {\it AstroSat}, during the rising phase of its 2017-18 outburst. We make both the spectral and the temporal analysis of the source during 2017 September 12-17 using data of {\it AstroSat}'s Large Area X-ray Proportional Counter (LAXPC) in the energy range of $3-40$~keV to infer the accretion flow properties of the source. Spectral analysis is done with the physical two-component advective flow (TCAF) solution-based {\it fits} file. From the nature of the variation of the TCAF model fitted physical flow parameters, we conclude and confirm that during our analysis period source was in the soft-intermediate spectral state. We observe sharp type-C quasi-periodic oscillations (QPOs) in the frequency range of $\sim 1.74-2.81$~Hz, which are not common in the soft-intermediate spectral state. For a better understanding of the nature and evolution of these type-C QPOs, a dynamic study of the power density spectra is done. We also investigate the origin of these QPOs from the shock oscillation model. We find that as is expected in the soft-Intermediate spectral states, non-satisfaction of Rankine-Hugoniot conditions for non-dissipative shocks and not their resonance oscillations may be the cause of the origin of type-C QPOs.

The class variable source GRS 1915+105 exhibits a wide range of time variabilities in time scales of a few seconds to a few days. Depending on the count rates in different energy bands and the nature of the conventional color-color diagram, the variabilities were classified into sixteen classes, which were later sequenced in ascending order of Comptonization Efficiency (CE). The CE is the ratio of the phenomenological power-law and disk blackbody model fitted spectral photon fluxes. For a better understanding of the flow dynamics in one of the variability classes, namely $\chi$, we fitted spectra with the physical two-component advective flow (TCAF) model. In $\chi$ class, X-ray flux is steady with no significant variation, however, various $\chi$ subclasses are observed at different X-ray flux and variations of count rates across different $\chi$ subclasses must be linked to the variation of flow parameters such as the accretion rates, be it the Keplerian disc rate and/or the low angular momentum halo rate. We find that in the $\chi_{2,4}$ classes, which are reportedly devoid of significant outflows, the spectra could be fitted well using TCAF solution alone. In the $\chi_{1,3}$ class, which are always linked with outflows, a cutoff power-law model is needed in addition to the TCAF solution. At the same time, the normalization required by this model along with the variation of photon index and exponential roll-off factor provides us with the information on the relative dominance of the outflow in the latter two classes. TCAF fit also provides us with the size and location of the Compton cloud along with its optical depth. Thus by fitting with TCAF, a physical understanding of the flow geometry in different $\chi$ classes of GRS 1915+105 has been obtained.

Surveys have shown radio-loud (RL) quasars constitute 10%-15% of the total quasar population and rest are radio-quiet (RQ). However, it is unknown if this radio-loud fraction (RLF) remains consistent among different parameter spaces. This study shows that RLF increases for increasing full width half maximum (FWHM) velocity of the H$\beta$ broad emission line (z < 0.75). To analyse the reason, we compared bolometric luminosity of RL and RQ quasars sample which have FWHM of H$\beta$ broad emission line greater than 15000km/s (High Broad Line or HBL) with which have FWHM of H$\beta$ emission line less than 2500km/s (Low Broad Line or LBL). From the distributions we can conclude for the HBL, RQ and RL quasars are peaking separately and RL quasars are having higher values whereas for the LBL the peaks are almost indistinguishable. We predicted selection effects could be the possible reason but to conclude anything more analysis is needed. Then we compared our result with Wills & Browne (1986) and have shown that some objects from our sample do not follow the pattern of the logR vs FWHM plot where R is the ratio of 5 GHz radio core flux density with the extended radio lobe flux density.

For the first time, we explore the tightly coupled interior-magnetosphere system of Mercury by employing a three-dimensional ten-moment multifluid model. This novel fluid model incorporates the non-ideal effects including the Hall effect, inertia, and tensorial pressures that are critical for collisionless magnetic reconnection; therefore, it is particularly well suited for investigating $collisionless$ magnetic reconnection in Mercury's magnetotail and at the planet's magnetopause. The model is able to reproduce the observed magnetic field vectors, field-aligned currents, and cross-tail current sheet asymmetry (beyond the MHD approach) and the simulation results are in good agreement with spacecraft observations. We also study the magnetospheric response of Mercury to a hypothetical extreme event with an enhanced solar wind dynamic pressure, which demonstrates the significance of induction effects resulting from the electromagnetically-coupled interior. More interestingly, plasmoids (or flux ropes) are formed in Mercury's magnetotail during the event, indicating the highly dynamic nature of Mercury's magnetosphere.

Analytic scaling relations are derived for a phenomenological model of the plasmoid instability in an evolving current sheet, including the effects of reconnection outflow. Two scenarios are considered, where the plasmoid instability can be triggered either by an injected initial perturbation or by the natural noise of the system (here referred to as the system noise). The two scenarios lead to different scaling relations because the initial noise decays when the linear growth of the plasmoid instability is not sufficiently fast to overcome the advection loss caused by the reconnection outflow, whereas the system noise represents the lowest level of fluctuations in the system. The leading order approximation for the current sheet width at disruption takes the form of a power-law multiplied by a logarithmic factor, and from that, the scaling relations for the wavenumber and the linear growth rate of the dominant mode are obtained. When the effects of the outflow are neglected, the scaling relations agree, up to the leading order approximation, with previously derived scaling relations based on a principle of least time. The analytic scaling relations are validated with numerical solutions of the model.

We use Smoothed Particle Hydrodynamics to study viscous accretion flows around a weakly magnetic neutron star. We show the formation of multiple ``boundary" layers in presence of both cooling and viscosity. We find that with the introduction of a small viscosity in a sub-Keplerian flow, much like the wind accretion in HMXBs such as Cir X-1, only a single Normal Boundary Layer (NBOL) forms to adjust the rotational velocity component. With the increase of viscosity, the region extends radially and beyond some critical value, a RAdiative KEplerian Disk/layer (RAKED) forms between the sub-Keplerian flow and the NBOL. When viscosity is increased further only NBOL and RAKED remain. In all such cases, the CENtrifugal pressure dominated BOundary Layer (CENBOL) is formed, away from the star, as in the case of black holes. This is the first self-consistent study where such a transition from sub-Keplerian flows has been reported for neutron stars. We also identify the connection between accretion and ejection of matter, following the Two-Component Advective Flow for black holes, for neutron stars. The results are crucial in the understanding of the formation of disks, boundary layers and outflows in wind dominated neutron star systems.

The galactic black hole candidate Cygnus X-1, one of the brightest sources in the sky, is the first ever black hole candidate to be discovered. Despite being a very well-studied object due to its persistent brightness in X-rays, there has been much difficulty in determining its mass since its discovery. In spite of continuous efforts since the early nineteen seventies, there is yet no concensus on its mass for nearly four decades. The uncertainties in the mass measurements are due to the high degree of error involved in the measurement of its distance. In 2011, Orosz et al. constrained the mass of this object to be M = 14.8 +/- 1.0 M_Sun using dynamical methods. In this paper, we obtained the mass of Cygnus X-1, using a completely independent method, namely, carrying out the spectral analysis using Two Component Advective Flow (TCAF) solution based fits file and the archival data of RXTE PCA instrument. Our result does not require the distance of the source or the information about the companion. Each spectral fit with the TCAF gives one best fitted mass. Averaging fitted masses of Cygnus X-1 over a span of five months of observation during its persistent hard phase, mass of the source comes out to be M_avg = 14.20 +/- 0.36 M_Sun, which is consistent with the dynamically estimated mass.

The Magnetospheric Multiscale (MMS) mission has given us unprecedented access to high cadence particle and field data of magnetic reconnection at Earth's magnetopause. MMS first passed very near an X-line on 16 October 2015, the Burch event, and has since observed multiple X-line crossings. Subsequent 3D particle-in-cell (PIC) modeling efforts of and comparison with the Burch event have revealed a host of novel physical insights concerning magnetic reconnection, turbulence induced particle mixing, and secondary instabilities. In this study, we employ the Gkeyll simulation framework to study the Burch event with different classes of extended, multi-fluid magnetohydrodynamics (MHD), including models that incorporate important kinetic effects, such as the electron pressure tensor, with physics-based closure relations designed to capture linear Landau damping. Such fluid modeling approaches are able to capture different levels of kinetic physics in global simulations and are generally less costly than fully kinetic PIC. We focus on the additional physics one can capture with increasing levels of fluid closure refinement via comparison with MMS data and existing PIC simulations.

We carry out the first robust numerical simulation of accretion flows on a weakly magnetized neutron star using Smoothed Particle Hydrodynamics (SPH). We follow the Two-Component Advective Flow (TCAF) paradigm for black holes, and focus only on the advective component for the case of a neutron star. This low viscosity sub-Keplerian flow will create a normal boundary layer (or, NBOL) right on the star surface in addition to the centrifugal pressure supported boundary layer (or, CENBOL) present in a black hole accretion. These density jumps could give rise to standing or oscillating shock fronts. During a hard spectral state, the incoming flow has a negligible viscosity causing more sub-Keplerian component as compared to the Keplerian disc component. We show that our simulation of flows with a cooling and a negligible viscosity produces precisely two shocks and strong supersonic winds from these boundary layers. We find that the specific angular momentum of matter dictates the locations and the nature of oscillations of these shocks. For low angular momentum flows, the radial oscillation appears to be preferred. For flows with higher angular momentum, the vertical oscillation appears to become dominant. In all the cases, asymmetries w.r.t. the Z=0 plane are seen and instabilities set in due to the interaction of inflow and outgoing strong winds. Our results capture both the low and high-frequency quasi-periodic oscillations without invoking magnetic fields or any precession mechanism. Most importantly, these solutions directly corroborate observed features of wind dominated high-mass X-ray binaries, such as Cir X-1.

We use data from the All Wavelength Extended Groth Strip International Survey to construct stacked X-ray maps of optically bright active galaxies (AGN) and an associated control sample of galaxies at high redshift (z less than 0.6). From our analysis of the surface brightness profiles obtained from these X-ray maps, we find evidence of feedback from the active nuclei. We find that excluding galaxies and AGN, residing in group environments, from our samples enhances the significance of our detection. Our results support the tentative findings of Chatterjee et al. who use X-ray selected AGN for their analysis. We discuss the implications of these results in the context of quantifying AGN feedback and show that the current method can be used to extract X-ray source population in high redshift galaxies.

Helioseismology has revealed the internal density and rotation profiles of the Sun. Yet, knowledge of its magnetic fields and meridional circulation is confined much closer to the surface, and latitudinal entropy gradients are below detectable limits. While numerical simulations can offer insight into the interior dynamics and help identify which ingredients are necessary to reproduce particular observations, some features of the Sun can be understood analytically from an equilibrium perspective. Examples of such features include: the 1D density profile arising from steady-state energy transport from the core to the surface, and the tilting of isorotation contours in the convection zone (CZ) due to baroclinic forcing. To help identify which features can be explained by equilibrium, we propose analyzing stationary axisymmetric ideal magnetohydrodynamic flows in the solar regime. By prescribing an appropriate entropy profile at the surface, we recover a rotation profile that reasonably matches observations in the bulk of the CZ. Additionally, by including the effects of poloidal flow, we reproduce a feature that is reminiscent of the near surface shear layer. However, no tachocline-like feature is seen in hydrodynamic equilibrium, suggesting the importance of either dynamics or magnetic fields in its description.

In this chapter, we present a brief and non-exhaustive review of the developments of theoretical models for accretion flows around neutron stars. A somewhat chronological summary of crucial observations and modelling of timing and spectral properties are given in sections 2 and 3. In section 4, we argue why and how the Two-Component Advective Flow (TCAF) solution can be applied to the cases of neutron stars when suitable modifications are made for the NSs. We showcase some of our findings from Monte Carlo and Smoothed Particle Hydrodynamic simulations which further strengthens the points raised in section 4. In summary, we remark on the possibility of future works using TCAF for both weakly magnetic and magnetic Neutron Stars.

The black hole transient H1743-322 exhibited several outbursts with temporal and spectral variability since its discovery in 1977. These outbursts occur at a quasi-regular recurrence period of around $0.5-2$ years, since its rediscovery in March 2003. We investigate accretion flow dynamics around the Low Mass X-ray Binary H1743-322 during its 2004 outburst using the RXTE/PCA archival data. We use Two Component Advective Flow (TCAF) solution to analyse the spectral data. From the fits with TCAF solution, we obtain day to day variation of physical accretion rates of Keplerian and sub-Keplerian components, size of the Compton cloud and its other properties. Analysis of the spectral properties of the 2004 outburst by keeping fitted normalization to be in a narrow range and its timing properties in terms of the presence and absence of QPOs, enable us to constrain the mass of the black hole in a range of $10.31 M_{\odot} - 14.07 M_{\odot}$ which is consistent with other estimates reported in the literature.

We explore spectral properties of a Two-Component Advective Flow (TCAF) around a neutron star. We compute the effects of thermal Comptonization of soft photons emitted from a Keplerian disc and the boundary layer of the neutron star by the post-shock region of a sub-Keplerian flow, formed due to the centrifugal barrier. The shock location $X_s$ is also the inner edge of the Keplerian disc. We compute a series of realistic spectra assuming a set of electron temperatures of the post-shock region $T_{CE}$, the temperature of the normal boundary layer (NBOL) $T_{NS}$ of the neutron star and the shock location $X_s$. These parameters depend on the disc and halo accretion rates ($\dot{m_d}$ and $\dot{m_h}$, respectively) which control the resultant spectra. We find that the spectrum becomes harder when $\dot{m}_h$ is increased. The spectrum is controlled strongly by $T_{NS}$ due to its proximity to the Comptonizing cloud since photons emitted from the NBOL cool down the post-shock region very effectively. We also show the evidence of spectral hardening as the inclination angle of the disc is increased.

The plasmoid instability in evolving current sheets has been widely studied due to its effects on the disruption of current sheets, the formation of plasmoids, and the resultant fast magnetic reconnection. In this Letter, we study the role of the plasmoid instability in two-dimensional magnetohydrodynamic (MHD) turbulence by means of high-resolution direct numerical simulations. At sufficiently large magnetic Reynolds number ($R_m=10^6$), the combined effects of dynamic alignment and turbulent intermittency lead to a copious formation of plasmoids in a multitude of intense current sheets. The disruption of current sheet structures facilitates the energy cascade towards small scales, leading to the breaking and steepening of the energy spectrum. In the plasmoid-mediated regime, the energy spectrum displays a scaling that is close to the spectral index $-2.2$ as proposed by recent analytic theories. We also demonstrate that the scale-dependent dynamic alignment exists in 2D MHD turbulence and the corresponding slope of the alignment angle is close to 0.25.

We present the first results from a reverberation-mapping campaign undertaken during the first half of 2012, with additional data on one AGN (NGC 3227) from a 2014 campaign. Our main goals are (1) to determine the black hole masses from continuum-Hbeta reverberation signatures, and (2) to look for velocity-dependent time delays that might be indicators of the gross kinematics of the broad-line region. We successfully measure Hbeta time delays and black hole masses for five AGNs, four of which have previous reverberation mass measurements. The values measured here are in agreement with earlier estimates, though there is some intrinsic scatter beyond the formal measurement errors. We observe velocity dependent Hbeta lags in each case, and find that the patterns have changed in the intervening five years for three AGNs that were also observed in 2007.

Magnetohydrodynamic turbulence and magnetic reconnection are ubiquitous in astrophysical environments. In most situations, these processes do not occur in isolation, but interact with each other. This renders a comprehensive theory of these processes highly challenging. Here, we propose a theory of magnetohydrodynamic turbulence driven at large scale that self-consistently accounts for the mutual interplay with magnetic reconnection occurring at smaller scales. Magnetic reconnection produces plasmoids that grow from turbulence-generated noise and eventually disrupt the sheet-like structures in which they are born. The disruption of these structures leads to a modification of the turbulent energy cascade, which, in turn, exerts a feedback effect on the plasmoid formation via the turbulence-generated noise. The energy spectrum in this plasmoid-mediated range steepens relative to the standard inertial range and does not follow a simple power law. As a result of the complex interplay between turbulence and reconnection, we also find that the length scale which marks the beginning of the plasmoid-mediated range and the dissipation length scale do not obey true power laws. The transitional magnetic Reynolds number above which the plasmoid formation becomes statistically significant enough to affect the turbulent cascade is fairly modest, implying that plasmoids are expected to modify the turbulent path to dissipation in many astrophysical systems.

We present optical continuum lags for two Seyfert 1 galaxies, MCG+08-11-011 and NGC 2617, using monitoring data from a reverberation mapping campaign carried out in 2014. Our light curves span the ugriz filters over four months, with median cadences of 1.0 and 0.6 days for MCG+08-11-011 and NGC\,2617, respectively, combined with roughly daily X-ray and near-UV data from Swift for NGC 2617. We find lags consistent with geometrically thin accretion-disk models that predict a lag-wavelength relation of $\tau \propto \lambda^{4/3}$. However, the observed lags are larger than predictions based on standard thin-disk theory by factors of 3.3 for MCG+08-11-011 and 2.3 for NGC\,2617. These differences can be explained if the mass accretion rates are larger than inferred from the optical luminosity by a factor of 4.3 in MCG+08-11-011 and a factor of 1.3 in NGC\,2617, although uncertainty in the SMBH masses determines the significance of this result. While the X-ray variability in NGC\,2617 precedes the UV/optical variability, the long 2.6 day lag is problematic for coronal reprocessing models.

We revisit Parker's conjecture of current singularity formation in 3D line-tied plasmas using a recently developed numerical method, variational integration for ideal magnetohydrodynamics in Lagrangian labeling. With the frozen-in equation built-in, the method is free of artificial reconnection, and hence it is arguably an optimal tool for studying current singularity formation. Using this method, the formation of current singularity has previously been confirmed in the Hahm--Kulsrud--Taylor problem in 2D. In this paper, we extend this problem to 3D line-tied geometry. The linear solution, which is singular in 2D, is found to be smooth for arbitrary system length. However, with finite amplitude, the linear solution can become pathological when the system is sufficiently long. The nonlinear solutions turn out to be smooth for short systems. Nonetheless, the scaling of peak current density versus system length suggests that the nonlinear solution may become singular at finite length. With the results in hand, we can neither confirm nor rule out this possibility conclusively, since we cannot obtain solutions with system length near the extrapolated critical value.

We present the first results from an optical reverberation mapping campaign executed in 2014, targeting the active galactic nuclei (AGN) MCG+08-11-011, NGC 2617, NGC 4051, 3C 382, and Mrk 374. Our targets have diverse and interesting observational properties, including a "changing look" AGN and a broad-line radio galaxy. Based on continuum-H$\beta$ lags, we measure black hole masses for all five targets. We also obtain H$\gamma$ and He{\sc ii}\,$\lambda 4686$ lags for all objects except 3C 382. The He{\sc ii}\,$\lambda 4686$ lags indicate radial stratification of the BLR, and the masses derived from different emission lines are in general agreement. The relative responsivities of these lines are also in qualitative agreement with photoionization models. These spectra have extremely high signal-to-noise ratios (100--300 per pixel) and there are excellent prospects for obtaining velocity-resolved reverberation signatures.

The scaling of plasmoid instability maximum linear growth rate with respect to Lundquist number $S$ in a Sweet-Parker current sheet, $\gamma_{max}\sim S^{1/4}$, indicates that at high $S$, the current sheet will break apart before it approaches the Sweet-Parker width. Therefore, a proper description for the onset of the plasmoid instability must incorporate the evolving process of the current sheet. We carry out a series of two-dimensional simulations and develop diagnostics to separate fluctuations from an evolving background. It is found that the fluctuation amplitude starts to grow only when the linear growth rate is sufficiently large ($\gamma_{max}\tau_{A}>O(1)$) to overcome convective losses. The linear growth rate continues to rise until the sizes of plasmoids become comparable to the inner layer width of the tearing mode. At this point the current sheet is disrupted and the instability enters the early nonlinear regime. The growth rate suddenly decreases, but the fluctuation amplitude continues to grow until it reaches nonlinear saturation. We identify important time scales of the instability development, as well as scalings for the linear growth rate, current sheet width, and dominant wavenumber at current sheet disruption. These scalings depend on not only the Lundquist number, but also the initial noise amplitude. A phenomenological model that reproduces scalings from simulation results is proposed. The model incorporates the effect of reconnection outflow, which is crucial for yielding a critical Lundquist number $S_{c}$ below which disruption does not occur. The critical Lundquist number $S_{c}$ is not a constant value but has a weak dependence on the noise amplitude.

We present a three-species multi-fluid MHD model (H$^+$, H$_2$O$^+$ and e$^-$), endowed with the requisite atmospheric chemistry, that is capable of accurately quantifying the magnitude of water ion losses from exoplanets. We apply this model to a water world with Earth-like parameters orbiting a Sun-like star for three cases: (i) current normal solar wind conditions, (ii) ancient normal solar wind conditions, and (iii) one extreme "Carrington-type" space weather event. We demonstrate that the ion escape rate for (ii), with a value of 6.0$\times$10$^{26}$ s$^{-1}$, is about an order of magnitude higher than the corresponding value of 6.7$\times$10$^{25}$ s$^{-1}$ for (i). Studies of ion losses induced by space weather events, where the ion escape rates can reach $\sim$ 10$^{28}$ s$^{-1}$, are crucial for understanding how an active, early solar-type star (e.g., with frequent coronal mass ejections) could have accelerated the depletion of the exoplanet's atmosphere. We briefly explore the ramifications arising from the loss of water ions, especially for planets orbiting M-dwarfs where such effects are likely to be significant.

The plasmoid instability has revolutionized our understanding of magnetic reconnection in astrophysical environments. By preventing the formation of highly elongated reconnection layers, it is crucial in enabling the rapid energy conversion rates that are characteristic of many astrophysical phenomena. Most of the previous studies have focused on Sweet-Parker current sheets, which, however, are unattainable in typical astrophysical systems. Here, we derive a general set of scaling laws for the plasmoid instability in resistive and visco-resistive current sheets that evolve over time. Our method relies on a principle of least time that enables us to determine the properties of the reconnecting current sheet (aspect ratio and elapsed time) and the plasmoid instability (growth rate, wavenumber, inner layer width) at the end of the linear phase. After this phase the reconnecting current sheet is disrupted and fast reconnection can occur. The scaling laws of the plasmoid instability are \emph{not} simple power laws, and depend on the Lundquist number ($S$), the magnetic Prandtl number ($P_m$), the noise of the system ($\psi_0$), the characteristic rate of current sheet evolution ($1/\tau$), as well as the thinning process. We also demonstrate that previous scalings are inapplicable to the vast majority of the astrophysical systems. We explore the implications of the new scaling relations in astrophysical systems such as the solar corona and the interstellar medium. In both these systems, we show that our scaling laws yield values for the growth rate, wavenumber, and aspect ratio that are much smaller than the Sweet-Parker based scalings.

A new formulation of the plasma resistivity that arises from the collisional momentum-transfer rate between electrons and ions is presented. The resistivity computed herein is shown to depend not only on the temperature and density but also on all other polynomial velocity-space moments of the distribution function, such as the pressure tensor and heat flux vector. The exact expression for the collisional momentum-transfer rate is determined, and is used to formulate the nonlinear anisotropic resistivity. The new formalism recovers the Spitzer resistivity, as well as the concept of thermal force if the heat flux is assumed to be proportional to a temperature gradient. Furthermore, if the pressure tensor is related to viscous stress, the latter enters the expression for the resistivity. The relative importance of the nonlinear term(s) with respect to the well-established electron inertia and Hall terms is also examined. The subtle implications of the nonlinear resistivity, and its dependence on the fluid variables, are discussed in the context of magnetized plasma environments and phenomena such as magnetic reconnection.

A general theory of the onset and development of the plasmoid instability is formulated by means of a principle of least time. The scaling relations for the final aspect ratio, transition time to rapid onset, growth rate, and number of plasmoids are derived, and shown to depend on the initial perturbation amplitude $\left({\hat w}_0\right)$, the characteristic rate of current sheet evolution $\left(1/\tau\right)$, and the Lundquist number $\left(S\right)$. They are not simple power laws, and are proportional to $S^{\alpha} \tau^{\beta} \left[\ln f(S,\tau,{\hat w}_0)\right]^\sigma$. The detailed dynamics of the instability is also elucidated, and shown to comprise of a period of quiescence followed by sudden growth over a short time scale.

Numerical simulations have consistently shown that the reconnection rate in certain collisionless regimes can be fast, on the order of $0.1 v_A B_u$, where $v_A$ and $B_u$ are the Alfv\'en speed and the reconnecting magnetic field upstream of the ion diffusion region. This particular value has been reported in myriad numerical simulations under disparate conditions. However, despite decades of research, the reasons underpinning this specific value remain mysterious. Here, we present an overview of this problem and discuss the conditions under which the "0.1 value" is attained. Furthermore, we explain why this problem should be interpreted in terms of the ion diffusion region length.

The role of the Hall term on large scale dynamo action is investigated by means of the First Order Smoothing Approximation. It is shown that the standard $\alpha$ coefficient is altered, and is zero when a specific double Beltrami state is attained, in contrast to the Alfv\'enic state for MHD dynamos. The $\beta$ coefficient is no longer positive definite, and thereby enables dynamo action even if $\alpha$-quenching were to operate. The similarities and differences with the (magnetic) shear-current effect are pointed out, and a mechanism that may be potentially responsible for $\beta < 0$ is advanced. The results are compared against previous studies, and their astrophysical relevance is also highlighted.

Although the Shakura-Sunyaev $\alpha$ viscosity prescription has been highly successful in characterizing myriad astrophysical environments, it has proven to be partly inadequate in modelling turbulent stresses driven by the MRI. Hence, we adopt the approach employed by \citet{GIO03}, but in the context of Hall magnetohydrodynamics (MHD), to study MRI turbulence. We utilize the exact evolution equations for the stresses, and the non-linear terms are closed through the invocation of dimensional analysis and physical considerations. We demonstrate that the inclusion of the Hall term leads to non-trivial results, including the modification of the Reynolds and Maxwell stresses, as well as the (asymptotic) non-equipartition between the kinetic and magnetic energies; the latter issue is also addressed via the analysis of non-linear waves. The asymptotic ratio of the kinetic and magnetic energies is shown to be \emph{independent} of the choice of initial conditions, but it is governed by the \emph{Hall parameter}. We contrast our model with the Kazantsev prescription from small-scale dynamo theory, and the Hall term does not contribute in the latter approach, illustrating the limitations of the Kazantsev formalism. We indicate potential astrophysical applications of our model, including the solar wind where a lack of equipartition has been observed.

A novel large-scale dynamo mechanism, the magnetic shear-current effect, is discussed and explored. The effect relies on the interaction of magnetic fluctuations with a mean shear flow, meaning the saturated state of the small-scale dynamo can drive a large-scale dynamo -- in some sense the inverse of dynamo quenching. The dynamo is nonhelical, with the mean-field $\alpha$ coefficient zero, and is caused by the interaction between an off-diagonal component of the turbulent resistivity and the stretching of the large-scale field by shear flow. Following up on previous numerical and analytic work, this paper presents further details of the numerical evidence for the effect, as well as an heuristic description of how magnetic fluctuations can interact with shear flow to produce the required electromotive force. The pressure response of the fluid is fundamental to this mechanism, which helps explain why the magnetic effect is stronger than its kinematic cousin, and the basic idea is related to the well-known lack of turbulent resistivity quenching by magnetic fluctuations. As well as being interesting for its applications to general high Reynolds number astrophysical turbulence, where strong small-scale magnetic fluctuations are expected to be prevalent, the magnetic shear-current effect is a likely candidate for large-scale dynamo in the unstratified regions of ionized accretion disks. Evidence for this is discussed, as well as future research directions and the challenges involved with understanding details of the effect in astrophysically relevant regimes.

Recently a variational integrator for ideal magnetohydrodynamics in Lagrangian labeling has been developed. Its built-in frozen-in equation makes it optimal for studying current sheet formation. We use this scheme to study the Hahm-Kulsrud-Taylor problem, which considers the response of a 2D plasma magnetized by a sheared field under sinusoidal boundary forcing. We obtain an equilibrium solution that preserves the magnetic topology of the initial field exactly, with a fluid mapping that is non-differentiable. Unlike previous studies that examine the current density output, we identify a singular current sheet from the fluid mapping. These results are benchmarked with a constrained Grad-Shafranov solver. The same signature of current singularity can be found in other cases with more complex magnetic topologies.

Our understanding of the process of fast reconnection has undergone a dramatic change in the last 10 years driven, in part, by the availability of high-resolution numerical simulations that have consistently demonstrated the break-up of current sheets into magnetic islands, with reconnection rates that become independent of Lundquist number, challenging the belief that fast magnetic reconnection in flares proceeds via the Petschek mechanism that invokes pairs of slow-mode shocks connected to a compact diffusion region. The reconnection sites are too small to be resolved with images but these reconnection mechanisms, Petschek and the plasmoid instability, have reconnection sites with very different density and velocity structures and so can be distinguished by high-resolution line-profiles observations. Using IRIS spectroscopic observations we obtain a survey of typical line profiles produced by small-scale events thought to be reconnection sites on the Sun. Slit-jaw images are used to investigate the plasma heating and re-configuration at the sites. A sample of 15 events from two active regions is presented. The line profiles are complex with bright cores and broad wings extending to over 300 km/s. The profiles can be reproduced with the multiple magnetic islands and acceleration sites that characterise the plasmoid instability but not by bi-directional jets that characterise the Petschek mechanism. This result suggests that if these small-scale events are reconnection sites, then fast reconnection proceeds via the plasmoid instability, rather than the Petschek mechanism during small-scale reconnection on the Sun.

We propose a new mechanism for turbulent mean-field dynamo in which the magnetic fluctuations resulting from a small-scale dynamo drive the generation of large-scale magnetic fields. This is in stark contrast to the common idea that small-scale magnetic fields should be harmful to large-scale dynamo action. These dynamos occur in the presence of large-scale velocity shear and do not require net helicity, resulting from off-diagonal components of the turbulent resistivity tensor as the magnetic analogue of the "shear-current" effect. Given the inevitable existence of non-helical small-scale magnetic fields in turbulent plasmas, as well as the generic nature of velocity shear, the suggested mechanism may help to explain generation of large-scale magnetic fields across a wide range of astrophysical objects.

This article presents a calculation of the mean electromotive force arising from general small-scale magnetohydrodynamical turbulence, within the framework of the second-order correlation approximation. With the goal of improving understanding of the accretion disk dynamo, effects arising through small-scale magnetic fluctuations, velocity gradients, density and turbulence stratification, and rotation, are included. The primary result, which supplements numerical findings, is that an off-diagonal turbulent resistivity due to magnetic fluctuations can produce large-scale dynamo action -- the magnetic analogue of the "shear-current" effect. In addition, consideration of $\alpha$ effects in the stratified regions of disks gives the puzzling result that there is no strong prediction for a sign of $\alpha$, since the effects due to kinetic and magnetic fluctuations, as well as those due to shear and rotation, are each of opposing signs and tend to cancel each other.

Nonhelical shear dynamos are studied with a particular focus on the possibility of coherent dynamo action. The primary results -- serving as a follow up to the results of Squire & Bhattacharjee [arXiv:1506.04109 (2015)] -- pertain to the "magnetic shear-current effect" as a viable mechanism to drive large-scale magnetic field generation. This effect raises the interesting possibility that the saturated state of the small-scale dynamo could drive large-scale dynamo action, and is likely to be important in the unstratified regions of accretion disk turbulence. In this paper, the effect is studied at low Reynolds numbers, removing the complications of small-scale dynamo excitation and aiding analysis by enabling the use of quasi-linear statistical simulation methods. In addition to the magnetically driven dynamo, new results on the kinematic nonhelical shear dynamo are presented. These illustrate the relationship between coherent and incoherent driving in such dynamos, demonstrating the importance of rotation in determining the relative dominance of each mechanism.

Recent result of BICEP2, revealing a larger value of tensor to scalar ratio (r), has opened up new investigations of the in ationary models to fit the experimental data. The experiment needs to reconfirm the results, specifically the consistency between Planck and BICEP2. On the other hand, the combined analysis of Planck and BICEP2 B, including the dust polarization uncertainty, brings down the upper limit on r. In this note, we reexamine the low energy effective theory of Chromo Natural In ation model and its generalization in view of such observational data. We find that the parameter space of the model admits a large value of r as well as other cosmological observables consistent with data.

Turbulence and dynamo induced by the magnetorotational instability (MRI) are analyzed using quasi-linear statistical simulation methods. It is found that homogenous turbulence is unstable to a large scale dynamo instability, which saturates to an inhomogenous equilibrium with a strong dependence on the magnetic Prandtl number (Pm). Despite its enormously reduced nonlinearity, the dependence of the angular momentum transport on Pm in the quasi-linear model is qualitatively similar to that of nonlinear MRI turbulence. This indicates that recent convergence problems may be related to large scale dynamo and suggests how dramatically simplified models may be used to gain insight into the astrophysically relevant regimes of very low or high Pm.

Supra-arcade fans are bright, irregular regions of emission that develop during eruptive flares, above flare arcades. The underlying flare arcades are thought to be a consequence of magnetic reconnection along a current sheet in the corona. At the same time, theory predicts plasma jets from the reconnection site which would be extremely difficult to observe directly because of their low density. It has been suggested that the dark supra-arcade downflows (SADs) seen falling through supra-arcade fans may be low density jet plasma. The head of a low density jet directed towards higher density plasma would be Rayleigh-Taylor unstable, and lead to the development of rapidly growing low and high density fingers along the interface. Using SDO/AIA 131A images, we show details of SADs seen from three different orientations with respect to the flare arcade and current sheet, and highlight features that have been previously unexplained, such as the splitting of SADs at their heads, but are a natural consequence of instabilities above the arcade. Comparison with 3-D magnetohydrodynamic simulations suggests that supra-arcade downflows are the result of secondary instabilities of the Rayleigh-Taylor type in the exhaust of reconnection jets.

Magnetic fields without a direction of continuous symmetry have the generic feature that neighboring field lines exponentiate away from each other and become stochastic, hence the ideal constraint of preserving magnetic field line connectivity becomes exponentially sensitive to small deviations from ideal Ohm's law. The idea of breaking field line connectivity by stochasticity as a mechanism for fast reconnection is tested with numerical simulations based on reduced magnetohydrodynamics equations with a strong guide field line-tied to two perfectly conducting end plates. Starting from an ideally stable force-free equilibrium, the system is allowed to undergo resistive relaxation. Two distinct phases are found in the process of resistive relaxation. During the quasi-static phase, rapid change of field line connectivity and strong induced flow are found in regions of high field line exponentiation. However, although the field line connectivity of individual field lines can change rapidly, the overall pattern of field line mapping appears to deform gradually. From this perspective, field line exponentiation appears to cause enhanced diffusion rather than reconnection. In some cases, resistive quasi-static evolution can cause the ideally stable initial equilibrium to cross a stability threshold, leading to formation of intense current filaments and rapid change of field line mapping into a qualitatively different pattern. It is in this onset phase that the change of field line connectivity is more appropriately designated as magnetic reconnection. Our results show that rapid change of field line connectivity appears to be a necessary, but not a sufficient condition for fast reconnection.

We study the magnetorotational instability (MRI) using nonmodal stability techniques. Despite the spectral instability of many forms of the MRI, this proves to be a natural method of analysis that is well-suited to deal with the non-self-adjoint nature of the linear MRI equations. We find that the fastest growing linear MRI structures on both local and global domains can look very different to the eigenmodes, invariably resembling waves shearing with the background flow (shear waves). In addition, such structures can grow many times faster than the least stable eigenmode over long time periods, and be localized in a completely different region of space. These ideas lead -- for both axisymmetric and non-axisymmetric modes -- to a natural connection between the global MRI and the local shearing box approximation. By illustrating that the fastest growing global structure is well described by the ordinary differential equations (ODEs) governing a single shear wave, we find that the shearing box is a very sensible approximation for the linear MRI, contrary to many previous claims. Since the shear wave ODEs are most naturally understood using nonmodal analysis techniques, we conclude by analyzing local MRI growth over finite time-scales using these methods. The strong growth over a wide range of wave-numbers suggests that nonmodal linear physics could be of fundamental importance in MRI turbulence.

We analyze the linear growth of the magnetorotational instability (MRI) in the short time limit using nonmodal methods. Our findings are quite different to standard results, illustrating that shearing wave energy can grow at the maximum MRI rate, $-d\Omega/d \ln r,$ for any choice of azimuthal and vertical wavelengths. In addition, by comparing the growth of shearing waves with static structures, we show that over short time-scales shearing waves will always be dynamically more important than static structures in the ideal limit. By demonstrating that fast linear growth is possible at all wavelengths, these results suggest that nonmodal linear physics could play a fundamental role in MRI turbulence.

We analyze orbital solutions for 48 massive multiple-star systems in the Cygnus OB2 Association, 23 of which are newly presented here, to find that the observed distribution of orbital periods is approximately uniform in log P for P<45 d, but it is not scale-free. Inflections in the cumulative distribution near 6 d, 14, d, and 45 d, suggest key physical scales of ~0.2, ~0.4, and ~1 A.U. where yet-to-be-identified phenomena create distinct features. No single power law provides a statistically compelling prescription, but if features are ignored, a power law with exponent beta = -0.22 provides a crude approximation over P=1.4 -- 2000 d, as does a piece-wise linear function with a break near 45 d. The cumulative period distribution flattens at P > 45 d, even after correction for completeness, indicating either a lower binary fraction or a shift toward low-mass companions. A high degree of similarity (91% likelihood) between the Cyg OB2 period distribution and that of other surveys suggests that the binary properties at P<25 d are determined by local physics of disk/clump fragmentation and are relatively insensitive to environmental and evolutionary factors. Fully 30% of the unbiased parent sample is a binary with period P < 45 d. Completeness corrections imply a binary fraction near 55% for P < 5000 d. The observed distribution of mass ratios 0.2<q<1 is consistent with uniform, while the observed distribution of eccentricities 0.1<e<0.6 is consistent with uniform plus an excess of e ~ 0 systems. We identify six stars, all supergiants, that exhibit aperiodic velocity variations of ~30 km/s attributed to atmospheric fluctuations.

Supra-arcade downflows (hereafter referred to as SADs) are low-emission, elongated, finger-like features usually observed in active-region coronae above post-eruption flare arcades. Observations exhibit downward moving SADs intertwined with bright upward moving spikes. Whereas SADs are dark voids, spikes are brighter, denser structures. Although SADs have been observed for decades, the mechanism of formation of SADs remains an open issue. In our three-dimensional resistive magnetohydrodynamic simulations, we demonstrate that secondary Rayleigh-Taylor type instabilities develop in the downstream region of a reconnecting current sheet. The instability results in the formation of low-density coherent structures that resemble SADs, and high-density structures that appear to be spike-like. Comparison between the simulation results and observations suggests that secondary Rayleigh-Taylor type instabilities in the exhaust of reconnecting current sheets provide a plausible mechanism for observed SADs and spikes.

The position of a mobile D-3 brane moving towards a stack of localized D-5 branes has been studied as a candidate driving inflation in the warm inflationary scenario. We compare the results obtained by considering the dissipation parameter $\Gamma$ as an arbitrary function of only the inflaton field and a particular form derived by Bastro-Gil et al \cite{Berera8}. We find that the observables remain well within the recent observational constraint for a wide range of model parameters for the first case whereas the spectral index in the later case is always predicted blue, other cosmological observables remaining well within bound for a wider range of parameters though. We also discuss the non-gaussianity generated during inflation in this model.

The constraint imposed by magnetic helicity conservation on the alpha effect is considered for both magnetically and flow dominated self-organizing plasmas. Direct numerical simulations are presented for a dominant contribution to the alpha effect, which can be cast in the functional form of a total divergence of an averaged helicity flux, called the helicity-flux-driven alpha ( H$\alpha$) effect. Direct numerical simulations of the H$\alpha$ effect are prese nted for two examples---the magnetically dominated toroidal plasma unstable to tearing modes, and the flow-dominated accretion disk.

The inflationary models developed in presence of a background radiation can be a solution to the reheating problem faced by common cold (isentropic) inflationary scenario. A D-brane system comprising of k Neuvo-Schwarz(NS) 5-branes with a transverse circle and BPS D3-branes with world volume parallel to the NS 5-branes, placed at a point on the transverse circle diametrically opposite to the NS 5-branes has a point of unstable equilibrium and the D3 brane has a geometric tachyonic mode associated with displacement of the brane along the circle. Cold inflationary scenario has been studied in connection with this geometric tachyon \cite{PA1} where it was found that one needs a background of minimum $10^4$ branes to realize a viable inflationary model. In this piece of work, we have tried to study a model of inflation driven by this geometric tachoyn in presence of radiation. We have found that compared to the isentropic scenario, to satisfy the observational bounds, the number of background branes required in this case reduces drastically and a viable model can be obtained with even six to seven NS 5-branes in the background.In this context, we have also analyzed the non-gaussianity associated with the model and observed that the concerned parameter lies well within the observation limit.

Magnetic reconnection mediated by the hyper-resistive plasmoid instability is studied with both linear analysis and nonlinear simulations. The linear growth rate is found to scale as $S_{H}^{1/6}$ with respect to the hyper-resistive Lundquist number $S_{H}\equiv L^{3}V_{A}/\eta_{H}$, where $L$ is the system size, $V_{A}$ is the Alfv\'en velocity, and $\eta_{H}$ is the hyper-resistivity. In the nonlinear regime, reconnection rate becomes nearly independent of $S_{H}$, the number of plasmoids scales as $S_{H}^{1/2}$, and the secondary current sheet length and width both scale as $S_{H}^{-1/2}$. These scalings are consistent with a heuristic argument assuming secondary current sheets are close to marginal stability. The distribution of plasmoids as a function of the enclosed flux $\psi$ is found to obey a $\psi^{-1}$ power law over an extended range, followed by a rapid fall off for large plasmoids. These results are compared with those from resistive magnetohydrodynamic studies.

The Hahm--Kulsrud (HK) [T. S. Hahm and R. M. Kulsrud, Phys. Fluids {\bf 28}, 2412 (1985)] solutions for a magnetically sheared plasma slab driven by a resonant periodic boundary perturbation illustrate fully shielded (current sheet) and fully reconnected (magnetic island) responses. On the global scale, reconnection involves solving a magnetohydrodynamic (MHD) equilibrium problem. In systems with a continuous symmetry such MHD equilibria are typically found by solving the Grad--Shafranov equation, and in slab geometry the elliptic operator in this equation is the 2-D Laplacian. Thus, assuming appropriate pressure and poloidal current profiles, a conformal mapping method can be used to transform one solution into another with different boundary conditions, giving a continuous sequence of solutions in the form of partially reconnected magnetic islands (plasmoids) separated by Syrovatsky current sheets. The two HK solutions appear as special cases.

Our understanding of magnetic reconnection in resistive magnetohydrodynamics has gone through a fundamental change in recent years. The conventional wisdom is that magnetic reconnection mediated by resistivity is slow in laminar high Lundquist ($S$) plasmas, constrained by the scaling of the reconnection rate predicted by Sweet-Parker theory. However, recent studies have shown that when $S$ exceeds a critical value $\sim10^{4}$, the Sweet-Parker current sheet is unstable to a super-Alfv\'enic plasmoid instability, with a linear growth rate that scales as $S^{1/4}$. In the fully developed statistical steady state of two-dimensional resistive magnetohydrodynamic simulations, the normalized average reconnection rate is approximately 0.01, nearly independent of $S$, and the distribution function $f(\psi)$ of plasmoid magnetic flux $\psi$ follows a power law $f(\psi)\sim\psi^{-1}$. When Hall effects are included, the plasmoid instability may trigger onset of Hall reconnection even when the conventional criterion for onset is not satisfied. The rich variety of possible reconnection dynamics is organized in the framework of a phase diagram.

In a recent numerical study [Ng et al., Astrophys. J. {\bf 747}, 109, 2012], with a three-dimensional model of coronal heating using reduced magnetohydrodynamics (RMHD), we have obtained scaling results of heating rate versus Lundquist number based on a series of runs in which random photospheric motions are imposed for hundreds to thousands of \al time in order to obtain converged statistical values. The heating rate found in these simulations saturate to a level that is independent of the Lundquist number. This scaling result was also supported by an analysis with the assumption of the Sweet-Parker scaling of the current sheets, as well as how the width, length and number of current sheets scale with Lundquist number. In order to test these assumptions, we have implemented an automated routine to analyze thousands of current sheets in these simulations and return statistical scalings for these quantities. It is found that the Sweet-Parker scaling is justified. However, some discrepancies are also found and require further study.

The distribution function $f(\psi)$ of magnetic flux $\psi$ in plasmoids formed in high-Lundquist-number current sheets is studied by means of an analytic phenomenological model and direct numerical simulations. The distribution function is shown to follow a power law $f(\psi)\sim\psi^{-1}$, which differs from other recent theoretical predictions. Physical explanations are given for the discrepant predictions of other theoretical models.

Parker's model is one of the most discussed mechanisms for coronal heating and has generated much debate. We have recently obtained new scaling results in a two-dimensional (2D) version of this problem suggesting that the heating rate becomes independent of resistivity in a statistical steady state [Ng and Bhattacharjee, Astrophys. J., 675, 899 (2008)]. Our numerical work has now been extended to 3D by means of large-scale numerical simulations. Random photospheric footpoint motion is applied for a time much longer than the correlation time of the motion to obtain converged average coronal heating rates. Simulations are done for different values of the Lundquist number to determine scaling. In the high-Lundquist number limit, the coronal heating rate obtained so far is consistent with a trend that is independent of the Lundquist number, as predicted by previous analysis as well as 2D simulations. In the same limit the average magnetic energy built up by the random footpoint motion tends to have a much weaker dependence on the Lundquist number than that in the 2D simulations, due to the formation of strong current layers and subsequent disruption when the equilibrium becomes unstable. We will present scaling analysis showing that when the dissipation time is comparable or larger than the correlation time of the random footpoint motion, the heating rate tends to become independent of Lundquist number, and that the magnetic energy production is also reduced significantly.

Mid-infrared images from the Spitzer Space Telescope Galactic Legacy Infrared MidPlane Survey Extraordinaire program reveal that the infrared source IRAS 03063+5735 is a bowshock nebula produced by an early B star, 2MASS 03101044+5747035. We present new optical spectra of this star, classify it as a B1.5 V, and determine a probable association with a molecular cloud complex at V_LSR=-38 -- -42 km/s in the outer Galaxy near l=140.59 degr, b=-0.250 degr. On the basis of spectroscopic parallax, we estimate a distance of 4.0 +/-1 kpc to both the bowshock nebula and the molecular complex. One plausible scenario is that this a high-velocity runaway star impinging upon a molecular cloud. We identify the HII region and stellar cluster associated with IRAS 03064+5638 at a projected distance of 64 pc as one plausible birth site. The spectrophotometric distance and linkage to a molecular feature provides another piece of data helping to secure the ill-determined rotation curve in the outer Galaxy. As a by-product of spectral typing this star, we present empirical spectral diagnostic diagrams suitable for approximate spectral classification of O and B stars using He lines in the little-used yellow-red portion of the optical spectrum.

We have recently carried out a computational campaign to investigate a model of coronal heating in three-dimensions using reduced magnetohydrodynamics (RMHD). Our code is built on a conventional scheme using the pseudo-spectral method, and is parallelized using MPI. The current investigation requires very long time integrations using high Lundquist numbers, where the formation of very fine current layers challenge the resolutions achievable even on massively parallel machines. We present here results of a port to Nvidia CUDA (Compute Unified Device Architecture) for hardware acceleration using graphics processing units (GPUs). In addition to a brief discussion of our general strategy, we will report code performance on several machines which span a variety of hardware configurations and capabilities. These include a desktop workstation with commodity hardware, a dedicated research workstation equipped with four Nvidia C2050 GPUs, as well as several large-scale GPU accelerated distributed memory machines: Lincoln/NCSA, Dirac/NERSC, and Keeneland/NICS.

An analytical and numerical treatment is given of a constrained version of the tectonics model developed by Priest, Heyvaerts, & Title [2002]. We begin with an initial uniform magnetic field ${\bf B} = B_0 \hat{\bf z}$ that is line-tied at the surfaces $z = 0$ and $z = L$. This initial configuration is twisted by photospheric footpoint motion that is assumed to depend on only one coordinate ($x$) transverse to the initial magnetic field. The geometric constraints imposed by our assumption precludes the occurrence of reconnection and secondary instabilities, but enables us to follow for long times the dissipation of energy due to the effects of resistivity and viscosity. In this limit, we demonstrate that when the coherence time of random photospheric footpoint motion is much smaller by several orders of magnitude compared with the resistive diffusion time, the heating due to Ohmic and viscous dissipation becomes independent of the resistivity of the plasma. Furthermore, we obtain scaling relations that suggest that even if reconnection and/or secondary instabilities were to limit the build-up of magnetic energy in such a model, the overall heating rate will still be independent of the resistivity.

The solar wind and the interstellar medium are permeated by large-scale magnetic fields that render magnetohydrodynamic (MHD) turbulence anisotropic. In the weak-turbulence limit in which three-wave interactions dominate, analytical and high-resolution numerical results based on random scattering of shear-Alfv\'en waves propagating parallel to a large-scale magnetic field, as well as direct simulations demonstrate rigorously an anisotropic energy spectrum that scales as $k^{-2}_\perp$, instead of the famous Iroshnikov-Kraichnan (IK) spectrum of $k^{-3/2}$ for the isotropic case. Even in the absence of a background magnetic field, anisotropy is found to develop with respect to the local magnetic field, although the energy spectrum is globally isotropic and is found to be consistent with a $k^{-3/2}$ scaling. It is also found in direct numerical simulations that the energy cascade rate is much closer to IK scaling than a Kolmogorov scaling. Recent observations in the solar wind on cascade rates (as functions of the proton temperature and solar wind speed at 1 AU) seem to support this result [Vasquez et al. 2007].

The role of a super-Alfv\'enic plasmoid instability in the onset of fast reconnection is studied by means of the largest Hall magnetohydrodynamics simulations to date, with system sizes up to $10^{4}$ ion skin depths ($d_{i}$). It is demonstrated that the plasmoid instability can facilitate the onset of rapid Hall reconnection, in a regime where the onset would otherwise be inaccessible because the Sweet-Parker width is significantly above $d_{i}$. However, the topology of Hall reconnection is not inevitably a single stable X-point. There exists an intermediate regime where the single X-point topology itself exhibits instability, causing the system to alternate between a single X-point geometry and an extended current sheet with multiple X-points produced by the plasmoid instability. Through a series of simulations with various system sizes relative to $d_{i}$, it is shown that system size affects the accessibility of the intermediate regime. The larger the system size is, the easier it is to realize the intermediate regime. Although our Hall MHD model lacks many important physical effects included in fully kinetic models, the fact that a single X-point geometry is not inevitable raises the interesting possibility for the first time that Hall MHD simulations may have the potential to realize reconnection with geometrical features similar to those seen in fully kinetic simulations, namely, extended current sheets and plasmoid formation.

Whether the phenomenology governing MHD turbulence is Kolmogorov or Iroshnikov-Kraichnan (IK) remains an open question, theoretically as well as observationally. The ion heating profile observed in the solar wind provides a quantitative, if indirect, observational constraint on the relevant phenomenology. Recently, a solar wind heating model based on Kolmogorov spectral scaling has produced reasonably good agreement with observations, provided the effect of turbulence generation due to pickup ions is included in the model. Without including the pickup ion contributions, the Kolmogorov scaling predicts a proton temperature profile that decays too rapidly beyond a radial distance of 15 AU. In the present study, we alter the heating model by applying an energy cascade rate based on IK scaling, and show that the model yields higher proton temperatures, within the range of observations, with or without the inclusion of the effect due to pickup ions. Furthermore, the turbulence correlation length based on IK scaling seems to follow the trend of observations better.

In recent years, a phenomenological solar wind heating model based on a turbulent energy cascade prescribed by the Kolmogorov theory has produced reasonably good agreement with observations on proton temperatures out to distances around 70 AU, provided the effect of turbulence generation due to pickup ions is included in the model. In a recent study [Ng et al., J. Geophys. Res., 115, A02101 (2010)], we have incorporated in the heating model the energy cascade rate based on Iroshnikov-Kraichnan (IK) scaling. We showed that the IK cascade rate can also produce good agreement with observations, with or without the inclusion of pickup ions. This effect was confirmed both by integrating the model using average boundary conditions at 1 AU, and by applying a method [Smith et al., Astrophys. J., 638, 508 (2006)] that uses directly observed values as boundary conditions. The effects due to pickup ions is found to be less important for the IK spectrum, which is shallower than the Kolmogorov spectrum. In this paper, we will present calculations of the pickup ions effect in more details, and discuss the physical reason why a shallower spectrum generates less waves and turbulence.

Understanding the physical processes responsible for accelerating the solar wind requires detailed measurements of the collisionless plasma in the extended solar corona. Some key clues about these processes have come from instruments that combine the power of an ultraviolet (UV) spectrometer with an occulted telescope. This combination enables measurements of ion emission lines far from the bright solar disk, where most of the solar wind acceleration occurs. Although the UVCS instrument on SOHO made several key discoveries, many questions remain unanswered because its capabilities were limited. This white paper summarizes these past achievements and also describes what can be accomplished with next-generation instrumentation of this kind.

An overview of some recent progress on magnetohydrodynamic stability and current sheet formation in a line-tied system is given. Key results on the linear stability of the ideal internal kink mode and resistive tearing mode are summarized. For nonlinear problems, a counterexample to the recent demonstration of current sheet formation by Low \emph{et al}. [B. C. Low and \AA. M. Janse, Astrophys. J. \textbf{696}, 821 (2009)] is presented, and the governing equations for quasi-static evolution of a boundary driven, line-tied magnetic field are derived. Some open questions and possible strategies to resolve them are discussed.

The Sweet-Parker layer in a system that exceeds a critical value of the Lundquist number ($S$) is unstable to the plasmoid instability. In this paper, a numerical scaling study has been done with an island coalescing system driven by a low level of random noise. In the early stage, a primary Sweet-Parker layer forms between the two coalescing islands. The primary Sweet-Parker layer breaks into multiple plasmoids and even thinner current sheets through multiple levels of cascading if the Lundquist number is greater than a critical value $S_{c}\simeq4\times10^{4}$. As a result of the plasmoid instability, the system realizes a fast nonlinear reconnection rate that is nearly independent of $S$, and is only weakly dependent on the level of noise. The number of plasmoids in the linear regime is found to scales as $S^{3/8}$, as predicted by an earlier asymptotic analysis (Loureiro \emph{et al.}, Phys. Plasmas \textbf{14}, 100703 (2007)). In the nonlinear regime, the number of plasmoids follows a steeper scaling, and is proportional to $S$. The thickness and length of current sheets are found to scale as $S^{-1}$, and the local current densities of current sheets scale as $S^{-1}$. Heuristic arguments are given in support of theses scaling relations.

(Abridged) An anisotropic theory of MHD turbulence with nonvanishing cross-helicity is constructed based on Boldyrev's alignment hypothesis and probabilities p and q for fluctuations v and b to be positively or negatively aligned. Guided by observations suggesting that the normalized cross-helicity and the probability p are approximately constant in the inertial range, a generalization of Boldyrev's theory is derived that depends on the three ratios w+/w-, epsilon+/epsilon-, and p/q. The theory reduces to Boldyrev's original theory when w+ = w-, epsilon+ = epsilon-, and p = q.

The recent demonstration of current singularity formation by Low et al. assumes that potential fields will remain potential under simple expansion or compression (Low 2006, 2007; Janse & Low 2009). An explicit counterexample to their key assumption is constructed. Our findings suggest that their results may need to be reconsidered.

(Abridged abstract) A theory of incompressible MHD turbulence recently developed by Boldyrev predicts the existence of a scale dependent angle of alignment between velocity and magnetic field fluctuations that is proportional to the lengthscale of the fluctuations to the power 1/4. In this study, plasma and magnetic field data from the Wind spacecraft are used to investigate the angle between velocity and magnetic field fluctuations in the solar wind as a function of the timescale of the fluctuations and to look for the power law scaling predicted by Boldyrev.

A recently developed spectral-element adaptive refinement incompressible magnetohydrodynamic (MHD) code [Rosenberg, Fournier, Fischer, Pouquet, J. Comp. Phys. 215, 59-80 (2006)] is applied to simulate the problem of MHD island coalescence instability (MICI) in two dimensions. MICI is a fundamental MHD process that can produce sharp current layers and subsequent reconnection and heating in a high-Lundquist number plasma such as the solar corona [Ng and Bhattacharjee, Phys. Plasmas, 5, 4028 (1998)]. Due to the formation of thin current layers, it is highly desirable to use adaptively or statically refined grids to resolve them, and to maintain accuracy at the same time. The output of the spectral-element static adaptive refinement simulations are compared with simulations using a finite difference method on the same refinement grids, and both methods are compared to pseudo-spectral simulations with uniform grids as baselines. It is shown that with the statically refined grids roughly scaling linearly with effective resolution, spectral element runs can maintain accuracy significantly higher than that of the finite difference runs, in some cases achieving close to full spectral accuracy.

The nature of current sheet formation in the vicinity of three-dimensional (3D) magnetic null points is investigated. The particular focus is upon the effect of the compressibility of the plasma on the qualitative and quantitative properties of the current sheet. An initially potential 3D null is subjected to shearing perturbations, as in a previous paper [Pontin et al., Phys. Plasmas, in press (2007)]. It is found that as the incompressible limit is approached, the collapse of the null point is suppressed, and an approximately planar current sheet aligned to the fan plane is present instead. This is the case regardless of whether the spine or fan of the null is sheared. Both the peak current and peak reconnection rate are reduced. The results have a bearing on previous analytical solutions for steady-state reconnection in incompressible plasmas, implying that fan current sheet solutions are dynamically accessible, while spine current sheet solutions are not.

The nature of the evolution of the magnetic field, and of current sheet formation, at three-dimensional (3D) magnetic null points is investigated. A kinematic example is presented which demonstrates that for certain evolutions of a 3D null (specifically those for which the ratios of the null point eigenvalues are time-dependent) there is no possible choice of boundary conditions which renders the evolution of the field at the null ideal. Resistive MHD simulations are described which demonstrate that such evolutions are generic. A 3D null is subjected to boundary driving by shearing motions, and it is shown that a current sheet localised at the null is formed. The qualitative and quantitative properties of the current sheet are discussed. Accompanying the sheet development is the growth of a localised parallel electric field, one of the signatures of magnetic reconnection. Finally, the relevance of the results to a recent theory of turbulent reconnection is discussed.