The onset of a small-scale turbulent dynamo at low magnetic Prandtl numbers

We study numerically the dependence of the critical magnetic Reynolds number Rmc for the turbulent small-scale dynamo on the hydrodynamic Reynolds number Re. The turbulence is statistically homogeneous, isotropic, and mirror--symmetric. We are interested in the regime of low magnetic Prandtl number Pm=Rm/Re<1, which is relevant for stellar convective zones, protostellar disks, and laboratory liquid-metal experiments. The two asymptotic possibilities are Rmc->const as Re->infinity (a small-scale dynamo exists at low Pm) or Rmc/Re=Pmc->const as Re->infinity (no small-scale dynamo exists at low Pm). Results obtained in two independent sets of simulations of MHD turbulence using grid and spectral codes are brought together and found to be in quantitative agreement. We find that at currently accessible resolutions, Rmc grows with Re with no sign of approaching a constant limit. We reach the maximum values of Rmc~500 for Re~3000. By comparing simulations with Laplacian viscosity, fourth-, sixth-, and eighth-order hyperviscosity and Smagorinsky large-eddy viscosity, we find that Rmc is not sensitive to the particular form of the viscous cutoff. This work represents a significant extension of the studies previously published by Schekochihin et al. 2004, PRL 92, 054502 and Haugen et al. 2004, PRE, 70, 016308 and the first detailed scan of the numerically accessible part of the stability curve Rmc(Re).Comment: 4 pages, emulateapj aastex, 2 figures; final version as published in ApJL (but with colour figures

Fourth Order Algorithms for Solving the Multivariable Langevin Equation and the Kramers Equation

We develop a fourth order simulation algorithm for solving the stochastic Langevin equation. The method consists of identifying solvable operators in the Fokker-Planck equation, factorizing the evolution operator for small time steps to fourth order and implementing the factorization process numerically. A key contribution of this work is to show how certain double commutators in the factorization process can be simulated in practice. The method is general, applicable to the multivariable case, and systematic, with known procedures for doing fourth order factorizations. The fourth order convergence of the resulting algorithm allowed very large time steps to be used. In simulating the Brownian dynamics of 121 Yukawa particles in two dimensions, the converged result of a first order algorithm can be obtained by using time steps 50 times as large. To further demostrate the versatility of our method, we derive two new classes of fourth order algorithms for solving the simpler Kramers equation without requiring the derivative of the force. The convergence of many fourth order algorithms for solving this equation are compared.Comment: 19 pages, 2 figure

Big Sagebrush Response to One-Way and Two-Way Chaining in Southeastern Utah

A decadent, mixed stand of Wyoming big sagebrush, Artemisia tridentata wyomingensis, and mountain big sagebrush, Artemisia tridentata vaseyana, located north of Cisco, Utah, was subjected to one-way and two-way chaining treatments in November 1987. The effect of the treatments on plant community characteristics and shrub vigor was documented over a 3-year period. Stand density was reduced 60 percent on sites chained two ways and 43 percent on sites chained over once. Shrubs on one-way chained sites produced more leader growth in 1989 and 1990 than those on untreated sites or sites chained two ways. Browse production on one-way chained sites surpassed that of untreated sites and two-way chained sites by 140 percent and 350 percent, respectively. Over the short term, a one-way chaining was shown to be an effective method for improving sagebrush vigor and production on a critical mule deer winter range

Quantum Monte Carlo Studies of Relativistic Effects in Light Nuclei

Relativistic Hamiltonians are defined as the sum of relativistic one-body kinetic energy, two- and three-body potentials and their boost corrections. In this work we use the variational Monte Carlo method to study two kinds of relativistic effects in the binding energy of 3H and 4He. The first is due to the nonlocalities in the relativistic kinetic energy and relativistic one-pion exchange potential (OPEP), and the second is from boost interaction. The OPEP contribution is reduced by about 15% by the relativistic nonlocality, which may also have significant effects on pion exchange currents. However, almost all of this reduction is canceled by changes in the kinetic energy and other interaction terms, and the total effect of the nonlocalities on the binding energy is very small. The boost interactions, on the other hand, give repulsive contributions of 0.4 (1.9) MeV in 3H (4He) and account for 37% of the phenomenological part of the three-nucleon interaction needed in the nonrelativistic Hamiltonians.Comment: 33 pages, RevTeX, 11 PostScript figures, submitted to Physical Review

Thermodynamics and Sound Speeds at the Chapman-Jouguet State

Some thermodynamic relations about an equilibrium Chapman-Jouguet (CJ) state are obtained. Relations for sound speeds in the wave velocity-particle velocity plane are derived. A relation between the slope of the sound speed in this plane and the asymptotic slope of the Hugoniot is suggested

On the construction of high-order force gradient algorithms for integration of motion in classical and quantum systems

A consequent approach is proposed to construct symplectic force-gradient algorithms of arbitrarily high orders in the time step for precise integration of motion in classical and quantum mechanics simulations. Within this approach the basic algorithms are first derived up to the eighth order by direct decompositions of exponential propagators and further collected using an advanced composition scheme to obtain the algorithms of higher orders. Contrary to the scheme by Chin and Kidwell [Phys. Rev. E 62, 8746 (2000)], where high-order algorithms are introduced by standard iterations of a force-gradient integrator of order four, the present method allows to reduce the total number of expensive force and its gradient evaluations to a minimum. At the same time, the precision of the integration increases significantly, especially with increasing the order of the generated schemes. The algorithms are tested in molecular dynamics and celestial mechanics simulations. It is shown, in particular, that the efficiency of the new fourth-order-based algorithms is better approximately in factors 5 to 1000 for orders 4 to 12, respectively. The results corresponding to sixth- and eighth-order-based composition schemes are also presented up to the sixteenth order. For orders 14 and 16, such highly precise schemes, at considerably smaller computational costs, allow to reduce unphysical deviations in the total energy up in 100 000 times with respect to those of the standard fourth-order-based iteration approach.Comment: 23 pages, 2 figures; submitted to Phys. Rev.

Cost effectiveness of adherence to IDSA/ATS guidelines in elderly patients hospitalized for Community-Aquired Pneumonia

A copy of the survey questionnaire given to the expert panel regarding patient utilities. (PDF 8.28 kb

Quantum Monte Carlo calculations of six-quark states

The variational Monte Carlo method is used to find the ground state of six quarks confined to a cavity of diameter R_c, interacting via an assumed non-relativistic constituent quark model (CQM) Hamiltonian. We use a flux-tube model augmented with one-gluon and one-pion exchange interactions, which has been successful in describing single hadron spectra. The variational wave function is written as a product of three-quark nucleon states with correlations between quarks in different nucleons. We study the role of quark exchange effects by allowing flux-tube configuration mixing. An accurate six-body variational wave function is obtained. It has only ~13% rms fluctuation in the total energy and yields a standard deviation of ~<.1%; small enough to be useful in discerning nuclear interaction effects from the large rest mass of the two nucleons. Results are presented for three values of the cavity diameter, R_c=2, 4, and 6 fm. They indicate that the flux-tube model Hamiltonian with gluon and pion exchange requires revisions in order to obtain agreement with the energies estimated from realistic two-nucleon interactions. We calculate the two-quark probability distribution functions and show how they may be used to study and adjust the model Hamiltonian.Comment: 49 pages, 13 figures, submitted to Phys. Rev.

Phenomenological Lambda-Nuclear Interactions

Variational Monte Carlo calculations for ${_{\Lambda}^4}H$ (ground and excited states) and ${_{\Lambda}^5}He$ are performed to decipher information on ${\Lambda}$-nuclear interactions. Appropriate operatorial nuclear and ${\Lambda}$-nuclear correlations have been incorporated to minimize the expectation values of the energies. We use the Argonne $\upsilon_{18}$ two-body NN along with the Urbana IX three-body NNN interactions. The study demonstrates that a large part of the splitting energy in ${_{\Lambda}^4}H$ ($0^+-1^+$) is due to the three-body ${\Lambda}$ NN forces. $_{\Lambda}^{17}O$ hypernucleus is analyzed using the {\it s}-shell results. $\Lambda$ binding to nuclear matter is calculated within the variational framework using the Fermi-Hypernetted-Chain technique. There is a need to correctly incorporate the three-body ${\Lambda}$ NN correlations for $\Lambda$ binding to nuclear matter.Comment: 18 pages (TeX), 2 figure