BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Werner Krauth

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

ResearcherUnclaimed source record
cond-mat.stat-mechcond-mat.softcond-mat.quant-gasphysics.comp-phcond-mat.dis-nncond-mat.str-elcond-mat.supr-conComputationhep-thquant-ph

Catalog footprint

What is connected

29works
10topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

29 published item(s)

preprint2022arXiv

Direction-sweep Markov chains

We discuss a non-reversible, lifted Markov-chain Monte Carlo (MCMC) algorithm for particle systems in which the direction of proposed displacements is changed deterministically. This algorithm sweeps through directions analogously to the popular MCMC sweep methods for particle or spin indices. Direction-sweep MCMC can be applied to a wide range of original reversible or non-reversible Markov chains, such as the Metropolis algorithm or the event-chain Monte Carlo algorithm. For a single two-dimensional dipole, we consider direction-sweep MCMC in the limit where restricted equilibrium is reached among the accessible configurations before changing the direction. We show rigorously that direction-sweep MCMC leaves the stationary probability distribution unchanged, and that it profoundly modifies the Markov-chain trajectory. Long excursions, with persistent rotation in one direction, alternate with long sequences of rapid zigzags resulting in persistent rotation in the opposite direction in the limit of small direction increments. The mapping to a Langevin equation then yields the exact scaling of excursions while the zigzags are described through a non-linear differential equation that is solved exactly. We show that the direction-sweep algorithm can have shorter mixing times than the algorithms with random updates of directions. We point out possible applications of direction-sweep MCMC in polymer physics and in molecular simulation.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Event-chain Monte Carlo: foundations, applications, and prospects

This review treats the mathematical and algorithmic foundations of non-reversible Markov chains in the context of event-chain Monte Carlo (ECMC), a continuous-time lifted Markov chain that employs the factorized Metropolis algorithm. It analyzes a number of model applications, and then reviews the formulation as well as the performance of ECMC in key models in statistical physics. Finally, the review reports on an ongoing initiative to apply the method to the sampling problem in molecular simulation, that is, to real-world models of peptides, proteins, and polymers in aqueous solution.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Hard-disk dipoles and non-reversible Markov chains

We benchmark event-chain Monte Carlo (ECMC) algorithms for tethered hard-disk dipoles in two dimensions in view of application of ECMC to water models in molecular simulation. We characterize the rotation dynamics of dipoles through the integrated autocorrelation times of the polarization. The non-reversible straight, reflective, forward, and Newtonian ECMC algorithms are all event-driven, and they differ only in their update rules at event times. They realize considerable speedups with respect to the local reversible Metropolis algorithm. We also find significant speed differences among the ECMC variants. Newtonian ECMC appears particularly well-suited for overcoming the dynamical arrest that has plagued straight ECMC for three-dimensional dipolar models with Coulomb interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Large-scale dynamics of event-chain Monte Carlo

Event-chain Monte Carlo (ECMC) accelerates the sampling of hard-sphere systems, and has been generalized to the potentials used in classical molecular simulation. Rather than imposing detailed balance on transition probabilities, the method enforces a weaker global-balance condition in order to guarantee convergence to equilibrium. In this paper we generalize the factor-field variant of ECMC to higher space dimensions. In the two-dimensional fluid phase, factor-field ECMC saturates the lower bound $z=0$ for the dynamical scaling exponent for local dynamics, whereas molecular dynamics is characterized by $z=1$ and local Metropolis Monte Carlo by $z=2$. In the presence of hexatic order, factor fields are not found to speed up the convergence. We indicate applications to the physics of glasses, and note that generalizations of factor fields could couple to orientational order.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Sparse hard-disk packings and local Markov chains

We propose locally stable sparse hard-disk packings, as introduced by Böröczky, as a model for the analysis and benchmarking of Markov-chain Monte Carlo (MCMC) algorithms. We first generate such packings in a square box with periodic boundary conditions and analyze their properties. We then study how local MCMC algorithms, namely the Metropolis algorithm and several versions of event-chain Monte Carlo (ECMC), escape from configurations that are obtained by slightly reducing all disk radii by a relaxation parameter. A scaling analysis is confirmed by simulation results. We obtain two classes of ECMC, one in which the escape time varies algebraically with the relaxation parameter (as for the local Metropolis algorithm) and another in which the escape time scales as the logarithm of the relaxation parameter. We discuss the connectivity of the hard-disk sample space, the ergodicity of local MCMC algorithms, as well as the meaning of packings in the context of the NPT ensemble. Our work is accompanied by open-source, arbitrary-precision software for Böröczky packings (in Python) and for straight, reflective, forward, and Newtonian ECMC (in Go).

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Multithreaded event-chain Monte Carlo with local times

