DRCMR Logo 300px Color
 

Quick start guide to the Bloch Simulator

This page provides step-by-step examples to practical use of the Bloch Simulator made for visualizing basic and advanced Magnetic Resonance techniques. The examples are spin-echo formation and the excitation process, which are also demonstrated in corresponding YouTube videos. The basic MR knowledge needed to understand the text below is described in an MRI tutorial targeted at a broad audience.

  • It is easiest to run the simulator directly in an internet browser. If you have not already started it, please do so by following this link starting the software in a new browser window: http://www.drcmr.dk/BlochSimulator
  • If you have Flash installed on your computer, you should now see the magnetization vector of an imaginary sample precessing in a magnetic field, B0, pointing upwards. You can always get back to this situation by choosing "Scene: Precession" in the pull-down menu on the lower left.
  • First try to adjust the B0 field strength by moving the slider named "B0". Notice how the precession frequency changes with the field. Think of the units as tesla, for example.

Spin-echo demonstration, step by step:

The spin echo sequence
The spin echo sequence
  • Now, let's introduce field inhomogeneity so that all nuclei do not precess at the same frequency. This can be done by choosing "Scene: Weak inhomogeneity" in the pull-down menu on the lower left. You now see the magnetization in thermal equilibrium. The floor rotates since you see this in the rotating frame of resonance which is only needed later. Therefore, shift back to the stationary frame by pressing the  "Change frame"-button.
  • Since weak inhomogeneity was selected above, the nuclei precess at different frequencies, but it does not show until the magnetization is rotated into the transversal plane. Do this by pressing the "90x hard" button (initiates 90 degree RF pulse rotating the magnetization quickly around the x-axis in the transversal plane). You should now see different components of the magnetization precessing at slightly different frequencies. Gradual dephasing leads to a loss of signal seen in the lower right graph. Press the "v" key on keyboard a few times - it toggles the viewing angle.
  • You may want to repeat the last few steps by pressing the "Scene: Weak inhomogeneity" button again and doing a new excitation. Proceed when you understand what you are seeing so far. Let's add a short 180 degree refocusing pulse to our sequence after some dephasing has occurred (press the "180y hard"-button). This will lead to the formation of a spin-echo after a period equal to the dephasing period. Notice how the nuclei get back in phase and how the lost signal is recovered.
  • This finishes the spin-echo demonstration, but you can continue experimenting with it, e.g. by playing out extra refocusing pulses or by adding a bit of relaxation using the T2-slider for example (choose a large value for T2 initially, e.g. 20 seconds). You can also watch the events in the rotating frame of resonance.

Off- and on-resonance excitation, step by step:

  • Go back to the initial situation by  choosing "Scene: Precession" in the pull-down menu on the lower left.
  • Adjust the B0 field strength by moving the slider named "B0". Notice how the precession frequency changes with the field. Set the value of B0 back to 3 by typing the value in the text field above the slider, followed by "enter". Think of the units as tesla, for example.
  • Radio waves applied off-resonance
    Radio waves applied off-resonance
    Now turn on a radio wave field. Do this by setting the RF amplitude to 0.1 using the appropriate slider or by typing the value in the corresponding text field, followed by "enter". The units are the same as for the B0 field. The red bar shows the push of the radio waves on the magnetization (the torque). Notice how the push is orthogonal to the magnetization, and how it rotates.
  • The push does not seem to have much of an effect on the magnetization. This is because the frequency of the radio waves is not matched to the Larmor frequency (the radiowaves are off-resonance). The reason for the missing effect may be easier to see in the rotating frame of reference. Press "Change frame" to toggle between viewing in the stationary and the rotating frames of reference as you wish. Notice how the magnetization is not pushed consistently by off-resonance radiowaves (in synchony with the precession). Therefore it just wiggles slightly in the radio wave field. Switch back to the stationary frame of reference.
  • Adjust the radio frequency to match the Larmor frequency simply by setting the RF frequency equal to B0 (3 in this case). This is valid since the gyromagnetic ratio in the simulator is set to 1 Hz/T (approximately). Hence the resonance frequency for a field of value 3 is simply 3. Notice how the magnetization is now gradually changed by the RF field. Watching this in the rotating frame makes it appear much simpler -- switch back and forth a few times, and notice how resonance radiowaves induces rotations in the rotating frame of reference.
  • Radio waves applied on resonance
    Radio waves applied on resonance
    Now you have seen how resonant radiowaves change the magnetization, and how off-resonant waves does not, lets start from equilibrium and do an excitation. Choose "Scene: Equilibrium". Add an RF field by setting the RF amplitude equal to 0.1. Notice how little happens since the RF frequency is not adjusted to the Larmor frequency. Adjust it as described above and see how the magnetization is rotated away from equilibirum, e.g. so it after a while points in the exact opposite directions. Similarly, on-resonance rotations away from equilibrium can be triggered by pressing the buttons labelled "90x hard", "90x selective", "180y hard" and "30x hard" that send bursts of radiowaves sufficiently long to cause the specified rotations.


