Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “MD simulation”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 703 records · Page 39Linked to original sources
PubMed2002

An improved OPLS-AA force field for carbohydrates.

This work describes an improved version of the original OPLS-all atom (OPLS-AA) force field for carbohydrates (Damm et al., J Comp Chem 1997, 18, 1955). The improvement is achieved by applying additional scaling factors for the electrostatic interactions between 1,5- and 1,6-interactions. This new model is tested first for improving the conformational energetics of 1,2-ethanediol, the smallest polyol. With a 1,5-scaling factor of 1.25 the force field calculated relative energies are in excellent agreement with the ab initio-derived data. Applying the new 1,5-scaling makes it also necessary to use a 1,6-scaling factor for the interactions between the C4 and C6 atoms in hexopyranoses. After torsional parameter fitting, this improves the conformational energetics in comparison to the OPLS-AA force field. The set of hexopyranoses included in the torsional parameter derivation consists of the two anomers of D-glucose, D-mannose, and D-galactose, as well as of the methyl-pyranosides of D-glucose, D-mannose. Rotational profiles for the rotation of the exocyclic group and of different hydroxyl groups are also compared for the two force fields and at the ab initio level of theory. The new force field reduces the overly high barriers calculated using the OPLS-AA force field. This leads to better sampling, which was shown to produce more realistic conformational behavior for hexopyranoses in liquid simulation. From 10-ns molecular dynamics (MD) simulations of alpha-D-glucose and alpha-D-galactose the ratios for the three different conformations of the hydroxymethylene group and the average (3)J(H,H) coupling constants are derived and compared to experimental values. The results obtained for OPLS-AA-SEI force field are in good agreement with experiment whereas the properties derived for the OPLS-AA force field suffer from sampling problems. The undertaken investigations show that the newly derived OPLS-AA-SEI force field will allow simulating larger carbohydrates or polysaccharides with improved sampling of the hydroxyl groups.

Carbohydrate Conformation↗
PubMed1996

Conformational studies on the selectin and natural killer cell receptor ligands sulfo- and sialyl-lacto-N-fucopentaoses (SuLNFPII and SLNFPII) using NMR spectroscopy and molecular dynamics simulations. Comparisons with the nonacidic parent molecule LNFPII.

