<div dir="ltr">I did not read the structure, but it is not a trivial problem. You have to choose among several different approaches, and may test all of them:<div>1) try to draw a molecule with "chemical sense" (may using "molecular" calculations") and only later try to fit in the solid state</div><div>2) substitute an atom with that of your interest in an already optimized structure and hope...</div><div>3) try to find a suitable "space" for an extra interstitial atom and hope...</div><div><br></div><div>None is expected to be straightforward, and all require a certain degree of "chemical experience"...</div><div>HTH,</div><div> Carlo</div></div><div class="gmail_extra"><br><div class="gmail_quote">2015-02-25 3:18 GMT+01:00 Chenyang Lee <span dir="ltr"><<a href="mailto:lee19921027@gmail.com" target="_blank">lee19921027@gmail.com</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Dear experts in PWforum,<br>
<br>
I am experiencing difficulties relaxing some of the structures of B-HAp. Because the Boron group substitutions are proposed (unknown), thus the initial structure may deviate a lot from the final optimized one. Quantum Espresso seems not to be very capable to do these jobs (no offense) and the SCF cycles usually exceed 500 steps and it still cannot converge.<br>
I am attaching the input file, and I’d appreciate it if someone who is able to identify the problem can help me solve it.<br>
<br>
Thank you very much,<br>
Chenyang<br>
__<br>
Chenyang Li<br>
Undergraduate Student, Materials Science and Engineering<br>
The Pennsylvania State University, University Park<br>
Phone: <a href="tel:%28814%29441-4312" value="+18144414312">(814)441-4312</a><br>
<br>
_______________________________________________<br>
<br>
&control<br>
calculation = 'vc-relax'<br>
restart_mode = 'from_scratch'<br>
prefix = ‘BHAp'<br>
pseudo_dir = '/gpfs/home/col5422/work/HAp/pseudo/'<br>
outdir = './scratch/'<br>
tprnfor = .true.<br>
forc_conv_thr = 1e-4<br>
etot_conv_thr = 1e-6<br>
nstep = 500<br>
/<br>
&system<br>
ibrav = 14<br>
A = 9.4166<br>
B = 9.4166<br>
C = 6.8745<br>
cosAC = 0.0<br>
cosAB = -0.5<br>
cosBC = 0.0<br>
nat = 52<br>
ntyp = 6<br>
ecutwfc = 40<br>
ecutrho = 400<br>
nosym = .true.<br>
occupations = 'smearing'<br>
smearing = 'gaussian'<br>
degauss = 0.02<br>
/<br>
&electrons<br>
electron_maxstep = 500<br>
diagonalization = 'david'<br>
mixing_mode = 'plain'<br>
mixing_beta = 0.3<br>
conv_thr = 1e-6<br>
/<br>
&ions<br>
ion_dynamics = ‘damp'<br>
/<br>
&cell<br>
cell_dynamics = ‘damp-pr'<br>
/<br>
<br>
ATOMIC_SPECIES<br>
H 1 H.pbe-kjpaw_psl.0.1.upf<br>
B 1 B.pbe-n-kjpaw_psl.0.1.upf<br>
C 1 C.pbe-n-kjpaw_psl.0.1.upf<br>
O 1 O.pbe-n-kjpaw_psl.0.1.upf<br>
P 1 P.pbe-n-kjpaw_psl.0.1.upf<br>
Ca 1 Ca.pbe-spn-kjpaw_psl.0.2.3.upf<br>
<br>
ATOMIC_POSITIONS crystal<br>
H 0.000000 0.000000 0.400000<br>
H 0.000000 0.000000 0.900000<br>
H 0.567369 0.551274 0.972591<br>
H 0.629383 0.851436 0.797858<br>
H 0.357078 0.434312 0.660517<br>
H 0.809237 0.663097 0.589841<br>
B 0.581602 0.618142 0.689068<br>
C 0.310984 0.129561 0.635342<br>
O 0.348000 0.259000 0.073000<br>
O 0.741000 0.089000 0.073000<br>
O 0.911000 0.652000 0.073000<br>
O 0.259000 0.911000 0.573000<br>
O 0.089000 0.348000 0.573000<br>
O 0.259000 0.911000 0.927000<br>
O 0.089000 0.348000 0.927000<br>
O 0.348000 0.259000 0.427000<br>
O 0.741000 0.089000 0.427000<br>
O 0.911000 0.652000 0.427000<br>
O 0.329000 0.484000 0.250000<br>
O 0.589000 0.466000 0.250000<br>
O 0.516000 -0.155000 0.250000<br>
O 0.534000 0.123000 0.250000<br>
O 1.155000 0.671000 0.250000<br>
O 0.877000 0.411000 0.250000<br>
O 0.484000 1.155000 0.750000<br>
O 0.466000 0.877000 0.750000<br>
O -0.155000 0.329000 0.750000<br>
O 0.123000 0.589000 0.750000<br>
O 0.000000 0.000000 0.250000<br>
O 0.000000 0.000000 0.750000<br>
O 0.603664 0.524999 0.843110<br>
O 0.718661 0.681136 0.557507<br>
O 0.546871 0.739153 0.774413<br>
O 0.441199 0.509454 0.566674<br>
O 0.464549 0.231819 0.634250<br>
O 0.235429 0.067896 0.472588<br>
O 0.237756 0.091652 0.804032<br>
P 0.400000 0.369000 0.250000<br>
P 0.631000 0.031000 0.250000<br>
P 0.969000 0.600000 0.250000<br>
P 0.369000 0.969000 0.750000<br>
P 0.031000 0.400000 0.750000<br>
Ca 0.333333 0.666667 0.001000<br>
Ca 0.666667 0.333333 0.501000<br>
Ca 0.666667 0.333333 0.999000<br>
Ca 0.333333 0.666667 0.499000<br>
Ca 0.246000 0.993000 0.250000<br>
Ca 0.007000 0.253000 0.250000<br>
Ca 0.747000 0.754000 0.250000<br>
Ca 0.754000 0.007000 0.750000<br>
Ca 0.993000 0.747000 0.750000<br>
Ca 0.253000 0.246000 0.750000<br>
<br>
K_POINTS automatic<br>
4 4 4 1 1 1<br>
_____________________________________________<br>
<br>
<br>
<br>
_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://pwscf.org/mailman/listinfo/pw_forum" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature"><div dir="ltr"><pre>------------------------------------------------------------
Prof. Carlo Nervi <a href="mailto:carlo.nervi@unito.it" target="_blank">carlo.nervi@unito.it</a> <a>Tel:+39</a> 0116707507/8
Fax: +39 0116707855 - Dipartimento di Chimica, via
P. Giuria 7, 10125 Torino, Italy. <a href="http://lem.ch.unito.it/" target="_blank">http://lem.ch.unito.it/</a><br></pre></div></div>
</div>