High-resolution motion- and phase-corrected functional MRI at 7 T using shuttered multishot echo-planar imaging.

PURPOSE: To achieve high-resolution multishot echo-planar imaging (EPI) for functional MRI (fMRI) with reduced sensitivity to in-plane motion and between-shot phase variations. METHODS: Two-dimensional radiofrequency pulses were incorporated in a multishot EPI sequence at 7T which selectively excited a set of in-plane bands (shutters) in the phase encoding direction, which moved between shots to cover the entire slice. A phase- and motion-corrected reconstruction was implemented for the acquisition. Brain imaging experiments were performed with instructed motion to evaluate image quality for conventional multishot and shuttered EPI. Temporal stability was assessed in three subjects by quantifying temporal SNR (tSNR) and artifact levels, and fMRI activation experiments using visual stimulation were performed to assess the strength and distribution of activation, using both conventional multishot and shuttered EPI. RESULTS: In the instructed motion experiment, ghosting was lower in shuttered EPI images without or with corrections and image quality metrics were improved with motion correction. tSNR was improved by phase correction in both conventional multishot and shuttered EPI and the acquisitions had similar tSNR without and with phase correction. However, while phase correction was necessary to maximize tSNR in conventional multishot EPI, it also increased intermittent ghosting, but did not increase intermittent ghosting in shuttered EPI. Phase correction increased activation strength in both conventional multishot and shuttered EPI, but caused increased spurious activation outside the brain and in frontal brain regions in conventional multishot EPI. CONCLUSION: Shuttered EPI supports multishot segmented EPI acquisitions with lower sensitivity to artifacts from motion for high-resolution fMRI.

Individually tailored spatial-spectral pulsed CEST MRI for ratiometric mapping of myocardial energetic species at 3T.

PURPOSE: CEST MRI has been used to probe changes in cardiac metabolism via assessment of CEST contrast from Cr. However, B1 variation across the myocardium leads to spatially variable Cr CEST contrast in healthy myocardium. METHODS: We developed a spatial-spectral (SPSP) saturation pulsed CEST protocol to compensate for B1 variation. Flip angle maps were used to individually tailor SPSP pulses comprised of a train of one-dimensional spatially selective subpulses selective along the principal B1 gradient dimension. Complete Z-spectra in the hearts of (n = 10) healthy individuals were acquired using conventional Gaussian saturation and SPSP schemes and supported by phantom studies. RESULTS: In simulations, the use of SPSP pulses reduced the average SD of the effective saturation B1 values within the myocardium (n = 10) from 0.12 ± 0.02 µT to 0.05 ± 0.01 µT (p < 0.01) and reduced the average SD of Cr CEST contrast in vivo from 10.0 ± 4.3% to 6.1 ± 3.5% (p < 0.05). Results from the hearts of human subjects showed a significant reduction of CEST contrast distribution at 2 ppm, as well as amplitude, when using SPSP saturation. Corresponding phantom experiments revealed PCr-specific contrast generation at body temperature when SPSP saturation was used but combined PCr and Cr contrast generation when Gaussian saturation was used. CONCLUSION: The use of SPSP saturation pulsed CEST resulted in PCr-specific contrast generation and enabled ratiometric mapping of PCr to total Cr CEST contrast in the human heart at 3T.

Detection of laser-associated heating in the brain during simultaneous fMRI and optogenetic stimulation.

PURPOSE: To calculate temperatures from T2 *-weighted images collected during optogenetic fMRI based on proton resonance frequency (PRF) shift thermometry, to monitor confounding heating effects and determine appropriate light parameters for optogenetic stimulation. METHODS: fMRI is mainly based on long-TE gradient-recalled echo acquisitions that are also suitable for measuring small temperature changes via the PRF shift. A motion- and respiration-robust processing pipeline was developed to calculate temperature changes based on the PRF shift directly from the T2 *-weighted images collected for fMRI with a two-shot 2D gradient-recalled echo-EPI sequence at 9.4T. Optogenetic fMRI protocols which differed in stimulation durations (3, 6 and 9 s) within a total block duration of 30 s were applied in a squirrel monkey to validate the methods with blue and green light (20 Hz, 30 mW) delivery interleaved between periods. General linear modeling was performed on the resulting time series temperature maps to verify if light delivery with each protocol resulted in significant heating in the brain around the optical fiber. RESULTS: The temperature SD was 0.05°C with the proposed imaging protocol and processing. Statistical analysis showed that the optogenetic stimulation protocol with a 3 s stimulation duration did not result in significant temperature rises. Significant temperature rises up to 0.13°C (p < 0. 05) were observed with 6 and 9 s stimulation durations for blue and green light. CONCLUSION: The proposed processing pipeline can be useful for the design of optogenetic stimulation protocols and for monitoring confounding heating effects.