This investigation is focused on the conformational behavior of the blood group Lewisa (Le(a)-active pentasaccharide lacto-N-fucopentaose II (LNFPII) and its sulfated and sialylated analogs, SuLNFPII and SLNFPII. The latter two are more potent oligosaccharide ligands for the animal lectins, E- and L-selectin, and the natural killer cell receptor, NKR-P1, than are the shorter chain analogs based on the trisaccharide Le(a) domain. We report here that the three oligosaccharides based on the fucopentasaccharide have very similar average solution conformations as determined from NMR spectroscopical parameters, in particular 13C chemical shift differences. From restrained simulated annealing and restrained molecular dynamics (MD) simulations performed in order to determine the most probable conformational distributions around the glycosidic linkages we derive models for these oligosaccharides that are in good agreement with experimental parameters, such as rotating-frame Overhauser effects (ROE's) and long-range 1H,13C coupling constants across the glycosidic linkages. In these model structures the Le(a) domain at the non-reducing end of the longer chain oligosaccharides approximates the same rigid structure as in the shorter analogs. The Gal beta 1-4Glc linkage at the reducing end is also rather rigid, showing only little more flexibility than the Le(a) domain. However, the NeuAc alpha 2-3Gal linkage in SLNFPII, and the GlcNAc beta 1-3Gal linkage in all three oligosaccharides are flexible, in each case fluctuating mainly between two minimum energy structures: (phi = -81 degrees, psi = 8 degrees) and (phi = -160 degrees, psi = -20 degrees) for the NeuAc alpha 2-3Gal linkage, as reported previously for the isomeric sequence 3'-sialyl Le(x), and (phi = -25 degrees, psi = -26 degrees) and (phi = 20 degrees, psi = 24 degrees) for the GlcNAc beta 1-3Gal linkage. The flexibility of the latter linkage may allow the lactosyl domain at the reducing end to fit with little strain into extended carbohydrate binding sites on the recognition proteins, and, for the purposes of drug designs, it will be important to establish which conformational distribution is assumed for the GlcNAc beta 1-3Gal linkage in these longer chain oligosaccharides in the bound state.

Animals↗
PubMed2001

Modeling platinum group metal complexes in aqueous solution.

We construct force fields suited for the study of three platinum group metals (PGM) as chloranions in aqueous solution from quantum chemical computations and report experimental data. Density functional theory (DFT) using the local density approximation (LDA), as well as extended basis sets that incorporate relativistic corrections for the transition metal atoms, has been used to obtain equilibrium geometries, harmonic vibrational frequencies, and atomic charges for the complexes. We found that DFT calculations of [PtCl(6)](2-).3H(2)O, [PdCl(4)](2-).2H(2)O, and [RhCl(6)](3-).3H(2)O water clusters compared well with molecular mechanics (MM) calculations using the specific force field developed here. The force field performed equally well in condensed phase simulations. A 500 ps molecular dynamics (MD) simulation of [PtCl(6)](2-) in water was used to study the structure of the solvation shell around the anion. The resulting data were compared to an experimental radial distribution function derived from X-ray diffraction experiments. We found the calculated pair correlation functions (PCF) for hexachloroplatinate to be in good agreement with experiment and were able to use the simulation results to identify and resolve two water-anion peaks in the experimental spectrum.

Journal Article↗
PubMed2000

The role of side chains in the interaction of new antitumor pyrimidoacridinetriones with DNA: molecular dynamics simulations.

Pyrimidoacridinetriones (PATs) are a new group of highly active antitumor compounds. It seems reasonable to assume that, like for some other acridine derivatives, intercalation into DNA is a necessary, however not a sufficient condition for antitumor activity of these compounds. Rational design of new compounds of this chemotype requires knowledge about the structure of the intercalation complex, as well as about interactions responsible for its stability. Computer simulation techniques such as molecular dynamics (MD) may provide valuable information about these problems. The results of MD simulations performed for three rationally selected PATs are presented in this paper. The compounds differ in the number and position of side chains. Each of the compounds was simulated in two systems: i) in water, and ii) in the intercalation complex with the dodecamer duplex d(GCGCGCGCGCGC)2. The orientation of the side chain in relation to the ring system is determined by the position of its attachment. Orientation of the ring system inside the intercalation cavity depends on the number and position of side chain(s). The conformations of the side chain(s) of all PATs studied in the intercalation complex were found to be very similar to those observed in water.

Acridines↗
PubMed1996

Characterization of lipid membrane dynamics by simulation: 3. Probing molecular transport across the phospholipid bilayer.

PURPOSE: The goal of this study is to elucidate the role of the motions of the hydrocarbon chains of a phospholipid bilayer in penetrant diffusion. Penetrant size, as well as its position in the hydrocarbon core of the lipid bilayer, has also been explored regarding impact on the diffusion rate in a phospholipid bilayer. METHOD: Molecular dynamics, MD, simulations were carried out on a model dimyristoyl phosphatidylcholine (DMPC) membrane bilayer with and without methanol and propanol as penetrants. The MD trajectories were analyzed in terms of estimating time and space properties. RESULTS: These simulations show that torsion angle kink shifts in the hydrocarbon chains of phospholipids are natural occurrences in a bilayer assembly. The diffusion coefficients of methanol and propanol in a DMPC lipid bilayer, as calculated from the MD simulations, agree with experimental measurements. Both methanol and propanol show different diffusion rates in different regions of the hydrocarbon chain matrix of the lipid bilayer. Solute size has more impact on diffusion rate in the bilayer regions with high torsion angle order parameters, as compared to the regions with low torsion angle order parameters. CONCLUSIONS: The simulated transport behavior suggests that a kink shift diffusion mechanism is more likely to occur in regions with high torsion angle order parameters, and a free volume transport mechanism is more likely operative in the region with low torsion angle order parameters, mainly the center core of the bilayer. A three zone diffusion model is proposed for transport of a penetrant across a bilayer.

1-Propanol↗
PubMed2005

Phase and interface behaviors in type-I and type-V Lennard-Jones mixtures: theory and simulations.

