ITRUSST consensus on standardised reporting for transcranial ultrasound stimulation

As transcranial ultrasound stimulation (TUS) advances as a precise, non-invasive neuromodulatory method, there is a need for consistent reporting standards to enable comparison and reproducibility across studies. To this end, the International Transcranial Ultrasonic Stimulation Safety and Standards Consortium (ITRUSST) formed a subcommittee of experts across several domains to review and suggest standardised reporting parameters for low intensity TUS, resulting in the guide presented here. The scope of the guide is limited to reporting the ultrasound aspects of a study. The guide and supplementary material provide a simple checklist covering the reporting of: (1) the transducer and drive system, (2) the drive system settings, (3) the free field acoustic parameters, (4) the pulse timing parameters, (5) in situ estimates of exposure parameters in the brain, and (6) intensity parameters. Detailed explanations for each of the parameters, including discussions on assumptions, measurements, and calculations, are also provided.


Introduction
Transcranial 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.
To 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 [1][2][3].
This 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 [4].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.
The 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 https://doi.org/10.1016/j.brs.2024.04.013Received 7 February 2024; Received in revised form 30 March 2024; Accepted 18 April 2024 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.

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

Transducer and drive system description
□ Transducer manufacturer and model number (Section 'Transducer Description') □ Transducer centre frequency (Section 'Transducer description') □ Transducer geometry (e.g., radius of curvature and aperture diameter) (Section 'Transducer description') □ Drive system components, including manufacturer and model number (e.g., signal generator and amplifier or integrated driving system) (Section 'Drive system description') Drive system settings (Section 'Drive system settings') Pulse timing parameters (Section 'Pulse timing parameters')

□ Pulse timing table
In situ estimates of exposure parameters □ Estimated in situ spatial-peak pressure amplitude (Section 'Estimated in situ pressure amplitude') □ Estimated in situ pressure amplitude at the target (Section 'Estimated in situ pressure amplitude') □ Estimated in situ mechanical index (Section 'Estimate of in situ mechanical index') □ One of the following thermal metrics: temperature rise, thermal index, or thermal dose (Section 'Thermal metrics') □ Description of how in situ estimates were obtained (Section 'Estimated in situ pressure amplitude')

System and free field acoustic pressure parameters
System 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 [5] (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.
Reporting 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.

Transducer description
A 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.

Drive system description
The 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.

Drive system settings
All 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.

Parameters to report
The 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.
The 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.

Methods of obtaining free field pressure parameters
The 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.
For 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 [5,8,9] with additional help from the literature (e.g.[10][11][12]).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.[5,[13][14][15]).

Spatial-peak pressure at study output levels
The 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.

Defining the spatial location of reported free field pressure parameters
The 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.

Uncertainty on free-field acoustic parameters
All 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.

Pulse timing parameters
Typical pulse timing parameters are laid out in Fig. 2. The period refers to the duration of one cycle of the operating frequency  0 .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   = 1   .In 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.
The 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.
A 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''.Fig. 2. Schematics of two ultrasound pressure waveforms.A rectangular ramp shape is used in A, while a Tukey window on the pressure waveform is used in B. The pressure amplitude is shown as , while the peak negative pressure is shown as   ; these should be the same when operating at low pressures within the linear regime.The duration of a single cycle of the operating frequency is the period (T).Pulses are often repeated in a pulse train.The time between two pulses in a pulse train is the pulse repetition interval (PRI) and is equal to 1 divided by the pulse repetition frequency (PRF).The pulse train has a duration which is the pulse train duration.The pulse train can be repeated, and if so has a structure similar to the pulse, as laid out in Tables 2.1-2.3.In this figure, the operating frequency  0 is lower than typically used in order to aid visualisation.

Table 2.1
Pulse timing parameters from two experiments in Johnstone et al. [6].The A stimulation is experiment 1, condition A2 and the B stimulation is experiment 4, condition A2.In both cases, the first line contains the pulse characteristics and the second line contains the pulse train characteristics.The 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.
For 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  pulse train .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  pulse train and  pulse train repeat .The overall  is then the product of each .For the patterned rTUS stimulation in Table 2.2 B,  pulse train = 32%,  pulse train repeat = 31.25%,and the product is the overall duty cycle,  = 10%.Note, 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'' [8].This difference is negligible for long pulses.
While the community is only starting to embrace ramping as one means to reduce the auditory confound [6,19,20], 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. [6].
It 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.

Table 2.2
Pulse timing parameters for the two stimulations from Zeng et al. [16].A was a single pulse train repeated every 200 ms for 80 s.B consisted of a short pulse train repeated every 1.6 s, again for a total time of 80 s.This is an example of patterned repetitive TUS.

Derived parameters: In situ estimates of exposure parameters
Exposure 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.

Estimated in situ pressure amplitude
The 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 [21]).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.
(a) Simulating the intracranial ultrasound field.Simulations 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.If 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.
An excellent reference to follow is the Reporting of Computational Modelling Studies in Medical Device Submissions FDA guidance document [13].Examples from the literature can be found in e.g., [22,23].(b) Derating.The 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. [24].
If 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.
If experimental measurements are used to determine derating factors, for example using a sample of human skulls, the methods, assumptions, and variability should be reported [4].This example adapted from Deffieux 2013 [20] 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.(c) Measuring the effect of the intracranial ultrasound field.
In 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 [25].
If 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.
In 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 [14].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.