The simulator can do much more. A few examples are given in YouTube videos and on the project homepage but you really need to experiment yourself to get full benefit. The "Challenges" section in the "Help" menu may also be of inspiration.

 

This page provides step-by-step examples to practical use of the Bloch Simulator made for visualizing basic and advanced Magnetic Resonance techniques. The examples are spin-echo formation and the excitation process, which are also demonstrated in corresponding YouTube videos. The basic MR knowledge needed to understand the text below is described in an MRI tutorial targeted at a broad audience.

  • It is easiest to run the simulator directly in an internet browser. If you have not already started it, please do so by following this link starting the software in a new browser window: http://www.drcmr.dk/BlochSimulator
  • If you have Flash installed on your computer, you should now see the magnetization vector of an imaginary sample precessing in a magnetic field, B0, pointing upwards. You can always get back to this situation by choosing "Scene: Precession" in the pull-down menu on the lower left.
  • First try to adjust the B0 field strength by moving the slider named "B0". Notice how the precession frequency changes with the field. Think of the units as tesla, for example.

Spin-echo demonstration, step by step:

The spin echo sequence
The spin echo sequence
  • Now, let's introduce field inhomogeneity so that all nuclei do not precess at the same frequency. This can be done by choosing "Scene: Weak inhomogeneity" in the pull-down menu on the lower left. You now see the magnetization in thermal equilibrium. The floor rotates since you see this in the rotating frame of resonance which is only needed later. Therefore, shift back to the stationary frame by pressing the  "Change frame"-button.
  • Since weak inhomogeneity was selected above, the nuclei precess at different frequencies, but it does not show until the magnetization is rotated into the transversal plane. Do this by pressing the "90x hard" button (initiates 90 degree RF pulse rotating the magnetization quickly around the x-axis in the transversal plane). You should now see different components of the magnetization precessing at slightly different frequencies. Gradual dephasing leads to a loss of signal seen in the lower right graph. Press the "v" key on keyboard a few times - it toggles the viewing angle.
  • You may want to repeat the last few steps by pressing the "Scene: Weak inhomogeneity" button again and doing a new excitation. Proceed when you understand what you are seeing so far. Let's add a short 180 degree refocusing pulse to our sequence after some dephasing has occurred (press the "180y hard"-button). This will lead to the formation of a spin-echo after a period equal to the dephasing period. Notice how the nuclei get back in phase and how the lost signal is recovered.
  • This finishes the spin-echo demonstration, but you can continue experimenting with it, e.g. by playing out extra refocusing pulses or by adding a bit of relaxation using the T2-slider for example (choose a large value for T2 initially, e.g. 20 seconds). You can also watch the events in the rotating frame of resonance.

