ITRUSST consensus on standardised reporting for transcranial ultrasound stimulation

1Transcranial ultrasound stimulation (TUS) is a non-invasive neuromodulation technique that employs focused ultrasound waves to modulate neuronal activity within the brain. TUS offers a promising avenue for therapeutic and research applications due to its spatial precision and ability to target deep neural structures. As the field transitions to more widespread human studies, the time is upon us to standardise the reporting of such studies to aid understanding and reproducibility.

2To this end, the International Transcranial Ultrasonic Stimulation Safety and Standards Consortium (ITRUSST) formed a subcommittee to review and suggest standardised reporting parameters for TUS, with multiple presentations to the ITRUSST group and multiple opportunities for feedback, resulting in the guide presented here. Similar guidelines have been produced for other analogous techniques and applications.

3This guide focuses on the ultrasound aspects of TUS experiments. It does not discuss reporting of additional elements such as EEG, fMRI, or behavioural readouts, which are important, but outside the scope of this work. It does not seek to establish safety values or thresholds, which are addressed in a separate ITRUSST consensus on TUS safety. While this guide endeavours to align with existing international standards related to diagnostic ultrasound imaging and focused ultrasound transducers, it may diverge from these standards when necessary. An important additional note is that we make the assumption that the ultrasound fields that this guide applies to are linear. Our rationale for doing so is provided in Section A of the Supplementary Material.

4The paper is structured as follows. Firstly, we provide a summary checklist of all required reporting parameters in Section ‘Checklist’. In Section ‘System and free field acoustic pressure parameters’, we discuss reporting of the transducer and drive system, and the free field acoustic pressure amplitude and spatial characteristics of the focal region, as measured by a hydrophone in free field in a water bath. These measurements, performed under standardised conditions, serve as a baseline for comparison between studies independent from the specific configuration of use during TUS studies. In Section ‘Pulse timing parameters’, we discuss the pulse timing parameters. Next, in Section ‘Derived parameters: in situ estimates of exposure parameters’, we discuss the reporting of exposure parameters, which describe the acoustic field inside the individual brain, after accounting for the skull bone, brain tissue, and any other acoustic distortions. Further details are provided in the Supplementary Material sections A–D.

5Report these parameters for all devices and settings used. Further details for each parameter are given in the referenced section.

6Transducer manufacturer and model number (Section ‘Transducer Description’)

7Transducer centre frequency (Section ‘Transducer description’)

8Transducer geometry (e.g., radius of curvature and aperture diameter) (Section ‘Transducer description’)

9Drive system components, including manufacturer and model number (e.g., signal generator and amplifier or integrated driving system) (Section ‘Drive system description’)

10Operating frequency

11Output level settings

12Focal position settings

13Description of transducer coupling method

14Reference position for measurements

15Spatial-peak pressure amplitude

16Position of spatial-peak pressure amplitude (relative to reference position)

17Size of focal volume (−3 dB and −6 dB axial and lateral widths)

18Position of centre of focal volume (centre of −3 dB relative to reference position)

19Description of how free field parameters were obtained (including details of measurement equipment)

20Pulse timing table

21Estimated in situ spatial-peak pressure amplitude (Section ‘Estimated in situ pressure amplitude’)

22Estimated in situ pressure amplitude at the target (Section ‘Estimated in situ pressure amplitude’)

23Estimated in situ mechanical index (Section ‘Estimate of in situ mechanical index’)

24One of the following thermal metrics: temperature rise, thermal index, or thermal dose (Section ‘Thermal metrics’)

25Description of how in situ estimates were obtained (Section ‘Estimated in situ pressure amplitude’)

26Spatial-peak pulse-average intensity (Section ‘Spatial-peak pulse-average intensity’)

27Spatial-peak time-average intensities (Section ‘Spatial-peak time-average intensities’)

28The acoustic impedance used for the conversion (Section ‘Spatial-peak pulse-average intensity’)

29System parameters describe the type and geometry of the transducer, and the signal chain or system used to drive the transducer in order to generate the acoustic field. The free field acoustic pressure parameters provide a description of the generated ultrasound field under reference conditions. The term ‘free field’ refers to the field generated when a transducer is radiating continuously (or with the comparatively long pulses used in TUS studies) into water without any obstruction by reflectors, scatterers, or aberrators. See Section B for details of how this is implemented during measurements.