Estimate of in situ mechanical index
For 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 [26].A detailed discussion can be found in [27].The in situ estimate of mechanical index is easily derived from the in situ estimate of spatial-peak pressure.
The mechanical index is formally defined as where  ,.3 is the peak-rarefactional pressure in MPa derated using an attenuation coefficient of 0.3 dB cm −1 MHz −1 , and  0 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.
In 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: This approach for using in situ, equivalent or effective MI calculations has been used before by several authors [28][29][30].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 [24].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.
The following example calculation uses values taken from Johnstone et al. [6]: The   was 0.44, calculated as follows.The spatialpeak 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 [21]), the derated spatial-peak pressure amplitude was 230 kPa.The operating frequency was 270 kHz.The  tc is then 0.23/sqrt(0.27)= 0.44.

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

(a) Temperature Rise
In 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'.(b) Thermal Index For 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 [32,33] as evidence of safety.Detailed discussions of these parameters can be found in [27].Following IEC 62359 [27], the relevant thermal index is the bone-at-surface or cranial   ( ), which is calculated as Here,  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.
The following example calculation uses values taken from Johnstone et al. [6]: 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 .

(c) Thermal Dose
The 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 [34] defined a 'thermal isoeffective dose' in terms of cumulative equivalent minutes (CEM) at 43 • C which allows conversion of any 2 From 60601-2-37 [31]: ''The TI gives a relative indication of the potential for temperature increase at a specific point along the ultrasound beam.The reason for the term ''relative'' is that the assumed conditions for heating in tissue are complex such that any single index or model cannot be expected to give the actual increase in temperature for all possible conditions and tissue types.Thus, for a particular beam shape, a TI of 2 represents a higher temperature rise than a TI of 1, but does not necessarily represent a rise of 2 • C''.
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  = 0.5 for  ≥ 43 • and  = 0.25 for  < 43 • .This formula accounts for tissue thermotolerance (e.g.mediated by heat-shock proteins) which occurs during exposure at mild hyperthermic temperatures.

Spatial-peak pulse-average intensity
Acoustic 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 [8], they are often used to calculate the spatial-peak intensity in a focused field, where the wave is approximately plane.
Under this assumption, the instantaneous intensity in the direction of the plane wave can be calculated by Here,  sp () 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 × 10 6 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  sp is the spatial-peak pressure amplitude (the amplitude of the sinusoidal pressure signal at the location of the spatial-peak).The  sppa 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.

Spatial-peak time-average intensities
The 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  spta, pulse train =  pulse train  sppa .(10) 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.
The 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.
Ultimately, the maximum continuous length of time where no TUS is administered included in the  spta 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. spta, pulse train and  spta, pulse train repeat .

Neuromodulation dose parameters
For 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 [35].

Summary
In 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.

Declaration of competing interest
The authors declare the following financial interests/personal relationships which may be considered as potential competing interests: EM has acted as a consultant for BrainBox Ltd.JFA holds 5 patents on transcranial ultrasound and had research support from Insightec within the last 3 years.MS is a consultant to Brainsonix Corp., co-inventor on several patents in tFUS, and is a member of the US Technical Advisory Group to the IEC responsible for diagnostic ultrasound output regulation.BET is a developer of and holds a financial interest in the commercially available k-Plan treatment planning software.LV is a member of the scientific board of the Brainbox Initiative.KBP has acted as a consultant for Attune Neurosciences and has had research support from MR Instruments and General Electric within the last 3 years.

□
Operating frequency □ Output level settings □ Focal position settings □ Description of transducer coupling method Free field acoustic parameters (Section 'Parameters to report') Reference position for measurements □ Spatial-peak pressure amplitude □ Position of spatial-peak pressure amplitude (relative to reference position) □ Size of focal volume (−3 dB and −6 dB axial and lateral widths) □ Position of centre of focal volume (centre of −3 dB relative to reference position) □ Description of how free field parameters were obtained (including details of measurement equipment)

Fig. 1 .
Fig. 1.Schematics of example transducers, geometric and field parameters, for information.(A) common transducer configurations including attached coupling media or internal lenses.The external transducer surface plane,   which may be coincident with the front surface of the transducer housing often serves as the reference plane,   from which distances are measured.For some transducers, e.g.complex multi-element devices,   may intersect with the centre of the surface on which the elements sit.(B) Aperture diameter, D, and radius of curvature, R, are used to describe the geometry of focusing transducers.(C) Free field pressure parameters shown with respect to the reference plane.  is the axial position of the location of spatial-peak pressure relative to the reference plane   ;   is the axial position of the centre of the −3 dB focal region relative to the reference plane.