Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userDiscrepancy between Monte-Carlo Results and Analytic Values for the Average Excluded Volume of Rectangular Prisms
126
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We perform Monte Carlo simulations to determine the average excluded volume <V_{ex}> of randomly oriented rectangular prisms, randomly oriented ellipsoids and randomly oriented capped cylinders in 3-D. There is agreement between the analytically obtained <V_{ex}> and the results of simulations for randomly oriented ellipsoids and randomly oriented capped cylinders. However, we find that the <V_{ex}> for randomly oriented prisms obtained from the simulations differs from the analytically obtained results. In particular, for cubes, the percentage difference is 3.92, far exceeding the bounds of statistical error in our simulation.{bf Added in Revision 2: We recently found the cause of the discrepancy between the simulation result and the analytic value of the excluded volume to be the effect of an error in our simulation code. Upon rectification of the simulation code, the simulation yields $ 11.00 pm 0.002 $ as the excluded volume of a pair of randomly oriented cubes of unit volume. The simulation also yields results as predicted by the analytic formula for all other cases of rectangular prisms that we study.}
rate research
Read More
We describe Janus, a massively parallel FPGA-based computer optimized for the simulation of spin glasses, theoretical models for the behavior of glassy materials. FPGAs (as compared to GPUs or many-core processors) provide a complementary approach to massively parallel computing. In particular, our model problem is formulated in terms of binary variables, and floating-point operations can be (almost) completely avoided. The FPGA architecture allows us to run many independent threads with almost no latencies in memory access, thus updating up to 1024 spins per cycle. We describe Janus in detail and we summarize the physics results obtained in four years of operation of this machine; we discuss two types of physics applications: long simulations on very large systems (which try to mimic and provide understanding about the experimental non-equilibrium dynamics), and low-temperature equilibrium simulations using an artificial parallel tempering dynamics. The time scale of our non-equilibrium simulations spans eleven orders of magnitude (from picoseconds to a tenth of a second). On the other hand, our equilibrium simulations are unprecedented both because of the low temperatures reached and for the large systems that we have brought to equilibrium. A finite-time scaling ansatz emerges from the detailed comparison of the two sets of simulations. Janus has made it possible to perform spin-glass simulations that would take several decades on more conventional architectures. The paper ends with an assessment of the potential of possible futu
We consider high dimensional random optimization problems where the dynamical variables are subjected to non-convex excluded volume constraints. We focus on the case in which the cost function is a simple quadratic cost and the excluded volume constraints are modeled by a perceptron constraint satisfaction problem. We show that depending on the density of constraints, one can have different situations. If the number of constraints is small, one typically has a phase where the ground state of the cost function is unique and sits on the boundary of the island of configurations allowed by the constraints. In this case there is an hypostatic number of constraints that are marginally satisfied. If the number of constraints is increased one enters in a glassy phase where the cost function has many local minima sitting again on the boundary of the regions of allowed configurations. At the phase transition point the total number of constraints that are marginally satisfied becomes equal to the number of degrees of freedom in the problem and therefore we say that these minima are isostatic. We conjecture that increasing further the constraints the system stays isostatic up to the point where the volume of available phase space shrinks to zero. We derive our results using the replica method and we also analyze a dynamical algorithm, the Karush-Kuhn-Tucker algorithm, through dynamical mean field theory and we show how to recover the results of the replica approach in the replica symmetric phase.
The key idea behind the renormalization group (RG) transformation is that properties of physical systems with very different microscopic makeups can be characterized by a few universal parameters. However, finding the optimal RG transformation remains difficult due to the many possible choices of the weight factors in the RG procedure. Here we show, by identifying the conditional distribution in the restricted Boltzmann machine (RBM) and the weight factor distribution in the RG procedure, an optimal real-space RG transformation can be learned without prior knowledge of the physical system. This neural Monte Carlo RG algorithm allows for direct computation of the RG flow and critical exponents. This scheme naturally generates a transformation that maximizes the real-space mutual information between the coarse-grained region and the environment. Our results establish a solid connection between the RG transformation in physics and the deep architecture in machine learning, paving the way to further interdisciplinary research.
Ising Monte Carlo simulations of the random-field Ising system Fe(0.80)Zn(0.20)F2 are presented for H=10T. The specific heat critical behavior is consistent with alpha approximately 0 and the staggered magnetization with beta approximately 0.25 +- 0.03.
We review the physics of the Bose-Hubbard model with disorder in the chemical potential focusing on recently published analytical arguments in combination with quantum Monte Carlo simulations. Apart from the superfluid and Mott insulator phases that can occur in this system without disorder, disorder allows for an additional phase, called the Bose glass phase. The topology of the phase diagram is subject to strong theorems proving that the Bose Glass phase must intervene between the superfluid and the Mott insulator and implying a Griffiths transition between the Mott insulator and the Bose glass. The full phase diagrams in 3d and 2d are discussed, and we zoom in on the insensitivity of the transition line between the superfluid and the Bose glass in the close vicinity of the tip of the Mott insulator lobe. We briefly comment on the established and remaining questions in the 1d case, and give a short overview of numerical work on related models.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