Off- and on-resonance excitation, step by step:

  • Go back to the initial situation by  choosing "Scene: Precession" in the pull-down menu on the lower left.
  • Adjust the B0 field strength by moving the slider named "B0". Notice how the precession frequency changes with the field. Set the value of B0 back to 3 by typing the value in the text field above the slider, followed by "enter". Think of the units as tesla, for example.
  • Radio waves applied off-resonance
    Radio waves applied off-resonance
    Now turn on a radio wave field. Do this by setting the RF amplitude to 0.1 using the appropriate slider or by typing the value in the corresponding text field, followed by "enter". The units are the same as for the B0 field. The red bar shows the push of the radio waves on the magnetization (the torque). Notice how the push is orthogonal to the magnetization, and how it rotates.
  • The push does not seem to have much of an effect on the magnetization. This is because the frequency of the radio waves is not matched to the Larmor frequency (the radiowaves are off-resonance). The reason for the missing effect may be easier to see in the rotating frame of reference. Press "Change frame" to toggle between viewing in the stationary and the rotating frames of reference as you wish. Notice how the magnetization is not pushed consistently by off-resonance radiowaves (in synchony with the precession). Therefore it just wiggles slightly in the radio wave field. Switch back to the stationary frame of reference.
  • Adjust the radio frequency to match the Larmor frequency simply by setting the RF frequency equal to B0 (3 in this case). This is valid since the gyromagnetic ratio in the simulator is set to 1 Hz/T (approximately). Hence the resonance frequency for a field of value 3 is simply 3. Notice how the magnetization is now gradually changed by the RF field. Watching this in the rotating frame makes it appear much simpler -- switch back and forth a few times, and notice how resonance radiowaves induces rotations in the rotating frame of reference.
  • Radio waves applied on resonance
    Radio waves applied on resonance
    Now you have seen how resonant radiowaves change the magnetization, and how off-resonant waves does not, lets start from equilibrium and do an excitation. Choose "Scene: Equilibrium". Add an RF field by setting the RF amplitude equal to 0.1. Notice how little happens since the RF frequency is not adjusted to the Larmor frequency. Adjust it as described above and see how the magnetization is rotated away from equilibirum, e.g. so it after a while points in the exact opposite directions. Similarly, on-resonance rotations away from equilibrium can be triggered by pressing the buttons labelled "90x hard", "90x selective", "180y hard" and "30x hard" that send bursts of radiowaves sufficiently long to cause the specified rotations.


The simulator can do much more. A few examples are given in YouTube videos and on the project homepage but you really need to experiment yourself to get full benefit. The "Challenges" section in the "Help" menu may also be of inspiration.

 

Selected Publications

Webb JL, Troise L, Hansen NW, Olsson C, Wojciechowski AM, Achard J, Brinza O, Staacke R, Kieschnick M, Meijer J, Thielscher A, Perrier J-F, Berg-Sørensen K, Huck A, Andersen UL. 2021. Detection of biological signals from a live mammalian muscle using an early stage diamond quantum sensor. Scientific Reports. 11(1):1-11. https://doi.org/10.1038/s41598-021-81828-x

von Conta J, Kasten FH, Ćurčić-Blake B, Aleman A, Thielscher A, Herrmann CS. 2021. Interindividual variability of electric fields during transcranial temporal interference stimulation (tTIS). Scientific Reports. 11(1):1-12. https://doi.org/10.1038/s41598-021-99749-0

Splittgerber M, Borzikowsky C, Salvador R, Puonti O, Papadimitriou K, Merschformann C, Biagi MC, Stenner T, Brauer H, Breitling-Ziegler C, Prehn-Kristensen A, Krauel K, Ruffini G, Pedersen A, Nees F, Thielscher A, Dempfle A, Siniatchkin M, Moliadze V. 2021. Multichannel anodal tDCS over the left dorsolateral prefrontal cortex in a paediatric population. Scientific Reports. 11(1):1-15. https://doi.org/10.1038/s41598-021-00933-z