Density gradient theory (DGT) and molecular-dynamics (MD) simulations have been used to predict subcritical phase and interface behaviors in type-I and type-V equal-size Lennard-Jones mixtures. Type-I mixtures exhibit a continuum critical line connecting their pure critical components, which implies that their subcritical phase equilibria are gas liquid. Type-V mixtures are characterized by two critical lines and a heteroazeotropic line. One of the two critical lines begins at the more volatile pure component critical point up to an upper critical end point and the other one comes from the less volatile pure component critical point ending at a lower critical end point. The heteroazeotropic line connects both critical end points and is characterized by gas-liquid-liquid equilibria. Therefore, subcritical states of this type exhibit gas-liquid and gas-liquid-liquid equilibria. In order to obtain a correct characterization of the phase and interface behaviors of these types of mixtures and to directly compare DGT and MD results, the global phase diagram of equal-size Lennard-Jones mixtures has been used to define the molecular parameters of these mixtures. According to our results, DGT and MD are two complementary methodologies able to obtain a complete and simultaneous prediction of phase equilibria and their interfacial properties. For the type of mixtures analyzed here, both approaches have shown excellent agreement in their phase equilibrium and interface properties in the full concentration range.

Journal Article↗
PubMed2003

Analysis of the effect of electrostatic energy truncation in molecular dynamics simulations of immunoglobulin G light chain dimer.

Molecular dynamics (MD) simulations of immunoglobulin G (IgG) light chain dimer using particle mesh Ewald (PME) and cutoff methods of treating electrostatic interactions were performed. The results indicate that structural parameters (RMSD, radius of gyration, solvent accessible surface) are very similar for both schemes; however, PME simulation shows increased mobility of side chains. This leads to larger fluctuations in the distance between the monomers in the dimer molecule, and, as a consequence, results in decreased number of interactions across the dimer interface. The wall clock time of the simulations was also compared. It was shown that the PME method is approximately 30% faster than the cutoff method for the system studied on a single processor.

Computer Simulation↗
PubMed2005

Supercooled water in PVA matrixes. II. A molecular dynamics simulation study and comparison with QENS results.

Molecular dynamics (MD) simulations were carried out to elucidate the dynamic behavior of water confined in poly(vinyl alcohol), PVA, hydrogels. Model topology is supported by experimental network parameters, and simulation results are compared to an incoherent quasielastic neutron scattering (QENS) investigation carried out on PVA hydrogels. From the QENS dynamic scattering law (part I), a random jump model was adopted for the description of water diffusion to extract a microscopic diffusion coefficient and a residence time between two "jumps". In the present work, consistently with this framework, water diffusion parameters as diffusion coefficients and residence times have been evaluated using the mean square displacement of water in a time window of 10 ps and the time autocorrelation function of water hydrogen bonds. The calculated parameters are in good agreement with the experimental ones, giving confidence to this approach. Further developments are in progress to take into account a more realistic description of hydrogel structure in the molecular dynamics simulations.

Journal Article↗
PubMed2006

OD vibrations and hydration structure in an Al3+ (aq) solution from a Car-Parrinello molecular-dynamics simulation.

The optimized geometry, energetics, and vibrational properties of Al(D2O)n3+ clusters, with n = 1,2,4, and 6, have been studied using plane waves, different local basis sets, different methodologies [density-functional theory, MP2, CCSD(T)], and different functionals (BLYP, PBE). Moreover, Car-Parrinello molecular-dynamics (MD) simulations using the BLYP functional, plane waves, and the Vanderbilt ultrasoft pseudopotentials have been performed for an aqueous Al3+ solution with 1 ion and 32 D2O molecules in a periodic box at room temperature, studied for 10 ps. The cluster calculations were performed to pinpoint possible shortcomings of the electronic structure description used in the Car-Parinello MD (CPMD) simulation. For the clusters, the hydration structure and interaction energies calculated with the 'BLYP/plane-wave' approach agree well with high-level ab initio methods but the exchange-correlation functional introduces errors in the OD stretching frequencies (both in the absolute values and in the ion-induced shifts). For the aqueous solution, the CPMD simulation yields structural properties in good agreement with experimental data. The CPMD-simulated OD stretching vibrational band for the first-shell water molecules around Al3+ is strongly downshifted by the influence of the ion and is compared with experimental data from the literature. To make such a comparison meaningful, the influences of a number of systematic effects have been addressed, such as the exchange-correlation functional, the fictitious electron mass, anharmonicity effects, and the small box size in the simulation. Each of these factors (except the last one) is found to affect the OD frequency by 100 cm(-1) or more. The final "corrected" frequencies agree with experiment within approximately 30 cm(-1) for bulk water but are too little downshifted for the first-shell Al3+ (aq) water molecules (by approximately 200 cm(-1)).

