<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>
</head>
<body dir="ltr">
<div id="divtagdefaultwrapper" dir="ltr" style="font-size:12pt; color:rgb(0,0,0); font-family:Calibri,Helvetica,sans-serif,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols">
<p></p>
<div>
<div style="margin-top:0; margin-bottom:0">Dear Quantum Espresso developers,</div>
<div style="margin-top:0; margin-bottom:0"><br>
</div>
<div style="margin-top:0; margin-bottom:0">as this is my first message posted here, let me first introduce myself: my name is Giovanni Marini, I am currently working as a postdoc at Istituto Italiano di Tecnologia, collaborating with Prof. Matteo Calandra.
It is a great pleasure for me to become part of this community.</div>
<div style="margin-top:0; margin-bottom:0"><br>
</div>
<div style="margin-top:0; margin-bottom:0">I am writing this post because we have been working on a new methodology to treat photoexcited insulating system in the ultrafast regime ( see Ref.
<a href="https://journals.aps.org/prb/abstract/10.1103/PhysRevB.104.144103" target="_blank" rel="noopener noreferrer" class="x_OWAAutoLink" id="LPlnk820480" previewremoved="true">
https://journals.aps.org/prb/abstract/10.1103/PhysRevB.104.144103</a> ), and we have been developing our implementation inside Quantum Espresso.
<br>
</div>
<div style="margin-top:0; margin-bottom:0"><br>
</div>
<div style="margin-top:0; margin-bottom:0">At present, we would like to publicly release our implementation, possibly within Quantum Espresso. Let me give some further details on what we have been doing.
<br>
</div>
<div style="margin-top:0; margin-bottom:0"><br>
</div>
<div style="margin-top:0; margin-bottom:0">Our work is divided in two parts: <br>
</div>
<div style="margin-top:0; margin-bottom:0"><br>
</div>
<div style="margin-top:0; margin-bottom:0">The first part operates inside pw.x and allows to perform self-consistent, non-self consistent, total energy, forces and stress tensor calculations in the presence of two chemical potentials, one for the holes and
one for the electrons, constraining some of the electrons in a subset of bands (conduction bands). Occupations in valence and conduction must be treated using smearing.
<br>
</div>
<div style="margin-top:0; margin-bottom:0">The second part operates inside ph.x (and hp.x) and is intended to allow linear response calculations and electron-phonon matrix elements calculations within the same framework.
<br>
</div>
<div style="margin-top:0; margin-bottom:0"><br>
</div>
<div style="margin-top:0; margin-bottom:0">As a first step, we would like to publish only the part relative to pw.x, and leave the second part (somewhat more complicated) as a future follow-up.
<br>
</div>
<div style="margin-top:0; margin-bottom:0">The current status of the project can be found in my fork of the develop branch,
<a href="https://gitlab.com/giovannimarini1993/q-e_two_chem_pot" class="OWAAutoLink" id="LPlnk495547" previewremoved="true">
https://gitlab.com/giovannimarini1993/q-e_two_chem_pot</a> . I did my best to follow the developers guidelines. I still didn't open any merge request.
<br>
</div>
<div dir="ltr" style="font-size:12pt; color:rgb(0,0,0); font-family:Calibri,Helvetica,sans-serif,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols">
<div style="margin-top:0; margin-bottom:0"><br>
</div>
<div style="margin-top:0; margin-bottom:0">You can find a brief description of the operated changes in the .txt file at the following link:
<a href="https://file.io/cO1lFd8pYiAW" class="OWAAutoLink" id="LPlnk190888" previewremoved="true">
https://file.io/cO1lFd8pYiAW</a><br>
</div>
<div style="margin-top:0; margin-bottom:0"><br>
</div>
Changes to the pw.x input include the addition of four new variables, one in the &control namelist, twochem, regulating the activation of a two-chemical potential calculation, and three other variables in &system. It would be useful for the future to add these
variables could be added to the .xml output, as they would be needed in the linear response calculation.
<br>
</div>
<br>
</div>
<div>The heart of this new implementation can be found in the new module file "PW/src/two_chem.f90", containing some new specific routines to treat the new case. The normal Quantum Espresso workflow would be modified in <span>"PW/src/weights.f90", where the
new subroutines are invoked when twochem=.true. (nothing changes when twochem=.false.).<br>
</span></div>
<div><span></span><span>Temporarily, I inserted a stop inside the PHonon code when it detects that pw.x ran with twochem=.true. (I let pw produce a certain empty file in this case.)</span></div>
<div><br>
</div>
<div>Let me also add that both magnetism (excluding fixed tot_magnetization case) and non-collinearity are implemented.
<br>
</div>
<div>Finally, I also put a new test folder (pw_twochem) in the test-suite folder, according to the procedure specified on gitlab, and I modified the documentation in PW/Doc/INPUT_PW.def to explain the new introduced variables.</div>
<div><br>
</div>
<div>Let me know what you think of this project and if you deem possible to integrate it inside Quantum Espresso.<br>
<div dir="ltr" style="font-size:12pt; color:rgb(0,0,0); font-family:Calibri,Helvetica,sans-serif,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols">
<br>
Thank you very much and</div>
<div dir="ltr" style="font-size:12pt; color:rgb(0,0,0); font-family:Calibri,Helvetica,sans-serif,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols">
Best Regards,</div>
<div dir="ltr" style="font-size:12pt; color:rgb(0,0,0); font-family:Calibri,Helvetica,sans-serif,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols">
<br>
</div>
Giovanni Marini</div>
<div>Postdoctoral Researcher</div>
<div>Istituto Italiano di Tecnologia</div>
<div><br>
</div>
<div>email: giovanni.marini@iit.it</div>
<div>phone number: +39 3888284164<br>
</div>
<br>
<p></p>
</div>
</body>
</html>