Selective excitation localized by the Bloch-Siegert shift and a B1+ gradient.

PURPOSE: To perform B1+$$ {B}_1^{+} $$ -selective excitation using the Bloch-Siegert shift for spatial localization. THEORY AND METHODS: A B1+$$ {B}_1^{+} $$ -selective excitation is produced by an radiofrequency (RF) pulse consisting of two summed component pulses: an off-resonant pulse that induces a B1+$$ {B}_1^{+} $$ -dependent Bloch-Siegert frequency shift and a frequency-selective excitation pulse. The passband of the pulse can be tailored by adjusting the frequency content of the frequency-selective pulse, as in conventional B0$$ {B}_0 $$ gradient-localized excitation. Fine magnetization profile control is achieved by using the Shinnar-Le Roux algorithm to design the frequency-selective excitation pulse. Simulations analyzed the pulses' robustness to off-resonance, their suitability for multi-echo spin echo pulse sequences, and how their performance compares to that of rotating-frame selective excitation pulses. The pulses were evaluated experimentally on a 47.5 mT MRI scanner using an RF gradient transmit coil. Multiphoton resonances produced by the pulses were characterized and their distribution across B1+$$ {B}_1^{+} $$ predicted. RESULTS: With correction for varying B1+$$ {B}_1^{+} $$ across the desired profile, the proposed pulses produced selective excitation with the specified profile characteristics. The pulses were robust against off-resonance and RF amplifier distortion, and suitable for multi-echo pulse sequences. Experimental profiles closely matched simulated patterns. CONCLUSION: The Bloch-Siegert shift can be used to perform B0$$ {B}_0 $$ -gradient-free selective excitation, enabling the excitation of slices or slabs in RF gradient-encoded MRI.

Reduced-field of view three-dimensional MR acoustic radiation force imaging with a low-rank reconstruction for targeting transcranial focused ultrasound.

PURPOSE: To rapidly image and localize the focus in MR-guided focused ultrasound (FUS) while maintaining a low ultrasound duty cycle to minimize tissue effects. METHODS: MR-acoustic radiation force imaging (ARFI) is key to targeting FUS procedures such as neuromodulation, and works by encoding ultrasound-induced displacements into the phase of MR images. However, it can require long scan times to cover a volume of tissue, especially when minimizing the FUS dose during targeting is paramount. To simultaneously minimize scan time and the FUS duty cycle, a 2-min three-dimensional (3D) reduced-FOV spin echo ARFI scan with two-dimensional undersampling was implemented at 3T with a FUS duty cycle of 0.85%. The 3D k-space sampling scheme incorporated uniform undersampling in one phase-encoded axis and partial Fourier (PF) sampling in the other. The scan interleaved FUS-on and FUS-off data collection to improve displacement map quality via a joint low-rank image reconstruction. Experiments in agarose and graphite phantoms and living macaque brains for neuromodulation and blood-brain barrier opening studied the effects of the sampling and reconstruction strategy on the acquisition, and evaluated its repeatability and accuracy. RESULTS: In the phantom, the distances between displacement centroids of 10 prospective reconstructions and a fully sampled reference were below 1 mm. In in vivo brain, the distances between centroids ranged from 1.3 to 2.1 mm. Results in phantom and in vivo brain both showed that the proposed method can recover the FUS focus compared to slower fully sampled scans. CONCLUSION: The proposed 3D MR-ARFI reduced-FOV method enables rapid imaging of the FUS focus while maintaining a low FUS duty cycle.

External Dynamic InTerference Estimation and Removal (EDITER) for low field MRI.

