SAMBA-ilum 1.0.0.400__tar.gz → 1.0.0.402__tar.gz

{SAMBA_ilum-1.0.0.400 → samba_ilum-1.0.0.402}/PKG-INFO

@@ -1,12 +1,20 @@
 Metadata-Version: 2.1
 Name: SAMBA_ilum
-Version: 1.0.0.400
+Version: 1.0.0.402
 Summary: ...
 Author: Augusto de Lelis Araujo
 Author-email: augusto-lelis@outlook.com
 License: Closed source
 Description-Content-Type: text/markdown
 License-File: LICENSE.txt
+Requires-Dist: matplotlib
+Requires-Dist: pyfiglet
+Requires-Dist: requests
+Requires-Dist: plotly
+Requires-Dist: scipy
+Requires-Dist: numpy
+Requires-Dist: uuid
+Requires-Dist: vasprocar
 
 # SAMBA_ilum Copyright (C) 2024 - Closed source
 

{SAMBA_ilum-1.0.0.400 → samba_ilum-1.0.0.402}/SAMBA_ilum.egg-info/PKG-INFO

@@ -1,12 +1,20 @@
 Metadata-Version: 2.1
-Name: SAMBA-ilum
-Version: 1.0.0.400
+Name: SAMBA_ilum
+Version: 1.0.0.402
 Summary: ...
 Author: Augusto de Lelis Araujo
 Author-email: augusto-lelis@outlook.com
 License: Closed source
 Description-Content-Type: text/markdown
 License-File: LICENSE.txt
+Requires-Dist: matplotlib
+Requires-Dist: pyfiglet
+Requires-Dist: requests
+Requires-Dist: plotly
+Requires-Dist: scipy
+Requires-Dist: numpy
+Requires-Dist: uuid
+Requires-Dist: vasprocar
 
 # SAMBA_ilum Copyright (C) 2024 - Closed source
 

{SAMBA_ilum-1.0.0.400 → samba_ilum-1.0.0.402}/SAMBA_ilum.egg-info/SOURCES.txt

@@ -32,9 +32,8 @@ samba_ilum/src/potcar.py
 samba_ilum/src/vdW_DF.py
 samba_ilum/src/xy-scan.py
 samba_ilum/src/xy-scan_analysis.py
-samba_ilum/src/xyz-scan - Copia.py
 samba_ilum/src/xyz-scan.py
-samba_ilum/src/xyz-scan_analysis - Copia.py
+samba_ilum/src/xyz-scan_analysis.py
 samba_ilum/src/z-scan.py
 samba_ilum/src/z-scan_analysis.py
 samba_ilum/src/INPUTS/INCAR_a-scan

{SAMBA_ilum-1.0.0.400 → samba_ilum-1.0.0.402}/samba_ilum/__main__.py

@@ -16,7 +16,7 @@ dir_samba = os.path.dirname(os.path.realpath(__file__))
 print(f'{dir_samba}')
 #--------------------
 
-version = '1.0.0.400'
+version = '1.0.0.402'
 
 print(" ")
 print("=============================================================")

{SAMBA_ilum-1.0.0.400 → samba_ilum-1.0.0.402}/samba_ilum/src/energy_scan.py

@@ -20,11 +20,16 @@ energia = VTemp.replace('=',' ').split()
 energy.write(f'{name_dir} {energia[4]} \n')
 
 temp_name = name_dir.replace('_', ' ').split()
+t_temp_name = temp_name
+#------------------------
 if (len(temp_name) == 2):
    if (temp_name[0] == '0.0'): temp_name[0] = '1.0'
    if (temp_name[1] == '0.0'): temp_name[1] = '1.0'
    if (temp_name[0] == '1.0' or temp_name[1] == '1.0'):
       new_name_dir = str(temp_name[0]) + '_' + str(temp_name[1])
       energy.write(f'{new_name_dir} {energia[4]} \n')
+   if (t_temp_name[0] == '0.0' and t_temp_name[1] == '0.0'):
+      energy.write(f'1.0_0.0 {energia[4]} \n')
+      energy.write(f'0.0_1.0 {energia[4]} \n')
 
 energy.close()

