Pulse Inversion Research Articles

Generalized Bloch model: A theory for pulsed magnetization transfer.

The paper introduces a classical model to describe the dynamics of large spin-1/2 ensembles associated with nuclei bound in large molecule structures, commonly referred to as the semi-solid spin pool, and their magnetization transfer (MT) to spins of nuclei in water. Like quantum-mechanical descriptions of spin dynamics and like the original Bloch equations, but unlike existing MT models, the proposed model is based on the algebra of angular momentum in the sense that it explicitly models the rotations induced by radiofrequency (RF) pulses. It generalizes the original Bloch model to non-exponential decays, which are, for example, observed for semi-solid spin pools. The combination of rotations with non-exponential decays is facilitated by describing the latter as Green's functions, comprised in an integro-differential equation. Our model describes the data of an inversion-recovery magnetization-transfer experiment with varying durations of the inversion pulse substantially better than established models. We made this observation for all measured data, but in particular for pulse durations smaller than 300μs. Furthermore, we provide a linear approximation of the generalized Bloch model that reduces the simulation time by approximately a factor 15,000, enabling simulation of the spin dynamics caused by a rectangular RF-pulse in roughly 2μs. The proposed theory unifies the original Bloch model, Henkelman's steady-state theory for MT, and the commonly assumed rotation induced by hard pulses (i.e., strong and infinitesimally short applications of RF-fields) and describes experimental data better than previous models.

Assessment of measurement precision in single‐voxel spectroscopy at 7 T: Toward minimal detectable changes of metabolite concentrations in the human brain in vivo

To introduce a study design and statistical analysis framework to assess the repeatability, reproducibility, and minimal detectable changes (MDCs) of metabolite concentrations determined by in vivo MRS. An unbalanced nested study design was chosen to acquire in vivo MRS data within different repeatability and reproducibility scenarios. A spin-echo, full-intensity acquired localized (SPECIAL) sequence was employed at 7 T utlizing three different inversion pulses: a hyperbolic secant (HS), a gradient offset independent adiabaticity (GOIA), and a wideband, uniform rate, smooth truncation (WURST) pulse. Metabolite concentrations, Cramér-Rao lower bounds (CRLBs) and coefficients of variation (CVs) were calculated. Both Bland-Altman analysis and a restricted maximum-likelihood estimation (REML) analysis were performed to estimate the different variance contributions of the repeatability and reproducibility of the measured concentration. A Bland-Altmann analysis of the spectral shape was performed to assess the variance of the spectral shape, independent of quantification model influences. For the used setup, minimal detectable changes of brain metabolite concentrations were found to be between 0.40 µmol/g and 2.23 µmol/g. CRLBs account for only 16 % to 74 % of the total variance of the metabolite concentrations. The application of gradient-modulated inversion pulses in SPECIAL led to slightly improved repeatability, but overall reproducibility appeared to be limited by differences in positioning, calibration, and other day-to-day variations throughout different sessions. A framework is introduced to estimate the precision of metabolite concentrations obtained by MRS in vivo, and the minimal detectable changes for 13 metabolite concentrations measured at 7 T using SPECIAL are obtained.

Evaluation of Nonlinear Contrast Pulse Sequencing for Use in Super-Resolution Ultrasound Imaging.