We present a multithreaded event-chain Monte Carlo algorithm (ECMC) for hard spheres. Threads synchronize at infrequent breakpoints and otherwise scan for local horizon violations. Using a mapping onto absorbing Markov chains, we rigorously prove the correctness of a sequential-consistency implementation for small test suites. On x86 and ARM processors, a C++ (OpenMP) implementation that uses compare-and-swap primitives for data access achieves considerable speed-up with respect to single-threaded code. The generalized birthday problem suggests that for the number of threads scaling as the square root of the number of spheres, the horizon-violation probability remains small for a fixed simulation time. We provide C++ and Python open-source code that reproduces all our results.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Cell-veto Monte Carlo algorithm for long-range systems

We present a rigorous efficient event-chain Monte Carlo algorithm for long-range interacting particle systems. Using a cell-veto scheme within the factorized Metropolis algorithm, we compute each single-particle move with a fixed number of operations. For slowly decaying potentials such as Coulomb interactions, screening line charges allow us to take into account periodic boundary conditions. We discuss the performance of the cell-veto Monte Carlo algorithm for general inverse-power-law potentials, and illustrate how it provides a new outlook on one of the prominent bottlenecks in large-scale atomistic Monte Carlo simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Momentum distribution in the unitary Bose gas from first principles

We consider a realistic bosonic N-particle model with unitary interactions relevant for Efimov physics. Using quantum Monte Carlo methods, we find that the critical temperature for Bose-Einstein condensation is decreased with respect to the ideal Bose gas. We also determine the full momentum distribution of the gas, including its universal asymptotic behavior, and compare this crucial observable to recent experimental data. Similar to the experiments with different atomic species, differentiated solely by a three-body length scale, our model only depends on a single parameter. We establish a weak influence of this parameter on physical observables. In current experiments, the thermodynamic instability of our model from the atomic gas towards an Efimov liquid could be masked by the dynamical instability due to three-body losses.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

2D Melting: From Liquid-Hexatic Coexistence to Continuous Transitions

The phase diagram of two-dimensional continuous particle systems is studied using Event-Chain Monte Carlo. For soft disks with repulsive power-law interactions $\propto r^{-n}$ with $n \gtrsim 6$, the recently established hard-disk melting scenario ($n \to \infty$) holds: a first-order liquid-hexatic and a continuous hexatic-solid transition are identified. Close to $n = 6$, the coexisting liquid exhibits very long orientational correlations, and positional correlations in the hexatic are extremely short. For $n\lesssim 6$, the liquid-hexatic transition is continuous, with correlations consistent with the Kosterlitz-Thouless-Halperin-Nelson-Yong (KTHNY) scenario. To illustrate the generality of these results, we demonstrate that Yukawa particles likewise may follow either the KTHNY or the hard-disk melting scenario, depending on the Debye-Hückel screening length as well as on the temperature.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Event-chain algorithm for the Heisenberg model: Evidence for $z \simeq 1$ dynamic scaling

We apply the event-chain Monte Carlo algorithm to the three-dimensional ferromagnetic Heisenberg model. The algorithm is rejection-free and also realizes an irreversible Markov chain that satisfies global balance. The autocorrelation functions of the magnetic susceptibility and the energy indicate a dynamical critical exponent $z \approx 1$ at the critical temperature, while that of the magnetization does not measure the performance of the algorithm. This seems to be the first report that the event-chain Monte Carlo algorithm substantially reduces the dynamical critical exponent from the conventional value of $z\simeq 2$.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Hard-sphere melting and crystallization with event-chain Monte Carlo

We simulate crystallization and melting with local Monte Carlo (LMC), event-chain Monte Carlo (ECMC), and with event-driven molecular dynamics (EDMD) in systems with up to one million three-dimensional hard spheres. We illustrate that our implementations of the three algorithms rigorously coincide in their equilibrium properties. We then study nucleation in the NVE ensemble from the fcc crystal into the homogeneous liquid phase and from the liquid into the homogeneous crystal. ECMC and EDMD both approach equilibrium orders of magnitude faster than LMC. ECMC is also notably faster than EDMD, especially for the equilibration into a crystal from a disordered initial condition at high density. ECMC can be trivially implemented for hard-sphere and for soft-sphere potentials, and we suggest possible applications of this algorithm for studying jamming and the physics of glasses, as well as disordered systems.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Liquid-solid transitions in the three-body hard-core model