{SAMBA_ilum-1.0.0.400 → samba_ilum-1.0.0.402}/samba_ilum/src/xy-scan_analysis.py

@@ -32,7 +32,18 @@ n_d = 301  # O dados do xy-scan serão interpolados para um grid de (n_d x n_d)
 # Obtendo a área no plano XY da Heteroestrutura =========
 #========================================================
 poscar = open('POSCAR.0', 'r')
-VTemp = poscar.readline().split()
+info = poscar.readline().split()
+#--------------------------------
+delta_A1_0 = 0.0
+delta_A2_0 = 0.0
+#---------------
+if (info[-7] == 'Shift_plane'):
+   delta_0 = info[-5].replace('_', ' ').split()
+   delta_A1_0 = float(delta_0[0])
+   delta_A2_0 = float(delta_0[1])
+   if (delta_A1_0 == 1.0): delta_A1_0 = 0.0
+   if (delta_A2_0 == 1.0): delta_A2_0 = 0.0
+#-----------------------------------------------
 VTemp = poscar.readline();  param = float(VTemp)
 #-----------------------------------------------
 A1 = poscar.readline().split();  A1x = float(A1[0])*param; A1y = float(A1[1])*param; A1z = float(A1[2])*param;  mA1 = np.linalg.norm(A1)
@@ -65,15 +76,27 @@ line     = np.argmin(date_E)
 delta    = date_shift[line]
 #------------------------------------------
 delta_min = delta.replace('_', ' ').split()
+#------------------------------------------
 a1_min = float(delta_min[0])
 a2_min = float(delta_min[1])
-x_min = (a1_min*A1x) + (a2_min*A2x)
-y_min = (a1_min*A1y) + (a2_min*A2y)
-#----------------------------------
-if (a1_min == -0.0): a1_min = 0.0
-if (a2_min == -0.0): a2_min = 0.0
+if (a1_min == -0.0 or a1_min == 1.0): a1_min = 0.0
+if (a2_min == -0.0 or a2_min == 1.0): a2_min = 0.0
+#-------------------------------------------------
+new_delta_A1 = delta_A1_0 + a1_min
+new_delta_A2 = delta_A2_0 + a2_min
+if (new_delta_A1 > 1.0): new_delta_A1 = new_delta_A1 -1.0
+if (new_delta_A2 > 1.0): new_delta_A2 = new_delta_A2 -1.0
+#-------------------------------------------------------------------------------------
+if (info[-7] == 'Shift_plane'): info[-5] = str(new_delta_A1) + '_' + str(new_delta_A2)
+new_line = ''
+for i in range(len(info)): new_line += info[i] + ' '
+new_line += '\n'
+#---------------------------------------------------
+x_min = (new_delta_A1*A1x) + (new_delta_A2*A2x)
+y_min = (new_delta_A1*A1y) + (new_delta_A2*A2y)
 if (x_min == -0.0):  x_min = 0.0
 if (y_min == -0.0):  y_min = 0.0
+#-------------------------------
 
 
 #-------------------------------------
@@ -152,9 +175,9 @@ T = np.linalg.inv(np.array([a, b, c]).T)  # Definindo a matriz de transformaçã
 #------------------------------------------------------------------------------
 r = np.array([delta_X, delta_Y, 0.0])        # Definindo o vetor posição cartesiano do átomo
 #------------------------------------
-f = np.dot(T, r)           # Calculando a correspondenre posição em coordenadas fracionárias    # ??????????????????????? Esta correto ???????????????????????
-for m in range(3):                                                                              # ?????????? Ao invés de 'f' não deveria ser 'f[m]' ??????????
-    f = np.where(f < 0, f + 1, f)                                                               # ??????????????????????? Esta correto ??????????????????????? 
+f = np.dot(T, r)           # Calculando a correspondenre posição em coordenadas fracionárias    # ??????????????????? Esta correto ???????????????????
+for m in range(3):                                                                              # ?????? Ao invés de 'f' não deveria ser 'f[m]' ??????
+    f = np.where(f < 0, f + 1, f)                                                               # ??????????????????? Esta correto ???????????????????
     f = np.where(f > 1, f - 1, f)
 #--------------------------------
 for m in range(3):