The use of super-resolution ultrasound (SR-US) imaging greatly improves visualization of microvascular structures, but clinical adoption is limited by long imaging times. This method depends on detecting and localizing isolated microbubbles (MBs), forcing the use of a dilute contrast agent concentration. Contrast-enhanced ultrasound (CEUS) image acquisition times as long as minutes arise as the localization of thousands of MBs are acquired to form a complete SR-US image. In this article, we explore the use of nonlinear CEUS strategies using nonlinear fundamental contrast pulse sequencing (CPS) to increase the contrast-to-tissue ratio (CTR) and compare MB detection effectiveness to linear B-mode CEUS imaging. The CPS compositions of amplitude modulation (AM), pulse inversion (PI), and a combination of the two (AMPI) were studied. A simulation study combined the Rayleigh-Plesset-Marmottant (RPM) model of MB characteristics and a nonlinear tissue model using the k-Wave toolbox for MATLAB (MathWorks Inc., Natick, MA, USA). Validation was conducted using an in vitro flow phantom and in vivo in the rat hind-limb. Imaging was performed with a programmable US scanner (Vantage 256, Verasonics Inc., Kirkland, WA, USA) and customized to transmit a set of basis US pulses from which both B-mode US (frame rate (FR) of 800 Hz) and multiple nonlinear CPS compositions (FR of 200 Hz) could be assessed from identical in vitro and in vivo datasets using a near simultaneous method. The simulations suggest that MB characteristics, such as diameter and motion, help to predict which US imaging strategy will enhance MB detection. The in vitro and in vivo US imaging studies revealed that different subpopulations of polydisperse MB contrast agents were detected by linear imaging and by each different nonlinear CPS composition. The most effective single imaging strategy at a 200-Hz FR was found to be B-mode US imaging. However, a combination of B-mode US imaging with a nonlinear CPS imaging strategy was more effective in detecting MBs in vivo at all depths and was shown to shorten image acquisition time by an average of 33.3%-76.7% when combining one or more CPS sequences.

Assessment of renal function using magnetic resonance quantitative histogram analysis based on spatial labeling with multiple inversion pulses.

BackgroundThe incidence of chronic kidney disease (CKD) is high, and is easy to develop into end-stage renal disease (ESRD), which requires kidney dialysis or kidney transplantation. Therefore, we want to explore the clinical value of magnetic resonance quantitative histogram analysis based on spatial labeling with multiple inversion pulses (SLEEK) in assessing renal function in the early stage.MethodsOne hundred and twenty-nine patients underwent abdominal MRI examination, including a coronal SLEEK sequence. The patients were divided into the control group [CG, 47 cases, estimated glomerular filtration rate (eGFR) >90], the mild renal function impairment (mRI) group (48 cases, eGFR =60–90), and the moderate to severe renal function impairment (m-sRI) group (34 cases, eGFR <60). Two experienced radiologists delineated cortex and medulla regions of interest (ROIs) on SLEEK images to obtain cortex and medulla quantitative histogram parameters [Mean, Median, Percentiles (5th, 10th, 25th, 75th, and 90th), Skewness, Kurtosis, and Entropy] using FireVoxel. These histogram parameters were compared by proper statistical methods such as one-way analysis of variance, the χ2 test, and receiver operating characteristic (ROC) curve analysis.ResultsFour histogram parameters (Inhomogeneitycortex, Skewnesscortex, Kurtosismedulla, and Entropymedulla) differed significantly between the CG and the mRI group. One medulla (Entropymedulla) and nine cortex (Meancortex, Mediancortex, Kurtosiscortex, Entropycortex, and 5th, 10th, 25th, 75th, and 90th Percentilecortex) histogram parameters were significantly different between the m-RI and m-sRI groups. The most relevant parameter to eGFR was Inhomogenitycortex (r=−0.450, P<0.001). Inhomogeneitycortex had the largest area under the curve (AUC) for differentiating the mRI group from the CG (AUC =0.718; 95% CI: 0.616–0.806), while 25th Percentilecortex generated the largest AUC (AUC =0.786; 95% CI: 0.681–0.869) for differentiating the mRI and m-sRI groups.ConclusionsQuantitative histogram parameters based on a SLEEK sequence can be used to supplement renal dysfunction assessment. Cortex histogram parameters are more valuable for evaluating renal function than medulla histogram parameters.

Efficient solvent suppression with adiabatic inversion for 1H-detected solid-state NMR.