Shirinpour S, Mantell K, Li X, Puonti O, Madsen K, Haigh Z, Casillo EC, Alekseichuk I, Hendrickson T, Xu T. 2021. New tools for computational modeling of non-invasive brain stimulation in SimNIBS. Brain Stimulation: Basic, Translational, and Clinical Research in Neuromodulation. 14(6):1644. https://doi.org/10.1016/j.brs.2021.10.180

Saturnino GB, Madsen KH, Thielscher A. 2021. Optimizing the electric field strength in multiple targets for multichannel transcranial electric stimulation. Journal of Neural Engineering. 18(1): Article 014001. https://doi.org/10.1088/1741-2552/abca15

Numssen O, Zier A-L, Thielscher A, Hartwigsen G, Knösche TR, Weise K. 2021. Efficient high-resolution TMS mapping of the human motor cortex by nonlinear regression. NeuroImage. 245:1-11. https://doi.org/10.1016/j.neuroimage.2021.118654

Montanaro H, Pasquinelli C, Lee HJ, Kim H, Siebner HR, Kuster N, Thielscher A, Neufeld E. 2021. The impact of CT image parameters and skull heterogeneity modeling on the accuracy of transcranial focused ultrasound simulations. Journal of Neural Engineering. 18(4):1-28. https://doi.org/10.1088/1741-2552/abf68d

Mezger E, Rauchmann B-S, Brunoni AR, Bulubas L, Thielscher A, Werle J, Mortazavi M, Karali T, Stöcklein S, Ertl-Wagner B, Goerigk S, Padberg F, Keeser D. 2021. Effects of bifrontal transcranial direct current stimulation on brain glutamate levels and resting state connectivity: multimodal MRI data for the cathodal stimulation site. European Archives of Psychiatry and Clinical Neuroscience. 271(1):111-122. https://doi.org/10.1007/s00406-020-01177-0

Karadas M, Olsson C, Winther Hansen N, Perrier J-F, Webb JL, Huck A, Andersen UL, Thielscher A. 2021. In-vitro Recordings of Neural Magnetic Activity From the Auditory Brainstem Using Color Centers in Diamond: A Simulation Study. Frontiers in Neuroscience. 15:1-17. https://doi.org/10.3389/fnins.2021.643614

Gregersen F, Göksu C, Schaefers G, Xue R, Thielscher A, Hanson LG. 2021. Safety Evaluation of a New Setup for Transcranial Electric Stimulation during Magnetic Resonance Imaging. Brain Stimulation. 14(3):488-497. https://doi.org/10.1016/j.brs.2021.02.019

Göksu C, Scheffler K, Gregersen F, Eroğlu HH, Heule R, Siebner HR, Hanson LG, Thielscher A. 2021. Sensitivity and resolution improvement for in vivo magnetic resonance current-density imaging of the human brain. Magnetic Resonance in Medicine. 86(6):3131-3146. https://doi.org/10.1002/mrm.28944

Eroğlu HH, Puonti O, Göksu C, Gregersen F, Siebner HR, Hanson LG, Thielscher A. 2021. On the reconstruction of magnetic resonance current density images of the human brain: Pitfalls and perspectives. NeuroImage. 243:1-15. https://doi.org/10.1016/j.neuroimage.2021.118517

Dubbioso R, Madsen KH, Thielscher A, Siebner HR. 2021. The myelin content of the human precentral hand knob reflects interindividual differences in manual motor control at the physiological and behavioral level. The Journal of Neuroscience: the official journal of the Society for Neuroscience. 41(14):3163-3179. https://doi.org/10.1523/JNEUROSCI.0390-20.2021

Antonenko D, Grittner U, Puonti O, Flöel A, Thielscher A. 2021. Estimation of individually induced e-field strength during transcranial electric stimulation using the head circumference. Brain Stimulation. 14(5):1055-1058. https://doi.org/10.1016/j.brs.2021.07.001

