The min-conflicts heuristic: Experimental and theoretical results

This paper describes a simple heuristic method for solving large-scale constraint satisfaction and scheduling problems. Given an initial assignment for the variables in a problem, the method operates by searching through the space of possible repairs. The search is guided by an ordering heuristic, the min-conflicts heuristic, that attempts to minimize the number of constraint violations after each step. We demonstrate empirically that the method performs orders of magnitude better than traditional backtracking techniques on certain standard problems. For example, the one million queens problem can be solved rapidly using our approach. We also describe practical scheduling applications where the method has been successfully applied. A theoretical analysis is presented to explain why the method works so well on certain types of problems and to predict when it is likely to be most effective

An extended abstract: A heuristic repair method for constraint-satisfaction and scheduling problems

The work described in this paper was inspired by a surprisingly effective neural network developed for scheduling astronomical observations on the Hubble Space Telescope. Our heuristic constraint satisfaction problem (CSP) method was distilled from an analysis of the network. In the process of carrying out the analysis, we discovered that the effectiveness of the network has little to do with its connectionist implementation. Furthermore, the ideas employed in the network can be implemented very efficiently within a symbolic CSP framework. The symbolic implementation is extremely simple. It also has the advantage that several different search strategies can be employed, although we have found that hill-climbing methods are particularly well-suited for the applications that we have investigated. We begin the paper with a brief review of the neural network. Following this, we describe our symbolic method for heuristic repair

Weighted-density approximation for general nonuniform fluid mixtures

In order to construct a general density-functional theory for nonuniform fluid mixtures, we propose an extension to multicomponent systems of the weighted-density approximation (WDA) of Curtin and Ashcroft [Phys. Rev. A 32, 2909 (1985)]. This extension corrects a deficiency in a similar extension proposed earlier by Denton and Ashcroft [Phys. Rev. A 42, 7312 (1990)], in that that functional cannot be applied to the multi-component nonuniform fluid systems with spatially varying composition, such as solid-fluid interfaces. As a test of the accuracy of our new functional, we apply it to the calculation of the freezing phase diagram of a binary hard-sphere fluid, and compare the results to simulation and the Denton-Ashcroft extension.Comment: 4 pages, 4 figures, to appear in Phys. Rev. E as Brief Repor

Particle Dark Energy

We explore the physics of a gas of particles interacting with a condensate that spontaneously breaks Lorentz invariance. The equation of state of this gas varies from 1/3 to less than -1 and can lead to the observed cosmic acceleration. The particles are always stable. In our particular class of models these particles are fermions with a chiral coupling to the condensate. They may behave as relativistic matter at early times, produce a brief period where they dominate the expansion with w<0 today, and behave as matter at late time. There are no small parameters in our models, which generically lead to dark energy clustering and, depending on the choice of parameters, smoothing of small scale power.Comment: 8 pages, 5 figures; minor update with added refs; version appearing in Phys. Rev.

Symplectic algorithm for constant-pressure molecular dynamics using a Nose-Poincare thermostat

We present a new algorithm for isothermal-isobaric molecular-dynamics simulation. The method uses an extended Hamiltonian with an Andersen piston combined with the Nos'e-Poincar'e thermostat, recently developed by Bond, Leimkuhler and Laird [J. Comp. Phys., 151, (1999)]. This Nos'e-Poincar'e-Andersen (NPA) formulation has advantages over the Nos'e-Hoover-Andersen approach in that the NPA is Hamiltonian and can take advantage of symplectic integration schemes, which lead to enhanced stability for long-time simulations. The equations of motion are integrated using a Generalized Leapfrog Algorithm and the method is easy to implement, symplectic, explicit and time reversible. To demonstrate the stability of the method we show results for test simulations using a model for aluminum.Comment: 7 page

Adjusting the melting point of a model system via Gibbs-Duhem integration: application to a model of Aluminum

Model interaction potentials for real materials are generally optimized with respect to only those experimental properties that are easily evaluated as mechanical averages (e.g., elastic constants (at T=0 K), static lattice energies and liquid structure). For such potentials, agreement with experiment for the non-mechanical properties, such as the melting point, is not guaranteed and such values can deviate significantly from experiment. We present a method for re-parameterizing any model interaction potential of a real material to adjust its melting temperature to a value that is closer to its experimental melting temperature. This is done without significantly affecting the mechanical properties for which the potential was modeled. This method is an application of Gibbs-Duhem integration [D. Kofke, Mol. Phys.78, 1331 (1993)]. As a test we apply the method to an embedded atom model of aluminum [J. Mei and J.W. Davenport, Phys. Rev. B 46, 21 (1992)] for which the melting temperature for the thermodynamic limit is 826.4 +/- 1.3K - somewhat below the experimental value of 933K. After re-parameterization, the melting temperature of the modified potential is found to be 931.5K +/- 1.5K.Comment: 9 pages, 5 figures, 4 table

Direct calculation of the hard-sphere crystal/melt interfacial free energy

We present a direct calculation by molecular-dynamics computer simulation of the crystal/melt interfacial free energy, $\gamma$, for a system of hard spheres of diameter $\sigma$. The calculation is performed by thermodynamic integration along a reversible path defined by cleaving, using specially constructed movable hard-sphere walls, separate bulk crystal and fluid systems, which are then merged to form an interface. We find the interfacial free energy to be slightly anisotropic with $\gamma$ = 0.62$\pm 0.01$, 0.64$\pm 0.01$ and 0.58$\pm 0.01 k_BT/\sigma^2$ for the (100), (110) and (111) fcc crystal/fluid interfaces, respectively. These values are consistent with earlier density functional calculations and recent experiments measuring the crystal nucleation rates from colloidal fluids of polystyrene spheres that have been interpreted [Marr and Gast, Langmuir {\bf 10}, 1348 (1994)] to give an estimate of $\gamma$ for the hard-sphere system of $0.55 \pm 0.02 k_BT/\sigma^2$, slightly lower than the directly determined value reported here.Comment: 4 pages, 4 figures, submitted to Physical Review Letter

Middle Atlantic Outer Continental Shelf Environmental Studies Volume II-B: Chemical and Biological Benchmark Studies

The Middle Atlantic Outer Continental Shelf Environmental Studies is comprised of three volumes. Volume I. Executive Summary. Volume IIA, IIB, IIC and IID. Chemical and Biological Benchmark Studies. Volume III. Geologic Studies. This is the second of four sections of the Chemical and Biological Benchmark Studies CHAPTER 5. BOTTOM SEDIMENTS AND SEDIMENTARY FRAMEWORK by Donald .F. Boesch CHAPTER 6. BENTHIC ECOLOGICAL STUDIES: MACROBENTHOS by Donald F. Boesch CHAPTER 7. BENTHIC ECOLOGICAL STUDIES: MEIOBENTHOS by D.J. Hartzband and Donald F. Boesch CHAPTER 8. BENTHIC ECOLOGICAL STUDIES: FORAMINIFERA by Robert L. Ellison Chapters of this report contain the Institutes\u27 Special Reports in Applied Marine Science and Ocean Engineering No.193,194,195,196