This study introduces a conceptually new solvent suppression scheme with adiabatic inversion pulses for 1H-detected multidimensional solid-state NMR (SSNMR) of biomolecules and other systems, which is termed "Solvent suppression of Liquid signal with Adiabatic Pulse" (SLAP). 1H-detected 2D 13C/1H SSNMR data of uniformly 13C- and 15N-labeled GB1 sample using ultra-fast magic angle spinning at a spinning rate of 60kHz demonstrated that the SLAP scheme showed up to 3.5-fold better solvent suppression performance over a traditional solvent-suppression scheme for SSNMR, MISSISSIPPI (Zhou and Rienstra, J Magn Reson 192:167-172, 2008) with 2/3 of the average RF power.

Characterizing hydrogen bonds in intact RNA from MS2 bacteriophage using magic angle spinning NMR

RNA is a polymer with pivotal functions in many biological processes. RNA structure determination is thus a vital step toward understanding its function. The secondary structure of RNA is stabilized by hydrogen bonds formed between nucleotide basepairs, and it defines the positions and shapes of functional stem-loops, internal loops, bulges, and other functional and structural elements. In this work, we present a methodology for studying large intact RNA biomolecules using homonuclear 15N solid-state NMR spectroscopy. We show that proton-driven spin-diffusion experiments with long mixing times, up to 16 s, improved by the incorporation of multiple rotor-synchronous 1H inversion pulses (termed radio-frequency dipolar recoupling pulses), reveal key hydrogen-bond contacts. In the full-length RNA isolated from MS2 phage, we observed strong and dominant contributions of guanine-cytosine Watson-Crick basepairs, and beyond these common interactions, we observe a significant contribution of the guanine-uracil wobble basepairs. Moreover, we can differentiate basepaired and non-basepaired nitrogen atoms. Using the improved technique facilitates characterization of hydrogen-bond types in intact large-scale RNA using solid-state NMR. It can be highly useful to guide secondary structure prediction techniques and possibly structure determination methods.

Dynamic Filtering of Adherent and Non-adherent Microbubble Signals Using Singular Value Thresholding and Normalized Singular Spectrum Area Techniques

Ultrasound molecular imaging techniques rely on the separation and identification of three types of signals: static tissue, adherent microbubbles and non-adherent microbubbles. In this study, the image filtering techniques of singular value thresholding (SVT) and normalized singular spectrum area (NSSA) were combined to isolate and identify vascular endothelial growth factor receptor 2-targeted microbubbles in a mouse hindlimb tumor model (n = 24). By use of a Verasonics Vantage 256 imaging system with an L12-5 transducer, a custom-programmed pulse inversion sequence employing synthetic aperture virtual source element imaging was used to collect contrast images of mouse tumors perfused with microbubbles. SVT was used to suppress static tissue signals by 9.6 dB while retaining adherent and non-adherent microbubble signals. NSSA was used to classify microbubble signals as adherent or non-adherent with high accuracy (receiver operating characteristic area under the curve [ROC AUC] = 0.97), matching the classification performance of differential targeted enhancement. The combined SVT + NSSA filtering method also outperformed differential targeted enhancement in differentiating MB signals from all other signals (ROC AUC = 0.89) without necessitating destruction of the contrast agent. The results from this study indicate that SVT and NSSA can be used to automatically segment and classify contrast signals. This filtering method with potential real-time capability could be used in future diagnostic settings to improve workflow and speed the clinical uptake of ultrasound molecular imaging techniques.

Xenon Diffusion in Ionic Liquids with Blurred Nanodomain Separation.

