ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/group/trunk/tengDissertation/dissertation.bib
(Generate patch)

Comparing trunk/tengDissertation/dissertation.bib (file contents):
Revision 2801 by tim, Tue Jun 6 14:56:36 2006 UTC vs.
Revision 2839 by tim, Fri Jun 9 02:41:58 2006 UTC

# Line 3 | Line 3 | Encoding: GBK
3  
4   @ARTICLE{Torre2003,
5    author = {J. G. {de la Torre} and H. E. Sanchez and A. Ortega and J. G. Hernandez
6 <    and M. X. Fernandes and F. G. Diaz and M. C. L. Martinez},
6 >        and M. X. Fernandes and F. G. Diaz and M. C. L. Martinez},
7    title = {Calculation of the solution properties of flexible macromolecules:
8 <    methods and applications},
8 >        methods and applications},
9    journal = {European Biophysics Journal with Biophysics Letters},
10    year = {2003},
11    volume = {32},
# Line 13 | Line 13 | Encoding: GBK
13    number = {5},
14    month = {Aug},
15    abstract = {While the prediction of hydrodynamic properties of rigid particles
16 <    is nowadays feasible using simple and efficient computer programs,
17 <    the calculation of such properties and, in general, the dynamic
18 <    behavior of flexible macromolecules has not reached a similar situation.
19 <    Although the theories are available, usually the computational work
20 <    is done using solutions specific for each problem. We intend to
21 <    develop computer programs that would greatly facilitate the task
22 <    of predicting solution behavior of flexible macromolecules. In this
23 <    paper, we first present an overview of the two approaches that are
24 <    most practical: the Monte Carlo rigid-body treatment, and the Brownian
25 <    dynamics simulation technique. The Monte Carlo procedure is based
26 <    on the calculation of properties for instantaneous conformations
27 <    of the macromolecule that are regarded as if they were instantaneously
28 <    rigid. We describe how a Monte Carlo program can be interfaced to
29 <    the programs in the HYDRO suite for rigid particles, and provide
30 <    an example of such calculation, for a hypothetical particle: a protein
31 <    with two domains connected by a flexible linker. We also describe
32 <    briefly the essentials of Brownian dynamics, and propose a general
33 <    mechanical model that includes several kinds of intramolecular interactions,
34 <    such as bending, internal rotation, excluded volume effects, etc.
35 <    We provide an example of the application of this methodology to
36 <    the dynamics of a semiflexible, wormlike DNA.},
16 >        is nowadays feasible using simple and efficient computer programs,
17 >        the calculation of such properties and, in general, the dynamic
18 >        behavior of flexible macromolecules has not reached a similar situation.
19 >        Although the theories are available, usually the computational work
20 >        is done using solutions specific for each problem. We intend to
21 >        develop computer programs that would greatly facilitate the task
22 >        of predicting solution behavior of flexible macromolecules. In this
23 >        paper, we first present an overview of the two approaches that are
24 >        most practical: the Monte Carlo rigid-body treatment, and the Brownian
25 >        dynamics simulation technique. The Monte Carlo procedure is based
26 >        on the calculation of properties for instantaneous conformations
27 >        of the macromolecule that are regarded as if they were instantaneously
28 >        rigid. We describe how a Monte Carlo program can be interfaced to
29 >        the programs in the HYDRO suite for rigid particles, and provide
30 >        an example of such calculation, for a hypothetical particle: a protein
31 >        with two domains connected by a flexible linker. We also describe
32 >        briefly the essentials of Brownian dynamics, and propose a general
33 >        mechanical model that includes several kinds of intramolecular interactions,
34 >        such as bending, internal rotation, excluded volume effects, etc.
35 >        We provide an example of the application of this methodology to
36 >        the dynamics of a semiflexible, wormlike DNA.},
37    annote = {724XK Times Cited:6 Cited References Count:64},
38    issn = {0175-7571},
39    uri = {<Go to ISI>://000185513400011},
# Line 42 | Line 42 | Encoding: GBK
42   @ARTICLE{Alakent2005,
43    author = {B. Alakent and M. C. Camurdan and P. Doruker},
44    title = {Hierarchical structure of the energy landscape of proteins revisited
45 <    by time series analysis. II. Investigation of explicit solvent effects},
45 >        by time series analysis. II. Investigation of explicit solvent effects},
46    journal = {Journal of Chemical Physics},
47    year = {2005},
48    volume = {123},
# Line 50 | Line 50 | Encoding: GBK
50    number = {14},
51    month = {Oct 8},
52    abstract = {Time series analysis tools are employed on the principal modes obtained
53 <    from the C-alpha trajectories from two independent molecular-dynamics
54 <    simulations of alpha-amylase inhibitor (tendamistat). Fluctuations
55 <    inside an energy minimum (intraminimum motions), transitions between
56 <    minima (interminimum motions), and relaxations in different hierarchical
57 <    energy levels are investigated and compared with those encountered
58 <    in vacuum by using different sampling window sizes and intervals.
59 <    The low-frequency low-indexed mode relationship, established in
60 <    vacuum, is also encountered in water, which shows the reliability
61 <    of the important dynamics information offered by principal components
62 <    analysis in water. It has been shown that examining a short data
63 <    collection period (100 ps) may result in a high population of overdamped
64 <    modes, while some of the low-frequency oscillations (< 10 cm(-1))
65 <    can be captured in water by using a longer data collection period
66 <    (1200 ps). Simultaneous analysis of short and long sampling window
67 <    sizes gives the following picture of the effect of water on protein
68 <    dynamics. Water makes the protein lose its memory: future conformations
69 <    are less dependent on previous conformations due to the lowering
70 <    of energy barriers in hierarchical levels of the energy landscape.
71 <    In short-time dynamics (< 10 ps), damping factors extracted from
72 <    time series model parameters are lowered. For tendamistat, the friction
73 <    coefficient in the Langevin equation is found to be around 40-60
74 <    cm(-1) for the low-indexed modes, compatible with literature. The
75 <    fact that water has increased the friction and that on the other
76 <    hand has lubrication effect at first sight contradicts. However,
77 <    this comes about because water enhances the transitions between
78 <    minima and forces the protein to reduce its already inherent inability
79 <    to maintain oscillations observed in vacuum. Some of the frequencies
80 <    lower than 10 cm(-1) are found to be overdamped, while those higher
81 <    than 20 cm(-1) are slightly increased. As for the long-time dynamics
82 <    in water, it is found that random-walk motion is maintained for
83 <    approximately 200 ps (about five times of that in vacuum) in the
84 <    low-indexed modes, showing the lowering of energy barriers between
85 <    the higher-level minima.},
53 >        from the C-alpha trajectories from two independent molecular-dynamics
54 >        simulations of alpha-amylase inhibitor (tendamistat). Fluctuations
55 >        inside an energy minimum (intraminimum motions), transitions between
56 >        minima (interminimum motions), and relaxations in different hierarchical
57 >        energy levels are investigated and compared with those encountered
58 >        in vacuum by using different sampling window sizes and intervals.
59 >        The low-frequency low-indexed mode relationship, established in
60 >        vacuum, is also encountered in water, which shows the reliability
61 >        of the important dynamics information offered by principal components
62 >        analysis in water. It has been shown that examining a short data
63 >        collection period (100 ps) may result in a high population of overdamped
64 >        modes, while some of the low-frequency oscillations (< 10 cm(-1))
65 >        can be captured in water by using a longer data collection period
66 >        (1200 ps). Simultaneous analysis of short and long sampling window
67 >        sizes gives the following picture of the effect of water on protein
68 >        dynamics. Water makes the protein lose its memory: future conformations
69 >        are less dependent on previous conformations due to the lowering
70 >        of energy barriers in hierarchical levels of the energy landscape.
71 >        In short-time dynamics (< 10 ps), damping factors extracted from
72 >        time series model parameters are lowered. For tendamistat, the friction
73 >        coefficient in the Langevin equation is found to be around 40-60
74 >        cm(-1) for the low-indexed modes, compatible with literature. The
75 >        fact that water has increased the friction and that on the other
76 >        hand has lubrication effect at first sight contradicts. However,
77 >        this comes about because water enhances the transitions between
78 >        minima and forces the protein to reduce its already inherent inability
79 >        to maintain oscillations observed in vacuum. Some of the frequencies
80 >        lower than 10 cm(-1) are found to be overdamped, while those higher
81 >        than 20 cm(-1) are slightly increased. As for the long-time dynamics
82 >        in water, it is found that random-walk motion is maintained for
83 >        approximately 200 ps (about five times of that in vacuum) in the
84 >        low-indexed modes, showing the lowering of energy barriers between
85 >        the higher-level minima.},
86    annote = {973OH Times Cited:1 Cited References Count:33},
87    issn = {0021-9606},
88    uri = {<Go to ISI>://000232532000064},
89   }
90  
91 @BOOK{Allen1987,
92  title = {Computer Simulations of Liquids},
93  publisher = {Oxford University Press},
94  year = {1987},
95  author = {M.~P. Allen and D.~J. Tildesley},
96  address = {New York},
97 }
98
91   @ARTICLE{Allison1991,
92    author = {S. A. Allison},
93    title = {A Brownian Dynamics Algorithm for Arbitrary Rigid Bodies - Application
94 <    to Polarized Dynamic Light-Scattering},
94 >        to Polarized Dynamic Light-Scattering},
95    journal = {Macromolecules},
96    year = {1991},
97    volume = {24},
# Line 107 | Line 99 | Encoding: GBK
99    number = {2},
100    month = {Jan 21},
101    abstract = {A Brownian dynamics algorithm is developed to simulate dynamics experiments
102 <    of rigid macromolecules. It is applied to polarized dynamic light
103 <    scattering from rodlike sturctures and from a model of a DNA fragment
104 <    (762 base pairs). A number of rod cases are examined in which the
105 <    translational anisotropy is increased form zero to a large value.
106 <    Simulated first cumulants as well as amplitudes and lifetimes of
107 <    the dynamic form factor are compared with predictions of analytic
108 <    theories and found to be in very good agreement with them. For DNA
109 <    fragments 762 base pairs in length or longer, translational anisotropy
110 <    does not contribute significantly to dynamic light scattering. In
111 <    a comparison of rigid and flexible simulations on semistiff models
112 <    of this fragment, it is shown directly that flexing contributes
113 <    to the faster decay processes probed by light scattering and that
114 <    the flexible model studies are in good agreement with experiment.},
102 >        of rigid macromolecules. It is applied to polarized dynamic light
103 >        scattering from rodlike sturctures and from a model of a DNA fragment
104 >        (762 base pairs). A number of rod cases are examined in which the
105 >        translational anisotropy is increased form zero to a large value.
106 >        Simulated first cumulants as well as amplitudes and lifetimes of
107 >        the dynamic form factor are compared with predictions of analytic
108 >        theories and found to be in very good agreement with them. For DNA
109 >        fragments 762 base pairs in length or longer, translational anisotropy
110 >        does not contribute significantly to dynamic light scattering. In
111 >        a comparison of rigid and flexible simulations on semistiff models
112 >        of this fragment, it is shown directly that flexing contributes
113 >        to the faster decay processes probed by light scattering and that
114 >        the flexible model studies are in good agreement with experiment.},
115    annote = {Eu814 Times Cited:8 Cited References Count:32},
116    issn = {0024-9297},
117    uri = {<Go to ISI>://A1991EU81400029},
# Line 128 | Line 120 | Encoding: GBK
120   @ARTICLE{Andersen1983,
121    author = {H. C. Andersen},
122    title = {Rattle - a Velocity Version of the Shake Algorithm for Molecular-Dynamics
123 <    Calculations},
123 >        Calculations},
124    journal = {Journal of Computational Physics},
125    year = {1983},
126    volume = {52},
# Line 142 | Line 134 | Encoding: GBK
134   @ARTICLE{Auerbach2005,
135    author = {A. Auerbach},
136    title = {Gating of acetylcholine receptor channels: Brownian motion across
137 <    a broad transition state},
137 >        a broad transition state},
138    journal = {Proceedings of the National Academy of Sciences of the United States
139 <    of America},
139 >        of America},
140    year = {2005},
141    volume = {102},
142    pages = {1408-1412},
143    number = {5},
144    month = {Feb 1},
145    abstract = {Acetylcholine receptor channels (AChRs) are proteins that switch between
146 <    stable #closed# and #open# conformations. In patch clamp recordings,
147 <    diliganded AChR gating appears to be a simple, two-state reaction.
148 <    However, mutagenesis studies indicate that during gating dozens
149 <    of residues across the protein move asynchronously and are organized
150 <    into rigid body gating domains (#blocks#). Moreover, there is an
151 <    upper limit to the apparent channel opening rate constant. These
152 <    observations suggest that the gating reaction has a broad, corrugated
153 <    transition state region, with the maximum opening rate reflecting,
154 <    in part, the mean first-passage time across this ensemble. Simulations
155 <    reveal that a flat, isotropic energy profile for the transition
156 <    state can account for many of the essential features of AChR gating.
157 <    With this mechanism, concerted, local structural transitions that
158 <    occur on the broad transition state ensemble give rise to fractional
159 <    measures of reaction progress (Phi values) determined by rate-equilibrium
160 <    free energy relationship analysis. The results suggest that the
161 <    coarse-grained AChR gating conformational change propagates through
162 <    the protein with dynamics that are governed by the Brownian motion
163 <    of individual gating blocks.},
146 >        stable #closed# and #open# conformations. In patch clamp recordings,
147 >        diliganded AChR gating appears to be a simple, two-state reaction.
148 >        However, mutagenesis studies indicate that during gating dozens
149 >        of residues across the protein move asynchronously and are organized
150 >        into rigid body gating domains (#blocks#). Moreover, there is an
151 >        upper limit to the apparent channel opening rate constant. These
152 >        observations suggest that the gating reaction has a broad, corrugated
153 >        transition state region, with the maximum opening rate reflecting,
154 >        in part, the mean first-passage time across this ensemble. Simulations
155 >        reveal that a flat, isotropic energy profile for the transition
156 >        state can account for many of the essential features of AChR gating.
157 >        With this mechanism, concerted, local structural transitions that
158 >        occur on the broad transition state ensemble give rise to fractional
159 >        measures of reaction progress (Phi values) determined by rate-equilibrium
160 >        free energy relationship analysis. The results suggest that the
161 >        coarse-grained AChR gating conformational change propagates through
162 >        the protein with dynamics that are governed by the Brownian motion
163 >        of individual gating blocks.},
164    annote = {895QF Times Cited:9 Cited References Count:33},
165    issn = {0027-8424},
166    uri = {<Go to ISI>://000226877300030},
# Line 177 | Line 169 | Encoding: GBK
169   @ARTICLE{Baber1995,
170    author = {J. Baber and J. F. Ellena and D. S. Cafiso},
171    title = {Distribution of General-Anesthetics in Phospholipid-Bilayers Determined
172 <    Using H-2 Nmr and H-1-H-1 Noe Spectroscopy},
172 >        Using H-2 Nmr and H-1-H-1 Noe Spectroscopy},
173    journal = {Biochemistry},
174    year = {1995},
175    volume = {34},
# Line 185 | Line 177 | Encoding: GBK
177    number = {19},
178    month = {May 16},
179    abstract = {The effect of the general anesthetics halothane, enflurane, and isoflurane
180 <    on hydrocarbon chain packing in palmitoyl(d(31))oleoylphosphatidylcholine
181 <    membranes in the liquid crystalline phase was investigated using
182 <    H-2 NMR. Upon the addition of the anesthetics, the first five methylene
183 <    units near the interface generally show a very small increase in
184 <    segmental order, while segments deeper within the bilayer show a
185 <    small decrease in segmental order. From the H-2 NMR results, the
186 <    chain length for the perdeuterated palmitoyl chain in the absence
187 <    of anesthetic was found to be 12.35 Angstrom. Upon the addition
188 <    of halothane enflurane, or isoflurane, the acyl chain undergoes
189 <    slight contractions of 0.11, 0.20, or 0.16 Angstrom, respectively,
190 <    at 50 mol % anesthetic. A simple model was used to estimate the
191 <    relative amounts of anesthetic located near the interface and deeper
192 <    in the bilayer hydrocarbon region, and only a slight preference
193 <    for an interfacial location was observed. Intermolecular H-1-H-1
194 <    nuclear Overhauser effects (NOEs) were measured between phospholipid
195 <    and halothane protons. These NOEs are consistent with the intramembrane
196 <    location of the anesthetics suggested by the H-2 NMR data. In addition,
197 <    the NOE data indicate that anesthetics prefer the interfacial and
198 <    hydrocarbon regions of the membrane and are not found in high concentrations
199 <    in the phospholipid headgroup.},
180 >        on hydrocarbon chain packing in palmitoyl(d(31))oleoylphosphatidylcholine
181 >        membranes in the liquid crystalline phase was investigated using
182 >        H-2 NMR. Upon the addition of the anesthetics, the first five methylene
183 >        units near the interface generally show a very small increase in
184 >        segmental order, while segments deeper within the bilayer show a
185 >        small decrease in segmental order. From the H-2 NMR results, the
186 >        chain length for the perdeuterated palmitoyl chain in the absence
187 >        of anesthetic was found to be 12.35 Angstrom. Upon the addition
188 >        of halothane enflurane, or isoflurane, the acyl chain undergoes
189 >        slight contractions of 0.11, 0.20, or 0.16 Angstrom, respectively,
190 >        at 50 mol % anesthetic. A simple model was used to estimate the
191 >        relative amounts of anesthetic located near the interface and deeper
192 >        in the bilayer hydrocarbon region, and only a slight preference
193 >        for an interfacial location was observed. Intermolecular H-1-H-1
194 >        nuclear Overhauser effects (NOEs) were measured between phospholipid
195 >        and halothane protons. These NOEs are consistent with the intramembrane
196 >        location of the anesthetics suggested by the H-2 NMR data. In addition,
197 >        the NOE data indicate that anesthetics prefer the interfacial and
198 >        hydrocarbon regions of the membrane and are not found in high concentrations
199 >        in the phospholipid headgroup.},
200    annote = {Qz716 Times Cited:38 Cited References Count:37},
201    issn = {0006-2960},
202    uri = {<Go to ISI>://A1995QZ71600035},
# Line 213 | Line 205 | Encoding: GBK
205   @ARTICLE{Banerjee2004,
206    author = {D. Banerjee and B. C. Bag and S. K. Banik and D. S. Ray},
207    title = {Solution of quantum Langevin equation: Approximations, theoretical
208 <    and numerical aspects},
208 >        and numerical aspects},
209    journal = {Journal of Chemical Physics},
210    year = {2004},
211    volume = {120},
# Line 221 | Line 213 | Encoding: GBK
213    number = {19},
214    month = {May 15},
215    abstract = {Based on a coherent state representation of noise operator and an
216 <    ensemble averaging procedure using Wigner canonical thermal distribution
217 <    for harmonic oscillators, a generalized quantum Langevin equation
218 <    has been recently developed [Phys. Rev. E 65, 021109 (2002); 66,
219 <    051106 (2002)] to derive the equations of motion for probability
220 <    distribution functions in c-number phase-space. We extend the treatment
221 <    to explore several systematic approximation schemes for the solutions
222 <    of the Langevin equation for nonlinear potentials for a wide range
223 <    of noise correlation, strength and temperature down to the vacuum
224 <    limit. The method is exemplified by an analytic application to harmonic
225 <    oscillator for arbitrary memory kernel and with the help of a numerical
226 <    calculation of barrier crossing, in a cubic potential to demonstrate
227 <    the quantum Kramers' turnover and the quantum Arrhenius plot. (C)
228 <    2004 American Institute of Physics.},
216 >        ensemble averaging procedure using Wigner canonical thermal distribution
217 >        for harmonic oscillators, a generalized quantum Langevin equation
218 >        has been recently developed [Phys. Rev. E 65, 021109 (2002); 66,
219 >        051106 (2002)] to derive the equations of motion for probability
220 >        distribution functions in c-number phase-space. We extend the treatment
221 >        to explore several systematic approximation schemes for the solutions
222 >        of the Langevin equation for nonlinear potentials for a wide range
223 >        of noise correlation, strength and temperature down to the vacuum
224 >        limit. The method is exemplified by an analytic application to harmonic
225 >        oscillator for arbitrary memory kernel and with the help of a numerical
226 >        calculation of barrier crossing, in a cubic potential to demonstrate
227 >        the quantum Kramers' turnover and the quantum Arrhenius plot. (C)
228 >        2004 American Institute of Physics.},
229    annote = {816YY Times Cited:8 Cited References Count:35},
230    issn = {0021-9606},
231    uri = {<Go to ISI>://000221146400009},
232   }
233  
242 @ARTICLE{Barojas1973,
243  author = {J. Barojas and D. Levesque},
244  title = {Simulation of Diatomic Homonuclear Liquids},
245  journal = {Phys. Rev. A},
246  year = {1973},
247  volume = {7},
248  pages = {1092-1105},
249 }
250
234   @ARTICLE{Barth1998,
235    author = {E. Barth and T. Schlick},
236    title = {Overcoming stability limitations in biomolecular dynamics. I. Combining
237 <    force splitting via extrapolation with Langevin dynamics in LN},
237 >        force splitting via extrapolation with Langevin dynamics in LN},
238    journal = {Journal of Chemical Physics},
239    year = {1998},
240    volume = {109},
# Line 259 | Line 242 | Encoding: GBK
242    number = {5},
243    month = {Aug 1},
244    abstract = {We present an efficient new method termed LN for propagating biomolecular
245 <    dynamics according to the Langevin equation that arose fortuitously
246 <    upon analysis of the range of harmonic validity of our normal-mode
247 <    scheme LIN. LN combines force linearization with force splitting
248 <    techniques and disposes of LIN'S computationally intensive minimization
249 <    (anharmonic correction) component. Unlike the competitive multiple-timestepping
250 <    (MTS) schemes today-formulated to be symplectic and time-reversible-LN
251 <    merges the slow and fast forces via extrapolation rather than impulses;
252 <    the Langevin heat bath prevents systematic energy drifts. This combination
253 <    succeeds in achieving more significant speedups than these MTS methods
254 <    which are Limited by resonance artifacts to an outer timestep less
255 <    than some integer multiple of half the period of the fastest motion
256 <    (around 4-5 fs for biomolecules). We show that LN achieves very
257 <    good agreement with small-timestep solutions of the Langevin equation
258 <    in terms of thermodynamics (energy means and variances), geometry,
259 <    and dynamics (spectral densities) for two proteins in vacuum and
260 <    a large water system. Significantly, the frequency of updating the
261 <    slow forces extends to 48 fs or more, resulting in speedup factors
262 <    exceeding 10. The implementation of LN in any program that employs
263 <    force-splitting computations is straightforward, with only partial
264 <    second-derivative information required, as well as sparse Hessian/vector
265 <    multiplication routines. The linearization part of LN could even
266 <    be replaced by direct evaluation of the fast components. The application
267 <    of LN to biomolecular dynamics is well suited for configurational
268 <    sampling, thermodynamic, and structural questions. (C) 1998 American
269 <    Institute of Physics.},
245 >        dynamics according to the Langevin equation that arose fortuitously
246 >        upon analysis of the range of harmonic validity of our normal-mode
247 >        scheme LIN. LN combines force linearization with force splitting
248 >        techniques and disposes of LIN'S computationally intensive minimization
249 >        (anharmonic correction) component. Unlike the competitive multiple-timestepping
250 >        (MTS) schemes today-formulated to be symplectic and time-reversible-LN
251 >        merges the slow and fast forces via extrapolation rather than impulses;
252 >        the Langevin heat bath prevents systematic energy drifts. This combination
253 >        succeeds in achieving more significant speedups than these MTS methods
254 >        which are Limited by resonance artifacts to an outer timestep less
255 >        than some integer multiple of half the period of the fastest motion
256 >        (around 4-5 fs for biomolecules). We show that LN achieves very
257 >        good agreement with small-timestep solutions of the Langevin equation
258 >        in terms of thermodynamics (energy means and variances), geometry,
259 >        and dynamics (spectral densities) for two proteins in vacuum and
260 >        a large water system. Significantly, the frequency of updating the
261 >        slow forces extends to 48 fs or more, resulting in speedup factors
262 >        exceeding 10. The implementation of LN in any program that employs
263 >        force-splitting computations is straightforward, with only partial
264 >        second-derivative information required, as well as sparse Hessian/vector
265 >        multiplication routines. The linearization part of LN could even
266 >        be replaced by direct evaluation of the fast components. The application
267 >        of LN to biomolecular dynamics is well suited for configurational
268 >        sampling, thermodynamic, and structural questions. (C) 1998 American
269 >        Institute of Physics.},
270    annote = {105HH Times Cited:29 Cited References Count:49},
271    issn = {0021-9606},
272    uri = {<Go to ISI>://000075066300006},
# Line 292 | Line 275 | Encoding: GBK
275   @ARTICLE{Batcho2001,
276    author = {P. F. Batcho and T. Schlick},
277    title = {Special stability advantages of position-Verlet over velocity-Verlet
278 <    in multiple-time step integration},
278 >        in multiple-time step integration},
279    journal = {Journal of Chemical Physics},
280    year = {2001},
281    volume = {115},
# Line 300 | Line 283 | Encoding: GBK
283    number = {9},
284    month = {Sep 1},
285    abstract = {We present an analysis for a simple two-component harmonic oscillator
286 <    that compares the use of position-Verlet to velocity-Verlet for
287 <    multiple-time step integration. The numerical stability analysis
288 <    based on the impulse-Verlet splitting shows that position-Verlet
289 <    has enhanced stability, in terms of the largest allowable time step,
290 <    for cases where an ample separation of time scales exists. Numerical
291 <    investigations confirm the advantages of the position-Verlet scheme
292 <    when used for the fastest time scales of the system. Applications
293 <    to a biomolecule. a solvated protein, for both Newtonian and Langevin
294 <    dynamics echo these trends over large outer time-step regimes. (C)
295 <    2001 American Institute of Physics.},
286 >        that compares the use of position-Verlet to velocity-Verlet for
287 >        multiple-time step integration. The numerical stability analysis
288 >        based on the impulse-Verlet splitting shows that position-Verlet
289 >        has enhanced stability, in terms of the largest allowable time step,
290 >        for cases where an ample separation of time scales exists. Numerical
291 >        investigations confirm the advantages of the position-Verlet scheme
292 >        when used for the fastest time scales of the system. Applications
293 >        to a biomolecule. a solvated protein, for both Newtonian and Langevin
294 >        dynamics echo these trends over large outer time-step regimes. (C)
295 >        2001 American Institute of Physics.},
296    annote = {469KV Times Cited:6 Cited References Count:30},
297    issn = {0021-9606},
298    uri = {<Go to ISI>://000170813800005},
# Line 318 | Line 301 | Encoding: GBK
301   @ARTICLE{Bates2005,
302    author = {M. A. Bates and G. R. Luckhurst},
303    title = {Biaxial nematic phases and V-shaped molecules: A Monte Carlo simulation
304 <    study},
304 >        study},
305    journal = {Physical Review E},
306    year = {2005},
307    volume = {72},
# Line 326 | Line 309 | Encoding: GBK
309    number = {5},
310    month = {Nov},
311    abstract = {Inspired by recent claims that compounds composed of V-shaped molecules
312 <    can exhibit the elusive biaxial nematic phase, we have developed
313 <    a generic simulation model for such systems. This contains the features
314 <    of the molecule that are essential to its liquid crystal behavior,
315 <    namely the anisotropies of the two arms and the angle between them.
316 <    The behavior of the model has been investigated using Monte Carlo
317 <    simulations for a wide range of these structural parameters. This
318 <    allows us to establish the relationship between the V-shaped molecule
319 <    and its ability to form a biaxial nematic phase. Of particular importance
320 <    are the criteria of geometry and the relative anisotropy necessary
321 <    for the system to exhibit a Landau point, at which the biaxial nematic
322 <    is formed directly from the isotropic phase. The simulations have
323 <    also been used to determine the orientational order parameters for
324 <    a selection of molecular axes. These are especially important because
325 <    they reveal the phase symmetry and are connected to the experimental
326 <    determination of this. The simulation results show that, whereas
327 <    some positions are extremely sensitive to the phase biaxiality,
328 <    others are totally blind to this.},
312 >        can exhibit the elusive biaxial nematic phase, we have developed
313 >        a generic simulation model for such systems. This contains the features
314 >        of the molecule that are essential to its liquid crystal behavior,
315 >        namely the anisotropies of the two arms and the angle between them.
316 >        The behavior of the model has been investigated using Monte Carlo
317 >        simulations for a wide range of these structural parameters. This
318 >        allows us to establish the relationship between the V-shaped molecule
319 >        and its ability to form a biaxial nematic phase. Of particular importance
320 >        are the criteria of geometry and the relative anisotropy necessary
321 >        for the system to exhibit a Landau point, at which the biaxial nematic
322 >        is formed directly from the isotropic phase. The simulations have
323 >        also been used to determine the orientational order parameters for
324 >        a selection of molecular axes. These are especially important because
325 >        they reveal the phase symmetry and are connected to the experimental
326 >        determination of this. The simulation results show that, whereas
327 >        some positions are extremely sensitive to the phase biaxiality,
328 >        others are totally blind to this.},
329    annote = {Part 1 988LQ Times Cited:0 Cited References Count:38},
330    issn = {1539-3755},
331    uri = {<Go to ISI>://000233603100030},
# Line 358 | Line 341 | Encoding: GBK
341    number = {5},
342    month = {Nov 1},
343    abstract = {We introduce an unbiased protocol for performing rotational moves
344 <    in rigid-body dynamics simulations. This approach - based on the
345 <    analytic solution for the rotational equations of motion for an
346 <    orthogonal coordinate system at constant angular velocity - removes
347 <    deficiencies that have been largely ignored in Brownian dynamics
348 <    simulations, namely errors for finite rotations that result from
349 <    applying the noncommuting rotational matrices in an arbitrary order.
350 <    Our algorithm should thus replace standard approaches to rotate
351 <    local coordinate frames in Langevin and Brownian dynamics simulations.},
344 >        in rigid-body dynamics simulations. This approach - based on the
345 >        analytic solution for the rotational equations of motion for an
346 >        orthogonal coordinate system at constant angular velocity - removes
347 >        deficiencies that have been largely ignored in Brownian dynamics
348 >        simulations, namely errors for finite rotations that result from
349 >        applying the noncommuting rotational matrices in an arbitrary order.
350 >        Our algorithm should thus replace standard approaches to rotate
351 >        local coordinate frames in Langevin and Brownian dynamics simulations.},
352    annote = {736UA Times Cited:0 Cited References Count:11},
353    issn = {0006-3495},
354    uri = {<Go to ISI>://000186190500018},
# Line 374 | Line 357 | Encoding: GBK
357   @ARTICLE{Beloborodov1998,
358    author = {I. S. Beloborodov and V. Y. Orekhov and A. S. Arseniev},
359    title = {Effect of coupling between rotational and translational Brownian
360 <    motions on NMR spin relaxation: Consideration using green function
361 <    of rigid body diffusion},
360 >        motions on NMR spin relaxation: Consideration using green function
361 >        of rigid body diffusion},
362    journal = {Journal of Magnetic Resonance},
363    year = {1998},
364    volume = {132},
# Line 383 | Line 366 | Encoding: GBK
366    number = {2},
367    month = {Jun},
368    abstract = {Using the Green function of arbitrary rigid Brownian diffusion (Goldstein,
369 <    Biopolymers 33, 409-436, 1993), it was analytically shown that coupling
370 <    between translation and rotation diffusion degrees of freedom does
371 <    not affect the correlation functions relevant to the NMR intramolecular
372 <    relaxation. It follows that spectral densities usually used for
373 <    the anisotropic rotation diffusion (Woessner, J. Chem. Phys. 37,
374 <    647-654, 1962) can be regarded as exact in respect to the rotation-translation
375 <    coupling for the spin system connected with a rigid body. (C) 1998
376 <    Academic Press.},
369 >        Biopolymers 33, 409-436, 1993), it was analytically shown that coupling
370 >        between translation and rotation diffusion degrees of freedom does
371 >        not affect the correlation functions relevant to the NMR intramolecular
372 >        relaxation. It follows that spectral densities usually used for
373 >        the anisotropic rotation diffusion (Woessner, J. Chem. Phys. 37,
374 >        647-654, 1962) can be regarded as exact in respect to the rotation-translation
375 >        coupling for the spin system connected with a rigid body. (C) 1998
376 >        Academic Press.},
377    annote = {Zu605 Times Cited:2 Cited References Count:6},
378    issn = {1090-7807},
379    uri = {<Go to ISI>://000074214800017},
# Line 399 | Line 382 | Encoding: GBK
382   @ARTICLE{Berardi1996,
383    author = {R. Berardi and S. Orlandi and C. Zannoni},
384    title = {Antiphase structures in polar smectic liquid crystals and their molecular
385 <    origin},
385 >        origin},
386    journal = {Chemical Physics Letters},
387    year = {1996},
388    volume = {261},
# Line 407 | Line 390 | Encoding: GBK
390    number = {3},
391    month = {Oct 18},
392    abstract = {We demonstrate that the overall molecular dipole organization in a
393 <    smectic liquid crystal formed of polar molecules can be strongly
394 <    influenced by the position of the dipole in the molecule. We study
395 <    by large scale Monte Carlo simulations systems of attractive-repulsive
396 <    ''Gay-Berne'' elongated ellipsoids with an axial dipole at the center
397 <    or near the end of the molecule and we show that monolayer smectic
398 <    liquid crystals and modulated antiferroelectric bilayer stripe domains
399 <    similar to the experimentally observed ''antiphase'' structures
400 <    are obtained in the two cases.},
393 >        smectic liquid crystal formed of polar molecules can be strongly
394 >        influenced by the position of the dipole in the molecule. We study
395 >        by large scale Monte Carlo simulations systems of attractive-repulsive
396 >        ''Gay-Berne'' elongated ellipsoids with an axial dipole at the center
397 >        or near the end of the molecule and we show that monolayer smectic
398 >        liquid crystals and modulated antiferroelectric bilayer stripe domains
399 >        similar to the experimentally observed ''antiphase'' structures
400 >        are obtained in the two cases.},
401    annote = {Vn637 Times Cited:49 Cited References Count:26},
402    issn = {0009-2614},
403    uri = {<Go to ISI>://A1996VN63700023},
# Line 422 | Line 405 | Encoding: GBK
405  
406   @ARTICLE{Berkov2005,
407    author = {D. V. Berkov and N. L. Gorn},
408 +  title = {Stochastic dynamic simulations of fast remagnetization processes:
409 +        recent advances and applications},
410 +  journal = {Journal of Magnetism and Magnetic Materials},
411 +  year = {2005},
412 +  volume = {290},
413 +  pages = {442-448},
414 +  month = {Apr},
415 +  abstract = {Numerical simulations of fast remagnetization processes using stochastic
416 +        dynamics are widely used to study various magnetic systems. In this
417 +        paper, we first address several crucial methodological problems
418 +        of such simulations: (i) the influence of finite-element discretization
419 +        on simulated dynamics, (ii) choice between Ito and Stratonovich
420 +        stochastic calculi by the solution of micromagnetic stochastic equations
421 +        of motion and (iii) non-trivial correlation properties of the random
422 +        (thermal) field. Next, we discuss several examples to demonstrate
423 +        the great potential of the Langevin dynamics for studying fast remagnetization
424 +        processes in technically relevant applications: we present numerical
425 +        analysis of equilibrium magnon spectra in patterned structures,
426 +        study thermal noise effects on the magnetization dynamics of nanoelements
427 +        in pulsed fields and show some results for a remagnetization dynamics
428 +        induced by a spin-polarized current. (c) 2004 Elsevier B.V. All
429 +        rights reserved.},
430 +  annote = {Part 1 Sp. Iss. SI 922KU Times Cited:2 Cited References Count:25},
431 +  issn = {0304-8853},
432 +  uri = {<Go to ISI>://000228837600109},
433 + }
434 +
435 + @ARTICLE{Berkov2005a,
436 +  author = {D. V. Berkov and N. L. Gorn},
437    title = {Magnetization precession due to a spin-polarized current in a thin
438 <    nanoelement: Numerical simulation study},
438 >        nanoelement: Numerical simulation study},
439    journal = {Physical Review B},
440    year = {2005},
441    volume = {72},
# Line 431 | Line 443 | Encoding: GBK
443    number = {9},
444    month = {Sep},
445    abstract = {In this paper a detailed numerical study (in frames of the Slonczewski
446 <    formalism) of magnetization oscillations driven by a spin-polarized
447 <    current through a thin elliptical nanoelement is presented. We show
448 <    that a sophisticated micromagnetic model, where a polycrystalline
449 <    structure of a nanoelement is taken into account, can explain qualitatively
450 <    all most important features of the magnetization oscillation spectra
451 <    recently observed experimentally [S. I. Kiselev , Nature 425, 380
452 <    (2003)], namely, existence of several equidistant spectral bands,
453 <    sharp onset and abrupt disappearance of magnetization oscillations
454 <    with increasing current, absence of the out-of-plane regime predicted
455 <    by a macrospin model, and the relation between frequencies of so-called
456 <    small-angle and quasichaotic oscillations. However, a quantitative
457 <    agreement with experimental results (especially concerning the frequency
458 <    of quasichaotic oscillations) could not be achieved in the region
459 <    of reasonable parameter values, indicating that further model refinement
460 <    is necessary for a complete understanding of the spin-driven magnetization
461 <    precession even in this relatively simple experimental situation.},
446 >        formalism) of magnetization oscillations driven by a spin-polarized
447 >        current through a thin elliptical nanoelement is presented. We show
448 >        that a sophisticated micromagnetic model, where a polycrystalline
449 >        structure of a nanoelement is taken into account, can explain qualitatively
450 >        all most important features of the magnetization oscillation spectra
451 >        recently observed experimentally [S. I. Kiselev , Nature 425, 380
452 >        (2003)], namely, existence of several equidistant spectral bands,
453 >        sharp onset and abrupt disappearance of magnetization oscillations
454 >        with increasing current, absence of the out-of-plane regime predicted
455 >        by a macrospin model, and the relation between frequencies of so-called
456 >        small-angle and quasichaotic oscillations. However, a quantitative
457 >        agreement with experimental results (especially concerning the frequency
458 >        of quasichaotic oscillations) could not be achieved in the region
459 >        of reasonable parameter values, indicating that further model refinement
460 >        is necessary for a complete understanding of the spin-driven magnetization
461 >        precession even in this relatively simple experimental situation.},
462    annote = {969IT Times Cited:2 Cited References Count:55},
463    issn = {1098-0121},
464    uri = {<Go to ISI>://000232228500058},
465   }
466  
455 @ARTICLE{Berkov2005a,
456  author = {D. V. Berkov and N. L. Gorn},
457  title = {Stochastic dynamic simulations of fast remagnetization processes:
458    recent advances and applications},
459  journal = {Journal of Magnetism and Magnetic Materials},
460  year = {2005},
461  volume = {290},
462  pages = {442-448},
463  month = {Apr},
464  abstract = {Numerical simulations of fast remagnetization processes using stochastic
465    dynamics are widely used to study various magnetic systems. In this
466    paper, we first address several crucial methodological problems
467    of such simulations: (i) the influence of finite-element discretization
468    on simulated dynamics, (ii) choice between Ito and Stratonovich
469    stochastic calculi by the solution of micromagnetic stochastic equations
470    of motion and (iii) non-trivial correlation properties of the random
471    (thermal) field. Next, we discuss several examples to demonstrate
472    the great potential of the Langevin dynamics for studying fast remagnetization
473    processes in technically relevant applications: we present numerical
474    analysis of equilibrium magnon spectra in patterned structures,
475    study thermal noise effects on the magnetization dynamics of nanoelements
476    in pulsed fields and show some results for a remagnetization dynamics
477    induced by a spin-polarized current. (c) 2004 Elsevier B.V. All
478    rights reserved.},
479  annote = {Part 1 Sp. Iss. SI 922KU Times Cited:2 Cited References Count:25},
480  issn = {0304-8853},
481  uri = {<Go to ISI>://000228837600109},
482 }
483
467   @ARTICLE{Berkov2002,
468    author = {D. V. Berkov and N. L. Gorn and P. Gornert},
469    title = {Magnetization dynamics in nanoparticle systems: Numerical simulation
470 <    using Langevin dynamics},
470 >        using Langevin dynamics},
471    journal = {Physica Status Solidi a-Applied Research},
472    year = {2002},
473    volume = {189},
# Line 492 | Line 475 | Encoding: GBK
475    number = {2},
476    month = {Feb 16},
477    abstract = {We report on recent progress achieved by the development of numerical
478 <    methods based on the stochastic (Langevin) dynamics applied to systems
479 <    of interacting magnetic nanoparticles. The method enables direct
480 <    simulations of the trajectories of magnetic moments taking into
481 <    account (i) all relevant interactions, (ii) precession dynamics,
482 <    and (iii) temperature fluctuations included via the random (thermal)
483 <    field. We present several novel results obtained using new methods
484 <    developed for the solution of the Langevin equations. In particular,
485 <    we have investigated magnetic nanodots and disordered granular systems
486 <    of single-domain magnetic particles. For the first case we have
487 <    calculated the spectrum and the spatial distribution of spin excitations.
488 <    For the second system the complex ac susceptibility chi(omega, T)
489 <    for various particle concentrations and particle anisotropies were
490 <    computed and compared with numerous experimental results.},
478 >        methods based on the stochastic (Langevin) dynamics applied to systems
479 >        of interacting magnetic nanoparticles. The method enables direct
480 >        simulations of the trajectories of magnetic moments taking into
481 >        account (i) all relevant interactions, (ii) precession dynamics,
482 >        and (iii) temperature fluctuations included via the random (thermal)
483 >        field. We present several novel results obtained using new methods
484 >        developed for the solution of the Langevin equations. In particular,
485 >        we have investigated magnetic nanodots and disordered granular systems
486 >        of single-domain magnetic particles. For the first case we have
487 >        calculated the spectrum and the spatial distribution of spin excitations.
488 >        For the second system the complex ac susceptibility chi(omega, T)
489 >        for various particle concentrations and particle anisotropies were
490 >        computed and compared with numerous experimental results.},
491    annote = {526TF Times Cited:4 Cited References Count:37},
492    issn = {0031-8965},
493    uri = {<Go to ISI>://000174145200026},
# Line 513 | Line 496 | Encoding: GBK
496   @ARTICLE{Bernal1980,
497    author = {J.M. Bernal and J. G. {de la Torre}},
498    title = {Transport Properties and Hydrodynamic Centers of Rigid Macromolecules
499 <    with Arbitrary Shape},
499 >        with Arbitrary Shape},
500    journal = {Biopolymers},
501    year = {1980},
502    volume = {19},
# Line 523 | Line 506 | Encoding: GBK
506   @ARTICLE{Brunger1984,
507    author = {A. Brunger and C. L. Brooks and M. Karplus},
508    title = {Stochastic Boundary-Conditions for Molecular-Dynamics Simulations
509 <    of St2 Water},
509 >        of St2 Water},
510    journal = {Chemical Physics Letters},
511    year = {1984},
512    volume = {105},
# Line 537 | Line 520 | Encoding: GBK
520   @ARTICLE{Budd1999,
521    author = {C. J. Budd and G. J. Collins and W. Z. Huang and R. D. Russell},
522    title = {Self-similar numerical solutions of the porous-medium equation using
523 <    moving mesh methods},
523 >        moving mesh methods},
524    journal = {Philosophical Transactions of the Royal Society of London Series
525 <    a-Mathematical Physical and Engineering Sciences},
525 >        a-Mathematical Physical and Engineering Sciences},
526    year = {1999},
527    volume = {357},
528    pages = {1047-1077},
529    number = {1754},
530    month = {Apr 15},
531    abstract = {This paper examines a synthesis of adaptive mesh methods with the
532 <    use of symmetry to study a partial differential equation. In particular,
533 <    it considers methods which admit discrete self-similar solutions,
534 <    examining the convergence of these to the true self-similar solution
535 <    as well as their stability. Special attention is given to the nonlinear
536 <    diffusion equation describing flow in a porous medium.},
532 >        use of symmetry to study a partial differential equation. In particular,
533 >        it considers methods which admit discrete self-similar solutions,
534 >        examining the convergence of these to the true self-similar solution
535 >        as well as their stability. Special attention is given to the nonlinear
536 >        diffusion equation describing flow in a porous medium.},
537    annote = {199EE Times Cited:4 Cited References Count:14},
538    issn = {1364-503X},
539    uri = {<Go to ISI>://000080466800005},
# Line 566 | Line 549 | Encoding: GBK
549    number = {21},
550    month = {Dec 1},
551    abstract = {Fluids of hard bent-core molecules have been studied using theory
552 <    and computer simulation. The molecules are composed of two hard
553 <    spherocylinders, with length-to-breadth ratio L/D, joined by their
554 <    ends at an angle 180 degrees - gamma. For L/D = 2 and gamma = 0,10,20
555 <    degrees, the simulations show isotropic, nematic, smectic, and solid
556 <    phases. For L/D = 2 and gamma = 30 degrees, only isotropic, nematic,
557 <    and solid phases are in evidence, which suggests that there is a
558 <    nematic-smectic-solid triple point at an angle in the range 20 degrees
559 <    < gamma < 30 degrees. In all of the orientationally ordered fluid
560 <    phases the order is purely uniaxial. For gamma = 10 degrees and
561 <    20 degrees, at the studied densities, the solid is also uniaxially
562 <    ordered, whilst for gamma = 30 degrees the solid layers are biaxially
563 <    ordered. For L/D = 2 and gamma = 60 degrees and 90 degrees we find
564 <    no spontaneous orientational ordering. This is shown to be due to
565 <    the interlocking of dimer pairs which precludes alignment. We find
566 <    similar results for L/D = 9.5 and gamma = 72 degrees, where an isotropic-biaxial
567 <    nematic transition is predicted by Onsager theory. Simulations in
568 <    the biaxial nematic phase show it to be at least mechanically stable
569 <    with respect to the isotropic phase, however. We have compared the
570 <    quasi-exact simulation results in the isotropic phase with the predicted
571 <    equations of state from three theories: the virial expansion containing
572 <    the second and third virial coefficients; the Parsons-Lee equation
573 <    of state; an application of Wertheim's theory of associating fluids
574 <    in the limit of infinite attractive association energy. For all
575 <    of the molecule elongations and geometries we have simulated, the
576 <    Wertheim theory proved to be the most accurate. Interestingly, the
577 <    isotropic equation of state is virtually independent of the dimer
578 <    bond angle-a feature that is also reflected in the lack of variation
579 <    with angle of the calculated second and third virial coefficients.
580 <    (C) 1999 American Institute of Physics. [S0021-9606(99)50445-5].},
552 >        and computer simulation. The molecules are composed of two hard
553 >        spherocylinders, with length-to-breadth ratio L/D, joined by their
554 >        ends at an angle 180 degrees - gamma. For L/D = 2 and gamma = 0,10,20
555 >        degrees, the simulations show isotropic, nematic, smectic, and solid
556 >        phases. For L/D = 2 and gamma = 30 degrees, only isotropic, nematic,
557 >        and solid phases are in evidence, which suggests that there is a
558 >        nematic-smectic-solid triple point at an angle in the range 20 degrees
559 >        < gamma < 30 degrees. In all of the orientationally ordered fluid
560 >        phases the order is purely uniaxial. For gamma = 10 degrees and
561 >        20 degrees, at the studied densities, the solid is also uniaxially
562 >        ordered, whilst for gamma = 30 degrees the solid layers are biaxially
563 >        ordered. For L/D = 2 and gamma = 60 degrees and 90 degrees we find
564 >        no spontaneous orientational ordering. This is shown to be due to
565 >        the interlocking of dimer pairs which precludes alignment. We find
566 >        similar results for L/D = 9.5 and gamma = 72 degrees, where an isotropic-biaxial
567 >        nematic transition is predicted by Onsager theory. Simulations in
568 >        the biaxial nematic phase show it to be at least mechanically stable
569 >        with respect to the isotropic phase, however. We have compared the
570 >        quasi-exact simulation results in the isotropic phase with the predicted
571 >        equations of state from three theories: the virial expansion containing
572 >        the second and third virial coefficients; the Parsons-Lee equation
573 >        of state; an application of Wertheim's theory of associating fluids
574 >        in the limit of infinite attractive association energy. For all
575 >        of the molecule elongations and geometries we have simulated, the
576 >        Wertheim theory proved to be the most accurate. Interestingly, the
577 >        isotropic equation of state is virtually independent of the dimer
578 >        bond angle-a feature that is also reflected in the lack of variation
579 >        with angle of the calculated second and third virial coefficients.
580 >        (C) 1999 American Institute of Physics. [S0021-9606(99)50445-5].},
581    annote = {255TC Times Cited:24 Cited References Count:38},
582    issn = {0021-9606},
583    uri = {<Go to ISI>://000083685400056},
# Line 610 | Line 593 | Encoding: GBK
593    number = {11},
594    month = {Nov},
595    abstract = {A review is presented of molecular and mesoscopic computer simulations
596 <    of liquid crystalline systems. Molecular simulation approaches applied
597 <    to such systems are described, and the key findings for bulk phase
598 <    behaviour are reported. Following this, recently developed lattice
599 <    Boltzmann approaches to the mesoscale modelling of nemato-dynanics
600 <    are reviewed. This paper concludes with a discussion of possible
601 <    areas for future development in this field.},
596 >        of liquid crystalline systems. Molecular simulation approaches applied
597 >        to such systems are described, and the key findings for bulk phase
598 >        behaviour are reported. Following this, recently developed lattice
599 >        Boltzmann approaches to the mesoscale modelling of nemato-dynanics
600 >        are reviewed. This paper concludes with a discussion of possible
601 >        areas for future development in this field.},
602    annote = {989TU Times Cited:2 Cited References Count:258},
603    issn = {0034-4885},
604    uri = {<Go to ISI>://000233697600004},
# Line 624 | Line 607 | Encoding: GBK
607   @ARTICLE{Carrasco1999,
608    author = {B. Carrasco and J. G. {de la Torre}},
609    title = {Hydrodynamic properties of rigid particles: Comparison of different
610 <    modeling and computational procedures},
610 >        modeling and computational procedures},
611    journal = {Biophysical Journal},
612    year = {1999},
613    volume = {76},
# Line 632 | Line 615 | Encoding: GBK
615    number = {6},
616    month = {Jun},
617    abstract = {The hydrodynamic properties of rigid particles are calculated from
618 <    models composed of spherical elements (beads) using theories developed
619 <    by Kirkwood, Bloomfield, and their coworkers. Bead models have usually
620 <    been built in such a way that the beads fill the volume occupied
621 <    by the particles. Sometimes the beads are few and of varying sizes
622 <    (bead models in the strict sense), and other times there are many
623 <    small beads (filling models). Because hydrodynamic friction takes
624 <    place at the molecular surface, another possibility is to use shell
625 <    models, as originally proposed by Bloomfield. In this work, we have
626 <    developed procedures to build models of the various kinds, and we
627 <    describe the theory and methods for calculating their hydrodynamic
628 <    properties, including approximate methods that may be needed to
629 <    treat models with a very large number of elements. By combining
630 <    the various possibilities of model building and hydrodynamic calculation,
631 <    several strategies can be designed. We have made a quantitative
632 <    comparison of the performance of the various strategies by applying
633 <    them to some test cases, for which the properties are known a priori.
634 <    We provide guidelines and computational tools for bead modeling.},
618 >        models composed of spherical elements (beads) using theories developed
619 >        by Kirkwood, Bloomfield, and their coworkers. Bead models have usually
620 >        been built in such a way that the beads fill the volume occupied
621 >        by the particles. Sometimes the beads are few and of varying sizes
622 >        (bead models in the strict sense), and other times there are many
623 >        small beads (filling models). Because hydrodynamic friction takes
624 >        place at the molecular surface, another possibility is to use shell
625 >        models, as originally proposed by Bloomfield. In this work, we have
626 >        developed procedures to build models of the various kinds, and we
627 >        describe the theory and methods for calculating their hydrodynamic
628 >        properties, including approximate methods that may be needed to
629 >        treat models with a very large number of elements. By combining
630 >        the various possibilities of model building and hydrodynamic calculation,
631 >        several strategies can be designed. We have made a quantitative
632 >        comparison of the performance of the various strategies by applying
633 >        them to some test cases, for which the properties are known a priori.
634 >        We provide guidelines and computational tools for bead modeling.},
635    annote = {200TT Times Cited:46 Cited References Count:57},
636    issn = {0006-3495},
637    uri = {<Go to ISI>://000080556700016},
# Line 657 | Line 640 | Encoding: GBK
640   @ARTICLE{Chandra1999,
641    author = {A. Chandra and T. Ichiye},
642    title = {Dynamical properties of the soft sticky dipole model of water: Molecular
643 <    dynamics simulations},
643 >        dynamics simulations},
644    journal = {Journal of Chemical Physics},
645    year = {1999},
646    volume = {111},
# Line 665 | Line 648 | Encoding: GBK
648    number = {6},
649    month = {Aug 8},
650    abstract = {Dynamical properties of the soft sticky dipole (SSD) model of water
651 <    are calculated by means of molecular dynamics simulations. Since
652 <    this is not a simple point model, the forces and torques arising
653 <    from the SSD potential are derived here. Simulations are carried
654 <    out in the microcanonical ensemble employing the Ewald method for
655 <    the electrostatic interactions. Various time correlation functions
656 <    and dynamical quantities associated with the translational and rotational
657 <    motion of water molecules are evaluated and compared with those
658 <    of two other commonly used models of liquid water, namely the transferable
659 <    intermolecular potential-three points (TIP3P) and simple point charge/extended
660 <    (SPC/E) models, and also with experiments. The dynamical properties
661 <    of the SSD water model are found to be in good agreement with the
662 <    experimental results and appear to be better than the TIP3P and
663 <    SPC/E models in most cases, as has been previously shown for its
664 <    thermodynamic, structural, and dielectric properties. Also, molecular
665 <    dynamics simulations of the SSD model are found to run much faster
666 <    than TIP3P, SPC/E, and other multisite models. (C) 1999 American
667 <    Institute of Physics. [S0021-9606(99)51430-X].},
651 >        are calculated by means of molecular dynamics simulations. Since
652 >        this is not a simple point model, the forces and torques arising
653 >        from the SSD potential are derived here. Simulations are carried
654 >        out in the microcanonical ensemble employing the Ewald method for
655 >        the electrostatic interactions. Various time correlation functions
656 >        and dynamical quantities associated with the translational and rotational
657 >        motion of water molecules are evaluated and compared with those
658 >        of two other commonly used models of liquid water, namely the transferable
659 >        intermolecular potential-three points (TIP3P) and simple point charge/extended
660 >        (SPC/E) models, and also with experiments. The dynamical properties
661 >        of the SSD water model are found to be in good agreement with the
662 >        experimental results and appear to be better than the TIP3P and
663 >        SPC/E models in most cases, as has been previously shown for its
664 >        thermodynamic, structural, and dielectric properties. Also, molecular
665 >        dynamics simulations of the SSD model are found to run much faster
666 >        than TIP3P, SPC/E, and other multisite models. (C) 1999 American
667 >        Institute of Physics. [S0021-9606(99)51430-X].},
668    annote = {221EN Times Cited:14 Cited References Count:66},
669    issn = {0021-9606},
670    uri = {<Go to ISI>://000081711200038},
# Line 711 | Line 694 | Encoding: GBK
694    number = {1-2},
695    month = {Jan},
696    abstract = {We investigate the asymptotic behavior of systems of nonlinear differential
697 <    equations and introduce a family of mixed methods from combinations
698 <    of explicit Runge-Kutta methods. These methods have better stability
699 <    behavior than traditional Runge-Kutta methods and generally extend
700 <    the range of validity of the calculated solutions. These methods
701 <    also give a way of determining if the numerical solutions are real
702 <    or spurious. Emphasis is put on examples coming from mathematical
703 <    models in ecology. (C) 2002 IMACS. Published by Elsevier Science
704 <    B.V. All rights reserved.},
697 >        equations and introduce a family of mixed methods from combinations
698 >        of explicit Runge-Kutta methods. These methods have better stability
699 >        behavior than traditional Runge-Kutta methods and generally extend
700 >        the range of validity of the calculated solutions. These methods
701 >        also give a way of determining if the numerical solutions are real
702 >        or spurious. Emphasis is put on examples coming from mathematical
703 >        models in ecology. (C) 2002 IMACS. Published by Elsevier Science
704 >        B.V. All rights reserved.},
705    annote = {633ZD Times Cited:0 Cited References Count:9},
706    issn = {0168-9274},
707    uri = {<Go to ISI>://000180314200002},
# Line 727 | Line 710 | Encoding: GBK
710   @ARTICLE{Cheung2004,
711    author = {D. L. Cheung and S. J. Clark and M. R. Wilson},
712    title = {Calculation of flexoelectric coefficients for a nematic liquid crystal
713 <    by atomistic simulation},
713 >        by atomistic simulation},
714    journal = {Journal of Chemical Physics},
715    year = {2004},
716    volume = {121},
# Line 735 | Line 718 | Encoding: GBK
718    number = {18},
719    month = {Nov 8},
720    abstract = {Equilibrium molecular dynamics calculations have been performed for
721 <    the liquid crystal molecule n-4-(trans-4-n-pentylcyclohexyl)benzonitrile
722 <    (PCH5) using a fully atomistic model. Simulation data have been
723 <    obtained for a series of temperatures in the nematic phase. The
724 <    simulation data have been used to calculate the flexoelectric coefficients
725 <    e(s) and e(b) using the linear response formalism of Osipov and
726 <    Nemtsov [M. A. Osipov and V. B. Nemtsov, Sov. Phys. Crstallogr.
727 <    31, 125 (1986)]. The temperature and order parameter dependence
728 <    of e(s) and e(b) are examined, as are separate contributions from
729 <    different intermolecular interactions. Values of e(s) and e(b) calculated
730 <    from simulation are consistent with those found from experiment.
731 <    (C) 2004 American Institute of Physics.},
721 >        the liquid crystal molecule n-4-(trans-4-n-pentylcyclohexyl)benzonitrile
722 >        (PCH5) using a fully atomistic model. Simulation data have been
723 >        obtained for a series of temperatures in the nematic phase. The
724 >        simulation data have been used to calculate the flexoelectric coefficients
725 >        e(s) and e(b) using the linear response formalism of Osipov and
726 >        Nemtsov [M. A. Osipov and V. B. Nemtsov, Sov. Phys. Crstallogr.
727 >        31, 125 (1986)]. The temperature and order parameter dependence
728 >        of e(s) and e(b) are examined, as are separate contributions from
729 >        different intermolecular interactions. Values of e(s) and e(b) calculated
730 >        from simulation are consistent with those found from experiment.
731 >        (C) 2004 American Institute of Physics.},
732    annote = {866UM Times Cited:4 Cited References Count:61},
733    issn = {0021-9606},
734    uri = {<Go to ISI>://000224798900053},
# Line 761 | Line 744 | Encoding: GBK
744    number = {1-2},
745    month = {Apr 15},
746    abstract = {Equilibrium molecular dynamics calculations have been performed for
747 <    the liquid crystal molecule n-4-(trans-4-npentylcyclohexyl)benzonitrile
748 <    (PCH5) using a fully atomistic model. Simulation data has been obtained
749 <    for a series of temperatures in the nematic phase. The rotational
750 <    viscosity co-efficient gamma(1), has been calculated using the angular
751 <    velocity correlation function of the nematic director, n, the mean
752 <    squared diffusion of n and statistical mechanical methods based
753 <    on the rotational diffusion co-efficient. We find good agreement
754 <    between the first two methods and experimental values. (C) 2002
755 <    Published by Elsevier Science B.V.},
747 >        the liquid crystal molecule n-4-(trans-4-npentylcyclohexyl)benzonitrile
748 >        (PCH5) using a fully atomistic model. Simulation data has been obtained
749 >        for a series of temperatures in the nematic phase. The rotational
750 >        viscosity co-efficient gamma(1), has been calculated using the angular
751 >        velocity correlation function of the nematic director, n, the mean
752 >        squared diffusion of n and statistical mechanical methods based
753 >        on the rotational diffusion co-efficient. We find good agreement
754 >        between the first two methods and experimental values. (C) 2002
755 >        Published by Elsevier Science B.V.},
756    annote = {547KF Times Cited:8 Cited References Count:31},
757    issn = {0009-2614},
758    uri = {<Go to ISI>://000175331000020},
# Line 778 | Line 761 | Encoding: GBK
761   @ARTICLE{Chin2004,
762    author = {S. A. Chin},
763    title = {Dynamical multiple-time stepping methods for overcoming resonance
764 <    instabilities},
764 >        instabilities},
765    journal = {Journal of Chemical Physics},
766    year = {2004},
767    volume = {120},
# Line 786 | Line 769 | Encoding: GBK
769    number = {1},
770    month = {Jan 1},
771    abstract = {Current molecular dynamics simulations of biomolecules using multiple
772 <    time steps to update the slowly changing force are hampered by instabilities
773 <    beginning at time steps near the half period of the fastest vibrating
774 <    mode. These #resonance# instabilities have became a critical barrier
775 <    preventing the long time simulation of biomolecular dynamics. Attempts
776 <    to tame these instabilities by altering the slowly changing force
777 <    and efforts to damp them out by Langevin dynamics do not address
778 <    the fundamental cause of these instabilities. In this work, we trace
779 <    the instability to the nonanalytic character of the underlying spectrum
780 <    and show that a correct splitting of the Hamiltonian, which renders
781 <    the spectrum analytic, restores stability. The resulting Hamiltonian
782 <    dictates that in addition to updating the momentum due to the slowly
783 <    changing force, one must also update the position with a modified
784 <    mass. Thus multiple-time stepping must be done dynamically. (C)
785 <    2004 American Institute of Physics.},
772 >        time steps to update the slowly changing force are hampered by instabilities
773 >        beginning at time steps near the half period of the fastest vibrating
774 >        mode. These #resonance# instabilities have became a critical barrier
775 >        preventing the long time simulation of biomolecular dynamics. Attempts
776 >        to tame these instabilities by altering the slowly changing force
777 >        and efforts to damp them out by Langevin dynamics do not address
778 >        the fundamental cause of these instabilities. In this work, we trace
779 >        the instability to the nonanalytic character of the underlying spectrum
780 >        and show that a correct splitting of the Hamiltonian, which renders
781 >        the spectrum analytic, restores stability. The resulting Hamiltonian
782 >        dictates that in addition to updating the momentum due to the slowly
783 >        changing force, one must also update the position with a modified
784 >        mass. Thus multiple-time stepping must be done dynamically. (C)
785 >        2004 American Institute of Physics.},
786    annote = {757TK Times Cited:1 Cited References Count:22},
787    issn = {0021-9606},
788    uri = {<Go to ISI>://000187577400003},
# Line 808 | Line 791 | Encoding: GBK
791   @ARTICLE{Cook2000,
792    author = {M. J. Cook and M. R. Wilson},
793    title = {Simulation studies of dipole correlation in the isotropic liquid
794 <    phase},
794 >        phase},
795    journal = {Liquid Crystals},
796    year = {2000},
797    volume = {27},
# Line 816 | Line 799 | Encoding: GBK
799    number = {12},
800    month = {Dec},
801    abstract = {The Kirkwood correlation factor g(1) determines the preference for
802 <    local parallel or antiparallel dipole association in the isotropic
803 <    phase. Calamitic mesogens with longitudinal dipole moments and Kirkwood
804 <    factors greater than 1 have an enhanced effective dipole moment
805 <    along the molecular long axis. This leads to higher values of Delta
806 <    epsilon in the nematic phase. This paper describes state-of-the-art
807 <    molecular dynamics simulations of two calamitic mesogens 4-(trans-4-n-pentylcyclohexyl)benzonitrile
808 <    (PCH5) and 4-(trans-4-n-pentylcyclohexyl) chlorobenzene (PCH5-Cl)
809 <    in the isotropic liquid phase using an all-atom force field and
810 <    taking long range electrostatics into account using an Ewald summation.
811 <    Using this methodology, PCH5 is seen to prefer antiparallel dipole
812 <    alignment with a negative g(1) and PCH5-Cl is seen to prefer parallel
813 <    dipole alignment with a positive g(1); this is in accordance with
814 <    experimental dielectric measurements. Analysis of the molecular
815 <    dynamics trajectories allows an assessment of why these molecules
816 <    behave differently.},
802 >        local parallel or antiparallel dipole association in the isotropic
803 >        phase. Calamitic mesogens with longitudinal dipole moments and Kirkwood
804 >        factors greater than 1 have an enhanced effective dipole moment
805 >        along the molecular long axis. This leads to higher values of Delta
806 >        epsilon in the nematic phase. This paper describes state-of-the-art
807 >        molecular dynamics simulations of two calamitic mesogens 4-(trans-4-n-pentylcyclohexyl)benzonitrile
808 >        (PCH5) and 4-(trans-4-n-pentylcyclohexyl) chlorobenzene (PCH5-Cl)
809 >        in the isotropic liquid phase using an all-atom force field and
810 >        taking long range electrostatics into account using an Ewald summation.
811 >        Using this methodology, PCH5 is seen to prefer antiparallel dipole
812 >        alignment with a negative g(1) and PCH5-Cl is seen to prefer parallel
813 >        dipole alignment with a positive g(1); this is in accordance with
814 >        experimental dielectric measurements. Analysis of the molecular
815 >        dynamics trajectories allows an assessment of why these molecules
816 >        behave differently.},
817    annote = {376BF Times Cited:10 Cited References Count:16},
818    issn = {0267-8292},
819    uri = {<Go to ISI>://000165437800002},
# Line 839 | Line 822 | Encoding: GBK
822   @ARTICLE{Cui2003,
823    author = {B. X. Cui and M. Y. Shen and K. F. Freed},
824    title = {Folding and misfolding of the papillomavirus E6 interacting peptide
825 <    E6ap},
825 >        E6ap},
826    journal = {Proceedings of the National Academy of Sciences of the United States
827 <    of America},
827 >        of America},
828    year = {2003},
829    volume = {100},
830    pages = {7087-7092},
831    number = {12},
832    month = {Jun 10},
833    abstract = {All-atom Langevin dynamics simulations have been performed to study
834 <    the folding pathways of the 18-residue binding domain fragment E6ap
835 <    of the human papillomavirus E6 interacting peptide. Six independent
836 <    folding trajectories, with a total duration of nearly 2 mus, all
837 <    lead to the same native state in which the E6ap adopts a fluctuating
838 <    a-helix structure in the central portion (Ser-4-Leu-13) but with
839 <    very flexible N and C termini. Simulations starting from different
840 <    core configurations exhibit the E6ap folding dynamics as either
841 <    a two- or three-state folder with an intermediate misfolded state.
842 <    The essential leucine hydrophobic core (Leu-9, Leu-12, and Leu-13)
843 <    is well conserved in the native-state structure but absent in the
844 <    intermediate structure, suggesting that the leucine core is not
845 <    only essential for the binding activity of E6ap but also important
846 <    for the stability of the native structure. The free energy landscape
847 <    reveals a significant barrier between the basins separating the
848 <    native and misfolded states. We also discuss the various underlying
849 <    forces that drive the peptide into its native state.},
834 >        the folding pathways of the 18-residue binding domain fragment E6ap
835 >        of the human papillomavirus E6 interacting peptide. Six independent
836 >        folding trajectories, with a total duration of nearly 2 mus, all
837 >        lead to the same native state in which the E6ap adopts a fluctuating
838 >        a-helix structure in the central portion (Ser-4-Leu-13) but with
839 >        very flexible N and C termini. Simulations starting from different
840 >        core configurations exhibit the E6ap folding dynamics as either
841 >        a two- or three-state folder with an intermediate misfolded state.
842 >        The essential leucine hydrophobic core (Leu-9, Leu-12, and Leu-13)
843 >        is well conserved in the native-state structure but absent in the
844 >        intermediate structure, suggesting that the leucine core is not
845 >        only essential for the binding activity of E6ap but also important
846 >        for the stability of the native structure. The free energy landscape
847 >        reveals a significant barrier between the basins separating the
848 >        native and misfolded states. We also discuss the various underlying
849 >        forces that drive the peptide into its native state.},
850    annote = {689LC Times Cited:3 Cited References Count:48},
851    issn = {0027-8424},
852    uri = {<Go to ISI>://000183493500037},
# Line 879 | Line 862 | Encoding: GBK
862    number = {1},
863    month = {Jan 1},
864    abstract = {We study the slow phase of thermally activated magnetic relaxation
865 <    in finite two-dimensional ensembles of dipolar interacting ferromagnetic
866 <    nanoparticles whose easy axes of magnetization are perpendicular
867 <    to the distribution plane. We develop a method to numerically simulate
868 <    the magnetic relaxation for the case that the smallest heights of
869 <    the potential barriers between the equilibrium directions of the
870 <    nanoparticle magnetic moments are much larger than the thermal energy.
871 <    Within this framework, we analyze in detail the role that the correlations
872 <    of the nanoparticle magnetic moments and the finite size of the
873 <    nanoparticle ensemble play in magnetic relaxation.},
865 >        in finite two-dimensional ensembles of dipolar interacting ferromagnetic
866 >        nanoparticles whose easy axes of magnetization are perpendicular
867 >        to the distribution plane. We develop a method to numerically simulate
868 >        the magnetic relaxation for the case that the smallest heights of
869 >        the potential barriers between the equilibrium directions of the
870 >        nanoparticle magnetic moments are much larger than the thermal energy.
871 >        Within this framework, we analyze in detail the role that the correlations
872 >        of the nanoparticle magnetic moments and the finite size of the
873 >        nanoparticle ensemble play in magnetic relaxation.},
874    annote = {642XH Times Cited:11 Cited References Count:31},
875    issn = {1098-0121},
876    uri = {<Go to ISI>://000180830400056},
# Line 903 | Line 886 | Encoding: GBK
886    number = {1},
887    month = {Jan},
888    abstract = {To explore the origin of the large-scale motion of triosephosphate
889 <    isomerase's flexible loop (residues 166 to 176) at the active site,
890 <    several simulation protocols are employed both for the free enzyme
891 <    in vacuo and for the free enzyme with some solvent modeling: high-temperature
892 <    Langevin dynamics simulations, sampling by a #dynamics##driver#
893 <    approach, and potential-energy surface calculations. Our focus is
894 <    on obtaining the energy barrier to the enzyme's motion and establishing
895 <    the nature of the loop movement. Previous calculations did not determine
896 <    this energy barrier and the effect of solvent on the barrier. High-temperature
897 <    molecular dynamics simulations and crystallographic studies have
898 <    suggested a rigid-body motion with two hinges located at both ends
899 <    of the loop; Brownian dynamics simulations at room temperature pointed
900 <    to a very flexible behavior. The present simulations and analyses
901 <    reveal that although solute/solvent hydrogen bonds play a crucial
902 <    role in lowering the energy along the pathway, there still remains
903 <    a high activation barrier, This finding clearly indicates that,
904 <    if the loop opens and closes in the absence of a substrate at standard
905 <    conditions (e.g., room temperature, appropriate concentration of
906 <    isomerase), the time scale for transition is not in the nanosecond
907 <    but rather the microsecond range. Our results also indicate that
908 <    in the context of spontaneous opening in the free enzyme, the motion
909 <    is of rigid-body type and that the specific interaction between
910 <    residues Ala(176) and Tyr(208) plays a crucial role in the loop
911 <    opening/closing mechanism.},
889 >        isomerase's flexible loop (residues 166 to 176) at the active site,
890 >        several simulation protocols are employed both for the free enzyme
891 >        in vacuo and for the free enzyme with some solvent modeling: high-temperature
892 >        Langevin dynamics simulations, sampling by a #dynamics##driver#
893 >        approach, and potential-energy surface calculations. Our focus is
894 >        on obtaining the energy barrier to the enzyme's motion and establishing
895 >        the nature of the loop movement. Previous calculations did not determine
896 >        this energy barrier and the effect of solvent on the barrier. High-temperature
897 >        molecular dynamics simulations and crystallographic studies have
898 >        suggested a rigid-body motion with two hinges located at both ends
899 >        of the loop; Brownian dynamics simulations at room temperature pointed
900 >        to a very flexible behavior. The present simulations and analyses
901 >        reveal that although solute/solvent hydrogen bonds play a crucial
902 >        role in lowering the energy along the pathway, there still remains
903 >        a high activation barrier, This finding clearly indicates that,
904 >        if the loop opens and closes in the absence of a substrate at standard
905 >        conditions (e.g., room temperature, appropriate concentration of
906 >        isomerase), the time scale for transition is not in the nanosecond
907 >        but rather the microsecond range. Our results also indicate that
908 >        in the context of spontaneous opening in the free enzyme, the motion
909 >        is of rigid-body type and that the specific interaction between
910 >        residues Ala(176) and Tyr(208) plays a crucial role in the loop
911 >        opening/closing mechanism.},
912    annote = {Zl046 Times Cited:30 Cited References Count:29},
913    issn = {0006-3495},
914    uri = {<Go to ISI>://000073393400009},
# Line 941 | Line 924 | Encoding: GBK
924    number = {15},
925    month = {Oct 15},
926    abstract = {Rigid body molecular models possess symplectic structure and time-reversal
927 <    symmetry. Standard numerical integration methods destroy both properties,
928 <    introducing nonphysical dynamical behavior such as numerically induced
929 <    dissipative states and drift in the energy during long term simulations.
930 <    This article describes the construction, implementation, and practical
931 <    application of fast explicit symplectic-reversible integrators for
932 <    multiple rigid body molecular simulations, These methods use a reduction
933 <    to Euler equations for the free rigid body, together with a symplectic
934 <    splitting technique. In every time step, the orientational dynamics
935 <    of each rigid body is integrated by a sequence of planar rotations.
936 <    Besides preserving the symplectic and reversible structures of the
937 <    flow, this scheme accurately conserves the total angular momentum
938 <    of a system of interacting rigid bodies. Excellent energy conservation
939 <    fan be obtained relative to traditional methods, especially in long-time
940 <    simulations. The method is implemented in a research code, ORIENT
941 <    and compared with a quaternion/extrapolation scheme for the TIP4P
942 <    model of water. Our experiments show that the symplectic-reversible
943 <    scheme is far superior to the more traditional quaternion method.
944 <    (C) 1997 American Institute of Physics.},
927 >        symmetry. Standard numerical integration methods destroy both properties,
928 >        introducing nonphysical dynamical behavior such as numerically induced
929 >        dissipative states and drift in the energy during long term simulations.
930 >        This article describes the construction, implementation, and practical
931 >        application of fast explicit symplectic-reversible integrators for
932 >        multiple rigid body molecular simulations, These methods use a reduction
933 >        to Euler equations for the free rigid body, together with a symplectic
934 >        splitting technique. In every time step, the orientational dynamics
935 >        of each rigid body is integrated by a sequence of planar rotations.
936 >        Besides preserving the symplectic and reversible structures of the
937 >        flow, this scheme accurately conserves the total angular momentum
938 >        of a system of interacting rigid bodies. Excellent energy conservation
939 >        fan be obtained relative to traditional methods, especially in long-time
940 >        simulations. The method is implemented in a research code, ORIENT
941 >        and compared with a quaternion/extrapolation scheme for the TIP4P
942 >        model of water. Our experiments show that the symplectic-reversible
943 >        scheme is far superior to the more traditional quaternion method.
944 >        (C) 1997 American Institute of Physics.},
945    annote = {Ya587 Times Cited:35 Cited References Count:32},
946    issn = {0021-9606},
947    uri = {<Go to ISI>://A1997YA58700024},
# Line 967 | Line 950 | Encoding: GBK
950   @ARTICLE{Edwards2005,
951    author = {S. A. Edwards and D. R. M. Williams},
952    title = {Stretching a single diblock copolymer in a selective solvent: Langevin
953 <    dynamics simulations},
953 >        dynamics simulations},
954    journal = {Macromolecules},
955    year = {2005},
956    volume = {38},
# Line 975 | Line 958 | Encoding: GBK
958    number = {25},
959    month = {Dec 13},
960    abstract = {Using the Langevin dynamics technique, we have carried out simulations
961 <    of a single-chain flexible diblock copolymer. The polymer consists
962 <    of two blocks of equal length, one very poorly solvated and the
963 <    other close to theta-conditions. We study what happens when such
964 <    a polymer is stretched, for a range of different stretching speeds,
965 <    and correlate our observations with features in the plot of force
966 <    vs extension. We find that at slow speeds this force profile does
967 <    not increase monotonically, in disagreement with earlier predictions,
968 <    and that at high speeds there is a strong dependence on which end
969 <    of the polymer is pulled, as well as a high level of hysteresis.},
961 >        of a single-chain flexible diblock copolymer. The polymer consists
962 >        of two blocks of equal length, one very poorly solvated and the
963 >        other close to theta-conditions. We study what happens when such
964 >        a polymer is stretched, for a range of different stretching speeds,
965 >        and correlate our observations with features in the plot of force
966 >        vs extension. We find that at slow speeds this force profile does
967 >        not increase monotonically, in disagreement with earlier predictions,
968 >        and that at high speeds there is a strong dependence on which end
969 >        of the polymer is pulled, as well as a high level of hysteresis.},
970    annote = {992EC Times Cited:0 Cited References Count:13},
971    issn = {0024-9297},
972    uri = {<Go to ISI>://000233866200035},
# Line 992 | Line 975 | Encoding: GBK
975   @ARTICLE{Egberts1988,
976    author = {E. Egberts and H. J. C. Berendsen},
977    title = {Molecular-Dynamics Simulation of a Smectic Liquid-Crystal with Atomic
978 <    Detail},
978 >        Detail},
979    journal = {Journal of Chemical Physics},
980    year = {1988},
981    volume = {89},
# Line 1017 | Line 1000 | Encoding: GBK
1000    uri = {<Go to ISI>://A1978FP21600004},
1001   }
1002  
1003 + @ARTICLE{Evans1977,
1004 +  author = {D. J. Evans},
1005 +  title = {Representation of Orientation Space},
1006 +  journal = {Molecular Physics},
1007 +  year = {1977},
1008 +  volume = {34},
1009 +  pages = {317-325},
1010 +  number = {2},
1011 +  annote = {Ds757 Times Cited:271 Cited References Count:18},
1012 +  issn = {0026-8976},
1013 +  uri = {<Go to ISI>://A1977DS75700002},
1014 + }
1015 +
1016   @ARTICLE{Fennell2004,
1017    author = {C. J. Fennell and J. D. Gezelter},
1018    title = {On the structural and transport properties of the soft sticky dipole
1019 <    and related single-point water models},
1019 >        and related single-point water models},
1020    journal = {Journal of Chemical Physics},
1021    year = {2004},
1022    volume = {120},
# Line 1028 | Line 1024 | Encoding: GBK
1024    number = {19},
1025    month = {May 15},
1026    abstract = {The density maximum and temperature dependence of the self-diffusion
1027 <    constant were investigated for the soft sticky dipole (SSD) water
1028 <    model and two related reparametrizations of this single-point model.
1029 <    A combination of microcanonical and isobaric-isothermal molecular
1030 <    dynamics simulations was used to calculate these properties, both
1031 <    with and without the use of reaction field to handle long-range
1032 <    electrostatics. The isobaric-isothermal simulations of the melting
1033 <    of both ice-I-h and ice-I-c showed a density maximum near 260 K.
1034 <    In most cases, the use of the reaction field resulted in calculated
1035 <    densities which were significantly lower than experimental densities.
1036 <    Analysis of self-diffusion constants shows that the original SSD
1037 <    model captures the transport properties of experimental water very
1038 <    well in both the normal and supercooled liquid regimes. We also
1039 <    present our reparametrized versions of SSD for use both with the
1040 <    reaction field or without any long-range electrostatic corrections.
1041 <    These are called the SSD/RF and SSD/E models, respectively. These
1042 <    modified models were shown to maintain or improve upon the experimental
1043 <    agreement with the structural and transport properties that can
1044 <    be obtained with either the original SSD or the density-corrected
1045 <    version of the original model (SSD1). Additionally, a novel low-density
1046 <    ice structure is presented which appears to be the most stable ice
1047 <    structure for the entire SSD family. (C) 2004 American Institute
1048 <    of Physics.},
1027 >        constant were investigated for the soft sticky dipole (SSD) water
1028 >        model and two related reparametrizations of this single-point model.
1029 >        A combination of microcanonical and isobaric-isothermal molecular
1030 >        dynamics simulations was used to calculate these properties, both
1031 >        with and without the use of reaction field to handle long-range
1032 >        electrostatics. The isobaric-isothermal simulations of the melting
1033 >        of both ice-I-h and ice-I-c showed a density maximum near 260 K.
1034 >        In most cases, the use of the reaction field resulted in calculated
1035 >        densities which were significantly lower than experimental densities.
1036 >        Analysis of self-diffusion constants shows that the original SSD
1037 >        model captures the transport properties of experimental water very
1038 >        well in both the normal and supercooled liquid regimes. We also
1039 >        present our reparametrized versions of SSD for use both with the
1040 >        reaction field or without any long-range electrostatic corrections.
1041 >        These are called the SSD/RF and SSD/E models, respectively. These
1042 >        modified models were shown to maintain or improve upon the experimental
1043 >        agreement with the structural and transport properties that can
1044 >        be obtained with either the original SSD or the density-corrected
1045 >        version of the original model (SSD1). Additionally, a novel low-density
1046 >        ice structure is presented which appears to be the most stable ice
1047 >        structure for the entire SSD family. (C) 2004 American Institute
1048 >        of Physics.},
1049    annote = {816YY Times Cited:5 Cited References Count:39},
1050    issn = {0021-9606},
1051    uri = {<Go to ISI>://000221146400032},
# Line 1058 | Line 1054 | Encoding: GBK
1054   @ARTICLE{Fernandes2002,
1055    author = {M. X. Fernandes and J. G. {de la Torre}},
1056    title = {Brownian dynamics simulation of rigid particles of arbitrary shape
1057 <    in external fields},
1057 >        in external fields},
1058    journal = {Biophysical Journal},
1059    year = {2002},
1060    volume = {83},
# Line 1066 | Line 1062 | Encoding: GBK
1062    number = {6},
1063    month = {Dec},
1064    abstract = {We have developed a Brownian dynamics simulation algorithm to generate
1065 <    Brownian trajectories of an isolated, rigid particle of arbitrary
1066 <    shape in the presence of electric fields or any other external agents.
1067 <    Starting from the generalized diffusion tensor, which can be calculated
1068 <    with the existing HYDRO software, the new program BROWNRIG (including
1069 <    a case-specific subprogram for the external agent) carries out a
1070 <    simulation that is analyzed later to extract the observable dynamic
1071 <    properties. We provide a variety of examples of utilization of this
1072 <    method, which serve as tests of its performance, and also illustrate
1073 <    its applicability. Examples include free diffusion, transport in
1074 <    an electric field, and diffusion in a restricting environment.},
1065 >        Brownian trajectories of an isolated, rigid particle of arbitrary
1066 >        shape in the presence of electric fields or any other external agents.
1067 >        Starting from the generalized diffusion tensor, which can be calculated
1068 >        with the existing HYDRO software, the new program BROWNRIG (including
1069 >        a case-specific subprogram for the external agent) carries out a
1070 >        simulation that is analyzed later to extract the observable dynamic
1071 >        properties. We provide a variety of examples of utilization of this
1072 >        method, which serve as tests of its performance, and also illustrate
1073 >        its applicability. Examples include free diffusion, transport in
1074 >        an electric field, and diffusion in a restricting environment.},
1075    annote = {633AD Times Cited:2 Cited References Count:43},
1076    issn = {0006-3495},
1077    uri = {<Go to ISI>://000180256300012},
# Line 1084 | Line 1080 | Encoding: GBK
1080   @ARTICLE{Gay1981,
1081    author = {J. G. Gay and B. J. Berne},
1082    title = {Modification of the Overlap Potential to Mimic a Linear Site-Site
1083 <    Potential},
1083 >        Potential},
1084    journal = {Journal of Chemical Physics},
1085    year = {1981},
1086    volume = {74},
# Line 1105 | Line 1101 | Encoding: GBK
1101    number = {6},
1102    month = {Nov},
1103    abstract = {To investigate the influence of inertial effects on the dynamics of
1104 <    an assembly of beads subjected to rigid constraints and placed in
1105 <    a buffer medium, a convenient method to introduce suitable generalized
1106 <    coordinates is presented. Without any restriction on the nature
1107 <    of the soft forces involved (both stochastic and deterministic),
1108 <    pertinent Langevin equations are derived. Provided that the Brownian
1109 <    forces are Gaussian and Markovian, the corresponding Fokker-Planck
1110 <    equation (FPE) is obtained in the complete phase space of generalized
1111 <    coordinates and momenta. The correct short time behavior for correlation
1112 <    functions (CFs) of generalized coordinates is established, and the
1113 <    diffusion equation with memory (DEM) is deduced from the FPE in
1114 <    the high friction Limit. The DEM is invoked to perform illustrative
1115 <    calculations in two dimensions of the orientational CFs for once
1116 <    broken nonrigid rods immobilized on a surface. These calculations
1117 <    reveal that the CFs under certain conditions exhibit an oscillatory
1118 <    behavior, which is irreproducible within the standard diffusion
1119 <    equation. Several methods are considered for the approximate solution
1120 <    of the DEM, and their application to three dimensional DEMs is discussed.},
1104 >        an assembly of beads subjected to rigid constraints and placed in
1105 >        a buffer medium, a convenient method to introduce suitable generalized
1106 >        coordinates is presented. Without any restriction on the nature
1107 >        of the soft forces involved (both stochastic and deterministic),
1108 >        pertinent Langevin equations are derived. Provided that the Brownian
1109 >        forces are Gaussian and Markovian, the corresponding Fokker-Planck
1110 >        equation (FPE) is obtained in the complete phase space of generalized
1111 >        coordinates and momenta. The correct short time behavior for correlation
1112 >        functions (CFs) of generalized coordinates is established, and the
1113 >        diffusion equation with memory (DEM) is deduced from the FPE in
1114 >        the high friction Limit. The DEM is invoked to perform illustrative
1115 >        calculations in two dimensions of the orientational CFs for once
1116 >        broken nonrigid rods immobilized on a surface. These calculations
1117 >        reveal that the CFs under certain conditions exhibit an oscillatory
1118 >        behavior, which is irreproducible within the standard diffusion
1119 >        equation. Several methods are considered for the approximate solution
1120 >        of the DEM, and their application to three dimensional DEMs is discussed.},
1121    annote = {257MM Times Cited:2 Cited References Count:82},
1122    issn = {1022-1344},
1123    uri = {<Go to ISI>://000083785700002},
1124   }
1125  
1130 @BOOK{Goldstein2001,
1131  title = {Classical Mechanics},
1132  publisher = {Addison Wesley},
1133  year = {2001},
1134  author = {H. Goldstein and C. Poole and J. Safko},
1135  address = {San Francisco},
1136  edition = {3rd},
1137 }
1138
1126   @ARTICLE{Gray2003,
1127    author = {J. J. Gray and S. Moughon and C. Wang and O. Schueler-Furman and
1128 <    B. Kuhlman and C. A. Rohl and D. Baker},
1128 >        B. Kuhlman and C. A. Rohl and D. Baker},
1129    title = {Protein-protein docking with simultaneous optimization of rigid-body
1130 <    displacement and side-chain conformations},
1130 >        displacement and side-chain conformations},
1131    journal = {Journal of Molecular Biology},
1132    year = {2003},
1133    volume = {331},
# Line 1148 | Line 1135 | Encoding: GBK
1135    number = {1},
1136    month = {Aug 1},
1137    abstract = {Protein-protein docking algorithms provide a means to elucidate structural
1138 <    details for presently unknown complexes. Here, we present and evaluate
1139 <    a new method to predict protein-protein complexes from the coordinates
1140 <    of the unbound monomer components. The method employs a low-resolution,
1141 <    rigid-body, Monte Carlo search followed by simultaneous optimization
1142 <    of backbone displacement and side-chain conformations using Monte
1143 <    Carlo minimization. Up to 10(5) independent simulations are carried
1144 <    out, and the resulting #decoys# are ranked using an energy function
1145 <    dominated by van der Waals interactions, an implicit solvation model,
1146 <    and an orientation-dependent hydrogen bonding potential. Top-ranking
1147 <    decoys are clustered to select the final predictions. Small-perturbation
1148 <    studies reveal the formation of binding funnels in 42 of 54 cases
1149 <    using coordinates derived from the bound complexes and in 32 of
1150 <    54 cases using independently determined coordinates of one or both
1151 <    monomers. Experimental binding affinities correlate with the calculated
1152 <    score function and explain the predictive success or failure of
1153 <    many targets. Global searches using one or both unbound components
1154 <    predict at least 25% of the native residue-residue contacts in 28
1155 <    of the 32 cases where binding funnels exist. The results suggest
1156 <    that the method may soon be useful for generating models of biologically
1157 <    important complexes from the structures of the isolated components,
1158 <    but they also highlight the challenges that must be met to achieve
1159 <    consistent and accurate prediction of protein-protein interactions.
1160 <    (C) 2003 Elsevier Ltd. All rights reserved.},
1138 >        details for presently unknown complexes. Here, we present and evaluate
1139 >        a new method to predict protein-protein complexes from the coordinates
1140 >        of the unbound monomer components. The method employs a low-resolution,
1141 >        rigid-body, Monte Carlo search followed by simultaneous optimization
1142 >        of backbone displacement and side-chain conformations using Monte
1143 >        Carlo minimization. Up to 10(5) independent simulations are carried
1144 >        out, and the resulting #decoys# are ranked using an energy function
1145 >        dominated by van der Waals interactions, an implicit solvation model,
1146 >        and an orientation-dependent hydrogen bonding potential. Top-ranking
1147 >        decoys are clustered to select the final predictions. Small-perturbation
1148 >        studies reveal the formation of binding funnels in 42 of 54 cases
1149 >        using coordinates derived from the bound complexes and in 32 of
1150 >        54 cases using independently determined coordinates of one or both
1151 >        monomers. Experimental binding affinities correlate with the calculated
1152 >        score function and explain the predictive success or failure of
1153 >        many targets. Global searches using one or both unbound components
1154 >        predict at least 25% of the native residue-residue contacts in 28
1155 >        of the 32 cases where binding funnels exist. The results suggest
1156 >        that the method may soon be useful for generating models of biologically
1157 >        important complexes from the structures of the isolated components,
1158 >        but they also highlight the challenges that must be met to achieve
1159 >        consistent and accurate prediction of protein-protein interactions.
1160 >        (C) 2003 Elsevier Ltd. All rights reserved.},
1161    annote = {704QL Times Cited:48 Cited References Count:60},
1162    issn = {0022-2836},
1163    uri = {<Go to ISI>://000184351300022},
# Line 1186 | Line 1173 | Encoding: GBK
1173    number = {5174},
1174    month = {Aug 12},
1175    abstract = {Some of the recently developed fast summation methods that have arisen
1176 <    in scientific computing are described. These methods require an
1177 <    amount of work proportional to N or N log N to evaluate all pairwise
1178 <    interactions in an ensemble of N particles. Traditional methods,
1179 <    by contrast, require an amount of work proportional to N-2. AS a
1180 <    result, large-scale simulations can be carried out using only modest
1181 <    computer resources. In combination with supercomputers, it is possible
1182 <    to address questions that were previously out of reach. Problems
1183 <    from diffusion, gravitation, and wave propagation are considered.},
1176 >        in scientific computing are described. These methods require an
1177 >        amount of work proportional to N or N log N to evaluate all pairwise
1178 >        interactions in an ensemble of N particles. Traditional methods,
1179 >        by contrast, require an amount of work proportional to N-2. AS a
1180 >        result, large-scale simulations can be carried out using only modest
1181 >        computer resources. In combination with supercomputers, it is possible
1182 >        to address questions that were previously out of reach. Problems
1183 >        from diffusion, gravitation, and wave propagation are considered.},
1184    annote = {Pb499 Times Cited:99 Cited References Count:44},
1185    issn = {0036-8075},
1186    uri = {<Go to ISI>://A1994PB49900031},
# Line 1223 | Line 1210 | Encoding: GBK
1210    number = {4},
1211    month = {Jun},
1212    abstract = {Backward error analysis is a useful tool for the study of numerical
1213 <    approximations to ordinary differential equations. The numerical
1214 <    solution is formally interpreted as the exact solution of a perturbed
1215 <    differential equation, given as a formal and usually divergent series
1216 <    in powers of the step size. For a rigorous analysis, this series
1217 <    has to be truncated. In this article we study the influence of this
1218 <    truncation to the difference between the numerical solution and
1219 <    the exact solution of the perturbed differential equation. Results
1220 <    on the long-time behaviour of numerical solutions are obtained in
1221 <    this way. We present applications to the numerical phase portrait
1222 <    near hyperbolic equilibrium points, to asymptotically stable periodic
1223 <    orbits and Hopf bifurcation, and to energy conservation and approximation
1224 <    of invariant tori in Hamiltonian systems.},
1213 >        approximations to ordinary differential equations. The numerical
1214 >        solution is formally interpreted as the exact solution of a perturbed
1215 >        differential equation, given as a formal and usually divergent series
1216 >        in powers of the step size. For a rigorous analysis, this series
1217 >        has to be truncated. In this article we study the influence of this
1218 >        truncation to the difference between the numerical solution and
1219 >        the exact solution of the perturbed differential equation. Results
1220 >        on the long-time behaviour of numerical solutions are obtained in
1221 >        this way. We present applications to the numerical phase portrait
1222 >        near hyperbolic equilibrium points, to asymptotically stable periodic
1223 >        orbits and Hopf bifurcation, and to energy conservation and approximation
1224 >        of invariant tori in Hamiltonian systems.},
1225    annote = {Xj488 Times Cited:50 Cited References Count:19},
1226    issn = {0029-599X},
1227    uri = {<Go to ISI>://A1997XJ48800002},
# Line 1243 | Line 1230 | Encoding: GBK
1230   @ARTICLE{Hao1993,
1231    author = {M. H. Hao and M. R. Pincus and S. Rackovsky and H. A. Scheraga},
1232    title = {Unfolding and Refolding of the Native Structure of Bovine Pancreatic
1233 <    Trypsin-Inhibitor Studied by Computer-Simulations},
1233 >        Trypsin-Inhibitor Studied by Computer-Simulations},
1234    journal = {Biochemistry},
1235    year = {1993},
1236    volume = {32},
# Line 1251 | Line 1238 | Encoding: GBK
1238    number = {37},
1239    month = {Sep 21},
1240    abstract = {A new procedure for studying the folding and unfolding of proteins,
1241 <    with an application to bovine pancreatic trypsin inhibitor (BPTI),
1242 <    is reported. The unfolding and refolding of the native structure
1243 <    of the protein are characterized by the dimensions of the protein,
1244 <    expressed in terms of the three principal radii of the structure
1245 <    considered as an ellipsoid. A dynamic equation, describing the variations
1246 <    of the principal radii on the unfolding path, and a numerical procedure
1247 <    to solve this equation are proposed. Expanded and distorted conformations
1248 <    are refolded to the native structure by a dimensional-constraint
1249 <    energy minimization procedure. A unique and reproducible unfolding
1250 <    pathway for an intermediate of BPTI lacking the [30,51] disulfide
1251 <    bond is obtained. The resulting unfolded conformations are extended;
1252 <    they contain near-native local structure, but their longest principal
1253 <    radii are more than 2.5 times greater than that of the native structure.
1254 <    The most interesting finding is that the majority of expanded conformations,
1255 <    generated under various conditions, can be refolded closely to the
1256 <    native structure, as measured by the correct overall chain fold,
1257 <    by the rms deviations from the native structure of only 1.9-3.1
1258 <    angstrom, and by the energy differences of about 10 kcal/mol from
1259 <    the native structure. Introduction of the [30,51] disulfide bond
1260 <    at this stage, followed by minimization, improves the closeness
1261 <    of the refolded structures to the native structure, reducing the
1262 <    rms deviations to 0.9-2.0 angstrom. The unique refolding of these
1263 <    expanded structures over such a large conformational space implies
1264 <    that the folding is strongly dictated by the interactions in the
1265 <    amino acid sequence of BPTI. The simulations indicate that, under
1266 <    conditions that favor a compact structure as mimicked by the volume
1267 <    constraints in our algorithm; the expanded conformations have a
1268 <    strong tendency to move toward the native structure; therefore,
1269 <    they probably would be favorable folding intermediates. The results
1270 <    presented here support a general model for protein folding, i.e.,
1271 <    progressive formation of partially folded structural units, followed
1272 <    by collapse to the compact native structure. The general applicability
1273 <    of the procedure is also discussed.},
1241 >        with an application to bovine pancreatic trypsin inhibitor (BPTI),
1242 >        is reported. The unfolding and refolding of the native structure
1243 >        of the protein are characterized by the dimensions of the protein,
1244 >        expressed in terms of the three principal radii of the structure
1245 >        considered as an ellipsoid. A dynamic equation, describing the variations
1246 >        of the principal radii on the unfolding path, and a numerical procedure
1247 >        to solve this equation are proposed. Expanded and distorted conformations
1248 >        are refolded to the native structure by a dimensional-constraint
1249 >        energy minimization procedure. A unique and reproducible unfolding
1250 >        pathway for an intermediate of BPTI lacking the [30,51] disulfide
1251 >        bond is obtained. The resulting unfolded conformations are extended;
1252 >        they contain near-native local structure, but their longest principal
1253 >        radii are more than 2.5 times greater than that of the native structure.
1254 >        The most interesting finding is that the majority of expanded conformations,
1255 >        generated under various conditions, can be refolded closely to the
1256 >        native structure, as measured by the correct overall chain fold,
1257 >        by the rms deviations from the native structure of only 1.9-3.1
1258 >        angstrom, and by the energy differences of about 10 kcal/mol from
1259 >        the native structure. Introduction of the [30,51] disulfide bond
1260 >        at this stage, followed by minimization, improves the closeness
1261 >        of the refolded structures to the native structure, reducing the
1262 >        rms deviations to 0.9-2.0 angstrom. The unique refolding of these
1263 >        expanded structures over such a large conformational space implies
1264 >        that the folding is strongly dictated by the interactions in the
1265 >        amino acid sequence of BPTI. The simulations indicate that, under
1266 >        conditions that favor a compact structure as mimicked by the volume
1267 >        constraints in our algorithm; the expanded conformations have a
1268 >        strong tendency to move toward the native structure; therefore,
1269 >        they probably would be favorable folding intermediates. The results
1270 >        presented here support a general model for protein folding, i.e.,
1271 >        progressive formation of partially folded structural units, followed
1272 >        by collapse to the compact native structure. The general applicability
1273 >        of the procedure is also discussed.},
1274    annote = {Ly294 Times Cited:27 Cited References Count:57},
1275    issn = {0006-2960},
1276    uri = {<Go to ISI>://A1993LY29400014},
# Line 1291 | Line 1278 | Encoding: GBK
1278  
1279   @ARTICLE{Hinsen2000,
1280    author = {K. Hinsen and A. J. Petrescu and S. Dellerue and M. C. Bellissent-Funel
1281 <    and G. R. Kneller},
1281 >        and G. R. Kneller},
1282    title = {Harmonicity in slow protein dynamics},
1283    journal = {Chemical Physics},
1284    year = {2000},
# Line 1300 | Line 1287 | Encoding: GBK
1287    number = {1-2},
1288    month = {Nov 1},
1289    abstract = {The slow dynamics of proteins around its native folded state is usually
1290 <    described by diffusion in a strongly anharmonic potential. In this
1291 <    paper, we try to understand the form and origin of the anharmonicities,
1292 <    with the principal aim of gaining a better understanding of the
1293 <    principal motion types, but also in order to develop more efficient
1294 <    numerical methods for simulating neutron scattering spectra of large
1295 <    proteins. First, we decompose a molecular dynamics (MD) trajectory
1296 <    of 1.5 ns for a C-phycocyanin dimer surrounded by a layer of water
1297 <    into three contributions that we expect to be independent: the global
1298 <    motion of the residues, the rigid-body motion of the sidechains
1299 <    relative to the backbone, and the internal deformations of the sidechains.
1300 <    We show that they are indeed almost independent by verifying the
1301 <    factorization of the incoherent intermediate scattering function.
1302 <    Then, we show that the global residue motions, which include all
1303 <    large-scale backbone motions, can be reproduced by a simple harmonic
1304 <    model which contains two contributions: a short-time vibrational
1305 <    term, described by a standard normal mode calculation in a local
1306 <    minimum, and a long-time diffusive term, described by Brownian motion
1307 <    in an effective harmonic potential. The potential and the friction
1308 <    constants were fitted to the MD data. The major anharmonic contribution
1309 <    to the incoherent intermediate scattering function comes from the
1310 <    rigid-body diffusion of the sidechains. This model can be used to
1311 <    calculate scattering functions for large proteins and for long-time
1312 <    scales very efficiently, and thus provides a useful complement to
1313 <    MD simulations, which are best suited for detailed studies on smaller
1314 <    systems or for shorter time scales. (C) 2000 Elsevier Science B.V.
1315 <    All rights reserved.},
1290 >        described by diffusion in a strongly anharmonic potential. In this
1291 >        paper, we try to understand the form and origin of the anharmonicities,
1292 >        with the principal aim of gaining a better understanding of the
1293 >        principal motion types, but also in order to develop more efficient
1294 >        numerical methods for simulating neutron scattering spectra of large
1295 >        proteins. First, we decompose a molecular dynamics (MD) trajectory
1296 >        of 1.5 ns for a C-phycocyanin dimer surrounded by a layer of water
1297 >        into three contributions that we expect to be independent: the global
1298 >        motion of the residues, the rigid-body motion of the sidechains
1299 >        relative to the backbone, and the internal deformations of the sidechains.
1300 >        We show that they are indeed almost independent by verifying the
1301 >        factorization of the incoherent intermediate scattering function.
1302 >        Then, we show that the global residue motions, which include all
1303 >        large-scale backbone motions, can be reproduced by a simple harmonic
1304 >        model which contains two contributions: a short-time vibrational
1305 >        term, described by a standard normal mode calculation in a local
1306 >        minimum, and a long-time diffusive term, described by Brownian motion
1307 >        in an effective harmonic potential. The potential and the friction
1308 >        constants were fitted to the MD data. The major anharmonic contribution
1309 >        to the incoherent intermediate scattering function comes from the
1310 >        rigid-body diffusion of the sidechains. This model can be used to
1311 >        calculate scattering functions for large proteins and for long-time
1312 >        scales very efficiently, and thus provides a useful complement to
1313 >        MD simulations, which are best suited for detailed studies on smaller
1314 >        systems or for shorter time scales. (C) 2000 Elsevier Science B.V.
1315 >        All rights reserved.},
1316    annote = {Sp. Iss. SI 368MT Times Cited:16 Cited References Count:31},
1317    issn = {0301-0104},
1318    uri = {<Go to ISI>://000090121700003},
# Line 1341 | Line 1328 | Encoding: GBK
1328    number = {4},
1329    month = {Oct},
1330    abstract = {Evidence has been found for the existence water at the protein-lipid
1331 <    hydrophobic interface ot the membrane proteins, gramicidin and apocytochrome
1332 <    C, using two related fluorescence spectroscopic approaches. The
1333 <    first approach exploited the fact that the presence of water in
1334 <    the excited state solvent cage of a fluorophore increases the rate
1335 <    of decay. For 1,6-diphenyl-1,3,5-hexatriene (DPH) and 1-palmitoyl-2-[[2-[4-(6-phenyl-trans-1,3,5-hexatrienyl)
1336 <    phenyl]ethyl]carbonyl]-3-sn-PC (DPH-PC), where the fluorophores
1337 <    are located in the hydrophobic core of the lipid bilayer, the introduction
1338 <    of gramicidin reduced the fluorescence lifetime, indicative of an
1339 <    increased presence of water in the bilayer. Since a high protein:lipid
1340 <    ratio was used, the fluorophores were forced to be adjacent to the
1341 <    protein hydrophobic surface, hence the presence of water in this
1342 <    region could be inferred. Cholesterol is known to reduce the water
1343 <    content of lipid bilayers and this effect was maintained at the
1344 <    protein-lipid interface with both gramicidin and apocytochrome C,
1345 <    again suggesting hydration in this region. The second approach was
1346 <    to use the fluorescence enhancement induced by exchanging deuterium
1347 <    oxide (D2O) for H2O. Both the fluorescence intensities of trimethylammonium-DPH,
1348 <    located in the lipid head group region, and of the gramicidin intrinsic
1349 <    tryptophans were greater in a D2O buffer compared with H2O, showing
1350 <    that the fluorophores were exposed to water in the bilayer at the
1351 <    protein-lipid interface. In the presence of cholesterol the fluorescence
1352 <    intensity ratio of D2O to H2O decreased, indicating a removal of
1353 <    water by the cholesterol, in keeping with the lifetime data. Altered
1354 <    hydration at the protein-lipid interface could affect conformation,
1355 <    thereby offering a new route by which membrane protein functioning
1356 <    may be modified.},
1331 >        hydrophobic interface ot the membrane proteins, gramicidin and apocytochrome
1332 >        C, using two related fluorescence spectroscopic approaches. The
1333 >        first approach exploited the fact that the presence of water in
1334 >        the excited state solvent cage of a fluorophore increases the rate
1335 >        of decay. For 1,6-diphenyl-1,3,5-hexatriene (DPH) and 1-palmitoyl-2-[[2-[4-(6-phenyl-trans-1,3,5-hexatrienyl)
1336 >        phenyl]ethyl]carbonyl]-3-sn-PC (DPH-PC), where the fluorophores
1337 >        are located in the hydrophobic core of the lipid bilayer, the introduction
1338 >        of gramicidin reduced the fluorescence lifetime, indicative of an
1339 >        increased presence of water in the bilayer. Since a high protein:lipid
1340 >        ratio was used, the fluorophores were forced to be adjacent to the
1341 >        protein hydrophobic surface, hence the presence of water in this
1342 >        region could be inferred. Cholesterol is known to reduce the water
1343 >        content of lipid bilayers and this effect was maintained at the
1344 >        protein-lipid interface with both gramicidin and apocytochrome C,
1345 >        again suggesting hydration in this region. The second approach was
1346 >        to use the fluorescence enhancement induced by exchanging deuterium
1347 >        oxide (D2O) for H2O. Both the fluorescence intensities of trimethylammonium-DPH,
1348 >        located in the lipid head group region, and of the gramicidin intrinsic
1349 >        tryptophans were greater in a D2O buffer compared with H2O, showing
1350 >        that the fluorophores were exposed to water in the bilayer at the
1351 >        protein-lipid interface. In the presence of cholesterol the fluorescence
1352 >        intensity ratio of D2O to H2O decreased, indicating a removal of
1353 >        water by the cholesterol, in keeping with the lifetime data. Altered
1354 >        hydration at the protein-lipid interface could affect conformation,
1355 >        thereby offering a new route by which membrane protein functioning
1356 >        may be modified.},
1357    annote = {Ju251 Times Cited:55 Cited References Count:44},
1358    issn = {0006-3495},
1359    uri = {<Go to ISI>://A1992JU25100002},
1360   }
1361  
1362 < @BOOK{Hockney1981,
1363 <  title = {Computer Simulation Using Particles},
1364 <  publisher = {McGraw-Hill},
1365 <  year = {1981},
1366 <  author = {R.W. Hockney and J.W. Eastwood},
1367 <  address = {New York},
1362 > @ARTICLE{Hoover1985,
1363 >  author = {W. G. Hoover},
1364 >  title = {Canonical Dynamics - Equilibrium Phase-Space Distributions},
1365 >  journal = {Physical Review A},
1366 >  year = {1985},
1367 >  volume = {31},
1368 >  pages = {1695-1697},
1369 >  number = {3},
1370 >  annote = {Acr30 Times Cited:1809 Cited References Count:11},
1371 >  issn = {1050-2947},
1372 >  uri = {<Go to ISI>://A1985ACR3000056},
1373   }
1374  
1375   @ARTICLE{Huh2004,
1376    author = {Y. Huh and N. M. Cann},
1377    title = {Discrimination in isotropic, nematic, and smectic phases of chiral
1378 <    calamitic molecules: A computer simulation study},
1378 >        calamitic molecules: A computer simulation study},
1379    journal = {Journal of Chemical Physics},
1380    year = {2004},
1381    volume = {121},
# Line 1391 | Line 1383 | Encoding: GBK
1383    number = {20},
1384    month = {Nov 22},
1385    abstract = {Racemic fluids of chiral calamitic molecules are investigated with
1386 <    molecular dynamics simulations. In particular, the phase behavior
1387 <    as a function of density is examined for eight racemates. The relationship
1388 <    between chiral discrimination and orientational order in the phase
1389 <    is explored. We find that the transition from the isotropic phase
1390 <    to a liquid crystal phase is accompanied by an increase in chiral
1391 <    discrimination, as measured by differences in radial distributions.
1392 <    Among ordered phases, discrimination is largest for smectic phases
1393 <    with a significant preference for heterochiral contact within the
1394 <    layers. (C) 2004 American Institute of Physics.},
1386 >        molecular dynamics simulations. In particular, the phase behavior
1387 >        as a function of density is examined for eight racemates. The relationship
1388 >        between chiral discrimination and orientational order in the phase
1389 >        is explored. We find that the transition from the isotropic phase
1390 >        to a liquid crystal phase is accompanied by an increase in chiral
1391 >        discrimination, as measured by differences in radial distributions.
1392 >        Among ordered phases, discrimination is largest for smectic phases
1393 >        with a significant preference for heterochiral contact within the
1394 >        layers. (C) 2004 American Institute of Physics.},
1395    annote = {870FJ Times Cited:0 Cited References Count:63},
1396    issn = {0021-9606},
1397    uri = {<Go to ISI>://000225042700059},
1398   }
1399  
1400 + @ARTICLE{Humphrey1996,
1401 +  author = {W. Humphrey and A. Dalke and K. Schulten},
1402 +  title = {VMD: Visual molecular dynamics},
1403 +  journal = {Journal of Molecular Graphics},
1404 +  year = {1996},
1405 +  volume = {14},
1406 +  pages = {33-\&},
1407 +  number = {1},
1408 +  month = {Feb},
1409 +  abstract = {VMD is a molecular graphics program designed for the display and analysis
1410 +        of molecular assemblies, in particular biopolymers such as proteins
1411 +        and nucleic acids. VMD can simultaneously display any number of
1412 +        structures using a wide variety of rendering styles and coloring
1413 +        methods. Molecules are displayed as one or more ''representations,''
1414 +        in which each representation embodies a particular rendering method
1415 +        and coloring scheme for a selected subset of atoms. The atoms displayed
1416 +        in each representation are chosen using an extensive atom selection
1417 +        syntax, which includes Boolean operators and regular expressions.
1418 +        VMD provides a complete graphical user interface for program control,
1419 +        as well as a text interface using the Tcl embeddable parser to allow
1420 +        for complex scripts with variable substitution, control loops, and
1421 +        function calls. Full session logging is supported, which produces
1422 +        a VMD command script for later playback. High-resolution raster
1423 +        images of displayed molecules may be produced by generating input
1424 +        scripts for use by a number of photorealistic image-rendering applications.
1425 +        VMD has also been expressly designed with the ability to animate
1426 +        molecular dynamics (MD) simulation trajectories, imported either
1427 +        from files or from a direct connection to a running MD simulation.
1428 +        VMD is the visualization component of MDScope, a set of tools for
1429 +        interactive problem solving in structural biology, which also includes
1430 +        the parallel MD program NAMD, and the MDCOMM software used to connect
1431 +        the visualization and simulation programs. VMD is written in C++,
1432 +        using an object-oriented design; the program, including source code
1433 +        and extensive documentation, is freely available via anonymous ftp
1434 +        and through the World Wide Web.},
1435 +  annote = {Uh515 Times Cited:1418 Cited References Count:19},
1436 +  issn = {0263-7855},
1437 +  uri = {<Go to ISI>://A1996UH51500005},
1438 + }
1439 +
1440   @ARTICLE{Izaguirre2001,
1441    author = {J. A. Izaguirre and D. P. Catarello and J. M. Wozniak and R. D. Skeel},
1442    title = {Langevin stabilization of molecular dynamics},
# Line 1415 | Line 1447 | Encoding: GBK
1447    number = {5},
1448    month = {Feb 1},
1449    abstract = {In this paper we show the possibility of using very mild stochastic
1450 <    damping to stabilize long time step integrators for Newtonian molecular
1451 <    dynamics. More specifically, stable and accurate integrations are
1452 <    obtained for damping coefficients that are only a few percent of
1453 <    the natural decay rate of processes of interest, such as the velocity
1454 <    autocorrelation function. Two new multiple time stepping integrators,
1455 <    Langevin Molly (LM) and Brunger-Brooks-Karplus-Molly (BBK-M), are
1456 <    introduced in this paper. Both use the mollified impulse method
1457 <    for the Newtonian term. LM uses a discretization of the Langevin
1458 <    equation that is exact for the constant force, and BBK-M uses the
1459 <    popular Brunger-Brooks-Karplus integrator (BBK). These integrators,
1460 <    along with an extrapolative method called LN, are evaluated across
1461 <    a wide range of damping coefficient values. When large damping coefficients
1462 <    are used, as one would for the implicit modeling of solvent molecules,
1463 <    the method LN is superior, with LM closely following. However, with
1464 <    mild damping of 0.2 ps(-1), LM produces the best results, allowing
1465 <    long time steps of 14 fs in simulations containing explicitly modeled
1466 <    flexible water. With BBK-M and the same damping coefficient, time
1467 <    steps of 12 fs are possible for the same system. Similar results
1468 <    are obtained for a solvated protein-DNA simulation of estrogen receptor
1469 <    ER with estrogen response element ERE. A parallel version of BBK-M
1470 <    runs nearly three times faster than the Verlet-I/r-RESPA (reversible
1471 <    reference system propagator algorithm) when using the largest stable
1472 <    time step on each one, and it also parallelizes well. The computation
1473 <    of diffusion coefficients for flexible water and ER/ERE shows that
1474 <    when mild damping of up to 0.2 ps-1 is used the dynamics are not
1475 <    significantly distorted. (C) 2001 American Institute of Physics.},
1450 >        damping to stabilize long time step integrators for Newtonian molecular
1451 >        dynamics. More specifically, stable and accurate integrations are
1452 >        obtained for damping coefficients that are only a few percent of
1453 >        the natural decay rate of processes of interest, such as the velocity
1454 >        autocorrelation function. Two new multiple time stepping integrators,
1455 >        Langevin Molly (LM) and Brunger-Brooks-Karplus-Molly (BBK-M), are
1456 >        introduced in this paper. Both use the mollified impulse method
1457 >        for the Newtonian term. LM uses a discretization of the Langevin
1458 >        equation that is exact for the constant force, and BBK-M uses the
1459 >        popular Brunger-Brooks-Karplus integrator (BBK). These integrators,
1460 >        along with an extrapolative method called LN, are evaluated across
1461 >        a wide range of damping coefficient values. When large damping coefficients
1462 >        are used, as one would for the implicit modeling of solvent molecules,
1463 >        the method LN is superior, with LM closely following. However, with
1464 >        mild damping of 0.2 ps(-1), LM produces the best results, allowing
1465 >        long time steps of 14 fs in simulations containing explicitly modeled
1466 >        flexible water. With BBK-M and the same damping coefficient, time
1467 >        steps of 12 fs are possible for the same system. Similar results
1468 >        are obtained for a solvated protein-DNA simulation of estrogen receptor
1469 >        ER with estrogen response element ERE. A parallel version of BBK-M
1470 >        runs nearly three times faster than the Verlet-I/r-RESPA (reversible
1471 >        reference system propagator algorithm) when using the largest stable
1472 >        time step on each one, and it also parallelizes well. The computation
1473 >        of diffusion coefficients for flexible water and ER/ERE shows that
1474 >        when mild damping of up to 0.2 ps-1 is used the dynamics are not
1475 >        significantly distorted. (C) 2001 American Institute of Physics.},
1476    annote = {397CQ Times Cited:14 Cited References Count:36},
1477    issn = {0021-9606},
1478    uri = {<Go to ISI>://000166676100020},
1479   }
1480  
1449 @ARTICLE{Torre1977,
1450  author = {Jose Garcia De La Torre, V.A. Bloomfield},
1451  title = {Hydrodynamic properties of macromolecular complexes. I. Translation},
1452  journal = {Biopolymers},
1453  year = {1977},
1454  volume = {16},
1455  pages = {1747-1763},
1456 }
1457
1481   @ARTICLE{Kane2000,
1482    author = {C. Kane and J. E. Marsden and M. Ortiz and M. West},
1483    title = {Variational integrators and the Newmark algorithm for conservative
1484 <    and dissipative mechanical systems},
1484 >        and dissipative mechanical systems},
1485    journal = {International Journal for Numerical Methods in Engineering},
1486    year = {2000},
1487    volume = {49},
# Line 1466 | Line 1489 | Encoding: GBK
1489    number = {10},
1490    month = {Dec 10},
1491    abstract = {The purpose of this work is twofold. First, we demonstrate analytically
1492 <    that the classical Newmark family as well as related integration
1493 <    algorithms are variational in the sense of the Veselov formulation
1494 <    of discrete mechanics. Such variational algorithms are well known
1495 <    to be symplectic and momentum preserving and to often have excellent
1496 <    global energy behaviour. This analytical result is verified through
1497 <    numerical examples and is believed to be one of the primary reasons
1498 <    that this class of algorithms performs so well. Second, we develop
1499 <    algorithms for mechanical systems with forcing, and in particular,
1500 <    for dissipative systems. In this case, we develop integrators that
1501 <    are based on a discretization of the Lagrange d'Alembert principle
1502 <    as well as on a variational formulation of dissipation. It is demonstrated
1503 <    that these types of structured integrators have good numerical behaviour
1504 <    in terms of obtaining the correct amounts by which the energy changes
1505 <    over the integration run. Copyright (C) 2000 John Wiley & Sons,
1506 <    Ltd.},
1492 >        that the classical Newmark family as well as related integration
1493 >        algorithms are variational in the sense of the Veselov formulation
1494 >        of discrete mechanics. Such variational algorithms are well known
1495 >        to be symplectic and momentum preserving and to often have excellent
1496 >        global energy behaviour. This analytical result is verified through
1497 >        numerical examples and is believed to be one of the primary reasons
1498 >        that this class of algorithms performs so well. Second, we develop
1499 >        algorithms for mechanical systems with forcing, and in particular,
1500 >        for dissipative systems. In this case, we develop integrators that
1501 >        are based on a discretization of the Lagrange d'Alembert principle
1502 >        as well as on a variational formulation of dissipation. It is demonstrated
1503 >        that these types of structured integrators have good numerical behaviour
1504 >        in terms of obtaining the correct amounts by which the energy changes
1505 >        over the integration run. Copyright (C) 2000 John Wiley & Sons,
1506 >        Ltd.},
1507    annote = {373CJ Times Cited:30 Cited References Count:41},
1508    issn = {0029-5981},
1509    uri = {<Go to ISI>://000165270600004},
# Line 1496 | Line 1519 | Encoding: GBK
1519    number = {2},
1520    month = {Jul 14},
1521    abstract = {The viscosity (eta) dependence of the folding rates for four sequences
1522 <    (the native state of three sequences is a beta sheet, while the
1523 <    fourth forms an alpha helix) is calculated for off-lattice models
1524 <    of proteins. Assuming that the dynamics is given by the Langevin
1525 <    equation, we show that the folding rates increase linearly at low
1526 <    viscosities eta, decrease as 1/eta at large eta, and have a maximum
1527 <    at intermediate values. The Kramers' theory of barrier crossing
1528 <    provides a quantitative fit of the numerical results. By mapping
1529 <    the simulation results to real proteins we estimate that for optimized
1530 <    sequences the time scale for forming a four turn alpha-helix topology
1531 <    is about 500 ns, whereas for beta sheet it is about 10 mu s.},
1522 >        (the native state of three sequences is a beta sheet, while the
1523 >        fourth forms an alpha helix) is calculated for off-lattice models
1524 >        of proteins. Assuming that the dynamics is given by the Langevin
1525 >        equation, we show that the folding rates increase linearly at low
1526 >        viscosities eta, decrease as 1/eta at large eta, and have a maximum
1527 >        at intermediate values. The Kramers' theory of barrier crossing
1528 >        provides a quantitative fit of the numerical results. By mapping
1529 >        the simulation results to real proteins we estimate that for optimized
1530 >        sequences the time scale for forming a four turn alpha-helix topology
1531 >        is about 500 ns, whereas for beta sheet it is about 10 mu s.},
1532    annote = {Xk293 Times Cited:77 Cited References Count:17},
1533    issn = {0031-9007},
1534    uri = {<Go to ISI>://A1997XK29300035},
# Line 1521 | Line 1544 | Encoding: GBK
1544    number = {7},
1545    month = {Aug 15},
1546    abstract = {Rigid-body molecular dynamics simulations typically are performed
1547 <    in a quaternion representation. The nonseparable form of the Hamiltonian
1548 <    in quaternions prevents the use of a standard leapfrog (Verlet)
1549 <    integrator, so nonsymplectic Runge-Kutta, multistep, or extrapolation
1550 <    methods are generally used, This is unfortunate since symplectic
1551 <    methods like Verlet exhibit superior energy conservation in long-time
1552 <    integrations. In this article, we describe an alternative method,
1553 <    which we call RSHAKE (for rotation-SHAKE), in which the entire rotation
1554 <    matrix is evolved (using the scheme of McLachlan and Scovel [J.
1555 <    Nonlin. Sci, 16 233 (1995)]) in tandem with the particle positions.
1556 <    We employ a fast approximate Newton solver to preserve the orthogonality
1557 <    of the rotation matrix. We test our method on a system of soft-sphere
1558 <    dipoles and compare with quaternion evolution using a 4th-order
1559 <    predictor-corrector integrator, Although the short-time error of
1560 <    the quaternion algorithm is smaller for fixed time step than that
1561 <    for RSHAKE, the quaternion scheme exhibits an energy drift which
1562 <    is not observed in simulations with RSHAKE, hence a fixed energy
1563 <    tolerance can be achieved by using a larger time step, The superiority
1564 <    of RSHAKE increases with system size. (C) 1997 American Institute
1565 <    of Physics.},
1547 >        in a quaternion representation. The nonseparable form of the Hamiltonian
1548 >        in quaternions prevents the use of a standard leapfrog (Verlet)
1549 >        integrator, so nonsymplectic Runge-Kutta, multistep, or extrapolation
1550 >        methods are generally used, This is unfortunate since symplectic
1551 >        methods like Verlet exhibit superior energy conservation in long-time
1552 >        integrations. In this article, we describe an alternative method,
1553 >        which we call RSHAKE (for rotation-SHAKE), in which the entire rotation
1554 >        matrix is evolved (using the scheme of McLachlan and Scovel [J.
1555 >        Nonlin. Sci, 16 233 (1995)]) in tandem with the particle positions.
1556 >        We employ a fast approximate Newton solver to preserve the orthogonality
1557 >        of the rotation matrix. We test our method on a system of soft-sphere
1558 >        dipoles and compare with quaternion evolution using a 4th-order
1559 >        predictor-corrector integrator, Although the short-time error of
1560 >        the quaternion algorithm is smaller for fixed time step than that
1561 >        for RSHAKE, the quaternion scheme exhibits an energy drift which
1562 >        is not observed in simulations with RSHAKE, hence a fixed energy
1563 >        tolerance can be achieved by using a larger time step, The superiority
1564 >        of RSHAKE increases with system size. (C) 1997 American Institute
1565 >        of Physics.},
1566    annote = {Xq332 Times Cited:11 Cited References Count:18},
1567    issn = {0021-9606},
1568    uri = {<Go to ISI>://A1997XQ33200046},
# Line 1548 | Line 1571 | Encoding: GBK
1571   @ARTICLE{Lansac2001,
1572    author = {Y. Lansac and M. A. Glaser and N. A. Clark},
1573    title = {Microscopic structure and dynamics of a partial bilayer smectic liquid
1574 <    crystal},
1574 >        crystal},
1575    journal = {Physical Review E},
1576    year = {2001},
1577    volume = {6405},
# Line 1556 | Line 1579 | Encoding: GBK
1579    number = {5},
1580    month = {Nov},
1581    abstract = {Cyanobiphenyls (nCB's) represent a useful and intensively studied
1582 <    class of mesogens. Many of the peculiar properties of nCB's (e.g.,
1583 <    the occurence of the partial bilayer smectic-A(d) phase) are thought
1584 <    to be a manifestation of short-range antiparallel association of
1585 <    neighboring molecules, resulting from strong dipole-dipole interactions
1586 <    between cyano groups. To test and extend existing models of microscopic
1587 <    ordering in nCB's, we carry out large-scale atomistic simulation
1588 <    studies of the microscopic structure and dynamics of the Sm-A(d)
1589 <    phase of 4-octyl-4'-cyanobiphenyl (8CB). We compute a variety of
1590 <    thermodynamic, structural, and dynamical properties for this material,
1591 <    and make a detailed comparison of our results with experimental
1592 <    measurements in order to validate our molecular model. Semiquantitative
1593 <    agreement with experiment is found: the smectic layer spacing and
1594 <    mass density are well reproduced, translational diffusion constants
1595 <    are similar to experiment, but the orientational ordering of alkyl
1596 <    chains is overestimated. This simulation provides a detailed picture
1597 <    of molecular conformation, smectic layer structure, and intermolecular
1598 <    correlations in Sm-A(d) 8CB, and demonstrates that pronounced short-range
1599 <    antiparallel association of molecules arising from dipole-dipole
1600 <    interactions plays a dominant role in determining the molecular-scale
1601 <    structure of 8CB.},
1582 >        class of mesogens. Many of the peculiar properties of nCB's (e.g.,
1583 >        the occurence of the partial bilayer smectic-A(d) phase) are thought
1584 >        to be a manifestation of short-range antiparallel association of
1585 >        neighboring molecules, resulting from strong dipole-dipole interactions
1586 >        between cyano groups. To test and extend existing models of microscopic
1587 >        ordering in nCB's, we carry out large-scale atomistic simulation
1588 >        studies of the microscopic structure and dynamics of the Sm-A(d)
1589 >        phase of 4-octyl-4'-cyanobiphenyl (8CB). We compute a variety of
1590 >        thermodynamic, structural, and dynamical properties for this material,
1591 >        and make a detailed comparison of our results with experimental
1592 >        measurements in order to validate our molecular model. Semiquantitative
1593 >        agreement with experiment is found: the smectic layer spacing and
1594 >        mass density are well reproduced, translational diffusion constants
1595 >        are similar to experiment, but the orientational ordering of alkyl
1596 >        chains is overestimated. This simulation provides a detailed picture
1597 >        of molecular conformation, smectic layer structure, and intermolecular
1598 >        correlations in Sm-A(d) 8CB, and demonstrates that pronounced short-range
1599 >        antiparallel association of molecules arising from dipole-dipole
1600 >        interactions plays a dominant role in determining the molecular-scale
1601 >        structure of 8CB.},
1602    annote = {Part 1 496QF Times Cited:10 Cited References Count:60},
1603    issn = {1063-651X},
1604    uri = {<Go to ISI>://000172406900063},
# Line 1591 | Line 1614 | Encoding: GBK
1614    number = {1},
1615    month = {Jan},
1616    abstract = {Recently, a new class of smectic liquid crystal phases characterized
1617 <    by the spontaneous formation of macroscopic chiral domains from
1618 <    achiral bent-core molecules has been discovered. We have carried
1619 <    out Monte Carlo simulations of a minimal hard spherocylinder dimer
1620 <    model to investigate the role of excluded volume interactions in
1621 <    determining the phase behavior of bent-core materials and to probe
1622 <    the molecular origins of polar and chiral symmetry breaking. We
1623 <    present the phase diagram of hard spherocylinder dimers of length-diameter
1624 <    ratio of 5 as a function of pressure or density and dimer opening
1625 <    angle psi. With decreasing psi, a transition from a nonpolar to
1626 <    a polar smectic A phase is observed near psi=167degrees, and the
1627 <    nematic phase becomes thermodynamically unstable for psi<135degrees.
1628 <    Free energy calculations indicate that the antipolar smectic A (SmAP(A))
1629 <    phase is more stable than the polar smectic A phase (SmAP(F)). No
1630 <    chiral smectic or biaxial nematic phases were found.},
1617 >        by the spontaneous formation of macroscopic chiral domains from
1618 >        achiral bent-core molecules has been discovered. We have carried
1619 >        out Monte Carlo simulations of a minimal hard spherocylinder dimer
1620 >        model to investigate the role of excluded volume interactions in
1621 >        determining the phase behavior of bent-core materials and to probe
1622 >        the molecular origins of polar and chiral symmetry breaking. We
1623 >        present the phase diagram of hard spherocylinder dimers of length-diameter
1624 >        ratio of 5 as a function of pressure or density and dimer opening
1625 >        angle psi. With decreasing psi, a transition from a nonpolar to
1626 >        a polar smectic A phase is observed near psi=167degrees, and the
1627 >        nematic phase becomes thermodynamically unstable for psi<135degrees.
1628 >        Free energy calculations indicate that the antipolar smectic A (SmAP(A))
1629 >        phase is more stable than the polar smectic A phase (SmAP(F)). No
1630 >        chiral smectic or biaxial nematic phases were found.},
1631    annote = {Part 1 646CM Times Cited:15 Cited References Count:38},
1632    issn = {1063-651X},
1633    uri = {<Go to ISI>://000181017300042},
1634   }
1635  
1613 @BOOK{Leach2001,
1614  title = {Molecular Modeling: Principles and Applications},
1615  publisher = {Pearson Educated Limited},
1616  year = {2001},
1617  author = {A. Leach},
1618  address = {Harlow, England},
1619  edition = {2nd},
1620 }
1621
1636   @ARTICLE{Leimkuhler1999,
1637    author = {B. Leimkuhler},
1638    title = {Reversible adaptive regularization: perturbed Kepler motion and classical
1639 <    atomic trajectories},
1639 >        atomic trajectories},
1640    journal = {Philosophical Transactions of the Royal Society of London Series
1641 <    a-Mathematical Physical and Engineering Sciences},
1641 >        a-Mathematical Physical and Engineering Sciences},
1642    year = {1999},
1643    volume = {357},
1644    pages = {1101-1133},
1645    number = {1754},
1646    month = {Apr 15},
1647    abstract = {Reversible and adaptive integration methods based on Kustaanheimo-Stiefel
1648 <    regularization and modified Sundman transformations are applied
1649 <    to simulate general perturbed Kepler motion and to compute classical
1650 <    trajectories of atomic systems (e.g. Rydberg atoms). The new family
1651 <    of reversible adaptive regularization methods also conserves angular
1652 <    momentum and exhibits superior energy conservation and numerical
1653 <    stability in long-time integrations. The schemes are appropriate
1654 <    for scattering, for astronomical calculations of escape time and
1655 <    long-term stability, and for classical and semiclassical studies
1656 <    of atomic dynamics. The components of an algorithm for trajectory
1657 <    calculations are described. Numerical experiments illustrate the
1658 <    effectiveness of the reversible approach.},
1648 >        regularization and modified Sundman transformations are applied
1649 >        to simulate general perturbed Kepler motion and to compute classical
1650 >        trajectories of atomic systems (e.g. Rydberg atoms). The new family
1651 >        of reversible adaptive regularization methods also conserves angular
1652 >        momentum and exhibits superior energy conservation and numerical
1653 >        stability in long-time integrations. The schemes are appropriate
1654 >        for scattering, for astronomical calculations of escape time and
1655 >        long-term stability, and for classical and semiclassical studies
1656 >        of atomic dynamics. The components of an algorithm for trajectory
1657 >        calculations are described. Numerical experiments illustrate the
1658 >        effectiveness of the reversible approach.},
1659    annote = {199EE Times Cited:11 Cited References Count:48},
1660    issn = {1364-503X},
1661    uri = {<Go to ISI>://000080466800007},
1662   }
1663  
1650 @BOOK{Leimkuhler2004,
1651  title = {Simulating Hamiltonian Dynamics},
1652  publisher = {Cambridge University Press},
1653  year = {2004},
1654  author = {B. Leimkuhler and S. Reich},
1655  address = {Cambridge},
1656 }
1657
1664   @ARTICLE{Levelut1981,
1665    author = {A. M. Levelut and R. J. Tarento and F. Hardouin and M. F. Achard
1666 <    and G. Sigaud},
1666 >        and G. Sigaud},
1667    title = {Number of Sa Phases},
1668    journal = {Physical Review A},
1669    year = {1981},
# Line 1672 | Line 1678 | Encoding: GBK
1678   @ARTICLE{Lieb1982,
1679    author = {W. R. Lieb and M. Kovalycsik and R. Mendelsohn},
1680    title = {Do Clinical-Levels of General-Anesthetics Affect Lipid Bilayers -
1681 <    Evidence from Raman-Scattering},
1681 >        Evidence from Raman-Scattering},
1682    journal = {Biochimica Et Biophysica Acta},
1683    year = {1982},
1684    volume = {688},
# Line 1685 | Line 1691 | Encoding: GBK
1691  
1692   @ARTICLE{Link1997,
1693    author = {D. R. Link and G. Natale and R. Shao and J. E. Maclennan and N. A.
1694 <    Clark and E. Korblova and D. M. Walba},
1694 >        Clark and E. Korblova and D. M. Walba},
1695    title = {Spontaneous formation of macroscopic chiral domains in a fluid smectic
1696 <    phase of achiral molecules},
1696 >        phase of achiral molecules},
1697    journal = {Science},
1698    year = {1997},
1699    volume = {278},
# Line 1695 | Line 1701 | Encoding: GBK
1701    number = {5345},
1702    month = {Dec 12},
1703    abstract = {A smectic liquid-crystal phase made from achiral molecules with bent
1704 <    cores was found to have fluid layers that exhibit two spontaneous
1705 <    symmetry-breaking instabilities: polar molecular orientational ordering
1706 <    about the layer normal and molecular tilt. These instabilities combine
1707 <    to form a chiral layer structure with a handedness that depends
1708 <    on the sign of the tilt. The bulk states are either antiferroelectric-racemic,
1709 <    with the layer polar direction and handedness alternating in sign
1710 <    from layer to layer, or antiferroelectric-chiral, which is of uniform
1711 <    layer handedness. Both states exhibit an electric field-induced
1712 <    transition from antiferroelectric to ferroelectric.},
1704 >        cores was found to have fluid layers that exhibit two spontaneous
1705 >        symmetry-breaking instabilities: polar molecular orientational ordering
1706 >        about the layer normal and molecular tilt. These instabilities combine
1707 >        to form a chiral layer structure with a handedness that depends
1708 >        on the sign of the tilt. The bulk states are either antiferroelectric-racemic,
1709 >        with the layer polar direction and handedness alternating in sign
1710 >        from layer to layer, or antiferroelectric-chiral, which is of uniform
1711 >        layer handedness. Both states exhibit an electric field-induced
1712 >        transition from antiferroelectric to ferroelectric.},
1713    annote = {Yl002 Times Cited:407 Cited References Count:25},
1714    issn = {0036-8075},
1715    uri = {<Go to ISI>://A1997YL00200028},
# Line 1712 | Line 1718 | Encoding: GBK
1718   @ARTICLE{Liwo2005,
1719    author = {A. Liwo and M. Khalili and H. A. Scheraga},
1720    title = {Ab initio simulations of protein folding pathways by molecular dynamics
1721 <    with the united-residue (UNRES) model of polypeptide chains},
1721 >        with the united-residue (UNRES) model of polypeptide chains},
1722    journal = {Febs Journal},
1723    year = {2005},
1724    volume = {272},
# Line 1726 | Line 1732 | Encoding: GBK
1732   @ARTICLE{Luty1994,
1733    author = {B. A. Luty and M. E. Davis and I. G. Tironi and W. F. Vangunsteren},
1734    title = {A Comparison of Particle-Particle, Particle-Mesh and Ewald Methods
1735 <    for Calculating Electrostatic Interactions in Periodic Molecular-Systems},
1735 >        for Calculating Electrostatic Interactions in Periodic Molecular-Systems},
1736    journal = {Molecular Simulation},
1737    year = {1994},
1738    volume = {14},
1739    pages = {11-20},
1740    number = {1},
1741    abstract = {We compare the Particle-Particle Particle-Mesh (PPPM) and Ewald methods
1742 <    for calculating electrostatic interactions in periodic molecular
1743 <    systems. A brief comparison of the theories shows that the methods
1744 <    are very similar differing mainly in the technique which is used
1745 <    to perform the ''k-space'' or mesh calculation. Because the PPPM
1746 <    utilizes the highly efficient numerical Fast Fourier Transform (FFT)
1747 <    method it requires significantly less computational effort than
1748 <    the Ewald method and scale's almost linearly with system size.},
1742 >        for calculating electrostatic interactions in periodic molecular
1743 >        systems. A brief comparison of the theories shows that the methods
1744 >        are very similar differing mainly in the technique which is used
1745 >        to perform the ''k-space'' or mesh calculation. Because the PPPM
1746 >        utilizes the highly efficient numerical Fast Fourier Transform (FFT)
1747 >        method it requires significantly less computational effort than
1748 >        the Ewald method and scale's almost linearly with system size.},
1749    annote = {Qf464 Times Cited:50 Cited References Count:20},
1750    issn = {0892-7022},
1751    uri = {<Go to ISI>://A1994QF46400002},
1752   }
1753  
1754 < @BOOK{Marion1990,
1755 <  title = {Classical Dynamics of Particles and Systems},
1756 <  publisher = {Academic Press},
1757 <  year = {1990},
1758 <  author = {J.~B. Marion},
1759 <  address = {New York},
1760 <  edition = {2rd},
1754 > @ARTICLE{Marrink1994,
1755 >  author = {S. J. Marrink and H. J. C. Berendsen},
1756 >  title = {Simulation of Water Transport through a Lipid-Membrane},
1757 >  journal = {Journal of Physical Chemistry},
1758 >  year = {1994},
1759 >  volume = {98},
1760 >  pages = {4155-4168},
1761 >  number = {15},
1762 >  month = {Apr 14},
1763 >  abstract = {To obtain insight in the process of water permeation through a lipid
1764 >        membrane, we performed molecular dynamics simulations on a phospholipid
1765 >        (DPPC)/water system with atomic detail. Since the actual process
1766 >        of permeation is too slow to be studied directly, we deduced the
1767 >        permeation rate indirectly via computation of the free energy and
1768 >        diffusion rate profiles of a water molecule across the bilayer.
1769 >        We conclude that the permeation of water through a lipid membrane
1770 >        cannot be described adequately by a simple homogeneous solubility-diffusion
1771 >        model. Both the excess free energy and the diffusion rate strongly
1772 >        depend on the position in the membrane, as a result from the inhomogeneous
1773 >        nature of the membrane. The calculated excess free energy profile
1774 >        has a shallow slope and a maximum height of 26 kJ/mol. The diffusion
1775 >        rate is highest in the middle of the membrane where the lipid density
1776 >        is low. In the interfacial region almost all water molecules are
1777 >        bound by the lipid headgroups, and the diffusion turns out to be
1778 >        1 order of magnitude smaller. The total transport process is essentially
1779 >        determined by the free energy barrier. The rate-limiting step is
1780 >        the permeation through the dense part of the lipid tails, where
1781 >        the resistance is highest. We found a permeation rate of 7(+/-3)
1782 >        x 10(-2) cm/s at 350 K, comparable to experimental values for DPPC
1783 >        membranes, if corrected for the temperature of the simulation. Taking
1784 >        the inhomogeneity of the membrane into account, we define a new
1785 >        ''four-region'' model which seems to be more realistic than the
1786 >        ''two-phase'' solubility-diffusion model.},
1787 >  annote = {Ng219 Times Cited:187 Cited References Count:25},
1788 >  issn = {0022-3654},
1789 >  uri = {<Go to ISI>://A1994NG21900040},
1790   }
1791  
1792   @ARTICLE{Marsden1998,
1793    author = {J. E. Marsden and G. W. Patrick and S. Shkoller},
1794    title = {Multisymplectic geometry, variational integrators, and nonlinear
1795 <    PDEs},
1795 >        PDEs},
1796    journal = {Communications in Mathematical Physics},
1797    year = {1998},
1798    volume = {199},
# Line 1765 | Line 1800 | Encoding: GBK
1800    number = {2},
1801    month = {Dec},
1802    abstract = {This paper presents a geometric-variational approach to continuous
1803 <    and discrete mechanics and field theories. Using multisymplectic
1804 <    geometry, we show that the existence of the fundamental geometric
1805 <    structures as well as their preservation along solutions can be
1806 <    obtained directly from the variational principle. In particular,
1807 <    we prove that a unique multisymplectic structure is obtained by
1808 <    taking the derivative of an action function, and use this structure
1809 <    to prove covariant generalizations of conservation of symplecticity
1810 <    and Noether's theorem. Natural discretization schemes for PDEs,
1811 <    which have these important preservation properties, then follow
1812 <    by choosing a discrete action functional. In the case of mechanics,
1813 <    we recover the variational symplectic integrators of Veselov type,
1814 <    while for PDEs we obtain covariant spacetime integrators which conserve
1815 <    the corresponding discrete multisymplectic form as well as the discrete
1816 <    momentum mappings corresponding to symmetries. We show that the
1817 <    usual notion of symplecticity along an infinite-dimensional space
1818 <    of fields can be naturally obtained by making a spacetime split.
1819 <    All of the aspects of our method are demonstrated with a nonlinear
1820 <    sine-Gordon equation, including computational results and a comparison
1821 <    with other discretization schemes.},
1803 >        and discrete mechanics and field theories. Using multisymplectic
1804 >        geometry, we show that the existence of the fundamental geometric
1805 >        structures as well as their preservation along solutions can be
1806 >        obtained directly from the variational principle. In particular,
1807 >        we prove that a unique multisymplectic structure is obtained by
1808 >        taking the derivative of an action function, and use this structure
1809 >        to prove covariant generalizations of conservation of symplecticity
1810 >        and Noether's theorem. Natural discretization schemes for PDEs,
1811 >        which have these important preservation properties, then follow
1812 >        by choosing a discrete action functional. In the case of mechanics,
1813 >        we recover the variational symplectic integrators of Veselov type,
1814 >        while for PDEs we obtain covariant spacetime integrators which conserve
1815 >        the corresponding discrete multisymplectic form as well as the discrete
1816 >        momentum mappings corresponding to symmetries. We show that the
1817 >        usual notion of symplecticity along an infinite-dimensional space
1818 >        of fields can be naturally obtained by making a spacetime split.
1819 >        All of the aspects of our method are demonstrated with a nonlinear
1820 >        sine-Gordon equation, including computational results and a comparison
1821 >        with other discretization schemes.},
1822    annote = {154RH Times Cited:88 Cited References Count:36},
1823    issn = {0010-3616},
1824    uri = {<Go to ISI>://000077902200006},
1825   }
1826  
1827 < @ARTICLE{McLachlan1993,
1828 <  author = {R.~I McLachlan},
1829 <  title = {Explicit Lie-Poisson integration and the Euler equations},
1830 <  journal = {prl},
1831 <  year = {1993},
1832 <  volume = {71},
1833 <  pages = {3043-3046},
1827 > @ARTICLE{Matthey2004,
1828 >  author = {T. Matthey and T. Cickovski and S. Hampton and A. Ko and Q. Ma and
1829 >        M. Nyerges and T. Raeder and T. Slabach and J. A. Izaguirre},
1830 >  title = {ProtoMol, an object-oriented framework for prototyping novel algorithms
1831 >        for molecular dynamics},
1832 >  journal = {Acm Transactions on Mathematical Software},
1833 >  year = {2004},
1834 >  volume = {30},
1835 >  pages = {237-265},
1836 >  number = {3},
1837 >  month = {Sep},
1838 >  abstract = {PROTOMOL is a high-performance framework in C++ for rapid prototyping
1839 >        of novel algorithms for molecular dynamics and related applications.
1840 >        Its flexibility is achieved primarily through the use of inheritance
1841 >        and design patterns (object-oriented programming): Performance is
1842 >        obtained by using templates that enable generation of efficient
1843 >        code for sections critical to performance (generic programming).
1844 >        The framework encapsulates important optimizations that can be used
1845 >        by developers, such as parallelism in the force computation. Its
1846 >        design is based on domain analysis of numerical integrators for
1847 >        molecular dynamics (MD) and of fast solvers for the force computation,
1848 >        particularly due to electrostatic interactions. Several new and
1849 >        efficient algorithms are implemented in PROTOMOL. Finally, it is
1850 >        shown that PROTOMOL'S sequential performance is excellent when compared
1851 >        to a leading MD program, and that it scales well for moderate number
1852 >        of processors. Binaries and source codes for Windows, Linux, Solaris,
1853 >        IRIX, HP-UX, and AIX platforms are available under open source license
1854 >        at http://protomol.sourceforge.net.},
1855 >  annote = {860EP Times Cited:2 Cited References Count:52},
1856 >  issn = {0098-3500},
1857 >  uri = {<Go to ISI>://000224325600001},
1858   }
1859  
1860 < @ARTICLE{McLachlan1998a,
1860 > @ARTICLE{McLachlan1998,
1861    author = {R. I. McLachlan and G. R. W. Quispel},
1862    title = {Generating functions for dynamical systems with symmetries, integrals,
1863 <    and differential invariants},
1863 >        and differential invariants},
1864    journal = {Physica D},
1865    year = {1998},
1866    volume = {112},
# Line 1809 | Line 1868 | Encoding: GBK
1868    number = {1-2},
1869    month = {Jan 15},
1870    abstract = {We give a survey and some new examples of generating functions for
1871 <    systems with symplectic structure, systems with a first integral,
1872 <    systems that preserve volume, and systems with symmetries and/or
1873 <    time-reversing symmetries. Both ODEs and maps are treated, and we
1874 <    discuss how generating functions may be used in the structure-preserving
1875 <    numerical integration of ODEs with the above properties.},
1871 >        systems with symplectic structure, systems with a first integral,
1872 >        systems that preserve volume, and systems with symmetries and/or
1873 >        time-reversing symmetries. Both ODEs and maps are treated, and we
1874 >        discuss how generating functions may be used in the structure-preserving
1875 >        numerical integration of ODEs with the above properties.},
1876    annote = {Yt049 Times Cited:7 Cited References Count:26},
1877    issn = {0167-2789},
1878    uri = {<Go to ISI>://000071558900021},
1879   }
1880  
1881 < @ARTICLE{McLachlan1998,
1881 > @ARTICLE{McLachlan1998a,
1882    author = {R. I. McLachlan and G. R. W. Quispel and G. S. Turner},
1883    title = {Numerical integrators that preserve symmetries and reversing symmetries},
1884    journal = {Siam Journal on Numerical Analysis},
# Line 1829 | Line 1888 | Encoding: GBK
1888    number = {2},
1889    month = {Apr},
1890    abstract = {We consider properties of flows, the relationships between them, and
1891 <    whether numerical integrators can be made to preserve these properties.
1892 <    This is done in the context of automorphisms and antiautomorphisms
1893 <    of a certain group generated by maps associated to vector fields.
1894 <    This new framework unifies several known constructions. We also
1895 <    use the concept of #covariance# of a numerical method with respect
1896 <    to a group of coordinate transformations. The main application is
1897 <    to explore the relationship between spatial symmetries, reversing
1898 <    symmetries, and time symmetry of flows and numerical integrators.},
1891 >        whether numerical integrators can be made to preserve these properties.
1892 >        This is done in the context of automorphisms and antiautomorphisms
1893 >        of a certain group generated by maps associated to vector fields.
1894 >        This new framework unifies several known constructions. We also
1895 >        use the concept of #covariance# of a numerical method with respect
1896 >        to a group of coordinate transformations. The main application is
1897 >        to explore the relationship between spatial symmetries, reversing
1898 >        symmetries, and time symmetry of flows and numerical integrators.},
1899    annote = {Zc449 Times Cited:14 Cited References Count:33},
1900    issn = {0036-1429},
1901    uri = {<Go to ISI>://000072580500010},
# Line 1852 | Line 1911 | Encoding: GBK
1911    number = {1},
1912    month = {Feb},
1913    abstract = {In this paper we revisit the Moser-Veselov description of the free
1914 <    rigid body in body coordinates, which, in the 3 x 3 case, can be
1915 <    implemented as an explicit, second-order, integrable approximation
1916 <    of the continuous solution. By backward error analysis, we study
1917 <    the modified vector field which is integrated exactly by the discrete
1918 <    algorithm. We deduce that the discrete Moser-Veselov (DMV) is well
1919 <    approximated to higher order by time reparametrizations of the continuous
1920 <    equations (modified vector field). We use the modified vector field
1921 <    to scale the initial data of the DMV to improve the order of the
1922 <    approximation and show the equivalence of the DMV and the RATTLE
1923 <    algorithm. Numerical integration with these preprocessed initial
1924 <    data is several orders of magnitude more accurate than the original
1925 <    DMV and RATTLE approach.},
1914 >        rigid body in body coordinates, which, in the 3 x 3 case, can be
1915 >        implemented as an explicit, second-order, integrable approximation
1916 >        of the continuous solution. By backward error analysis, we study
1917 >        the modified vector field which is integrated exactly by the discrete
1918 >        algorithm. We deduce that the discrete Moser-Veselov (DMV) is well
1919 >        approximated to higher order by time reparametrizations of the continuous
1920 >        equations (modified vector field). We use the modified vector field
1921 >        to scale the initial data of the DMV to improve the order of the
1922 >        approximation and show the equivalence of the DMV and the RATTLE
1923 >        algorithm. Numerical integration with these preprocessed initial
1924 >        data is several orders of magnitude more accurate than the original
1925 >        DMV and RATTLE approach.},
1926    annote = {911NS Times Cited:0 Cited References Count:14},
1927    issn = {1615-3375},
1928    uri = {<Go to ISI>://000228011900003},
1929   }
1930  
1931 + @ARTICLE{Meineke2005,
1932 +  author = {M. A. Meineke and C. F. Vardeman and T. Lin and C. J. Fennell and
1933 +        J. D. Gezelter},
1934 +  title = {OOPSE: An object-oriented parallel simulation engine for molecular
1935 +        dynamics},
1936 +  journal = {Journal of Computational Chemistry},
1937 +  year = {2005},
1938 +  volume = {26},
1939 +  pages = {252-271},
1940 +  number = {3},
1941 +  month = {Feb},
1942 +  abstract = {OOPSE is a new molecular dynamics simulation program that is capable
1943 +        of efficiently integrating equations of motion for atom types with
1944 +        orientational degrees of freedom (e.g. #sticky# atoms and point
1945 +        dipoles). Transition metals can also be simulated using the embedded
1946 +        atom method (EAM) potential included in the code. Parallel simulations
1947 +        are carried out using the force-based decomposition method. Simulations
1948 +        are specified using a very simple C-based meta-data language. A
1949 +        number of advanced integrators are included, and the basic integrator
1950 +        for orientational dynamics provides substantial improvements over
1951 +        older quaternion-based schemes. (C) 2004 Wiley Periodicals, Inc.},
1952 +  annote = {891CF Times Cited:1 Cited References Count:56},
1953 +  issn = {0192-8651},
1954 +  uri = {<Go to ISI>://000226558200006},
1955 + }
1956 +
1957 + @ARTICLE{Melchionna1993,
1958 +  author = {S. Melchionna and G. Ciccotti and B. L. Holian},
1959 +  title = {Hoover Npt Dynamics for Systems Varying in Shape and Size},
1960 +  journal = {Molecular Physics},
1961 +  year = {1993},
1962 +  volume = {78},
1963 +  pages = {533-544},
1964 +  number = {3},
1965 +  month = {Feb 20},
1966 +  abstract = {In this paper we write down equations of motion (following the approach
1967 +        pioneered by Hoover) for an exact isothermal-isobaric molecular
1968 +        dynamics simulation, and we extend them to multiple thermostating
1969 +        rates, to a shape-varying cell and to molecular systems, coherently
1970 +        with the previous 'extended system method'. An integration scheme
1971 +        is proposed together with a numerical illustration of the method.},
1972 +  annote = {Kq355 Times Cited:172 Cited References Count:17},
1973 +  issn = {0026-8976},
1974 +  uri = {<Go to ISI>://A1993KQ35500002},
1975 + }
1976 +
1977   @ARTICLE{Memmer2002,
1978    author = {R. Memmer},
1979    title = {Liquid crystal phases of achiral banana-shaped molecules: a computer
1980 <    simulation study},
1980 >        simulation study},
1981    journal = {Liquid Crystals},
1982    year = {2002},
1983    volume = {29},
# Line 1880 | Line 1985 | Encoding: GBK
1985    number = {4},
1986    month = {Apr},
1987    abstract = {The phase behaviour of achiral banana-shaped molecules was studied
1988 <    by computer simulation. The banana-shaped molecules were described
1989 <    by model intermolecular interactions based on the Gay-Berne potential.
1990 <    The characteristic molecular structure was considered by joining
1991 <    two calamitic Gay-Berne particles through a bond to form a biaxial
1992 <    molecule of point symmetry group C-2v with a suitable bending angle.
1993 <    The dependence on temperature of systems of N=1024 rigid banana-shaped
1994 <    molecules with bending angle phi=140degrees has been studied by
1995 <    means of Monte Carlo simulations in the isobaric-isothermal ensemble
1996 <    (NpT). On cooling an isotropic system, two phase transitions characterized
1997 <    by phase transition enthalpy, entropy and relative volume change
1998 <    have been observed. For the first time by computer simulation of
1999 <    a many-particle system of banana-shaped molecules, at low temperature
2000 <    an untilted smectic phase showing a global phase biaxiality and
2001 <    a spontaneous local polarization in the layers, i.e. a local polar
2002 <    arrangement of the steric dipoles, with an antiferroelectric-like
2003 <    superstructure could be proven, a phase structure which recently
2004 <    has been discovered experimentally. Additionally, at intermediate
2005 <    temperature a nematic-like phase has been proved, whereas close
2006 <    to the transition to the smectic phase hints of a spontaneous achiral
2007 <    symmetry breaking have been determined. Here, in the absence of
2008 <    a layered structure a helical superstructure has been formed. All
2009 <    phases have been characterized by visual representations of selected
2010 <    configurations, scalar and pseudoscalar correlation functions, and
2011 <    order parameters.},
1988 >        by computer simulation. The banana-shaped molecules were described
1989 >        by model intermolecular interactions based on the Gay-Berne potential.
1990 >        The characteristic molecular structure was considered by joining
1991 >        two calamitic Gay-Berne particles through a bond to form a biaxial
1992 >        molecule of point symmetry group C-2v with a suitable bending angle.
1993 >        The dependence on temperature of systems of N=1024 rigid banana-shaped
1994 >        molecules with bending angle phi=140degrees has been studied by
1995 >        means of Monte Carlo simulations in the isobaric-isothermal ensemble
1996 >        (NpT). On cooling an isotropic system, two phase transitions characterized
1997 >        by phase transition enthalpy, entropy and relative volume change
1998 >        have been observed. For the first time by computer simulation of
1999 >        a many-particle system of banana-shaped molecules, at low temperature
2000 >        an untilted smectic phase showing a global phase biaxiality and
2001 >        a spontaneous local polarization in the layers, i.e. a local polar
2002 >        arrangement of the steric dipoles, with an antiferroelectric-like
2003 >        superstructure could be proven, a phase structure which recently
2004 >        has been discovered experimentally. Additionally, at intermediate
2005 >        temperature a nematic-like phase has been proved, whereas close
2006 >        to the transition to the smectic phase hints of a spontaneous achiral
2007 >        symmetry breaking have been determined. Here, in the absence of
2008 >        a layered structure a helical superstructure has been formed. All
2009 >        phases have been characterized by visual representations of selected
2010 >        configurations, scalar and pseudoscalar correlation functions, and
2011 >        order parameters.},
2012    annote = {531HT Times Cited:12 Cited References Count:37},
2013    issn = {0267-8292},
2014    uri = {<Go to ISI>://000174410500001},
2015   }
2016  
1912 @ARTICLE{Metropolis1949,
1913  author = {N. Metropolis and S. Ulam},
1914  title = {The $\mbox{Monte Carlo}$ Method},
1915  journal = {J. Am. Stat. Ass.},
1916  year = {1949},
1917  volume = {44},
1918  pages = {335-341},
1919 }
1920
2017   @ARTICLE{Mielke2004,
2018    author = {S. P. Mielke and W. H. Fink and V. V. Krishnan and N. Gronbech-Jensen
2019 <    and C. J. Benham},
2019 >        and C. J. Benham},
2020    title = {Transcription-driven twin supercoiling of a DNA loop: A Brownian
2021 <    dynamics study},
2021 >        dynamics study},
2022    journal = {Journal of Chemical Physics},
2023    year = {2004},
2024    volume = {121},
# Line 1930 | Line 2026 | Encoding: GBK
2026    number = {16},
2027    month = {Oct 22},
2028    abstract = {The torque generated by RNA polymerase as it tracks along double-stranded
2029 <    DNA can potentially induce long-range structural deformations integral
2030 <    to mechanisms of biological significance in both prokaryotes and
2031 <    eukaryotes. In this paper, we introduce a dynamic computer model
2032 <    for investigating this phenomenon. Duplex DNA is represented as
2033 <    a chain of hydrodynamic beads interacting through potentials of
2034 <    linearly elastic stretching, bending, and twisting, as well as excluded
2035 <    volume. The chain, linear when relaxed, is looped to form two open
2036 <    but topologically constrained subdomains. This permits the dynamic
2037 <    introduction of torsional stress via a centrally applied torque.
2038 <    We simulate by Brownian dynamics the 100 mus response of a 477-base
2039 <    pair B-DNA template to the localized torque generated by the prokaryotic
2040 <    transcription ensemble. Following a sharp rise at early times, the
2041 <    distributed twist assumes a nearly constant value in both subdomains,
2042 <    and a succession of supercoiling deformations occurs as superhelical
2043 <    stress is increasingly partitioned to writhe. The magnitude of writhe
2044 <    surpasses that of twist before also leveling off when the structure
2045 <    reaches mechanical equilibrium with the torsional load. Superhelicity
2046 <    is simultaneously right handed in one subdomain and left handed
2047 <    in the other, as predicted by the #transcription-induced##twin-supercoiled-domain#
2048 <    model [L. F. Liu and J. C. Wang, Proc. Natl. Acad. Sci. U.S.A. 84,
2049 <    7024 (1987)]. The properties of the chain at the onset of writhing
2050 <    agree well with predictions from theory, and the generated stress
2051 <    is ample for driving secondary structural transitions in physiological
2052 <    DNA. (C) 2004 American Institute of Physics.},
2029 >        DNA can potentially induce long-range structural deformations integral
2030 >        to mechanisms of biological significance in both prokaryotes and
2031 >        eukaryotes. In this paper, we introduce a dynamic computer model
2032 >        for investigating this phenomenon. Duplex DNA is represented as
2033 >        a chain of hydrodynamic beads interacting through potentials of
2034 >        linearly elastic stretching, bending, and twisting, as well as excluded
2035 >        volume. The chain, linear when relaxed, is looped to form two open
2036 >        but topologically constrained subdomains. This permits the dynamic
2037 >        introduction of torsional stress via a centrally applied torque.
2038 >        We simulate by Brownian dynamics the 100 mus response of a 477-base
2039 >        pair B-DNA template to the localized torque generated by the prokaryotic
2040 >        transcription ensemble. Following a sharp rise at early times, the
2041 >        distributed twist assumes a nearly constant value in both subdomains,
2042 >        and a succession of supercoiling deformations occurs as superhelical
2043 >        stress is increasingly partitioned to writhe. The magnitude of writhe
2044 >        surpasses that of twist before also leveling off when the structure
2045 >        reaches mechanical equilibrium with the torsional load. Superhelicity
2046 >        is simultaneously right handed in one subdomain and left handed
2047 >        in the other, as predicted by the #transcription-induced##twin-supercoiled-domain#
2048 >        model [L. F. Liu and J. C. Wang, Proc. Natl. Acad. Sci. U.S.A. 84,
2049 >        7024 (1987)]. The properties of the chain at the onset of writhing
2050 >        agree well with predictions from theory, and the generated stress
2051 >        is ample for driving secondary structural transitions in physiological
2052 >        DNA. (C) 2004 American Institute of Physics.},
2053    annote = {861ZF Times Cited:3 Cited References Count:34},
2054    issn = {0021-9606},
2055    uri = {<Go to ISI>://000224456500064},
# Line 1962 | Line 2058 | Encoding: GBK
2058   @ARTICLE{Naess2001,
2059    author = {S. N. Naess and H. M. Adland and A. Mikkelsen and A. Elgsaeter},
2060    title = {Brownian dynamics simulation of rigid bodies and segmented polymer
2061 <    chains. Use of Cartesian rotation vectors as the generalized coordinates
2062 <    describing angular orientations},
2063 <  journal = {Physica A},
2064 <  year = {2001},
2065 <  volume = {294},
2066 <  pages = {323-339},
2067 <  number = {3-4},
2068 <  month = {May 15},
2069 <  abstract = {The three Eulerian angles constitute the classical choice of generalized
2070 <    coordinates used to describe the three degrees of rotational freedom
2071 <    of a rigid body, but it has long been known that this choice yields
2072 <    singular equations of motion. The latter is also true when Eulerian
2073 <    angles are used in Brownian dynamics analyses of the angular orientation
2074 <    of single rigid bodies and segmented polymer chains. Starting from
2075 <    kinetic theory we here show that by instead employing the three
2076 <    components of Cartesian rotation vectors as the generalized coordinates
2077 <    describing angular orientation, no singularity appears in the configuration
2078 <    space diffusion equation and the associated Brownian dynamics algorithm.
2079 <    The suitability of Cartesian rotation vectors in Brownian dynamics
2080 <    simulations of segmented polymer chains with spring-like or ball-socket
2081 <    joints is discussed. (C) 2001 Elsevier Science B.V. All rights reserved.},
2061 >        chains. Use of Cartesian rotation vectors as the generalized coordinates
2062 >        describing angular orientations},
2063 >  journal = {Physica A},
2064 >  year = {2001},
2065 >  volume = {294},
2066 >  pages = {323-339},
2067 >  number = {3-4},
2068 >  month = {May 15},
2069 >  abstract = {The three Eulerian angles constitute the classical choice of generalized
2070 >        coordinates used to describe the three degrees of rotational freedom
2071 >        of a rigid body, but it has long been known that this choice yields
2072 >        singular equations of motion. The latter is also true when Eulerian
2073 >        angles are used in Brownian dynamics analyses of the angular orientation
2074 >        of single rigid bodies and segmented polymer chains. Starting from
2075 >        kinetic theory we here show that by instead employing the three
2076 >        components of Cartesian rotation vectors as the generalized coordinates
2077 >        describing angular orientation, no singularity appears in the configuration
2078 >        space diffusion equation and the associated Brownian dynamics algorithm.
2079 >        The suitability of Cartesian rotation vectors in Brownian dynamics
2080 >        simulations of segmented polymer chains with spring-like or ball-socket
2081 >        joints is discussed. (C) 2001 Elsevier Science B.V. All rights reserved.},
2082    annote = {433TA Times Cited:7 Cited References Count:19},
2083    issn = {0378-4371},
2084    uri = {<Go to ISI>://000168774800005},
# Line 1991 | Line 2087 | Encoding: GBK
2087   @ARTICLE{Niori1996,
2088    author = {T. Niori and T. Sekine and J. Watanabe and T. Furukawa and H. Takezoe},
2089    title = {Distinct ferroelectric smectic liquid crystals consisting of banana
2090 <    shaped achiral molecules},
2090 >        shaped achiral molecules},
2091    journal = {Journal of Materials Chemistry},
2092    year = {1996},
2093    volume = {6},
# Line 1999 | Line 2095 | Encoding: GBK
2095    number = {7},
2096    month = {Jul},
2097    abstract = {The synthesis of a banana-shaped molecule is reported and it is found
2098 <    that the smectic phase which it forms is biaxial with the molecules
2099 <    packed in the best,direction into a layer. Because of this characteristic
2100 <    packing, spontaneous polarization appears parallel to the layer
2101 <    and switches on reversal of an applied electric field. This is the
2102 <    first obvious example of ferroelectricity in an achiral smectic
2103 <    phase and is ascribed to the C-2v symmetry of the molecular packing.},
2098 >        that the smectic phase which it forms is biaxial with the molecules
2099 >        packed in the best,direction into a layer. Because of this characteristic
2100 >        packing, spontaneous polarization appears parallel to the layer
2101 >        and switches on reversal of an applied electric field. This is the
2102 >        first obvious example of ferroelectricity in an achiral smectic
2103 >        phase and is ascribed to the C-2v symmetry of the molecular packing.},
2104    annote = {Ux855 Times Cited:447 Cited References Count:18},
2105    issn = {0959-9428},
2106    uri = {<Go to ISI>://A1996UX85500025},
# Line 2020 | Line 2116 | Encoding: GBK
2116    number = {5},
2117    month = {may},
2118    abstract = {We Studied the structural changes of bilayer vesicles induced by mechanical
2119 <    forces using a Brownian dynamics simulation. Two nanoparticles,
2120 <    which interact repulsively with amphiphilic molecules, are put inside
2121 <    a vesicle. The position of one nanoparticle is fixed, and the other
2122 <    is moved by a constant force as in optical-trapping experiments.
2123 <    First, the pulled vesicle stretches into a pear or tube shape. Then
2124 <    the inner monolayer in the tube-shaped region is deformed, and a
2125 <    cylindrical structure is formed between two vesicles. After stretching
2126 <    the cylindrical region, fission occurs near the moved vesicle. Soon
2127 <    after this the cylindrical region shrinks. The trapping force similar
2128 <    to 100 pN is needed to induce the formation of the cylindrical structure
2129 <    and fission.},
2119 >        forces using a Brownian dynamics simulation. Two nanoparticles,
2120 >        which interact repulsively with amphiphilic molecules, are put inside
2121 >        a vesicle. The position of one nanoparticle is fixed, and the other
2122 >        is moved by a constant force as in optical-trapping experiments.
2123 >        First, the pulled vesicle stretches into a pear or tube shape. Then
2124 >        the inner monolayer in the tube-shaped region is deformed, and a
2125 >        cylindrical structure is formed between two vesicles. After stretching
2126 >        the cylindrical region, fission occurs near the moved vesicle. Soon
2127 >        after this the cylindrical region shrinks. The trapping force similar
2128 >        to 100 pN is needed to induce the formation of the cylindrical structure
2129 >        and fission.},
2130    annote = {Part 1 568PX Times Cited:5 Cited References Count:39},
2131    issn = {1063-651X},
2132    uri = {<Go to ISI>://000176552300084},
# Line 2046 | Line 2142 | Encoding: GBK
2142    number = {20},
2143    month = {Nov 22},
2144    abstract = {We studied the fusion dynamics of vesicles using a Brownian dynamics
2145 <    simulation. Amphiphilic molecules spontaneously form vesicles with
2146 <    a bilayer structure. Two vesicles come into contact and form a stalk
2147 <    intermediate, in which a necklike structure only connects the outer
2148 <    monolayers, as predicted by the stalk hypothesis. We have found
2149 <    a new pathway of pore opening from stalks at high temperature: the
2150 <    elliptic stalk bends and contact between the ends of the arc-shaped
2151 <    stalk leads to pore opening. On the other hand, we have clarified
2152 <    that the pore-opening process at low temperature agrees with the
2153 <    modified stalk model: a pore is induced by contact between the inner
2154 <    monolayers inside the stalk. (C) 2001 American Institute of Physics.},
2145 >        simulation. Amphiphilic molecules spontaneously form vesicles with
2146 >        a bilayer structure. Two vesicles come into contact and form a stalk
2147 >        intermediate, in which a necklike structure only connects the outer
2148 >        monolayers, as predicted by the stalk hypothesis. We have found
2149 >        a new pathway of pore opening from stalks at high temperature: the
2150 >        elliptic stalk bends and contact between the ends of the arc-shaped
2151 >        stalk leads to pore opening. On the other hand, we have clarified
2152 >        that the pore-opening process at low temperature agrees with the
2153 >        modified stalk model: a pore is induced by contact between the inner
2154 >        monolayers inside the stalk. (C) 2001 American Institute of Physics.},
2155    annote = {491UW Times Cited:48 Cited References Count:25},
2156    issn = {0021-9606},
2157    uri = {<Go to ISI>://000172129300049},
2158   }
2159  
2064 @BOOK{Olver1986,
2065  title = {Applications of Lie groups to differential equatitons},
2066  publisher = {Springer},
2067  year = {1986},
2068  author = {P.J. Olver},
2069  address = {New York},
2070 }
2071
2160   @ARTICLE{Omelyan1998,
2161    author = {I. P. Omelyan},
2162    title = {On the numerical integration of motion for rigid polyatomics: The
2163 <    modified quaternion approach},
2163 >        modified quaternion approach},
2164    journal = {Computers in Physics},
2165    year = {1998},
2166    volume = {12},
# Line 2080 | Line 2168 | Encoding: GBK
2168    number = {1},
2169    month = {Jan-Feb},
2170    abstract = {A revised version of the quaternion approach for numerical integration
2171 <    of the equations of motion for rigid polyatomic molecules is proposed.
2172 <    The modified approach is based on a formulation of the quaternion
2173 <    dynamics with constraints. This allows one to resolve the rigidity
2174 <    problem rigorously using constraint forces. It is shown that the
2175 <    procedure for preservation of molecular rigidity can be realized
2176 <    particularly simply within the Verlet algorithm in velocity form.
2177 <    We demonstrate that the method presented leads to an improved numerical
2178 <    stability with respect to the usual quaternion rescaling scheme
2179 <    and it is roughly as good as the cumbersome atomic-constraint technique.
2180 <    (C) 1998 American Institute of Physics.},
2171 >        of the equations of motion for rigid polyatomic molecules is proposed.
2172 >        The modified approach is based on a formulation of the quaternion
2173 >        dynamics with constraints. This allows one to resolve the rigidity
2174 >        problem rigorously using constraint forces. It is shown that the
2175 >        procedure for preservation of molecular rigidity can be realized
2176 >        particularly simply within the Verlet algorithm in velocity form.
2177 >        We demonstrate that the method presented leads to an improved numerical
2178 >        stability with respect to the usual quaternion rescaling scheme
2179 >        and it is roughly as good as the cumbersome atomic-constraint technique.
2180 >        (C) 1998 American Institute of Physics.},
2181    annote = {Yx279 Times Cited:12 Cited References Count:28},
2182    issn = {0894-1866},
2183    uri = {<Go to ISI>://000072024300025},
# Line 2098 | Line 2186 | Encoding: GBK
2186   @ARTICLE{Omelyan1998a,
2187    author = {I. P. Omelyan},
2188    title = {Algorithm for numerical integration of the rigid-body equations of
2189 <    motion},
2189 >        motion},
2190    journal = {Physical Review E},
2191    year = {1998},
2192    volume = {58},
# Line 2106 | Line 2194 | Encoding: GBK
2194    number = {1},
2195    month = {Jul},
2196    abstract = {An algorithm for numerical integration of the rigid-body equations
2197 <    of motion is proposed. The algorithm uses the leapfrog scheme and
2198 <    the quantities involved are angular velocities and orientational
2199 <    variables that can be expressed in terms of either principal axes
2200 <    or quaternions. Due to specific features of the algorithm, orthonormality
2201 <    and unit norms of the orientational variables are integrals of motion,
2202 <    despite an approximate character of the produced trajectories. It
2203 <    is shown that the method presented appears to be the most efficient
2204 <    among all such algorithms known.},
2197 >        of motion is proposed. The algorithm uses the leapfrog scheme and
2198 >        the quantities involved are angular velocities and orientational
2199 >        variables that can be expressed in terms of either principal axes
2200 >        or quaternions. Due to specific features of the algorithm, orthonormality
2201 >        and unit norms of the orientational variables are integrals of motion,
2202 >        despite an approximate character of the produced trajectories. It
2203 >        is shown that the method presented appears to be the most efficient
2204 >        among all such algorithms known.},
2205    annote = {101XL Times Cited:8 Cited References Count:22},
2206    issn = {1063-651X},
2207    uri = {<Go to ISI>://000074893400151},
# Line 2122 | Line 2210 | Encoding: GBK
2210   @ARTICLE{Orlandi2006,
2211    author = {S. Orlandi and R. Berardi and J. Steltzer and C. Zannoni},
2212    title = {A Monte Carlo study of the mesophases formed by polar bent-shaped
2213 <    molecules},
2213 >        molecules},
2214    journal = {Journal of Chemical Physics},
2215    year = {2006},
2216    volume = {124},
# Line 2130 | Line 2218 | Encoding: GBK
2218    number = {12},
2219    month = {Mar 28},
2220    abstract = {Liquid crystal phases formed by bent-shaped (or #banana#) molecules
2221 <    are currently of great interest. Here we investigate by Monte Carlo
2222 <    computer simulations the phases formed by rigid banana molecules
2223 <    modeled combining three Gay-Berne sites and containing either one
2224 <    central or two lateral and transversal dipoles. We show that changing
2225 <    the dipole position and orientation has a profound effect on the
2226 <    mesophase stability and molecular organization. In particular, we
2227 <    find a uniaxial nematic phase only for off-center dipolar models
2228 <    and tilted phases only for the one with terminal dipoles. (c) 2006
2229 <    American Institute of Physics.},
2221 >        are currently of great interest. Here we investigate by Monte Carlo
2222 >        computer simulations the phases formed by rigid banana molecules
2223 >        modeled combining three Gay-Berne sites and containing either one
2224 >        central or two lateral and transversal dipoles. We show that changing
2225 >        the dipole position and orientation has a profound effect on the
2226 >        mesophase stability and molecular organization. In particular, we
2227 >        find a uniaxial nematic phase only for off-center dipolar models
2228 >        and tilted phases only for the one with terminal dipoles. (c) 2006
2229 >        American Institute of Physics.},
2230    annote = {028CP Times Cited:0 Cited References Count:42},
2231    issn = {0021-9606},
2232    uri = {<Go to ISI>://000236464000072},
# Line 2154 | Line 2242 | Encoding: GBK
2242    number = {6},
2243    month = {Nov},
2244    abstract = {Continuous, explicit Runge-Kutta methods with the minimal number of
2245 <    stages are considered. These methods are continuously differentiable
2246 <    if and only if one of the stages is the FSAL evaluation. A characterization
2247 <    of a subclass of these methods is developed for orders 3, 4, and
2248 <    5. It is shown how the free parameters of these methods can be used
2249 <    either to minimize the continuous truncation error coefficients
2250 <    or to maximize the stability region. As a representative for these
2251 <    methods the fifth-order method with minimized error coefficients
2252 <    is chosen, supplied with an error estimation method, and analysed
2253 <    by using the DETEST software. The results are compared with a similar
2254 <    implementation of the Dormand-Prince 5(4) pair with interpolant,
2255 <    showing a significant advantage in the new method for the chosen
2256 <    problems.},
2245 >        stages are considered. These methods are continuously differentiable
2246 >        if and only if one of the stages is the FSAL evaluation. A characterization
2247 >        of a subclass of these methods is developed for orders 3, 4, and
2248 >        5. It is shown how the free parameters of these methods can be used
2249 >        either to minimize the continuous truncation error coefficients
2250 >        or to maximize the stability region. As a representative for these
2251 >        methods the fifth-order method with minimized error coefficients
2252 >        is chosen, supplied with an error estimation method, and analysed
2253 >        by using the DETEST software. The results are compared with a similar
2254 >        implementation of the Dormand-Prince 5(4) pair with interpolant,
2255 >        showing a significant advantage in the new method for the chosen
2256 >        problems.},
2257    annote = {Ju936 Times Cited:25 Cited References Count:20},
2258    issn = {0196-5204},
2259    uri = {<Go to ISI>://A1992JU93600013},
# Line 2174 | Line 2262 | Encoding: GBK
2262   @ARTICLE{Palacios1998,
2263    author = {J. L. Garcia-Palacios and F. J. Lazaro},
2264    title = {Langevin-dynamics study of the dynamical properties of small magnetic
2265 <    particles},
2265 >        particles},
2266    journal = {Physical Review B},
2267    year = {1998},
2268    volume = {58},
# Line 2182 | Line 2270 | Encoding: GBK
2270    number = {22},
2271    month = {Dec 1},
2272    abstract = {The stochastic Landau-Lifshitz-Gilbert equation of motion for a classical
2273 <    magnetic moment is numerically solved (properly observing the customary
2274 <    interpretation of it as a Stratonovich stochastic differential equation),
2275 <    in order to study the dynamics of magnetic nanoparticles. The corresponding
2276 <    Langevin-dynamics approach allows for the study of the fluctuating
2277 <    trajectories of individual magnetic moments, where we have encountered
2278 <    remarkable phenomena in the overbarrier rotation process, such as
2279 <    crossing-back or multiple crossing of the potential barrier, rooted
2280 <    in the gyromagnetic nature of the system. Concerning averaged quantities,
2281 <    we study the linear dynamic response of the archetypal ensemble
2282 <    of noninteracting classical magnetic moments with axially symmetric
2283 <    magnetic anisotropy. The results are compared with different analytical
2284 <    expressions used to model the relaxation of nanoparticle ensembles,
2285 <    assessing their accuracy. It has been found that, among a number
2286 <    of heuristic expressions for the linear dynamic susceptibility,
2287 <    only the simple formula proposed by Shliomis and Stepanov matches
2288 <    the coarse features of the susceptibility reasonably. By comparing
2289 <    the numerical results with the asymptotic formula of Storonkin {Sov.
2290 <    Phys. Crystallogr. 30, 489 (1985) [Kristallografiya 30, 841 (1985)]},
2291 <    the effects of the intra-potential-well relaxation modes on the
2292 <    low-temperature longitudinal dynamic response have been assessed,
2293 <    showing their relatively small reflection in the susceptibility
2294 <    curves but their dramatic influence on the phase shifts. Comparison
2295 <    of the numerical results with the exact zero-damping expression
2296 <    for the transverse susceptibility by Garanin, Ishchenko, and Panina
2297 <    {Theor. Math. Phys. (USSR) 82, 169 (1990) [Teor. Mat. Fit. 82, 242
2298 <    (1990)]}, reveals a sizable contribution of the spread of the precession
2299 <    frequencies of the magnetic moment in the anisotropy field to the
2300 <    dynamic response at intermediate-to-high temperatures. [S0163-1829
2301 <    (98)00446-9].},
2273 >        magnetic moment is numerically solved (properly observing the customary
2274 >        interpretation of it as a Stratonovich stochastic differential equation),
2275 >        in order to study the dynamics of magnetic nanoparticles. The corresponding
2276 >        Langevin-dynamics approach allows for the study of the fluctuating
2277 >        trajectories of individual magnetic moments, where we have encountered
2278 >        remarkable phenomena in the overbarrier rotation process, such as
2279 >        crossing-back or multiple crossing of the potential barrier, rooted
2280 >        in the gyromagnetic nature of the system. Concerning averaged quantities,
2281 >        we study the linear dynamic response of the archetypal ensemble
2282 >        of noninteracting classical magnetic moments with axially symmetric
2283 >        magnetic anisotropy. The results are compared with different analytical
2284 >        expressions used to model the relaxation of nanoparticle ensembles,
2285 >        assessing their accuracy. It has been found that, among a number
2286 >        of heuristic expressions for the linear dynamic susceptibility,
2287 >        only the simple formula proposed by Shliomis and Stepanov matches
2288 >        the coarse features of the susceptibility reasonably. By comparing
2289 >        the numerical results with the asymptotic formula of Storonkin {Sov.
2290 >        Phys. Crystallogr. 30, 489 (1985) [Kristallografiya 30, 841 (1985)]},
2291 >        the effects of the intra-potential-well relaxation modes on the
2292 >        low-temperature longitudinal dynamic response have been assessed,
2293 >        showing their relatively small reflection in the susceptibility
2294 >        curves but their dramatic influence on the phase shifts. Comparison
2295 >        of the numerical results with the exact zero-damping expression
2296 >        for the transverse susceptibility by Garanin, Ishchenko, and Panina
2297 >        {Theor. Math. Phys. (USSR) 82, 169 (1990) [Teor. Mat. Fit. 82, 242
2298 >        (1990)]}, reveals a sizable contribution of the spread of the precession
2299 >        frequencies of the magnetic moment in the anisotropy field to the
2300 >        dynamic response at intermediate-to-high temperatures. [S0163-1829
2301 >        (98)00446-9].},
2302    annote = {146XW Times Cited:66 Cited References Count:45},
2303    issn = {0163-1829},
2304    uri = {<Go to ISI>://000077460000052},
2305 + }
2306 +
2307 + @ARTICLE{Parr1995,
2308 +  author = {T. J. Parr and R. W. Quong},
2309 +  title = {Antlr - a Predicated-Ll(K) Parser Generator},
2310 +  journal = {Software-Practice \& Experience},
2311 +  year = {1995},
2312 +  volume = {25},
2313 +  pages = {789-810},
2314 +  number = {7},
2315 +  month = {Jul},
2316 +  abstract = {Despite the parsing power of LR/LALR algorithms, e.g. YACC, programmers
2317 +        often choose to write recursive-descent parsers by hand to obtain
2318 +        increased flexibility, better error handling, and ease of debugging.
2319 +        We introduce ANTLR, a public-domain parser generator that combines
2320 +        the flexibility of hand-coded parsing with the convenience of a
2321 +        parser generator, which is a component of PCCTS. ANTLR has many
2322 +        features that make it easier to use than other language tools. Most
2323 +        important, ANTLR provides predicates which let the programmer systematically
2324 +        direct the parse via arbitrary expressions using semantic and syntactic
2325 +        context; in practice, the use of predicates eliminates the need
2326 +        to hand-tweak the ANTLR output, even for difficult parsing problems.
2327 +        ANTLR also integrates the description of lexical and syntactic analysis,
2328 +        accepts LL(k) grammars for k > 1 with extended BNF notation, and
2329 +        can automatically generate abstract syntax trees. ANTLR is widely
2330 +        used, with over 1000 registered industrial and academic users in
2331 +        37 countries. It has been ported to many popular systems such as
2332 +        the PC, Macintosh, and a variety of UNIX platforms; a commercial
2333 +        C++ front-end has been developed as a result of one of our industrial
2334 +        collaborations.},
2335 +  annote = {Rk104 Times Cited:19 Cited References Count:10},
2336 +  issn = {0038-0644},
2337 +  uri = {<Go to ISI>://A1995RK10400004},
2338   }
2339  
2340   @ARTICLE{Pastor1988,
# Line 2247 | Line 2368 | Encoding: GBK
2368   @ARTICLE{Perram1985,
2369    author = {J. W. Perram and M. S. Wertheim},
2370    title = {Statistical-Mechanics of Hard Ellipsoids .1. Overlap Algorithm and
2371 <    the Contact Function},
2371 >        the Contact Function},
2372    journal = {Journal of Computational Physics},
2373    year = {1985},
2374    volume = {58},
# Line 2258 | Line 2379 | Encoding: GBK
2379    uri = {<Go to ISI>://A1985AKB9300008},
2380   }
2381  
2261 @ARTICLE{Rotne1969,
2262  author = {F. Perrin},
2263  title = {Variational treatment of hydrodynamic interaction in polymers},
2264  journal = {J. Chem. Phys.},
2265  year = {1969},
2266  volume = {50},
2267  pages = {4831¨C4837},
2268 }
2269
2270 @ARTICLE{Perrin1936,
2271  author = {F. Perrin},
2272  title = {Mouvement brownien d'un ellipsoid(II). Rotation libre et depolarisation
2273    des fluorescences. Translation et diffusion de moleculese ellipsoidales},
2274  journal = {J. Phys. Radium},
2275  year = {1936},
2276  volume = {7},
2277  pages = {1-11},
2278 }
2279
2280 @ARTICLE{Perrin1934,
2281  author = {F. Perrin},
2282  title = {Mouvement brownien d'un ellipsoid(I). Dispersion dielectrique pour
2283    des molecules ellipsoidales},
2284  journal = {J. Phys. Radium},
2285  year = {1934},
2286  volume = {5},
2287  pages = {497-511},
2288 }
2289
2382   @ARTICLE{Petrache1998,
2383    author = {H. I. Petrache and S. Tristram-Nagle and J. F. Nagle},
2384    title = {Fluid phase structure of EPC and DMPC bilayers},
# Line 2297 | Line 2389 | Encoding: GBK
2389    number = {1},
2390    month = {Sep},
2391    abstract = {X-ray diffraction data taken at high instrumental resolution were
2392 <    obtained for EPC and DMPC under various osmotic pressures, primarily
2393 <    at T = 30 degrees C. The headgroup thickness D-HH was obtained from
2394 <    relative electron density profiles. By using volumetric results
2395 <    and by comparing to gel phase DPPC we obtain areas A(EPC)(F) = 69.4
2396 <    +/- 1.1 Angstrom(2) and A(DMPC)(F) = 59.7 +/- 0.2 Angstrom(2). The
2397 <    analysis also gives estimates for the areal compressibility K-A.
2398 <    The A(F) results lead to other structural results regarding membrane
2399 <    thickness and associated waters. Using the recently determined absolute
2400 <    electrons density profile of DPPC, the AF results also lead to absolute
2401 <    electron density profiles and absolute continuous transforms \F(q)\
2402 <    for EPC and DMPC, Limited measurements of temperature dependence
2403 <    show directly that fluctuations increase with increasing temperature
2404 <    and that a small decrease in bending modulus K-c accounts for the
2405 <    increased water spacing reported by Simon et al. (1995) Biophys.
2406 <    J. 69, 1473-1483. (C) 1998 Elsevier Science Ireland Ltd. All rights
2407 <    reserved.},
2392 >        obtained for EPC and DMPC under various osmotic pressures, primarily
2393 >        at T = 30 degrees C. The headgroup thickness D-HH was obtained from
2394 >        relative electron density profiles. By using volumetric results
2395 >        and by comparing to gel phase DPPC we obtain areas A(EPC)(F) = 69.4
2396 >        +/- 1.1 Angstrom(2) and A(DMPC)(F) = 59.7 +/- 0.2 Angstrom(2). The
2397 >        analysis also gives estimates for the areal compressibility K-A.
2398 >        The A(F) results lead to other structural results regarding membrane
2399 >        thickness and associated waters. Using the recently determined absolute
2400 >        electrons density profile of DPPC, the AF results also lead to absolute
2401 >        electron density profiles and absolute continuous transforms \F(q)\
2402 >        for EPC and DMPC, Limited measurements of temperature dependence
2403 >        show directly that fluctuations increase with increasing temperature
2404 >        and that a small decrease in bending modulus K-c accounts for the
2405 >        increased water spacing reported by Simon et al. (1995) Biophys.
2406 >        J. 69, 1473-1483. (C) 1998 Elsevier Science Ireland Ltd. All rights
2407 >        reserved.},
2408    annote = {130AT Times Cited:98 Cited References Count:39},
2409    issn = {0009-3084},
2410    uri = {<Go to ISI>://000076497600007},
2411   }
2412  
2321 @ARTICLE{Powles1973,
2322  author = {J.~G. Powles},
2323  title = {A general ellipsoid can not always serve as a modle for the rotational
2324    diffusion properties of arbitrary shaped rigid molecules},
2325  journal = {Advan. Phys.},
2326  year = {1973},
2327  volume = {22},
2328  pages = {1-56},
2329 }
2330
2413   @ARTICLE{Recio2004,
2414    author = {J. Fernandez-Recio and M. Totrov and R. Abagyan},
2415    title = {Identification of protein-protein interaction sites from docking
2416 <    energy landscapes},
2416 >        energy landscapes},
2417    journal = {Journal of Molecular Biology},
2418    year = {2004},
2419    volume = {335},
# Line 2339 | Line 2421 | Encoding: GBK
2421    number = {3},
2422    month = {Jan 16},
2423    abstract = {Protein recognition is one of the most challenging and intriguing
2424 <    problems in structural biology. Despite all the available structural,
2425 <    sequence and biophysical information about protein-protein complexes,
2426 <    the physico-chemical patterns, if any, that make a protein surface
2427 <    likely to be involved in protein-protein interactions, remain elusive.
2428 <    Here, we apply protein docking simulations and analysis of the interaction
2429 <    energy landscapes to identify protein-protein interaction sites.
2430 <    The new protocol for global docking based on multi-start global
2431 <    energy optimization of an allatom model of the ligand, with detailed
2432 <    receptor potentials and atomic solvation parameters optimized in
2433 <    a training set of 24 complexes, explores the conformational space
2434 <    around the whole receptor without restrictions. The ensembles of
2435 <    the rigid-body docking solutions generated by the simulations were
2436 <    subsequently used to project the docking energy landscapes onto
2437 <    the protein surfaces. We found that highly populated low-energy
2438 <    regions consistently corresponded to actual binding sites. The procedure
2439 <    was validated on a test set of 21 known protein-protein complexes
2440 <    not used in the training set. As much as 81% of the predicted high-propensity
2441 <    patch residues were located correctly in the native interfaces.
2442 <    This approach can guide the design of mutations on the surfaces
2443 <    of proteins, provide geometrical details of a possible interaction,
2444 <    and help to annotate protein surfaces in structural proteomics.
2445 <    (C) 2003 Elsevier Ltd. All rights reserved.},
2424 >        problems in structural biology. Despite all the available structural,
2425 >        sequence and biophysical information about protein-protein complexes,
2426 >        the physico-chemical patterns, if any, that make a protein surface
2427 >        likely to be involved in protein-protein interactions, remain elusive.
2428 >        Here, we apply protein docking simulations and analysis of the interaction
2429 >        energy landscapes to identify protein-protein interaction sites.
2430 >        The new protocol for global docking based on multi-start global
2431 >        energy optimization of an allatom model of the ligand, with detailed
2432 >        receptor potentials and atomic solvation parameters optimized in
2433 >        a training set of 24 complexes, explores the conformational space
2434 >        around the whole receptor without restrictions. The ensembles of
2435 >        the rigid-body docking solutions generated by the simulations were
2436 >        subsequently used to project the docking energy landscapes onto
2437 >        the protein surfaces. We found that highly populated low-energy
2438 >        regions consistently corresponded to actual binding sites. The procedure
2439 >        was validated on a test set of 21 known protein-protein complexes
2440 >        not used in the training set. As much as 81% of the predicted high-propensity
2441 >        patch residues were located correctly in the native interfaces.
2442 >        This approach can guide the design of mutations on the surfaces
2443 >        of proteins, provide geometrical details of a possible interaction,
2444 >        and help to annotate protein surfaces in structural proteomics.
2445 >        (C) 2003 Elsevier Ltd. All rights reserved.},
2446    annote = {763GQ Times Cited:21 Cited References Count:59},
2447    issn = {0022-2836},
2448    uri = {<Go to ISI>://000188066900016},
# Line 2369 | Line 2451 | Encoding: GBK
2451   @ARTICLE{Reddy2006,
2452    author = {R. A. Reddy and C. Tschierske},
2453    title = {Bent-core liquid crystals: polar order, superstructural chirality
2454 <    and spontaneous desymmetrisation in soft matter systems},
2454 >        and spontaneous desymmetrisation in soft matter systems},
2455    journal = {Journal of Materials Chemistry},
2456    year = {2006},
2457    volume = {16},
2458    pages = {907-961},
2459    number = {10},
2460    abstract = {An overview on the recent developments in the field of liquid crystalline
2461 <    bent-core molecules (so-called banana liquid crystals) is given.
2462 <    After some basic issues, dealing with general aspects of the systematisation
2463 <    of the mesophases, development of polar order and chirality in this
2464 <    class of LC systems and explaining some general structure-property
2465 <    relationships, we focus on fascinating new developments in this
2466 <    field, such as modulated, undulated and columnar phases, so-called
2467 <    B7 phases, phase biaxiality, ferroelectric and antiferroelectric
2468 <    polar order in smectic and columnar phases, amplification and switching
2469 <    of chirality and the spontaneous formation of superstructural and
2470 <    supramolecular chirality.},
2461 >        bent-core molecules (so-called banana liquid crystals) is given.
2462 >        After some basic issues, dealing with general aspects of the systematisation
2463 >        of the mesophases, development of polar order and chirality in this
2464 >        class of LC systems and explaining some general structure-property
2465 >        relationships, we focus on fascinating new developments in this
2466 >        field, such as modulated, undulated and columnar phases, so-called
2467 >        B7 phases, phase biaxiality, ferroelectric and antiferroelectric
2468 >        polar order in smectic and columnar phases, amplification and switching
2469 >        of chirality and the spontaneous formation of superstructural and
2470 >        supramolecular chirality.},
2471    annote = {021NS Times Cited:2 Cited References Count:316},
2472    issn = {0959-9428},
2473    uri = {<Go to ISI>://000235990500001},
# Line 2401 | Line 2483 | Encoding: GBK
2483    number = {5},
2484    month = {Sep 8},
2485    abstract = {Backward error analysis has become an important tool for understanding
2486 <    the long time behavior of numerical integration methods. This is
2487 <    true in particular for the integration of Hamiltonian systems where
2488 <    backward error analysis can be used to show that a symplectic method
2489 <    will conserve energy over exponentially long periods of time. Such
2490 <    results are typically based on two aspects of backward error analysis:
2491 <    (i) It can be shown that the modified vector fields have some qualitative
2492 <    properties which they share with the given problem and (ii) an estimate
2493 <    is given for the difference between the best interpolating vector
2494 <    field and the numerical method. These aspects have been investigated
2495 <    recently, for example, by Benettin and Giorgilli in [J. Statist.
2496 <    Phys., 74 (1994), pp. 1117-1143], by Hairer in [Ann. Numer. Math.,
2497 <    1 (1994), pp. 107-132], and by Hairer and Lubich in [Numer. Math.,
2498 <    76 (1997), pp. 441-462]. In this paper we aim at providing a unifying
2499 <    framework and a simplification of the existing results and corresponding
2500 <    proofs. Our approach to backward error analysis is based on a simple
2501 <    recursive definition of the modified vector fields that does not
2502 <    require explicit Taylor series expansion of the numerical method
2503 <    and the corresponding flow maps as in the above-cited works. As
2504 <    an application we discuss the long time integration of chaotic Hamiltonian
2505 <    systems and the approximation of time averages along numerically
2506 <    computed trajectories.},
2486 >        the long time behavior of numerical integration methods. This is
2487 >        true in particular for the integration of Hamiltonian systems where
2488 >        backward error analysis can be used to show that a symplectic method
2489 >        will conserve energy over exponentially long periods of time. Such
2490 >        results are typically based on two aspects of backward error analysis:
2491 >        (i) It can be shown that the modified vector fields have some qualitative
2492 >        properties which they share with the given problem and (ii) an estimate
2493 >        is given for the difference between the best interpolating vector
2494 >        field and the numerical method. These aspects have been investigated
2495 >        recently, for example, by Benettin and Giorgilli in [J. Statist.
2496 >        Phys., 74 (1994), pp. 1117-1143], by Hairer in [Ann. Numer. Math.,
2497 >        1 (1994), pp. 107-132], and by Hairer and Lubich in [Numer. Math.,
2498 >        76 (1997), pp. 441-462]. In this paper we aim at providing a unifying
2499 >        framework and a simplification of the existing results and corresponding
2500 >        proofs. Our approach to backward error analysis is based on a simple
2501 >        recursive definition of the modified vector fields that does not
2502 >        require explicit Taylor series expansion of the numerical method
2503 >        and the corresponding flow maps as in the above-cited works. As
2504 >        an application we discuss the long time integration of chaotic Hamiltonian
2505 >        systems and the approximation of time averages along numerically
2506 >        computed trajectories.},
2507    annote = {237HV Times Cited:43 Cited References Count:41},
2508    issn = {0036-1429},
2509    uri = {<Go to ISI>://000082650600010},
# Line 2436 | Line 2518 | Encoding: GBK
2518    pages = {5093-5098},
2519    number = {48},
2520    abstract = {The recent literature in the field of liquid crystals shows that banana-shaped
2521 <    mesogenic materials represent a bewitching and stimulating field
2522 <    of research that is interesting both academically and in terms of
2523 <    applications. Numerous topics are open to investigation in this
2524 <    area because of the rich phenomenology and new possibilities that
2525 <    these materials offer. The principal concepts in this area are reviewed
2526 <    along with recent results. In addition, new directions to stimulate
2527 <    further research activities are highlighted.},
2521 >        mesogenic materials represent a bewitching and stimulating field
2522 >        of research that is interesting both academically and in terms of
2523 >        applications. Numerous topics are open to investigation in this
2524 >        area because of the rich phenomenology and new possibilities that
2525 >        these materials offer. The principal concepts in this area are reviewed
2526 >        along with recent results. In addition, new directions to stimulate
2527 >        further research activities are highlighted.},
2528    annote = {990XA Times Cited:3 Cited References Count:72},
2529    issn = {0959-9428},
2530    uri = {<Go to ISI>://000233775500001},
2531   }
2532  
2533 + @ARTICLE{Roux1991,
2534 +  author = {B. Roux and M. Karplus},
2535 +  title = {Ion-Transport in a Gramicidin-Like Channel - Dynamics and Mobility},
2536 +  journal = {Journal of Physical Chemistry},
2537 +  year = {1991},
2538 +  volume = {95},
2539 +  pages = {4856-4868},
2540 +  number = {12},
2541 +  month = {Jun 13},
2542 +  abstract = {The mobility of water, Na+. and K+ has been calculated inside a periodic
2543 +        poly-(L,D)-alanine beta-helix, a model for the interior of the gramicidin
2544 +        channel. Because of the different dynamical regimes for the three
2545 +        species (high barrier for Na+, low barrier for K+, almost free diffusion
2546 +        for water), different methods are used to calculate the mobilities.
2547 +        By use of activated dynamics and a potential of mean force determined
2548 +        previously (Roux, B.; Karplus, M. Biophys. J. 1991, 59, 961), the
2549 +        barrier crossing rate of Na+ ion is determined. The motion of Na+
2550 +        at the transition state is controlled by local interactions and
2551 +        collisions with the neighboring carbonyls and the two nearest water
2552 +        molecules. There are significant deviations from transition-state
2553 +        theory; the transmission coefficient is equal to 0.11. The water
2554 +        and K+ motions are found to be well described by a diffusive model;
2555 +        the motion of K+ appears to be controlled by the diffusion of water.
2556 +        The time-dependent friction functions of Na+ and K+ ions in the
2557 +        periodic beta-helix are calculated and analyzed by using a generalized
2558 +        Langevin equation approach. Both Na+ and K+ suffer many rapid collisions,
2559 +        and their dynamics is overdamped and noninertial. Thus, the selectivity
2560 +        sequence of ions in the beta-helix is not influenced strongly by
2561 +        their masses.},
2562 +  annote = {Fr756 Times Cited:97 Cited References Count:65},
2563 +  issn = {0022-3654},
2564 +  uri = {<Go to ISI>://A1991FR75600049},
2565 + }
2566 +
2567   @ARTICLE{Roy2005,
2568    author = {A. Roy and N. V. Madhusudana},
2569    title = {A frustrated packing model for the B-6-B-1-SmAP(A) sequence of phases
2570 <    in banana shaped molecules},
2570 >        in banana shaped molecules},
2571    journal = {European Physical Journal E},
2572    year = {2005},
2573    volume = {18},
# Line 2459 | Line 2575 | Encoding: GBK
2575    number = {3},
2576    month = {Nov},
2577    abstract = {A vast majority of compounds with bent core or banana shaped molecules
2578 <    exhibit the phase sequence B-6-B-1-B-2 as the chain length is increased
2579 <    in a homologous series. The B-6 phase has an intercalated fluid
2580 <    lamellar structure with a layer spacing of half the molecular length.
2581 <    The B-1 phase has a two dimensionally periodic rectangular columnar
2582 <    structure. The B-2 phase has a monolayer fluid lamellar structure
2583 <    with molecules tilted with respect to the layer normal. Neglecting
2584 <    the tilt order of the molecules in the B-2 phase, we have developed
2585 <    a frustrated packing model to describe this phase sequence qualitatively.
2586 <    The model has some analogy with that of the frustrated smectics
2587 <    exhibited by highly polar rod like molecules.},
2578 >        exhibit the phase sequence B-6-B-1-B-2 as the chain length is increased
2579 >        in a homologous series. The B-6 phase has an intercalated fluid
2580 >        lamellar structure with a layer spacing of half the molecular length.
2581 >        The B-1 phase has a two dimensionally periodic rectangular columnar
2582 >        structure. The B-2 phase has a monolayer fluid lamellar structure
2583 >        with molecules tilted with respect to the layer normal. Neglecting
2584 >        the tilt order of the molecules in the B-2 phase, we have developed
2585 >        a frustrated packing model to describe this phase sequence qualitatively.
2586 >        The model has some analogy with that of the frustrated smectics
2587 >        exhibited by highly polar rod like molecules.},
2588    annote = {985FW Times Cited:0 Cited References Count:30},
2589    issn = {1292-8941},
2590    uri = {<Go to ISI>://000233363300002},
# Line 2477 | Line 2593 | Encoding: GBK
2593   @ARTICLE{Ryckaert1977,
2594    author = {J. P. Ryckaert and G. Ciccotti and H. J. C. Berendsen},
2595    title = {Numerical-Integration of Cartesian Equations of Motion of a System
2596 <    with Constraints - Molecular-Dynamics of N-Alkanes},
2596 >        with Constraints - Molecular-Dynamics of N-Alkanes},
2597    journal = {Journal of Computational Physics},
2598    year = {1977},
2599    volume = {23},
# Line 2491 | Line 2607 | Encoding: GBK
2607   @ARTICLE{Sagui1999,
2608    author = {C. Sagui and T. A. Darden},
2609    title = {Molecular dynamics simulations of biomolecules: Long-range electrostatic
2610 <    effects},
2610 >        effects},
2611    journal = {Annual Review of Biophysics and Biomolecular Structure},
2612    year = {1999},
2613    volume = {28},
2614    pages = {155-179},
2615    abstract = {Current computer simulations of biomolecules typically make use of
2616 <    classical molecular dynamics methods, as a very large number (tens
2617 <    to hundreds of thousands) of atoms are involved over timescales
2618 <    of many nanoseconds. The methodology for treating short-range bonded
2619 <    and van der Waals interactions has matured. However, long-range
2620 <    electrostatic interactions still represent a bottleneck in simulations.
2621 <    In this article, we introduce the basic issues for an accurate representation
2622 <    of the relevant electrostatic interactions. In spite of the huge
2623 <    computational time demanded by most biomolecular systems, it is
2624 <    no longer necessary to resort to uncontrolled approximations such
2625 <    as the use of cutoffs. In particular, we discuss the Ewald summation
2626 <    methods, the fast particle mesh methods, and the fast multipole
2627 <    methods. We also review recent efforts to understand the role of
2628 <    boundary conditions in systems with long-range interactions, and
2629 <    conclude with a short perspective on future trends.},
2616 >        classical molecular dynamics methods, as a very large number (tens
2617 >        to hundreds of thousands) of atoms are involved over timescales
2618 >        of many nanoseconds. The methodology for treating short-range bonded
2619 >        and van der Waals interactions has matured. However, long-range
2620 >        electrostatic interactions still represent a bottleneck in simulations.
2621 >        In this article, we introduce the basic issues for an accurate representation
2622 >        of the relevant electrostatic interactions. In spite of the huge
2623 >        computational time demanded by most biomolecular systems, it is
2624 >        no longer necessary to resort to uncontrolled approximations such
2625 >        as the use of cutoffs. In particular, we discuss the Ewald summation
2626 >        methods, the fast particle mesh methods, and the fast multipole
2627 >        methods. We also review recent efforts to understand the role of
2628 >        boundary conditions in systems with long-range interactions, and
2629 >        conclude with a short perspective on future trends.},
2630    annote = {213KJ Times Cited:126 Cited References Count:73},
2631    issn = {1056-8700},
2632    uri = {<Go to ISI>://000081271400008},
# Line 2519 | Line 2635 | Encoding: GBK
2635   @ARTICLE{Sandu1999,
2636    author = {A. Sandu and T. Schlick},
2637    title = {Masking resonance artifacts in force-splitting methods for biomolecular
2638 <    simulations by extrapolative Langevin dynamics},
2638 >        simulations by extrapolative Langevin dynamics},
2639    journal = {Journal of Computational Physics},
2640    year = {1999},
2641    volume = {151},
# Line 2527 | Line 2643 | Encoding: GBK
2643    number = {1},
2644    month = {May 1},
2645    abstract = {Numerical resonance artifacts have become recognized recently as a
2646 <    limiting factor to increasing the timestep in multiple-timestep
2647 <    (MTS) biomolecular dynamics simulations. At certain timesteps correlated
2648 <    to internal motions (e.g., 5 fs, around half the period of the fastest
2649 <    bond stretch, T-min), visible inaccuracies or instabilities can
2650 <    occur. Impulse-MTS schemes are vulnerable to these resonance errors
2651 <    since large energy pulses are introduced to the governing dynamics
2652 <    equations when the slow forces are evaluated. We recently showed
2653 <    that such resonance artifacts can be masked significantly by applying
2654 <    extrapolative splitting to stochastic dynamics. Theoretical and
2655 <    numerical analyses of force-splitting integrators based on the Verlet
2656 <    discretization are reported here for linear models to explain these
2657 <    observations and to suggest how to construct effective integrators
2658 <    for biomolecular dynamics that balance stability with accuracy.
2659 <    Analyses for Newtonian dynamics demonstrate the severe resonance
2660 <    patterns of the Impulse splitting, with this severity worsening
2661 <    with the outer timestep. Delta t: Constant Extrapolation is generally
2662 <    unstable, but the disturbances do not grow with Delta t. Thus. the
2663 <    stochastic extrapolative combination can counteract generic instabilities
2664 <    and largely alleviate resonances with a sufficiently strong Langevin
2665 <    heat-bath coupling (gamma), estimates for which are derived here
2666 <    based on the fastest and slowest motion periods. These resonance
2667 <    results generally hold for nonlinear test systems: a water tetramer
2668 <    and solvated protein. Proposed related approaches such as Extrapolation/Correction
2669 <    and Midpoint Extrapolation work better than Constant Extrapolation
2670 <    only for timesteps less than T-min/2. An effective extrapolative
2671 <    stochastic approach for biomolecules that balances long-timestep
2672 <    stability with good accuracy for the fast subsystem is then applied
2673 <    to a biomolecule using a three-class partitioning: the medium forces
2674 <    are treated by Midpoint Extrapolation via position Verlet, and the
2675 <    slow forces are incorporated by Constant Extrapolation. The resulting
2676 <    algorithm (LN) performs well on a solvated protein system in terms
2677 <    of thermodynamic properties and yields an order of magnitude speedup
2678 <    with respect to single-timestep Langevin trajectories. Computed
2679 <    spectral density functions also show how the Newtonian modes can
2680 <    be approximated by using a small gamma in the range Of 5-20 ps(-1).
2681 <    (C) 1999 Academic Press.},
2646 >        limiting factor to increasing the timestep in multiple-timestep
2647 >        (MTS) biomolecular dynamics simulations. At certain timesteps correlated
2648 >        to internal motions (e.g., 5 fs, around half the period of the fastest
2649 >        bond stretch, T-min), visible inaccuracies or instabilities can
2650 >        occur. Impulse-MTS schemes are vulnerable to these resonance errors
2651 >        since large energy pulses are introduced to the governing dynamics
2652 >        equations when the slow forces are evaluated. We recently showed
2653 >        that such resonance artifacts can be masked significantly by applying
2654 >        extrapolative splitting to stochastic dynamics. Theoretical and
2655 >        numerical analyses of force-splitting integrators based on the Verlet
2656 >        discretization are reported here for linear models to explain these
2657 >        observations and to suggest how to construct effective integrators
2658 >        for biomolecular dynamics that balance stability with accuracy.
2659 >        Analyses for Newtonian dynamics demonstrate the severe resonance
2660 >        patterns of the Impulse splitting, with this severity worsening
2661 >        with the outer timestep. Delta t: Constant Extrapolation is generally
2662 >        unstable, but the disturbances do not grow with Delta t. Thus. the
2663 >        stochastic extrapolative combination can counteract generic instabilities
2664 >        and largely alleviate resonances with a sufficiently strong Langevin
2665 >        heat-bath coupling (gamma), estimates for which are derived here
2666 >        based on the fastest and slowest motion periods. These resonance
2667 >        results generally hold for nonlinear test systems: a water tetramer
2668 >        and solvated protein. Proposed related approaches such as Extrapolation/Correction
2669 >        and Midpoint Extrapolation work better than Constant Extrapolation
2670 >        only for timesteps less than T-min/2. An effective extrapolative
2671 >        stochastic approach for biomolecules that balances long-timestep
2672 >        stability with good accuracy for the fast subsystem is then applied
2673 >        to a biomolecule using a three-class partitioning: the medium forces
2674 >        are treated by Midpoint Extrapolation via position Verlet, and the
2675 >        slow forces are incorporated by Constant Extrapolation. The resulting
2676 >        algorithm (LN) performs well on a solvated protein system in terms
2677 >        of thermodynamic properties and yields an order of magnitude speedup
2678 >        with respect to single-timestep Langevin trajectories. Computed
2679 >        spectral density functions also show how the Newtonian modes can
2680 >        be approximated by using a small gamma in the range Of 5-20 ps(-1).
2681 >        (C) 1999 Academic Press.},
2682    annote = {194FM Times Cited:14 Cited References Count:32},
2683    issn = {0021-9991},
2684    uri = {<Go to ISI>://000080181500004},
2685 + }
2686 +
2687 + @ARTICLE{Sasaki2004,
2688 +  author = {Y. Sasaki and R. Shukla and B. D. Smith},
2689 +  title = {Facilitated phosphatidylserine flip-flop across vesicle and cell
2690 +        membranes using urea-derived synthetic translocases},
2691 +  journal = {Organic \& Biomolecular Chemistry},
2692 +  year = {2004},
2693 +  volume = {2},
2694 +  pages = {214-219},
2695 +  number = {2},
2696 +  abstract = {Tris(2-aminoethyl) amine derivatives with appended urea and sulfonamide
2697 +        groups are shown to facilitate the translocation of fluorescent
2698 +        phospholipid probes and endogenous phosphatidylserine across vesicle
2699 +        and erythrocyte cell membranes. The synthetic translocases appear
2700 +        to operate by binding to the phospholipid head groups and forming
2701 +        lipophilic supramolecular complexes which diffuse through the non-polar
2702 +        interior of the bilayer membrane.},
2703 +  annote = {760PX Times Cited:8 Cited References Count:25},
2704 +  issn = {1477-0520},
2705 +  uri = {<Go to ISI>://000187843800012},
2706   }
2707  
2708   @ARTICLE{Satoh1996,
2709    author = {K. Satoh and S. Mita and S. Kondo},
2710    title = {Monte Carlo simulations using the dipolar Gay-Berne model: Effect
2711 <    of terminal dipole moment on mesophase formation},
2711 >        of terminal dipole moment on mesophase formation},
2712    journal = {Chemical Physics Letters},
2713    year = {1996},
2714    volume = {255},
# Line 2579 | Line 2716 | Encoding: GBK
2716    number = {1-3},
2717    month = {Jun 7},
2718    abstract = {The effects of dipole-dipole interaction on mesophase formation are
2719 <    investigated with a Monte Carlo simulation using the dipolar Gay-Berne
2720 <    potential. It is shown that the dipole moment at the end of a molecule
2721 <    causes a shift in the nematic-isotropic transition toward higher
2722 <    temperature and a spread of the temperature range of the nematic
2723 <    phase and that layer structures with various interdigitations are
2724 <    formed in the smectic phase.},
2719 >        investigated with a Monte Carlo simulation using the dipolar Gay-Berne
2720 >        potential. It is shown that the dipole moment at the end of a molecule
2721 >        causes a shift in the nematic-isotropic transition toward higher
2722 >        temperature and a spread of the temperature range of the nematic
2723 >        phase and that layer structures with various interdigitations are
2724 >        formed in the smectic phase.},
2725    annote = {Uq975 Times Cited:32 Cited References Count:33},
2726    issn = {0009-2614},
2727    uri = {<Go to ISI>://A1996UQ97500017},
2728   }
2729  
2730 + @ARTICLE{Schaps1999,
2731 +  author = {G. L. Schaps},
2732 +  title = {Compiler construction with ANTLR and Java - Tools for building tools},
2733 +  journal = {Dr Dobbs Journal},
2734 +  year = {1999},
2735 +  volume = {24},
2736 +  pages = {84-+},
2737 +  number = {3},
2738 +  month = {Mar},
2739 +  annote = {163EC Times Cited:0 Cited References Count:0},
2740 +  issn = {1044-789X},
2741 +  uri = {<Go to ISI>://000078389200023},
2742 + }
2743 +
2744   @ARTICLE{Shen2002,
2745    author = {M. Y. Shen and K. F. Freed},
2746    title = {Long time dynamics of met-enkephalin: Comparison of explicit and
2747 <    implicit solvent models},
2747 >        implicit solvent models},
2748    journal = {Biophysical Journal},
2749    year = {2002},
2750    volume = {82},
# Line 2601 | Line 2752 | Encoding: GBK
2752    number = {4},
2753    month = {Apr},
2754    abstract = {Met-enkephalin is one of the smallest opiate peptides. Yet, its dynamical
2755 <    structure and receptor docking mechanism are still not well understood.
2756 <    The conformational dynamics of this neuron peptide in liquid water
2757 <    are studied here by using all-atom molecular dynamics (MID) and
2758 <    implicit water Langevin dynamics (LD) simulations with AMBER potential
2759 <    functions and the three-site transferable intermolecular potential
2760 <    (TIP3P) model for water. To achieve the same simulation length in
2761 <    physical time, the full MID simulations require 200 times as much
2762 <    CPU time as the implicit water LID simulations. The solvent hydrophobicity
2763 <    and dielectric behavior are treated in the implicit solvent LD simulations
2764 <    by using a macroscopic solvation potential, a single dielectric
2765 <    constant, and atomic friction coefficients computed using the accessible
2766 <    surface area method with the TIP3P model water viscosity as determined
2767 <    here from MID simulations for pure TIP3P water. Both the local and
2768 <    the global dynamics obtained from the implicit solvent LD simulations
2769 <    agree very well with those from the explicit solvent MD simulations.
2770 <    The simulations provide insights into the conformational restrictions
2771 <    that are associated with the bioactivity of the opiate peptide dermorphin
2772 <    for the delta-receptor.},
2755 >        structure and receptor docking mechanism are still not well understood.
2756 >        The conformational dynamics of this neuron peptide in liquid water
2757 >        are studied here by using all-atom molecular dynamics (MID) and
2758 >        implicit water Langevin dynamics (LD) simulations with AMBER potential
2759 >        functions and the three-site transferable intermolecular potential
2760 >        (TIP3P) model for water. To achieve the same simulation length in
2761 >        physical time, the full MID simulations require 200 times as much
2762 >        CPU time as the implicit water LID simulations. The solvent hydrophobicity
2763 >        and dielectric behavior are treated in the implicit solvent LD simulations
2764 >        by using a macroscopic solvation potential, a single dielectric
2765 >        constant, and atomic friction coefficients computed using the accessible
2766 >        surface area method with the TIP3P model water viscosity as determined
2767 >        here from MID simulations for pure TIP3P water. Both the local and
2768 >        the global dynamics obtained from the implicit solvent LD simulations
2769 >        agree very well with those from the explicit solvent MD simulations.
2770 >        The simulations provide insights into the conformational restrictions
2771 >        that are associated with the bioactivity of the opiate peptide dermorphin
2772 >        for the delta-receptor.},
2773    annote = {540MH Times Cited:36 Cited References Count:45},
2774    issn = {0006-3495},
2775    uri = {<Go to ISI>://000174932400010},
# Line 2641 | Line 2792 | Encoding: GBK
2792   @ARTICLE{Shimada1993,
2793    author = {J. Shimada and H. Kaneko and T. Takada},
2794    title = {Efficient Calculations of Coulombic Interactions in Biomolecular
2795 <    Simulations with Periodic Boundary-Conditions},
2795 >        Simulations with Periodic Boundary-Conditions},
2796    journal = {Journal of Computational Chemistry},
2797    year = {1993},
2798    volume = {14},
# Line 2649 | Line 2800 | Encoding: GBK
2800    number = {7},
2801    month = {Jul},
2802    abstract = {To make improved treatments of electrostatic interactions in biomacromolecular
2803 <    simulations, two possibilities are considered. The first is the
2804 <    famous particle-particle and particle-mesh (PPPM) method developed
2805 <    by Hockney and Eastwood, and the second is a new one developed here
2806 <    in their spirit but by the use of the multipole expansion technique
2807 <    suggested by Ladd. It is then numerically found that the new PPPM
2808 <    method gives more accurate results for a two-particle system at
2809 <    small separation of particles. Preliminary numerical examination
2810 <    of the various computational methods for a single configuration
2811 <    of a model BPTI-water system containing about 24,000 particles indicates
2812 <    that both of the PPPM methods give far more accurate values with
2813 <    reasonable computational cost than do the conventional truncation
2814 <    methods. It is concluded the two PPPM methods are nearly comparable
2815 <    in overall performance for the many-particle systems, although the
2816 <    first method has the drawback that the accuracy in the total electrostatic
2817 <    energy is not high for configurations of charged particles randomly
2818 <    generated.},
2803 >        simulations, two possibilities are considered. The first is the
2804 >        famous particle-particle and particle-mesh (PPPM) method developed
2805 >        by Hockney and Eastwood, and the second is a new one developed here
2806 >        in their spirit but by the use of the multipole expansion technique
2807 >        suggested by Ladd. It is then numerically found that the new PPPM
2808 >        method gives more accurate results for a two-particle system at
2809 >        small separation of particles. Preliminary numerical examination
2810 >        of the various computational methods for a single configuration
2811 >        of a model BPTI-water system containing about 24,000 particles indicates
2812 >        that both of the PPPM methods give far more accurate values with
2813 >        reasonable computational cost than do the conventional truncation
2814 >        methods. It is concluded the two PPPM methods are nearly comparable
2815 >        in overall performance for the many-particle systems, although the
2816 >        first method has the drawback that the accuracy in the total electrostatic
2817 >        energy is not high for configurations of charged particles randomly
2818 >        generated.},
2819    annote = {Lh164 Times Cited:27 Cited References Count:47},
2820    issn = {0192-8651},
2821    uri = {<Go to ISI>://A1993LH16400011},
# Line 2680 | Line 2831 | Encoding: GBK
2831    number = {24},
2832    month = {Dec 20},
2833    abstract = {The best simple method for Newtonian molecular dynamics is indisputably
2834 <    the leapfrog Stormer-Verlet method. The appropriate generalization
2835 <    to simple Langevin dynamics is unclear. An analysis is presented
2836 <    comparing an 'impulse method' (kick; fluctuate; kick), the 1982
2837 <    method of van Gunsteren and Berendsen, and the Brunger-Brooks-Karplus
2838 <    (BBK) method. It is shown how the impulse method and the van Gunsteren-Berendsen
2839 <    methods can be implemented as efficiently as the BBK method. Other
2840 <    considerations suggest that the impulse method is the best basic
2841 <    method for simple Langevin dynamics, with the van Gunsteren-Berendsen
2842 <    method a close contender.},
2834 >        the leapfrog Stormer-Verlet method. The appropriate generalization
2835 >        to simple Langevin dynamics is unclear. An analysis is presented
2836 >        comparing an 'impulse method' (kick; fluctuate; kick), the 1982
2837 >        method of van Gunsteren and Berendsen, and the Brunger-Brooks-Karplus
2838 >        (BBK) method. It is shown how the impulse method and the van Gunsteren-Berendsen
2839 >        methods can be implemented as efficiently as the BBK method. Other
2840 >        considerations suggest that the impulse method is the best basic
2841 >        method for simple Langevin dynamics, with the van Gunsteren-Berendsen
2842 >        method a close contender.},
2843    annote = {633RX Times Cited:8 Cited References Count:22},
2844    issn = {0026-8976},
2845    uri = {<Go to ISI>://000180297200014},
# Line 2697 | Line 2848 | Encoding: GBK
2848   @ARTICLE{Skeel1997,
2849    author = {R. D. Skeel and G. H. Zhang and T. Schlick},
2850    title = {A family of symplectic integrators: Stability, accuracy, and molecular
2851 <    dynamics applications},
2851 >        dynamics applications},
2852    journal = {Siam Journal on Scientific Computing},
2853    year = {1997},
2854    volume = {18},
# Line 2705 | Line 2856 | Encoding: GBK
2856    number = {1},
2857    month = {Jan},
2858    abstract = {The following integration methods for special second-order ordinary
2859 <    differential equations are studied: leapfrog, implicit midpoint,
2860 <    trapezoid, Stormer-Verlet, and Cowell-Numerov. We show that all
2861 <    are members, or equivalent to members, of a one-parameter family
2862 <    of schemes. Some methods have more than one common form, and we
2863 <    discuss a systematic enumeration of these forms. We also present
2864 <    a stability and accuracy analysis based on the idea of ''modified
2865 <    equations'' and a proof of symplecticness. It follows that Cowell-Numerov
2866 <    and ''LIM2'' (a method proposed by Zhang and Schlick) are symplectic.
2867 <    A different interpretation of the values used by these integrators
2868 <    leads to higher accuracy and better energy conservation. Hence,
2869 <    we suggest that the straightforward analysis of energy conservation
2870 <    is misleading.},
2859 >        differential equations are studied: leapfrog, implicit midpoint,
2860 >        trapezoid, Stormer-Verlet, and Cowell-Numerov. We show that all
2861 >        are members, or equivalent to members, of a one-parameter family
2862 >        of schemes. Some methods have more than one common form, and we
2863 >        discuss a systematic enumeration of these forms. We also present
2864 >        a stability and accuracy analysis based on the idea of ''modified
2865 >        equations'' and a proof of symplecticness. It follows that Cowell-Numerov
2866 >        and ''LIM2'' (a method proposed by Zhang and Schlick) are symplectic.
2867 >        A different interpretation of the values used by these integrators
2868 >        leads to higher accuracy and better energy conservation. Hence,
2869 >        we suggest that the straightforward analysis of energy conservation
2870 >        is misleading.},
2871    annote = {We981 Times Cited:30 Cited References Count:35},
2872    issn = {1064-8275},
2873    uri = {<Go to ISI>://A1997WE98100012},
# Line 2724 | Line 2875 | Encoding: GBK
2875  
2876   @ARTICLE{Tao2005,
2877    author = {Y. G. Tao and W. K. {den Otter} and J. T. Padding and J. K. G. Dhont
2878 <    and W. J. Briels},
2878 >        and W. J. Briels},
2879    title = {Brownian dynamics simulations of the self- and collective rotational
2880 <    diffusion coefficients of rigid long thin rods},
2880 >        diffusion coefficients of rigid long thin rods},
2881    journal = {Journal of Chemical Physics},
2882    year = {2005},
2883    volume = {122},
# Line 2734 | Line 2885 | Encoding: GBK
2885    number = {24},
2886    month = {Jun 22},
2887    abstract = {Recently a microscopic theory for the dynamics of suspensions of long
2888 <    thin rigid rods was presented, confirming and expanding the well-known
2889 <    theory by Doi and Edwards [The Theory of Polymer Dynamics (Clarendon,
2890 <    Oxford, 1986)] and Kuzuu [J. Phys. Soc. Jpn. 52, 3486 (1983)]. Here
2891 <    this theory is put to the test by comparing it against computer
2892 <    simulations. A Brownian dynamics simulation program was developed
2893 <    to follow the dynamics of the rods, with a length over a diameter
2894 <    ratio of 60, on the Smoluchowski time scale. The model accounts
2895 <    for excluded volume interactions between rods, but neglects hydrodynamic
2896 <    interactions. The self-rotational diffusion coefficients D-r(phi)
2897 <    of the rods were calculated by standard methods and by a new, more
2898 <    efficient method based on calculating average restoring torques.
2899 <    Collective decay of orientational order was calculated by means
2900 <    of equilibrium and nonequilibrium simulations. Our results show
2901 <    that, for the currently accessible volume fractions, the decay times
2902 <    in both cases are virtually identical. Moreover, the observed decay
2903 <    of diffusion coefficients with volume fraction is much quicker than
2904 <    predicted by the theory, which is attributed to an oversimplification
2905 <    of dynamic correlations in the theory. (c) 2005 American Institute
2906 <    of Physics.},
2888 >        thin rigid rods was presented, confirming and expanding the well-known
2889 >        theory by Doi and Edwards [The Theory of Polymer Dynamics (Clarendon,
2890 >        Oxford, 1986)] and Kuzuu [J. Phys. Soc. Jpn. 52, 3486 (1983)]. Here
2891 >        this theory is put to the test by comparing it against computer
2892 >        simulations. A Brownian dynamics simulation program was developed
2893 >        to follow the dynamics of the rods, with a length over a diameter
2894 >        ratio of 60, on the Smoluchowski time scale. The model accounts
2895 >        for excluded volume interactions between rods, but neglects hydrodynamic
2896 >        interactions. The self-rotational diffusion coefficients D-r(phi)
2897 >        of the rods were calculated by standard methods and by a new, more
2898 >        efficient method based on calculating average restoring torques.
2899 >        Collective decay of orientational order was calculated by means
2900 >        of equilibrium and nonequilibrium simulations. Our results show
2901 >        that, for the currently accessible volume fractions, the decay times
2902 >        in both cases are virtually identical. Moreover, the observed decay
2903 >        of diffusion coefficients with volume fraction is much quicker than
2904 >        predicted by the theory, which is attributed to an oversimplification
2905 >        of dynamic correlations in the theory. (c) 2005 American Institute
2906 >        of Physics.},
2907    annote = {943DN Times Cited:3 Cited References Count:26},
2908    issn = {0021-9606},
2909    uri = {<Go to ISI>://000230332400077},
2910   }
2911  
2761 @BOOK{Tolman1979,
2762  title = {The Principles of Statistical Mechanics},
2763  publisher = {Dover Publications, Inc.},
2764  year = {1979},
2765  author = {R.~C. Tolman},
2766  address = {New York},
2767  chapter = {2},
2768  pages = {19-22},
2769 }
2770
2912   @ARTICLE{Tu1995,
2913    author = {K. Tu and D. J. Tobias and M. L. Klein},
2914    title = {Constant pressure and temperature molecular dynamics simulation of
2915 <    a fully hydrated liquid crystal phase dipalmitoylphosphatidylcholine
2916 <    bilayer},
2915 >        a fully hydrated liquid crystal phase dipalmitoylphosphatidylcholine
2916 >        bilayer},
2917    journal = {Biophysical Journal},
2918    year = {1995},
2919    volume = {69},
# Line 2780 | Line 2921 | Encoding: GBK
2921    number = {6},
2922    month = {Dec},
2923    abstract = {We report a constant pressure and temperature molecular dynamics simulation
2924 <    of a fully hydrated liquid crystal (L(alpha) phase bilayer of dipalmitoylphosphatidylcholine
2925 <    at 50 degrees C and 28 water molecules/lipid. We have shown that
2926 <    the bilayer is stable throughout the 1550-ps simulation and have
2927 <    demonstrated convergence of the system dimensions. Several important
2928 <    aspects of the bilayer structure have been investigated and compared
2929 <    favorably with experimental results. For example, the average positions
2930 <    of specific carbon atoms along the bilayer normal agree well with
2931 <    neutron diffraction data, and the electron density profile is in
2932 <    accord with x-ray diffraction results. The hydrocarbon chain deuterium
2933 <    order parameters agree reasonably well with NMR results for the
2934 <    middles of the chains, but the simulation predicts too much order
2935 <    at the chain ends. In spite of the deviations in the order parameters,
2936 <    the hydrocarbon chain packing density appears to be essentially
2937 <    correct, inasmuch as the area/lipid and bilayer thickness are in
2938 <    agreement with the most refined experimental estimates. The deuterium
2939 <    order parameters for the glycerol and choline groups, as well as
2940 <    the phosphorus chemical shift anisotropy, are in qualitative agreement
2941 <    with those extracted from NMR measurements.},
2924 >        of a fully hydrated liquid crystal (L(alpha) phase bilayer of dipalmitoylphosphatidylcholine
2925 >        at 50 degrees C and 28 water molecules/lipid. We have shown that
2926 >        the bilayer is stable throughout the 1550-ps simulation and have
2927 >        demonstrated convergence of the system dimensions. Several important
2928 >        aspects of the bilayer structure have been investigated and compared
2929 >        favorably with experimental results. For example, the average positions
2930 >        of specific carbon atoms along the bilayer normal agree well with
2931 >        neutron diffraction data, and the electron density profile is in
2932 >        accord with x-ray diffraction results. The hydrocarbon chain deuterium
2933 >        order parameters agree reasonably well with NMR results for the
2934 >        middles of the chains, but the simulation predicts too much order
2935 >        at the chain ends. In spite of the deviations in the order parameters,
2936 >        the hydrocarbon chain packing density appears to be essentially
2937 >        correct, inasmuch as the area/lipid and bilayer thickness are in
2938 >        agreement with the most refined experimental estimates. The deuterium
2939 >        order parameters for the glycerol and choline groups, as well as
2940 >        the phosphorus chemical shift anisotropy, are in qualitative agreement
2941 >        with those extracted from NMR measurements.},
2942    annote = {Tv018 Times Cited:108 Cited References Count:34},
2943    issn = {0006-3495},
2944    uri = {<Go to ISI>://A1995TV01800037},
# Line 2813 | Line 2954 | Encoding: GBK
2954    number = {3},
2955    month = {Aug 1},
2956    abstract = {The Trotter factorization of the Liouville propagator is used to generate
2957 <    new reversible molecular dynamics integrators. This strategy is
2958 <    applied to derive reversible reference system propagator algorithms
2959 <    (RESPA) that greatly accelerate simulations of systems with a separation
2960 <    of time scales or with long range forces. The new algorithms have
2961 <    all of the advantages of previous RESPA integrators but are reversible,
2962 <    and more stable than those methods. These methods are applied to
2963 <    a set of paradigmatic systems and are shown to be superior to earlier
2964 <    methods. It is shown how the new RESPA methods are related to predictor-corrector
2965 <    integrators. Finally, we show how these methods can be used to accelerate
2966 <    the integration of the equations of motion of systems with Nose
2967 <    thermostats.},
2957 >        new reversible molecular dynamics integrators. This strategy is
2958 >        applied to derive reversible reference system propagator algorithms
2959 >        (RESPA) that greatly accelerate simulations of systems with a separation
2960 >        of time scales or with long range forces. The new algorithms have
2961 >        all of the advantages of previous RESPA integrators but are reversible,
2962 >        and more stable than those methods. These methods are applied to
2963 >        a set of paradigmatic systems and are shown to be superior to earlier
2964 >        methods. It is shown how the new RESPA methods are related to predictor-corrector
2965 >        integrators. Finally, we show how these methods can be used to accelerate
2966 >        the integration of the equations of motion of systems with Nose
2967 >        thermostats.},
2968    annote = {Je891 Times Cited:680 Cited References Count:19},
2969    issn = {0021-9606},
2970    uri = {<Go to ISI>://A1992JE89100044},
2971   }
2972  
2832 @BOOK{Varadarajan1974,
2833  title = {Lie groups, Lie algebras, and their representations},
2834  publisher = {Prentice-Hall},
2835  year = {1974},
2836  author = {V.S. Varadarajan},
2837  address = {New York},
2838 }
2839
2840 @ARTICLE{Wegener1979,
2841  author = {W.~A. Wegener, V.~J. Koester and R.~M. Dowben},
2842  title = {A general ellipsoid can not always serve as a modle for the rotational
2843    diffusion properties of arbitrary shaped rigid molecules},
2844  journal = {Proc. Natl. Acad. Sci.},
2845  year = {1979},
2846  volume = {76},
2847  pages = {6356-6360},
2848  number = {12},
2849 }
2850
2973   @ARTICLE{Withers2003,
2974    author = {I. M. Withers},
2975    title = {Effects of longitudinal quadrupoles on the phase behavior of a Gay-Berne
2976 <    fluid},
2976 >        fluid},
2977    journal = {Journal of Chemical Physics},
2978    year = {2003},
2979    volume = {119},
# Line 2859 | Line 2981 | Encoding: GBK
2981    number = {19},
2982    month = {Nov 15},
2983    abstract = {The effects of longitudinal quadrupole moments on the formation of
2984 <    liquid crystalline phases are studied by means of constant NPT Monte
2985 <    Carlo simulation methods. The popular Gay-Berne model mesogen is
2986 <    used as the reference fluid, which displays the phase sequences
2987 <    isotropic-smectic A-smectic B and isotropic-smectic B at high (T*=2.0)
2988 <    and low (T*=1.5) temperatures, respectively. With increasing quadrupole
2989 <    magnitude the smectic phases are observed to be stabilized with
2990 <    respect to the isotropic liquid, while the smectic B is destabilized
2991 <    with respect to the smectic A. At the lower temperature, a sufficiently
2992 <    large quadrupole magnitude results in the injection of the smectic
2993 <    A phase into the phase sequence and the replacement of the smectic
2994 <    B phase by the tilted smectic J phase. The nematic phase is also
2995 <    injected into the phase sequence at both temperatures considered,
2996 <    and ultimately for sufficiently large quadrupole magnitudes no coherent
2997 <    layered structures were observed. The stabilization of the smectic
2998 <    A phase supports the commonly held belief that, while the inclusion
2999 <    of polar groups is not a prerequisite for the formation of the smectic
3000 <    A phase, quadrupolar interactions help to increase the temperature
3001 <    and pressure range for which the smectic A phase is observed. The
3002 <    quality of the layered structure is worsened with increasing quadrupole
3003 <    magnitude. This behavior, along with the injection of the nematic
3004 <    phase into the phase sequence, indicate that the general tendency
3005 <    of the quadrupolar interactions is to destabilize the layered structure.
3006 <    A pressure dependence upon the smectic layer spacing is observed.
3007 <    This behavior is in much closer agreement with experimental findings
3008 <    than has been observed previously for nonpolar Gay-Berne and hard
3009 <    spherocylinder models. (C) 2003 American Institute of Physics.},
2984 >        liquid crystalline phases are studied by means of constant NPT Monte
2985 >        Carlo simulation methods. The popular Gay-Berne model mesogen is
2986 >        used as the reference fluid, which displays the phase sequences
2987 >        isotropic-smectic A-smectic B and isotropic-smectic B at high (T*=2.0)
2988 >        and low (T*=1.5) temperatures, respectively. With increasing quadrupole
2989 >        magnitude the smectic phases are observed to be stabilized with
2990 >        respect to the isotropic liquid, while the smectic B is destabilized
2991 >        with respect to the smectic A. At the lower temperature, a sufficiently
2992 >        large quadrupole magnitude results in the injection of the smectic
2993 >        A phase into the phase sequence and the replacement of the smectic
2994 >        B phase by the tilted smectic J phase. The nematic phase is also
2995 >        injected into the phase sequence at both temperatures considered,
2996 >        and ultimately for sufficiently large quadrupole magnitudes no coherent
2997 >        layered structures were observed. The stabilization of the smectic
2998 >        A phase supports the commonly held belief that, while the inclusion
2999 >        of polar groups is not a prerequisite for the formation of the smectic
3000 >        A phase, quadrupolar interactions help to increase the temperature
3001 >        and pressure range for which the smectic A phase is observed. The
3002 >        quality of the layered structure is worsened with increasing quadrupole
3003 >        magnitude. This behavior, along with the injection of the nematic
3004 >        phase into the phase sequence, indicate that the general tendency
3005 >        of the quadrupolar interactions is to destabilize the layered structure.
3006 >        A pressure dependence upon the smectic layer spacing is observed.
3007 >        This behavior is in much closer agreement with experimental findings
3008 >        than has been observed previously for nonpolar Gay-Berne and hard
3009 >        spherocylinder models. (C) 2003 American Institute of Physics.},
3010    annote = {738EF Times Cited:3 Cited References Count:43},
3011    issn = {0021-9606},
3012    uri = {<Go to ISI>://000186273200027},
# Line 2893 | Line 3015 | Encoding: GBK
3015   @ARTICLE{Wolf1999,
3016    author = {D. Wolf and P. Keblinski and S. R. Phillpot and J. Eggebrecht},
3017    title = {Exact method for the simulation of Coulombic systems by spherically
3018 <    truncated, pairwise r(-1) summation},
3018 >        truncated, pairwise r(-1) summation},
3019    journal = {Journal of Chemical Physics},
3020    year = {1999},
3021    volume = {110},
# Line 2901 | Line 3023 | Encoding: GBK
3023    number = {17},
3024    month = {May 1},
3025    abstract = {Based on a recent result showing that the net Coulomb potential in
3026 <    condensed ionic systems is rather short ranged, an exact and physically
3027 <    transparent method permitting the evaluation of the Coulomb potential
3028 <    by direct summation over the r(-1) Coulomb pair potential is presented.
3029 <    The key observation is that the problems encountered in determining
3030 <    the Coulomb energy by pairwise, spherically truncated r(-1) summation
3031 <    are a direct consequence of the fact that the system summed over
3032 <    is practically never neutral. A simple method is developed that
3033 <    achieves charge neutralization wherever the r(-1) pair potential
3034 <    is truncated. This enables the extraction of the Coulomb energy,
3035 <    forces, and stresses from a spherically truncated, usually charged
3036 <    environment in a manner that is independent of the grouping of the
3037 <    pair terms. The close connection of our approach with the Ewald
3038 <    method is demonstrated and exploited, providing an efficient method
3039 <    for the simulation of even highly disordered ionic systems by direct,
3040 <    pairwise r(-1) summation with spherical truncation at rather short
3041 <    range, i.e., a method which fully exploits the short-ranged nature
3042 <    of the interactions in ionic systems. The method is validated by
3043 <    simulations of crystals, liquids, and interfacial systems, such
3044 <    as free surfaces and grain boundaries. (C) 1999 American Institute
3045 <    of Physics. [S0021-9606(99)51517-1].},
3026 >        condensed ionic systems is rather short ranged, an exact and physically
3027 >        transparent method permitting the evaluation of the Coulomb potential
3028 >        by direct summation over the r(-1) Coulomb pair potential is presented.
3029 >        The key observation is that the problems encountered in determining
3030 >        the Coulomb energy by pairwise, spherically truncated r(-1) summation
3031 >        are a direct consequence of the fact that the system summed over
3032 >        is practically never neutral. A simple method is developed that
3033 >        achieves charge neutralization wherever the r(-1) pair potential
3034 >        is truncated. This enables the extraction of the Coulomb energy,
3035 >        forces, and stresses from a spherically truncated, usually charged
3036 >        environment in a manner that is independent of the grouping of the
3037 >        pair terms. The close connection of our approach with the Ewald
3038 >        method is demonstrated and exploited, providing an efficient method
3039 >        for the simulation of even highly disordered ionic systems by direct,
3040 >        pairwise r(-1) summation with spherical truncation at rather short
3041 >        range, i.e., a method which fully exploits the short-ranged nature
3042 >        of the interactions in ionic systems. The method is validated by
3043 >        simulations of crystals, liquids, and interfacial systems, such
3044 >        as free surfaces and grain boundaries. (C) 1999 American Institute
3045 >        of Physics. [S0021-9606(99)51517-1].},
3046    annote = {189PD Times Cited:70 Cited References Count:34},
3047    issn = {0021-9606},
3048    uri = {<Go to ISI>://000079913000008},
# Line 2940 | Line 3062 | Encoding: GBK
3062    uri = {<Go to ISI>://A1990EJ79800009},
3063   }
3064  
2943 @Book{Frenkel1996,
2944  author =   {D. Frenkel and B. Smit},
2945  title =    {Understanding Molecular Simulation : From Algorithms
2946                  to Applications},
2947  publisher =    {Academic Press},
2948  year =     1996,
2949  address =  {New York}
2950 }

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines