The magnetohydrodynamics (MHD) of protoplanetary disks are strongly subject to the non-ideal MHD effects arising from the low ionization fraction of the disk gas. A strong electric field induced by gas motions can heat ionized gas particles and can thereby affect the ionization balance in the disks. Our previous studies revealed that in dusty protoplanetary disks, the Ohmic conductivity decreases with increasing electric field strength until the electrical breakdown of the disk gas occurs. In this study, we extend our previous work to more general cases where both electric and magnetic fields affect the motion of plasma particles, allowing us to study the impacts of plasma heating on all non-ideal MHD effects: Ohmic, Hall, and ambipolar diffusion. We find that the upper limit on the electric current we previously derived applies even in the presence of magnetic fields. Although the Hall and ambipolar resistivities can either increase or decrease with electric field strength depending on the abundance of charged dust grains, the Ohmic resistivity always increases with electric field strength. An order-of-magnitude estimate suggests that a large-scale electric current generated by gas motions in the inner part of protoplanetary disks could exceed the upper limit. This implies that MHD motions of the inner disk, such as the motion driven by the Hall-shear instability, could either get suppressed or trigger electrical breakdown (lightning discharge). This may have important implications for gas accretion and chondrule formation in the inner part of protoplanetary disks.

The gas temperature in protoplanetary disks (PPDs) is determined by a combination of irradiation heating and accretion heating, with the latter conventionally attributed to turbulent dissipation. However, recent studies have suggested that the inner disk (a few AU) is largely laminar, with accretion primarily driven by magnetized disk winds, as a result of nonideal magnetohydrodynamic (MHD) effects from weakly ionized gas, suggesting an alternative heating mechanism by Joule dissipation. We perform local stratified MHD simulations including all three nonideal MHD effects (ohmic, Hall, and ambipolar diffusion) and investigate the role of Joule heating and the resulting disk vertical temperature profiles. We find that in the inner disk, as ohmic and ambipolar diffusion strongly suppress electrical current around the midplane, Joule heating primarily occurs at several scale heights above the midplane, making the midplane temperature much lower than that with the conventional viscous heating model. Including the Hall effect, Joule heating is enhanced/reduced when the magnetic fields threading the disks are aligned/anti-aligned with the disk rotation, but it is overall ineffective. Our results further suggest that the midplane temperature in the inner PPDs is almost entirely determined by irradiation heating, unless viscous heating can trigger thermal ionization in the disk innermost region to self-sustain magnetorotational instability turbulence.

Magnetorotational instability (MRI) has a potential to generate the vigorous turbulence in protoplanetary disks, although its turbulence strength and accretion stress remains debatable because of the uncertainty of MRI with low ionization fraction. We focus on the heating of electrons by strong electric fields which amplifies nonideal magnetohydrodynamic effects. The heated electrons frequently collide with and stick to dust grains, which in turn decreases the ionization fraction and is expected to weaken the turbulent motion driven by MRI. In order to quantitatively investigate the nonlinear evolution of MRI including the electron heating, we perform magnetohydrodynamical simulation with the unstratified shearing box. We introduce a simple analytic resistivity model depending on the current density by mimicking resistivity given by the calculation of ionization. Our simulation confirms that the electron heating suppresses magnetic turbulence when the electron heating occurs with low current density. We find a clear correlation between magnetic stress and its current density, which means that the magnetic stress is proportional to the squared current density. When the turbulent motion is completely suppressed, laminar accretion flow is caused by ordered magnetic field. We give an analytical description of the laminar state by using a solution of linear perturbation equations with resistivity. We also propose a formula that successfully predicts the accretion stress in the presence of the electron heating.