The dynamics of xenon gas, loaded in a series of 1-alkyl-3-methylimidazolium based ionic liquids, probes the formation of increasingly blurred polar/apolar nanodomains as a function of the anion type and the cation chain length. Exploiting 129 Xe NMR spectroscopy techniques, like Pulse Gradient Spin Echo (PGSE) and inversion recovery (IR), the diffusion motion and relaxation times are determined for 1-alkyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide [Cn C1 im][TFSI]. A correlation between the ILs nano-structure and both xenon diffusivity and relaxation times, as well as chemical shifts, is outlined. Interestingly, comparison with previous results of the same properties in the homologous imidazolium chlorides and hexafluorophospate shows an opposite trend with the alkyl chain length. Classical molecular dynamics (MD) simulations are used to calculate the xenon and cation and anion diffusion coefficients in the same systems, including imidazolium cations with longer chains (n=4, 6, 8 … 20). An almost quantitative agreement with the experiments validates the MD simulations and, at the same time, provides the necessary structural and dynamic microscopic insights on the nano-segregation and diffusion of xenon in bistriflimide, chloride and hexafluorphosphate salts allowing to observe and rationalize the shaping effect of the cation in the nanostructure.

Cerebrospinal fluid dynamics

The classical cerebrospinal fluid (CSF) circulation theory has been accepted as an established theory of CSF physiology. It describes bulk CSF flow from production site to absorption site. However, much controversy remains regarding the basic CSF physiology and the mechanisms behind the development of hydrocephalus. In the recent observations made using advanced magnetic resonance imaging (MRI) technique, namely, the time spatial inversion pulse (Time-SLIP) method, CSF was used as internal CSF tracer to trace true CSF movement. Observation of the CSF dynamics using this method reveals aspects of CSF dynamics that are different from those of classical CSF circulation theory. Cerebrospinal fluid shows pulsation but does not show bulk flow from production site to absorption site, a theory that was built upon externally injected tracer studies. Observation of the exogeneous tracer studies were true but misinterpreted. Causes of misinterpretations are the differences between results obtained using the true CSF tracer and exogenous tracers. A better understanding of the real CSF physiology can be significant for the advancement of medical sciences in the future. Revisiting CSF flow physiology is a necessary step toward this goal.

Evaluation of metallic bonded plates with nonlinear ultrasound and comparison with destructive testing

In the last decades, the use of structural adhesives has increased. Indeed, they are able to replace advantageously traditional assembly techniques such as riveting or bolting, and generally tend to reduce the global weight of the structure and/or improve the strength of the structure thanks to more homogeneous stress distribution. In addition, it allows the assembly of composite materials or materials having different nature, which is of great interest in the aeronautic industry. However, to be used for structural joining and critical application, reliable non-destructive testing techniques are compulsory for evident safety reasons, be it after fabrication or during the whole life of the structure. While linear ultrasounds have demonstrated their capability to detect easily decohesion or voids in a structure, their use for bond strength inspection is less straightforward. Another approach relies on the analysis of nonlinear signature of a bond defect inspected by high amplitude ultrasound. In the present paper, a Chaotic Cavity Transducer is used to inspect various metallic bonded plates (titanium or aluminum) with several bond defects. The defects were introduced by depositing localized surface pollution (introduction of PTFE spray, release agent, or fingerprints) before the adhesive is deposited on top of the surface. Combined with the pulse inversion technique, high amplitude plane waves were sent in the structure, leading to harmonic components observed in the defect region. Mechanical destructive tests were performed and display a good agreement with nondestructive tests.

Can 7T MPRAGE match MP2RAGE for gray-white matter contrast?