PURPOSE: Point-of-care MRI requires operation outside of Faraday shielded rooms normally used to block image-degrading electromagnetic interference (EMI). To address this, we introduce the EDITER method (External Dynamic InTerference Estimation and Removal), an external sensor-based method to retrospectively remove image artifacts from time-varying external interference sources. THEORY AND METHODS: The method acquires data from multiple EMI detectors (tuned receive coils as well as untuned electrodes placed on the body) simultaneously with the primary MR coil during and between image data acquisition. We calculate impulse response functions dynamically that map the data from the detectors to the time varying artifacts then remove the transformed detected EMI from the MR data. Performance of the EDITER algorithm was assessed in phantom and in vivo imaging experiments in an 80 mT portable brain MRI in a controlled EMI environment and with an open 47.5 mT MRI scanner in an uncontrolled EMI setting. RESULTS: In the controlled setting, the effectiveness of the EDITER technique was demonstrated for specific types of introduced EMI sources with up to a 97% reduction of structured EMI and up to 76% reduction of broadband EMI in phantom experiments. In the uncontrolled EMI experiments, we demonstrate EMI reductions of up to 99% using an electrode and pick-up coil in vivo. We demonstrate up to a nine-fold improvement in image SNR with the method. CONCLUSION: The EDITER technique is a flexible and robust method to improve image quality in portable MRI systems with minimal passive shielding and could reduce the reliance of MRI on shielded rooms and allow for truly portable MRI.

On the design and manufacturing of miniaturized microstripline power splitters for driving multicoil transmit arrays with arbitrary ratios at 7 T.

The purpose of the current study was to implement unequal microstrip power splitters for parallel transmission at 7 T that are optimized for size and loss and that can be configured for a wide range of power ratios. The splitters will enable the use of more transmit coils without a corresponding increase in the number of transmit channels or amplifiers to control specific absorption rate, shorten RF pulses, and shim inhomogeneous RF fields. Wilkinson unequal power splitters based on a novel microstrip network design were optimized to minimize their size under 8 cm in length and 9 cm in width, enabling them to be included in coil housing or cascaded in multiple stages. Splitters were designed and constructed for a wide range of output power ratios at 298 MHz. Simulations and bench tests were performed for each ratio, and a methodology was established to adapt the designs to other ratios and frequencies. The designs and code are open source and can be reproduced as is or reconfigured. The single-stage designs achieved good matches and isolations between output ports (worst isolation -15.9 dB, worst match -15.1 dB). A two-stage cascaded (one input to four outputs) power splitter with 1:2.5, 1:10, 1:3, and 1:6 ratio outputs was constructed. The worst isolation between output ports was -19.7 dB in simulation and the worst match of the three ports was -17.8 dB. The measured ratios for one- and two-stage boards were within 10% of the theoretical ratios. The power-handling capability of the smallest trace was approximately 70 W. Power loss for the one- and two-stage boards ranged from 1% to 3% in simulation compared with 5.1% to 7.2% on the bench. It was concluded that Wilkinson unequal microstrip power splitters can be implemented with a small board size (low height) and low loss, and across a wide range of output power ratios. The splitters can be cascaded in multiple stages while maintaining the expected ratios and low loss. This will enable the construction of large fixed transmit array-compression matrices with low loss.

Hybrid-pair ratio adjustable power splitters for add-on RF shimming and array-compressed parallel transmission.

PURPOSE: A ratio adjustable power splitter (RAPS) circuit was recently proposed for add-on RF shimming and array-compressed parallel transmission. Here we propose a new RAPS circuit design based on off-the-shelf components for improved performance and manufacturability. THEORY AND METHODS: The original RAPS used a pair of home-built Wilkinson splitter and hybrid coupler connected by a pair of connectorized coaxial cables. Here we propose a new hybrid-pair RAPS (or HP-RAPS) circuit that replaces the home-built circuits with two commercially available hybrid couplers and replaces connectorized cables with interchangeable microstrip lines. We derive the relation between the desired splitting ratio and the required phase shifts for HP-RAPS and investigate how to generate arbitrary splitting ratios using paired meandering and straight lines. Several HP-RAPSs with different splitting ratios were fabricated and tested on the workbench and MRI experiments. RESULTS: The splitting ratio of an HP-RAPS circuit has a tan or cot dependence on the meandering line's additional length compared to the straight line. The fabricated HP-RAPSs exhibit accurate splitting ratios as expected (<4% deviations) and generate transmit fields that well agree with predicted fields. They also demonstrated a low insertion loss of 0.33 dB, high output isolation of -26 dB, and acceptable impedance matching of -16 dB. CONCLUSION: A novel HP-RAPS circuit was developed and implemented. It is easy-to-fabricate/reproduce with minimal expertise. It also preserves the features of the original RAPS circuit (ratio-adjustable, small footprint, etc.) with lower insertion loss.