We determine the phase diagram for a generalisation of two-and three-dimensional hard spheres: a classical system with three-body interactions realised as a hard cut-off on the mean-square distance for each triplet of particles. Quantum versions of this model are important in the context of the unitary Bose gas, which is currently under close theoretical and experimental scrutiny. In two dimensions, the three-body hard-core model possesses a conventional atomic liquid phase and a peculiar solid phase formed by dimers. These dimers interact effectively as hard disks. In three dimensions, the solid phase consists of isolated atoms that arrange in a simple-hexagonal lattice.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Efimov-driven phase transitions of the unitary Bose gas

In quantum physics, Efimov trimers are bound states of three particles that fall apart like Borromean rings when one of them is removed. Initially predicted in nuclear physics, these striking bosonic states are hard to observe, but the "unitary" interactions at which they form is commonly realized in current cold atoms experiments. There, they set the stage for a new class of universal physics: Two-body interactions are all but invisible, but three-body effects allow the emergence of a largely uncharted new world of many-particle bound states. Three-particle systems were characterized theoretically, and the ground-state properties of small unitary clusters computed numerically, but the macroscopic many-body behaviour has remained unknown. Here we show, using a Path-Integral Monte Carlo algorithm backed up by theoretical arguments, that the unitary Bose gas presents a first-order phase transition from a normal gas to a superfluid Efimov liquid. The normal gas is very well described by the available virial coefficients. At unitarity, the phase diagram of the bosonic system is universal in rescaled pressure and temperature. A triple point separates the normal gas, the superfluid Efimov liquid, and a third phase, the conventional superfluid gas. These two superfluid phases are separated by a critical line that ends in a critical point at high temperature. This rich phase diagram should allow for a number of experimental protocols that would probe these universal transitions between the normal gas, the superfluid gas, and the superfluid Efimov liquid.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Generalized event-chain Monte Carlo: Constructing rejection-free global-balance algorithms from infinitesimal steps

In this article, we present an event-driven algorithm that generalizes the recent hard-sphere event-chain Monte Carlo method without introducing discretizations in time or in space. A factorization of the Metropolis filter and the concept of infinitesimal Monte Carlo moves are used to design a rejection-free Markov-chain Monte Carlo algorithm for particle systems with arbitrary pairwise interactions. The algorithm breaks detailed balance, but satisfies maximal global balance and performs better than the classic, local Metropolis algorithm in large systems. The new algorithm generates a continuum of samples of the stationary probability density. This allows us to compute the pressure and stress tensor as a byproduct of the simulation without any additional computations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Hard-disk equation of state: First-order liquid-hexatic transition in two dimensions with three simulation methods

We report large-scale computer simulations of the hard-disk system at high densities in the region of the melting transition. Our simulations reproduce the equation of state, previously obtained using the event-chain Monte Carlo algorithm, with a massively parallel implementation of the local Monte Carlo method and with event-driven molecular dynamics. We analyze the relative performance of these simulation methods to sample configuration space and approach equilibrium. Our results confirm the first-order nature of the melting phase transition in hard disks. Phase coexistence is visualized for individual configurations via the orientational order parameter field. The analysis of positional order confirms the existence of the hexatic phase.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Sampling from a polytope and hard-disk Monte Carlo

The hard-disk problem, the statics and the dynamics of equal two-dimensional hard spheres in a periodic box, has had a profound influence on statistical and computational physics. Markov-chain Monte Carlo and molecular dynamics were first discussed for this model. Here we reformulate hard-disk Monte Carlo algorithms in terms of another classic problem, namely the sampling from a polytope. Local Markov-chain Monte Carlo, as proposed by Metropolis et al. in 1953, appears as a sequence of random walks in high-dimensional polytopes, while the moves of the more powerful event-chain algorithm correspond to molecular dynamics evolution. We determine the convergence properties of Monte Carlo methods in a special invariant polytope associated with hard-disk configurations, and the implications for convergence of hard-disk sampling. Finally, we discuss parallelization strategies for event-chain Monte Carlo and present results for a multicore implementation.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Supplementary information to "Efimov-driven phase transitions of the unitary Bose gas"

