Source author record

David P. Landau

David P. Landau appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

20works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

20 published item(s)

preprint2022arXiv

Quantum algorithm for Wang-Landau sampling

It has been shown that the Metropolis algorithm can be implemented on quantum computers in a way that avoids the sign problem. However, flat histogram techniques are often preferred as they don't suffer from the same limitations that afflict Metropolis for problems of real-world interest and provide a host of other benefits. In particular, the Wang-Landau method is known for its efficiency and accuracy. In this work we design, implement, and validate a quantum algorithm for Wang-Landau sampling, greatly expanding the range of quantum many body problems solvable by Monte Carlo simulation.

preprint2020arXiv

High-resolution Monte Carlo study of the order-parameter distribution of the three-dimensional Ising model

We apply extensive Monte Carlo simulations to study the probability distribution $P(m)$ of the order parameter $m$ for the simple cubic Ising model with periodic boundary condition at the transition point. Sampling is performed with the Wolff cluster flipping algorithm, and histogram reweighting together with finite-size scaling analyses are then used to extract a precise functional form for the probability distribution of the magnetization, $P(m)$, in the thermodynamic limit. This form should serve as a benchmark for other models in the three-dimensional Ising Universality class.

preprint2016arXiv

Magnetic phase transition in coupled spin-lattice systems: A replica-exchange Wang-Landau study

Coupled, dynamical spin-lattice models provide a unique test ground for simulations investigating the finite-temperature magnetic properties of materials under the direct influence of the lattice vibrations. These models are constructed by combining a coordinate-dependent interatomic potential with a Heisenberg-like spin Hamiltonian, facilitating the treatment of both the atomic coordinates and spins as explicit phase variables. Using a model parameterized for bcc iron, we study the magnetic phase transition in these complex systems via the recently introduced, massively parallel replica-exchange Wang-Landau Monte Carlo method. Comparison with the results obtained from rigid lattice (spin only) simulations show that the transition temperature as well as the amplitude of the peak in the specific heat curve is marginally affected by the lattice vibrations. Moreover, the results were found to be sensitive to the particular choice of the interatomic potential.

preprint2016arXiv

Reinventing atomistic magnetic simulations with spin-orbit coupling

We propose a powerful extension to combined molecular and spin dynamics that fully captures the coupling between the atomic and spin subsystems via spin-orbit interactions. Its foundation is the inclusion of the local magnetic anisotropies that arise as a consequence of the lattice symmetry breaking due to phonons or defects. We demonstrate that our extension enables the exchange of angular momentum between the atomic and spin subsystems, which is critical to the challenges arising in the study of fluctuations and non-equilibrium processes in complex, natural, and engineered magnetic materials.

preprint2015arXiv

Exploring Replica-Exchange Wang-Landau sampling in higher-dimensional parameter space

We considered a higher-dimensional extension for the replica-exchange Wang-Landau algorithm to perform a random walk in the energy and magnetization space of the two-dimensional Ising model. This hybrid scheme combines the advantages of Wang-Landau and Replica-Exchange algorithms, and the one-dimensional version of this approach has been shown to be very efficient and to scale well, up to several thousands of computing cores. This approach allows us to split the parameter space of the system to be simulated into several pieces and still perform a random walk over the entire parameter range, ensuring the ergodicity of the simulation. Previous work, in which a similar scheme of parallel simulation was implemented without using replica exchange and with a different way to combine the result from the pieces, led to discontinuities in the final density of states over the entire range of parameters. From our simulations, it appears that the replica-exchange Wang-Landau algorithm is able to overcome this difficulty, allowing exploration of higher parameter phase space by keeping track of the joint density of states.

preprint2015arXiv

Versatile approach to access the low temperature thermodynamics of lattice polymers and proteins

We show that Wang-Landau sampling, combined with suitable Monte Carlo trial moves, provides a powerful method for both the ground state search and the determination of the density of states for the hydrophobic-polar (HP) protein model and the interacting self-avoiding walk (ISAW) model for homopolymers. We obtained accurate estimates of thermodynamic quantities for HP sequences with $>100$ monomers and for ISAWs up to $>500$ monomers. Our procedure possesses an intrinsic simplicity and overcomes the limitations inherent in more tailored approaches making it interesting for a broad range of protein and polymer models.

preprint2014arXiv

A new paradigm for petascale Monte Carlo simulation: Replica exchange Wang-Landau sampling

We introduce a generic, parallel Wang-Landau method that is naturally suited to implementation on massively parallel, petaflop supercomputers. The approach introduces a replica-exchange framework in which densities of states for overlapping sub-windows in energy space are determined iteratively by traditional Wang-Landau sampling. The advantages and general applicability of the method are demonstrated for several distinct systems that possess discrete or continuous degrees of freedom, including those with complex free energy landscapes and topological constraints.

preprint2014arXiv

Identifying transitions in finite systems by means of partition function zeros and microcanonical inflection-point analysis: A comparison for elastic flexible polymers