Antonenko D, Grittner U, Saturnino G, Nierhaus T, Thielscher A, Flöel A. 2021. Inter-individual and age-dependent variability in simulated electric fields induced by conventional transcranial electrical stimulation. NeuroImage. 224:1-9. https://doi.org/10.1016/j.neuroimage.2020.117413

Weise K, Numssen O, Thielscher A, Hartwigsen G, Knösche TR. 2020. A novel approach to localize cortical TMS effects. NeuroImage. 209:1-17. Available from: 10.1016/j.neuroimage.2019.116486

Puonti O, Van Leemput K, Saturnino GB, Siebner HR, Madsen KH, Thielscher A. 2020. Accurate and robust whole-head segmentation from magnetic resonance images for individualized head modeling. NeuroImage. 219:1-17. Available from: 10.1016/j.neuroimage.2020.117044

Puonti O, Saturnino GB, Madsen KH, Thielscher A. 2020. Value and limitations of intracranial recordings for validating electric field modeling for transcranial brain stimulation. NeuroImage. 208:1-14. Available from: 10.1016/j.neuroimage.2019.116431

Pasquinelli C, Montanaro H, Lee HJ, Hanson LG, Kim H, Kuster N, Siebner HR, Neufeld E, Thielscher A. 2020. Transducer modeling for accurate acoustic simulations of transcranial focused ultrasound stimulation. Journal of Neural Engineering. 17(4):1-22. Available from: 10.1088/1741-2552/ab98dc

Habich A, Fehér KD, Antonenko D, Boraxbekk C-J, Flöel A, Nissen C, Siebner HR, Thielscher A, Klöppel S. 2020. Stimulating aged brains with transcranial direct current stimulation: Opportunities and challenges: Opportunities and challenges. Psychiatry Research - Neuroimaging. 306:1-9. Available from: 10.1016/j.pscychresns.2020.111179

Boayue NM, Csifcsák G, Aslaksen P, Turi Z, Antal A, Groot J, Hawkins GE, Forstmann B, Opitz A, Thielscher A, Mittner M. 2020. Increasing propensity to mind-wander by transcranial direct current stimulation? A registered report. European Journal of Neuroscience. 51(3):755-780. Available from: 10.1111/ejn.14347

Bikson M, Hanlon CA, Woods AJ, Gillick BT, Charvet L, Lamm C, Madeo G, Holczer A, Almeida J, Antal A, Ay MR, Baeken C, Blumberger DM, Campanella S, Camprodon J, Christiansen L, Colleen L, Crinion J, Fitzgerald P, Gallimberti L, Ghobadi-Azbari P, Ghodratitoostani I, Grabner R, Hartwigsen G, Hirata A, Kirton A, Knotkova H, Krupitsky E, Marangolo P, Nakamura-Palacios EM, Potok W, Praharaj SK, Ruff CC, Schlaug G, Siebner HR, Stagg CJ, Thielscher A, Wenderoth N, Yuan T-F, Zhang X, Ekhtiari H. 2020. Guidelines for TMS/tES Clinical Services and Research through the COVID-19 Pandemic. Brain Stimulation. 13(4):1124-1149. Available from: 10.1016/j.brs.2020.05.010

Saturnino GB, Siebner HR, Thielscher A, Madsen KH Accessibility of cortical regions to focal TES: Dependence on spatial position, safety, and practical constraints. Neuroimage. 2019 doi: 10.1016/j.neuroimage.2019.116183

Saturnino GB, Madsen KH, Thielscher A Electric field simulations for transcranial brain stimulation using FEM: an efficient implementation and error analysis. J Neural Eng. doi: 10.1088/1741-2552/ab41ba, 2019

Pasquinelli C, Hanson LG, Siebner HR, Lee HJ, Thielscher A Safety of transcranial focused ultrasound stimulation: A systematic review of the state of knowledge from both human and animal studies Brain Stimul. doi: 10.1016/j.brs.2019.07.024, 2019

