Dual Geometric Worm Algorithm for Two-Dimensional Discrete Classical Lattice Models

We present a dual geometrical worm algorithm for two-dimensional Ising models. The existence of such dual algorithms was first pointed out by Prokof'ev and Svistunov \cite{ProkofevClassical}. The algorithm is defined on the dual lattice and is formulated in terms of bond-variables and can therefore be generalized to other two-dimensional models that can be formulated in terms of bond-variables. We also discuss two related algorithms formulated on the direct lattice, applicable in any dimension. These latter algorithms turn out to be less efficient but of considerable intrinsic interest. We show how such algorithms quite generally can be "directed" by minimizing the probability for the worms to erase themselves. Explicit proofs of detailed balance are given for all the algorithms. In terms of computational efficiency the dual geometrical worm algorithm is comparable to well known cluster algorithms such as the Swendsen-Wang and Wolff algorithms, however, it is quite different in structure and allows for a very simple and efficient implementation. The dual algorithm also allows for a very elegant way of calculating the domain wall free energy.Comment: 12 pages, 6 figures, Revtex

Dissipative Transport of a Bose-Einstein Condensate

We investigate the effects of impurities, either correlated disorder or a single Gaussian defect, on the collective dipole motion of a Bose-Einstein condensate of $^7$Li in an optical trap. We find that this motion is damped at a rate dependent on the impurity strength, condensate center-of-mass velocity, and interatomic interactions. Damping in the Thomas-Fermi regime depends universally on the disordered potential strength scaled to the condensate chemical potential and the condensate velocity scaled to the peak speed of sound. The damping rate is comparatively small in the weakly interacting regime, and the damping in this case is accompanied by strong condensate fragmentation. \textit{In situ} and time-of-flight images of the atomic cloud provide evidence that this fragmentation is driven by dark soliton formation.Comment: 14 pages, 20 figure

Photoemission Spectroscopy and the Unusually Robust One Dimensional Physics of Lithium Purple Bronze

Temperature dependent photoemission spectroscopy in Li0.9Mo6O17 contributes to evidence for one dimensional physics that is unusually robust. Three generic characteristics of the Luttinger liquid are observed, power law behavior of the k-integrated spectral function down to temperatures just above the superconducting transition, k-resolved lineshapes that show holon and spinon features, and quantum critical (QC) scaling in the lineshapes. Departures of the lineshapes and the scaling from expectations in the Tomonaga Luttinger model can be partially described by a phenomenological momentum broadening that is presented and discussed. The possibility that some form of 1d physics obtains even down to the superconducting transition temperature is assessed.Comment: submitted to JPCM, Special issue article "Physics in one dimension

Application of SWAT Hydrologic Model for TMDL Development on Chapel Branch Creek Watershed, SC

2008 S.C. Water Resources Conference - Addressing Water Challenges Facing the State and Regio

Bayesian model averaging with fixed and flexible priors: theory, concepts, and calibration experiments for rainfall-runoff modeling

This paper introduces for the first time the concept of Bayesian Model Averaging (BMA) with multiple prior structures, for rainfall‐runoff modeling applications. The original BMA model proposed by Raftery et al. (2005) assumes that the prior probability density function (pdf) is adequately described by a mixture of Gamma and Gaussian distributions. Here we discuss the advantages of using BMA with fixed and flexible prior distributions. Uniform, Binomial, Binomial‐Beta, Benchmark, and Global Empirical Bayes priors along with Informative Prior Inclusion and Combined Prior Probabilities were applied to calibrate daily streamflow records of a coastal plain watershed in the South‐East USA. Various specifications for Zellner's g prior including Hyper, Fixed, and Empirical Bayes Local (EBL) g priors were also employed to account for the sensitivity of BMA and derive the conditional pdf of each constituent ensemble member. These priors were examined using the simulation results of conceptual and semi‐distributed rainfall‐runoff models. The hydrologic simulations were first coupled with a new sensitivity analysis model and a parameter uncertainty algorithm to assess the sensitivity and uncertainty associated with each model. BMA was then used to subsequently combine the simulations of the posterior pdf of each constituent hydrological model. Analysis suggests that a BMA based on combined fixed and flexible priors provides a coherent mechanism and promising results for calculating a weighted posterior probability compared to individual model calibration. Furthermore, the probability of Uniform and Informative Prior Inclusion priors received significantly lower predictive error whereas more uncertainty resulted from a fixed g prior (i.e. EBL)

Strain Relaxation Mechanisms and Local Structural Changes in Si_{1-x}$Ge_{x} Alloys

In this work, we address issues pertinent to the understanding of the structural and electronic properties of Si_{1-x} Ge_{x}alloys, namely, (i) how does the lattice constant mismatch between bulk Si and bulk Ge manifests itself in the alloy system? and (ii) what are the relevant strain release mechanisms? To provide answers to these questions, we have carried out an in-depth study of the changes in the local geometric and electronic structures arising from the strain relaxation in Si_{1-x} Ge_{x} alloys using an ab initio molecular dynamics scheme. The optimized lattice constant, while exhibiting a general trend of linear dependence on the composition (Vegard's law), shows a negative deviation from Vegard's law in the vicinity of x=0.5. We delineate the mechanisms responsible for each one of the above features. We show that the radial-strain relaxation through bond stretching is responsible for the overall trend of linear dependence of the lattice constant on the composition. On the other hand, the negative deviation from Vegard's law is shown to arise from the angular-strain relaxation.Comment: 21 pages, 7 figure

Creating Digital Coastal Watersheds: The Remote Data Acquisition Network at Bannockburn Plantation, Georgetown County, SC

2008 S.C. Water Resources Conference - Addressing Water Challenges Facing the State and Regio

A DSEL for Studying and Explaining Causation

We present a domain-specific embedded language (DSEL) in Haskell that supports the philosophical study and practical explanation of causation. The language provides constructs for modeling situations comprised of events and functions for reliably determining the complex causal relationships that emerge between these events. It enables the creation of visual explanations of these causal relationships and a means to systematically generate alternative, related scenarios, along with corresponding outcomes and causes. The DSEL is based on neuron diagrams, a visual notation that is well established in practice and has been successfully employed for causation explanation and research. In addition to its immediate applicability by users of neuron diagrams, the DSEL is extensible, allowing causation experts to extend the notation to introduce special-purpose causation constructs. The DSEL also extends the notation of neuron diagrams to operate over non-boolean values, improving its expressiveness and offering new possibilities for causation research and its applications.Comment: In Proceedings DSL 2011, arXiv:1109.032

Nonequilibrium molecular dynamics simulation of rapid directional solidification

We present the results of non-equilibrium molecular dynamics simulations for the growth of a solid binary alloy from its liquid phase. The regime of high pulling velocities, $V$, for which there is a progressive transition from solute segregation to solute trapping, is considered. In the segregation regime, we recover the exponential form of the concentration profile within the liquid phase. Solute trapping is shown to settle in progressively as $V$ is increased and our results are in good agreement with the theoretical predictions of Aziz [J. Appl. Phys. {\bf 53}, 1158 (1981)]. In addition, the fluid advection velocity is shown to remain directly proportional to $V$, even at the highest velocities considered here ($V\simeq10$ms$^{-1}$).Comment: Submitted to Phys. Rev.