Ultra-High Field (UHF) MRI provides a significant increase in Signal-to-Noise Ratio (SNR) and gains in contrast weighting in several functional and structural acquisitions. Unfortunately, an increase in field strength also induces non-uniformities in the transmit field (B1+) that can compromise image contrast non-uniformly. The MPRAGE is one of the most common T1 weighted (T1w) image acquisitions for structural imaging. It provides excellent contrast between gray and white matter and is widely used for brain segmentation. At 7T, the signal non-uniformities tend to complicate this and therefore, the self-bias-field corrected MP2RAGE is often used there. In both MPRAGE and MP2RAGE, more homogeneous image contrast can be achieved with adiabatic pulses, like the TR-FOCI inversion pulse, or special pulse design on parallel transmission systems, like Universal Pulses (UP). In the present study, we investigate different strategies to improve the bias-field for MPRAGE at 7T, comparing the contrast and GM/WM segmentability against MP2RAGE. The higher temporal efficiency of MPRAGE combined with the potential of the user-friendly UPs was the primary motivation for this MPRAGE-MP2RAGE comparison. We acquired MPRAGE data in six volunteers, adding a k-space shutter to reduce scan time, a kt-point UP approach for homogeneous signal excitation, and a TR-FOCI pulse for homogeneous inversion. Our results show remarkable signal contrast improvement throughout the brain, including regions of low B1+ such as the cerebellum. The improvements in the MPRAGE were largest following the introduction of the UPs. In addition to the CNR, both SNR and GM/WM segmentability were also assessed. Among the MPRAGEs, the combined strategy (UP + TR-FOCI) yielded highest SNR and showed highest spatial similarity between GM segments to the MP2RAGE. Interestingly, the distance between gray and white matter peaks in the intensity histograms did not increase, as better pulses and higher SNR especially benefitted the (cerebellar) gray matter. Overall, the gray-white matter contrast from MP2RAGE is higher, with higher CNR and higher intensity peak distances, even when scaled to scan time. Hence, the extra acquisition time for MP2RAGE is justified by the improved segmentability.

Evaluation of cerebrospinal fluid flow by phase-contrast technique

Objective The objective of this study was to evaluate cerebrospinal fluid (CSF) flow by phase-contrast technique. Background Multiple CSF flow disorders are present like normal pressure hydrocephalus (NPH), atrophic brain changes, communicating hydrocephalus, aqueductal stenosis, and others. Routine MRI as well as new techniques like steady-state free precession sequence, 2D phase-contrast cine MRI, and 2D time spatial labeling inversion pulse can be used to evaluate these CSF flow disorders. Patients and methods This was a case–control study and included 50 patients and 10 healthy volunteers as a control group. All were referred to the radiodiagnosis departments at Menoufia University Hospital and Tanta University Hospital from July 2017 to August 2019. Results This study revealed that CSF flowmetry can differentiate between NPH and atrophic brain changes by the value of stroke volume, where the normal range was 10–21 μl, and all cases of atrophic brain changes showed stroke volume less than 10 μl, where all cases of NPH showed stroke volume more than 21 μl. This study can also confirm the diagnosis of aqueductal stenosis and differentiate it from communicating hydrocephalus (cases of aqueductal stenosis showed stroke volume less than 10 μl, where communicating hydrocephalus showed stroke volume is between 10 and 1 μl). It was also useful in follow-up of patients with endoscopic third ventriculostomy by detection of stroke volume through the third ventricular floor stoma. Conclusion The study showed usefulness of this investigation for confirmation of diagnosis of various CSF disorders. Cine phase-contrast MRI studies, however, should be interpreted in conjunction with the clinical manifestations and the conventional MRI findings in cases of ventriculomegaly. The technique enables measuring aqueductal flow parameters in a reliable and reproducible way and predicts response to shunting and also allows accurate follow-up results.

Assessment of cortical bone fatigue using coded nonlinear ultrasound* *Project supported by the China Postdoctoral Science Foundation (Grant No. 2021M690709), the National Natural Science Foundation of China (Grant Nos. 11827808, 11874289, 11804056, and 12034005), the Program of Shanghai Academic Research Leader (Grant No. 19XD1400500), and the Project of Shanghai Science and Technology Innovation Plan (Grant No. 19441903400).