Korshøj AR, Sørensen JCH, von Oettingen G, Poulsen FR, Thielscher A Optimization of tumor treating fields using singular value decomposition and minimization of field anisotropy. Phys Med Biol. 64(4):04NT03. 2019

Saturnino GB, Thielscher A, Madsen KH, Knösche TR, Weise K. A principled approach to conductivity uncertainty analysis in electric field calculations. Neuroimage. 188:821-834, 2019

Karadas M, Wojciechowski AM, Huck A, Dalby NO, Andersen UL, Thielscher A. Feasibility and resolution limits of opto-magnetic imaging of neural network activity in brain slices using color centers in diamond. Sci Rep. 8(1):4503, 2018.

Nielsen JD, Madsen KH, Puonti O, Siebner HR, Bauer C, Madsen CG, Saturnino GB, Thielscher A Automatic skull segmentation from MR images for realistic volume conductor models of the head: Assessment of the state-of-the-art. Neuroimage. 174:587-598, 2018.

Göksu, C., Hanson, L. G., Siebner, H. R., Ehses, P., Scheffler, K. & Thielscher, A. Human in-vivo brain magnetic resonance current density imaging (MRCDI). NeuroImage. 171, p. 26-39, 2018.

Göksu C, Scheffler K, Ehses P, Hanson L.G, Thielscher A. Sensitivity Analysis of Magnetic Field Measurements for Magnetic Resonance Electrical Impedance Tomography (MREIT), Magnetic Resonance in Medicine. 79, p. 748-760, 2018.

Bungert, A., Antunes, A., Espenhahn, S. & Thielscher, A. Where does TMS Stimulate the Motor Cortex? Combining Electrophysiological Measurements and Realistic Field Estimates to Reveal the Affected Cortex Position. Cerebral cortex, 27(11):5083-5094, 2017.

Minjoli, S., Saturnino, G. B., Blicher, J. U., Stagg, C. J., Siebner, H. R., Antunes, A. & Thielscher, A. The impact of large structural brain changes in chronic stroke patients on the electric field caused by transcranial brain stimulation. NeuroImage. Clinical. 15, p. 106-117, 2017.

Saturnino, G. B., Madsen, K. H., Siebner, H. R. & Thielscher, A. How to target inter-regional phase synchronization with dual-site Transcranial Alternating Current Stimulation.
NeuroImage. 163, p. 68-80, 2017.

Madsen, K.H., Ewald, L., Siebner, H.R., Thielscher, A. Transcranial Magnetic Stimulation: An Automated Procedure to Obtain Coil-specific Models for Field Calculations. Brain Stimulation 8, 1205-1208, 2015

Saturnino, G.B., Antunes, A., Thielscher, A. On the importance of electrode parameters for shaping electric field patterns generated by tDCS. Neuroimage 120, 25-35, 2015.

Moisa, M., Siebner, H.R., Pohmann, R., Thielscher, A. Uncovering a context-specific connectional fingerprint of human dorsal premotor cortex. J Neurosci 32, 7244-7252, 2012.

Thielscher, A., Opitz, A., Windhoff, M. Impact of the gyral geometry on the electric field induced by transcranial magnetic stimulation. Neuroimage 54, 234-243, 2011

Group Members

Axel Thielscher

Group Leader

Oula Puonti

Jesper Duemose Nielsen

Show all group members (8)

External Collaborators

Professor Ulrik Lund Andersen 

DTU Physics


Assoc. Professor Alexander Huck

DTU Physics


Assoc. Prof. Koen van Lempuut

Harvard University & DTU Compute


Assoc. Prof. Gottfried Schlaug

Beth Israel Deaconess Medical Center


Prof. Dr. Klaus Scheffler

Max Planck Institut for Biological Cybernetics


Assistant Prof. Hyunjoo Jenny Lee

KAIST


Professor Michael Siniatchkin

University of Bielefeld