1 |
|
%\documentclass[prb,aps,twocolumn,tabularx]{revtex4} |
2 |
|
%\documentclass[aps,prb,preprint]{revtex4} |
3 |
< |
\documentclass[10pt]{article} |
4 |
< |
%\usepackage{endfloat} |
3 |
> |
\documentclass[11pt]{article} |
4 |
> |
\usepackage{endfloat} |
5 |
|
\usepackage{amsmath,bm} |
6 |
|
\usepackage{amssymb} |
7 |
|
\usepackage{epsf} |
27 |
|
|
28 |
|
\title{Is the Ewald summation still necessary? \\ |
29 |
|
Pairwise alternatives to the accepted standard for \\ |
30 |
< |
long-range electrostatics} |
30 |
> |
long-range electrostatics in molecular simulations} |
31 |
|
|
32 |
|
\author{Christopher J. Fennell and J. Daniel Gezelter\footnote{Corresponding author. \ Electronic mail: |
33 |
|
gezelter@nd.edu} \\ |
38 |
|
\date{\today} |
39 |
|
|
40 |
|
\maketitle |
41 |
< |
%\doublespacing |
41 |
> |
\doublespacing |
42 |
|
|
43 |
|
\begin{abstract} |
44 |
|
We investigate pairwise electrostatic interaction methods and show |
531 |
|
|
532 |
|
\begin{figure} |
533 |
|
\centering |
534 |
< |
\includegraphics[width = 3.25in]{./dualLinear.pdf} |
534 |
> |
\includegraphics[width = \linewidth]{./dualLinear.pdf} |
535 |
|
\caption{Example least squares regressions of the configuration energy |
536 |
|
differences for SPC/E water systems. The upper plot shows a data set |
537 |
|
with a poor correlation coefficient ($R^2$), while the lower plot |
717 |
|
|
718 |
|
\begin{figure} |
719 |
|
\centering |
720 |
< |
\includegraphics[width=3.25in]{./delEplot.pdf} |
720 |
> |
\includegraphics[width=5.5in]{./delEplot.pdf} |
721 |
|
\caption{Statistical analysis of the quality of configurational energy |
722 |
|
differences for a given electrostatic method compared with the |
723 |
|
reference Ewald sum. Results with a value equal to 1 (dashed line) |
777 |
|
|
778 |
|
\begin{figure} |
779 |
|
\centering |
780 |
< |
\includegraphics[width=3.25in]{./frcMagplot.pdf} |
780 |
> |
\includegraphics[width=5.5in]{./frcMagplot.pdf} |
781 |
|
\caption{Statistical analysis of the quality of the force vector |
782 |
|
magnitudes for a given electrostatic method compared with the |
783 |
|
reference Ewald sum. Results with a value equal to 1 (dashed line) |
816 |
|
|
817 |
|
\begin{figure} |
818 |
|
\centering |
819 |
< |
\includegraphics[width=3.25in]{./trqMagplot.pdf} |
819 |
> |
\includegraphics[width=5.5in]{./trqMagplot.pdf} |
820 |
|
\caption{Statistical analysis of the quality of the torque vector |
821 |
|
magnitudes for a given electrostatic method compared with the |
822 |
|
reference Ewald sum. Results with a value equal to 1 (dashed line) |
856 |
|
|
857 |
|
\begin{figure} |
858 |
|
\centering |
859 |
< |
\includegraphics[width=3.25in]{./frcTrqAngplot.pdf} |
859 |
> |
\includegraphics[width=5.5in]{./frcTrqAngplot.pdf} |
860 |
|
\caption{Statistical analysis of the width of the angular distribution |
861 |
|
that the force and torque vectors from a given electrostatic method |
862 |
|
make with their counterparts obtained using the reference Ewald sum. |
991 |
|
|
992 |
|
\begin{figure} |
993 |
|
\centering |
994 |
< |
\includegraphics[width = 3.25in]{./vCorrPlot.pdf} |
994 |
> |
\includegraphics[width = \linewidth]{./vCorrPlot.pdf} |
995 |
|
\caption{Velocity autocorrelation functions of NaCl crystals at |
996 |
|
1000 K using {\sc spme}, {\sc sf} ($\alpha$ = 0.0, 0.1, \& 0.2), and {\sc |
997 |
|
sp} ($\alpha$ = 0.2). The inset is a magnification of the area around |
1031 |
|
|
1032 |
|
\begin{figure} |
1033 |
|
\centering |
1034 |
< |
\includegraphics[width = 3.25in]{./spectraSquare.pdf} |
1034 |
> |
\includegraphics[width = \linewidth]{./spectraSquare.pdf} |
1035 |
|
\caption{Power spectra obtained from the velocity auto-correlation |
1036 |
|
functions of NaCl crystals at 1000 K while using {\sc spme}, {\sc sf} |
1037 |
|
($\alpha$ = 0, 0.1, \& 0.2), and {\sc sp} ($\alpha$ = 0.2). The inset |
1070 |
|
|
1071 |
|
\begin{figure} |
1072 |
|
\centering |
1073 |
< |
\includegraphics[width = 3.25in]{./increasedDamping.pdf} |
1073 |
> |
\includegraphics[width = \linewidth]{./increasedDamping.pdf} |
1074 |
|
\caption{Effect of damping on the two lowest-frequency phonon modes in |
1075 |
|
the NaCl crystal at 1000~K. The undamped shifted force ({\sc sf}) |
1076 |
|
method is off by less than 10 cm$^{-1}$, and increasing the |