In quantum physics, Efimov trimers are bound states of three particles that fall apart like Borromean rings when one of them is removed. Initially predicted in nuclear physics, these striking bosonic states are hard to observe, but the "unitary" interactions at which they form is commonly realized in current cold atoms experiments. There, they set the stage for a new class of universal physics: Two-body interactions are all but invisible, but three-body effects allow the emergence of a largely uncharted new world of many-particle bound states. Three-particle systems were characterized theoretically, and the ground-state properties of small unitary clusters computed numerically, but the macroscopic many-body behaviour has remained unknown. Here we show, using a Path-Integral Monte Carlo algorithm backed up by theoretical arguments, that the unitary Bose gas presents a first-order phase transition from a normal gas to a superfluid Efimov liquid. The normal gas is very well described by the available virial coefficients. At unitarity, the phase diagram of the bosonic system is universal in rescaled pressure and temperature. A triple point separates the normal gas, the superfluid Efimov liquid, and a third phase, the conventional superfluid gas. These two superfluid phases are separated by a critical line that ends in a critical point at high temperature. This rich phase diagram should allow for a number of experimental protocols that would probe these universal transitions between the normal gas, the superfluid gas, and the superfluid Efimov liquid.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Damage spreading and coupling in Markov chains

In this paper, we relate the coupling of Markov chains, at the basis of perfect sampling methods, with damage spreading, which captures the chaotic nature of stochastic dynamics. For two-dimensional spin glasses and hard spheres we point out that the obstacle to the application of perfect-sampling schemes is posed by damage spreading rather than by the survey problem of the entire configuration space. We find dynamical damage-spreading transitions deeply inside the paramagnetic and liquid phases, and show that critical values of the transition temperatures and densities depend on the coupling scheme. We discuss our findings in the light of a classic proof that for arbitrary Monte Carlo algorithms damage spreading can be avoided through non-Markovian coupling schemes.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Dynamics and Thermodynamics of the Low-Temperature Strongly Interacting Bose Gas

We measure the zero-temperature equation of state of a homogeneous Bose gas of $^7$Li atoms by analyzing the \emph{in-situ} density distributions of trapped samples. For increasing repulsive interactions our data shows a clear departure from mean-field theory and provides a quantitative test of the many-body corrections first predicted in 1957 by Lee, Huang and Yang. We further probe the dynamic response of the Bose gas to a varying interaction strength and compare it to simple theoretical models. We deduce a lower bound for the value of the universal constant $ξ>0.44(8)$ that would characterize the universal Bose gas at the unitary limit.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Event-driven Monte Carlo algorithm for general potentials

We extend the event-chain Monte Carlo algorithm from hard-sphere interactions to the micro-canonical ensemble (constant potential energy) for general potentials. This event-driven Monte Carlo algorithm is non-local, rejection-free, and allows for the breaking of detailed balance. The algorithm uses a discretized potential, but its running speed is asymptotically independent of the discretization. We implement the algorithm for the cut-off linear potential, and discuss its possible implementation directly in the continuum limit.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Two-step melting in two dimensions: First-order liquid-hexatic transition

Melting in two spatial dimensions, as realized in thin films or at interfaces, represents one of the most fascinating phase transitions in nature, but it remains poorly understood. Even for the fundamental hard-disk model, the melting mechanism has not been agreed on after fifty years of studies. A recent Monte Carlo algorithm allows us to thermalize systems large enough to access the thermodynamic regime. We show that melting in hard disks proceeds in two steps with a liquid phase, a hexatic phase, and a solid. The hexatic-solid transition is continuous while, surprisingly, the liquid-hexatic transition is of first-order. This melting scenario solves one of the fundamental statistical-physics models, which is at the root of a large body of theoretical, computational and experimental research.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Event-chain Monte Carlo algorithms for hard-sphere systems

In this paper we present the event-chain algorithms, which are fast Markov-chain Monte Carlo methods for hard spheres and related systems. In a single move of these rejection-free methods, an arbitrarily long chain of particles is displaced, and long-range coherent motion can be induced. Numerical simulations show that event-chain algorithms clearly outperform the conventional Metropolis method. Irreversible versions of the algorithms, which violate detailed balance, improve the speed of the method even further. We also compare our method with a recent implementations of the molecular-dynamics algorithm.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The equilibrium state of a trapped two-dimensional Bose gas

We study experimentally and numerically the equilibrium density profiles of a trapped two-dimensional $^{87}$Rb Bose gas, and investigate the equation of state of the homogeneous system using the local density approximation. We find a clear discrepancy between in-situ measurements and Quantum Monte Carlo simulations, which we attribute to a non-linear variation of the optical density of the atomic cloud with its spatial density. However, good agreement between experiment and theory is recovered for the density profiles measured after time-of-flight, taking advantage of their self-similarity in a two-dimensional expansion.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Convergence and coupling for spin glasses and hard spheres