k-Space Domain Parallel Transmit Pulse Design.

PURPOSE: To accelerate the design of (under- or oversampled) multidimensional parallel transmission pulses. METHODS: A k-space domain parallel transmission pulse design algorithm was proposed that produces a sparse matrix relating a complex-valued target excitation pattern to the pulses that produce it, and can be finely parallelized. The algorithm was applied in simulations to the design of 3D SPINS pulses for inner volume excitation in the brain at 7 Tesla. It was characterized in terms of the dependence of computation time, excitation error, and required memory on algorithm parameters, and it was compared to an iterative spatial domain pulse design method in terms of computation time, excitation error, Gibbs ringing, and ability to compensate off-resonance. RESULTS: The proposed algorithm achieved approximately 80% faster pulse design compared to the spatial domain method with the same number of parallel threads, with the tradeoff of increased excitation error and RMS RF amplitude. It reduced the memory required to store the design matrix by 99% compared to a full matrix solution. Even with a coarse design grid, the algorithm produced patterns that were free of Gibbs ringing. It was similarly sensitive to k-space undersampling as the spatial domain method, and was similarly capable of compensating for off-resonance. CONCLUSIONS: The proposed k-space domain algorithm accelerates and finely parallelizes parallel transmission pulse design, with a modest tradeoff of excitation error and RMS RF amplitude.

Ultra-high spatial resolution BOLD fMRI in humans using combined segmented-accelerated VFA-FLEET with a recursive RF pulse design.

PURPOSE: To alleviate the spatial encoding limitations of single-shot echo-planar imaging (EPI) by developing multi-shot segmented EPI for ultra-high-resolution functional MRI (fMRI) with reduced ghosting artifacts from subject motion and respiration. THEORY AND METHODS: Segmented EPI can reduce readout duration and reduce acceleration factors, however, the time elapsed between segment acquisitions (on the order of seconds) can result in intermittent ghosting, limiting its use for fMRI. Here, "FLEET" segment ordering, where segments are looped over before slices, was combined with a variable flip angle progression (VFA-FLEET) to improve inter-segment fidelity and maximize signal for fMRI. Scaling a sinc pulse's flip angle for each segment (VFA-FLEET-Sinc) produced inconsistent slice profiles and ghosting, therefore, a recursive Shinnar-Le Roux (SLR) radiofrequency (RF) pulse design was developed (VFA-FLEET-SLR) to generate unique pulses for every segment that together produce consistent slice profiles and signals. RESULTS: The temporal stability of VFA-FLEET-SLR was compared against conventional-segmented EPI and VFA-FLEET-Sinc at 3T and 7T. VFA-FLEET-SLR showed reductions in both intermittent and stable ghosting compared to conventional-segmented and VFA-FLEET-Sinc, resulting in improved image quality with a minor trade-off in temporal SNR. Combining VFA-FLEET-SLR with acceleration, we achieved a 0.6-mm isotropic acquisition at 7T, without zoomed imaging or partial Fourier, demonstrating reliable detection of blood oxygenation level-dependent (BOLD) responses to a visual stimulus. To counteract the increased repetition time from segmentation, simultaneous multi-slice VFA-FLEET-SLR was demonstrated using RF-encoded controlled aliasing. CONCLUSIONS: VFA-FLEET with a recursive RF pulse design supports acquisitions with low levels of artifact and spatial blur, enabling fMRI at previously inaccessible spatial resolutions with a "full-brain" field of view.

Reducing temperature errors in transcranial MR-guided focused ultrasound using a reduced-field-of-view sequence.

