Search results for "Algorithm"
showing 10 items of 4887 documents
Monte Carlo simulation of phase separation and clustering in the ABV model
1991
As a model for a binary alloy undergoing an unmixing phase transition, we consider a square lattice where each site can be either taken by an A atom, a B atom, or a vacancy (V), and there exists a repulsive interaction between AB nearest neighbor pairs. Starting from a random initial configuration, unmixing proceeds via random jumps of A atoms or B atoms to nearest neighbor vacant sites. In the absence of any interaction, these jumps occur at jump ratesΓ A andΓ B, respectively. For a small concentration of vacancies (c v=0.04) the dynamics of the structure factorS(k,t) and its first two momentsk 1(t),k 2 2 (t) is studied during the early stages of phase separation, for several choices of co…
Multicanonical multigrid Monte Carlo method.
1994
To further improve the performance of Monte Carlo simulations of first-order phase transitions we propose to combine the multicanonical approach with multigrid techniques. We report tests of this proposition for the d-dimensional ${\mathrm{\ensuremath{\Phi}}}^{4}$ field theory in two different situations. First, we study quantum tunneling for d=1 in the continuum limit, and second, we investigate first-order phase transitions for d=2 in the infinite volume limit. Compared with standard multicanonical simulations we obtain improvement factors of several, and of about one order of magnitude, respectively.
Finite-Size Scaling Study of the Simple Cubic Three-State Potts Glass
1991
During the last few years the Potts glass model has attracted more and more attention. It is considered as a first step towards modelling the phase transition of structural and orientational glasses. A mean-field approach /1/ predicts a low temperature behavior completely different from what is known from Ising spin glasses /2/. But short range models differ markedly from mean-field-predictions. So it is natural to ask, how the short range Potts glass behaves. Especially the question of the lower critical dimension d l is important, below which a finite temperature transition ceases to occur. We tried to answer this by combining Monte-Carlo simulations with a finite-size scaling analysis. T…
Comparison between transmission and scattering spectrum reconstruction methods based on EPID images.
2013
Numerous improved physics-based methods for Linac photon spectra reconstruction have been published; some of them are based on transmission data analysis and others on scattering data. In this work, the two spectrum unfolding approaches are compared in order to experimentally validate its robustness and to determine which is the optimal methodology for application on a clinical quality assurance routine. Both studied methods are based on EPID images generated when the incident photon beam impinges onto plastic blocks. The distribution of transmitted/scatter radiation produced by this object centered at the beam field size was measured. Measurements were performed using a 6 MeV photon beam p…
Modelling systems of classical/quantum identical particles by focusing on algorithms
2012
A procedure modelling ideal classical and quantum gases is discussed. The proposed approach is mainly based on the idea that modelling and algorithm analysis can provide a deeper understanding of particularly complex physical systems. Appropriate representations and physical models able to mimic possible pseudo-mechanisms of functioning and having predictive validity are developed.
Trajectory Statistics of Confined L\'evy Flights and Boltzmann-type Equilibria
2013
We analyze a specific class of random systems that are driven by a symmetric L\'{e}vy stable noise, where Langevin representation is absent. In view of the L\'{e}vy noise sensitivity to environmental inhomogeneities, the pertinent random motion asymptotically sets down at the Boltzmann-type equilibrium, represented by a probability density function (pdf) $\rho_*(x) \sim \exp [-\Phi (x)]$. Here, we infer pdf $\rho (x,t)$ based on numerical path-wise simulation of the underlying jump-type process. A priori given data are jump transition rates entering the master equation for $\rho (x,t)$ and its target pdf $\rho_*(x)$. To simulate the above processes, we construct a suitable modification of t…
A Fast and Very Accurate Approach to the Computation of Microlensing Magnification Patterns Based on Inverse Polygon Mapping
2006
A new method of calculating microlensing magnification patterns is proposed that is based on the properties of the backward gravitational lens mapping of a lattice of polygonal cells defined at the image plane. To a first-order approximation, the local linearity of the transformation allows us to compute the contribution of each image-plane cell to the magnification by apportioning the area of the inverse image of the cell (transformed cell) among the source-plane pixels covered by it. Numerical studies in the κ = 0.1-0.8 range of mass surface densities demonstrate that this method (provided with an exact algorithm for distributing the area of the transformed cells among the source-plane pi…
Probabilistic Fault-Tolerant Universal Quantum Computation and Sampling Problems in Continuous Variables
2019
Continuous-Variable (CV) devices are a promising platform for demonstrating large-scale quantum information protocols. In this framework, we define a general quantum computational model based on a CV hardware. It consists of vacuum input states, a finite set of gates - including non-Gaussian elements - and homodyne detection. We show that this model incorporates encodings sufficient for probabilistic fault-tolerant universal quantum computing. Furthermore, we show that this model can be adapted to yield sampling problems that cannot be simulated efficiently with a classical computer, unless the polynomial hierarchy collapses. This allows us to provide a simple paradigm for short-term experi…
Monte Carlo calculation of dose rate distributions around the Walstam CDC.K-type137Cs sources
2001
Basic dosimetric data for the Walstam CDC.K-type low dose rate 137Cs sources in water have been calculated using Monte Carlo techniques. These sources, CDC.K1 -K3 and CDC.K4, are widely used in a range of applicators and moulds for the treatment of intracavitary and superficial cancers. Our purpose is to improve existing data about these sources using the Monte Carlo simulation code GEANT3. Absolute dose rate distributions in water have been calculated around these sources and are presented as conventional 2D Cartesian look-up tables. Also the AAPM Task Group 43 formalism for dose calculation has been applied. The calculated dose rate constant for the CDC.K1-K3 source is A = 1.106 +/- 0.001…
Method specific Cholesky decomposition : Coulomb and exchange energies
2008
We present a novel approach to the calculation of the Coulomb and exchange contributions to the total electronic energy in self consistent field and density functional theory. The numerical procedure is based on the Cholesky decomposition and involves decomposition of specific Hadamard product matrices that enter the energy expression. In this way, we determine an auxiliary basis and obtain a dramatic reduction in size as compared to the resolution of identity (RI) method. Although the auxiliary basis is determined from the energy expression, we have complete control of the errors in the gradient or Fock matrix. Another important advantage of this method specific Cholesky decomposition is t…