2004
DOI: 10.1063/1.1802672
|View full text |Cite
|
Sign up to set email alerts
|

Some issues on the calculation of interfacial properties by molecular simulation

Abstract: Some of the pitfalls that may befall molecular simulations of interfaces are discussed. They are all related to the calculation of the pressure tensor profiles, which are needed in order to compute surface tensions. We focus on three controversial points: (1) the calculation of the pressure tensor profiles for polyatomic systems, in particular, when the SHAKE algorithm is employed, (2) the addition of long-range corrections to compensate the truncation of the potential, and (3) the importance of including adeq… Show more

Help me understand this report

Search citation statements

Order By: Relevance

Paper Sections

Select...
3
1
1

Citation Types

4
73
0
1

Year Published

2006
2006
2016
2016

Publication Types

Select...
8

Relationship

1
7

Authors

Journals

citations
Cited by 82 publications
(78 citation statements)
references
References 22 publications
4
73
0
1
Order By: Relevance
“…We also recall that the evaluation of the pressure tensor for molecular fluids can be tricky, particularly in the case of molecules subject to geometrical constraints, as pointed out by Duque and Vega. 71 Obviously, the test-area method can be easily incorporated into a molecular dynamics simulation program or, as is done here, used for analyzing a set of trajectories generated from molecular dynamics. As for the estimates of the errors in the computed averages, both methods appear to provide a similar level of accuracy for all models of water at all the values of temperature considered here.…”
Section: Discussionmentioning
confidence: 99%
“…We also recall that the evaluation of the pressure tensor for molecular fluids can be tricky, particularly in the case of molecules subject to geometrical constraints, as pointed out by Duque and Vega. 71 Obviously, the test-area method can be easily incorporated into a molecular dynamics simulation program or, as is done here, used for analyzing a set of trajectories generated from molecular dynamics. As for the estimates of the errors in the computed averages, both methods appear to provide a similar level of accuracy for all models of water at all the values of temperature considered here.…”
Section: Discussionmentioning
confidence: 99%
“…The GCMC simulations for the cylindrical pore were done employing: (1) a simulation cell consisting of a cylinder with a diameter ranging from 1.0 to 4.0 nm and a fixed length of 10.0 nm, (2) a cutoff radius of at least 6 times the collision diameter of the fluid molecules (r LJ ), 46 (3) 1.5 Â 10 5 MC steps for equilibrating the system and 2.0 Â 10 6 MC steps for data collection, (4) an equal a priori probability to displace, insert or remove a molecule in each simulation step and (5) periodic boundary conditions in the z-direction.…”
Section: Simulation Detailsmentioning
confidence: 99%
“…These values are chosen in order to reduce the truncation and system size effects involved in the phase equilibrium and interfacial tension calculations (see Refs. [60][61][62][63][64] for a complete discussion of truncation effects).…”
Section: Simulation Detailsmentioning
confidence: 99%