2,524 research outputs found
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 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, , 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 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 , even
at the highest velocities considered here (ms).Comment: Submitted to Phys. Rev.
- …