Hi Quentin, thank you so much for the recommendation - I will look into packmol!!<br /><div>Yes, the idea would be to introduce a highly pressurized flow of CO2 into our system - to be honest, my advisor and I have not really worked out how pressurized it should be, so I will look into it.</div><div><br /></div><div>Michela</div><br /><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Friday, February 14, 2025 at 4:50:57 AM UTC-5 Quentin Pessemesse wrote:<br/></div><blockquote class="gmail_quote" style="margin: 0 0 0 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;"><div>Hello,</div><div>About building initial configurations, packmol is a great free option to generate "disordered" inital structures. It's very easy to use.</div><div>Regarding your CO2 cell, you have 36 CO2 molecules in a cell with a volume of around ~
2,0 × 10^-27 m³. This is an ideal gas pressure of around 3000 bars, so you are indeed looking at a very dense phase of CO2 and not CO2 gas (of course this is no longer an ideal gas and the 3000 bar figure is not really meaningful, but your initial density is ~1300 kg/
m^-3

, which is between the densities of liquid and solid CO2). Are you sure this is what you mean to look at?</div><div>I hope someone can chime in regarding the barostat as it's not something I am very familiar with.</div><div>Take care :)</div><div>Quentin</div><div class="gmail_quote"><div dir="auto" class="gmail_attr">Le jeudi 13 février 2025 à 15:08:02 UTC+1, Michela Benazzi a écrit :<br></div><blockquote class="gmail_quote" style="margin:0 0 0 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div>Hello everyone,</div><div><br></div><div>I have posted here before about my simulations involving liquid metal + CO2. I realized that my extremely high T, P were arising once I added CO2 instead of vacuum on top of a previously converged 2D periodic slab of liquid metal.</div><div><br></div><div>I am now trying to run DFT on CO2 alone at 1200 K. I attached the files I used before and after adding the Grimme D3 correction: I thought it would "fix it" by accounting for long-range vdW forces that I was not accounting for before, but nothing has changed. Is there anything else I am missing?<br><br>I have gotten comments before on how my initial CO2 system is too ordered (I generate it with Python): I am just curious about how others prepare an initial gas structure before using DFT to get optimized geometry. Should I run GEO_OPT first?</div><div><br></div><div>Is that even the main problem with my work?<br><br>Thank you so much! Best,</div><div><br></div><div>Michela</div></blockquote></div></blockquote></div>

<p></p>

-- <br />
You received this message because you are subscribed to the Google Groups "cp2k" group.<br />
To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp2k+unsubscribe@googlegroups.com">cp2k+unsubscribe@googlegroups.com</a>.<br />
To view this discussion visit <a href="https://groups.google.com/d/msgid/cp2k/309a7441-6e54-4beb-a231-a5457a7bae90n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/309a7441-6e54-4beb-a231-a5457a7bae90n%40googlegroups.com</a>.<br />