Bone fatigue accumulation is a factor leading to bone fracture, which is a progressive process of microdamage deteriorating under long-term and repeated stress. Since the microdamage of the early stage in bone is difficult to be investigated by linear ultrasound, the second harmonic generation method in nonlinear ultrasound technique is employed in this paper, which is proved to be more sensitive to microdamage. To solve the deficiency that the second harmonic component is easily submerged by noise in traditional nonlinear measurement, a weighted chirp coded sinusoidal signal was applied as the ultrasonic excitation, while pulse inversion is implemented at the receiving side. The effectiveness of this combination to improve the signal-to-noise ratio has been demonstrated by in vitro experiment. Progressive fatigue loading experiments were conducted on the cortical bone plate in vitro for microdamage generation. There was a significant increase in the slope of the acoustic nonlinearity parameter with the propagation distance (increased by 8% and 24% respectively) when the bone specimen was at a progressive level of microdamage. These results indicate that the coded nonlinear ultrasonic method might have the potential in diagnosing bone fatigue.

The Impact of Vendor-Specific Ultrasound Beam-Forming and Processing Techniques on the Visualization of In Vitro Experimental “Scar”: Implications for Myocardial Scar Imaging Using Two-Dimensional and Three-Dimensional Echocardiography

Myocardial scar appears brighter compared with normal myocardium on echocardiography because of differences in tissue characteristics. The aim of this study was to test how different ultrasound pulse characteristics affect the brightness contrast (i.e., contrast ratio [CR]) between tissues of different acoustic properties, as well as the accuracy of assessing tissue volume. An experimental invitro "scar" model was created using overheated and raw pieces of commercially available bovine muscle. Two-dimensional and three-dimensional ultrasound scanning of the model was performed using combinations of ultrasound pulse characteristics: ultrasound frequency, harmonics, pulse amplitude, steady pulse (SP) emission, power modulation (PM), and pulse inversion modalities. On both two-dimensional and three-dimensional imaging, the CR between the "scar" and its adjacent tissue was higher when PM was used. PM, as well as SP ultrasound imaging, provided good "scar" volume quantification. When tested on 10 "scars" of different size and shape, PM resulted in lower bias (-9.7 vs 54.2mm3) and narrower limits of agreement (-168.6 to 149.2mm3 vs -296.0 to 404.4mm3, P=.03). The interobserver variability for "scar" volume was better with PM (intraclass correlation coefficient=0.901 vs 0.815). Two-dimensional and three-dimensional echocardiography with PM and SP was performed on 15 individuals with myocardial scar secondary to infarction. The CR was higher on PM imaging. Using cardiac magnetic resonance as a reference, quantification of myocardial scar volume showed better agreement when PM was used (bias, -645mm3; limits of agreement, -3,158 to 1,868mm3) as opposed to SP (bias, -1,138mm3; limits of agreement, -5,510 to 3,233mm3). The PM modality increased the CR between tissues with different acoustic properties in an experimental invitro "scar" model while allowing accurate quantification of "scar" volume. By applying the invitro findings to humans, PM resulted in higher CR between scarred and healthy myocardium, providing better scar volume quantification than SP compared with cardiac magnetic resonance.

Joint Design of RF and Gradient Waveforms via Auto-differentiation for 3D Tailored Excitation in MRI.

This paper proposes a new method for joint design of radiofrequency (RF) and gradient waveforms in Magnetic Resonance Imaging (MRI), and applies it to the design of 3D spatially tailored saturation and inversion pulses. The joint design of both waveforms is characterized by the ODE Bloch equations, to which there is no known direct solution. Existing approaches therefore typically rely on simplified problem formulations based on, e.g., the small-tip approximation or constraining the gradient waveforms to particular shapes, and often apply only to specific objective functions for a narrow set of design goals (e.g., ignoring hardware constraints). This paper develops and exploits an auto-differentiable Bloch simulator to directly compute Jacobians of the (Bloch-simulated) excitation pattern with respect to RF and gradient waveforms. This approach is compatible with arbitrary sub-differentiable loss functions, and optimizes the RF and gradients directly without restricting the waveform shapes. For computational efficiency, we derive and implement explicit Bloch simulator Jacobians (approximately halving computation time and memory usage). To enforce hardware limits (peak RF, gradient, and slew rate), we use a change of variables that makes the 3D pulse design problem effectively unconstrained; we then optimize the resulting problem directly using the proposed auto-differentiation framework. We demonstrate our approach with two kinds of 3D excitation pulses that cannot be easily designed with conventional approaches: Outer-volume saturation (90° flip angle), and inner-volume inversion.

