Comparative Monte Carlo Efficiency by Monte Carlo Analysis

We propose a modified power method for computing the subdominant eigenvalue $\lambda_2$ of a matrix or continuous operator. Here we focus on defining simple Monte Carlo methods for its application. The methods presented use random walkers of mixed signs to represent the subdominant eigenfuction. Accordingly, the methods must cancel these signs properly in order to sample this eigenfunction faithfully. We present a simple procedure to solve this sign problem and then test our Monte Carlo methods by computing the $\lambda_2$ of various Markov chain transition matrices. We first computed ${\lambda_2}$ for several one and two dimensional Ising models, which have a discrete phase space, and compared the relative efficiencies of the Metropolis and heat-bath algorithms as a function of temperature and applied magnetic field. Next, we computed $\lambda_2$ for a model of an interacting gas trapped by a harmonic potential, which has a mutidimensional continuous phase space, and studied the efficiency of the Metropolis algorithm as a function of temperature and the maximum allowable step size $\Delta$. Based on the $\lambda_2$ criterion, we found for the Ising models that small lattices appear to give an adequate picture of comparative efficiency and that the heat-bath algorithm is more efficient than the Metropolis algorithm only at low temperatures where both algorithms are inefficient. For the harmonic trap problem, we found that the traditional rule-of-thumb of adjusting $\Delta$ so the Metropolis acceptance rate is around 50% range is often sub-optimal. In general, as a function of temperature or $\Delta$, $\lambda_2$ for this model displayed trends defining optimal efficiency that the acceptance ratio does not. The cases studied also suggested that Monte Carlo simulations for a continuum model are likely more efficient than those for a discretized version of the model.Comment: 23 pages, 8 figure

New techniques for experimental generation of two-dimensional blade-vortex interaction at low Reynolds numbers

An experimental investigation of two dimensional blade vortex interaction was held at NASA Langley Research Center. The first phase was a flow visualization study to document the approach process of a two dimensional vortex as it encountered a loaded blade model. To accomplish the flow visualization study, a method for generating two dimensional vortex filaments was required. The numerical study used to define a new vortex generation process and the use of this process in the flow visualization study were documented. Additionally, photographic techniques and data analysis methods used in the flow visualization study are examined

Coloured mulch as a weed control technology and yield booster for summer savory

An investigation into the effect of coloured mulch technology as a technique to control weeds when growing the essential oil plant, summer savory (Satureja hortensis) was made. As well as weed control, the effects on the production of crop biomass and essential oil content and quality were also considered. The mulch treatments produced significantly more biomass than either of the control treatments (which used no mulch either with or without herbicide). The white mulch treatment produced the greatest biomass, closely followed by the red mulch treatment. The blue mulch treatment was third in ranking, although not significantly greater than the black mulch. Estimates of the quantity of essential oil produced by each treatment followed a similar trend to that shown by biomass production

Randomised positive control trial of NSAID and antimicrobial treatment for calf fever caused by pneumonia

One hundred and fifty-four preweaning calves were followed between May and October 2015. Calves were fitted with continuous monitoring temperature probes (TempVerified FeverTag), programmed so a flashing light emitting diode (LED) light was triggered following six hours of a sustained ear canal temperature of ≥39.7°C. A total of 83 calves (61.9 per cent) developed undifferentiated fever, with a presumptive diagnosis of pneumonia through exclusion of other calf diseases. Once fever was detected, calves were randomly allocated to treatment groups. Calves in group 1 (NSAID) received 2 mg/kg flunixin meglumine (Allevinix, Merial) for three consecutive days and group 2 (antimicrobial) received 6 mg/kg gamithromycin (Zactran, Merial). If fever persisted for 72 hours after the initial treatment, calves were given further treatment (group 1 received antimicrobial and group 2 received NSAID). Calves in group 1 (NSAID) were five times more likely (P=0.002) to require a second treatment (the antimicrobial) after 72 hours to resolve the fever compared with the need to give group 2 (antimicrobial) calves a second treatment (NSAID). This demonstrates the importance of ongoing monitoring and follow-up of calves with respiratory disease. However, of calves with fever in group 1 (NSAID), 25.7 per cent showed resolution following NSAID-only treatment with no detrimental effect on the development of repeated fever or daily live weight gain. This suggests that NSAID alone may be a useful first-line treatment, provided adequate attention is given to ongoing monitoring to identify those cases that require additional antimicrobial treatment

Laser-Doppler gas-velocity instrument

Three-D instrument using a laser light source measures both turbulence and mean velocity of subsonic and supersonic gas flows. This instrument is based on the measurement of the Doppler frequency shift of light waves scattered by moving particles in the gas stream

Monte Carlo Determination of Multiple Extremal Eigenpairs

We present a Monte Carlo algorithm that allows the simultaneous determination of a few extremal eigenpairs of a very large matrix without the need to compute the inner product of two vectors or store all the components of any one vector. The new algorithm, a Monte Carlo implementation of a deterministic one we recently benchmarked, is an extension of the power method. In the implementation presented, we used a basic Monte Carlo splitting and termination method called the comb, incorporated the weight cancellation method of Arnow {\it et al.}, and exploited a new sampling method, the sewing method, that does a large state space sampling as a succession of small state space samplings. We illustrate the effectiveness of the algorithm by its determination of the two largest eigenvalues of the transfer matrices for variously-sized two-dimensional, zero field Ising models. While very likely useful for other transfer matrix problems, the algorithm is however quite general and should find application to a larger variety of problems requiring a few dominant eigenvalues of a matrix.Comment: 22 pages, no figure

Is U3Ni3Sn4 best described as near a quantum critical point?

Although most known non-Fermi liquid (NFL) materials are structurally or chemically disordered, the role of this disorder remains unclear. In particular, very few systems have been discovered that may be stoichiometric and well ordered. To test whether U3Ni3Sn4 belongs in this latter class, we present measurements of the x-ray absorption fine structure (XAFS) of polycrystalline and single-crystal U3Ni3Sn4 samples that are consistent with no measurable local structural disorder. We also present temperature-dependent specific heat data in applied magnetic fields as high as 8 T that show features that are inconsistent with the antiferromagnetic Griffiths' phase model, but do support the conclusion that a Fermi liquid/NFL crossover temperature increases with applied field. These results are inconsistent with theoretical explanations that require strong disorder effects, but do support the view that U3Ni3Sn4 is a stoichiometric, ordered material that exhibits NFL behavior, and is best described as being near an antiferromagnetic quantum critical point.Comment: 9 pages, 8 figures, in press with PR

Flexible story generation with Norms and Preferences in computer role playing games

Interactive storytelling is a strength of table-top role playing games as they are facilitated by a game master (GM) who directs the narrative and devises game scenarios. One difficulty with the implementation is the large amount of time