Simulating ESR echo transients

Topics related to Spinach package
Post Reply
zikriehb
Posts: 2
Joined: Tue Apr 16, 2024 8:46 pm

Simulating ESR echo transients

Post by zikriehb »

Hi,

I have a question about simulating ESR transient signals. If I simulate a Davies ENDOR experiment, is it possible for me to extract the net magnetization in the x-direction as a function of time? I'm particularly interested in seeing how the echo shape changes as I change the pulses. Additionally, I'm curious about the unwanted echoes to understand how phase-cycling affects the summed transient.

Part of me thinks that this is simple since I can just modify the davies_endor.m file such that I can extract trajectories from both shaped_pulse_af() and evolution() in the m file. However, I feel like I'm missing something else because of a passage in the examples/tutorials/time_domain_epr.pdf (page 5) that says 'it is not usually possible to predict what the echo is going to look like'.

Any insight is appreciated.
kuprov
Posts: 201
Joined: Mon Mar 29, 2021 4:26 pm

Re: Simulating ESR echo transients

Post by kuprov »

Yes, this is quite easy. Instead of the "final" option in the evolution calls inside endor_davies.m, use "trajectory" option with some suitable time step and an appropriate number of steps. Then put those trajectories into the "answer" variable so that the sequence returns them.

Optionally, in the main script, set parameters.sum_up=0 to tell the powder wrapper that it should not do the powder averaging, but return the trajectory for each individual orientation. Or maybe set the powder grid to 'single_crystal' so that the powder averaging operation des not get in the way.

Essentially, the dynamics you want is being computed anyway, you just need to tell Spinach to return it.
zikriehb
Posts: 2
Joined: Tue Apr 16, 2024 8:46 pm

Re: Simulating ESR echo transients

Post by zikriehb »

Thank you for the response Prof. Kuprov.

As a follow up, I tried to write a script to predict an echo of a nitroxide from a two-pulse Hahn echo. Attached to this post is the main script that uses the powder wrapper, a script for the Hahn echo pulse sequence, and a the signal as a function of time after the second pulse. I believe I am seeing the expected echo at the 400 ns mark.
Hahn_echo_test.png
Hahn_echo_test.png (25.45 KiB) Viewed 90762 times
However, it seems to be smaller than what I expected. I ensured that the pulses provided a correct pi/2 and pi pulse using the holeburn experiment (also included as a section in the main script).

So, I am currently struggling to understand why the echo is small (which I'm assuming is due to a mistake in my implementation of the trajectory).
Attachments
hahn_echo.m
(6.33 KiB) Downloaded 5225 times
Hahn_Echo_trajectory_no.m
(2.56 KiB) Downloaded 5253 times
kuprov
Posts: 201
Joined: Mon Mar 29, 2021 4:26 pm

Re: Simulating ESR echo transients

Post by kuprov »

I do get a better one using standard functionality in Spinach.
untitled.png
untitled.png (17.95 KiB) Viewed 90757 times
The enclosed echo test pulse sequence is a part of Spinach 2.9, the current developer version. It will be released within the next couple of weeks.
Attachments
endor_mims_echo_nitroxide.m
(1.11 KiB) Downloaded 5307 times
endor_mims_echo.m
(4.22 KiB) Downloaded 5165 times
Post Reply