For the estimation of transition points of finite elastic, flexible polymers with chain lengths from $13$ to $309$ monomers, we compare systematically transition temperatures obtained by the Fisher partition function zeros approach with recent results from microcanonical inflection-point analysis. These methods rely on accurate numerical estimates of the density of states, which have been obtained by advanced multicanonical Monte Carlo sampling techniques. Both the Fisher zeros method and microcanonical inflection-point analysis yield very similar results and enable the unique identification of transition points in finite systems, which is typically impossible in the conventional canonical analysis of thermodynamic quantities.

preprint2014arXiv

Replica-exchange Wang-Landau sampling: Pushing the limits of Monte Carlo simulations in materials sciences

We describe the study of thermodynamics of materials using replica-exchange Wang-Landau (REWL) sampling, a generic framework for massively parallel implementations of the Wang-Landau Monte Carlo method. To evaluate the performance and scalability of the method, we investigate the magnetic phase transition in body-centered cubic (bcc) iron using the classical Heisenberg model parametrized with first principles calculations. We demonstrate that our framework leads to a significant speedup without compromising the accuracy and precision and facilitates the study of much larger systems than is possible with its serial counterpart.

preprint2014arXiv

Scalable replica-exchange framework for Wang-Landau sampling

We investigate a generic, parallel replica-exchange framework for Monte Carlo simulations based on the Wang-Landau method. To demonstrate its advantages and general applicability for massively parallel simulations of complex systems, we apply it to lattice spin models, the self-assembly process in amphiphilic solutions, and the adsorption of molecules on surfaces. While of general, current interest, the latter phenomena are challenging to study computationally because of multiple structural transitions occurring over a broad temperature range. We show how the parallel framework facilitates simulations of such processes and, without any loss of accuracy or precision, gives a significant speedup and allows for the study of much larger systems and much wider temperature ranges than possible with single-walker methods.

preprint2013arXiv

A generic, hierarchical framework for massively parallel Wang-Landau sampling

We introduce a parallel Wang-Landau method based on the replica-exchange framework for Monte Carlo simulations. To demonstrate its advantages and general applicability for simulations of complex systems, we apply it to different spin models including spin glasses, the Ising model and the Potts model, lattice protein adsorption, and the self-assembly process in amphiphilic solutions. Without loss of accuracy, the method gives significant speed-up and potentially scales up to petaflop machines.

preprint2013arXiv

Exploring new frontiers in statistical physics with a new, parallel Wang-Landau framework

Combining traditional Wang-Landau sampling for multiple replica systems with an exchange of densities of states between replicas, we describe a general framework for simulations on massively parallel Petaflop supercomputers. The advantages and general applicability of the method for simulations of complex systems are demonstrated for the classical 2D Potts spin model featuring a strong first-order transition and the self-assembly of lipid bilayers in amphiphilic solutions in a continuous model.

preprint2013arXiv

Generic folding and transition hierarchies for surface adsorption of hydrophobic-polar lattice model proteins

The thermodynamic behavior and structural properties of hydrophobic-polar (HP) lattice proteins interacting with attractive surfaces are studied by means of Wang-Landau sampling. Three benchmark HP sequences (48mer, 67mer, and 103mer) are considered with different types of surfaces, each of which attract either all monomers, only hydrophobic (H) monomers, or only polar (P) monomers, respectively. The diversity of folding behavior in dependence of surface strength is discussed. Analyzing the combined patterns of various structural observables, such as, e.g., the derivatives of the numbers of surface contacts, together with the specific heat, we are able to identify generic categories of folding and transition hierarchies. We also infer a connection between these transition categories and the relative surface strengths, i.e., the ratio of the surface attractive strength to the interchain attraction among H monomers. The validity of our proposed classification scheme is reinforced by the analysis of additional benchmark sequences. We thus believe that the folding hierarchies and identification scheme are generic for HP proteins interacting with attractive surfaces, regardless of chain length, sequence, or surface attraction.

preprint2013arXiv

Surface adsorption of lattice HP proteins: Thermodynamics and structural transitions using Wang-Landau sampling

Wang-Landau sampling has been applied to investigate the thermodynamics and structural properties of a lattice hydrophobic-polar heteropolymer (the HP protein model) interacting with an attractive substrate. For simplicity, we consider a short HP sequence consisting of only 36 monomers interacting with a substrate which attracts all monomers in the sequence. The conformational "phase transitions" have been identified by a canonical analysis of the specific heat and suitable structural observables. Three major "transitions", namely, adsorption, hydrophobic core formation and "flattening" of adsorbed structures, are observed. Depending on the surface attractive strength relative to the intra-protein attraction among the H monomers, these processes take place in different sequences upon cooling.

preprint2013arXiv

Unraveling the beautiful complexity of simple lattice model polymers and proteins using Wang-Landau sampling

