Comments (3)
Hi
After looking into post update scripts for pressure correction, I found out
that we dont have to do much to enable pressure correction with relative
entropy and we can use (almost) the same post_update scripts for relative
entropy. Only thing, that is needed to take care is that instead of adding
correction to .dpot, the correction should be added to .pot.new, the potential
file written by csg_reupdate.
Now, I am a bit confused about how the pressure correction is done in the case
of mixtures. For mixtures, user would specify how much each interaction should
contribute to the total pressure of a system, i.e., the target pressure for an
interaction. Then votca should determine the actual pressure contribution of
each specific interaction and add corresponding correction based on error in
each case.
However, from the script, I understand that gromacs' g_energy is used to
compute CG pressure from current CG-MD simulations. Then this pressure is used
to determine the difference between user specified target pressure. However,
g_energy is computing the total pressure and not the contribution from a
specified CG interaction. Therefore, in the case of mixture, wouldnt it cause
excessive correction (multiple counting of error) in each interaction?
I think, one solution to this issue would be to specify pressure correction
frequency for each interaction such that only one interaction is corrected at
any given iteration. Am I correct? Am I just re-iterating what you guys have
already thought about :)
Thanks
Sikandar
Original comment by [email protected]
on 13 Nov 2013 at 5:29
from votca.
Well, yes, for multi-component systems, one would need to calculate the partial
pressures. Espresso is able to do that, but Gromacs isn't. We always had the
plan to add this feature to Gromacs, but as pressure correction got less and
less popular with the years, so it never happened. (We should add a warning to
csg_inverse though!)
Considering that all components contribute equally according to the number of
pairs, one can calculate an estimate for the partial pressures. Technically the
ratio, N_pairs_XY/N_pairs_total can then be used as the scale of the pressure
correction. However, this approximation is pretty crude.
Original comment by [email protected]
on 13 Nov 2013 at 5:56
from votca.
one more option to get partial pressures from gromacs would be to do mdrun
-rerun with energy exclusion group. With this option one can exclude all
interactions except for a given pair of particles and then the pressure
from g_energy will have virial component due to interaction of interest.
Original comment by [email protected]
on 13 Nov 2013 at 6:14
from votca.
Related Issues (20)
- Reduce warnings XML options not set in pyxtp HOT 5
- 2 tests fail HOT 1
- COHSEX approximation in GW
- Aux-AO matrix elements of the MM reaction field
- libint2 warning on ubuntu:devel HOT 6
- Votca fails to build with Boost 1.83 HOT 2
- Failing CI builds HOT 7
- boost/progress.hpp:23:3: error: This header is deprecated and will be removed HOT 2
- [2023] configure first says that it found Python, and then fails because it says that it didn't find Python HOT 2
- [2023] Several tests fail HOT 7
- Forbid segment types in KMC
- duplicated names make csg_stat segfault
- xtp-tutorials fail on rawhide (f41) HOT 8
- Can't set source for GROMACS binaries HOT 4
- Changes for the v2024.1 JOSS paper HOT 2
- Update CITATION.cff
- Gromacs Relative entropy HOT 5
- IBI with Lammps HOT 4
- installation problem on ubuntu HOT 4
- CI fails on FreeBSD HOT 5
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from votca.