To initialize the geometry

[tirthbha]

This is the one way to initialize the geometry.

[reshnashrestha]

Can you elaborate the code you have used like what does it do?

[pratibha77118]

@tirthbh how is this method different and effective from the one given in Com-Phy ? Could you explain how can your code be implemented for different structure?

[reshnashrestha]

Yes, exactly we can find the similar in compphys. But how can we be sure to choose for different configurations?

[tirthbha]

@reshnashrestha @pratibha77118 I have added details in the code shown below. This is just an example for the tetrahedron structure. You can take this as a reference for your chosen sructure.

import numpy as np

# Define the edge length of the tetrahedron based on expected ion separation
# Edge length is set to 0.3 nanometers, an arbitrary starting point that might represent the average ion distance
edge_length = 0.3  

# Define the vertices of a regular tetrahedron
# These vertices are in a normalized form and represent the corners of a regular tetrahedron
# Each vertex is a 3D coordinate [x, y, z]
vertices = np.array([[1, 1, 1],
                     [-1, -1, 1],
                     [-1, 1, -1],
                     [1, -1, -1]])

# Scale the vertices to the actual edge length of the tetrahedron
# Division by sqrt(3) scales the unit-length edges of the normalized tetrahedron to the desired edge length
tetrahedron_positions = vertices * (edge_length / np.sqrt(3))

# Assign Na+ and Cl- ions to the vertices alternately
# This loop creates pairs of ion type (Na or Cl) and their positions at the tetrahedron vertices
ideal_geometry = []
for i, pos in enumerate(tetrahedron_positions):
    ion = 'Na' if i % 2 == 0 else 'Cl'  # Alternating assignment: even index gets Na, odd index gets Cl
    ideal_geometry.append((ion, pos))   # Append a tuple of ion type and position

# Convert the list of ion positions and types into a structured NumPy array
# This array now represents the ideal geometric arrangement of Na+ and Cl- ions in a tetrahedral structure
ideal_geometry = np.array(ideal_geometry)