Selection on Golay complementary sequences in binary pulse compression for microbubble detection

In medical ultrasound imaging using microbubbles (MBs), the nonlinear echoes from the MBs are used for contrast-specific image construction. Techniques such as pulse inversion, amplitude modulation, and the combined method (PIAM) are employed to increase nonlinear components in the echoes. In addition, employment of pulse compression using binary-coded ultrasound can potentially increase the components. In the case of the nonlinear echo, however, a nonlinear sidelobe occurs around the compressed pulse (correlation peak), and nonlinear components for PIAM in the correlation peak are reduced. The shape of the nonlinear sidelobe and the nonlinear-component reduction in the correlation peak can be estimated from the binary codes used. In this study, the optimal binary codes for PIAM are determined from all patterns of 10, 16, and 20 bit Golay codes. Then, the performance of PIAM with pulse compression using each code is evaluated via computer simulations and experiments using the SonoVue MBs.

Three-dimensional whole-brain mapping of cerebral blood volume and venous cerebral blood volume using Fourier transform-based velocity-selective pulse trains.

To develop 3D MRI methods for cerebral blood volume (CBV) and venous cerebral blood volume (vCBV) estimation with whole-brain coverage using Fourier transform-based velocity-selective (FT-VS) pulse trains. For CBV measurement, FT-VS saturation pulse trains were used to suppress static tissue, whereas CSF contamination was corrected voxel-by-voxel using a multi-readout acquisition and a fast CSF T2 scan. The vCBV mapping was achieved by inserting an arterial-nulling module that included a FT-VS inversion pulse train. Using these methods, CBV and vCBV maps were obtained on 6 healthy volunteers at 3 T. The mean CBV and vCBV values in gray matter and white matter in different areas of the brain showed high correlation (r = 0.95 and P < .0001). The averaged CBV and vCBV values of the whole brain were 5.4 ± 0.6 mL/100 g and 2.5 ± 0.3 mL/100 g in gray matter, and 2.6 ± 0.5 mL/100 g and 1.5 ± 0.2 mL/100 g in white matter, respectively, comparable to the literature. The feasibility of FT-VS-based CBV and vCBV estimation was demonstrated for 3D acquisition with large spatial coverage.

Magnetic resonance angiography and perfusion mapping by arterial spin labeling using Fourier transform-based velocity-selective pulse trains: Examination on a commercial perfusion phantom.

Purpose:Benchmarking of flow and perfusion MR techniques on standardized phantoms can facilitate the use of advanced angiography and perfusion-mapping techniques across multiple sites, field strength, and vendors. Here, MRA and perfusion mapping by arterial spin labeling (ASL) using Fourier transform (FT)–based velocity-selective saturation and inversion pulse trains were evaluated on a commercial perfusion phantom.Methods:The FT velocity-selective saturation–based MRA and FT velocity-selective inversion–based ASL perfusion imaging were compared with time-of-flight and pseudo-continuous ASL at 3 T on the perfusion phantom at two controlled flow rates, 175 mL/min and 350 mL/min. Velocity-selective MRA (VSMRA) and velocity-selective ASL (VSASL) were each performed with three velocity-encoding directions: foot–head, left–right, and oblique 45°. The contrast-to-noise ratio for MRA scans and perfusion-weighted signal, as well as labeling efficiency for ASL methods, were quantified.Results:On this phantom with feeding tubes having only vertical and transverse flow directions, VSMRA and VSASL exhibited the dependence of velocity-encoding directions. The foot–head-encoded VSMRA and VSASL generated similar signal contrasts as time of flight and pseudo-continuous ASL for the two flow rates, respectively. The oblique 45°–encoded VSMRA yielded more uniform contrast-to-noise ratio across slices than foot–head and left–right-encoded VSMRA scans. The oblique 45°–encoded VSASL elevated labeling efficiency from 0.22–0.68 to 0.82–0.90 through more uniform labeling of the entire feeding tubes.Conclusion:Both FT velocity-selective saturation–based VSMRA and FT velocity-selective inversion–based VSASL were characterized on a commercial perfusion phantom. Careful selection of velocity-encoding directions along the major vessels is recommended for their applications in various organs.