Journal Article↗
PubMed2004

Unveiling functional protein motions with picosecond x-ray crystallography and molecular dynamics simulations.

A joint analysis of all-atom molecular dynamics (MD) calculations and picosecond time-resolved x-ray structures was performed to gain single-molecule insights into mechanisms of protein function. Ensemble-averaged MD simulations of the L29F mutant of myoglobin after ligand dissociation reproduce the direction, amplitude, and time scales of crystallographically determined structural changes. This close agreement with experiments at comparable resolution in space and time validates the individual MD trajectories. From 1,700 single-molecule trajectories, we identified and structurally characterized a conformational switch that directs dissociated ligands to one of two nearby protein cavities. Subsequent ligand migration proceeds through a network of transiently interconnected internal cavities, with passage between them involving correlated protein-ligand motions. The simulations also suggest how picosecond protein motions modulate the functional dissociation of oxygen and suppress the geminate recombination of toxic carbon monoxide.

Carbon Monoxide↗
PubMed2025

Computational Insights on the Assembly of the Dengue Virus Membrane-Capsid-RNA Complex.

Dengue virus, an arbovirus from the genus Flavivirus in the family Flaviviridae, forms a nucleocapsid structure through interactions between its genome and multiple copies of the capsid protein. Experimental studies have confirmed the interaction between the viral capsid protein and lipid droplets, indicating a protein-lipid interaction. Cryo-EM studies show that in immature viruses, the nucleocapsid is located close to the viral membrane. This study uses multiple MD simulations to explore the orientation of the capsid protein relative to the lipid membrane, focusing on how the protein's hydrophobic pocket interacts with the membrane. We also investigated the interaction between the capsid protein and RNA, considering the effects of sequence length and identity. Finally, we construct a model of the lipid-protein-RNA complex, demonstrating that the capsid protein's hydrophobic pocket interacts with the membrane, while the positively charged H4 helix interacts with the negatively charged RNA. This research may identify crucial interactions for immature virus particle formation and provide insights for future therapeutic interventions.

Dengue Virus↗
PubMed2006

Molecular dynamics simulations of sarcin-ricin rRNA motif.

Explicit solvent molecular dynamics (MD) simulations were carried out for sarcin-ricin domain (SRD) motifs from 23S (Escherichia coli) and 28S (rat) rRNAs. The SRD motif consists of GAGA tetraloop, G-bulged cross-strand A-stack, flexible region and duplex part. Detailed analysis of the overall dynamics, base pairing, hydration, cation binding and other SRD features is presented. The SRD is surprisingly static in multiple 25 ns long simulations and lacks any non-local motions, with root mean square deviation (r.m.s.d.) values between averaged MD and high-resolution X-ray structures of 1-1.4 A. Modest dynamics is observed in the tetraloop, namely, rotation of adenine in its apex and subtle reversible shift of the tetraloop with respect to the adjacent base pair. The deformed flexible region in low-resolution rat X-ray structure is repaired by simulations. The simulations reveal few backbone flips, which do not affect positions of bases and do not indicate a force field imbalance. Non-Watson-Crick base pairs are rigid and mediated by long-residency water molecules while there are several modest cation-binding sites around SRD. In summary, SRD is an unusually stiff rRNA building block. Its intrinsic structural and dynamical signatures seen in simulations are strikingly distinct from other rRNA motifs such as Loop E and Kink-turns.

Animals↗
PubMed2006

Equation of state of nitrogen (N2) at high pressures and high temperatures: molecular dynamics simulation.

