diff --git a/nexus_infinity_core/qen.py b/nexus_infinity_core/qen.py
new file mode 100644
index 0000000..449d641
--- /dev/null
+++ b/nexus_infinity_core/qen.py
@@ -0,0 +1,146 @@
+import numpy as np
+from scipy.linalg import expm
+from qiskit import QuantumCircuit, execute
+from qiskit.providers.aer import AerSimulator
+
+class QEN:
+    """
+    Quantum Entanglement Network (QEN) for secure communication.
+
+    Attributes:
+    -----------
+    nodes : list
+        List of nodes in the QEN.
+    entanglements : dict
+        Dictionary of entanglements between nodes.
+    simulator : AerSimulator
+        Qiskit Aer simulator for simulating quantum circuits.
+    """
+
+    def __init__(self):
+        self.nodes = []
+        self.entanglements = {}
+        self.simulator = AerSimulator()
+
+    def add_node(self, node):
+        """
+        Add a node to the QEN.
+
+        Parameters:
+        -----------
+        node : str
+            Node identifier.
+        """
+        self.nodes.append(node)
+
+    def entangle_nodes(self, node1, node2):
+        """
+        Entangle two nodes in the QEN.
+
+        Parameters:
+        -----------
+        node1 : str
+            First node identifier.
+        node2 : str
+            Second node identifier.
+        """
+        # Create a quantum circuit for entanglement
+        qc = QuantumCircuit(2)
+        qc.h(0)
+        qc.cx(0, 1)
+
+        # Simulate the circuit
+        job = execute(qc, self.simulator)
+        result = job.result()
+        entangled_state = result.get_statevector()
+
+        # Store the entanglement
+        self.entanglements[(node1, node2)] = entangled_state
+
+    def measure_entanglement(self, node1, node2):
+        """
+        Measure the entanglement between two nodes.
+
+        Parameters:
+        -----------
+        node1 : str
+            First node identifier.
+        node2 : str
+            Second node identifier.
+
+        Returns:
+        -------
+        measurement : str
+            Measurement outcome (0 or 1).
+        """
+        entangled_state = self.entanglements[(node1, node2)]
+        measurement = np.random.choice([0, 1], p=[0.5, 0.5])
+        return measurement
+
+    def teleport_qubit(self, node1, node2, qubit):
+        """
+        Teleport a qubit from one node to another.
+
+        Parameters:
+        -----------
+        node1 : str
+            Source node identifier.
+        node2 : str
+            Destination node identifier.
+        qubit : str
+            Qubit to teleport (0 or 1).
+        """
+        # Create a quantum circuit for teleportation
+        qc = QuantumCircuit(3)
+        qc.x(qubit)
+        qc.cx(0, 1)
+        qc.cx(1, 2)
+        qc.measure(0, 1)
+
+        # Simulate the circuit
+        job = execute(qc, self.simulator)
+        result = job.result()
+        measurement = result.get_counts()
+
+        # Apply corrections based on measurement
+        if measurement[0] == '00':
+            correction = 'I'
+        elif measurement[0] == '01':
+            correction = 'X'
+        elif measurement[0] == '10':
+            correction = 'Z'
+        else:
+            correction = 'ZX'
+
+        # Apply correction to the destination node
+        if correction == 'X':
+            self.nodes[node2] = '1' if self.nodes[node2] == '0' else '0'
+        elif correction == 'Z':
+            self.nodes[node2] = '1' if self.nodes[node2] == '0' else '0'
+
+    def generate_shared_secret(self, node1, node2):
+        """
+        Generate a shared secret between two nodes.
+
+        Parameters:
+        -----------
+        node1 : str
+            First node identifier.
+        node2 : str
+            Second node identifier.
+
+        Returns:
+        -------
+        shared_secret : str
+            Shared secret key.
+        """
+        # Measure the entanglement
+        measurement1 = self.measure_entanglement(node1, node2)
+        measurement2 = self.measure_entanglement(node1, node2)
+
+        # Generate the shared secret
+        shared_secret = ''
+        for i in range(256):  # 256-bit shared secret
+            shared_secret += str(measurement1 ^ measurement2)
+
+        return shared_secret