@@ -265,95 +288,20 @@ plt.savefig('xy-scan_direct.png', dpi = 600, bbox_inches='tight', pad_inches = 0
 """
 
 
-#==========================================================
-# Obtendo os vetores de rede A1 e A2 da Heteroestrutura ===
-#==========================================================
-poscar = open('POSCAR.0', "r")
-#-----------------------------
-VTemp = poscar.readline()
-VTemp = poscar.readline();  param = float(VTemp)
-VTemp = poscar.readline().split();  A1x = float(VTemp[0])*param;  A1y = float(VTemp[1])*param;  A1z = float(VTemp[2])*param
-VTemp = poscar.readline().split();  A2x = float(VTemp[0])*param;  A2y = float(VTemp[1])*param;  A2z = float(VTemp[2])*param
-VTemp = poscar.readline().split();  A3x = float(VTemp[0])*param;  A3y = float(VTemp[1])*param;  A3z = float(VTemp[2])*param
-#-------------
-poscar.close()
-#-------------
-
-
-#================================================
-# Gerando o arquivo POSCAR deslocado no plano ===
-#================================================
-poscar = open('CONTCAR.0', "r")
-poscar_new = open('CONTCAR', "w")
-#-------------------------------
-VTemp = poscar.readline()
-poscar_new.write(f'{VTemp}')
-VTemp = VTemp.split()
-nions1 = int(VTemp[2]);  nions2 = int(VTemp[3])
-#----------------------------------------------
-for k in range(7):
-    VTemp = poscar.readline()
-    poscar_new.write(f'{VTemp}')
-#-------------------------------
-for k in range(nions1):
-    VTemp = poscar.readline().split()
-    poscar_new.write(f'{VTemp[0]} {VTemp[1]} {VTemp[2]} \n')
-#-----------------------------------------------------------
-for k in range(nions2):
-    VTemp = poscar.readline().split()
-    #---------------------------------
-    temp_c1 = float(VTemp[0]) + a1_min
-    temp_c2 = float(VTemp[1]) + a2_min
-    #---------------------------------
-    # temp_c1 = float(VTemp[0]) + delta_A1
-    # temp_c2 = float(VTemp[1]) + delta_A2
-    #-------------------------------------
-    for i in range(3):
-        if (temp_c1 > 1.0): temp_c1 = temp_c1 -1.0
-        if (temp_c2 > 1.0): temp_c2 = temp_c2 -1.0
-    #---------------------------------------------
-    poscar_new.write(f'{temp_c1} {temp_c2} {VTemp[2]} \n')
-#---------------------------------------------------------
-poscar.close()
-poscar_new.close()
-#-----------------
-
-
-#===========================================================================
-# Gerando arquivo POSCAR com coordenadas cartesianas =======================
-#===========================================================================
-a = np.array([A1x, A1y, A1z])
-b = np.array([A2x, A2y, A2z])
-c = np.array([A3x, A3y, A3z])
-T = np.linalg.inv(np.array([a, b, c]).T)  # Definindo a matriz de transformação
-#------------------------------------------------------------------------------
-poscar = open('CONTCAR', "r")
-poscar_new = open('POSCAR', "w")
-#--------------------------------
-for k in range(7):
-    VTemp = poscar.readline()
-    poscar_new.write(f'{VTemp}')
-#------------------------
-VTemp = poscar.readline()
-poscar_new.write(f'Cartesian \n')
-
-#-------------------------------------------------------------------------------------------------------
-# Convertendo as posições atomicas diretas de todos os átomos da Supercélula para a forma cartesiana ---
-#-------------------------------------------------------------------------------------------------------
-for k in range(nions1 + nions2):
-    VTemp = poscar.readline().split()
-    #--------------------------------
-    k1 = float(VTemp[0]); k2 = float(VTemp[1]); k3 = float(VTemp[2])
-    #---------------------------------------------------------------
-    coord_x = ((k1*A1x) + (k2*A2x) + (k3*A3x))
-    coord_y = ((k1*A1y) + (k2*A2y) + (k3*A3y))
-    coord_z = ((k1*A1z) + (k2*A2z) + (k3*A3z))
-    poscar_new.write(f'{coord_x:>28,.21f} {coord_y:>28,.21f} {coord_z:>28,.21f} \n')
-#-------------
-poscar.close()
-poscar_new.close()
-#-----------------
-
+#------------------------------------------------
+shutil.copyfile(str(delta) + '/POSCAR', 'POSCAR')
+shutil.copyfile(str(delta) + '/CONTCAR', 'CONTCAR')
+#--------------------------------------------------
+if (info[-7] == 'Shift_plane'):
+   #---------------------------------------------------------
+   with open('CONTCAR', 'r') as file: line = file.readlines()
+   line[0] = new_line
+   with open('CONTCAR', 'w') as file: file.writelines(line)
+   #--------------------------------------------------------
+   with open('POSCAR', 'r') as file: line = file.readlines()
+   line[0] = new_line
+   with open('POSCAR', 'w') as file: file.writelines(line)
+   #------------------------------------------------------
 
 #=====================================================
 info = open('info_xy-scan.dat', "w", encoding='utf-8')
@@ -361,8 +309,8 @@ info.write(f'====================================================== \n')
 info.write(f'Displacement carried out over the 2nd material lattice   \n')
 #------------------------------------------------------------------
 info.write(f'Displacement_XY = ({x_min}, {y_min}) in Å \n')
-info.write(f'Displacement_XY = ({a1_min}*A1, {a2_min}*A2) \n')
-#----------------------------------------------------------------
+info.write(f'Displacement_XY = ({new_delta_A1}*A1, {new_delta_A2}*A2) \n')
+#-------------------------------------------------------------------------
 # info.write(f'Displacement_XY = ({delta_X}, {delta_Y}) in Å \n')
 # info.write(f'Displacement_XY = ({delta_A1}*A1, {delta_A2}*A2) \n')
 info.write(f'------------------------------------------------------ \n')
@@ -370,4 +318,3 @@ info.write(f'ΔE = {Delta_E_meV:.12f} meV/Å^2  or  {Delta_E_J:.12f} J/m^2 \n')
 info.write(f'====================================================== \n')
 info.close()
 #===========
-                                                                                                                                     

{SAMBA_ilum-1.0.0.400 → samba_ilum-1.0.0.402}/setup.py

@@ -6,7 +6,7 @@ import json
 
 setup(
     name = "SAMBA_ilum",
-    version = "1.0.0.400",
+    version = "1.0.0.402",
     entry_points={'console_scripts': ['samba_ilum = samba_ilum:main']},
     description = "...",
     author = "Augusto de Lelis Araujo", 

SAMBA_ilum-1.0.0.400/samba_ilum/src/xyz-scan - Copia.py

@@ -1,300 +0,0 @@
-# SAMBA_ilum Copyright (C) 2024 - Closed source
-
-
-import numpy as np
-import shutil
-import os
-
-
-# Atenção: ================================================================
-# O código foi escrito pensado em uma Heteroestrutura com n_Lattice = 2 ===
-# Para redes com n_Lattice > 2 testes e generalizações devem ser feitas === 
-#==========================================================================
-
-# Vacuo mínimo aplicado a Heteroestrutura
-vacuo = replace_vacuo
-
-# Separação no eixo-z a ser aplicada entre os materiais da Heteroestrutura
-z_scan = replace_zscan
-
-# Deslocamentos a ser aplicado no 2º material, referente ao vetor A1 da rede
-displacement_A1 = replace_displacement_xyz_A1
-
-# Deslocamentos a ser aplicado no 2º material, referente ao vetor A2 da rede 
-displacement_A2 = replace_displacement_xyz_A2
-
-
-"""
-#----------------------------------------------------------------
-# Testando a compatibilidade do arquivo POSCAR ------------------
-#----------------------------------------------------------------
-poscar = open('POSCAR', "r")
-VTemp = poscar.readline().split()
-poscar.close()
-#-------------
-crit = 0
-for k in range(len(VTemp)):
-    try:
-       inteiro = int(VTemp[k])
-       if (k > 0 and k < 3): crit += 1
-    except ValueError:
-       if (k == 0):  crit += 1
-    #------------------------------
-    if (len(VTemp) < 3 or crit < 3):
-       print(f' ')
-       print(f'========================================')
-       print(f'Verifique o arquivo POSCAR utilizado!   ')
-       print(f'INCOMPATIBILIDADE com o código detectada')
-       print(f'========================================')
-       print(f' ')
-       #==========
-       sys.exit()   
-       #=========
-"""
-
-
-#==========================================================================
-# Obtendo os vetores de rede da Heteroestrutura ===========================
-#==========================================================================
-contcar = open('CONTCAR', "r")
-#-----------------------------
-VTemp = contcar.readline()
-VTemp = contcar.readline();  param = float(VTemp)
-VTemp = contcar.readline().split();  A1x = float(VTemp[0]);  A1y = float(VTemp[1]);  A1z = float(VTemp[2])
-VTemp = contcar.readline().split();  A2x = float(VTemp[0]);  A2y = float(VTemp[1]);  A2z = float(VTemp[2])
-VTemp = contcar.readline().split();  A3x = float(VTemp[0]);  A3y = float(VTemp[1]);  A3z = float(VTemp[2])
-#--------------
-contcar.close()
-#--------------
-
-
-#==========================================================================
-# Obtendo a altura no eixo-z dos diferentes materiais =====================
-#==========================================================================
-contcar = open('CONTCAR', "r")
-#--------------------------------
-VTemp = contcar.readline().split()
-n_Lattice = len(VTemp[1].replace('+', ' ').split())
-nions_Lattice = []
-for m in range(n_Lattice):  nions_Lattice.append(int(VTemp[m+2]))  
-#----------------------------------------------------------------
-VTemp = contcar.readline();  param = float(VTemp)
-#----------------------------------------------------
-for k in range(3): VTemp = contcar.readline().split()
-fator_Z = float(VTemp[2])*param
-#--------------------------------------------
-for k in range(3): VTemp = contcar.readline()
-#--------------------------------------------------------------
-minZ = [0]*n_Lattice;  dZ = [0]*(n_Lattice +1);  dZ_total = 0.0
-#--------------------------------------------------------------
-for k in range(n_Lattice):
-    vZ = []
-    for m in range(nions_Lattice[k]):
-        VTemp = contcar.readline().split()
-        vZ.append(float(VTemp[2]))
-    #-----------------------------
-    dZ[k+1] = (max(vZ) - min(vZ))
-    dZ_total += dZ[k+1]*fator_Z
-    minZ[k] = min(vZ)
-#--------------------
-contcar.close()
-#--------------
-
-
-#==========================================================================
-# Deslocando os materiais para Z = 0.0 ====================================
-#==========================================================================
-contcar = open('CONTCAR', "r")
-poscar_new = open('POSCAR_temp', "w")
-#------------------------------------
-for k in range(8):
-    VTemp = contcar.readline()
-    poscar_new.write(f'{VTemp}')
-for k in range(n_Lattice):
-    for m in range(nions_Lattice[k]):
-        VTemp = contcar.readline().split()
-        temp_z = float(VTemp[2]) -minZ[k] +dZ[k]
-        if (temp_z < 0.0):  temp_z = 0.0
-        poscar_new.write(f'{float(VTemp[0])} {float(VTemp[1])} {temp_z} \n')
-#-------------
-contcar.close()
-poscar_new.close()
-#-----------------
-
-
-#===========================================================
-# Convertendo as coordenadas para a forma cartesiana =======
-#===========================================================
-poscar = open('POSCAR_temp', "r")
-poscar_new = open('POSCAR_cart', "w")
-VTemp = poscar.readline();  poscar_new.write(f'{VTemp}')
-VTemp = poscar.readline();  poscar_new.write(f'{VTemp}')
-VTemp = poscar.readline();  poscar_new.write(f'{VTemp}');  VTemp = VTemp.split();  A = [float(VTemp[0]), float(VTemp[1]), float(VTemp[2])]  
-VTemp = poscar.readline();  poscar_new.write(f'{VTemp}');  VTemp = VTemp.split();  B = [float(VTemp[0]), float(VTemp[1]), float(VTemp[2])]
-VTemp = poscar.readline();  poscar_new.write(f'{VTemp}');  VTemp = VTemp.split();  C = [float(VTemp[0]), float(VTemp[1]), float(VTemp[2])]
-VTemp = poscar.readline();  poscar_new.write(f'{VTemp}')
-VTemp = poscar.readline();  poscar_new.write(f'{VTemp}')
-VTemp = poscar.readline();  poscar_new.write(f'Cartesian \n')
-#-----------------------------------------------------------
-# Escrita das coordenadas cartesianas ----------------------
-#-----------------------------------------------------------
-for k in range(n_Lattice):
-    for m in range(nions_Lattice[k]):
-        VTemp = poscar.readline().split()
-        k1 = float(VTemp[0]); k2 = float(VTemp[1]); k3 = float(VTemp[2])
-        coord_x = ((k1*A[0]) + (k2*B[0]) + (k3*C[0]))
-        coord_y = ((k1*A[1]) + (k2*B[1]) + (k3*C[1]))
-        coord_z = ((k1*A[2]) + (k2*B[2]) + (k3*C[2]))
-        poscar_new.write(f'{coord_x:>28,.21f} {coord_y:>28,.21f} {coord_z:>28,.21f} \n')
-#-------------
-poscar.close()   
-poscar_new.close()
-#-----------------
-
-
-for z in range(len(z_scan)):
-    #-----------------
-    deltaZ = z_scan[z]
-    #-----------------
-
-    #===========================================================
-    # Gerando os arquivos POSCAR para cada separação em Z ======
-    #===========================================================
-    poscar = open('POSCAR_cart', "r")
-    poscar_new = open('POSCAR_deltaZ', "w")
-    VTemp = poscar.readline();  poscar_new.write(f'{VTemp}')
-    VTemp = poscar.readline();  poscar_new.write(f'{VTemp}')
-    VTemp = poscar.readline();  poscar_new.write(f'{VTemp}');  VTemp = VTemp.split();  A = [float(VTemp[0]), float(VTemp[1]), float(VTemp[2])]
-    VTemp = poscar.readline();  poscar_new.write(f'{VTemp}');  VTemp = VTemp.split();  B = [float(VTemp[0]), float(VTemp[1]), float(VTemp[2])]
-    VTemp = poscar.readline().split();  C = [float(VTemp[0]), float(VTemp[1]), float(VTemp[2])]
-    #---------------------------------------------
-    # temp_Z = (dZ_total + deltaZ_f + vacuo)/param
-    temp_Z = (dZ_total + deltaZ + vacuo)/param;  A3z = temp_Z
-    poscar_new.write(f'{float(VTemp[0]):>28,.21f} {float(VTemp[1]):>28,.21f} {float(temp_Z):>28,.21f} \n')
-    #-----------------------------------------------------------------------------------------------------
-    VTemp = poscar.readline();  poscar_new.write(f'{VTemp}')
-    VTemp = poscar.readline();  poscar_new.write(f'{VTemp}')
-    VTemp = poscar.readline();  poscar_new.write(f'Cartesian \n')
-    #------------------------------------------------------------
-    for k in range(n_Lattice):
-        for m in range(nions_Lattice[k]):
-            VTemp = poscar.readline().split()
-            coord_x = float(VTemp[0]); coord_y = float(VTemp[1]); coord_z = float(VTemp[2])
-            #-------------------------------------------------------------------------------
-            coord_z = coord_z + (vacuo/2)
-            if (k > 0):  coord_z = coord_z + deltaZ
-            poscar_new.write(f'{coord_x:>28,.21f} {coord_y:>28,.21f} {coord_z:>28,.21f} \n')
-            #-------------------------------------------------------------------------------
-    #-------------
-    poscar.close()   
-    poscar_new.close()
-    #-----------------
-
-
-    #==========================================================================
-    # Gerando e Armazenando os arquivos POSCAR deslocados no plano ============
-    #==========================================================================
-    os.mkdir('POSCAR_temp_cart')
-    #---------------------------
-    number = 0
-    #---------
-    for ii in displacement_A1:
-        for jj in displacement_A2:
-            #-----------------------------------------------
-            displacement_X = (ii*A1x*param) + (jj*A2x*param)        
-            displacement_Y = (ii*A1y*param) + (jj*A2y*param)   
-            #-----------------------------------------------
-            number += 1
-            #---------------------------------
-            poscar = open('POSCAR_deltaZ', "r")
-            poscar_new = open('POSCAR_temp_cart/POSCAR_' + str(number), "w") 
-            #---------------------------------------------------------------------------
-            VTemp = poscar.readline()
-            poscar_new.write(f'{VTemp}')
-            VTemp = VTemp.split()
-            nions1 = int(VTemp[2]);  nions2 = int(VTemp[3])
-            #----------------------------------------------
-            for k in range(7 + nions1):
-                VTemp = poscar.readline()
-                poscar_new.write(f'{VTemp}')
-            #-------------------------------
-            for k in range(nions2):
-                VTemp = poscar.readline().split()
-                poscar_new.write(f'{float(VTemp[0]) + displacement_X} {float(VTemp[1]) + displacement_Y} {VTemp[2]} \n')
-            #-----------------------------------------------------------------------------------------------------------
-            poscar.close()
-            poscar_new.close()
-            #-----------------
-
-
-    #=============================================================================
-    # Convertendo as coordenadas dos arquivos POSCAR de cartesiano para direto ===
-    #=============================================================================
-
-    #---------
-    number = 0
-    #----------
-    a = np.array([A1x*param, A1y*param, A1z*param])
-    b = np.array([A2x*param, A2y*param, A2z*param])
-    c = np.array([A3x*param, A3y*param, A3z*param])
-    T = np.linalg.inv(np.array([a, b, c]).T)  # Definindo a matriz de transformação
-    #------------------------------------------------------------------------------
-
-    for ii in displacement_A1:
-        for jj in displacement_A2: 
-            #----------
-            number += 1
-            #---------------------------------
-            dir_temp = str(ii) + '_' + str(jj) + '_' + str(deltaZ)
-            os.mkdir(dir_temp)
-            if os.path.isfile('vdw_kernel.bindat'): shutil.copyfile('vdw_kernel.bindat', dir_temp + '/vdw_kernel.bindat')
-            shutil.copyfile('contcar_update.py', dir_temp + '/contcar_update.py')
-            shutil.copyfile('energy_scan.py', dir_temp + '/energy_scan.py')
-            shutil.copyfile('KPOINTS', dir_temp + '/KPOINTS')
-            shutil.copyfile('POTCAR', dir_temp + '/POTCAR')
-            shutil.copyfile('INCAR', dir_temp + '/INCAR')
-            #-----------------------------------------------------------
-            poscar = open('POSCAR_temp_cart/POSCAR_' + str(number), "r")
-            poscar_new = open(dir_temp + '/POSCAR', "w") 
-            #-------------------------------------------
-            for k in range(7):
-                VTemp = poscar.readline()
-                poscar_new.write(f'{VTemp}')
-            #------------------------
-            VTemp = poscar.readline()
-            poscar_new.write(f'Direct \n')
-
-            #----------------------------------------------------------------------------------------------------
-            # Convertendo as posições atomicas cartesianas de todos os átomos da Supercélula para a forma direta,
-            # e ajustando as posições dos átomos que se encontram fora da célula.
-            #--------------------------------------------------------------------
-            for k in range(nions1 + nions2):
-                VTemp = poscar.readline().split()
-                x = float(VTemp[0])
-                y = float(VTemp[1])
-                z = float(VTemp[2])    
-                #----------------------
-                r = np.array([x, y, z])        # Definindo o vetor posição cartesiano do átomo  
-                #----------------------           
-                f = np.dot(T, r)               # Calculando a correspondenre posição em coordenadas fracionárias
-                for m in range(3):
-                    f = np.where(f < 0, f + 1, f)
-                    f = np.where(f > 1, f - 1, f)
-                #-------------------------------- 
-                for m in range(3):
-                    # f[m] = round(f[m], 6)
-                    if (f[m] > 0.9999 or f[m] < 0.0001):
-                       f[m] = 0.0
-                poscar_new.write(f'{f[0]} {f[1]} {f[2]} \n')
-
-
-            #-------------
-            poscar.close()
-            poscar_new.close()
-    #------------------------
-    os.remove('POSCAR_deltaZ')
-    shutil.rmtree('POSCAR_temp_cart')
-#------------------------------------
-os.remove('POSCAR_temp')
-os.remove('POSCAR_cart')
-#-----------------------