Nitrogen equation of state at pressures up to 30 GPa (300 kbars) and temperatures above 800 K was studied by molecular dynamics (MD) simulations. The dynamics of the N(2) molecules is treated in hard rotor approximation, i.e., it accounts both translational and rotational degrees of freedom. The rotational motion of the N(2) molecule is treated assuming constant moment of inertia of the nitrogen molecule. The new MD program fully accounts anisotropic molecular nitrogen interaction. The N(2)-N(2) interaction potential has been derived by van der Avoird et al. [J. Chem. Phys. 84, 1629 (1986)] using the results of high precision Hartree-Fock ab initio quantum mechanical calculations. The potential, fully accounts rotational symmetry of the N(2)-N(2) system, by employing 6-j Wigner symbols, i.e., preserving full rotational symmetry of the system. Various numerical algorithms were tested, in order to achieve the energy preservation during the simulation. It has been demonstrated that the standard Verlet algorithm was not preserving the energy for the standard MD time step, equal to 5x10(-16) s. Runge-Kutta fourth order method was able to preserve the energy within 10(-4) relative error, but it requires calculation of the force four times for each time step and therefore it is highly inefficient. A predictor-corrector method of the fifth order (PC5) was found to be efficient and precise and was therefore adopted for the simulation of the molecular nitrogen properties at high pressure. Singer and Fincham algorithms were tested and were found to be as precise as PC5 algorithm and they were also used in the simulation of the equation of state. Results of MD simulations are in very good agreement with the experimental data on nitrogen equation of state at pressures below 1 GPa (10 kbars). For higher pressures, up to 30 GPa (300 kbars), i.e., close to molecular nitrogen stability limit, determined by Nellis et al. [Phys. Rev. Lett. 85, 1262 (1984)], the obtained numerical results provide new data of the experimentally unexplored region. These data were formulated in the analytical form of pressure-density-temperature equation of state.

Journal Article↗
PubMed2005

Onset of patterns in an oscillated granular layer: continuum and molecular dynamics simulations.

We study the onset of patterns in vertically oscillated layers of frictionless dissipative particles. Using both numerical solutions of continuum equations to Navier-Stokes order and molecular dynamics (MD) simulations, we find that standing waves form stripe patterns above a critical acceleration of the cell. Changing the frequency of oscillation of the cell changes the wavelength of the resulting pattern; MD and continuum simulations both yield wavelengths in accord with previous experimental results. The value of the critical acceleration for ordered standing waves is approximately 10% higher in molecular dynamics simulations than in the continuum simulations, and the amplitude of the waves differs significantly between the models. The delay in the onset of order in molecular dynamics simulations and the amplitude of noise below this onset are consistent with the presence of fluctuations which are absent in the continuum theory. The strength of the noise obtained by fit to Swift-Hohenberg theory is orders of magnitude larger than the thermal noise in fluid convection experiments, and is comparable to the noise found in experiments with oscillated granular layers and in recent fluid experiments on fluids near the critical point. Good agreement is found between the mean field value of onset from the Swift-Hohenberg fit and the onset in continuum simulations. Patterns are compared in cells oscillated at two different frequencies in MD; the layer with larger wavelength patterns has less noise than the layer with smaller wavelength patterns.

Journal Article↗
PubMed2003

Membrane structure of the human immunodeficiency virus gp41 fusion peptide by molecular dynamics simulation. II. The glycine mutants.

In this work, molecular dynamics (MD) simulation of the interaction of three mutants, G3V, G5V and G10V, of the human immunodeficiency virus (HIV) gp41 16-residue fusion peptide (FP) with an explicit palmitoyloleoylphosphatidyl-ethanolamine (POPE) lipid bilayer was performed. The goals of this work are to study the correlation of the fusogenic activity of the FPs with the mode of their interaction with the bilayer and to examine the roles of the many glycine residues in the FP in the fusion process. The results of this work corroborate the main conclusion of our earlier MD work of the WT FP and several mutants with polar substitution. These two studies provide correlation between the mode of insertion and the fusogenic activity of these peptides and support the hypothesis that an oblique insertion of the fusion domain of the viral protein is required for fusogenic activity. Inactive mutants interact with the bilayer by a surface-binding mode. The results of this work, combined with the results of our earlier work, show that, while the secondary structures of the wild-type FP and its mutants do not affect the fusogenic activities, the conformational flexibility appears to be an important factor. The active WT FP and its partially active mutants, G3V and G5V, all have significant conformational transitions at one of the glycine sites. They occur at Gly(5) in FP-wt, at Gly(10) in FP-G5V and at Gly(13) in FP-G3V. Thus, a glycine site in each of these active (or partially active) FPs provides conformational flexibility. On the other hand, the inactive mutants FP-G10V, FP-L9R and FP-V2E do not have any conformational transitions except at either terminus and thus possess no conformational flexibility. Thus, the results of this work support the suggestion that the role of glycine residues in the fusion domain is to provide the necessary conformational flexibility for fusion activity. The glycines also form a "glycine strip" in the FP that locates on one (the less hydrophobic) face of the helix (the "sided helix"). However, whether this "glycine strip" is disrupted or not does not seem to correlate with the retention of fusogenic activities. Finally, although the FLGFL (8-12) motif is absolutely conserved in the HIV fusion domain, a well-structured motif stabilized by hydrogen bonding does not appear to be required for activity. In fact, hydrogen bonding in this motif was found to be missing in FP-G3V and FP-G5V. Both of these mutants are partially active.