30Reporting free field acoustic parameters allows comparison of the transducer output at the chosen output level and focal settings under standardised conditions, serving as a baseline for comparison between studies. In this guide, the term output level refers to the amplitude of the acoustic pressure output, which is governed by settings on the drive system such as voltage, power or intensity. The focal position setting refers to the position of the focal region where this can be steered to different positions. It may be fixed for a particular transducer, or set by the operator, for example, by selecting a distance setting, or by setting the relative phases applied to the transducer elements. The free field parameters described in this section should be reported at the output and focal settings used during exposure of participants during studies.

31A description of the ultrasound transducer should be given, including the manufacturer and model number, and the operating frequency of the transducer. For spherically focusing transducers, the geometry is described by the radius of curvature and aperture diameter (see Fig. 1). In some cases, the exact geometry may not be known, for example, if the transducer construction includes a lens or permanently attached coupling medium, in which case, any other information such as the nominal position of the focus relative to the transducer face should be reported. For multi-element transducers, the number of elements and information about their shape, size, and positions should be given if available. The exact size and position of elements may not be available for commercial transducers, so a brief description accompanied by the model number is sufficient. Any other features of the transducer such as permanently attached coupling media, or lenses should be described, including a description of the material, its thickness and shape/geometry, and properties of the material such as the sound speed, density and attenuation coefficient where known.

32The input signal used to drive the transducer may be generated by an integrated drive system, or using a signal generator and radio frequency amplifier. All components of the system, including manufacturer and model numbers and any external electrical impedance matching networks used to couple the drive signal to the transducer should be reported.

33All drive system settings used during studies should be described, including the operating frequency, output level settings (e.g. displayed electrical power, focal intensity, voltage etc.), and for multi-element transducers, focal distance or position settings applied. The method of coupling the transducer to the participant forms part of the in situ acoustic transmission path, so the materials and methods used should be described.

34The free field acoustic pressure field parameters are illustrated in Fig. 1. Table 1.1 shows the quantities to be reported, which include spatial-peak pressure amplitude, focal dimensions, and location. An explanation of how to calculate focal dimensions from the axial pressure profile is given in Section B.8. Two illustrative reporting examples are provided in Table 1.1, and annotated example axial pressure profiles are shown in Fig B.2.

35The amplitude of the spatial-peak acoustic pressure should be reported at each of the output level settings, and all quantities should be reported for each of the focal position settings used during studies, where practical. This will provide a reference spatial-peak pressure in water for each study exposure. For multi-element array devices where steering and aberration correction are performed, the focal position and amplitude may vary between participants and target locations. Therefore, reporting of free field parameters for each condition may not be useful. Instead information about the pressure amplitude and lateral and axial width of the focal region over the focal steering range utilised during the study should be reported. See Section C for recommendations for reporting average or ranges of parameters.

36The free field pressure parameters can be obtained by investigators in one of several ways, for example, directly from acoustic field measurements using a hydrophone in a water bath, from a test report provided by the manufacturer or other calibration/characterisation provider, or via a combination of simulation and measurement. In any case, the source and/or method used to obtain the field parameters should be reported, and where measurement or simulation were used, the procedures should be described, as well as the equipment used (including hydrophone manufacturer, model, and element size, and calibration parameters). Further guidance can be found in Section B.

37For those performing their own measurements to determine the free field acoustic pressure parameters, it is recommended that information on measurement best practice be sought from hydrophone measurement standards with additional help from the literature (e.g.). While a full description is out of scope of this document, some further details and discussion are provided in Section B. Similarly, where simulation is used to obtain some or all of the field parameters, the literature should be consulted for methods and best practice (e.g.).

38The free field spatial-peak pressure amplitude at each of the study output level settings and focal position settings may be obtained either directly from a measurement with the transducer operated at the study output level and focal setting, or by scaling a measurement made with the same focal setting, but at a different output level as described in Section B.7. Typically, a lower output level can be used for free field measurements than would be used to obtain a similar pressure amplitude in the brain after propagation through the skull. Measurements with sufficient signal to noise ratio can be made, while reducing the possibility of causing damage to the hydrophone (especially at the lowest frequencies). This is discussed in more detail in Section B.