We discuss convergence and coupling of Markov chains, and present general relations between the transfer matrices describing these two processes. We then analyze a recently developed local-patch algorithm, which computes rigorous upper bound for the coupling time of a Markov chain for non-trivial statistical-mechanics models. Using the coupling from the past protocol, this allows one to exactly sample the underlying equilibrium distribution. For spin glasses in two and three spatial dimensions, the local-patch algorithm works at lower temperatures than previous exact-sampling methods. We discuss variants of the algorithm which might allow one to reach, in three dimensions, the spin-glass transition temperature. The algorithm can be adapted to hard-sphere models. For two-dimensional hard disks, the algorithm allows us to draw exact samples at higher densities than previously possible.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Universal correlations and coherence in quasi-two-dimensional trapped Bose gases

We study the quasi-two-dimensional Bose gas in harmonic traps at temperatures above the Kosterlitz-Thouless transition, where the gas is in the normal phase. We show that mean-field theory takes into account the dominant interaction effects for experimentally relevant trap geometries. Comparing with Quantum Monte Carlo calculations, we quantify the onset of the fluctuation regime, where correlations beyond mean-field become important. Although the density profile depends on the microscopic parameters of the system, we show that the correlation density (the difference between the exact and the mean-field density) is accurately described by a universal expression, obtained from classical-field calculations of the homogeneous strictly two-dimensional gas. Deviations from universality, due to the finite value of the interaction or to the trap geometry, are shown to be small for current experiments. We further study coherence and pair correlations on a microscopic scale. Finite-size effects in the off-diagonal density matrix allows us to characterize the cross-over from Kosterlitz-Thouless to Bose-Einstein behavior for small particle numbers. Bose-Einstein condensation occurs below a characteristic number of particles which rapidly diverges with vanishing interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Geometry of Gaussian signals

We consider Gaussian signals, i.e. random functions $u(t)$ ($t/L \in [0,1]$) with independent Gaussian Fourier modes of variance $\sim 1/q^α$, and compute their statistical properties in small windows $[x, x+δ]$. We determine moments of the probability distribution of the mean square width of $u(t)$ in powers of the window size $δ$. We show that the moments, in the small-window limit $δ\ll 1$, become universal, whereas they strongly depend on the boundary conditions of $u(t)$ for larger $δ$. For $α> 3$, the probability distribution is computed in the small-window limit and shown to be independent of $α$.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

Coulomb and Liquid Dimer Models in Three Dimensions

We study classical hard-core dimer models on three-dimensional lattices using analytical approaches and Monte Carlo simulations. On the bipartite cubic lattice, a local gauge field generalization of the height representation used on the square lattice predicts that the dimers are in a critical Coulomb phase with algebraic, dipolar, correlations, in excellent agreement with our large-scale Monte Carlo simulations. The non-bipartite FCC and Fisher lattices lack such a representation, and we find that these models have both confined and exponentially deconfined but no critical phases. We conjecture that extended critical phases are realized only on bipartite lattices, even in higher dimensions.

0 citations0 reviewsNo code linkedOpen access
preprint2002arXiv

Pocket Monte Carlo algorithm for classical doped dimer models

We study the correlations of classical hardcore dimer models doped with monomers by Monte Carlo simulation. We introduce an efficient cluster algorithm, which is applicable in any dimension, for different lattices and arbitrary doping. We use this algorithm for the dimer model on the square lattice, where a finite density of monomers destroys the critical confinement of the two-monomer problem. The monomers form a two-component plasma located in its high-temperature phase, with the Coulomb interaction screened at finite densities. On the triangular lattice, a single pair of monomers is not confined. The monomer correlations are extremely short-ranged and hardly change with doping.

0 citations0 reviewsNo code linkedOpen access
preprint1999arXiv

Yang-Mills Integrals

SU(N) Yang-Mills integrals form a new class of matrix models which, in their maximally supersymmetric version, are relevant to recent non-perturbative definitions of 10-dimensional IIB superstring theory and 11-dimensional M-theory. We demonstrate how Monte Carlo methods may be used to establish important properties of these models. In particular we consider the partition functions as well as the matrix eigenvalue distributions. For the latter we derive a number of new exact results for SU(2). We also report preliminary computations of Wilson loops. (Based on talk presented by M. Staudacher at Strings '99, Potsdam, July 19-24 1999)

0 citations0 reviewsNo code linkedOpen access