<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" style="font-size:12pt;color:#000000;font-family:Calibri,Arial,Helvetica,sans-serif;" dir="ltr">
<p></p>
<div>
<div dir="ltr">
<div id="x_divtagdefaultwrapper" dir="ltr" style="font-size:12pt; color:rgb(0,0,0); font-family:Calibri,Arial,Helvetica,sans-serif,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols">
<p>Hi there,</p>
<p><br>
</p>
<p> I am trying to understand the ESM example with bc2 (metal-slab-metal) distributed with QE 6.0 which consists of an Al(001) slab between two electrodes, and at the same time I am trying to replicate the results of Otani and Sugino in their paper (Phys.
Rev. B, 73, 115407, 2006) which consists of an Al/Si(111)/Al slab between two screening boundaries.</p>
<p><br>
</p>
<p>I have a few questions:<br>
</p>
<p><br>
</p>
<p>** 1. ** In the QE example the following cell params and atomic positions are specified:</p>
<p><br>
</p>
<div>CELL_PARAMETERS bohr<br>
10.82227686 0.00000000 0.00000000<br>
0.00000000 10.82227686 0.00000000<br>
0.00000000 0.00000000 22.67672253<br>
ATOMIC_POSITIONS bohr<br>
Al 0.00000000 0.00000000 0.00000000<br>
Al 5.41113843 0.00000000 0.00000000 <br>
Al 0.00000000 5.41113843 0.00000000 <br>
Al 5.41113843 5.41113843 0.00000000 </div>
<p><br>
</p>
<p>My understanding is that the cell will repeat without overlapping atoms along X and Y. What I'm not clear on is why the 3rd lattice vector is so large?</p>
<p><br>
</p>
<p>The parameter 'esm_w' defaults to zero and is not specified in the example, so does this large value mean that there is a region of vacuum? Or does it just repeat the atoms such that in the Z-direction there are 5 Al layers? I'm guessing that this 3rd cell
parameter Z-value corresponds to L_z, which in the paper corresponds to z_0?</p>
<p><br>
</p>
<p>** 2. ** I have reproduced the 'surface unit cell' from Otani/Sugino's paper using VESTA.</p>
<p><br>
</p>
<p>My concern here is with the correct repetition of atoms in X and Y, and what exactly happens with the Z direction. I realize I must first rotate the unit cell such that the repeated Al atoms lie perpendicular to the Z Cartesian axis which I have not done
yet.</p>
<p><br>
</p>
<p>I have tested the repetition of this cell within VESTA along vectors perpendicular to the Al plane. So my most important question is: if the height of the cell along Z is ~14.4281 Angstrom (measured from Al center plane), then the third 'CELL_PARAMETERS'
entry (after rotation) should be this value exactly + the 'decay length' of the wavefunction in Z? Then the parameter 'esm_w' is used to set z_1 from the paper I suppose?<br>
</p>
<p><br>
</p>
<p>I will generate an xyz file of the surface cell, but if I repeat by 3.816 A along a surface lattice vector, it is a matter of just deleting the entries that cause an overlap of the atoms I suppose?</p>
<p><br>
</p>
<p>Thank you for your time and sorry for the very long email,</p>
<p>Kindest regards,</p>
<p>Louis<br>
</p>
</div>
</div>
</div>
<span class="PersonaPaneLauncher">
<div tabindex="-1" class="_pe_d _pe_12"></div>
</span><br>
<p></p>
</div>
</body>
</html>