39The spatial location of the spatial-peak pressure amplitude can be defined relative to a variety of different reference locations, for example the radiating surface or the external transducer surface plane, also known as the exit plane (see Fig. 1). In some cases the full construction of the transducer or the distance from the radiating surface to the outer surface of any permanently attached or integrated lens or coupling medium may not be known. In this case, a useful reference location would be the outer surface of these layers, the external transducer surface plane. The chosen reference plane location should be reported along with the parameters.

40All hydrophone measurements have an associated measurement uncertainty which arises from a combination of systematic and random uncertainties. One major source of systematic uncertainty is hydrophone sensitivity, which can be up to 20% depending on frequency and calibration method. Additionally, uncertainty on pressure measurements will propagate to derived quantities. Understanding of uncertainty on pressure measurements is essential when reporting and comparing measurements, therefore uncertainties should be reported alongside hydrophone measurements. Further discussion of this can be found in Section B.3.

41Typical pulse timing parameters are laid out in Fig. 2. The period refers to the duration of one cycle of the operating frequency. A pulse is a single continuous sonication and has a duration referred to as the pulse duration (PD). If the pulse is repeated, the pulse repetition interval (PRI) is the time between successive pulses. The pulse repetition frequency is given by.

42In Tables 2.1–2.3 we suggest a table reporting structure. In the rare case of only a single continuous pulse, only the first line of a table need be provided, with the pulse duration, and if a ramp is used, the pulse ramp duration and the pulse ramp shape. The pulse repetition interval is not applicable. When the pulse is repeated, the PRI/PRF are entered in the first line, and the second line provides information on the pulse train including the pulse train duration, the pulse train ramp duration, and the pulse train ramp shape.

43The examples chosen from the literature illustrate a number of variations in TUS pulse parameters in order to illustrate increasing levels of patterning, as well as rectangular and ramped waveforms. While figures are helpful and illustrative, we feel that tables of the pulse parameters would be most clear.

44A number of phrases are commonly used to refer to ultrasound pulse trains. We recommend avoiding ambiguous phrases while retaining clear conventions of the community. The use of the word “burst” is ambiguous as it has different meanings in different contexts. For example, in diagnostic ultrasound, burst refers to a single long pulse, while on function generators and in TMS, burst is often used to refer to pulse trains. We suggest avoiding the use of the word “burst” to remove this ambiguity and instead to use the terms “pulse” and “pulse trains”.

45The phrase “repetitive TUS” (rTUS) is conventionally used to refer to pulse trains intended to elicit cumulative or delayed effects, as opposed to only acute effects. The distinction between acute and delayed effects is critical for study design (‘online’ vs. ‘offline’ designs), the underlying neurophysiology (acute modulatory vs. early-phase plasticity mechanisms), and the safety assessment. Distinguishing characteristics of an rTUS protocol, for example, the pulse repetition frequency, can be included in the label, as in the phrase ‘10 Hz rTUS’. It is recommended to avoid using unique labels when an (extended) rTUS label is also appropriate. By definition, all rTUS protocols include at least the pulse and pulse train timing parameters, as in the 5 Hz rTUS example A of Table 2.2. Patterned rTUS protocols are a subset of rTUS where the pulse trains are repeated in a pattern and, thus, include more than two rows, as in the patterned rTUS example B of Table 2.2.

46For rectangular pulses in a single pulse train, duty cycle (DC) is defined as the percentage of time that a pulse is on, Duty cycle can be optionally reported since it can be easily derived from the timing tables and is only defined for rectangular pulses. However, if provided, we recommend referring to it in a particular way. The rTUS stimulation protocol of Table 2.2 A has a single duty cycle of 10%. In this case, the overall DC is equal to the. However, when rTUS is patterned, such as in the stimulation of Table 2.2 B, it would be clearer to refer to the duty cycle of each level of patterning, such as and. The overall is then the product of each. For the patterned rTUS stimulation in Table 2.2 B,,, and the product is the overall duty cycle,.

47Note, our definition of the pulse duration differs slightly from the International Electrotechnical Commission definition, which defines it as “1.25 times the interval between the time when the time integral of the square of the instantaneous acoustic pressure reaches 10% and 90% of its final value”. This difference is negligible for long pulses.

48While the community is only starting to embrace ramping as one means to reduce the auditory confound, it is likely that this will continue and potentially become a source of confusion unless reported in a standardised way. The shape and duration of the ramp must be reported as already described, as well as whether that shape is applied to the pressure (voltage) waveform or the intensity (power) waveform. In Table 2.1, we specify how ramping is applied in Johnstone et al..