PURPOSE: To reduce temperature errors due to water motion in transcranial MR-guided focused ultrasound (tcMRgFUS) ablation. THEORY AND METHODS: In tcMRgFUS, water is circulated in the transducer bowl around the patient's head for acoustic coupling and heat removal. The water moves during sonications that are monitored by MR thermometry, which causes it to alias into the brain and create temperature errors. To reduce these errors, a two-dimensional excitation pulse was implemented in a gradient-recalled echo thermometry sequence. The pulse suppresses water signal by selectively exciting the brain only, which reduces the imaging FOV. Improvements in temperature precision compared to the conventional full-FOV scan were evaluated in healthy subject scans outside the tcMRgFUS system, gel phantom scans in the system with heating, and in 2×-accelerated head phantom scans in the system without heating. RESULTS: In vivo temperature precision (standard deviation of temperature errors) outside the tcMRgFUS system was improved 43% on average, due to the longer TR and TE of the reduced-FOV sequence. In the phantom heating experiments, the hot spot was less distorted in the reduced-FOV scans, and background temperature precision was improved 59% on average. In the accelerated head phantom temperature reconstructions, temperature precision was improved 89% using the reduced-FOV sequence. CONCLUSIONS: Reduced-FOV temperature imaging alleviates temperature errors due to water bath motion in tcMRgFUS, and enables accelerated temperature mapping with greater precision.

Simultaneous multislice MRI thermometry with a single coil using incoherent blipped-controlled aliasing.

PURPOSE: To increase volume coverage in real-time MR thermometry for transcranial MR-guided focused ultrasound (tcMRgFUS) ablation, without multiple receive coils. THEORY AND METHODS: Multiband excitation and incoherent blipped-controlled aliasing were implemented in a 2DFT pulse sequence used clinically for tcMRgFUS, and an extended k-space hybrid reconstruction was developed that recovers slice-separated temperature maps assuming that heating is focal, given slice-separated pretreatment images. Simulations were performed to characterize slice leakage, the number of slices that can be simultaneously imaged with low-temperature error, and robustness across random slice-phase k-space permutations. In vivo experiments were performed using a single receive coil without heating to measure temperature precision, and gel phantom FUS experiments were performed to test the method with heating and with a water bath. RESULTS: Simulations showed that with large hot spots and identical magnitude images on each slice, up to three slices can be simultaneously imaged with less than 1∘ C temperature root-mean-square error. They also showed that hot spots do not alias coherently between slices, and that an average 86% of random slice-phase k-space permutations yielded less than 1∘ C temperature error. Temperature precision was not degraded compared to single-slice imaging in the in vivo SMS scans, and the gel phantom SMS temperature maps closely tracked single-slice temperature in the hot spot, with no coherent aliasing to other slices. CONCLUSIONS: Incoherent controlled aliasing SMS enables accurate reconstruction of focal heating maps from two or three slices simultaneously, using a single receive coil and a sparsity-promoting temperature reconstruction.

Designing parallel transmit head coil arrays based on radiofrequency pulse performance.

PURPOSE: A new approach to design parallel transmit (pTx) head arrays is proposed that integrates transmit radiofrequency pulse designs with electromagnetic modeling of array coil elements. THEORY AND METHODS: An approach to design pTx head arrays is proposed that finds optimal groupings of a large number of coils into a small number of channels. An algorithm is proposed to extend array-compressed parallel transmit pulse design by adding the ability to optimally select and prune coil elements, in addition to optimizing compression weights. The performance of the method is demonstrated in simulations of dynamic multislice shimming of the human brain in axial, coronal, and sagittal directions, and of reduced field-of-view excitation targeting the human occipital lobe, with simulated electromagnetic field maps from a group of 5 human head models at 7T. RESULTS: For both dynamic multislice shimming and reduced field-of-view excitation, the method successfully designed pTx arrays that simultaneously achieved in general 15% lower mean excitation errors with 20% lower SDs, along with 20% lower mean global averaged specific absorption rate and 50% lower SD than previously reported pTx head array designs. CONCLUSION: With the proposed optimal coil element selection algorithm, the array-compressed parallel transmit pulse design can be extended to design pTx transmit head arrays with joint consideration of the fields within the sample and the radiofrequency pulse. The pTx arrays from such an approach achieved higher transmit excitation accuracy, lower radiofrequency heating in subjects, and more robust performance across subjects compared with previously reported pTx head arrays with the same number of channels.

