<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
</head>
<body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">
Dear Laurent,
<div class=""><br class="">
</div>
<div class="">In Environ you can introduce a pressure term defined on a variety of interfaces. In particular, one can specify an electronic interface, similar to the one of Cococcioni et al., but based off a reformulation (from our J Chem Phys paper of 2012).
 But you can also specify an ionic interface, defined following the soft-sphere approach of Fisicaro et al, JTCT doi: 10.1021/acs.jctc.7b00375. The first case would introduce a ‘pressure' on the electrons via a contribution to the Kohn-Sham potential, the second
 case would introduce explicit forces on the ions. </div>
<div class=""><br class="">
</div>
<div class="">There is a third interface, dubbed system, which corresponds to analytic functions (a plane, a cylinder, a sphere). Right now this latter kind of interface is rigid, has no direct effect on forces or Kohn sham potential, but it could be modified
 easily. </div>
<div class=""><br class="">
</div>
<div class="">Best,</div>
<div class=""><br class="">
</div>
<div class="">Oliviero<br class="">
<div><br class="">
<blockquote type="cite" class="">
<div class="">On May 27, 2021, at 3:59 PM, laurent pizzagalli <<a href="mailto:laurent.pizzagalli@univ-poitiers.fr" class="">laurent.pizzagalli@univ-poitiers.fr</a>> wrote:</div>
<br class="Apple-interchange-newline">
<div class="">
<div class="">Dear Oliviero, <br class="">
<br class="">
thanks for your email. I know the paper you mentioned, but I thought (wrongly apparently) that this approach was only implemented in cp.x (the &PRESS-AI namelist). Also I was not really sure this could be used for applying a uniaxial stress on an isolated system.
 The pressure is accounted for by considering the PV term in the enthalpy, with a volume computed from rho. In the case of a uniaxial compression, one should instead define a surface where the stress should be present? I did not see how it could be done. This
 is why I considered forces on ions instead of forces on electrons. <br class="">
It would be interesting to study whether these different approaches could imply different mechanical behaviors !
<br class="">
<br class="">
Best regards<br class="">
<br class="">
Laurent <br class="">
  <br class="">
<br class="">
<div class="moz-cite-prefix">On 27/05/2021 22:34, Andreussi, Oliviero wrote:<br class="">
</div>
<blockquote type="cite" cite="mid:D439A443-44A9-4D97-812D-A62310C0CC9A@unt.edu" class="">
Dear Laurent,
<div class=""><br class="">
</div>
<div class="">I am not fully sure about the details of your implementation and I am not commenting on the best place to put your keywords (as Paolo said, Environ adopts an external file for its own input). I mostly wanted to mention that Environ allows you
 to introduce an external pressure into the system. With the right settings this could be used to apply a compressive pressure on rigid planes above or below a slab or, using a cylindrical symmetry, on a nanowire (but also to a isolated system, as per Cococcioni
 et al, PRL 10.1103/PhysRevLett.94.145501).</div>
<div class=""><br class="">
</div>
<div class="">Best,</div>
<div class=""><br class="">
</div>
<div class="">Oliviero</div>
<div class=""> <br class="">
<div class=""><br class="">
<div class="">
<div dir="auto" style="caret-color: rgb(0, 0, 0); letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">
<div dir="auto" style="caret-color: rgb(0, 0, 0); letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">
<div dir="auto" style="word-wrap: break-word;
                  -webkit-nbsp-mode: space; line-break:
                  after-white-space;" class="">
<div style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;" class="">
Oliviero Andreussi<br class="">
--<br class="">
Assistant Professor<br class="">
Department of Physics<br class="">
University of North Texas<br class="">
<a href="mailto:oliviero.andreussi@unt.edu" class="" moz-do-not-send="true">Email: oliviero.andreussi@unt.edu</a> </div>
<div style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;" class="">
Phone: +1-(940)-369-5316<br class="">
Skype: olivieroandreussi<br class="">
Web: <a href="https://nam04.safelinks.protection.outlook.com/?url=https%3A%2F%2Fwww.materialab.org%2F&data=04%7C01%7COliviero.Andreussi%40unt.edu%7C7bf7277f229c4050612508d9215248fb%7C70de199207c6480fa318a1afcba03983%7C0%7C1%7C637577459574757892%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&sdata=fK%2FXSMuliEpygyjMalUGn5XIwsgUGsFWTaRG2T7%2FcWE%3D&reserved=0" originalsrc="https://www.materialab.org/" shash="o5KhOUlCV3zbEwNvW4dcXNzQIm15LqkgOdxSSuYjLaD5liq6jIkteshtWMCw4gLdaAYdxuld4mVdCXJtbDOlNu3VsPK/WYm4a6r7oyCzqM5ZwTSmW/h1PJ7CUrPy/nfVT9LsP5ZLGyxG/MN1woEzhC/sx8GWdXQxbAl2qJzEq7Y=" class="" moz-do-not-send="true">https://www.materialab.org</a></div>
</div>
</div>
</div>
</div>
<div class=""><br class="">
<blockquote type="cite" class="">
<div class="">On May 27, 2021, at 3:04 PM, Paolo Giannozzi <<a href="mailto:paolo.giannozzi@uniud.it" class="" moz-do-not-send="true">paolo.giannozzi@uniud.it</a>> wrote:</div>
<br class="Apple-interchange-newline">
<div class="">
<div class="">Dear Laurent<br class="">
<br class="">
the proliferation of input variables in QE is an old and known problem. I am not sure what the best solution should be, though. One possibility I would consider is to read the needed variables from a separate file.<br class="">
If I remember correctly this is what the "Environ" library does.<br class="">
Another one is to extend the existing ATOMIC_FORCES card, that currently supports only a list of forces, to cover your case as well.<br class="">
<br class="">
Paolo<br class="">
<br class="">
On 27/05/2021 10:00, Laurent Pizzagalli wrote:<br class="">
<blockquote type="cite" class="">Dear Dr. Giannozzi,<br class="">
I write to you concerning the implementation of an external force field in quantum espresso. The motivation is the investigation of the mechanical properties of nanoparticles, especially in the plastic regime, obtained by uniaxial compression. The implementation
 is not complicated, since you (the developers) had already planned this potentiality in the file 'plugin_ext_forces.f90'.<br class="">