Refractory organic compounds formed in molecular clouds are among the building blocks of the solar system objects and could be the precursors of organic matter found in primitive meteorites and cometary materials. However, little is known about the evolutionary pathways of molecular cloud organics from dense molecular clouds to planetary systems. In this study, we focus on the evolution of the morphological and viscoelastic properties of molecular cloud refractory organic matter. We found that the organic residue, experimentally synthesized at about 10 K from UV-irradiated H2O-CH3OH-NH3 ice, changed significantly in terms of its nanometer- to micrometer-scale morphology and viscoelastic properties after UV irradiation at room temperature. The dose of this irradiation was equivalent to that experienced after short residence in diffuse clouds (equal or less than 10,000 years) or irradiation in outer protoplanetary disks. The irradiated organic residues became highly porous and more rigid and formed amorphous nanospherules. These nanospherules are morphologically similar to organic nanoglobules observed in the least-altered chondrites, chondritic porous interplanetary dust particles, and cometary samples, suggesting that irradiation of refractory organics could be a possible formation pathway for such nanoglobules. The storage modulus (elasticity) of photo-irradiated organic residues is about 100 MPa irrespective of vibrational frequency, a value that is lower than the storage moduli of minerals and ice. Dust grains coated with such irradiated organics would therefore stick together efficiently, but growth to larger grains might be suppressed due to an increase in aggregate brittleness caused by the strong connections between grains.

The magnetorotational instability (MRI) drives vigorous turbulence in a region of protoplanetary disks where the ionization fraction is sufficiently high. It has recently been shown that the electric field induced by the MRI can heat up electrons and thereby affect the ionization balance in the gas. In particular, in a disk where abundant dust grains are present, the electron heating causes a reduction of the electron abundance, thereby preventing further growth of the MRI. By using the nonlinear Ohm's law that takes into account electron heating, we investigate where in protoplanetary disks this negative feedback between the MRI and ionization chemistry becomes important. We find that the "e-heating zone," the region where the electron heating limits the saturation of the MRI, extends out up to 80 AU in the minimum-mass solar nebula with abundant submicron-sized grains. This region is considerably larger than the conventional dead zone whose radial extent is $\sim20$ AU in the same disk model. Scaling arguments show that the MRI turbulence in the e-heating zone should have a significantly lower saturation level. Submicron-sized grains in the e-heating zone are so negatively charged that their collisional growth is unlikely to occur. Our present model neglects ambipolar and Hall diffusion, but our estimate shows that ambipolar diffusion would also affect the MRI in the e-heating zone.

We analyze the temporal variation of the diurnal anisotropy of sub-TeV cosmic ray intensity observed with the Matsushiro (Japan) underground muon detector over two full solar activity cycles in 1985-2008. The average sidereal amplitude over the entire period is 0.034+-0.003 %, which is roughly one third of the amplitude reported from AS and deep-underground muon experiments monitoring multi-TeV GCR intensity suggesting a significant attenuation of the anisotropy due to the solar modulation. We find, on the other hand, only weak correlations either with the solar activity- or magnetic-cycles. We examine the temporal variation of the "single-band valley depth" (SBVD) quoted by the Milagro experiment and, by contrast with recent Milagro reports, we find no steady increase in the Matsushiro observations in a 7-year period between 2000 and 2007. We suggest, therefore, that the steady increase of the SBVD reported by the Milagro experiment is not caused by the decreasing solar modulation in the declining phase of the 23rd solar activity cycle.

We analyze the sidereal daily variations observed between 1985 and 2006 at Matsushiro, Japan (MAT) and between 1993 and 2005 at Liapootah, Tasmania (LPT). These stations comprise the two hemisphere network (THN) of underground muon detectors in Japan and Australia. Yearly mean harmonic vectors at MAT and LPT are more or less stable without any significant change in phase and amplitude in correlation with either the solar activity- or magnetic-cycles. In this paper, therefore, we analyze the average anisotropy over the entire observation periods, i.e. 1985-2006 for MAT and 1993-2005 for LPT. We apply to the THN data a best-fitting analysis based on a model anisotropy in space identical to that adopted by Amenomori et al. (2007) for Tibet III data. The median energies of primary cosmic rays recorded are ~0.5 TeV for THN and ~5 TeV for the Tibet III experiment. It is shown that the intensity distribution of the best-fit anisotropy is quite similar to that derived from Tibet III data, regardless of the order of magnitude difference in energies of primary particles. This, together with the THN observations, confirms that the analysis by Amenomori et al. (2007) based on the Tibet III experiment in the northern hemisphere is not seriously biased. The best-fit amplitudes of the anisotropy, on the other hand, are only one third or less of those reported by the Tibet III experiment, indicating attenuation due to solar modulation. The rigidity dependence of the anisotropy amplitude in the sub-TeV region is consistent with the spectrum reported by Hall et al. (1999), smoothly extending to the Tibet III result in the multi-TeV region. The amplitude at higher energies appears almost constant or gradually decreasing with increasing rigidity.