Accelerated spin-echo functional MRI using multisection excitation by simultaneous spin-echo interleaving (MESSI) with complex-encoded generalized slice dithered enhanced resolution (cgSlider) simultaneous multislice echo-planar imaging.

PURPOSE: Spin-echo functional MRI (SE-fMRI) has the potential to improve spatial specificity when compared with gradient-echo fMRI. However, high spatiotemporal resolution SE-fMRI with large slice-coverage is challenging as SE-fMRI requires a long echo time to generate blood oxygenation level-dependent (BOLD) contrast, leading to long repetition times. The aim of this work is to develop an acquisition method that enhances the slice-coverage of SE-fMRI at high spatiotemporal resolution. THEORY AND METHODS: An acquisition scheme was developed entitled multisection excitation by simultaneous spin-echo interleaving (MESSI) with complex-encoded generalized slice dithered enhanced resolution (cgSlider). MESSI uses the dead-time during the long echo time by interleaving the excitation and readout of 2 slices to enable 2× slice-acceleration, while cgSlider uses the stable temporal background phase in SE-fMRI to encode/decode 2 adjacent slices simultaneously with a "phase-constrained" reconstruction method. The proposed cgSlider-MESSI was also combined with simultaneous multislice (SMS) to achieve further slice-acceleration. This combined approach was used to achieve 1.5-mm isotropic whole-brain SE-fMRI with a temporal resolution of 1.5 s and was evaluated using sensory stimulation and breath-hold tasks at 3T. RESULTS: Compared with conventional SE-SMS, cgSlider-MESSI-SMS provides 4-fold increase in slice-coverage for the same repetition time, with comparable temporal signal-to-noise ratio. Corresponding fMRI activation from cgSlider-MESSI-SMS for both fMRI tasks were consistent with those from conventional SE-SMS. Overall, cgSlider-MESSI-SMS achieved a 32× encoding-acceleration by combining Rinplane × MB × cgSlider × MESSI = 4 × 2 × 2 × 2. CONCLUSION: High-quality, high-resolution whole-brain SE-fMRI was acquired at a short repetition time using cgSlider-MESSI-SMS. This method should be beneficial for high spatiotemporal resolution SE-fMRI studies requiring whole-brain coverage.

High-fidelity, high-isotropic-resolution diffusion imaging through gSlider acquisition with B1+ and T1 corrections and integrated ΔB0 /Rx shim array.

PURPOSE: B1+ and T1 corrections and dynamic multicoil shimming approaches were proposed to improve the fidelity of high-isotropic-resolution generalized slice-dithered enhanced resolution (gSlider) diffusion imaging. METHODS: An extended reconstruction incorporating B1+ inhomogeneity and T1 recovery information was developed to mitigate slab-boundary artifacts in short-repetition time (TR) gSlider acquisitions. Slab-by-slab dynamic B0 shimming using a multicoil integrated ΔB0 /Rx shim array and high in-plane acceleration (Rinplane = 4) achieved with virtual-coil GRAPPA were also incorporated into a 1-mm isotropic resolution gSlider acquisition/reconstruction framework to achieve a significant reduction in geometric distortion compared to single-shot echo planar imaging (EPI). RESULTS: The slab-boundary artifacts were alleviated by the proposed B1+ and T1 corrections compared to the standard gSlider reconstruction pipeline for short-TR acquisitions. Dynamic shimming provided >50% reduction in geometric distortion compared to conventional global second-order shimming. One-millimeter isotropic resolution diffusion data show that the typically problematic temporal and frontal lobes of the brain can be imaged with high geometric fidelity using dynamic shimming. CONCLUSIONS: The proposed B1+ and T1 corrections and local-field control substantially improved the fidelity of high-isotropic-resolution diffusion imaging, with reduced slab-boundary artifacts and geometric distortion compared to conventional gSlider acquisition and reconstruction. This enabled high-fidelity whole-brain 1-mm isotropic diffusion imaging with 64 diffusion directions in 20 min using a 3T clinical scanner.

Identifying the Role of Block Length in Neural Heat Block to Reduce Temperatures During Infrared Neural Inhibition.