We describe a class of "bare bones" models of homopolymers which undergo coil-globule collapse and proteins which fold into their native states in free space or into denatured states when captured by an attractive substrate as the temperature is lowered. We then show how, with the use of a properly chosen trial move set, Wang-Landau Monte Carlo sampling can be used to study the rough free energy landscape and ground (native) states of these intriguingly simple systems and thus elucidate their thermodynamic complexity.

preprint2012arXiv

Fictitious excitations in the classical Heisenberg antiferromagnet on the kagome lattice

Using an advanced Monte Carlo algorithm and high precision spin dynamics simulation, we investigated the dynamical behavior of the Heisenberg antiferromagnet on a kagome lattice at extremely low temperatures. We demonstrate that the detection and correct identification of propagating excitations depends crucially on the choice of coordinates and we show how modes are displaced in Fourier space if a single, global Cartesian coordinate system is chosen.

preprint2012arXiv

Optimized Wang-Landau sampling of lattice polymers: Ground state search and folding thermodynamics of HP model proteins

Coarse-grained (lattice-) models have a long tradition in aiding efforts to decipher the physical or biological complexity of proteins. Despite the simplicity of these models, however, numerical simulations are often computationally very demanding and the quest for efficient algorithms is as old as the models themselves. Expanding on our previous work [T. Wüst and D. P. Landau, Phys. Rev. Lett. 102, 178101 (2009)], we present a complete picture of a Monte Carlo method based on Wang-Landau sampling in combination with efficient trial moves (pull, bond-rebridging and pivot moves) which is particularly suited to the study of models such as the hydrophobic-polar (HP) lattice model of protein folding. With this generic and fully blind Monte Carlo procedure, all currently known putative ground states for the most difficult benchmark HP sequences could be found. For most sequences we could also determine the entire energy density of states and, together with suitably designed structural observables, explore the thermodynamics and intricate folding behavior in the virtually inaccessible low-temperature regime. We analyze the differences between random and protein-like heteropolymers for sequence lengths up to 500 residues. Our approach is powerful both in terms of robustness and speed, yet flexible and simple enough for the study of many related problems in protein folding.

preprint2011arXiv

Microcanonical entropy inflection points: Key to systematic understanding of transitions in finite systems

We introduce a systematic classification method for the analogs of phase transitions in finite systems. This completely general analysis, which is applicable to any physical system and extends towards the thermodynamic limit, is based on the microcanonical entropy and its energetic derivative, the inverse caloric temperature. Inflection points of this quantity signal cooperative activity and thus serve as distinct indicators of transitions. We demonstrate the power of this method through application to the long-standing problem of liquid-solid transitions in elastic, flexible homopolymers.

preprint2011arXiv

Monte Carlo simulations of the HP model (the "Ising model" of protein folding)

Using Wang-Landau sampling with suitable Monte Carlo trial moves (pull moves and bond-rebridging moves combined) we have determined the density of states and thermodynamic properties for a short sequence of the HP protein model. For free chains these proteins are known to first undergo a collapse "transition" to a globule state followed by a second "transition" into a native state. When placed in the proximity of an attractive surface, there is a competition between surface adsorption and folding that leads to an intriguing sequence of "transitions". These transitions depend upon the relative interaction strengths and are largely inaccessible to "standard" Monte Carlo methods.

preprint2010arXiv

The Dynamic Scaling Study of Vapor Deposition Polymerization: A Monte Carlo Approach

The morphological scaling properties of linear polymer films grown by vapor deposition polymerization (VDP) are studied by 1+1D Monte Carlo simulations. The model implements the basic processes of random angle ballistic deposition ($F$), free-monomer diffusion ($D$) and monomer adsorption along with the dynamical processes of polymer chain initiation, extension, and merger. The ratio $G=D/F$ is found to have a strong influence on the polymer film morphology. Spatial and temporal behavior of kinetic roughening has been extensively studied using finite-length scaling and height-height correlations $H(r,t)$. The scaling analysis has been performed within the no-overhang approximation and the scaling behaviors at local and global length scales were found to be very different. The global and local scaling exponents for morphological evolution have been evaluated for varying free-monomer diffusion by growing the films at $G$ = $10$, $10^2$, $10^3$, and $10^4$ and fixing the deposition flux $F$. With an increase in $G$ from $10$ to $10^4$, the average growth exponent $β\approx 0.50$ was found to be invariant, where as the global roughness exponent $α_{g}$ decreased from $0.87(1)$ to $0.73(1)$ along with a corresponding decrease in the global dynamic exponent $z_g$ from $1.71(1)$ to $1.38(2)$. The global scaling exponents were observed to follow the dynamic scaling hypothesis, $z_g=α_{g}/β$. With a similar increase in $G$ however, the average local roughness exponent ${α_{l}}$ remained close to $0.46$ and the anomalous growth exponent $β_{\ast}$ decreased from 0.23(4) to 0.18(8). The interfaces display anomalous scaling and multiscaling in the relevant height-height correlations. The variation of $H(r,t)$ with deposition time $t$ indicates non-stationary growth.