A New Godunov Scheme for MHD, with Application to the MRI in disks

We describe a new numerical scheme for MHD which combines a higher order Godunov method (PPM) with Constrained Transport. The results from a selection of multidimensional test problems are presented. The complete test suite used to validate the method, as well as implementations of the algorithm in both F90 and C, are available from the web. A fully three-dimensional version of the algorithm has been developed, and is being applied to a variety of astrophysical problems including the decay of supersonic MHD turbulence, the nonlinear evolution of the MHD Rayleigh-Taylor instability, and the saturation of the magnetorotational instability in the shearing box. Our new simulations of the MRI represent the first time that a higher-order Godunov scheme has been applied to this problem, providing a quantitative check on the accuracy of previous results computed with ZEUS; the latter are found to be reliable.Comment: 11 pages, style files included, Conference Proceedings: "Magnetic Fields in the Universe: from Laboratory and Stars to Primordial Structures", More information on Athena can be found at http://www.astro.princeton.edu/~jstone/athena.htm

Research in cosmic and gamma ray astrophysics

Discussed here is research in cosmic ray and gamma ray astrophysics at the Space Radiation Laboratory (SRL) of the California Institute of Technology. The primary activities discussed involve the development of new instrumentation and techniques for future space flight. In many cases these instrumentation developments were tested in balloon flight instruments designed to conduct new investigations in cosmic ray and gamma ray astrophysics. The results of these investigations are briefly summarized. Specific topics include a quantitative investigation of the solar modulation of cosmic ray protons and helium nuclei, a study of cosmic ray positron and electron spectra in interplanetary and interstellar space, the solar modulation of cosmic rays, an investigation of techniques for the measurement and interpretation of cosmic ray isotopic abundances, and a balloon measurement of the isotopic composition of galactic cosmic ray boron, carbon, and nitrogen

Propagation of a Topological Transition: the Rayleigh Instability

The Rayleigh capillary instability of a cylindrical interface between two immiscible fluids is one of the most fundamental in fluid dynamics. As Plateau observed from energetic considerations and Rayleigh clarified through hydrodynamics, such an interface is linearly unstable to fission due to surface tension. In traditional descriptions of this instability it occurs everywhere along the cylinder at once, triggered by infinitesimal perturbations. Here we explore in detail a recently conjectured alternate scenario for this instability: front propagation. Using boundary integral techniques for Stokes flow, we provide numerical evidence that the viscous Rayleigh instability can indeed spread behind a front moving at constant velocity, in some cases leading to a periodic sequence of pinching events. These basic results are in quantitative agreement with the marginal stability criterion, yet there are important qualitative differences associated with the discontinuous nature of droplet fission. A number of experiments immediately suggest themselves in light of these results.Comment: 15 pages, 7 figures, Te

Finite Nuclei in the Quark-Meson Coupling (QMC) Model

We report the first use of the effective QMC energy density functional (EDF), derived from a quark model of hadron structure, to study a broad range of ground state properties of even-even nuclei across the periodic table in the non-relativistic Hartree-Fock+BCS framework. The novelty of the QMC model is that the nuclear medium effects are treated through modification of the internal structure of the nucleon. The density dependence is microscopically derived and the spin-orbit term arises naturally. The QMC EDF depends on a single set of four adjustable parameters having clear physical basis. When applied to diverse ground state data the QMC EDF already produces, in its present simple form, overall agreement with experiment of a quality comparable to a representative Skyrme EDF. There exist however multiple Skyrme paramater sets, frequently tailored to describe selected nuclear phenomena. The QMC EDF set of fewer parameters, as derived in this work, is not open to such variation, chosen set being applied, without adjustment, to both the properties of finite nuclei and nuclear matter.Comment: 9 pages, 1 table, 4 figures; in print in Phys. Rev. Letters. A minor change in the abstract, a few typos corrected and some small technical adjustments made to comply with the journal regulation

Quasi-Chemical and Structural Analysis of Polarizable Anion Hydration

Quasi-chemical theory is utilized to analyze the roles of solute polarization and size in determining the structure and thermodynamics of bulk anion hydration for the Hofmeister series Cl$^-$, Br$^-$, and I$^-$. Excellent agreement with experiment is obtained for whole salt hydration free energies using the polarizable AMOEBA force field. The quasi-chemical approach exactly partitions the solvation free energy into inner-shell, outer-shell packing, and outer-shell long-ranged contributions by means of a hard-sphere condition. Small conditioning radii, even well inside the first maximum of the ion-water(oxygen) radial distribution function, result in Gaussian behavior for the long-ranged contribution that dominates the ion hydration free energy. The spatial partitioning allows for a mean-field treatment of the long-ranged contribution, leading to a natural division into first-order electrostatic, induction, and van der Waals terms. The induction piece exhibits the strongest ion polarizability dependence, while the larger-magnitude first-order electrostatic piece yields an opposing but weaker polarizability dependence. In addition, a structural analysis is performed to examine the solvation anisotropy around the anions. As opposed to the hydration free energies, the solvation anisotropy depends more on ion polarizability than on ion size: increased polarizability leads to increased anisotropy. The water dipole moments near the ion are similar in magnitude to bulk water, while the ion dipole moments are found to be significantly larger than those observed in quantum mechanical studies. Possible impacts of the observed over-polarization of the ions on simulated anion surface segregation are discussed.Comment: slight revision, in press at J. Chem. Phy

Energetics in MRI driven Turbulence

In these proceedings we present recent efforts to understand the energetics of magnetohydrodynamic (MHD) turbulence driven by the magneto-rotational instability (MRI). These studies are carried out in the local (shearing box) approximation using the Athena simulation code. Athena is a higher order Godunov algorithm based on the piecewise parabolic method (PPM), the corner transport upwind (CTU) integration algorithm, and the constrained transport (CT) algorithm for evolving the magnetic field. This algorithm is particularly suited for these studies owing to the conservation properties of a Godunov scheme and the particular implementation of the shearing box source terms used here. We present a variety of calculations which may be compared directly to previously published results and discuss them in some detail. The only significant discrepancy found between the results presented here and in the published literature involves the turbulent heating rate. We observe the presence of recurrent channel solutions in calculations involving a mean vertical magnetic field and the associated time lag between the energy injection and thermalization rate. We also present the results of a shearing box calculation which includes an optically thin radiative term with a cooling rate selected to match the turbulent heating rate. Some properties which we find uniformly present in all of the calculations presented here are compressible fluctuations, spiral waves and weak shocks. It is found that these compressible modes dominate the temporal fluctuations in the probability distribution functions for most of the thermodynamic variables; only the specific entropy is relatively immune to their effects.Comment: 14 pages, Conference Proceedings: "Magnetic Fields in the Universe: from Laboratory and Stars to Primordial Structures