BACKGROUND AND OBJECTIVES: The objective of this study is to assess the hypothesis that the length of axon heated, defined here as block length (BL), affects the temperature required for thermal inhibition of action potential propagation applied using laser heating. The presence of such a phenomenon has implications for how this technique, called infrared neural inhibition (INI), may be applied in a clinically safe manner since it suggests that temperatures required for therapy may be reduced through the proper spatial application of light. Here, we validate the presence of this phenomenon by assessing how the peak temperatures during INI are reduced when two different BLs are applied using irradiation from either one or two adjacent optical fibers. STUDY DESIGN/MATERIALS AND METHODS: Assessment of the role of BL was carried out over two phases. First, a computational proof of concept was performed in the neural conduction simulation environment, NEURON, simulating the response of action potentials to increased temperatures applied at different full-width at half-maxima (FWHM) along axons. Second, ex vivo validation of these predictions was performed by measuring the radiant exposure, peak temperature rise, and FWHM of heat distributions associated with INI from one or two adjacent optical fibers. Electrophysiological assessment of radiant exposures at inhibition threshold were carried out in ex vivo Aplysia californica (sea slug) pleural abdominal nerves ( n = 6), an invertebrate with unmyelinated axons. Measurement of the maximum temperature rise required for induced heat block was performed in a water bath using a fine wire thermocouple. Finally, magnetic resonance thermometry (MRT) was performed on a nerve immersed in saline to assess the elevated temperature distribution at these radiant exposures. RESULTS: Computational modeling in NEURON provided a theoretical proof of concept that the BL is an important factor contributing to the peak temperature required during neural heat block, predicting a 11.7% reduction in temperature rise when the FWHM along an axon is increased by 42.9%. Experimental validation showed that, when using two adjacent fibers instead of one, a 38.5 ± 2.2% (mean ± standard error of the mean) reduction in radiant exposure per pulse per fiber threshold at the fiber output (P = 7.3E-6) is measured, resulting in a reduction in peak temperature rise under each fiber of 23.5 ± 2.1% ( P = 9.3E-5) and 15.0 ± 2.4% ( P = 1.4E-3) and an increase in the FWHM of heating by 37.7 ± 6.4% ( P = 1E-3), 68.4 ± 5.2% ( P = 2.4E-5), and 51.9 ± 9.9% ( P = 1.7E-3) in three MRT slices. CONCLUSIONS: This study provides the first experimental evidence for a phenomenon during the heat block in which the temperature for inhibition is dependent on the BL. While more work is needed to further reduce the temperature during INI, the results highlight that spatial application of the temperature rise during INI must be considered. Optimized implementation of INI may leverage this cellular response to provide optical modulation of neural signals with lower temperatures over greater time periods, which may increase the utility of the technique for laboratory and clinical use. Lasers Surg. Med. © 2019 Wiley Periodicals, Inc.

Rapid quantitative imaging of high intensity ultrasonic pressure fields.

High intensity focused ultrasound (FUS) is a noninvasive technique for treatment of tissues that can lie deep within the body. There is a need for methods to rapidly and quantitatively map FUS pressure beams for quality assurance and accelerate development of FUS systems and techniques. However, conventional ultrasound pressure beam mapping instruments, including hydrophones and optical techniques, are slow, not portable, and expensive, and most cannot map beams at actual therapeutic pressure levels. Here, a rapid projection imaging method to quantitatively map FUS pressure beams based on continuous-wave background-oriented schlieren (CW-BOS) imaging is reported. The method requires only a water tank, a background pattern, and a camera and uses a multi-layer deep neural network to reconstruct two-dimensional root-mean-square (RMS) projected pressure maps that resolve the ultrasound propagation dimension and one lateral dimension. In this work, the method was applied to collect beam maps over a 3 × 1 cm2 field-of-view with 0.425 mm resolution for focal pressures up to 9 MPa. Results at two frequencies and comparisons to hydrophone measurements show that CW-BOS imaging produces high-resolution quantitative RMS projected FUS pressure maps in under 10 s, the technique is linear and robust to beam rotations and translations, and it can map aberrated beams.

Low-rank plus sparse compressed sensing for accelerated proton resonance frequency shift MR temperature imaging.