49It is likely that the field will continue to evolve in ways that we cannot now anticipate. The table structure allows for evolution by adding columns or rows as needed.

50Exposure parameters describe the properties of the acoustic field inside the individual brain, after accounting for the skull bone, brain tissue, and any other acoustic distortions. The exact ultrasound exposure is generally difficult to directly measure in vivo. Instead, exposure parameters are typically estimated using knowledge of the free field parameters along with a derating procedure, simulations, or indirect measurements. Here, in situ estimate is used as a general term to describe the procedure used to calculate the exposure parameters, i.e. the in vivo or in situ properties of the acoustic field. The recommended in situ parameters and some examples are provided here.

51The skull bone has very different acoustic material properties (sound speed, density, and attenuation) to the surrounding soft tissues. As ultrasound traverses the skull, these differences cause the waves to become distorted and lose energy. Brain tissue also has a higher attenuation coefficient compared to water. Consequently, the acoustic focus inside the brain will generally have a much lower intensity (sometimes by a factor of 10). In some cases, the shape and position of the focus can also become distorted. To complicate matters, the shape and properties of the skull vary significantly both within and between subjects. The same system and pulse timing parameters can therefore generate very different acoustic fields in the brain depending on the subject and the position of the transducer. An estimate of in situ pressure amplitude should thus be reported. Several different approaches for obtaining this are outlined below. In all cases, the method used should be reported, along with details of any calculations, parameters, and assumptions.

52Simulating the intracranial ultrasound field.

53Simulations are increasingly used to calculate the ultrasound and temperature fields inside the skull and brain. The material properties for the simulations might be based on individualised images for each subject or a CT or MR template image.

54If using simulations to calculate the exposure parameters, the following should be reported: Details of the simulation tool.Details of how the acoustic and thermal material properties are assigned.How the transducer modelling and positioning was performed within the simulation.Details of relevant simulation input parameters, program settings including mesh or grid parameters, and processing steps.

55Details of the simulation tool.

56Details of how the acoustic and thermal material properties are assigned.

57How the transducer modelling and positioning was performed within the simulation.

58Details of relevant simulation input parameters, program settings including mesh or grid parameters, and processing steps.

59An excellent reference to follow is the Reporting of Computational Modelling Studies in Medical Device Submissions FDA guidance document. Examples from the literature can be found in e.g.,.

60Derating.

61The simplest approach to calculating exposure parameters is to derate the pressure by applying acoustic attenuation coefficients to the free field parameters. However, one must be careful as a simple derating may not account for skull resonance effects. This is described with suggested skull derating factors for several experimental conditions in Attali et al..

62If using attenuation factors to calculate the in situ exposure parameters, the following should be reported: The individual acoustic attenuation coefficients in dB.(MHz.cm)−1 or dB.cm−1 (for example, for skull, scalp and brain), and the values used for frequency and distance.The total attenuation applied in dB.

63The individual acoustic attenuation coefficients in dB.(MHz.cm)−1 or dB.cm−1 (for example, for skull, scalp and brain), and the values used for frequency and distance.

64The total attenuation applied in dB.

65If experimental measurements are used to determine derating factors, for example using a sample of human skulls, the methods, assumptions, and variability should be reported.

66This example adapted from Deffieux 2013 reports attenuation from only skull in the following way: The pressure amplitude at focus was set to 0.6 MPa, as measured in free water with a heterodyne interferometer. Skull transmission was estimated on a clean and degassed primate skull specimen (Macaca mulatta skull) at seven different locations and was found to be 58% ±8% (derating factor of −4.7 dB). This allowed us to estimate the derated spatial-peak pressure at 0.35 MPa in the brain of the monkeys.

67Measuring the effect of the intracranial ultrasound field.

68In some cases, indirect measurement of the physical effects of the ultrasound field on the brain can be made. For example, a temperature rise induced by the absorption of ultrasound energy could be measured using MR thermometry. Similarly, for some beam shapes, the bulk displacement induced by an acoustic radiation force could be measured using MR-ARFI. With some restrictions, these measurements could be correlated with the ultrasound exposure parameters in the brain.

69If using measurements to calculate the exposure parameters, the following should be reported: Details of the measured quantity and the measurement method.Details of how the exposure parameters are calculated from the measured quantity.