Hemodynamic study of unenhanced magnetic resonance angiography using spatial labeling with multiple inversion pulses sequence: a phantom study.

This study sought to explore the functional relationship between displayed vascular length and blood suppression inversion time (BSP TI) and flow velocity in a phantom, and to provide a theoretical basis for quantitatively assessing vascular hemodynamic responses using unenhanced magnetic resonance angiography (MRA) and spatial labeling with multiple inversion pulses sequence (SLEEK). A polyethylene catheter was laid in a long rectangular container filled with pork fat. The entrance of the catheter into the container was connected to a high-pressure syringe filled with normal saline. The high-pressure injector flow rates were set at 0.0, 0.2, 0.4, 0.8, 1.2, 1.6, 2.0, and 2.4 mL/s. SLEEK was performed 19 times for each flow rate with parameter BSP TI values of 50, 75, 100, 150, 200, 300, 400, 500, 600, 700, 800, 900, 1,000, 1,100, 1,200, 1,300, 1,400, 1,600, and 1,800 ms. Maximum intensity projection was employed to reconstruct all SLEEK original images to determine the measurements of the displayed vascular lengths. A regression analysis was undertaken to assess the relationship between the displayed vascular lengths and BSP TI values for each flow rate, and to assess the relationship between the displayed vascular lengths and flow rates at each BSP TI. The displayed vascular length had a linear relationship with BSP TI for each flow rate (P<0.05) (R2=0.754, 0.941, 0.988, 0.988, 0.977, 0.966, and 0.982 for flow rates of 0.0, 0.2, 0.4, 0.8, 1.2, 1.6, and 2.0 mL/s, respectively). The displayed vascular length also had a linear relationship with flow rate for each BSP TI value (P<0.05) (R2 =0.914, 0.912, 0.834, 0.989, 0.980, 0.996, 0.992, 0.960, 0.975, 0.979, 0.982, 0.981, 0.976, and 0.993 for BSP TI 50, 75, 100, 150, 200, 300, 400, 500, 600, 700, 800, 900, 1,000, and 1,100 ms, respectively). No significant linear relationship was found between displayed vascular length and flow rate when the BSP TI value was 1,200 ms (P>0.05). Vascular displayed length has a linear relationship to BSP TI for flow ranges from 0.0 to 2.0 mL/s. Vascular displayed length has a linear relationship to flow rate for BSP TI values of 50 to 1,100 ms. Flow rate can be assessed in relation to vascular displayed length.

Numerical Engineering of Robust Adiabatic Operations

Adiabatic operations are powerful tools for robust quantum control in numerous fields of physics, chemistry and quantum information science. The inherent robustness due to adiabaticity can, however, be impaired in applications requiring short evolution times. We present a single versatile gradient-based optimization protocol that combines adiabatic control with effective Hamiltonian engineering in order to design adiabatic operations tailored to the specific imperfections and resources of an experimental setup. The practicality of the protocol is demonstrated by engineering a fast, 2.3 Rabi cycle-long adiabatic inversion pulse for magnetic resonance with built-in robustness to Rabi field inhomogeneities and resonance offsets. The performance and robustness of the pulse is validated in a nanoscale force-detected magnetic resonance experiment on a solid-state sample, indicating an ensemble-averaged inversion accuracy of $\sim 99.997\%$. We further showcase the utility of our protocol by providing examples of adiabatic pulses robust to spin-spin interactions, parameter-selective operations and operations connecting arbitrary states, each motivated by experiments.