PURPOSE: To improve multichannel compressed sensing (CS) reconstruction for MR proton resonance frequency (PRF) shift thermography, with application to MRI-induced RF heating evaluation and MR guided high intensity focused ultrasound (MRgFUS) temperature monitoring. METHODS: A new compressed sensing reconstruction is proposed that enforces joint low rank and sparsity of complex difference domain PRF data between post heating and baseline images. Validations were performed on 4 retrospectively undersampled dynamic data sets in PRF applications, by comparing the proposed method to a previously described L1 and total variation- (TV-) based CS approach that also operates on complex difference domain data, and to a conventional low rank plus sparse (L+S) separation-based CS reconstruction applied to the original domain data. RESULTS: In all 4 retrospective validations, the proposed reconstruction method outperformed the conventional L+S and L1 +TV CS reconstruction methods with a 3.6× acceleration ratio in terms of temperature accuracy with respect to fully sampled data. For RF heating evaluation, the proposed method achieved RMS error of 12%, compared to 19% for the L+S method and 17% for the L1 +TV method. For in vivo MRgFUS thalamotomy, the peak temperature reconstruction errors were 19%, 31%, and 35%, respectively. CONCLUSION: The complex difference-based low rank and sparse model enhances compressibility for dynamic PRF temperature imaging applications. The proposed multichannel CS reconstruction method enables high acceleration factors for PRF applications including RF heating evaluation and MRgFUS sonication.

Multi-echo MR thermometry using iterative separation of baseline water and fat images.

PURPOSE: To perform multi-echo water/fat separated proton resonance frequency (PRF)-shift temperature mapping. METHODS: State-of-the-art, iterative multi-echo water/fat separation algorithms produce high-quality water and fat images in the absence of heating but are not suitable for real-time imaging due to their long compute times and potential errors in heated regions. Existing fat-referenced PRF-shift temperature reconstruction methods partially address these limitations but do not address motion or large time-varying and spatially inhomogeneous B0 shifts. We describe a model-based temperature reconstruction method that overcomes these limitations by fitting a library of separated water and fat images measured before heating directly to multi-echo data measured during heating, while accounting for the PRF shift with temperature. RESULTS: Simulations in a mixed water/fat phantom with focal heating showed that the proposed algorithm reconstructed more accurate temperature maps in mixed tissues compared to a fat-referenced thermometry method. In a porcine phantom experiment with focused ultrasound heating at 1.5 Tesla, temperature maps were accurate to within 1∘ C of fiber optic probe temperature measurements and were calculated in 0.47 s per time point. Free-breathing breast and liver imaging experiments demonstrated motion and off-resonance compensation. The algorithm can also accurately reconstruct water/fat separated temperature maps from a single echo during heating. CONCLUSIONS: The proposed model-based water/fat separated algorithm produces accurate PRF-shift temperature maps in mixed water and fat tissues in the presence of spatiotemporally varying off-resonance and motion.

Volumetric MRI thermometry using a three-dimensional stack-of-stars echo-planar imaging pulse sequence.

PURPOSE: To measure temperature over a large brain volume with fine spatiotemporal resolution. METHODS: A three-dimensional stack-of-stars echo-planar imaging sequence combining echo-planar imaging and radial sampling with golden angle spacing was implemented at 3T for proton resonance frequency-shift temperature imaging. The sequence acquires a 188x188x43 image matrix with 1.5x1.5x2.75 mm3 spatial resolution. Temperature maps were reconstructed using sensitivity encoding (SENSE) image reconstruction followed by the image domain hybrid method, and using the k-space hybrid method. In vivo temperature maps were acquired without heating to measure temperature precision in the brain, and in a phantom during high-intensity focused ultrasound sonication. RESULTS: In vivo temperature standard deviation was less than 1°C at dynamic scan times down to 0.75 s. For a given frame rate, scanning at a minimum repetition time (TR) with minimum acceleration yielded the lowest standard deviation. With frame rates around 3 s, the scan was tolerant to a small number of receive coils, and temperature standard deviation was 48% higher than a standard two-dimensional Fourier transform temperature mapping scan, but provided whole-brain coverage. Phantom temperature maps with no visible aliasing were produced for dynamic scan times as short as 0.38 s. k-Space hybrid reconstructions were more tolerant to acceleration. CONCLUSION: Three-dimensional stack-of-stars echo-planar imaging temperature mapping provides volumetric brain coverage and fine spatiotemporal resolution. Magn Reson Med 79:2003-2013, 2018. © 2017 International Society for Magnetic Resonance in Medicine.