70Details of the measured quantity and the measurement method.

71Details of how the exposure parameters are calculated from the measured quantity.

72In principle, all of the spatial acoustic parameters described in Section ‘System and free field acoustic pressure parameters’ could be reported after accounting for the skull and brain. In particular, if using simulations to calculate exposure parameters, it is possible to extract the focal characteristics as well as pressure values at the spatial-peak and target locations. As a minimum we recommend the following parameters are reported: (1) the in situ estimate for the spatial-peak pressure amplitude and its location, and (2) the in situ estimate for the pressure amplitude at the target, if it is known and different from the spatial-peak pressure amplitude. The former values are related to safety, the latter values are related to the study efficacy.

73For diagnostic ultrasound imaging, the mechanical index provides a standardised indicator related to the potential for mechanical bioeffects, specifically cavitation. In lieu of a safety standard specifically related to TUS, is also often reported in TUS studies in relation to regulatory limits given in the FDA guidance document. A detailed discussion can be found in. The in situ estimate of mechanical index is easily derived from the in situ estimate of spatial-peak pressure.

74The mechanical index is formally defined as where is the peak-rarefactional pressure in MPa derated using an attenuation coefficient of 0.3 dB cm−1 MHz−1, and is the operating frequency in MHz (note units). This index and derating factor is intended to apply to ultrasound propagating in soft tissue. Note, in this expression, is used, as this is the formal definition of MI. Assuming operation in the linear regime, the peak positive and negative pressures are equal, and equal to the pressure amplitude, i.e.,, as used in elsewhere in this guide.

75In TUS applications where ultrasound propagates through the skull, a more suitable derating factor may be used to define an alternative index, specific to this application, called the, where the subscript “tc” refers to transcranial applications:

76This approach for using in situ, equivalent or effective MI calculations has been used before by several authors. In this case, the method used to calculate the in situ or derated peak-rarefactional pressure, (for example, a simulation) should be reported. If derating a free field value,, the attenuation coefficient, insertion loss of the skull, and any methods and parameters used for calculation should be reported, including whether resonances are taken into account as in. If simulation indicates that the target pressure is lower than the peak pressure, then the MI for both the spatial-peak and the target should be given.

77The following example calculation uses values taken from Johnstone et al.: The was 0.44, calculated as follows. The spatial-peak pressure amplitude measured in water was 700 kPa. Assuming a constant derating factor of −9.8 dB (the average insertion loss of the skull at 270 kHz measured in), the derated spatial-peak pressure amplitude was 230 kPa. The operating frequency was 270 kHz. The is then 0.23/sqrt(0.27) = 0.44.

78At least one of the following metrics should be reported.

79Temperature Rise

80In some cases, it may be possible to directly measure the temperature rise in soft tissue due to the applied TUS (for example, using MR thermometry or thermocouples), or to estimate the temperature rise in soft tissue using a thermal simulation. If reporting the temperature rise in soft tissue, the measurement or simulation method should be described as discussed in Section ‘Estimated in situ pressure amplitude’.

81Thermal Index

82For diagnostic ultrasound imaging, the thermal index provides a standardised indicator related to the potential for thermal bioeffects. is a unitless quantity intended to indicate a potential temperature rise in degrees Celsius. While scales with temperature rise, it does not equal the temperature rise.2 In lieu of a safety standard specifically related to TUS, it is also often reported in TUS studies in relation to limits given by the British Medical Ultrasound Society and the American Institute of Ultrasound in Medicine as evidence of safety. Detailed discussions of these parameters can be found in. Following IEC 62359, the relevant thermal index is the bone-at-surface or cranial, which is calculated as

83Here, is the time-averaged acoustic power of the transducer in free field in mW, and is the equivalent aperture diameter of the transducer in cm. For transducers in direct contact with the scalp, the equivalent aperture diameter would be taken to be the nominal aperture diameter of the transducer. However, when TUS transducers are applied with some stand off (such as a coupling pad) between the transducer surface and the scalp, the equivalent aperture diameter should be taken as the beam diameter at the scalp. This can be estimated from the distance between the transducer surface and the scalp, and the geometry of the transducer. Therefore, the transducer to scalp distance should be reported. If reporting TIC, the values for and (and the methods used to measure or calculate them) should also be reported.