Glycine↗
PubMed1999

An alamethicin channel in a lipid bilayer: molecular dynamics simulations.

We present the results of 2-ns molecular dynamics (MD) simulations of a hexameric bundle of Alm helices in a 1-palmitoyl-2-oleoylphosphatidylcholine bilayer. These simulations explore the dynamic properties of a model of a helix bundle channel in a complete phospholipid bilayer in an aqueous environment. We explore the stability and conformational dynamics of the bundle in a phospholipid bilayer. We also investigate the effect on bundle stability of the ionization state of the ring of Glu18 side chains. If all of the Glu18 side chains are ionised, the bundle is unstable; if none of the Glu18 side chains are ionized, the bundle is stable. pKA calculations suggest that either zero or one ionized Glu18 is present at neutral pH, correlating with the stable form of the helix bundle. The structural and dynamic properties of water in this model channel were examined. As in earlier in vacuo simulations (Breed et al., 1996 .Biophys. J. 70:1643-1661), the dipole moments of water molecules within the pore were aligned antiparallel to the helix dipoles. This contributes to the stability of the helix bundle.

Alamethicin↗
PubMed2004

Estimating entropies from molecular dynamics simulations.

While the determination of free-energy differences by MD simulation has become a standard procedure for which many techniques have been developed, total entropies and entropy differences are still hardly ever computed. An overview of techniques to determine entropy differences is given, and the accuracy and convergence behavior of five methods based on thermodynamic integration and perturbation techniques was evaluated using liquid water as a test system. Reasonably accurate entropy differences are obtained through thermodynamic integration in which many copies of a solute are desolvated. When only one solute molecule is involved, only two methods seem to yield useful results, the calculation of solute-solvent entropy through thermodynamic integration, and the calculation of solvation entropy through the temperature derivative of the corresponding free-energy difference. One-step perturbation methods seem unsuitable to obtain entropy estimates.

Journal Article↗
PubMed2005

Molecular dynamics in the isothermal-isobaric ensemble: the requirement of a "shell" molecule. I. Theory and phase-space analysis.

Current constant pressure molecular-dynamics (MD) algorithms are not consistent with the recent reformulation of the isothermal-isobaric (NpT) ensemble. The NpT ensemble partition function requires the use of a "shell" molecule to identify uniquely the volume of the system, thereby avoiding the redundant counting of configurations [e.g., G. J. M. Koper and H. Reiss, J. Phys. Chem. 100, 422 (1996); D. S. Corti, Phys. Rev. E, 64, 016128 (2001)]. So far, only the NpT Monte Carlo method has been updated to allow the system volume to be defined by a shell particle [D. S. Corti, Mol. Phys. 100, 1887 (2002)]. A shell particle has yet to be incorporated into MD simulations. The proper modification of the NpT MD algorithm is therefore the subject of this paper. Unlike Andersen's method [H. C. Andersen, J. Chem. Phys. 72, 2384 (1980)] where a piston of unknown mass serves to control the response time of volume fluctuations, the newly proposed equations of motion impose a constant external pressure via the introduction of a shell particle of known mass. Hence, the system itself sets the time scales for pressure and volume fluctuations. The new algorithm is subject to a number of fundamentally rigorous tests to ensure that the equations of motion sample phase space correctly. We also show that the Hoover NpT algorithm [W. G. Hoover, Phys. Rev. A. 31, 1695 (1985); 34, 2499 (1986)] does sample phase correctly, but only when periodic boundary conditions are employed.

Journal Article↗