diff --git a/.gitignore b/.gitignore
index c96c3829d..e91ccac81 100644
--- a/.gitignore
+++ b/.gitignore
@@ -39,6 +39,8 @@ stack-local.yaml
 /result-*
 cardano-crypto-class/output
 cardano-crypto-praos/output
+cardano-crypto-tests/sig-pop-test-vector/target
+cardano-crypto-tests/sig-pop-test-vector/Cargo.lock
 
 # ghcid
 **/.ghcid
diff --git a/cardano-crypto-tests/sig-pop-test-vectors/Cargo.toml b/cardano-crypto-tests/sig-pop-test-vectors/Cargo.toml
new file mode 100644
index 000000000..aaf5ca0cd
--- /dev/null
+++ b/cardano-crypto-tests/sig-pop-test-vectors/Cargo.toml
@@ -0,0 +1,15 @@
+[package]
+name = "sig-pop-test-vectors"
+version = "0.1.0"
+edition = "2021"
+
+[dependencies]
+bls12_381 = { version = "0.8.0", features = ["default", "experimental"] }
+rand_chacha = "0.3.1"
+ff = "0.13.0"
+group = "0.13.0"
+hex = "0.4.3"
+blst = "0.3.10"
+sha2 = "0.9"
+mithril = "1.0.0"
+mithril-stm = "0.5.0"
diff --git a/cardano-crypto-tests/sig-pop-test-vectors/README.md b/cardano-crypto-tests/sig-pop-test-vectors/README.md
new file mode 100644
index 000000000..977fb81d8
--- /dev/null
+++ b/cardano-crypto-tests/sig-pop-test-vectors/README.md
@@ -0,0 +1,135 @@
+## Test vectors for signatures and Proof of possession
+
+TODO: Modify the README
+
+This is a rust script to generate test vectors for the following:
+- Using [bls12_381](https://github.com/zkcrypto/bls12_381)
+  - Pairing properties
+  - Elliptic curve operations
+  - Deserialization/decompression
+- Using [blst](https://github.com/supranational/blst) bindings from [bls12_381](https://github.com/zkcrypto/bls12_381) BLS signature with `aug` and `dst`.
+  The results are in hex encoding and stored under the folder `test_vectors`.
+
+### 1- Test vectors for pairing properties
+The properties to be tested:
+- `e([a]P, Q) = e(P, [a]Q)`
+- `e([a]P, [b]Q) = e([b]P, [a]Q)`
+- `e([a]P, [b]Q) = e([a * b]P, Q)`
+- `e([a]P, Q) * e([b]P, Q) = e([a + b]P, Q)`
+- `e([a]P, [b]Q) = e(P, [a * b]Q)`
+- `e(P, [a]Q) * e(P, [b]Q) = e(P, [a + b]Q)`
+
+The values used to generate test vectors:
+```
+a = 0x0e51216fa879b2ce727b596d065dd9b7fd8a84d94ffacf9ca30ad114304272d3 // scalar
+b = 0x437c2d7d852637c2ef23645a5abcbb308d6150bfcccbf3a8fdbc9daaa91496ef // scalar
+aplusb = 0x51cd4eed2d9fea91619ebdc7611a94e88aebd5991cc6c345a0c76ebed95709c2 // scalar
+atimesb = 0x2d70bbc706812d56e805ae67934b3275ff67f304a76ea9b3c96d31b9c0d607ba // scalar
+```
+
+Order of the values printed on `pairing_test_vectors`:
+- `P`
+- `[a]P`
+- `[b]P`
+- `[a + b]P`
+- `[a * b]P`
+- `Q`
+- `[a]Q`
+- `[b]Q`
+- `[a + b]Q`
+- `[a * b]Q`
+
+
+
+### 2- Test vectors for elliptic curve operations
+Operations to be tested:
+- Addition
+- Subtraction
+- Scalar multiplication
+- Negation
+
+The scalar used in scalar multiplication:
+```
+0x40df499974f62e2f268cd5096b0d952073900054122ffce0a27c9d96932891a5
+```
+
+Order of the values printed on `ec_operations_test_vectors`:
+
+- `G1_P` - random point on `G_1`
+- `G1_Q` - random point on `G_1`
+- `G1_ADD = G1_P + G1_Q`
+- `G1_SUB = G1_P - G1_Q`
+- `G1_MULL = [scalar]G1_Q`
+- `G1_NEG = -G1_P`
+- `G2_P` - random point on `G_2`
+- `G2_Q` - random point on `G_2`
+- `G2_ADD = G2_P + G2_Q`
+- `G2_SUB = G2_P - G2_Q`
+- `G2_MULL = [scalar]G2_Q`
+- `G2_NEG = -G2_P`
+
+
+### 3- Test vectors for deserialization/decompression
+- Point not in curve should fail deserialisation
+- Point not in group should fail deserialisation
+- Point not in curve should fail decompression
+- Point not in group should fail decompression
+
+Order of the values printed on `serde_test_vectors`:
+- `G1_uncomp_not_on_curve`
+- `G1_comp_not_on_curve`
+- `G1_comp_not_on_group`
+- `G1_uncomp_not_on_group`
+- `G2_uncomp_not_on_curve`
+- `G2_comp_not_on_curve`
+- `G2_comp_not_on_group`
+- `G2_uncomp_not_on_group`
+
+### 4- BLS Signature
+Test vectors for BLS signature, using `blst` bindings.
+
+The explicit usage of `aug` is not allowed in Cardano-base bindings. Therefore, before verification, they need to be appended into the message by following [hash-to-curve](https://datatracker.ietf.org/doc/html/draft-irtf-cfrg-hash-to-curve#name-expand_message) spec.
+
+`DST` and `msg` values used to generate test vectors:
+
+```rust
+let dst = b"BLS_SIG_BLS12381G2_XMD:SHA-256_SSWU_RO_NUL_";
+let msg = b"blst is such a blast";
+let aug = b"Random value for test aug";
+```
+Order of the values printed on the files `bls_sig_aug_test_vectors` (BLS signature with `aug`):
+
+- `sig`
+- `pk`
+
+To validate these test vectors, one needs to proceed as follows:
+```
+let hashed_msg = HashToG1Curve(aug || msg, dst);
+
+assert!(pairing(sig, G2Generator) ==  pairing(hashed_msg, pk)) 
+```
+
+### 5- Hash to curve with large DST
+The plutus bindings bound the DST to be at most 255 bytes, following the standard draft specification. If
+applications require a domain separation tag that is longer than 255 bytes, they should convert it to a smaller
+DST following the instructions of the standard draft (see [section 5.3.3](https://datatracker.ietf.org/doc/html/draft-irtf-cfrg-hash-to-curve#name-using-dsts-longer-than-255-)).
+
+We create test vectors to ensure that hashing to a curve by first hashing a large DST with SHA256, and then
+hashing to the curve works as expected. The test vectors of this file are stored in `h2c_large_dst` with the
+following order: 
+
+- 'msg'
+- 'large_dst'
+- Compressed G1 'output'
+
+To validate these test vectors, one needs to proceed as follows: 
+
+```
+let hashed_dst = Sha256(b"H2C-OVERSIZE-DST-" | large_dst);
+
+let hashed_output = HashToG1Curve(msg, hashed_dst);
+
+let expected_output = G1FromCompressed(output);
+
+assert!(expected_oputput == hashed_output); 
+```
diff --git a/cardano-crypto-tests/sig-pop-test-vectors/src/main.rs b/cardano-crypto-tests/sig-pop-test-vectors/src/main.rs
new file mode 100644
index 000000000..edbca8d4d
--- /dev/null
+++ b/cardano-crypto-tests/sig-pop-test-vectors/src/main.rs
@@ -0,0 +1,76 @@
+#![warn(missing_docs)]
+#![doc = include_str!("../README.md")]
+#![allow(non_snake_case)]
+
+use blst::{blst_p1, blst_scalar, blst_sk_to_pk_in_g1};
+use blst::min_sig::SecretKey as BlstSk;
+use rand_chacha::rand_core::{RngCore, SeedableRng};
+use rand_chacha::ChaCha20Rng;
+use std::fs::File;
+use std::io::prelude::*;
+
+
+
+fn generate_sk_and_pk<R: RngCore>(mut rng: R) -> std::io::Result<()>  {
+    let mut ikm = [0u8; 32];
+    rng.fill_bytes(&mut ikm);
+
+    let sk = BlstSk::key_gen(&ikm, &[])
+        .expect("Error occurs when the length of ikm < 32. This will not happen here.");
+
+    
+    let pk = sk.sk_to_pk();
+
+    let dst = b"BLS_SIG_BLS12381G2_XMD:SHA-256_SSWU_RO_NUL_";
+    let aug = b"Random value for test aug. ";
+    let msg = b"blst is such a blast";
+
+    let sig = sk.sign(msg, dst, aug);
+
+    // Should we keep the verify signature here?
+    let _verif = sig.verify(false, msg, dst, aug, &pk, false);
+
+    // Proof of Possession
+    const POP: &[u8] = &[80, 111, 80];
+    let k1 = sk.sign(POP, &[], &[]);
+    let k2 = unsafe {
+        let sk_scalar = std::mem::transmute::<&BlstSk, &blst_scalar>(&sk);
+        let mut out = blst_p1::default();
+        blst_sk_to_pk_in_g1(&mut out, sk_scalar);
+        out
+    };
+    let _pop = (k1, k2);
+
+    // Verify PoP
+    // let _ = pk.validate().unwrap();
+
+    // Write values to file to create test vectors
+    // I believe we need to save the following values:
+    // - ikm
+    // - sk
+    // - pk
+    // - sig
+    // - k1
+    // - k2
+    // - What else?
+
+
+    Ok(())
+
+}
+
+
+fn write_hex_to_file(file_name: &str, hex_strings: &[String]) -> std::io::Result<()> {
+    let mut file = File::create(file_name)?;
+
+    for string in hex_strings {
+        file.write_all(string.as_ref())?;
+        file.write_all(b"\n")?;
+    }
+    Ok(())
+}
+
+fn main() {
+    let mut rng = ChaCha20Rng::from_seed([0u8; 32]);
+    generate_sk_and_pk(&mut rng).expect("Failed to create large dst test vectors!");
+}