84The following example calculation uses values taken from Johnstone et al.: The TIC was 0.48, calculated as follows. The electrical power was 4.8 W. Assuming a nominal electrical efficiency of 85%, this gives an acoustic power of 4.1 W. The minimum pulse train repetition interval was 5 s, giving an overall duty cycle of 3% (150 ms on every 5 s). This gives a time-averaged power of 0.12 W. The nominal aperture diameter was 64 mm. The TIC is then (1000 * 4.8 * 0.85 * 0.03)/(40 * 6.4) = 0.48.

85Thermal Dose

86The extent of biological changes in tissue resulting from thermal exposure is correlated with the amount of energy absorbed in tissue. For thermal energy, it is the temperature to which the tissue is raised, and the duration of the heating that play the predominant biological roles. Sapareto and Dewey defined a ‘thermal isoeffective dose’ in terms of cumulative equivalent minutes (CEM) at 43 °C which allows conversion of any temperature–time (T–t) combination to the equivalent time for which the reference temperature of 43 °C must be applied to obtain the same level of thermal damage, given by where for and for. This formula accounts for tissue thermotolerance (e.g. mediated by heat-shock proteins) which occurs during exposure at mild hyperthermic temperatures.

87Acoustic intensity is a measure of the flow of energy from one point in an acoustic medium to another. The instantaneous acoustic intensity is a vector quantity (it depends on direction) and is given by the product of the acoustic pressure and acoustic particle velocity. For a plane wave, the pressure and particle velocity are related by the characteristic acoustic impedance of the medium (the product of sound speed and density). While plane wave expressions have only limited validity and cannot generally be applied throughout a focused acoustic field, they are often used to calculate the spatial-peak intensity in a focused field, where the wave is approximately plane.

88Under this assumption, the instantaneous intensity in the direction of the plane wave can be calculated by Here, is the time varying acoustic pressure at the location of the spatial-peak, and is the characteristic acoustic impedance of the medium, which is approximately 1.5 × 106 Rayls for soft tissue. For a harmonically-varying acoustic wave, assuming that the pulse duration is an integer multiple of the acoustic period (or otherwise long compared to the acoustic period), the intensity time-averaged over the pulse duration (generally called the pulse average intensity) can be calculated using which, for long pulses without ramping, reduces to where is the spatial-peak pressure amplitude (the amplitude of the sinusoidal pressure signal at the location of the spatial-peak). The should be easily derivable by the reader from the estimated in situ pressure and timing, and is therefore an optional reporting parameter. If reporting intensity values, the characteristic acoustic impedance used for the calculation must be reported.

89The intensity time-averaged over the pulse train or pulse train repetition interval (generally called the time average intensity) can be calculated by where is the time period over which the average is taken. In the case of a pulse train without ramping, this reduces to Intensity parameters are widely reported in TUS studies. However, as mentioned above, it is important to note that in general, acoustic intensity is a vector quantity. Thus, intensity values calculated in this manner should always be reported alongside (not instead of) acoustic pressure values.

90The pulse timing parameters used in TUS frequently include intermittent pulsing, often with different repetition intervals at different pulse levels (see Section ‘Pulse timing parameters’). This causes ambiguity in the definition of time-averaged intensity, specifically regarding which time period should be used for averaging. Clearly, if 5 s pulse trains are administered every hour, taking the time-average over the latter period does not give a meaningful quantity.

91Ultimately, the maximum continuous length of time where no TUS is administered included in the calculation should be determined based on knowledge of characteristic diffusion times relevant to TUS. However, this is still a topic of active research, thus a recommendation cannot yet be made. For standardised reporting, our current recommendation is therefore to report the time period over which the average is taken. This can be done succinctly using the pulse level nomenclature introduced in Section ‘Pulse timing parameters’, e.g. and, pulse train repeat.

92For TUS, the precise mechanisms through which ultrasound affects brain function are only just beginning to be understood. This means it is not yet possible to define or adopt appropriate dose parameters.

93In this paper, we have suggested a minimal set of reporting guidelines for TUS parameters. We provide some examples of how to estimate In Situ parameters. We also provide a checklist of parameters to be reported. In the supplementary material, we provide more details on how to obtain free field measurements. We briefly list other aspects of the TUS experiment that fall outside of the scope of this paper but are also important in fully describing a study.

94Supplementary material related to this article can be found online at https://doi.org/10.1016/j.brs.2024.04.013.