Actually, i already made the implementation in cp.x, but in a quick way, in a local copy of the code. Since it is successful, I would like to do it properly, merging it with the development version using git. I got almost all the information I need in the developer
 manual. But I still have one question, concerning the general organization of the input. Right now, I implemented only one kind of force field, fully repulsive. It allows for compressing the nanoparticle between two planes. Six keywords are needed, one for
 the activation of the force field, four to define the position of the planes and the increment for their displacements at each timestep, and a last one to define the strength of the repulsion. My question is: what would be the best way to organize these keywords?
 Originally, I put all of these into the &IONS namelist. Another option would be in the &SYSTEM namelist. And finally, maybe it would be better to create a specific NAMELIST to include these keywords (something like &EXTFF)? I also plan to implement another
 force field (including an attractive interaction), which would require other keywords. So the creation of a specific NAMELIST seems the best option in my opinion. But I would like to have yours before proceeding.<br class="">
Best regards<br class="">
Laurent<br class="">
</blockquote>
<br class="">
-- <br class="">
Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,<br class="">
Univ. Udine, via delle Scienze 208, 33100 Udine, Italy<br class="">
Phone +39-0432-558216, fax +39-0432-558222<br class="">
_______________________________________________<br class="">
developers mailing list<br class="">
<a href="mailto:developers@lists.quantum-espresso.org" class="" moz-do-not-send="true">developers@lists.quantum-espresso.org</a><br class="">
<a class="moz-txt-link-freetext" href="https://nam04.safelinks.protection.outlook.com/?url=https%3A%2F%2Flists.quantum-espresso.org%2Fmailman%2Flistinfo%2Fdevelopers&data=04%7C01%7COliviero.Andreussi%40unt.edu%7C7bf7277f229c4050612508d9215248fb%7C70de199207c6480fa318a1afcba03983%7C0%7C1%7C637577459574757892%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&sdata=xvjE2%2BkEdI7QaXQq3hqqszv1aia%2FtmS2uNzV3oA%2F57Q%3D&reserved=0" originalsrc="https://lists.quantum-espresso.org/mailman/listinfo/developers" shash="nMXop6bL6NWGEN15dR6cCaLWHBerInMvHfnSK3rEtg1UTclRSqgn9CZwnUqf5vsp3aeoNwQWKbZm5z2+a7a9lNE7ObLFFsNeiKqLhSqiqlkw7DHzsBlM11aoNjYH/XIQnHfycMGnTwlCT3kt/q8aSawhuN5Co0lTNd3DoEfpgLM=">https://nam04.safelinks.protection.outlook.com/?url=https%3A%2F%2Flists.quantum-espresso.org%2Fmailman%2Flistinfo%2Fdevelopers&amp;data=04%7C01%7Coliviero.andreussi%40unt.edu%7Cb9b4175eea884b881cce08d9214ac4da%7C70de199207c6480fa318a1afcba03983%7C0%7C0%7C637577427286515016%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C1000&amp;sdata=X49jcc%2FVhCsjeq%2BzOrROZJn3TLla6KVzog0P2JE9TLs%3D&amp;reserved=0</a><br class="">
</div>
</div>
</blockquote>
</div>
<br class="">
</div>
</div>
</blockquote>
<br class="">
<pre class="moz-signature" cols="72">-- 
                                                      ,,,     __,
                                                     /'^'\   |__|
                                                    ( o o )  |
--------------------------------------------------oOOO--(_)--OO|o------
<a class="moz-txt-link-rfc2396E" href="mailto:Laurent.Pizzagalli@univ-poitiers.fr"><Laurent.Pizzagalli@univ-poitiers.fr></a>
<a class="moz-txt-link-freetext" href="https://nam04.safelinks.protection.outlook.com/?url=http%3A%2F%2Flaurent.pizzagalli.free.fr%2F&data=04%7C01%7COliviero.Andreussi%40unt.edu%7C7bf7277f229c4050612508d9215248fb%7C70de199207c6480fa318a1afcba03983%7C0%7C1%7C637577459574767888%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&sdata=iUTq%2Bww%2Fvlr4ehxRe%2BTSwjzupaquB0orV%2FmU9h%2Brqgw%3D&reserved=0" originalsrc="http://laurent.pizzagalli.free.fr/" shash="BKQnCePTfRozRvrDgV3EDdzridWCSAh/V5PQSy25ly8qSIraTZG1Sq2Lvqgza3elccLUWdCu7sD7eyuifYPWTR0uzYPWA4Mo5IId0x/wchiyNyqky1sYqjWGoGUotQiCx//voPENGBMXVYzqdmgYbf6IPVWCshN2LmFlnaJgdPI=">http://laurent.pizzagalli.free.fr/</a>            Tel +33 549 49 74 99
------------------------------------------    Fax +33 549 49 66 92
Institut P'
Departement de Physique et de Mécanique des Matériaux
CNRS UPR 3346
Université de Poitiers
SP2MI
TSA 41123                                          .oooO
86073 Poitiers Cedex 9, FRANCE                     (   )   Oooo.
----------------------------------------------------\ (----(   )-------
                                                    \_)    ) /
                                                          (_/

</pre>
</div>
</div>
</blockquote>
</div>
<br class="">
</div>
</body>
</html>